Phân tích hiệu quả giảm chấn của hệ cản MR nối giữa hai kết cấu

TRƯỜNG ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA - - LÊ THANH CƯỜNG PHÂN TÍCH HIỆU QUẢ GIẢM CHẤN CỦA HỆ CẢN MR NỐI GIỮA HAI KẾT CẤU Chuyên ngành : Xây Dựng Dân Dụng – Công Nghiệp Mã số ngành : 60 58 20 LUẬN VĂN THẠC SĨ Thành phố Hồ Chí Minh tháng 12 năm 2012 CƠNG TRÌNH ĐƯỢC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA – ĐHQG – TP.HCM Cán hướng dẫn khoa học: TS NGUYỄN TRỌNG PHƯỚC Cán chấm nhận xét 1: Cán chấm nhận xét 2: Luận văn thạc sĩ bảo vệ Trường Đại học Bách Khoa, ĐHQG, TP HCM Ngày tháng năm Thành phần đánh giá luận văn thạc sĩ gồm: Xác nhận Chủ tịch Hội đồng đánh giá LV Trưởng Khoa quản lý chuyên nghành sau luận văn sửa chữ (nếu có) CHỦ TỊCH HỘI ĐỒNG TRƯỞNG KHOA TRƯỜNG ĐẠI HỌC BÁCH KHOA KHOA XÂY DỰNG CỘNG HOÀ XÃ HỘI CHỦ NGHIÃ VIỆT NAM Độc Lập - Tự Do - Hạnh Phúc -o0o Tp HCM, ngày 02 tháng 07 năm 2012 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: LÊ THANH CƯỜNG Phái: NAM Ngày, tháng, năm sinh: 11/05/1981 Nơi sinh: HÀ TĨNH Chuyên ngành: Xây Dựng Dân Dụng – Cơng Nghiệp MSHV: 11211006 Khóa: 2011 1- TÊN ĐỀ TÀI: PHÂN TÍCH HIỆU QUẢ GIẢM CHẤN CỦA HỆ CẢN MR KHI NỐI GIỮA HAI KẾT CẤU 2- NHIỆM VỤ LUẬN VĂN: - Tìm hiểu hệ cản lưu biến từ MR: vật liệu, đặc tính mơ hình học - Áp dụng hệ cản MR để nối hai kết cấu phân tích hiệu tồn hệ chịu động đất: Thiết lập phương trình chủ đạo có mơ hình MR, giải tốn động lực học đánh giá kết - Viết chương trình máy tính ngơn ngữ MATLAB giải tốn động lực học trên: Trong bước thời gian mô ứng xử MR phương pháp Runge Kutta tồn miền thời gian phương pháp tích phân Newmark; lựa chọn liệu động đất phân tích phổ lượng theo chuỗi Fourier (FFT) - Đánh giá hiệu MR nối hai kết cấu 3- NGÀY GIAO NHIỆM VỤ : 02/07/2012 4- NGÀY HOÀN THÀNH NHIỆM VỤ : 25 - 12 - 2012 5- HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN: TS NGUYỄN TRỌNG PHƯỚC CÁN BỘ HƯỚNG DẪN (Họ tên chữ ký) TS Nguyễn Trọng Phước BAN QL CHUYÊN NGÀNH (Họ tên chữ ký) LỜI CẢM ƠN Trước tiên, xin gởi lời cảm ơn chân thành sâu sắc đến người hướng dẫn khoa học TS Nguyễn Trọng Phước Thầy tận tâm hướng dẫn, động viên tạo điều kiện cho học tập, nghiên cứu thực đề tài Trường Đại Học Bách Khoa TP.HCM Những dẫn Thầy không kiến thức khoa học quý báu giúp hồn thành luận văn mà Thầy cịn giúp tơi nhiều khả tư khoa học Tôi xin gởi lời cảm ơn chân thành đến Thầy, Cơ giảng dạy chương trình Sau đại học nghành Xây dựng dân dụng khóa 2011 truyền đạt cho kiến thức khoa học quý báu Tôi xin cảm ơn anh chị quản lý thư viện Trường Đại Học Học Bách Khoa TP.HCM tạo điều kiện thuận lợi cho tham khảo nguồn tài liệu quý giá trường Và cảm ơn đến bạn học viên ln có chia sẻ trao đổi kiến thức bổ ích suốt thời gian qua Cuối cùng, xin chân thành cảm ơn tất người thân, gia đình, thầy cô, bè bạn tri kỷ gắn bó tơi, khơng ngừng khuyến khích tơi suốt trình học tập, nghiên cứu khoa học thực đề tài Tôi xin chân thành cảm ơn ! TĨM TẮT “Phân tích hiệu giảm chấn hệ cản MR nối hai kết cấu” Lê Thanh Cường Luận văn phân tích hiệu giảm chấn hệ cản MR (MagnetoRheological, lưu biến từ) bố trí nối hai kết cấu chịu gia tốc động đất Thiết bị cản MR mô lị xo cản nhớt, lực cản sinh hàm phụ thuộc vào điện thơng số khác thiết bị Phương trình chuyển động kết cấu hệ cản chịu tác dụng gia tốc động đất thiết lập dựa nguyên lý cân động giải phương pháp tích phân trực tiếp bước theo thời gian Sự đáp ứng hệ cản MR bước thời gian mơ tả phương trình vi phân cấp giải phương pháp số Runge-Kutta Một chương trình máy tính viết ngơn ngữ MATLAB để phân tích phản ứng động hệ kết cấu Kết số từ phản ứng động gồm có chuyển vị, vận tốc lượng tiêu tán cho thấy hiệu hệ cản MR với điện áp điều khiển khác nối hai kết cấu ABSTRACT “Analyses effective for dynamic response reduction of MR damper – connected coubled structures” Le Thanh Cuong This thesis analysed dynamic reduction of MR damper (Magneto-Rheological) which set between two structures under ground motion MR damper is simulated by springs , viscous damper and dashpot This equipment produce a force which is a function depend on voltage and parameter The governing equations of motion of the coupled system is established base on principle of dynamic balancing, this equations are solved using Newmark’s step-by-step method The response of MR damper over a small time interval is simulated by 4th order Runge-Kutta method A program was establish base on Mathlap to analyse dynamic response of structures The result obtained from the numerical study is drawn effective of MR damper which is connected coupled structure MỤC LỤC Danh sách hình vẽ Danh sách bảng biểu Chữ viết tắt Chương GIỚI THIỆU 1.1 ĐẶT VẤN ĐỀ 1.2 MỤC TIÊU CỦA LUẬN VĂN 1.3 CẤU TRÚC CỦA LUẬN VĂN Chương TỔNG QUAN 2.1 GIỚI THIỆU 2.2 TỔNG QUAN ĐIỀU KHIỂN KẾT CẤU 2.2.1 Điều khiển bị động 2.2.2 Điều khiển chủ động 10 2.2.3 Điều khiển kết hợp 13 2.2.4 Điều khiển bán chủ động 13 2.3 CHẤT LƯU MR 16 2.3.1 Thành phần cấu tạo 16 2.3.2 Cơ chế hoạt động 17 2.4 CẤU TẠO HỆ CẢN MR 18 2.4.1 Cấu tạo chung 18 2.4.2 Giá thành số loại hệ cản MR 19 2.4.3 Mơ hình học hệ cản MR 19 2.5 LIÊN KẾT HAI CƠNG TRÌNH 25 2.6 KẾT LUẬN CHƯƠNG 28 Chương CƠ SỞ LÝ THUYẾT 29 3.1 GIỚI THIỆU 29 3.2 MƠ HÌNH KẾT CẤU 29 3.3 MƠ HÌNH CƠ HỌC 30 3.4 THIẾT LẬP PHƯƠNG TRÌNH CHUYỀN ĐỘNG 33 3.5 TÍNH TỐN LỰC ĐIỀU KHIỂN MR 38 3.6 CÁC PHƯƠNG PHÁP SỐ GIẢI HỆ PHƯƠNG TRÌNH VI PHÂN ĐỘNG LỰC HỌC 39 3.6.1 Phương pháp tích phân Newmark 40 3.6.2 Thuật toán giải phương trình vi phân chuyển động 41 3.7 NĂNG LƯỢNG 44 3.8 KẾT LUẬN CHƯƠNG 45 Chương VÍ DỤ SỐ 45 4.1 GIỚI THIỆU 45 4.2 MÔ TẢ KẾT CẤU GHÉP NỐI 45 4.3 SỐ LIỆU TRẬN ĐỘNG ĐẤT SỬ DỤNG ĐỂ PHÂN TÍCH KẾT CẤU 46 4.4 PHÂN TÍCH ĐÁP ỨNG CỦA KẾT CẤU DƯỚI TẢI ĐỘNG ĐẤT ELCENTRO 51 4.4.1 Phân tích hiệu giảm đáp ứng hệ cản MR 51 4.4.2 Phân tích ảnh hưởng số lượng hệ cản MR 62 4.4.3 Nhận xét 64 4.5 PHÂN TÍCH ĐÁP ỨNG CỦA KẾT CẤU DƯỚI TẢI ĐỘNG ĐẤT SUPERSTITION 65 4.5.1 Phân tích hiệu giảm đáp ứng hệ cản MR 65 4.5.2 Phân tích ảnh hưởng số lượng hệ cản MR 76 4.5.3 Nhận xét 77 4.6 PHÂN TÍCH ĐÁP ỨNG CỦA KẾT CẤU DƯỚI TẢI ĐỘNG ĐẤT NORTHRIDGE 79 4.6.1 Phân tích hiệu giảm đáp ứng hệ cản MR 79 4.6.2 Phân tích ảnh hưởng số lượng hệ cản MR 84 4.6.3 Nhận xét 85 4.7 KẾT LUẬN CHƯƠNG 86 Chương KẾT LUẬN 88 5.1 KẾT LUẬN 88 5.2 HƯỚNG PHÁT TRIỂN ĐỀ TÀI 89 TÀI LIỆU THAM KHẢO 90 PHỤ LỤC 94 DANH MỤC CÁC HÌNH VẼ Hình 1.1 Động đất Northridge, California – 1994 Hình 1.2 Động đất Kobe, Nhật Bản – 1995 Hình 1.3 Động đất Tứ Xuyên, Trung Quốc – 2008 Hình 1.4 Động đất Haiti – 2010 Hình 1.5 Động đất Đơng Bắc, Nhật Bản – 2012 Hình 2.1Phân loại điều khiển kết cấu Hình 2.2 Hệ điều khiển bị động – Liquid Tuned Mass Damper Hình 2.3 Hệ điều khiển bị động – Viscous Fluid Dampers 10 Hình 2.4 Hệ điều khiển bị động –Base Isolation 10 Hình 2.5 Điều khiển chủ động –Tịa nhàKyobashi Seiwa 12 Hình 2.6Điều khiển chủ động –Tòa nhàApplause 12 Hình 2.7Cơ cấu điều khiển kết hợp khối lượng 13 Hình 2.8Viện bảo tàng quốc gia Tokyo thiết bị MR-30T 15 Hình 2.9Cầu Dongting Lake thiết bị MR Lord SD-1005 15 Hình 2.10Cầu Binzhou Yellow River 16 103 bar(Y,1); colormap hot legend('Passive-off','Passive-on'); xlabel('Tang')%,ylabel('Do giam (%)'); title('BUILDING 2'); %% CHUYEN VI TRUNG BINH figure(4) subplot(1,2,1); hold on box on axis([0 16.5 -5 5]) tam1=mean_amt(u_building1_sg.*100); tam2=mean_amt(u_building1_db_Vmin); tam3=mean_amt(u_building1_db_Vmax); Y1=[tam1(:,1) tam2(:,1) tam3(:,1)]; Y2=[tam1(:,2) tam2(:,2) tam3(:,2)]; bar(Y1,1); bar(Y2,1); legend('uncontrol','Passive-off','Passive-on'); xlabel('Tang'),ylabel('Chuyen Vi Trung Binh(cm)'); title('BUILDING 1'); subplot(1,2,2); hold on box on grid off axis([0 8.5 -1 1]) tam1=mean_amt(u_building2_sg.*100); tam2=mean_amt(u_building2_db_Vmin); tam3=mean_amt(u_building2_db_Vmax); Y1=[tam1(:,1) tam2(:,1) tam3(:,1)]; Y2=[tam1(:,2) tam2(:,2) tam3(:,2)]; bar(Y1,1); bar(Y2,1); legend('uncontrol','Passive-off','Passive-on'); xlabel('Tang'), ylabel('Chuyen Vi Trung Binh(cm)'); title('BUILDING 2'); %% DO GIAM CHUYEN VI TRUNG BINH figure(5) subplot(1,2,1); hold on box on axis([0 17.5 70]) Pa1=max((abs(mean_amt(u_building1_sg.*100)))') ; Pa2=max((abs(mean_amt(u_building1_db_Vmin)))') ; Pa3=max((abs(mean_amt(u_building1_db_Vmax)))') ; 104 Reduce1=(Pa1'-Pa2').*100./Pa1' ; Reduce2=(Pa1'Pa3').*100./Pa1'; Y=[Reduce1 Reduce2]; bar(Y,1.5); colormap hot legend('Passive-off','Passive-on') ; xlabel('Tang'),ylabel('Do giam (%)'); title('BUILDING 1'); subplot(1,2,2); hold on box on axis([0 8.5 70]) Pa1=max((abs(mean_amt(u_building2_sg.*100)))') ; Pa2=max((abs(mean_amt(u_building2_db_Vmin)))'); Pa3=max((abs(mean_amt(u_building2_db_Vmax)))') ; Reduce1=(Pa1'-Pa2').*100./Pa1' ; Reduce2=(Pa1'Pa3').*100./Pa1'; Y=[Reduce1 Reduce2]; bar(Y,1.5); colormap hot legend('Passive-off','Passive-on'); xlabel('Tang'),ylabel('Do giam (%)'); title('BUILDING 2'); %% CHUYEN VI BINH PHUONG TRUNG BINH figure(6) subplot(1,2,1); hold on box on tam1=rms(u_building1_sg.*100); tam2=rms(u_building1_db_Vmin); tam3=rms(u_building1_db_Vmax); axis([0 16.5 6]) Y=[tam1 tam2 tam3]; bar(Y,1); legend('uncontrol','Passive-off','Passive-on'); xlabel('Tang'),ylabel('Chuyen Vi (cm)'); title('BUILDING 1'); subplot(1,2,2); hold on box on tam1=rms(u_building2_sg.*100); tam2=rms(u_building2_db_Vmin); tam3=rms(u_building2_db_Vmax); axis([0 8.5 1.5]) Y=[tam1 tam2 tam3]; bar(Y,1); xlabel('Tang') 105 title('BUILDING 2'); %% DO GIAM CHUYEN VI BINH PHUONG TRUNG BINH figure(7) subplot(1,2,1); hold on box on tam1=rms(u_building1_sg.*100); tam2=rms(u_building1_db_Vmin); tam3=rms(u_building1_db_Vmax); Reduce1=(tam1-tam2).*100./tam1 ; Reduce2=(tam1tam3).*100./tam1; Y=[Reduce1 Reduce2]; bar(Y,1); colormap hot legend('Passive-off','Passive-on'); xlabel('Tang'),ylabel('Do giam (%)'); title('BUILDING 1'); subplot(1,2,2); hold on box on tam1=rms(u_building2_sg.*100); tam2=rms(u_building2_db_Vmin); tam3=rms(u_building2_db_Vmax); Reduce1=(tam1-tam2).*100./tam1 ; Reduce2=(tam1tam3).*100./tam1; Y=[Reduce1 Reduce2]; bar(Y,1); colormap hot legend('Passive-off','Passive-on'); xlabel('Tang'),ylabel('Do giam (%)'); title('BUILDING 2'); Giatoc.m clc clear all %% Passive-on load Vmax_all udot2_building1_db_Vmax10=udot2_building1_db(16,:); udot2_building2_db_Vmax10=udot2_building2_db(8,:); udot2_building1_db_Vmax=udot2_building1_db; udot2_building2_db_Vmax=udot2_building2_db; %% Passive-off 106 ks_1ks_2u_doubleu_building1_sgu_building2_sg udot2_building1_db_Vmin10=udot2_building1_db(16,:); udot2_building2_db_Vmin10=udot2_building2_db(8,:); udot2_building1_db_Vmin=udot2_building1_db; udot2_building2_db_Vmin=udot2_building2_db; %% GIA TOC TANG DINH figure(1) subplot(1,2,1); axis([0 25 -6 6]) hold on grid off box on u_sg=udot2_building1_sg(16,:); plot(t,udot2_building1_sg(16,:),':k','LineWidth',1.5); plot(t,udot2_building1_db_Vmin10,' b','LineWidth',1.5); plot(t,udot2_building1_db_Vmax10,'-r','LineWidth',1.5); legend('uncontrol','Passive-off','Passive-on'); xlabel('Thoi gian(s)'),ylabel('Gia Toc (m/s2)'); title('BUILDING 1'); u_inter=abs(u_sg); index=find(u_inter==max(u_inter)); mainPeriodStr=num2str(u_sg(index)); plot(t(index),u_sg(index),'r.', 'MarkerSize',25); text(t(index)+1,u_sg(index),mainPeriodStr); subplot(1,2,2); axis([0 25 -8 8]) hold on grid off box on u_sg=udot2_building2_sg(8,:); plot(t,udot2_building2_sg(8,:),':k','LineWidth',1.5); plot(t,udot2_building2_db_Vmin10,' b','LineWidth',1.5); plot(t,udot2_building2_db_Vmax10,'-r','LineWidth',1.5); legend('uncontrol','Passive-off','Passive-on'); xlabel('Thoi gian(s)')%,ylabel('Gia Toc (m/s2)'); title('BUILDING 2'); u_inter=abs(u_sg); index=find(u_inter==max(u_inter)); mainPeriodStr=num2str(u_sg(index)); plot(t(index),u_sg(index),'r.', 'MarkerSize',25); text(t(index)+1,u_sg(index),mainPeriodStr); %% GIA TOC MAX figure(2) subplot(1,2,1); hold on box on 107 grid off axis([1 16.5 6]) plot(1:n_1,max(abs(udot2_building1_sg')),'s','LineWidth',1.5,'Color','k'); plot(1:n_1,max(abs(udot2_building1_db_Vmin')),'o','LineWidth',1.5,'Color','b'); plot(1:n_1,max(abs(udot2_building1_db_Vmax')),'s','LineWidth',1.5,'Color','r'); legend('uncontrol','Passive-off','Passive-on'); xlabel('Tang'),ylabel('Gia Toc (m/s2)'); title('BUILDING 1'); subplot(1,2,2); hold on box on grid off axis([1 8.5 8]) plot(1:n_2,max(abs(udot2_building2_sg')),'s','LineWidth',1.5,'Color','k'); plot(1:n_2,max(abs(udot2_building2_db_Vmin')),'o','LineWidth',1.5,'Color','b'); plot(1:n_2,max(abs(udot2_building2_db_Vmax')),'p','LineWidth',1.5,'Color','r'); legend('uncontrol','Passive-off','Passive-on'); xlabel('Tang')%,ylabel('Gia Toc (m/s2)'); title('BUILDING 2'); %% DO GIAM GIA TOC MAX figure(3) subplot(1,2,1); hold on box on axis([0 16.5 -1 15]) Pa1=max(abs(udot2_building1_sg'))' ; Pa2=max(abs(udot2_building1_db_Vmin'))'; Pa3=max(abs(udot2_building1_db_Vmax'))' ; Reduce1=(Pa1-Pa2).*100./Pa1 ; Reduce2=(Pa1Pa3).*100./Pa1; Y=[Reduce1 Reduce2]; bar(Y,1); legend('Passive-off','Passive-on'); xlabel('Tang'),ylabel('Do giam (%)'); title('BUILDING 1'); subplot(1,2,2); hold on box on axis([0 8.5 -1 15]) Pa1=max(abs(udot2_building2_sg'))' ; Pa2=max(abs(udot2_building2_db_Vmin'))'; 108 Pa3=max(abs(udot2_building2_db_Vmax'))' ; Reduce1=(Pa1-Pa2).*100./Pa1 ; Reduce2=(Pa1Pa3).*100./Pa1; Y=[Reduce1 Reduce2]; bar(Y,1); colormap hot legend('Passive-off','Passive-on'); xlabel('Tang')%,ylabel('Do giam (%)'); title('BUILDING 2'); %% GIA TOC TRUNG BINH figure(4) subplot(1,2,1); hold on box on axis([0 16.5 -1 1]) tam1=mean_amt(udot2_building1_sg); tam2=mean_amt(udot2_building1_db_Vmin); tam3=mean_amt(udot2_building1_db_Vmax); Y1=[tam1(:,1) tam2(:,1) tam3(:,1)]; Y2=[tam1(:,2) tam2(:,2) tam3(:,2)]; bar(Y1,1); bar(Y2,1); legend('uncontrol','Passive-off','Passive-on'); xlabel('Tang'),ylabel('Gia Toc Trung Binh(m/s2)'); title('BUILDING 1'); subplot(1,2,2); hold on box on grid off axis([0 8.5 -1 1]) tam1=mean_amt(udot2_building2_sg); tam2=mean_amt(udot2_building2_db_Vmin); tam3=mean_amt(udot2_building2_db_Vmax); Y1=[tam1(:,1) tam2(:,1) tam3(:,1)]; Y2=[tam1(:,2) tam2(:,2) tam3(:,2)]; bar(Y1,1); bar(Y2,1); legend('uncontrol','Passive-off','Passive-on'); xlabel('Tang'), ylabel('Gia Toc Trung Binh(m/s2)'); title('BUILDING 2'); %% DO GIAM GIA TOC TRUNG BINH figure(5) subplot(1,2,1); hold on box on axis([0 17.5 50]) 109 Pa1=max((abs(mean_amt(udot2_building1_sg)))') ; Pa2=max((abs(mean_amt(udot2_building1_db_Vmin)))'); Pa3=max((abs(mean_amt(udot2_building1_db_Vmax)))') ; Reduce1=(Pa1'-Pa2').*100./Pa1' ; Reduce2=(Pa1'Pa3').*100./Pa1'; Y=[Reduce1 Reduce2]; bar(Y,1); colormap hot legend('Passive-off','Passive-on'); xlabel('Tang'),ylabel('Do giam (%)'); title('BUILDING 1'); subplot(1,2,2); hold on box on axis([0 8.5 50]) Pa1=max((abs(mean_amt(udot2_building2_sg)))') ; Pa2=max((abs(mean_amt(udot2_building2_db_Vmin)))'); Pa3=max((abs(mean_amt(udot2_building2_db_Vmax)))') ; Reduce1=(Pa1'-Pa2').*100./Pa1' ; Reduce2=(Pa1'Pa3').*100./Pa1'; Y=[Reduce1 Reduce2]; bar(Y,1); colormap hot legend('Passive-off','Passive-on'); xlabel('Tang'),ylabel('Do giam (%)'); title('BUILDING 2'); %% CHUYEN VI BINH PHUONG TRUNG BINH figure(6) subplot(1,2,1); hold on box on tam1=rms(udot2_building1_sg); tam2=rms(udot2_building1_db_Vmin); tam3=rms(udot2_building1_db_Vmax); %axis([0 16.5 6]) Y=[tam1 tam2 tam3]; bar(Y,1); legend('uncontrol','Passive-off','Passive-on'); xlabel('Tang'),ylabel('Gia Toc (m/s2)'); title('BUILDING 1'); subplot(1,2,2); hold on box on tam1=rms(udot2_building2_sg); tam2=rms(udot2_building2_db_Vmin); tam3=rms(udot2_building2_db_Vmax); %axis([0 8.5 1.5]) 110 Y=[tam1 tam2 tam3]; bar(Y,1); xlabel('Tang') title('BUILDING 2'); %% DO GIAM GIA TOC TRUNG BINH figure(7) subplot(1,2,1); hold on box on %axis([0 17.5 50]) tam1=rms(udot2_building1_sg); tam2=rms(udot2_building1_db_Vmin); tam3=rms(udot2_building1_db_Vmax); Reduce1=(tam1-tam2).*100./tam1 tam3).*100./tam1; Y=[Reduce1 Reduce2]; bar(Y,1); colormap hot legend('Passive-off','Passive-on'); xlabel('Tang'),ylabel('Do giam (%)'); title('BUILDING 1'); subplot(1,2,2); hold on box on tam1=rms(udot2_building2_sg); tam2=rms(udot2_building2_db_Vmin); tam3=rms(udot2_building2_db_Vmax); Reduce1=(tam1-tam2).*100./tam1 tam3).*100./tam1; Y=[Reduce1 Reduce2]; bar(Y,1); colormap hot legend('Passive-off','Passive-on'); xlabel('Tang'),ylabel('Do giam (%)'); title('BUILDING 2'); Luccat.m clc clear all %% Passive-on load Vmax_all % Tinh luc cat cong trinh Luccat1_db_Vmax=zeros(n_1,length(t)); Luccat1_sg=zeros(n_1,length(t)); ; Reduce2=(tam1- ; Reduce2=(tam1- 111 Luccat1_sg(1,:)=u_building1_sg(1,:).*ks_1(1); Luccat1_db_Vmax(1,:)=u_building1_db(1,:).*ks_1(1); for i=2:n_1 Luccat1_db_Vmax(i,:)=ks_1(i).*(u_building1_db(i,:)u_building1_db(i-1,:)); Luccat1_sg(i,:)=ks_1(i).*(u_building1_sg(i,:)u_building1_sg(i-1,:)); end Luccat1_sg=Luccat1_sg./1000; % Chuyen ve KN Luccat1_db_Vmax=Luccat1_db_Vmax./1000; % Chuyen ve KN % Tinh luc cat cong trinh Luccat2_db_Vmax=zeros(n_2,length(t)); Luccat2_sg=zeros(n_2,length(t)); Luccat2_sg(1,:)=u_building2_sg(1,:).*ks_2(1); Luccat2_db_Vmax(1,:)=u_building2_db(1,:).*ks_2(1); for i=2:n_2 Luccat2_db_Vmax(i,:)=ks_2(i).*(u_building2_db(i,:)u_building2_db(i-1,:)); Luccat2_sg(i,:)=ks_2(i).*(u_building2_sg(i,:)u_building2_sg(i-1,:)); end Luccat2_sg=Luccat2_sg./1000; % Chuyen ve KN Luccat2_db_Vmax=Luccat2_db_Vmax./1000; % Chuyen ve KN %% Passive-off load Vmin_all % Tinh luc cat cong trinh Luccat1_db_Vmin=zeros(n_1,length(t)); Luccat1_db_Vmin(1,:)=u_building1_db(1,:).*ks_1(1); for i=2:n_1 Luccat1_db_Vmin(i,:)=ks_1(i).*(u_building1_db(i,:)u_building1_db(i-1,:)); end Luccat1_db_Vmin=Luccat1_db_Vmin./1000; % Chuyen ve KN % Tinh luc cat cong trinh Luccat2_db_Vmin=zeros(n_2,length(t)); Luccat2_db_Vmin(1,:)=u_building2_db(1,:).*ks_2(1); for i=2:n_2 Luccat2_db_Vmin(i,:)=ks_2(i).*(u_building2_db(i,:)u_building2_db(i-1,:)); end Luccat2_db_Vmin=Luccat2_db_Vmin./1000; % Chuyen ve KN %LUC CAT MAX figure(1) subplot(1,2,1); hold on 112 box on %axis([0 25 -6 6]) stairs(1:n_1,max(abs(Luccat1_sg')),'-k','LineWidth',2); stairs(1:n_1,max(abs(Luccat1_db_Vmin')),'-b','LineWidth',2.5); stairs(1:n_1,max(abs(Luccat1_db_Vmax')),':r','LineWidth',2.5) ; legend('uncontrol','Passive-off','Passive-on'); xlabel('Tang'),ylabel('Luc cat (KN)'); title('BUILDING 1'); subplot(1,2,2); hold on box on stairs(1:n_2,max(abs(Luccat2_sg')),'-k','LineWidth',2); stairs(1:n_2,max(abs(Luccat2_db_Vmin')),'-b','LineWidth',2.5); stairs(1:n_2,max(abs(Luccat2_db_Vmax')),':r','LineWidth',2.5) ; legend('uncontrol','Passive-off','Passive-on'); xlabel('Tang'),ylabel('Luc cat (KN)'); title('BUILDING 2'); %% DO GIAM LUC CAT MAX figure(2) subplot(1,2,1); hold on box on axis([0 16.5 40]) Pa1=max(abs(Luccat1_sg'))'; Pa2=max(abs(Luccat1_db_Vmin'))'; Pa3=max(abs(Luccat1_db_Vmax'))'; Reduce1=(Pa1-Pa2).*100./Pa1; Reduce2=(Pa1-Pa3).*100./Pa1; Y=[Reduce1 Reduce2]; bar(Y,1); legend('Passive-off','Passive-on'); xlabel('Tang'),ylabel('Do giam (%)'); title('BUILDING 1'); subplot(1,2,2); hold on box on axis([0 8.5 40]) Pa1=max(abs(Luccat2_sg'))'; Pa2=max(abs(Luccat2_db_Vmin'))'; Pa3=max(abs(Luccat2_db_Vmax'))'; Reduce1=(Pa1-Pa2).*100./Pa1; Reduce2=(Pa1-Pa3).*100./Pa1; Y=[Reduce1 Reduce2]; bar(Y,1); legend('Passive-off','Passive-on'); 113 xlabel('Tang'),ylabel('Do giam (%)'); title('BUILDING 2'); colormap hot %% LUC CAT TANG figure(3) subplot(1,2,1); hold on box on u_sg=Luccat1_sg(1,:); grid off axis([0 25 -4000 4000]) plot(t,Luccat1_sg(1,:),':k','LineWidth',1.5); plot(t,Luccat1_db_Vmin(1,:),' b','LineWidth',1.5); plot(t,Luccat1_db_Vmax(1,:),'-r','LineWidth',1.5); legend('uncontrol','Passive-off','Passive-on'); xlabel('Thoi gian(s)'),ylabel('Luc cat (KN)'); title('BUILDING 1'); u_inter=abs(u_sg); index=find(u_inter==max(u_inter)); mainPeriodStr=num2str(u_sg(index)); plot(t(index),u_sg(index),'r.', 'MarkerSize',25); text(t(index)+1,u_sg(index),mainPeriodStr); subplot(1,2,2); hold on box on u_sg=Luccat2_sg(1,:); grid off axis([0 25 -3000 3000]) plot(t,Luccat2_sg(1,:),':k','LineWidth',1.5); plot(t,Luccat2_db_Vmin(1,:),' b','LineWidth',1.5); plot(t,Luccat2_db_Vmax(1,:),'-r','LineWidth',1.5); legend('uncontrol','Passive-off','Passive-on'); xlabel('Thoi gian(s)'),ylabel('Luc cat (KN)'); title('BUILDING 2'); u_inter=abs(u_sg); index=find(u_inter==max(u_inter)); mainPeriodStr=num2str(u_sg(index)); plot(t(index),u_sg(index),'r.', 'MarkerSize',25); text(t(index)+1,u_sg(index),mainPeriodStr); %% LUC CAT TRUNG BINH figure(4) subplot(1,2,1); hold on box on axis([0 16.5 -1500 1500]) tam1=mean_amt(Luccat1_sg); 114 tam2=mean_amt(Luccat1_db_Vmin); tam3=mean_amt(Luccat1_db_Vmax); Y1=[tam1(:,1) tam2(:,1) tam3(:,1)]; Y2=[tam1(:,2) tam2(:,2) tam3(:,2)]; bar(Y1,1); bar(Y2,1); legend('uncontrol','Passive-off','Passive-on'); xlabel('Tang'),ylabel('Luc Cat Trung Binh(KN)'); title('BUILDING 1'); subplot(1,2,2); hold on box on grid off axis([0 8.5 -600 600]) tam1=mean_amt(Luccat2_sg); tam2=mean_amt(Luccat2_db_Vmin); tam3=mean_amt(Luccat2_db_Vmax); Y1=[tam1(:,1) tam2(:,1) tam3(:,1)]; Y2=[tam1(:,2) tam2(:,2) tam3(:,2)]; bar(Y1,1); bar(Y2,1); legend('uncontrol','Passive-off','Passive-on'); xlabel('Tang'), ylabel('Luc Cat Trung Binh(KN)'); title('BUILDING 2'); %% DO GIAM LUC CAT TRUNG BINH figure(5) subplot(1,2,1); hold on box on axis([0 17.5 70]) Pa1=max((abs(mean_amt(Luccat1_sg)))') ; Pa2=max((abs(mean_amt(Luccat1_db_Vmin)))'); Pa3=max((abs(mean_amt(Luccat1_db_Vmax)))') ; Reduce1=(Pa1'-Pa2').*100./Pa1' ; Reduce2=(Pa1'Pa3').*100./Pa1'; Y=[Reduce1 Reduce2]; bar(Y,1.5); colormap hot legend('Passive-off','Passive-on'); xlabel('Tang'),ylabel('Do giam (%)'); title('BUILDING 1'); subplot(1,2,2); hold on box on axis([0 8.5 70]) Pa1=max((abs(mean_amt(Luccat2_sg)))') ; Pa2=max((abs(mean_amt(Luccat2_db_Vmin)))'); Pa3=max((abs(mean_amt(Luccat2_db_Vmax)))') ; 115 Reduce1=(Pa1'-Pa2').*100./Pa1' Pa3').*100./Pa1'; Y=[Reduce1 Reduce2]; bar(Y,1.5); colormap hot legend('Passive-off','Passive-on'); xlabel('Tang'),ylabel('Do giam (%)'); title('BUILDING 2'); ; Reduce2=(Pa1'- %% LUC CAT BINH PHUONG TRUNG BINH figure(6) subplot(1,2,1); hold on box on tam1=rms(Luccat1_sg); tam2=rms(Luccat1_db_Vmin); tam3=rms(Luccat1_db_Vmax); Y=[tam1 tam2 tam3]; bar(Y,1); legend('uncontrol','Passive-off','Passive-on'); xlabel('Tang'),ylabel('Luc Cat (KN)'); title('BUILDING 1'); subplot(1,2,2); hold on box on tam1=rms(Luccat2_sg); tam2=rms(Luccat2_db_Vmin); tam3=rms(Luccat2_db_Vmax); Y=[tam1 tam2 tam3]; bar(Y,1); xlabel('Tang') title('BUILDING 2'); %% DO GIAM LUC CAT BINH PHUONG TRUNG BINH figure(7) subplot(1,2,1); hold on box on tam1=rms(Luccat1_sg); tam2=rms(Luccat1_db_Vmin); tam3=rms(Luccat1_db_Vmax); Reduce1=(tam1-tam2).*100./tam1 ; Reduce2=(tam1tam3).*100./tam1; Y=[Reduce1 Reduce2]; bar(Y,1); colormap hot legend('Passive-off','Passive-on'); xlabel('Tang'),ylabel('Do giam (%)'); 116 title('BUILDING 1'); subplot(1,2,2); hold on box on tam1=rms(Luccat2_sg); tam2=rms(Luccat2_db_Vmin); tam3=rms(Luccat2_db_Vmax); Reduce1=(tam1-tam2).*100./tam1 tam3).*100./tam1; Y=[Reduce1 Reduce2]; bar(Y,1); colormap hot legend('Passive-off','Passive-on'); xlabel('Tang'),ylabel('Do giam (%)'); title('BUILDING 2'); ; Reduce2=(tam1- mean_amt.m function [mean_amt_x]=mean_amt(x) [row_x,column_x]=size(x); for j=1:row_x tam_pos=0; count_pos=0; tam_neg=0; count_neg=0; for k=1:column_x if x(j,k)>0 count_pos=count_pos+1; tam_pos=tam_pos+x(j,k); elseif x(j,k)