BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƢỜNG ĐẠI HỌC CƠNG NGHỆ TP.HCM  NGUYỄN VĂN THI ẢNH HƢỞNG CỦA BỀ RỘNG TIẾT DIỆN DẦM GIẢM YẾU ĐẾN PHẢN ỨNG VỚI ĐỘNG ĐẤT CỦA MỘT KHUNG THÉP ĐẶC BIỆT BA TẦNG LUẬN VĂN THẠC SĨ Chun ngành: KTXD Cơng trình Dân dụng Cơng nghiệp Mã ngành: 60 58 02 08 TP Hồ Chí Minh, tháng 09 năm 2015 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƢỜNG ĐẠI HỌC CƠNG NGHỆ TP.HCM  NGUYỄN VĂN THI ẢNH HƢỞNGCỦA BỀ RỘNG TIẾT DIỆN DẦM GIẢM YẾU ĐẾN PHẢN ỨNG VỚI ĐỘNG ĐẤT CỦA MỘT KHUNG THÉP ĐẶC BIỆT BA TẦNG LUẬN VĂN THẠC SĨ Chun ngành: KTXD Cơng trình Dân dụng Cơng nghiệp Mã ngành: 60 58 02 08 CÁN BỘ HDKH: TS ĐÀO ĐÌNH NHÂN TP Hồ Chí Minh, tháng 09 năm 2015 i CƠNG TRÌNH ĐƢỢC HỒN THÀNH TẠI TRƢỜNG ĐẠI HỌC CƠNG NGHỆ TP.HCM Cán hƣớng dẫn khoa học: TS ĐÀO ĐÌNH NHÂN Luận văn Thạc sĩ đƣợc bảo vệ Trƣờng Đại học Kỹ thuật Cơng nghệ TP.HCM ngày … tháng … năm 2015 Thành phần Hội đồng đánh giá Luận văn Thạc sĩ gồm: Họ tên TT Chức danh Hội đồng Chủ tịch PGS.TS Võ Phán PGS.TS Lƣơng Văn Hải Phản biện TS Nguyễn Hồng Ân Phản biện PGS.TS Nguyễn Xn Hùng Ủy viên PGS.TS Dƣơng Hồng Thẩm Ủy viên, Thƣ ký Xác nhận Chủ tịch Hội đồng đánh giá Luận sau Luận văn đƣợc sửa chữa (nếu có) Chủ tịch Hội đồng đánh giá Luận văn ii BỘ GIÁO DỤC VÀ ĐÀO TẠO CỘNG HÕA XÃ HỘI CHỦ NGHĨA VIỆT NAM TRƢỜNG ĐH CƠNG NGHỆ TP HCM Độ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: NGUYỄN VĂN THI Giới tính: Nam Ngày, tháng, năm sinh: 08/10/1989 Nơi sinh: Bình Định Chun ngành: Kỹ thuật xây dựng MSHV: 1341870050 cơng trình dân dụng cơng nghiệp I Tên đề tài Ảnh hƣởngcủa bề rộng tiết diện dầm giảm yếu đến phản ứng với động đất khung thép đặc biệt ba tầng II Nhiệm vụ nội dung Nhiệm vụ Nghiên cứu phản ứng động đất tòa nhà tầng nhà có khung moment thép đặc biệt (Special Moment Resisting Frame - SMRF) với dầm có nhiều trƣờng hợp mặt cắt giảm yếu khác (Reduced Beam Section – RBS) Nội dung 1) Khảo sát mơ hình phƣơng pháp phân tích động, phạm vi khảo sát ngồi miền đàn hồi với tổng số lần khảo sát 660 lần 2) Từ trƣờng hợp giảm yếu ta tìm đặc trƣng dao động (gia tốc, chuyển vị) tầng khác ứng với tác động dao động trận động đất khác nhau, đƣợc mơ phần mềm chun dụng OpenSees từ đánh giá phân tích 3) Xác định phân tích lực cắt tầng theo tuyến tính phi tuyến sau so sánh kết tính tốn với tiêu chuẩn TCVN 9368:2012, ASCE 7-10 4) Tìm lực tới hạn khung moment thép đặc biệt (Special Moment Resisting Frame - SMRF) qua lần giảm yếu mặt cắt dầm iii 5) Kết luận kiến nghị giới hạn nghiên cứu III Ngày giao nhiệm vụ: 17/3/2015 IV Ngày hồn thành nhiệm vụ: 17/9/2015 V Cán hƣớng dẫn: Tiến sĩ Đào Đình Nhân CÁN BỘ HƢỚNG DẪN KHOA QUẢN LÝ CHUN NGÀNH (Họ tên chữ ký) (Họ tên chữ ký) iv LỜI CAM ĐOAN Tơi cam đoan cơng trình nghiên cứu riêng tơi Các số liệu, kết nêu Luận văn trung thực chƣa đƣợc cơng bố cơng trình khác Đồng thời, thơng tin trích dẫn Luận văn đƣợc tơn trọng đƣợc rõ nguồn gốc Tác giả NGUYỄN VĂN THI v LỜI CẢM ƠN Cơng trình nghiên cứu đƣợc thực mƣời tám tháng, q trình thực đề tài, khơng có giúp đỡ thầy Đào Đình Nhân bạn nhóm khơng có kết Em xin chân thành gửi lời cảm ơn đến thầy Khoa xây dựng nhƣ thầy dạy em mơn học vừa qua Đặc biệt, luận văn em xin chân thành cảm ơn thầy Đào Đình Nhân Cảm ơn thầy với lòng nhiệt huyết, thầy truyền đạt kiến thức kinh nghiệm bổ ích để em hồn thành viết Em mong góp ý chân thành đến từ Hội đồng khoa học để luận đƣợc hiệu chỉnh hồn thiện NGUYỄN VĂN THI vi TĨM TẮT LUẬN VĂN Đề tài : Ảnh hưởng bề rộng tiết diện dầm giảm yếu đến phản ứng với động đất khung thép đặc biệt ba tầng Từ khóa: Dầm giảm yếu tiết diện, động đất, khung mơ men, phi tuyến Nội dung nghiên cứu Nghiên cứu phản ứng động đất tòa nhà tầng nhà có khung mơ men thép đặc biệt có độ lệch tâm khối lƣợng 0,05 theo phƣơng cạnh dài Phản ứng kết cấu ứng với nhiều trƣờng hợp giảm yếu tiết diện dầm khác đƣợc khảo sát Mơ hình kết cấu đƣợc phân tích phƣơng pháp phân tích động miền thời gian, phạm vi khảo sát ngồi miền đàn hồi Tổng số lần phân tích đáp ứng động kết cấu với băng gia tốc 1320 lần Cơng trình dùng cho khảo sát đƣợc thiết kế theo tiêu chuẩn Mỹ Mơ hình 3D kết cấu đƣợc xây dựng phần mềm mơ chun dụng cho ngành Kỹ Thuật Động Đất, phần mềm OpenSees Mơ hình đƣợc phân tích cho 60 băng gia tốc đại diện cho cấp độ động đất khác Các đại lƣợng phản ứng đƣợc xem xét gồm gia tốc sàn, chuyển vị tầng tƣơng đối lực cắt tầng Kết khảo sát cho thấy việc giảm yếu bề rộng tiết diện dầm đến 50% khơng ảnh hƣởng đáng kể đến phản ứng động cơng trình tất cấp động đất khảo sát Cụ thể, đáp ứng với cấp động đất có chu kỳ lặp 2475 năm, gia tốc sàn trung bình giảm 15%, 12.7% 9% sàn 2, mái; chuyển vị tầng tƣơng đối giảm tầng 14%, tăng dần tầng tầng lần lƣợt 7%, 2% Khi tiết diện dầm giảm yếu xuống 50%, cƣờng độ cơng trình giảm xuống 80% so với khơng giảm tiết diện Kết khảo sát cho thấy giảm yếu tiết diện bé 25%, cƣờng độ kết cấu giảm chậm Ngồi ra, phần khảo sát luận văn cho thấy vii cơng thức phân phối lực tĩnh ngang tƣơng đƣơng động đất lên sàn tiêu chuẩn Mỹ tiêu chuẩn Việt Nam phù hợp tốt với phản ứng động kết cấu làm việc tuyến tính Khi kết cấu làm việc phi tuyến, việc phân phối lực tĩnh ngang theo tiêu chuẩn có phần sai lệch so với kết phản ứng động phân tích theo lịch sử thời gian Tuy nhiên, việc sai lệch tiêu chuẩn thiên an tồn viii ABSTRACT Thesis’s name: Effect of flange width of reduced beam section to the seismic responses of a 3-story steel special moment frame building Keywords: Reduced beam section, earthquake, moment frame, nonliear Contents: The seismic responses of a 3-story steel special moment frame building with 5% mass eccentricity was studied The structural responses corresponding to several reduced flange width of the I-shape beam was investigated Both linear and nonliear responses of the structural model were dynamically analyzed in time domain A total number of 1320 trials was dynamically analyzed The analyzed building was code-designed per USA standards The 3D model of the building was developed in OpenSees, a software for earthquake engineering simulation 60 pairs of ground motions representing earthquake levels was applied to the model The structural responses in consideration include floor acceleration, story drift and story shear The analysis results show that the reduction of flange width of beams up to 50% does not significantly affect the structural responses Specifically, in the responses to 2475 years return period earthquakes, floor acceleration decreases 15%, 12,7% and 9% at floors 2, and roof, respectively The story drift decreases 14% at story 1, increases 7% and 2% at story and When the flange width reduces to 50%, the lateral strength of the structure reduces to 80% of the strength of the full beam section structure It was observed that the lateral strength reduces slowly when the reduction of beam width is smaller than 25% Another important part of the investigation was to compare the equivalent lateral earthquake load distribution computed from codes to the effective lateral load analyzed from time history responses of the model The analysis indicates that the equivalent lateral load distribution computed from Vietnamese code and the ASCE - 109 - for m=1:1:length(Level) StAcc=MeanAcc(:,(m-1)*length(BHinge)+1:1:m*length(BHinge)); plot(RBS,StAcc,LineColorBHinge{m}); end legend('Tang 0','Tang 1','Tang 2','Tang 3','Location','SouthEast'); xlabel('Hinge'); ylabel('StAcc'); title('GIA TOC TANG CHO MOI LAN GIAM TIET DIEN'); Print_meanAcc = ['B',num2str(B),'_Hinge',num2str(BHinge(m)),'_StAcc.eps']; print('-depsc2',Print_meanAcc) hold off - 110 - PHỤ LỤC MÃ LẬP TRÌNH MATLAB: TÍNH TỐN VÀ VẼ BIỂU ĐỒ CHUYỂN VỊ clear all; clc; %% Cac cap dong dat B = ; % Nhap cap dong dat = [1 3] if B == YearNum = 72; elseif B ==2 YearNum = 475; else YearNum = 2475; end %% Nhap thong tin %GM={'BO0910','NGA1045','BO1516','NGA0884','NGA1084','NGA1054','NG A1063','NGA0173','NGA1116','NGA1752','NGA0728','NGA0169','NGA1044' ,'NGA0721','NGA0313','NGA1082','NGA0367','NGA0322','NGA0368','NGA0 558',}; %GM={'NGA1084','NGA1106','NGA0723','NGA1054','NGA0569','NGA0721',' NGA1054','NGA0821','NGA1141','NGA1116','la01','la02','la03','la04' ,'la05','la06','la07','la08','la09','la10'} GM={'NGA1084','NGA0527','NGA1106','NGA1054','NGA1044','NGA1120','N GA0821','NGA1605','NGA1141','bo31','se01','se02','se03','se04','se 05','se06','se07','se08','se09','se10'} BHinge = [7.75 8.525 9.3 10.075 10.85 11.625 12.4 13.175 13.95 14.725 15.5]; RBS = [50 45 40 35 30 25 20 15 10 ]; RBS1={'50%','45%','40%','35%','30%','25%','20%','15%','10%','5%',' 0%'}; Level = [0 3]; for i=1:1:length(BHinge) for j=1:1:length(GM) FAFile = ['F:\Output\B3_',GM{j},'_Hinge',num2str(BHinge(i)),'_Disp.OUT']; dFL = load(FAFile); RelDisp = dFL; RelDisp(:,1)=[]; PeakDriftL1_x(j,i)= max(max(abs(RelDisp(:,1:3:18)))); PeakDriftL2_x(j,i)= max(max(abs(RelDisp(:,19:3:36)RelDisp(:,1:3:18)))); PeakDriftL3_x(j,i)= max(max(abs(RelDisp(:,37:3:54)RelDisp(:,19:3:36)))); PeakDrift_x = [PeakDriftL1_x*0, PeakDriftL1_x, PeakDriftL2_x, PeakDriftL3_x]; - 111 - MeanDriftL1_x= mean(PeakDriftL1_x); % Gia tri trung binh MeanDriftL2_x= mean(PeakDriftL2_x); MeanDriftL3_x= mean(PeakDriftL3_x); MeanDrift_x = [MeanDriftL1_x*0, MeanDriftL1_x, MeanDriftL2_x, MeanDriftL3_x]; StdDriftL1_x StdDriftL2_x StdDriftL3_x StdDrift_x = StdDriftL3_x]; = std(PeakDriftL1_x); % Dolech chuan = std(PeakDriftL2_x); = std(PeakDriftL3_x); [MeanDriftL1_x*0, StdDriftL1_x, StdDriftL2_x, SigMeanDrift_x = StdDrift_x + MeanDrift_x; %Gia tri trung binh + lech chuan %% Create matrix of drift along y global axis in every story PeakDriftL1_y(j,i)= max(max(abs(RelDisp(:,2:3:18)))); PeakDriftL2_y(j,i)= max(max(abs(RelDisp(:,20:3:36)RelDisp(:,2:3:18)))); PeakDriftL3_y(j,i)= max(max(abs(RelDisp(:,38:3:54)RelDisp(:,20:3:36)))); PeakDrift_y = [PeakDriftL1_y*0, PeakDriftL1_y, PeakDriftL2_y, PeakDriftL3_y]; MeanDriftL1_y = mean(PeakDriftL1_y); % Gia tri trung binh MeanDriftL2_y = mean(PeakDriftL2_y); MeanDriftL3_y = mean(PeakDriftL3_y); MeanDrift_y = [ MeanDriftL1_y*0, MeanDriftL1_y, MeanDriftL2_y, MeanDriftL3_y]; StdDriftL1_y = std(PeakDriftL1_y); % Dolech chuan StdDriftL2_y = std(PeakDriftL2_y); StdDriftL3_y = std(PeakDriftL3_y); StdDrift_y = [StdDriftL1_y*0, StdDriftL1_y, StdDriftL2_y, StdDriftL3_y]; SigMeanDrift_y = StdDrift_y + MeanDrift_y; %Gia tri trung binh + lech chuan %% Create matrix of torque about z global axis in every story PeakTorL1_z(j,i)= max(max(abs(RelDisp(:,3:3:18)))); PeakTorL2_z(j,i)= max(max(abs(RelDisp(:,21:3:36)RelDisp(:,3:3:18)))); PeakTorL3_z(j,i)= max(max(abs(RelDisp(:,39:3:54)RelDisp(:,21:3:36)))); PeakTor_z = [PeakTorL1_z*0, PeakTorL1_z, PeakTorL2_z, PeakTorL3_z]; MeanTorL1_z = mean(PeakTorL1_z); % Gia tri trung binh MeanTorL2_z = mean(PeakTorL2_z); MeanTorL3_z = mean(PeakTorL3_z); MeanTor_z = [MeanTorL1_z*0, MeanTorL1_z, MeanTorL2_z, MeanTorL3_z]; StdTorL1_z = std(PeakTorL1_z); % Dolech chuan - 112 - StdTorL2_z = std(PeakTorL2_z); StdTorL3_z = std(PeakTorL3_z); StdTor_z = [StdTorL1_z*0, StdTorL1_z, StdTorL2_z, StdTorL3_z]; SigMeanTor_z = StdTor_z + MeanTor_z; %Gia tri trung binh + lech chuan end end %% Plot LineColorDash ={'g-','b-','c-','m-','y-','k-','r ','g ','b-','c ','m ','y ','k ','r-o','g-o','b-o','c-o','m-o','m-','m-'}; LineWidth = [1.0 1.2 1.4 1.6 1.8 2.0 0.8 1.0 1.2 1.4 1.6 1.8 2.0 0.8 1.0 1.2 1.4 1.6 1.8 2.0 ]; %% Plot elastic responses along x global axis for m=1:1:length(BHinge) figure hold on peakDrift_x=PeakDrift_x(:,m:length(BHinge):length(PeakDrift_x)); for nGM=1:1:length(LineColorDash) plot(peakDrift_x(nGM,:),Level,LineColorDash{nGM}); title(sprintf('CAP DO RBS X %.2f',BHinge(m))); xlabel('Story Drift (in)'); ylabel('Story Level'); legend(GM,'Location','SouthEast'); end Print_Drift_x = ['B',num2str(B),'_Hinge',num2str(BHinge(m)),'_Drift_x.eps']; print('-depsc',Print_Drift_x) hold off end %% Plot elastic responses along y global axis for m=1:1:length(BHinge) figure hold on peakDrift_y=PeakDrift_y(:,m:length(BHinge):length(PeakDrift_y)); for nGM=1:1:length(LineColorDash) plot(peakDrift_y(nGM,:),Level,LineColorDash{nGM}); title(sprintf('CAP DO RBS Y %.2f',BHinge(m))); xlabel('Story Drift (in)'); ylabel('Story Level'); legend(GM,'Location','SouthEast'); end Print_Drift_y = ['B',num2str(B),'_Hinge',num2str(BHinge(m)),'_eDrift_y.eps']; print('-depsc',Print_Drift_y) hold off end %% Plot elastic torque about z global axis - 113 - for m=1:1:length(BHinge) figure hold on peakTor_z=PeakTor_z(:,m:length(BHinge):length(PeakTor_z)); for nGM=1:1:length(LineColorDash) plot( peakTor_z(nGM,:),Level,LineColorDash{nGM}); title(sprintf('CAP DO RBS Z %.2f',BHinge(m))); xlabel('Story Drift (in)'); ylabel('Story Level'); legend(GM,'Location','SouthEast'); end Print_Tor_z = ['B',num2str(B),'_Hinge',num2str(BHinge(m)),'_Tor_z.eps']; print('-depsc2',Print_Tor_z) hold off end %% Plot the average value of responses - x global axis for m=1:1:length(BHinge) figure hold on meanDrift_x=MeanDrift_x(:,m:length(BHinge):length(MeanDrift_x)); sigMeanDrift_x=SigMeanDrift_x(:,m:length(BHinge):length(SigMeanDri ft_x)); for nBH=1:1:length(LineColorDash) plot(meanDrift_x,Level,LineColorDash{1}), grid, ; plot(sigMeanDrift_x,Level,LineColorDash{2}); xlabel('Story Drift (in)'); ylabel('Story Level'); title(sprintf('CAP DO RBS X %.2f',BHinge(m))); legend('meanDrift_x','sigMeanDrift_x','Location','SouthEast'); end hold off Print_meanDrift_x = ['B',num2str(B),'_Hinge',num2str(BHinge(m)),'_MeanDrift_SigMeanDru ff_x.eps']; print('-depsc2',Print_meanDrift_x) end %% Plot the average value of responses - y global axis for m=1:1:length(BHinge) figure hold on meanDrift_y=MeanDrift_y(:,m:length(BHinge):length(MeanDrift_y)); sigMeanDrift_y=SigMeanDrift_y(:,m:length(BHinge):length(SigMeanDri ft_y)); for nBH=1:1:length(LineColorDash) plot(meanDrift_y,Level,LineColorDash{1}), grid, ; plot(sigMeanDrift_y,Level,LineColorDash{2}); - 114 - xlabel('Story Drift (in)'); ylabel('Story Level'); title(sprintf('CAP DO RBS Y %.2f',BHinge(m))); legend('meanDrift_y','sigMeanDrift_y','Location','SouthEast'); end hold off Print_meanDrift_y = ['B',num2str(B),'_Hinge',num2str(BHinge(m)),'_MeanDrift_SigMeanDru ff_y.eps']; print('-depsc2',Print_meanDrift_y) end %% Plot the average value of responses - about z global axis figure for m=1:1:length(BHinge) figure hold on meanTor_z=MeanTor_z(:,m:length(BHinge):length(MeanTor_z)); sigMeanTor_z=SigMeanTor_z(:,m:length(BHinge):length(SigMeanTor_z)) ; for nBH=1:1:length(LineColorDash) plot(meanDrift_x,Level,LineColorDash{1}), grid, ; plot(sigMeanDrift_x,Level,LineColorDash{2}); xlabel('Story Drift (in)'); ylabel('Story Level'); title(sprintf('CAP DO RBS Z %.2f',BHinge(m))); legend('meanDrift_z','sigMeanDrift_z','Location','SouthEast'); end hold off Print_meanTor_z = ['B',num2str(B),'_Hinge',num2str(BHinge(m)),'_MeanDrift_SigMeanDru ff_z.eps']; print('-depsc2',Print_meanTor_z) end %% Plot the average value of responses - about nstep global axis figure figure hold on for m=1:1:length(BHinge) Bhingh meanDrift_x=MeanDrift_x(:,m:length(BHinge):length(MeanDrift_x)); plot(meanDrift_x,Level,LineColorDash{m}); end xlabel('Story Drift (in)'); ylabel('Story Level'); title('CAP DO DONG DAT TUONG UNG VOI CAC LAN GIAM TIET DIEN'); legend(RBS1,'Location','SouthEast'); - 115 - Print_meanTor_z = ['B',num2str(B),'_Hinge',num2str(BHinge(m)),'_MeanDrift_x.eps']; print('-depsc2',Print_meanDrift_x) hold off %% Plot the average value of responses - about nstep Bhingh global axis figure hold on for m=2:1:length(Level) StrDrift_x= MeanDrift_x(:,(m1)*length(BHinge)+1:1:m*length(BHinge)); plot(RBS,StrDrift_x,LineColorDash{m}); end legend('Tang 1','Tang 2','Tang 3','Location','SouthEast'); xlabel('Hinge'); ylabel('StrDrift'); title('CHUYEN VI TANG CHO MOI LAN GIAM TIET DIEN'); hold off %% figure hold on for m=1:1:length(BHinge) meanDrift_y=MeanDrift_y(:,m:length(BHinge):length(MeanDrift_y)); plot(meanDrift_y,Level,LineColorDash{m}); end xlabel('Story Drift (in)'); ylabel('Story Level'); title('CAP DO DONG DAT TUONG UNG VOI CAC LAN GIAM TIET DIEN'); legend(RBS1,'Location','SouthEast'); Print_meanTor_z = ['B',num2str(B),'_Hinge',num2str(BHinge(m)),'_MeanDrift_y.eps']; print('-depsc2',Print_meanDrift_x) hold off %% Plot the average value of responses - about nstep Bhingh global axis figure hold on for m=2:1:length(Level) StrDrift_y=MeanDrift_y(:,(m1)*length(BHinge)+1:1:m*length(BHinge)); plot(RBS,StrDrift_y,LineColorDash{m}); end legend('Tang 1','Tang 2','Tang 3','Location','SouthEast'); xlabel('Hinge'); ylabel('StrDrift'); title('CHUYEN VI TANG CHO MOI LAN GIAM TIET DIEN'); hold off %% figure hold on - 116 - for m=1:1:length(BHinge) meanDrift_z=MeanTor_z(:,m:length(BHinge):length(MeanTor_z)); plot(meanDrift_z,Level,LineColorDash{m}); end xlabel('Story Drift (in)'); ylabel('Story Level'); title('CAP DO DONG DAT TUONG UNG VOI CAC LAN GIAM TIET DIEN'); legend(RBS1,'Location','SouthEast'); hold off %% Plot the average value of responses - about nstep Bhingh global axis figure hold on for m=2:1:length(Level) StrDrift_z=MeanTor_z(:,(m1)*length(BHinge)+1:1:m*length(BHinge)); plot(RBS,StrDrift_z,LineColorDash{m}); end legend('Tang 1','Tang 2','Tang 3','Location','SouthEast'); xlabel('Hinge'); ylabel('StrDrift'); title('CHUYEN VI TANG CHO MOI LAN GIAM TIET DIEN'); hold off - 117 - PHỤ LỤC MÃ LẬP TRÌNH MATLAB: VẼ BIỂU LỰC CẮT TẦNG PP SỐ, TCVN, ASCE Level = [1,2,3]; %Luc cat, luc tinh x cap dong %% PHI TUYEN DataIne = [DigMethod; TCVN; ASCE] %Ve luc cat % Phuong X DataXIne = [2744 2306 1346; 438 960 1346; 4919 4093 2420; 826 1673 2420; 5924 4677 2382; 1246 2295 2382; 2744 2295 1402; 448 894 1402; 4919 4115 2513; 804 1602 2513; 5924 4956 3026; 968 1930 3026; 2744 2338 1456; 406 882 1456; 4919 4191 2610; 728 1582 2610; 5924 5047 3142; 876 1905 3142]; DataYIne = [2669 2281 1302; 387 980 1302; 4444 3765 2181; 679 1585 2181; 5514 4373 2420; 1141 1952 2420; 2669 2232 1363; 436 869 1363; 4444 3718 2270; 726 1448 2270; 5514 4613 2817; 901 1796 2817; 2669 2271 1411; 398 859 1411; 4444 3781 2350; 663 1431 2350; 5514 4692 2916; 822 1775 2916]; figure; subplot(1,3,1); barh(Level,DataXIne(1:6:18,:)','grouped'); % groups by row title('72 năm') subplot(1,3,2); barh(Level,DataXIne(3:6:18,:)','grouped'); % groups by row title('475 năm') subplot(1,3,3); barh(Level,DataXIne(5:6:18,:)','grouped'); % groups by row title('2.475 năm') legend('PP số', 'TCVN', 'ASCE','Location','SouthEast'); % Ve luc tinh - Phuong X figure; subplot(1,3,1); barh(Level,DataXIne(2:6:18,:)','grouped'); % groups by row title('72 năm') subplot(1,3,2); barh(Level,DataXIne(4:6:18,:)','grouped'); % groups by row title('475 năm') subplot(1,3,3); barh(Level,DataXIne(6:6:18,:)','grouped'); % groups by row title('2.475 năm') legend('PP số', 'TCVN', 'ASCE','Location','SouthEast'); %% Ve luc cat % Phuong Y figure; subplot(1,3,1); barh(Level,DataYIne(1:6:18,:)','grouped'); % groups by row title('72 năm') - 118 - subplot(1,3,2); barh(Level,DataYIne(3:6:18,:)','grouped'); % groups by row title('475 năm') subplot(1,3,3); barh(Level,DataYIne(5:6:18,:)','grouped'); % groups by row title('2.475 năm') legend('PP số', 'TCVN', 'ASCE','Location','SouthEast'); %% Ve luc tinh - Phuong Y figure; subplot(1,3,1); barh(Level,DataYIne(2:6:18,:)','grouped'); % groups by row title('72 năm'); subplot(1,3,2); barh(Level,DataYIne(4:6:18,:)','grouped'); % groups by row title('475 năm'); subplot(1,3,3); barh(Level,DataYIne(6:6:18,:)','grouped'); % groups by row title('2.475 năm') legend('PP số', 'TCVN', 'ASCE','Location','SouthEast'); %% TUYEN TINH, DataEla = [DigMethod; TCVN; ASCE] %Ve luc cat % Phuong X DataXEla = [2885 2447 1463;437 984 1463;5578 4719 2823;859 1896 2823;11190 9434 5640;1756 3794 5640;2885 2413 1474;471 940 1474;5578 4667 2850;911 1817 2850;11190 9361 5716;1828 3645 5716;2885 2458 1530;427 928 1530;5578 4753 2959;825 1794 2959;11190 9535 5936;1655 3599 5936;]; DataYEla = [2723 2325 1321;398 1004 1321;5244 4497 2579;747 1917 2579;8620 7387 4252;1233 3135 4252;2723 2278 1391;445 887 1391;5244 4387 2679;857 1708 2679;8620 7211 4403;1408 2808 4403;2723 2317 1440;406 877 1440;5244 4462 2774;782 1688 2774;8620 7334 4559;1285 2775 4559]; figure; subplot(1,3,1); barh(Level,DataXEla(1:6:18,:)','grouped'); % groups by row title('72 năm') subplot(1,3,2); barh(Level,DataXEla(3:6:18,:)','grouped'); % groups by row title('475 năm') subplot(1,3,3); barh(Level,DataXEla(5:6:18,:)','grouped'); % groups by row title('2.475 năm') legend('PP số', 'TCVN', 'ASCE','Location','SouthEast'); % Ve luc tinh - Phuong X figure; subplot(1,3,1); barh(Level,DataXEla(2:6:18,:)','grouped'); % groups by row title('72 năm') subplot(1,3,2); barh(Level,DataXEla(4:6:18,:)','grouped'); % groups by row - 119 - title('475 năm') subplot(1,3,3); barh(Level,DataXEla(6:6:18,:)','grouped'); % groups by row title('2.475 năm') legend('PP số', 'TCVN', 'ASCE','Location','SouthEast'); %% Ve luc cat % Phuong Y figure; subplot(1,3,1); barh(Level,DataYEla(1:6:18,:)','grouped'); % groups by row title('72 năm') subplot(1,3,2); barh(Level,DataYEla(3:6:18,:)','grouped'); % groups by row title('475 năm') subplot(1,3,3); barh(Level,DataYEla(5:6:18,:)','grouped'); % groups by row title('2.475 năm') legend('PP số', 'TCVN', 'ASCE','Location','SouthEast'); %% Ve luc tinh - Phuong Y figure; subplot(1,3,1); barh(Level,DataYEla(2:6:18,:)','grouped'); % groups by row title('72 năm'); subplot(1,3,2); barh(Level,DataYEla(4:6:18,:)','grouped'); % groups by row title('475 năm'); subplot(1,3,3); barh(Level,DataYEla(6:6:18,:)','grouped'); % groups by row title('2.475 năm') legend('PP số', 'TCVN', 'ASCE','Location','SouthEast'); - 120 - PHỤ LỤC 4: MÃ LẬP TRÌNH OPENSEES KHẢO SÁT HIỆN TƢỢNG ĐẨY DẦN (PUSHOVER) # EARTHQUAKE RESPONSE OF THREE-STORY BUILDING # CREATED BY: NGUYEN VAN THI # PURPOSE: DOING RESEARCH FOR THESIS OF MASTER OF SCIENCE IN CIVIL CONSTRUCTION # UNIT: [F] = kip, [L] = in ############################################################################### wipe; # clear previous model set BHinge1 "14.75"; # Tiet dien giam yeu Girder set BHinge2 "14.75"; # Tiet dien giam yeu Girder set BHinge3 "9.5"; # Tiet dien giam yeu Girder set e 0.05; # Define loops to analyze for {set iHinge 0} {$iHinge < [llength $BHinge1]} {incr iHinge} { set b1 [lindex $BHinge1 $iHinge]; set b2 [lindex $BHinge2 $iHinge]; set b3 [lindex $BHinge3 $iHinge]; # initialize the model model BasicBuilder -ndm -ndf 6; # INPUT INFORMATION set BayX 360.; # Kich thuoc nhip theo phuong truc X set BayY 360.; # Kich thuoc nhip theo phuong truc Y set HStory 180.; # Chieu cao tang set NFloor 4.; # level (3 san) set NBayX 6.; # co nhip theo phuong X set NBayY 4.; # co nhip theo phuong Y set LCol 17; # Chieu dai doan khop o cot khung moment tuong ung voi chieu cao tiet dien cot set LHinge1 20; # Chieu dai doan khop Girder tang set LHinge2 20; # Chieu dai doan khop Girder tang set LHinge3 15; # Chieu dai doan khop Girder tang set Lij1 7; # Chieu dai dau I J cua Girder tang set Lij2 7; # Chieu dai dau I J cua Girder tang set Lij3 6; # Chieu dai dau I J cua Girder tang #set LSmall 2; # Chieu dai khop dam phu # -source Procedure.tcl; - 121 - source Node.tcl; source Material.tcl; source Section.tcl; source Element.tcl; source Boundary.tcl; # -# Static Analysis puts "Static Analysis "; #file mkdir GravityAnalysis source Load.tcl; system UmfPack -lvalueFact 50; # method of solving system of linear equations constraints Transformation; # Create the constraint handler numberer Plain; # method of numbering test NormDispIncr 1e-5 50; # criterial for checking convergence algorithm Newton; # algorithm for solving nonlinear problem integrator LoadControl 1.; # method for predicting next step analysis Static; # type of problem | Transient/VariableTransient/Static/ analyze 1; # perform analysis | analyze $numIncr loadConst -time 0.; # Define recoder file mkdir Output recorder Node -file Output/MasterNodeDisp.out -time -node -dof disp puts "Pushover Analysis " pattern Plain Linear { load [expr 1246./(1246.+2295.+2382.)] 0 0 load [expr 2295./(1246.+2295.+2382.)] 0 0 load [expr 2382./(1246.+2295.+2382.)] 0 0 } # Define analysis object system UmfPack -lvalueFact 50; # method of solving system of linear equations test NormDispIncr 1e-5 20; # criterial for checking convergence algorithm Newton; # algorithm for solving nonlinear problem integrator DisplacementControl 0.1 ; # method for predicting next step analysis Static; # type of problem | Transient/VariableTransient/Static/ analyze 300; # perform analysis | analyze $numIncr wipe; }; # End Main - 122 - PHỤ LỤC MÃ LẬP TRÌNH MATLAB: BIỂU DIỄN LỰC GÂY SỤP ĐỔ VÀ CHUYỂN VỊ TẦNG clear all clc %% Nhap thong tin %GM={'BO0910','NGA1045','BO1516','NGA0884','NGA1084','NGA1054','NG A1063','NGA0173','NGA1116','NGA1752','NGA0728','NGA0169','NGA1044' ,'NGA0721','NGA0313','NGA1082','NGA0367','NGA0322','NGA0368'}; %GM={'NGA1084','NGA1106','NGA0723','NGA1054','NGA0569','NGA0721',' NGA1054','NGA0821','NGA1141','NGA1116','la01','la02','la03','la04' ,'la05','la06','la07','la08','la09','la10'} %GM={'NGA1084','NGA0527','NGA1106','NGA1054','NGA1044','NGA1120',' NGA0821','NGA1605','NGA1141','bo31','se01','se02','se03','se04','s e05','se06','se07','se08','se09','se10'} BHinge = [7.75 8.525 9.3 10.075 10.85 11.625 12.4 13.175 13.95 14.725 15.5]; RBS={'50%','45%','40%','35%','30%','25%','20%','15%','10%','5%','0 %'}; g = 386.22; Level = [1 3]; LineColorDash={'g-','b-','c-','m-','y-','k-','r ','g ','b ','c-','m ','y ','k ','r-o','g-o','b-o','c-o','m-o','m-','m '}; LineWidth = [1.0 1.2 1.4 1.6 1.8 2.0 0.8 1.0 1.2 1.4 1.6 1.8 2.0 0.8 1.0 1.2 1.4 1.6 1.8 2.0 ]; DatMatrix=zeros(300,length(BHinge)*4); for iBHinge=1:1:length(BHinge) FAFile=['I:\HOCTAP\luanvan\baithaysua\StoreyMomentFrame Hinges20 15.09.12\Output\MasterNodeDisp_Hinge',num2str(BHinge(iBHinge)),'.O UT']; Disp = load(FAFile); DatMatrix(:,(iBHinge-1)*4+1)=Disp(:,1); DatMatrix(:,(iBHinge-1)*4+2)=abs(Disp(:,2)); DatMatrix(:,(iBHinge-1)*4+3)=abs(Disp(:,5)-Disp(:,2)); DatMatrix(:,(iBHinge-1)*4+4)=abs(Disp(:,8)-Disp(:,5)); end for iSt=1:1:length(Level) - 123 - figure hold on for iCol=1:1:11 plot(DatMatrix(:,(iCol1)*4+iSt+1),DatMatrix(:,(iCo1)*4+1)*100/(12*15),LineColorDash{iCol},'linewidth,Line Width(iCol)); end title(sprintf('Tang %d',iSt)); xlabel('Story Drift (%)'); ylabel('Force'); legend(RBS,'Location','SouthEast'); hold off end RBS1= [0 10 15 20 25 30 35 40 45 50]; for iSt=1:1:4 figure hold on plot(RBS1,P(iSt,:),'linewidth',LineWidth(iCol)); title(sprintf('Tang %d',iSt)); xlabel('Story Drift (%)'); ylabel('Force'); legend('tầng1','tầng2','tầng3','Location','SouthEast'); end hold off