Phân tích hiệu quả giảm chấn của nhiều hệ cản điều chỉnh cột chất lỏng trong kết cấu chịu động đất

I H C QU C GIA TP H CHÍ MINH TR NG I H C BÁCH KHOA TR NG CH QUANG HI U PHÂN TÍCH HI U QU GI M CH N C A NHI U H C N I U CH NH C T CH T L NG TRONG K T C U CH U NG T Chuyên ngành : Xây d ng cơng trình dân d ng công nghi p Mã s ngành : 60 58 02 08 LU N V N TH C S Thành ph H Chí Minh, tháng 01 n m 2020 CƠNG TRÌNH C HỒN THÀNH T I TR NG I H C QU C GIA TP H Cán b h ng d n khoa h c: Cán b h ng d n: PGS.TS Nguy n Tr ng Ph I H C BÁCH KHOA CHÍ MINH c Cán b ch m nh n xét 1: TS Kh ng Tr ng Toàn Cán b ch m nh n xét 2: PGS TS L Lu n v n th c s c b o v t i Tr ng V n H i ng i h c Bách Khoa, HQG Tp.HCM, ngày 09 tháng 01 n m 2020 Thành ph n H i ng ánh giá c ng lu n v n th c s g m: PGS TS Chu Qu c Th ng - Ch t ch H i TS Nguy n Thái Bình - Th ký TS Kh ng Tr ng Toàn - y viên (Ph n bi n 1) PGS TS L - y viên (Ph n bi n 2) ng V n H i - y viên TS Liêu Xuân Quí CH T CH H I ng NG TR NG KHOA K THU T XÂY D NG i ĐẠI HỌC QUỐC GIA TP.HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập - Tự - Hạnh phúc NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: TRƯƠNG CHẾ QUANG HIỆU MSHV: 1770076 Ngày, tháng, năm sinh: 20/01/1994 Nơi sinh: Quảng Ngãi Chuyên ngành: Kỹ Thuật Xây Dựng Cơng Trình Dân Dụng Và Cơng Nghiệp Mã số : 60.58.02.08 I TÊN ĐỀ TÀI: PHÂN TÍCH HIỆU QUẢ GIẢM CHẤN CỦA NHIỀU HỆ CẢN ĐIỀU CHỈNH CỘT CHẤT LỎNG TRONG KẾT CẤU CHỊU ĐỘNG ĐẤT II NHIỆM VỤ VÀ NỘI DUNG: Tìm hiểu mơ hình kết cấu có lắp đặt hệ MTLCD (Multiple Tuned Liquid Column Dampers), lập phương trình chuyển động chủ đạo, viết chương trình tính Phân tích hiệu giảm dao động với việc lắp đặt hệ cản nhiều cột chất lỏng bị động MTLCD không lắp đặt hệ giảm dao động kết cấu chịu tải điều hòa gia tốc động đất dựa chương trình tính tốn viết ngơn ngữ lập trình MATLAB Nhận xét hệ giảm dao động nhiều cột chất lỏng MTLCD III NGÀY GIAO NHIỆM VỤ: 11/02/2019 IV NGÀY HOÀN THÀNH NHIỆM VỤ: 08/12/2019 V CÁN BỘ HƯỚNG DẪN: PGS TS NGUYỄN TRỌNG PHƯỚC TP.HCM, ngày tháng năm CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM BỘ MÔN TRƯỞNG KHOA KỸ THUẬT XÂY DỰNG ii LỜI CẢM ƠN Trước tiên xin bày tỏ lòng biết ơn đến PGS.TS Nguyễn Trọng Phước Với tận tụy nhiệt tình, Thầy giúp phát triển ý tưởng giúp đỡ tơi nhiều việc tìm tài liệu định hướng nghiên cứu khoa học đến lúc hoàn thành luận văn thạc sĩ Những tài liệu tham khảo kiến thức quý báu Thầy mang lại giúp có cách nhận định đắn vấn đề nghiên cứu Bên cạnh tơi xin chân thành cảm ơn quý Thầy Cô Khoa Kỹ Thuật Xây dựng, trường Đại học Bách Khoa TP.HCM truyền dạy kiến thức q giá cho tơi, kiến thức thiếu đường nghiên cứu khoa học nghiệp sau Đồng thời cảm ơn anh chị đồng khóa giúp đỡ tơi q trình học tập suốt thời gian qua Cuối xin cảm ơn tất người thân, gia đình, bạn bè ln bên cạnh động viên khuyến khích tơi suốt thời gian học tập nghiên cứu thực đề tài Luận văn thạc sĩ hoàn thành thời gian quy định với nỗ lực thân, nhiên tránh khỏi thiếu sót Kính mong q Thầy Cô dẫn thêm để bổ sung kiến thức hồn thiện thân Tơi xin trân trọng cảm ơn! iii TÓM TẮT LUẬN VĂN THẠC SĨ Giảm chấn cho kết cấu hệ cản cột cản chất lỏng (Tuned Liquid Column Damper-TLCD) ngày quan tâm năm gần Một mơ hình hệ cản nhiều cột chất lỏng (Multiple Tuned Liquid Column Dampers-MTLCD) đề xuất gần để khắc phục số nhược điểm hệ đơn TLCD Mô hình kết cấu để phân tích luận văn bao gồm nhiều bể chứa chất lỏng hình chữ U đặt tầng chịu tác động tải điều hòa gia tốc động đất Kết cấu mơ tả hệ hữu hạn bậc tự chuyển vị ngang sàn tầng Thông số đặc trưng thiết bị MTLCD phân tích luận văn Phương trình chuyển động hệ chịu tải trọng động thiết lập dựa nguyên lý cân động giải phương pháp tích phân Newmark tồn miền thời gian Một chương trình máy tính viết ngơn ngữ lập trình MATLAB để phân tích ứng xử động lực học hệ Từ kết số khảo sát luận văn đưa kết luận hiệu giảm dao động mơ hình MTLCD so với mơ hình đơn TLCD hệ kết cấu khơng gắn hệ cản Cuối cùng, luận văn xét ảnh hưởng thông số đặc trưng hệ cản đến giảm giao động cho hệ kết cấu nhằm tìm thơng số phù hợp cho làm việc hiệu thiết bị MTLCD iv ABSTRACT Tuned Liquid Column Damper (TLCD) has been increasingly concerned in recent years A model of Multiple Tuned Liquid Column Dampers (MTLCD) has been proposed to overcome some disadvantages of single TLCD The structural model for analysis in this thesis includes many U-shaped tubes located on the top floor that are affected by the harmonic load and earthquake acceleration The main structure described as a multi-degrees of freedom is the horizontal displacement at the floor level Typical parameters of MTLCD equipment are analyzed in the thesis The equation of motion of the dynamic loadbearing system is established based on the dynamic equilibrium theory which is solved by using the Newmark integration method in the time domain A computer program is written in the MATLAB programming language to analyze the dynamic behavior of the whole system From the numerical results surveyed in the thesis, it is concluded that the effect of vibration reduction of MTLCD model compared to a single model TLCD or structural system without barrier Finally, the thesis examines the influence of each characteristic of the barrier system to reduce the vibration of the structural system in order to find suitable parameters for the efficient working of MTLCD equipment v LỜI CAM ĐOAN Tơi xin cam đoan cơng việc thực hướng dẫn thầy PGS.TS Nguyễn Trọng Phước Các kết Luận văn thật chưa công bố nghiên cứu khác Tôi xin chịu trách nhiệm công việc thực Tp HCM, ngày … tháng … năm … Trương Chế Quang Hiệu vi MỤC LỤC NHIỆM VỤ LUẬN VĂN THẠC SĨ i LỜI CẢM ƠN ii TÓM TẮT LUẬN VĂN THẠC SĨ iii ABSTRACT iv LỜI CAM ĐOAN .v MỤC LỤC vi DANH MỤC BẢNG BIỂU ix DANH MỤC HÌNH VẼ x DANH MỤC KÍ HIỆU xiii MỘT SỐ KÍ HIỆU VIẾT TẮT xv GIỚI THIỆU 1.1 Đặt vấn đề 1.2 Mục tiêu luận văn 10 1.3 Phương pháp thực 11 1.4 Cấu trúc luận văn .11 TỔNG QUAN .13 2.1 Giới thiệu chương 13 2.2 Phân loại thiết bị giảm chấn sử dụng chất lỏng 13 Hệ cản điều chỉnh chuyển động sóng chất lỏng bề mặt 14 Hệ cản cột chất lỏng (TLCD) 16 2.3 Ứng dụng hệ cản cột chất lỏng hệ cản bể chất lỏng 18 2.4 Tổng quan tài liệu giảm dao động sử dụng chất lỏng 22 vii Tình hình nghiên cứu nước 22 Tình hình nghiên cứu nước 26 2.5 Kết luận chương 27 CƠ SỞ LÝ THUYẾT 28 3.1 Giới thiệu chương 28 3.2 Nguyên lý hoạt động hệ cản cột chất lỏng 28 Mơ hình hệ bậc tự có gắn TLCD 29 Mơ hình hệ bậc tự có gắn TLCD cải tiến .30 3.3 Nguyên lý hoạt động hệ cản nhiều cột chất lỏng 31 Mơ hình hệ bậc tự có gắn MTLCD 32 Xác định hàm truyền hệ MTLCD 37 Các thông số thiết kế hệ MTLCD 40 Mơ hình hệ nhiều bậc tự có gắn MTLCD 43 3.4 Phương pháp giải thuật toán 47 Phương pháp giải 48 Thuật toán 50 3.5 Kết luận chương 52 VÍ DỤ SỐ 53 4.1 Giới thiệu chương 53 4.2 Kiểm chứng chương trình tính 53 Bài toán chu kỳ dao động dạng dao động riêng .53 Bài toán kiểm tra độ tin cậy thuật toán Newmark 54 Bài toán kết cấu chịu động đất 56 Bài toán kết cấu có gắn TLCD 57 4.3 Phản ứng hệ kết cấu – MTLCD tác động tải trọng điều hòa 59 viii Khảo sát tỷ số khối lượng µi .60 Khảo sát tỷ số tần số β 63 Khảo sát tỷ số điều chỉnh γi 65 Khảo sát tỷ số khối lượng bể τ 68 Khảo sát hệ số kích thước bể α 70 Khảo sát nhiều hệ cản cột chất lỏng có tần số tự nhiên khác .71 4.4 Phản ứng hệ kết cấu – MTLCD chịu động đất El Centro 74 Khảo sát tỷ số khối lượng µi .75 Khảo sát tỷ số khối lượng bể τ 78 Khảo sát hệ số kích thước bể α 80 Khảo sát bước thời gian Δt trận động đất .81 Khảo sát nhiều hệ cản cột chất lỏng có tần số tự nhiên khác .82 4.5 Phản ứng hệ kết cấu – MTLCD chịu động đất ParkField 84 Khảo sát tỷ số khối lượng µi .85 Khảo sát tỷ số khối lượng bể τ 87 Khảo sát hệ số kích thước bể α 89 Khảo sát bước thời gian Δt trận động đất .90 Khảo sát nhiều hệ cản cột chất lỏng có tần số tự nhiên khác .91 4.6 Kết luận chương 92 KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN 94 5.1 KẾT LUẬN 94 5.2 HƯỚNG PHÁT TRIỂN ĐỀ TÀI 96 TÀI LIỆU THAM KHẢO 97 PHỤ LỤC .101 LÝ LỊCH TRÍCH NGANG .116 102 beta=1/6; end %% Phan tich he BTD tim So cua tai dT=stepT(2)-stepT(1); % buoc thoi gian P1=-1*accelerate; % Tai % Dieu kien ban dau U=zeros(length(accelerate),1); % Vecto chuyen vi ban dau U(1)=0; Udot=zeros(length(accelerate),1); % Vecto van toc ban dau Udot(1)=0; U2dot=zeros(length(accelerate),1); % Vecto gia toc ban dau U2dot(1)=0; for i=1:length(accelerate)-1; K_eff1=1/(beta*dT^2)+gama*2*khi1*omega1/(beta*dT)+omega1^2; P_eff1=P1(i+1)+1*(U(i)/(beta*dT^2)+Udot(i)/(beta*dT)+U2dot(i)* (1/2/beta-1))+ 2*khi1*omega1*(U(i)*gama/(beta*dT)+Udot(i)* (gama/beta-1)+dT/2*(gama/beta-2)*U2dot(i)); % Chuyen vi u thu i U(i+1)=P_eff1/K_eff1; Udot(i+1)=Udot(i)*(1-gama/beta)+dT*U2dot(i)*(1-gama/2/beta)+ (U(i+1)-U(i))*gama/beta/dT; U2dot(i+1)=(U(i+1)-U(i))/(beta*dT^2)-Udot(i)/(beta*dT)- (1/2/beta-1)*U2dot(i); end So=rms(U)^2*4*khi1*omega1^3/pi; %% Khong co he can (NoTLCD) dT=stepT(2)-stepT(1); % Buoc thoi gian P=accelerate; % Tai tac dong % Dieu kien ban dau X_noTLCD=zeros(dofs,length(accelerate)); % Vecto chuyen vi ban dau X_noTLCD(:,1)=0; Xdot_noTLCD=zeros(dofs,length(accelerate)); % Vecto van toc ban dau Xdot_noTLCD(:,1)=0; X2dot_noTLCD=zeros(dofs,length(accelerate)); % Vecto gia toc ban dau X2dot_noTLCD(:,1)=0; A=zeros(dofs,1); B=zeros(dofs,1); for i=1:length(accelerate)-1; for k=1:dofs A(k,1)=(-X_noTLCD(k,i)/beta/dT^2-Xdot_noTLCD(k,i)/beta/dT- (1/2/beta-1)*X2dot_noTLCD(k,i)); B(k,1)=(-gama*X_noTLCD(k,i)/beta/dT+(1-gama/beta)* Xdot_noTLCD(k,i)+(1-gama/2/beta)*dT*X2dot_noTLCD(k,i)); end % Chuyen vi u thu i K_eff=(1/beta/dT^2*Ms+gama/beta/dT*Cs+Ks); P_eff=(-ms'*P(i+1)-Ms*A-Cs*B); X_noTLCD(:,i+1)=K_eff\P_eff; Xdot_noTLCD(:,i+1)=Xdot_noTLCD(:,i)*(1-gama/beta)+ dT*X2dot_noTLCD(:,i)*(1-gama/2/beta)+ (X_noTLCD(:,i+1)-X_noTLCD(:,i))*gama/beta/dT; X2dot_noTLCD(:,i+1)=(X_noTLCD(:,i+1)-X_noTLCD(:,i))/ (beta*dT^2)-Xdot_noTLCD(:,i)/(beta*dT)- (1/2/beta-1)*X2dot_noTLCD(:,i); end X17_MAX_noTLCD=max(X_noTLCD(nstory,:))*100 figure('color',[1 1]); 103 plot(stepT,X_noTLCD(dofs,:),':k','linewidth',0.1); %% He co TLCD (1_TLCD) -muy=0.01; % Ti so khoi luong (he can/ket cau chinh) gama1=1/(1+muy); % Ty so dieu chinh (tan so he can/tan so ket cau chinh) anpha=0.9; % Ti so giua B/L cua he can tau=1; % Ti so khoi luong be/khoi luong chat long omega2=omega1*gama1; % Tan so he can hold on khi2=0.01; % Ti so can cua he can omegaL=omega1*gama1; % Tan so tu nhien cua chat long xi_L=100; [xichmaucham]=timcp7_cho_xiL_tim_xichma_u_cham(khi1,muy,anpha,Biendo,omeg a_0,tau,khi2,xi_L,omega1,omega2,omegaL); xichma_u_cham=xichmaucham; Cp=xichma_u_cham*xi_L/sqrt(2*pi); Cp_bd=Cp; % -M1=zeros(dofs+2); M1(dofs+1,dofs)=1; M1(dofs+1,dofs+1)=1; M1(dofs+1,dofs+2)=anpha/(1+tau); M1(dofs+2,dofs)=anpha; M1(dofs+2,dofs+1)=anpha; M1(dofs+2,dofs+2)=1; M2=[Ms zeros(dofs,2)]; M2=[M2;zeros(2,dofs+2)]; M_TLCD=M1+M2; % -C1=zeros(dofs+2); C1(dofs,dofs+1)=-2*khi2*omega2*muy*M_all; C1(dofs+1,dofs+1)=2*khi2*omega2; C1(dofs+2,dofs+2)=2*Cp; C2=[Cs zeros(dofs,2)]; C2=[C2;zeros(2,dofs+2)]; C_TLCD=C1+C2; % -K1=zeros(dofs+2); K1(dofs,dofs+1)=-omega2^2*muy*M_all; K1(dofs+1,dofs+1)=omega2^2; K1(dofs+2,dofs+2)=omegaL^2; K2=[Ks zeros(dofs,2)]; K2=[K2;zeros(2,dofs+2)]; K_TLCD=K1+K2; % -ms_TLCD=[ms';1;anpha]; % -% Initial conditions X_TLCD=zeros(dofs+2,length(accelerate)); % Vecto chuyen vi ban dau X_TLCD(:,1)=0; Xdot_TLCD=zeros(dofs+2,length(accelerate)); % Vecto van toc ban dau Xdot_TLCD(:,1)=0; X2dot_TLCD=zeros(dofs+2,length(accelerate)); % Vecto gia toc ban dau X2dot_TLCD(:,1)=0; A=zeros(dofs+2,1); B=zeros(dofs+2,1); for i=1:length(accelerate)-1; 104 for r=1:dofs+2 A(r,1)=(-X_TLCD(r,i)/beta/dT^2-Xdot_TLCD(r,i)/beta/dT- (1/2/beta-1)*X2dot_TLCD(r,i)); B(r,1)=(-gama*X_TLCD(r,i)/beta/dT+(1-gama/beta)* Xdot_TLCD(r,i)+(1-gama/2/beta)*dT*X2dot_TLCD(r,i)); end % Chuyen vi u thu i K_eff=(1/beta/dT^2*M_TLCD+gama/beta/dT*C_TLCD+K_TLCD); P_eff=(-ms_TLCD*P(i+1)-M_TLCD*A-C_TLCD*B); X_TLCD(:,i+1)=K_eff\P_eff; Xdot_TLCD(:,i+1)=Xdot_TLCD(:,i)*(1-gama/beta)+ dT*X2dot_TLCD(:,i)*(1-gama/2/beta)+ (X_TLCD(:,i+1)-X_TLCD(:,i))*gama/beta/dT; X2dot_TLCD(:,i+1)=(X_TLCD(:,i+1)-X_TLCD(:,i))/(beta*dT^2)- Xdot_TLCD(:,i)/(beta*dT)-(1/2/beta-1)* X2dot_TLCD(:,i); end X17_MAX_TLCD=max(X_TLCD(nstory,:))*100 plot(stepT,X_TLCD(dofs,:),'-g','linewidth',1.5); %% He co TLCD -N=3; % So luong he can muy=0.01; % Ti so khoi luong (he can/ket cau chinh) gama1=1/(1+muy); % Ty so dieu chinh (tan so he can/tan so ket cau chinh) anpha=0.9; % Ti so giua B/L cua he can tau=1; % Ti so khoi luong be/khoi luong chat long omega2=omega1*gama1; % Tan so he can hold on khi2=0.01; % Ti so can cua he can omegaL=omega1*gama1; % Tan so tu nhien cua chat long xi_L=100; [xichmaucham]=f_tim_xichmaucham(N,khi1,muy,anpha,Biendo,omega_0,tau,khi2, xi_L,omega1,omega2,omegaL); xichma_u_cham=xichmaucham; Cp=xichma_u_cham*xi_L/sqrt(2*pi); Cp_bd=Cp; % Ma tran khoi luong M1=zeros(dofs+2*N); for i=dofs+1:dofs+N M1(i,dofs)=muy*M_all; end for i=dofs+N+1:dofs+2*N M1(i,dofs)=anpha; end for i=dofs+1:dofs+N M1(i,i)=muy*M_all; end for i=dofs+N+1:dofs+2*N M1(i,i-N)=anpha; end for i=dofs+1:dofs+N M1(i,i+N)=muy*anpha*M_all/(1+tau); end for i=dofs+N+1:dofs+2*N M1(i,i)=1; end M2=[Ms zeros(dofs,2*N)]; M2=[M2;zeros(2*N,dofs+2*N)]; M_MTLCD=M1+M2; 105 % Ma tran can -C1=zeros(dofs+2*N); for i=dofs+1:dofs+N C1(dofs,i)=-2*khi2*omega2*muy*M_all; end for i=dofs+1:dofs+N C1(i,i)=2*khi2*omega2*muy*M_all; end for i=dofs+N+1:dofs+2*N C1(i,i)=2*Cp; end C2=[Cs zeros(dofs,2*N)]; C2=[C2;zeros(2*N,dofs+2*N)]; C_MTLCD=C1+C2; % Ma tran cung -K1=zeros(dofs+2*N); for i=dofs+1:dofs+N K1(dofs,i)=-omega2^2*muy*M_all; end for i=dofs+1:dofs+N K1(i,i)=omega2^2*muy*M_all; end for i=dofs+N+1:dofs+2*N K1(i,i)=omegaL^2; end K2=[Ks zeros(dofs,2*N)]; K2=[K2;zeros(2*N,dofs+2*N)]; K_MTLCD=K1+K2; % -ms_MTLCD_1=zeros(dofs+2*N,1); for i=dofs+1:dofs+N ms_MTLCD_1(i,1)=muy*M_all; end for i=dofs+N+1:dofs+2*N ms_MTLCD_1(i,1)=anpha; end ms_MTLCD_2=[ms';zeros(2*N,1)]; ms_MTLCD=ms_MTLCD_1+ms_MTLCD_2; % % Dieu kien ban dau X_MTLCD=zeros(dofs+2*N,length(accelerate)); % Vecto chuyen vi ban dau X_MTLCD(:,1)=0; Xdot_MTLCD=zeros(dofs+2*N,length(accelerate)); % Vecto van toc ban dau Xdot_MTLCD(:,1)=0; X2dot_MTLCD=zeros(dofs+2*N,length(accelerate)); % Vecto gia toc ban dau X2dot_MTLCD(:,1)=0; A=zeros(dofs+2*N,1); B=zeros(dofs+2*N,1); for i=1:length(accelerate)-1; for r=1:dofs+2*N A(r,1)=(-X_MTLCD(r,i)/beta/dT^2-Xdot_MTLCD(r,i)/beta/dT- (1/2/beta-1)*X2dot_MTLCD(r,i)); B(r,1)=(-gama*X_MTLCD(r,i)/beta/dT+(1-gama/beta)* Xdot_MTLCD(r,i)+(1-gama/2/beta)*dT*X2dot_MTLCD(r,i)); 106 end % Chuyen vi u thu i K_eff=(1/beta/dT^2*M_MTLCD+gama/beta/dT*C_MTLCD+K_MTLCD); P_eff=(-ms_MTLCD*P(i+1)-M_MTLCD*A-C_MTLCD*B); X_MTLCD(:,i+1)=K_eff\P_eff; Xdot_MTLCD(:,i+1)=Xdot_MTLCD(:,i)*(1-gama/beta)+ dT*X2dot_MTLCD(:,i)*(1-gama/2/beta)+ (X_MTLCD(:,i+1)-X_MTLCD(:,i))*gama/beta/dT; X2dot_MTLCD(:,i+1)=(X_MTLCD(:,i+1)-X_MTLCD(:,i))/(beta*dT^2)- Xdot_MTLCD(:,i)/(beta*dT)-(1/2/beta-1)* X2dot_MTLCD(:,i); end X17_MAX_MTLCD=max(X_MTLCD(nstory,:))*100 plot(stepT,X_MTLCD(dofs,:),'-r','linewidth',1.5); %% He co TLCD -N=5; % So luong he can muy=0.01; % Ti so khoi luong (he can/ket cau chinh) gama1=1/(1+muy); % Ty so dieu chinh (tan so he can/tan so ket cau chinh) anpha=0.9; % Ti so giua B/L cua he can tau=1; % Ti so khoi luong be/khoi luong chat long omega2=omega1*gama1; % Tan so he can hold on khi2=0.01; % Ti so can cua he can omegaL=omega1*gama1; % Tan so tu nhien cua chat long xi_L=100; [xichmaucham]=f_tim_xichmaucham(N,khi1,muy,anpha,Biendo,omega_0,tau,khi2, xi_L,omega1,omega2,omegaL); xichma_u_cham=xichmaucham; Cp=xichma_u_cham*xi_L/sqrt(2*pi); Cp_bd=Cp; % Ma tran khoi luong M1=zeros(dofs+2*N); for i=dofs+1:dofs+N M1(i,dofs)=muy*M_all; end for i=dofs+N+1:dofs+2*N M1(i,dofs)=anpha; end for i=dofs+1:dofs+N M1(i,i)=muy*M_all; end for i=dofs+N+1:dofs+2*N M1(i,i-N)=anpha; end for i=dofs+1:dofs+N M1(i,i+N)=muy*anpha*M_all/(1+tau); end for i=dofs+N+1:dofs+2*N M1(i,i)=1; end M2=[Ms zeros(dofs,2*N)]; M2=[M2;zeros(2*N,dofs+2*N)]; M_MTLCD=M1+M2; % Ma tran can -C1=zeros(dofs+2*N); for i=dofs+1:dofs+N C1(dofs,i)=-2*khi2*omega2*muy*M_all; end 107 for i=dofs+1:dofs+N C1(i,i)=2*khi2*omega2*muy*M_all; end for i=dofs+N+1:dofs+2*N C1(i,i)=2*Cp; end C2=[Cs zeros(dofs,2*N)]; C2=[C2;zeros(2*N,dofs+2*N)]; C_MTLCD=C1+C2; % Ma tran cung -K1=zeros(dofs+2*N); for i=dofs+1:dofs+N K1(dofs,i)=-omega2^2*muy*M_all; end for i=dofs+1:dofs+N K1(i,i)=omega2^2*muy*M_all; end for i=dofs+N+1:dofs+2*N K1(i,i)=omegaL^2; end K2=[Ks zeros(dofs,2*N)]; K2=[K2;zeros(2*N,dofs+2*N)]; K_MTLCD=K1+K2; % -ms_MTLCD_1=zeros(dofs+2*N,1); for i=dofs+1:dofs+N ms_MTLCD_1(i,1)=muy*M_all; end for i=dofs+N+1:dofs+2*N ms_MTLCD_1(i,1)=anpha; end ms_MTLCD_2=[ms';zeros(2*N,1)]; ms_MTLCD=ms_MTLCD_1+ms_MTLCD_2; % -% Dieu kien ban dau X_MTLCD=zeros(dofs+2*N,length(accelerate)); % Vecto chuyen vi ban dau X_MTLCD(:,1)=0; Xdot_MTLCD=zeros(dofs+2*N,length(accelerate)); % Vecto van toc ban dau Xdot_MTLCD(:,1)=0; X2dot_MTLCD=zeros(dofs+2*N,length(accelerate)); % Vecto gia toc ban dau X2dot_MTLCD(:,1)=0; A=zeros(dofs+2*N,1); B=zeros(dofs+2*N,1); for i=1:length(accelerate)-1; for r=1:dofs+2*N A(r,1)=(-X_MTLCD(r,i)/beta/dT^2-Xdot_MTLCD(r,i)/beta/dT- (1/2/beta-1)*X2dot_MTLCD(r,i)); B(r,1)=(-gama*X_MTLCD(r,i)/beta/dT+(1-gama/beta)* Xdot_MTLCD(r,i)+(1-gama/2/beta)*dT*X2dot_MTLCD(r,i)); end % Chuyen vi u thu i K_eff=(1/beta/dT^2*M_MTLCD+gama/beta/dT*C_MTLCD+K_MTLCD); P_eff=(-ms_MTLCD*P(i+1)-M_MTLCD*A-C_MTLCD*B); X_MTLCD(:,i+1)=K_eff\P_eff; Xdot_MTLCD(:,i+1)=Xdot_MTLCD(:,i)*(1-gama/beta)+ 108 dT*X2dot_MTLCD(:,i)*(1-gama/2/beta)+ (X_MTLCD(:,i+1)-X_MTLCD(:,i))*gama/beta/dT; X2dot_MTLCD(:,i+1)=(X_MTLCD(:,i+1)-X_MTLCD(:,i))/(beta*dT^2)- Xdot_MTLCD(:,i)/(beta*dT)-(1/2/beta-1)* X2dot_MTLCD(:,i); end X17_MAX_MTLCD=max(X_MTLCD(nstory,:))*100 plot(stepT,X_MTLCD(dofs,:),'-b','linewidth',1.5); %% VE DO THI xlabel('Thoi gian (s)'); ylabel('Chuyen vi lon nhat tai tang 17 (m)'); % -legend('no-TLCD','1-TLCD','3-TLCD','5-TLCD'); hold off  Hệ kết cấu – MTLCD chịu tác động tải điều hòa clear all clc nstory = 17; % Bac tu cua ket cau chinh ms = ones(1,nstory)*500*10^3; % Khoi luong cua ket cau chinh M_all=0; for m=1:nstory; M_all=M_all+ms(1,m); end ks = ones(1,nstory)*7*10^8;% Do cung cua ket cau chinh dofs = length(ms); cs = zeros(dofs,1); % Luc can cua ket cau chinh khi1 = 0.01; % Ty so can cua ket cau chinh %% Ket noi ma tran khoi luong, cung cua ket cau chinh [Ms,Ks] = f_connection_matrices(ms,cs,ks,dofs); %% Phan tich dang dao dong [Omega,Phi,Frequency] = f_ModalAnalysis(Ms,Ks); omega1=Omega(1); %rad/s T1=2*pi/Omega(1) % %% Determing the effective mass [Meff] = f_EffectiveMass(Ms,Ks); %% Rayleigh damper pRay = inv([1/Omega(1) Omega(1);1/Omega(2) Omega(2)])*2*khi1*ones(2,1); Cs = pRay(1)*Ms + pRay(2)*Ks; X_noTLCD=zeros(dofs,1); %% disp('type = 1: Average Acceleration Method') % disp('type = 2: Linear Acceleration method') method={'Average Acceleration Method';'Linear Acceleration method'}; % type=input('Nhap loai noi suy, type ='); type=1; % switch type case disp('Average Acceleration Method') gama=1/2; beta=1/4; case disp('Linear Acceleration Method') 109 gama=1/2; beta=1/6; end %% Thong so tai tran dong dat % Ten tran dong dat disp('STT cac tran dong dat') disp(' El Centro S00E') disp(' Northridge Earthquake') disp(' ParkField California N85E') disp(' PETROLIA CHAN 2') EQ=input('Nhap tran dong dat, so thu tu ='); % Nhap du lieu buoc thoi gian if EQ==1 stepT = xlsread('Earthquake Datas.xls','Sheet2','A2:A2690');%ElCentro disp('tran dong dat da chon:El Centro S00E') elseif EQ==2 stepT = xlsread('Earthquake Datas.xls','Sheet2','F2:F3002');%Northridge disp('tran dong dat da chon:Northridge Earthquake') elseif EQ==3 stepT = xlsread('Earthquake Datas.xls','Sheet2','O2:O2200');%ParkField disp('tran dong dat da chon:ParkField California N85E') elseif EQ==4 stepT = xlsread('Earthquake Datas.xls','Sheet2','T2:T3002');%PETROLIA disp('tran dong dat da chon:PETROLIA CHAN 2') else error('nhap lai so thu tu tran dong dat') end % Nhap du lieu gia toc nen if EQ==1 accelerate = xlsread('Earthquake Datas.xls','Sheet2','B2:B2690');%ElCentro elseif EQ==2 accelerate = xlsread('Earthquake Datas.xls','Sheet2','G2:G3002');%Northridge elseif EQ==3 accelerate = xlsread('Earthquake Datas.xls','Sheet2','P2:P2200');%ParkField elseif EQ==4 accelerate = xlsread('Earthquake Datas.xls','Sheet2','U2:U3002');%PETROLIA else error('nhap lai so thu tu tran dong dat') end % timestep=0.031; %results=zeros(3,8); for e=1:length(stepT)-1 stepT(e+1)=stepT(e)+timestep; end %% Phan tich he BTD tim So cua tai dT=stepT(2)-stepT(1); % buoc thoi gian P1=-1*accelerate; % Tai % Dieu kien ban dau U=zeros(length(accelerate),1); % Vecto chuyen vi ban dau U(1)=0; 110 Udot=zeros(length(accelerate),1); % Vecto van toc ban dau Udot(1)=0; U2dot=zeros(length(accelerate),1); % Vecto gia toc ban dau U2dot(1)=0; for i=1:length(accelerate)-1; K_eff1=1/(beta*dT^2)+gama*2*khi1*omega1/(beta*dT)+omega1^2; P_eff1=P1(i+1)+1*(U(i)/(beta*dT^2)+Udot(i)/(beta*dT)+U2dot(i)* (1/2/beta-1))+ 2*khi1*omega1*(U(i)*gama/(beta*dT)+Udot(i)* (gama/beta-1)+dT/2*(gama/beta-2)*U2dot(i)); % Chuyen vi u thu i U(i+1)=P_eff1/K_eff1; Udot(i+1)=Udot(i)*(1-gama/beta)+dT*U2dot(i)*(1-gama/2/beta)+ (U(i+1)-U(i))*gama/beta/dT; U2dot(i+1)=(U(i+1)-U(i))/(beta*dT^2)-Udot(i)/(beta*dT)- (1/2/beta-1)*U2dot(i); end So=rms(U)^2*4*khi1*omega1^3/pi; %% Khong co he can (NoTLCD) -dT=stepT(2)-stepT(1); % Buoc thoi gian P=accelerate; % Tai tac dong % Dieu kien ban dau X_noTLCD=zeros(dofs,length(accelerate)); % Vecto chuyen vi ban dau X_noTLCD(:,1)=0; Xdot_noTLCD=zeros(dofs,length(accelerate)); % Vecto van toc ban dau Xdot_noTLCD(:,1)=0; X2dot_noTLCD=zeros(dofs,length(accelerate)); % Vecto gia toc ban dau X2dot_noTLCD(:,1)=0; A=zeros(dofs,1); B=zeros(dofs,1); for i=1:length(accelerate)-1; for k=1:dofs A(k,1)=(-X_noTLCD(k,i)/beta/dT^2-Xdot_noTLCD(k,i)/beta/dT- (1/2/beta-1)*X2dot_noTLCD(k,i)); B(k,1)=(-gama*X_noTLCD(k,i)/beta/dT+(1-gama/beta)* Xdot_noTLCD(k,i)+(1-gama/2/beta)*dT*X2dot_noTLCD(k,i)); end % Chuyen vi u thu i K_eff=(1/beta/dT^2*Ms+gama/beta/dT*Cs+Ks); P_eff=(-ms'*P(i+1)-Ms*A-Cs*B); X_noTLCD(:,i+1)=K_eff\P_eff; Xdot_noTLCD(:,i+1)=Xdot_noTLCD(:,i)*(1-gama/beta)+ dT*X2dot_noTLCD(:,i)*(1-gama/2/beta)+ (X_noTLCD(:,i+1)-X_noTLCD(:,i))*gama/beta/dT; X2dot_noTLCD(:,i+1)=(X_noTLCD(:,i+1)-X_noTLCD(:,i))/ (beta*dT^2)-Xdot_noTLCD(:,i)/(beta*dT)- (1/2/beta-1)*X2dot_noTLCD(:,i); end X17_MAX_noTLCD=max(X_noTLCD(nstory,:))*100 figure('color',[1 1]); plot(stepT,X_noTLCD(dofs,:),':k','linewidth',0.1); %% He co TLCD (1_TLCD) N=1; muy=0.01; % Ti so khoi luong (he can/ket cau chinh) gama1=1/(1+muy); % Ty so dieu chinh (tan so he can/tan so ket cau chinh) anpha=0.9; % Ti so giua B/L cua he can tau=1; % Ti so khoi luong be/khoi luong chat long omega2=omega1*gama1; % Tan so he can hold on 111 khi2=0.01; % Ti so can cua he can omegaL=omega1*gama1; % Tan so tu nhien cua chat long xi_L=100; [xichmaucham]=f_tim_xichmaucham_TDD(N,khi1,muy,anpha,tau,khi2,xi_L,omega1 ,omega2,omegaL,So); xichma_u_cham=xichmaucham; Cp=xichma_u_cham*xi_L/sqrt(2*pi); Cp_bd=Cp; % -M1=zeros(dofs+2); M1(dofs+1,dofs)=1; M1(dofs+1,dofs+1)=1; M1(dofs+1,dofs+2)=anpha/(1+tau); M1(dofs+2,dofs)=anpha; M1(dofs+2,dofs+1)=anpha; M1(dofs+2,dofs+2)=1; M2=[Ms zeros(dofs,2)]; M2=[M2;zeros(2,dofs+2)]; M_TLCD=M1+M2; % -C1=zeros(dofs+2); C1(dofs,dofs+1)=-2*khi2*omega2*muy*M_all; C1(dofs+1,dofs+1)=2*khi2*omega2; C1(dofs+2,dofs+2)=2*Cp; C2=[Cs zeros(dofs,2)]; C2=[C2;zeros(2,dofs+2)]; C_TLCD=C1+C2; % -K1=zeros(dofs+2); K1(dofs,dofs+1)=-omega2^2*muy*M_all; K1(dofs+1,dofs+1)=omega2^2; K1(dofs+2,dofs+2)=omegaL^2; K2=[Ks zeros(dofs,2)]; K2=[K2;zeros(2,dofs+2)]; K_TLCD=K1+K2; % -ms_TLCD=[ms';1;anpha]; % -% Initial conditions X_TLCD=zeros(dofs+2,length(accelerate)); % Vecto chuyen vi ban dau X_TLCD(:,1)=0; Xdot_TLCD=zeros(dofs+2,length(accelerate)); % Vecto van toc ban dau Xdot_TLCD(:,1)=0; X2dot_TLCD=zeros(dofs+2,length(accelerate)); % Vecto gia toc ban dau X2dot_TLCD(:,1)=0; A=zeros(dofs+2,1); B=zeros(dofs+2,1); for i=1:length(accelerate)-1; for r=1:dofs+2 A(r,1)=(-X_TLCD(r,i)/beta/dT^2-Xdot_TLCD(r,i)/beta/dT- (1/2/beta-1)*X2dot_TLCD(r,i)); B(r,1)=(-gama*X_TLCD(r,i)/beta/dT+(1-gama/beta)* Xdot_TLCD(r,i)+(1-gama/2/beta)*dT*X2dot_TLCD(r,i)); end % Chuyen vi u thu i K_eff=(1/beta/dT^2*M_TLCD+gama/beta/dT*C_TLCD+K_TLCD); P_eff=(-ms_TLCD*P(i+1)-M_TLCD*A-C_TLCD*B); 112 X_TLCD(:,i+1)=K_eff\P_eff; Xdot_TLCD(:,i+1)=Xdot_TLCD(:,i)*(1-gama/beta)+ dT*X2dot_TLCD(:,i)*(1-gama/2/beta)+ (X_TLCD(:,i+1)-X_TLCD(:,i))*gama/beta/dT; X2dot_TLCD(:,i+1)=(X_TLCD(:,i+1)-X_TLCD(:,i))/(beta*dT^2)- Xdot_TLCD(:,i)/(beta*dT)-(1/2/beta-1)* X2dot_TLCD(:,i); end X17_MAX_TLCD=max(X_TLCD(nstory,:))*100 plot(stepT,X_TLCD(dofs,:),'-g','linewidth',1.5); %% He co TLCD -N=3; % So luong he can muy=0.01; % Ti so khoi luong (he can/ket cau chinh) gama1=1/(1+muy); % Ty so dieu chinh (tan so he can/tan so ket cau chinh) anpha=0.9; % Ti so giua B/L cua he can tau=1; % Ti so khoi luong be/khoi luong chat long omega2=omega1*gama1; % Tan so he can hold on khi2=0.01; % Ti so can cua he can omegaL=omega1*gama1; % Tan so tu nhien cua chat long xi_L=100; [xichmaucham]=f_tim_xichmaucham_TDD(N,khi1,muy,anpha,tau,khi2,xi_L,omega1 ,omega2,omegaL,So); xichma_u_cham=xichmaucham; Cp=xichma_u_cham*xi_L/sqrt(2*pi); Cp_bd=Cp; % Ma tran khoi luong M1=zeros(dofs+2*N); for i=dofs+1:dofs+N M1(i,dofs)=muy*M_all; end for i=dofs+N+1:dofs+2*N M1(i,dofs)=anpha; end for i=dofs+1:dofs+N M1(i,i)=muy*M_all; end for i=dofs+N+1:dofs+2*N M1(i,i-N)=anpha; end for i=dofs+1:dofs+N M1(i,i+N)=muy*anpha*M_all/(1+tau); end for i=dofs+N+1:dofs+2*N M1(i,i)=1; end M2=[Ms zeros(dofs,2*N)]; M2=[M2;zeros(2*N,dofs+2*N)]; M_MTLCD=M1+M2; % Ma tran can -C1=zeros(dofs+2*N); for i=dofs+1:dofs+N C1(dofs,i)=-2*khi2*omega2*muy*M_all; end for i=dofs+1:dofs+N C1(i,i)=2*khi2*omega2*muy*M_all; end for i=dofs+N+1:dofs+2*N 113 C1(i,i)=2*Cp; end C2=[Cs zeros(dofs,2*N)]; C2=[C2;zeros(2*N,dofs+2*N)]; C_MTLCD=C1+C2; % Ma tran cung -K1=zeros(dofs+2*N); for i=dofs+1:dofs+N K1(dofs,i)=-omega2^2*muy*M_all; end for i=dofs+1:dofs+N K1(i,i)=omega2^2*muy*M_all; end for i=dofs+N+1:dofs+2*N K1(i,i)=omegaL^2; end K2=[Ks zeros(dofs,2*N)]; K2=[K2;zeros(2*N,dofs+2*N)]; K_MTLCD=K1+K2; % -ms_MTLCD_1=zeros(dofs+2*N,1); for i=dofs+1:dofs+N ms_MTLCD_1(i,1)=muy*M_all; end for i=dofs+N+1:dofs+2*N ms_MTLCD_1(i,1)=anpha; end ms_MTLCD_2=[ms';zeros(2*N,1)]; ms_MTLCD=ms_MTLCD_1+ms_MTLCD_2; % -% Dieu kien ban dau X_MTLCD=zeros(dofs+2*N,length(accelerate)); % Vecto chuyen vi ban dau X_MTLCD(:,1)=0; Xdot_MTLCD=zeros(dofs+2*N,length(accelerate)); % Vecto van toc ban dau Xdot_MTLCD(:,1)=0; X2dot_MTLCD=zeros(dofs+2*N,length(accelerate)); % Vecto gia toc ban dau X2dot_MTLCD(:,1)=0; A=zeros(dofs+2*N,1); B=zeros(dofs+2*N,1); for i=1:length(accelerate)-1; for r=1:dofs+2*N A(r,1)=(-X_MTLCD(r,i)/beta/dT^2-Xdot_MTLCD(r,i)/beta/dT- (1/2/beta-1)*X2dot_MTLCD(r,i)); B(r,1)=(-gama*X_MTLCD(r,i)/beta/dT+(1-gama/beta)* Xdot_MTLCD(r,i)+(1-gama/2/beta)*dT*X2dot_MTLCD(r,i)); end % Chuyen vi u thu i K_eff=(1/beta/dT^2*M_MTLCD+gama/beta/dT*C_MTLCD+K_MTLCD); P_eff=(-ms_MTLCD*P(i+1)-M_MTLCD*A-C_MTLCD*B); X_MTLCD(:,i+1)=K_eff\P_eff; Xdot_MTLCD(:,i+1)=Xdot_MTLCD(:,i)*(1-gama/beta)+ dT*X2dot_MTLCD(:,i)*(1-gama/2/beta)+ (X_MTLCD(:,i+1)-X_MTLCD(:,i))*gama/beta/dT; X2dot_MTLCD(:,i+1)=(X_MTLCD(:,i+1)-X_MTLCD(:,i))/(beta*dT^2)- Xdot_MTLCD(:,i)/(beta*dT)-(1/2/beta-1)* 114 X2dot_MTLCD(:,i); end X17_MAX_MTLCD=max(X_MTLCD(nstory,:))*100 plot(stepT,X_MTLCD(dofs,:),'-r','linewidth',1.5); %% He co TLCD -N=5; % So luong he can muy=0.01; % Ti so khoi luong (he can/ket cau chinh) gama1=1/(1+muy); % Ty so dieu chinh (tan so he can/tan so ket cau chinh) anpha=0.9; % Ti so giua B/L cua he can tau=1; % Ti so khoi luong be/khoi luong chat long omega2=omega1*gama1; % Tan so he can hold on khi2=0.01; % Ti so can cua he can omegaL=omega1*gama1; % Tan so tu nhien cua chat long xi_L=100; [xichmaucham]=f_tim_xichmaucham_TDD(N,khi1,muy,anpha,tau,khi2,xi_L,omega1 ,omega2,omegaL,So); xichma_u_cham=xichmaucham; Cp=xichma_u_cham*xi_L/sqrt(2*pi); Cp_bd=Cp; % Ma tran khoi luong M1=zeros(dofs+2*N); for i=dofs+1:dofs+N M1(i,dofs)=muy*M_all; end for i=dofs+N+1:dofs+2*N M1(i,dofs)=anpha; end for i=dofs+1:dofs+N M1(i,i)=muy*M_all; end for i=dofs+N+1:dofs+2*N M1(i,i-N)=anpha; end for i=dofs+1:dofs+N M1(i,i+N)=muy*anpha*M_all/(1+tau); end for i=dofs+N+1:dofs+2*N M1(i,i)=1; end M2=[Ms zeros(dofs,2*N)]; M2=[M2;zeros(2*N,dofs+2*N)]; M_MTLCD=M1+M2; % Ma tran can -C1=zeros(dofs+2*N); for i=dofs+1:dofs+N C1(dofs,i)=-2*khi2*omega2*muy*M_all; end for i=dofs+1:dofs+N C1(i,i)=2*khi2*omega2*muy*M_all; end for i=dofs+N+1:dofs+2*N C1(i,i)=2*Cp; end C2=[Cs zeros(dofs,2*N)]; C2=[C2;zeros(2*N,dofs+2*N)]; C_MTLCD=C1+C2; % Ma tran cung 115 K1=zeros(dofs+2*N); for i=dofs+1:dofs+N K1(dofs,i)=-omega2^2*muy*M_all; end for i=dofs+1:dofs+N K1(i,i)=omega2^2*muy*M_all; end for i=dofs+N+1:dofs+2*N K1(i,i)=omegaL^2; end K2=[Ks zeros(dofs,2*N)]; K2=[K2;zeros(2*N,dofs+2*N)]; K_MTLCD=K1+K2; % -ms_MTLCD_1=zeros(dofs+2*N,1); for i=dofs+1:dofs+N ms_MTLCD_1(i,1)=muy*M_all; end for i=dofs+N+1:dofs+2*N ms_MTLCD_1(i,1)=anpha; end ms_MTLCD_2=[ms';zeros(2*N,1)]; ms_MTLCD=ms_MTLCD_1+ms_MTLCD_2; % -% Dieu kien ban dau X_MTLCD=zeros(dofs+2*N,length(accelerate)); % Vecto chuyen vi ban dau X_MTLCD(:,1)=0; Xdot_MTLCD=zeros(dofs+2*N,length(accelerate)); % Vecto van toc ban dau Xdot_MTLCD(:,1)=0; X2dot_MTLCD=zeros(dofs+2*N,length(accelerate)); % Vecto gia toc ban dau X2dot_MTLCD(:,1)=0; A=zeros(dofs+2*N,1); B=zeros(dofs+2*N,1); for i=1:length(accelerate)-1; for r=1:dofs+2*N A(r,1)=(-X_MTLCD(r,i)/beta/dT^2-Xdot_MTLCD(r,i)/beta/dT- (1/2/beta-1)*X2dot_MTLCD(r,i)); B(r,1)=(-gama*X_MTLCD(r,i)/beta/dT+(1-gama/beta)* Xdot_MTLCD(r,i)+(1-gama/2/beta)*dT*X2dot_MTLCD(r,i)); end % Chuyen vi u thu i K_eff=(1/beta/dT^2*M_MTLCD+gama/beta/dT*C_MTLCD+K_MTLCD); P_eff=(-ms_MTLCD*P(i+1)-M_MTLCD*A-C_MTLCD*B); X_MTLCD(:,i+1)=K_eff\P_eff; Xdot_MTLCD(:,i+1)=Xdot_MTLCD(:,i)*(1-gama/beta)+ dT*X2dot_MTLCD(:,i)*(1-gama/2/beta)+ (X_MTLCD(:,i+1)-X_MTLCD(:,i))*gama/beta/dT; X2dot_MTLCD(:,i+1)=(X_MTLCD(:,i+1)-X_MTLCD(:,i))/(beta*dT^2)- Xdot_MTLCD(:,i)/(beta*dT)-(1/2/beta-1)* X2dot_MTLCD(:,i); end X17_MAX_MTLCD=max(X_MTLCD(nstory,:))*100 plot(stepT,X_MTLCD(dofs,:),'-b','linewidth',1.5); 116 LÝ LỊCH TRÍCH NGANG Họ tên: Trương Chế Quang Hiệu Ngày, tháng, năm sinh: 20/01/1994 Nơi sinh: Quảng Ngãi Địa liên lạc: 245 Nguyễn Trãi, Phương Nguyễn Cư Trinh, Quận 1, Tp.HCM ĐTDĐ: 0974 858 607 Email: truongchequanghieu@gmail.com QUÁ TRÌNH ĐÀO TẠO 2012 – 2016: Sinh viên đại học chun ngành Cảng – Cơng trình biển, Trường Đại Học Bách Khoa Tp Hồ Chí Minh 2017 - 2019: Học viên cao học chuyên ngành Xây dựng cơng trình dân dụng cơng nghiệp, Trường Đại học Bách Khoa Tp HCM ... cản cột chất lỏng Các hệ cản cột chất lỏng gắn tầng mái thông số xác định dựa đặc tính động lực học kết cấu hệ cản  Xây dựng mơ hình kết cấu có gắn hệ cản cột chất lỏng nhiều hệ cản cột chất lỏng: ... sử dụng nhiều cản điều chỉnh cột chất lỏng (MTLCD) Hình 2.5: Mơ hình mơ hệ cản nhiều cột chất lỏng [8] 2.3 Ứng dụng hệ cản cột chất lỏng hệ cản bể chất lỏng Trên giới ứng dụng nhiều hệ cản TLD,... giải, thuật tốn dùng để phân tích hiệu giảm dao động MTLCD kết cấu chịu tác động tải điều hòa tải động đất 3.2 Nguyên lý hoạt động hệ cản cột chất lỏng Hệ cản cột chất lỏng (TLCD) gồm có bể chứa