Phân tích khung phẳng liên hợp thép bê tông cốt thép liên kết nửa cứng chịu tải trọng tỉnh và động

ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA YZ ĐỖ VIỆT HƯNG PHÂN TÍCH KHUNG PHẲNG LIÊN HP THÉP – BÊ TÔNG CỐT THÉP LIÊN KẾT NỬA CỨNG CHỊU TẢI TRỌNG TĨNH VÀ ĐỘNG CHUYÊN NGÀNH : XÂY DỰNG DÂN DỤNG VÀ CÔNG NGHIỆP MÃ SỐ NGÀNH : 23.04.10 LUẬN VĂN THẠC SĨ TP.HỒ CHÍ MINH – THÁNG 3.2007 ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA YZ ĐỖ VIỆT HƯNG PHÂN TÍCH KHUNG PHẲNG LIÊN HP THÉP – BÊ TÔNG CỐT THÉP LIÊN KẾT NỬA CỨNG CHỊU TẢI TRỌNG TĨNH VÀ ĐỘNG CHUYÊN NGÀNH : XÂY DỰNG DÂN DỤNG VÀ CÔNG NGHIỆP MÃ SỐ NGÀNH : 23.04.10 LUẬN VĂN THẠC SĨ NGƯỜI HƯỚNG DẪN KHOA HỌC PGS.TS BÙI CÔNG THÀNH TP.HỒ CHÍ MINH – THÁNG 3.2007 i CÔNG TRÌNH ĐƯC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA - ĐẠI HỌC QUỐC GIA TP.HCM Cán hướng dẫn khoa học : PGS.TS BÙI CÔNG THÀNH 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ệ : HỘI ĐỒNG CHẤM BẢO VỆ LUẬN VĂN THẠC SĨ TRƯỜNG ĐẠI HỌC BÁCH KHOA, ngày ……………………tháng ……………………năm 2007 ii ĐẠI HỌC QUỐC GIA TP.HỒ CHÍ MINH CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập – Tự – Hạnh phúc TRƯỜNG ĐẠI HỌC BÁCH KHOA NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên Ngày tháng năm sinh Chuyên ngành : ĐỖ VIỆT HƯNG : 27/04/1978 : XÂY DỰNG DD&CN Phái: Nam Nôi sinh : Gia Lai MSHV : 02104547 I - TÊN ĐỀ TÀI : PHÂN TÍCH KHUNG PHẲNG LIÊN HP THÉP – BÊ TÔNG CỐT THÉP LIÊN KẾT NỬA CỨNG CHỊU TẢI TRỌNG TĨNH VÀ ĐỘNG II - NHIỆM VỤ VÀ NỘI DUNG - Trình bày mô hình liên kết nửa cứng tuyến tính, Richard-Abbott, Eurocode 4, cách xác định độ cứng liên kết - Nghiên cứu, áp dụng khái niệm hệ số ngàm-đầu mút ma trận độ cứng hiệu chỉnh phần tử dầm-cột cứng thành phần tử dầm-cột nửa cứng sử dụng cho khung thép phẳng để phân tích khung liên hợp có liên kết nửa cứng Trình bày lý thuyết động lực học kết cấu lựa chọn phương pháp giải toán động lực học cho kết cấu khung phẳng liên hợp thép – bê tông cốt thép - Xây dựng chương trình tính toán ngôn ngữ lập trình Matlab 7.0 Sử dụng chương trình để khảo sát ví dụ đồng thời so sánh với nghiên cứu có trước nhằm kiểm tra tính đắn hiệu chương trình, từ rút kết luận kiến nghị III- NGÀY GIAO NHIỆM VỤ: IV- NGÀY HOÀN THÀNH NHIỆM VỤ: V- HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN CÁN BỘ HƯỚNG DẪN Ngày 03 tháng 07 năm 2003 Ngày 05 tháng 03 năm 2007 CHỦ NHIỆM NGÀNH BỘ MÔN QUẢN LÝ NGÀNH PGS.TS BÙI CÔNG THÀNH Nội dung đề cương luận văn thạc só Hội Đồng Chuyên Ngành thông qua Ngày …… tháng…… năm 2007 PHÒNG ĐÀO TẠO SAU ĐẠI HỌC KHOA QUẢN LÝ NGÀNH iii LỜI CẢM ƠN Tôi xin chân thành cám ơn Ban Giám hiệu, Phòng Đào tạo Sau Đại học, Thầy, Cô giảng dạy sau đại học Trường đại học Bách khoa Tp Hồ Chí Minh kiến thức mà có ngày hôm Đặc biệt, xin bày tỏ lòng biết ơn sâu sắc đến Thầy PGS.TS Bùi Công Thành, người Thầy trực tiếp hướng dẫn nhiệt tình, tận tâm cho lời khuyên quý báu, định hướng giúp vượt qua khó khăn để hoàn thành luận văn Cám ơn quan công tác tạo điều kiện cho trình học tập Cuối cùng, xin gởi lòng biết ơn đến Cha, Mẹ, người thân bạn bè tôi, người quan tâm động viên suốt trình học tập thực luận văn iv Chương TOÅNG QUAN .4 1.1 Liên kết nửa cứng 1.2 Tình hình nghiên cứu kết cấu liên hợp thép - bê tông cốt thép 1.3 Sự cần thiết việc nghiên cứu phản ứng động lực học kết cấu thép - bê tông cốt thép liên hợp 1.4 Mục tiêu phạm vi đề tài 1.5 So saùnh nội dung luận văn số công trình nghiên cứu trước kết cấu khung liên hợp Chương MÔ HÌNH LIÊN KẾT NỬA CỨNG 10 2.1 Các mô hình tuyến tính 10 2.2 Mô hình thông số Richard-Abbott 11 2.3 Mô hình liên kết nửa cứng Eurocode .12 2.3.1 Phân loại liên kết liên hợp 13 2.3.2 Phương pháp thành phần .15 2.3.3 Mô hình liên kết liên hợp nửa cứng theo tiêu chuẩn Eurocode .15 2.3.4 Lý tưởng hóa đường cong phi tuyến M-ϕ 16 2.4 Độ cứng liên kết liên hợp: 17 2.4.1 Mô hình 18 2.4.2 Độ cứng thành phần liên hợp 18 2.4.2.1 Cốt thép gia cường .18 2.4.2.2 Chốt chống cắt .19 2.5 Đặc trưng tiết diện dầm liên hợp 19 2.5.1 Chiều rộng làm việc sàn .19 2.5.2 Mô men quán tính tiết diện dầm liên hợp 20 2.6 Kết luận 21 Chương PHÂN TÍCH KHUNG PHẲNG 22 3.1 Giới thiệu 22 3.2 Khái niệm hệ số ngàm-đầu mút 22 3.3 Phần tử dầm-cột nửa cứng 24 3.4 Ma trận độ cứng phần tử nửa cứng 25 3.5 Phân tích khung phẳng với tải trọng động 26 3.5.1 Cơ sở lý thuyết phân tích động lực học kết cấu 26 3.5.2 Trường chuyển vị phần tử nửa cứng 27 3.5.3 Xác định ma trận tính chất kết cấu .29 3.5.3.1 Ma trận khối lượng tương thích 29 3.5.3.2 Ma trận cản 30 3.5.4 Phản ứng hệ khung phẳng tác dụng tải trọng động .30 3.5.4.1 Dao động tự 30 3.5.4.2 Tải trọng điều hòa 30 3.5.4.3 Tải trọng xung 31 3.5.5 Phương pháp tích phân số Newmark 31 3.6 Phân tích khung không xét phi tuyến hình học 32 3.7 Phân tích khung có xét phi tuyến hình học 33 3.7.1 Ma trận độ cứng xét đến ảnh hưởng bậc phần tử dầm-cột nửa cứng .33 3.7.2 Phân tích bậc 34 3.8 Sơ đồ thuật toán phân tích khung phẳng liên hợp 35 3.9 Kết luận 37 Chương CHƯƠNG TRÌNH TÍNH TOÁN .38 4.1 Giới thiệu 38 4.2 Môđun tính toán liên kết 38 4.3 Phân tích khung phẳng .40 4.3.1 Phân tích với tải trọng tónh 40 4.3.2 Phân tích với tải trọng động 45 Chương KHẢO SÁT CÁC BÀI TOÁN VÍ DỤ 47 5.1 Tính toán liên kết liên hợp 47 5.2 Phân tích khung thép chịu tải trọng tónh 49 5.2.1 Phân tích khung thép liên kết cứng bậc 49 5.2.2 Phân tích bậc hai khung thép phẳng liên kết nửa cứng 52 5.3 Phân tích khung phẳng liên hợp chịu tải trọng tónh 55 5.3.1 Khung tầng nhịp 55 5.3.2 Khung tầng nhịp 60 5.4 Khảo sát dao động tự khung phẳng 63 5.4.1 Dao động tự khung thép phẳng .63 5.4.2 Dao động tự khung liên hợp 66 5.4.3 Khảo sát hội tụ toán dao động tự 69 5.4.3.1 Khung liên kết cứng 69 5.4.3.2 Khung liên kết nửa cứng 74 5.4.4 Ảnh hưởng tỉ số cản đến dao động tự khung phẳng 78 5.5 Phân tích khung liên hợp chịu tải trọng điều hòa 82 5.6 Phân tích khung liên hợp chịu tải trọng xung 86 Chương KẾT LUẬN VÀ KIẾN NGHỊ 90 6.1 Kết luận 90 6.2 Kiến nghị 91 TÀI LIỆU THAM KHẢO 92 PHUÏ LUÏC 94 composite.m Chuong trinh chinh - Composite 94 Phần A Tính toán liên kết 95 Phaàn A.1 Tính liên kết extended end-plate connections 95 Con1.m - tinh lien ket extended end-plate connections 95 Phần A.2 Tính liên kết partial-depth end-plate connections 108 Con2.m - tinh lien ket partial-depth end-plate connections 108 10 11 12 13 14 15 16 17 18 19 20 21 22 23 Phần B Tính toán khung phẳng 116 opendata.m - Nap du lieu bai toan tu mot file co san 116 newdata.m - tao mot bai toan moi 116 savedata_as.m - Luu toan bai toan va dat ten moi .116 savedata.m - luu du lieu cua bai toan dang thuc hien 117 frame.m - nhap so lieu khung .117 loaddata.m - nhap tai 121 dis_graphic.m - ve so hinh hoc khung phang 124 loadvecto.m tinh tai nut khung theo he so ngam dau mut r 131 solvestatic.m - tinh noi luc khung voi tai tinh .132 stiff1.m - tinh cung bac khung co lien ket cung hoac nua cung 134 stiff2.m - tinh cung bac khung co lien ket cung hoac nua cung 134 sematrix.m - matran tinh momen 135 dynamic_com.m - tinh khung phang voi tai dong 136 rigid_dynamic.m - Tinh khung cung voi tai dong 139 semi_dynamic.m - Tinh khung nua cung voi tai dong 142 integration.m - Giai phuong trinh vi phan dao dong 146 rayleigh.m - thiet lap ma tran can tong the 147 vibration.m - tinh tan so va chu ky dao dong rieng 147 mass.m - thiet lap ma tran khoi luong tong the 147 ouput.m - xuat ket qua tinh toan 148 Chương TỔNG QUAN 1.1 Liên kết nửa cứng Liên kết đóng vai trò quan trọng hệ thống kết cấu khung Liên kết nửa cứng loại liên kết mô gần trạng thái làm việc thực kết cấu Loại liên kết diễn tả trạng thái làm việc liên kết, tùy theo cách cấu tạo đường cong quan hệ mômen (M) góc xoay (θ) liên kết nửa cứng nằm đường cong M-ϕ hai loại liên kết thông thường khớp lý tưởng (không có ràng buộc góc xoay, mômen liên kết 0) liên kết cứng tuyệt đối (không xoay, mômen dầm truyền hoàn toàn sang cột) Wilson Moore người đưa nghiên cứu liên kết nửa cứng vào năm 1917 [2] Trong thời gian dài sau đó, nhiều nghiên cứu tập trung vào ứng xử liên kết nửa cứng Có thể kể đến khảo sát sau đây: G.R Monforton T.S Wu đưa cách phân tích khung có liên kết nửa cứng cách hiệu chỉnh ma trận độ cứng (1963) Frye Morris (1975) đề xuất mô hình đa thức cho số loại liên kết thủ tục phân tích phi tuyến liên kết nhằm xác định biến dạng khung tác dụng tải trọng, Lei Xu Donale E.Grierson (1992) giới thiệu phương pháp tự động hóa máy tính để thiết kế tối ưu khung thép có xét đến ảnh hưởng liên kết nửa cứng Liew J.Y.R cộng (1997) tiến hành thí nghiệm khung không giằng có liên kết nửa cứng làm sở so sánh tiến hành phân tích dẻo bậc hai; Yosuk Kim Wai - Fah Chen (1998) đề nghị phương pháp thực hành khác để thiết kế khung thép có liên kết nửa cứng sở quy phạm AISC-LRFD; G.C Clifton J.W Butterworth (1999), G.C Clifton, J.W Butterworth M Pantke (2001) đề nghị mô hình liên kết nửa cứng hệ khung chống động đất Hệ khung thép liên kết nửa cứng quy định quy phạm thiết kế Mỹ AISC-ASD (AISC, 1989), AISC-LRFD (1999), châu Âu Eurocode Trong phân tích khung phẳng, liên kết truyền nội lực từ dầm sang cột Khi tác dụng lên liên kết mô men M liên kết bị xoay góc θr Góc xoay biểu diễn thay đổi hình dạng ban đầu liên kết dầm cột θ M Hình 1.1 Biến dạng xoay liên kết Đường cong quan hệ mômen (M) góc xoay (θ) liên kết nửa cứng đường cong phi tuyến, mô thông qua hàm số mô men tác dụng biến dạng xoay Có nhiều mô hình liên kết nửa cứng đề xuất, việc lựa chọn mô hình ứng xử liên kết có ảnh hưởng lớn đến kết tính toán Ở Việt Nam khoảng thời gian gần đây, có nhiều nghiên cứu liên kết nửa cứng Tác giả Vũ Quốc Anh (ĐH Kiến trúc Hà Nội) trình bày báo “Stability analysis of steel frames with semi-rigid connections and semi-rigid connections with rigidzones” (2001) Một số luận văn cao học vấn đề liên kết nửa cứng thực ĐH Bách khoa Tp.HCM năm gần như: - Trần Tuấn Kiệt: “Phân tích khung thép phẳng có liên kết nửa cứng phương pháp nâng cao” (2002) - Nguyễn Trung: “Phân tích khung thép không gian có liên kết nửa cứng phương pháp nâng cao” (2003) - Trần Chí Hoàng: “Phân tích thiết kế tối ưu khung thép phẳng có liên kết nửa cứng với tiết diện chữ I” (2003) - Bùi Lâm: “Phân tích phản ứng khung thép phẳng liên kết nửa cứng chịu tải trọng động” (2004)… 1.2 Tình hình nghiên cứu kết cấu liên hợp thép - bê tông cốt thép Kết cấu liên hợp thép - bê tông cốt thép gần sử dụng rộng rãi xây dựng công trình cao tầng ngành cầu đường giới Sự hình thành kết cấu liên hợp dựa sở kết hợp ưu điểm vốn có loại vật liệu thép (bền, nhẹ, dễ lắp ráp, chịu kéo tốt, tốc độ thi công nhanh) bê tông (rẻ tiền, nặng, cứng, có khả chịu nén lớn) [20] Ngoài việc sử dụng kết cấu liên hợp thép - bê tông cốt thép liên hợp để tạo công trình cao tầng, kết cấu khung nhà có nhịp rộng, loại kết cấu thích hợp cho công trình xây dựng vùng động đất có khả chống cháy tốt [25] Như vậy, kết cấu liên hợp loại kết cấu có vị trí quan trọng ngành công nghiệp xây dựng Ở Việt Nam, thiết kế kết cấu liên hợp điều mẻ kỹ sư chưa có quy trình phân tích hoàn chỉnh quy phạm xây dựng Việt Nam chưa đề cập đến loại kết cấu Để áp dụng loại kết cấu có nhiều ưu điểm xây dựng công trình Việt Nam, việc nghiên cứu điều cần thiết Kết cấu thép - bê tông cốt thép liên hợp thật sử dụng cách hàng trăm năm Vào năm 1894, công trình sử dụng loại kết cấu cầu Rock Rapids, bang Iowa tòa nhà Methodist Pittsburgh, USA Những cấu kiện dầm liên hợp kiểm tra vào năm 1922 Canada, chốt hàn chống cắt thử nghiệm lần đầu trường đại học Illinois vào năm 1954 [13] Gi=N(i)/L(i)*[0 0 0 0; 6/5 L(i)/10 -6/5 L(i)/10; L(i)/10 2*L(i)^2/15 -L(i)/10 -L(i)^2/30; 0 0 0; -6/5 -L(i)/10 6/5 -L(i)/10; L(i)/10 -L(i)^2/30 -L(i)/10 2*L(i)^2/15]; %Xac dinh ma tran hieu chinh de ke den tinh nua cung cho ma tran cung hinh hoc cua phan tu dam-cot thu i Ci_g=[0 0 0 0; 0 0; -4/(5*L(i)*(4-r(i,1)*r(i,2))^2)*(8*r(i,1)^2*r(i,2)-13*r(i,1)*r(i,2)^2-32*r(i,1)^28*r(i,2)^2+25*r(i,1)*r(i,2)+20) r(i,1)/(5*(4-r(i,1)*r(i,2))^2)*(16*r(i,2)^2+25*r(i,1)*r(i,2)^296*r(i,1)*r(i,2)+128*r(i,1)-28*r(i,2)) 4/(5*L(i)*(4-r(i,1)*r(i,2))^2)*(8*r(i,1)^2*r(i,2)13*r(i,1)*r(i,2)^2-32*r(i,1)^2-8*r(i,2)^2+25*r(i,1)*r(i,2)+20) 4*r(i,2)/(5*(4r(i,1)*r(i,2))^2)*(16*r(i,1)^2-5*r(i,1)^2*r(i,2)+9*r(i,1)*r(i,2)-28*r(i,1)+8*r(i,2)); 0 0 0; 0 0 0; -4/(5*L(i)*(4-r(i,1)*r(i,2))^2)*(8*r(i,1)*r(i,2)^2-13*r(i,1)^2*r(i,2)-32*r(i,2)^28*r(i,1)^2+25*r(i,1)*r(i,2)+20) 4*r(i,1)/(5*(4-r(i,1)*r(i,2))^2)*(16*r(i,2)^25*r(i,1)*r(i,2)^2+9*r(i,1)*r(i,2)-28*r(i,2)+8*r(i,1)) 4/(5*L(i)*(4-r(i,1)*r(i,2))^2)*(8*r(i,1)*r(i,2)^213*r(i,1)^2*r(i,2)-32*r(i,2)^2-8*r(i,1)^2+25*r(i,1)*r(i,2)+20) r(i,2)/(5*(4r(i,1)*r(i,2))^2)*(16*r(i,1)^2+25*r(i,1)^2*r(i,2)-96*r(i,1)*r(i,2)+128*r(i,2)-28*r(i,1))]; Gi_sr=Gi*Ci_g; %Xac dinh ma tran cung cho phan tu dam-cot nua cung thu i ke den anh huong P-delta if Bac==1 Ki=Ti.'*(Si*Ci_e)*Ti; elseif Bac==2 Ki=Ti.'*(Si*Ci_e+Gi_sr)*Ti; end Se=[Se -Ki(3,:); Ki(6,:)]; 16 %dynamic_com.m - tinh khung phang voi tai dong commandwindow; syms tv; disp('THONG SO VAT LIEU THEP'); prompt={'Ti so can cua thep: ','Khoi luong rieng (kN/m3)', 'He so no hong'}; tisocan_def=0.05; MperV_def=78.255; hesonohong_def=0.3; title='Thong so vat lieu'; lines=1; def={num2str(tisocan_def),num2str(MperV_def),num2str(hesonohong_def),''}; answer=inputdlg(prompt,title,lines,def); j={'j1' 'j2' 'j3'}; if ~isempty(answer) answer=cell2struct(answer,j,1); tisocan=str2num(answer.j1); MperV=str2num(answer.j2); hesonohong=str2num(answer.j3); else break; end disp('NHAP CAC SO LIEU DIEU KIEN BAN DAU') fprintf('\n'); t0_def=0; tf_def=30; Dt_def=0.01; u0_def=0; v0_def=0.1; prompt={'Thoi diem bat dau t0 (s):','Thoi diem ket thuc tf (s):', 'Buoc thoi gian Delta t:','Chuyen vi ban dau u0:','Van toc ban dau v0:'}; title= 'Nhap cac dieu kien ban dau'; lines=1; 136 def={num2str(t0_def),num2str(tf_def),num2str(Dt_def),num2str(u0_def),num2str(v0_def),''}; answer=inputdlg(prompt,title,lines,def,'on'); if ~isempty(answer) j={'j1' 'j2' 'j3' 'j4' 'j5'}; answer=cell2struct(answer,j,1); t0=str2num(answer.j1); tf=str2num(answer.j2); Dt=str2num(answer.j3); u0=str2num(answer.j4); v0=str2num(answer.j5); else return end clear prompt title lines def ans answer j t0_def tf_def Dt_def u0_def v0_def %============================================================= loadtype = questdlg('Phan tich khung voi loai tai trong?','Load Type','Dao dong tu do','Tai dieu hoa','Tai xung','Dao dong tu do'); if strcmp(loadtype,'Dao dong tu do') disp('Phan tich dao dong tu do'); load_type=1; Fn=zeros(1,3*Nnode); elseif strcmp(loadtype,'Tai dieu hoa') load_type=2; disp('Phan tich dao dong voi tai dieu hoa'); loadfunc=input(['Nhap ham tai dieu hoa-Vi du: sin(5*tv) ']); elseif strcmp(loadtype,'Tai xung') pulse_type = questdlg('Chon loai tai xung?','Impulse Load Type','Xung chu nhat','Xung hinh sin','Xung chu nhat'); if strcmp(pulse_type,'Xung chu nhat') load_type=3; disp('Phan tich dao dong voi tai xung chu nhat'); Po_def=20; tp_def=2; prompt={'Gia tri luc xung (kN)','Thoi gian tac dung (s)'}; title=['Thong so tai xung chu nhat']; lines=1; def={num2str(Po_def),num2str(tp_def),''}; answer=inputdlg(prompt,title,lines,def,'on'); if ~isempty(answer) j={'j1' 'j2'}; answer=cell2struct(answer,j,1); Po=str2num(answer.j1); tp=str2num(answer.j2); else return end elseif strcmp(pulse_type,'Xung hinh sin') load_type=4; loadfunc=input(['Nhap ham tai xung-Vi du: P0*sin(5*tv) ']); tp=input('Thoi gian tac dung (s): '); end end %loop=0; %Gan chi so lap ban dau la zero %delta=10^-2; N=zeros(Nel,2); frametype = questdlg('Chon loai khung can phan tich?','Frame Type','Khung cung','Khung nua cung','Khung cung'); if strcmp(frametype,'Khung cung') Lkhung=1; Bac=1; rigid_dynamic; elseif strcmp(frametype,'Khung nua cung') Lkhung=2; Bac=1; escap1=1; r=ones(Nel,2); 137 Mu=zeros(Nel,2); Re=ones(Nel,2)*inf; n=zeros(Nel,2); Rp=zeros(Nel,2); while escap1==1 ques1=questdlg('Tinh lien ket bang mo hinh lien ket nua cung nao?','Semi-connection Model','Tuyen tinh','Eurocode','Richard-Abbott','Eurocode'); if strcmp(ques1,'Tuyen tinh') % Quan he M-theta la tuyen tinh (nut nua cung) Lienket=1; for i=1:Ldam str=num2str(i); rd_def=0.5; prompt={'He so ngam-dau mut (end_fixity factor)'}; title=['He so ngam dau mut loai dam thu ', str]; lines=1; def={num2str(rd_def),''}; answer=inputdlg(prompt,title,lines,def,'on'); if ~isempty(answer) j={'j1'}; answer=cell2struct(answer,j,1); rd(i)=str2num(answer.j1); else return end Red(i)=rd(i)*3*Es*Jd(i)/L_span/(1-rd(i)); for j=1:SLdam(i) Re(Sdamg(i,j),1)=Red(i); Re(Sdamg(i,j),2)=Red(i); end end elseif strcmp(ques1,'Eurocode') Lienket=2; for i=1:Ldam %disp('Connection - Thong so lien ket nua cung') str=num2str(i); prompt={'Do cung ban dau Sj_ini (kN.m/rad)', 'Moment cuc han Mj_Rd (kN.m)'}; title=['Nhap thong so cua lien ket loai dam thu ', str]; lines=1; answer=inputdlg(prompt,title,lines); if ~isempty(answer) j={'j1' 'j2'}; answer=cell2struct(answer,j,1); Red(i)=str2num(answer.j1); Mud(i)=str2num(answer.j2); else return end for j=1:SLdam(i) Re(Sdamg(i,j),1)=Red(i); Re(Sdamg(i,j),2)=Red(i); Mu(Sdamg(i,j),1)=Mud(i); Mu(Sdamg(i,j),2)=Mud(i); end end elseif strcmp(ques1,'Richard-Abbott') Lienket=3; for i=1:Ldam %disp('Connection - Thong so lien ket nua cung') str=num2str(i); prompt={'Thong so hinh dang n cua lien ket','Do cung khoi dau Re cua lien ket', 'Moment tham chieu Mu cua lien ket','Do cung bien dang-cung Rp cua lien ket'}; title=['Nhap thong so cua lien ket loai dam thu ', str]; 138 lines=1; answer=inputdlg(prompt,title,lines); if ~isempty(answer) j={'j1' 'j2' 'j3' 'j4'}; answer=cell2struct(answer,j,1); nd(i)=str2num(answer.j1); Red(i)=str2num(answer.j2); Mud(i)=str2num(answer.j3); Rpd(i)=str2num(answer.j4); else return end for j=1:SLdam(i) n(Sdamg(i,j),1)=nd(i); n(Sdamg(i,j),2)=nd(i); Re(Sdamg(i,j),1)=Red(i); Re(Sdamg(i,j),2)=Red(i); Mu(Sdamg(i,j),1)=Mud(i); Mu(Sdamg(i,j),2)=Mud(i); Rp(Sdamg(i,j),1)=Rpd(i); Rp(Sdamg(i,j),2)=Rpd(i); end end end ques1=questdlg('Co muon sua lai cac thong da nhap khong?','Connection','Khong','Co','Khong'); if strcmp(ques1,'Co') escap1=1; else escap1=0; end end %end vong while chon lien ket for i=1:Nel for j=1:2 r0(i,j)=1/(1+3*Es*J(i)/(Re(i,j)*L(i))); %He so ngam dau mut cua cac phan tu end end semi_dynamic; end hay 17 %rigid_dynamic.m - Tinh khung cung voi tai dong ti=[t0:Dt:tf]; omega=[]; %Gan he so ngam dau mut va cung lien ket cho khung cung r=ones(Nel,2); R=ones(Nel,2)*inf; %Thoi gian tinh close all hidden; timebegin=cputime; disp(' ANALYZING '); %================================================== %Bat dau qua trinh tinh toan stiff1; %Matran tong the khung cung mass; %Matran khoi luong tong the P=zeros(1,3*Nnode); %Veto tai tong the ban dau cua bai toan dong bang 0; boundary; %Nen cac ma tran vibration; %Tinh tan so dao dong rayleigh; %Xac dinh ma tran can Fn=zeros(1,size(PP,2)); if (load_type==2|load_type==4) Fn(3*N_story-2)=subs(loadfunc,tv,ti(1)); elseif load_type==3 Fn(3*N_story-2)=Po; end ut=u0*ones(3*NnodeR,1); %Vecto chuyen vi ban dau vt=v0*ones(3*NnodeR,1); %Vecto van toc ban dau 139 at=inv(MM)*(Fn'-KK*ut-CC*vt); %Vecto gia toc ban dau %================================================== alpha=0.25; delta=0.5; a0=1/(alpha*Dt^2); a1=delta/(alpha*Dt); a2=1/(alpha*Dt); a3=1/(2*alpha)-1; a4=delta/alpha-1; a5=Dt/2*(delta/alpha-2); a6=Dt*(1-delta); a7=delta*Dt; %========================================= Uxtop=[];Uxtop(1)=ut(N_story*3-2); Uxtop_max=Uxtop(1); RAmax=0; nt=(tf-t0)/Dt; if (load_type==3|load_type==4) %tai xung nt_xung=tp/Dt; end %Tinh lap theo thoi gian for it=2:(nt+1) if load_type==2 %tai Dieu hoa Fn(3*N_story-2)=subs(loadfunc,tv,ti(it)); %chi dat len nut dinh khung elseif load_type==3 %xung chu nhat if it*Dt