Phân tích khung nhà công nghiệp thấp tầng (2 tầng hoặc 3 tầng) chịu tải trọng động đất với liên kết chân cột nửa cứng

Đại Học Quốc Gia Tp HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA HOÀNG QUỐC THANH ĐỀ TÀI: PHÂN TÍCH KHUNG NHÀ CÔNG NGHIỆP THẤP TẦNG (2 H0ẶC TẦNG) CHỊU TẢI TRỌNG ĐỘNG ĐẤT VỚI LIÊN KẾT CHÂN CỘT NỬA CỨNG Chuyên ngành : XÂY DỰNG DÂN DỤNG VÀ CÔNG NGHIỆP LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH, tháng 07 năm 2009 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 HỒ CHÍ MINH Cán hướng dẫn khoa học: PGS PHAN NGỌC CHÂU (Ghi rõ họ, tên, học hàm, học vị chữ ký) Cán chấm nhận xét 1: TS NGÔ HỮU CƯỜNG (Ghi rõ họ, tên, học hàm, học vị chữ ký) Cán chấm nhận xét 2: PGS.TS BÙI CÔNG THÀNH (Ghi rõ họ, tên, học hàm, học vị chữ ký) 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 …29…tháng……08…năm 2009 TRƯỜNG ĐẠI HỌC BÁCH KHOA PHÒNG ĐÀO TẠO SĐH CỘNG HOÀ XÃ HỘI CHỦ NGHĨA VIỆT NAM ĐỘC LẬP - TỰ DO - HẠNH PHÚC oOo - oOo -Tp Hồ Chí Minh, Ngày … tháng … năm 2009 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên : HOÀNG QUỐC THANH Ngày, tháng, năm sinh: 07/11/1982 Chuyên ngành : XÂY DỰNG DÂN DỤNG & C.NGHIỆP Phái : Nam Nơi sinh : BÌNH DƯƠNG MSHV : 02105517 I/ ĐỀ TÀI : PHÂN TÍCH KHUNG NHÀ CÔNG NGHIỆP THẤP TẦNG (2 HOẶC TẦNG) CHỊU TẢI TRỌNG ĐỘNG ĐẤT VỚI LIÊN KẾT CHÂN CỘT NỬA CỨNG II/ NHIỆM VỤ VÀ NỘI DUNG : Phân tích động đất ảnh hưởng lên khung nhà công nghiệp thấp tầng với ứng xử phi tuyến liên kết nửa cứng dầm cột chân cột Phân tích sở lý thuyết tính tóan phương pháp tính tóan khung nhà công nghiệp chịu tải trọng động đất Vận dụng lý thuyết tính tóan vào tính kết cấu khung nhà công nghiệp thấp tầng chịu tải động đất Kết luận – Kiến nghị III/ NGÀY GIAO NHIỆM VỤ (Ngày bắt đầu thực LV ghi Quyết định giao đề tài): Ngày / 06 /2007 IV/ NGÀY HOÀN THÀNH NHIỆM VỤ : V/ CÁN BỘ HƯỚNG DẪN : CÁN BỘ HƯỚNG DẪN Ngày 02/07/2009 PGS PHAN NGỌC CHÂU CN BỘ MÔN QL CHUYÊN NGÀNH PGS PHAN NGỌC CHÂU 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 2009 TRƯỞNG PHÒNG ĐT-SĐH TRƯỞNG KHOA QL NGÀNH LỜI CẢM ƠN Em xin chân thành cảm ơn Thầy PGS PHAN NGỌC CHÂU tận tình giúp đỡ đóng góp ý kiến cho em trình thực Luận văn Em xin bày tỏ lòng biết ơn đến Thầy Ths BÙI VĂN CHÚNG góp ý tận tâm vô q giá , tài liệu, kiến thức kinh nghiệm khoa học thực tế thầy giúp đỡ em nhiều trình thực Luận văn Em xin chân thành cảm ơn Ban Giám hiệu, thầy cô phòng Đào tạo Sau đại học, thầy cô Bộ môn Công trình, Khoa kỹ thuật Xây dựng tận tâm truyền đạt cho em kiến thức kỹ nghiên cứu khoa học Xin chân thành cám ơn thầy cô giảng dạy Trường ĐẠI HỌC BÁCH KHOA TP HỒ CHÍ MINH truyền đạt cho em kiến thức qúi báu Tôi xin gởi lời cám ơn đến Ban Giám Hiệu trường Đại học Bình Dương tạo nhiều điều kiện thuận lợi cho thực luận văn, xin gửi lời cảm ơn đến bạn đồng nghiệp tôi, người nhiệt tình ủng hộ thời gian thực Luận văn Cuối xin cám ơn các người thân gia dình, anh em bạn bè, người tìm cách ủng hộ, động viên giúp đỡ trình thực Luận văn Xin chân thành cảm ơn Luận Văn Tôt Nghiệp GVHD: PGS Phan Ngọc Châu MỤC LỤC NHIỆM VỤ LUẬN VĂN THẠC SĨ TRANG LỜI CÁM ƠN MỞ ĐẦU MỤC LỤC CHƯƠNG I : TỔNG QUAN I.1.Các Khái Niệm 1 I.1.1 Sự hình thành động đất I.1.2 Chấn tâm (tâm chấn), chấn tiêu I.1.3 Sóng địa chấn I.1.4 Thang đo cường độ động đất I.1.5 Nguyên nhân phá hoại công trình động đất xảy I.1.6 Liên kết nửa cứng I.2.Tình hình nghiên cứu I.3.Mục tiêu phạm vi đề tài CHƯƠNG II : CÁC MÔ HÌNH LIÊN KẾT NỬA CỨNG 12 II.1.Giới thiệu 12 II.2.ng xử liên kết 12 II.3.Các kiểu cấu tạo liên kết nửa cứng 15 II.4.Các mô hình liên kết nửa cứng 19 II.5.Chọn mô hình liên kết 28 II.6.Chọn kiểu liên kết nửa cứng 30 II.6.1 Liên kết dầm cột thép góc nối cánh cánh (TSA) 31 II.6.2 Liên kết dầm cột thép góc liên kết cánh trên, cánh thép góc nối bụng dầm với cột (TSAW) 33 II.6.3 Liên kết chân cột 35 CHƯƠNG III: MÔ HÌNH PHẦN TỬ LIÊN KẾT NỬA CỨNG 37 III.1 Mô hình phần tử liên kết nửa cứng dầm cột 37 III.2 Hệ số ngàm liên kết 38 III.3 Khảo sát dầm nửa cứng chịu tải trọng phân bố 40 III.4 Thiết lập trường chuyển vị cho phần tử nửa cưng 42 CHƯƠNG IV: PHÂN TÍCH ĐỘNG LỰC HỌC KHUNG NỬA CỨNG VÀ TÍNH TOÁN TẢI TRỌNG ĐỘNG 47 HV: Hoàng Quốc Thanh i Luận Văn Tôt Nghiệp IV.1 IV.2 IV.3 IV.4 GVHD: PGS Phan Ngọc Châu Phân tích động lực học khung nửa cứng 47 IV.1.1 Phương trình cân dao động 47 IV.1.2 Các đặc trưng động lực học khung nửa cứng 50 Tính toán tải trọng động 59 IV.2.1 Tính toán tải trọng động đất 59 IV.2.2 Phương pháp phổ đáp ứng 73 IV.2.3 Phương pháp phân tích theo thời gian lịch sử 82 Giải hệ phương trình động lực học 86 IV.3.1 Các phương pháp giải 86 IV.3.2 Phân tích khung không xét phi tuyến hình học 87 Tác động động đất theo Tiêu chuẩn Việt Nam TCXDVN 375:2006 89 IV.4.1 Tổng quát 89 IV.4.2 Phổ phản ứng đàn hồi theo phương nằm ngang 90 IV.4.3 Phổ thiết kế dùng cho phân tích đàn hồi 93 IV.4.4 Sơ đồ thuật toán phân tích khung có liên kết nửa cứng 94-95 CHƯƠNG V: CHƯƠNG TRÌNH PHÂN TÍCH KHUNG CÓ LIÊN KẾT NỬA CỨNG CHỊU TẢI TRỌNG ĐỘNG 96 V.1.Giới thiệu chương trình máy tính ứng dụng 96 V.2.Giới thiệu hàm chương trình 98 V.2.1 Phần nhập liệu 98 V.2.2 Phần giải toán 101 V.2.3 Phần xuất kết 105 CHƯƠNG VI: THÍ DỤ MINH HOẠ, KẾT LUẬN, KIẾN NGHỊ 106 VI.1 106 Thí dụ minh họa VI.1.1 Tìm tần số góc khung ứng với mode dao động 108 VI.1.2 Phân tích khung chịu tải trọng động đất 122 VI.2 Kết luận, kiến nghị 129 VI.2.1 Kết luận 129 VI.2.2 Kiến nghị 131 PHỤ LUC CHƯƠNG TRÌNH PL1 PL.1 PHẦN NHẬP DỮ LIỆU PL1 PL.2 PHẦN GIẢI BÀI TÓAN PL9 PL.3 PHẦN XUẤT KẾT QUẢ PL58 TÀI LIỆU KHAM THẢO HV: Hoàng Quốc Thanh ii Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu MỞ ĐẦU Các trận động đất liên tiếp năm gần Việt Nam thu hút quan tâm công luận nhà nghiên cứu khoa học Đã có nhiều phát biểu, tranh luận khoa học cần thiết hay không cần thiết phải thiết kế công trình chịu tải trọng động đất, mức độ ảnh hưởng động đất đến công trình xây dựng Trong bối cảnh kinh tế Việt Nam có bước phát triển mạnh mẽ, nói riêng, ngành xây dựng có nhiều bước phát triển vượt bậc Nhiều công trình xây dựng lónh vực Cảng công trình biển, Cầu đường, Thủy lợi đặc biệt ngành Dân dụng Công nghiệp xây dựng phát triển số lượng qui mô tầm cỡ công trình Bên cạnh vấn đề chất lượng công trình đòi hỏi ngày cao hơn, vấn đề quan trọng đưa công trình Việt Nam có bảo đãm khả chịu lực xuất động đất Tại nhiều nước giới, chủ yếu nước thường xảy động đất, vấn đề quan tâm nghiên cứu từ năm đầu kỷ 20 Nhiều nước giới, việc thiết kế công trình chịu tải trọng động đất vấn đề bắt buộc Ở Việt Nam, nhiều hoàn cảnh khách quan, vấn đề động đất chưa có nhiều nghiên cứu ứng dụng rộng rãi Tại trường Đại học Bách khoa thành phố Hồ Chí Minh, chưa có nhiều Luận Văn Cao học nghiên cứu đề tài này, đề tài nghiên cứu mức độ áp dụng công thức phần mềm tính toán nước mà chưa nghiên cứu chất lực động đất tác dụng vào công trình HV: Hoàng Quốc Thanh Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu Qua báo đài, thấy thiệt hại người tài sản trận động đất gây ra, em nhận thức tầm quan trọng mức độ cần thiết công tác thiết kế công trình chịu tải trọng động đất Có thể nói vấn đề đặc biệt quan tâm xã hội nhiều nhà khoa học Với tất lòng ham mê tìm hiểu vấn đề kết cấu công trình, em cố gắng tìm hiểu chất lực động đất tác động lên kết cấu công trình, phương pháp tính toán áp dụng thực tế thiết kế phương pháp sử dụng công tác nghiên cứu khoa học Mục tiêu đưa giải vấn đề thiết kế công trình dân dụng công nghiệp chịu tải trọng động đất sở lý luận khoa học dễ dàng vận dụng phương pháp tính thực tế thiết kế Do đó, Luận Văn bao gồm phần trình bày sở lý luận cần thiết, nhấn mạnh vào công thức ứng dụng thực tế Phần cuối Luận Văn chương trình hàm tính toán kết cấu chịu tải trọng động đất áp dụng thực tế Do tài liệu thu thập chưa phong phú hoàn cảnh thuật ngữ chưa thống nhất, trình độ thân nhiều hạn chế, nên phần trình bày thuật ngữ sử dụng luận văn không tránh khỏi sai sót định Kính mong nhận nhiều ý kiến giúp đỡ thầy cô nhằm giúp em nâng cao hiểu biết Xin chân thành cảm ơn HV: Hoàng Quốc Thanh Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu CHƯƠNG I: TỔNG QUAN I.1 Các Khái niệm:[5, 6] I.1.1 Sự hình thành động đất: Động đất liên quan đến toả khối lượng lớn từ nơi định, nơi nằm sâu lòng đất Có nhiều nguyên nhân dẫn đến phát sinh khối lương gây động đất nguyên nhân chuyển động tương hỗ không ngừng khối vật chất nằm sâu lòng đất để thiết lập cân mới, gọi vận động kiến tạo động đất kết vận động kiến tạo I.1.2 Chấn tâm (tâm chấn), chấn tiêu Điểm phát lượng trân động đất gọi chấn tiêu Điểm chiếu chấn tiêu theo đường thẳng lên mặt đất gọi chấn tâm Hình 1.1 CHẤN TÂM VÀ CHẤN TIÊU I.1.3 Sóng địa chấn Khi động đất xảy ra, lượng giải phóng từ chấn tiêu truyền môi trường xung quanh dạng sóng đàn hồi vật lý: sóng dọc, sóng ngang sóng mặt Tất sóng động đất gây gọi la sóng địa chấn a Sóng dọc (hay sóng sơ cấp P) truyền nhờ thay đổi thể tích môi trường, gây biến dạng kéo nén lòng đất Sóng dọc truyến theo HV: Hoàng Quốc Thanh   1  Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu phương từ chấn tiêu đến trạm quan sát với vận tốc lớn (VP= – km/s), nên chúng máy đo địa chấn ghi nhận sớm b Sóng ngang (sóng thứ cấp S) truyền di theo phương vuông góc với sóng dọc, có vận tốc nhỏ sóng dọc (VS= – Km/s) Sóng ngang không gây thay đổ thể tích mà gây tượng xoắn cắt môi trường (nên gọi sóng cắt) Các sóng dọc sóng ngang tới mặt đất chuyển thành sóng mặt gây chuyển động độ sâu nhỏ Sóng mặt (gần giống sóng biển) gây kéo, nén, cắt mặt đất Vận tốc truyền sóng mặt phụ thuộc vào tính chất lý tần đất phủ Nền đất cứng (vững chắc) vận tốc truyền sóng tăng (2.5 đến 15Km/s), đất yếu, vận tốc truyền sóng giảm (0.5 – 1.5 Km/s) Sóng mặt xuất môi trường đồng gọi sóng Reyleigh Trên sở lý thuyết học môi trường liên tục, đàn hồi lý tưởng, đồng chất đẳng hướng, vật tốc truyền sóng dọc VP sóng ngang VS tính theo công thức: ⎛ ⎞ E (1 − ν ) ⎟⎟ V P = ⎜⎜ ⎝ ρ (1 + ν )(1 − 2ν ) ⎠ ⎛G⎞ VS = ⎜ ⎟ ⎝P⎠ 1/ 1/ ⎛ ⎞ E ⎟⎟ = ⎜⎜ ⎝ ρ (1 + ν ) ⎠ (1.1) 1/ (1.2) Trong đó: E = Modul đàn hồi đất G = Modul chống cắt đất ρ = Mật độ môi trường sóng truyền ν = Hệ số Poisson môi trường truyền sóng HV: Hoàng Quốc Thanh   2  Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu global a0 a1 a2 a3 a4 a5 a6 a7 ut vt at global Xjoint Seb global PnQ Pnq it Connection global Nspan Lspan Hstory Zjoint Reb rb Rec rc Rbas Mub nb Uz Ry Q global NSbeam SBeam X q NScolumn SColumn global Rki Mu n global rluu Rluu global Se time groundtype ag Seb=[]; for j=1:2 W(i,j)=E*J(i)/(L(i)*R(i,j)); end DD(i)=(1+4*W(i,1))*(1+4*W(i,2))-4*W(i,1)*W(i,2); Seb=E*J(i)/(L(i)^2*DD(i))* [6*(1+2*W(i,2)), 4*L(i)*(1+3*W(i,2)), -6*(1+2*W(i,2)), 2*L(i); 6*(1+2*W(i,1)), 2*L(i), -6*(1+2*W(i,1)), 4*L(i)*(1+3*W(i,1))]; function SemiConn %************************************************************************** % SemiConn.m % CHUONG TRINH XAC DINH THONG SO CUA LIEN KET NUA CUNG MO HINH KISHI-CHEN % CHO LOAI LIEN KET NUA CUNG SAU: % 1) LIEN KET BANG THEP GOC NOI BAN CANH TREN VA DUOI CUA DAM VOI COT % (TOP-AND SEAT-ANGLE CONNECTION: TSA) % 2) LIEN KET BANG THEP GOC NOI BAN CANH TREN VA DUOI KET HOP VOI THEP % GOC BAN BUNG (TOP-AND SEAT-ANGLE WITH WEB-ANGLE CONNECTION: TSAW) % Nguoi lap trinh : HOANG QUOC THANH % % Thay huong dan : PGS.PHAN NGOC CHAU % % Ngay hoan : 20/06/2008 % %************************************************************************** global Nelement Njoint C S L E A J Ac Jc Ab Jb Re Ree r0 M_theta Hc Hb tc tb Hcd Bcd tbd tcd global DynaLoad global Ncolumn Nbeam Nstory Nbeamtype Ncolumntype global MperV WperV muy Fy Xuy global Pn LoadFunc global iFrame jFrame global r N global KK MM global K M R Pdelta global qt Qt ax global P Ux Uxx NjointR global Omega T Factor Phi Mnorm global PP t Giot CC PPi global a0 a1 a2 a3 a4 a5 a6 a7 ut vt at global Xjoint Seb global PnQ Pnq it Connection global Nspan Lspan Hstory Zjoint Reb rb Rec rc Rbas Mub nb Uz Ry Q global NSbeam SBeam X q NScolumn SColumn HV: Hoàng Quốc Thanh  PL 50  Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu global Rki Mu n global rluu Rluu global Se time groundtype ag disp(' ***** Don vi la kN-m *****'); disp('*=========================================================*'); disp('* NHAP LOAI LIEN KET: *'); disp('*=========================================================*'); SemiType=input('Lien ket TSA nhap 1, Lien ket TSAW nhap 2, (1/2): '); Fy=248200; %Gioi han chay cua thep goc (kN/m2) E=200E+6; %Modul dan hoi cua thep goc (kN/m2) % Gan kich thuoc ban dau cho thep goc lien ket canh tren Lt=Bcd+0.02; Tt=0.016; Kt=0.024; Gt=0.08; Woc=0.04; D=Hcd; % Gan kich thuoc ban dau cho thep goc lien ket ban bung Da=Hcd-2*tcd-0.02; Ta=0.008; Ka=0.018; Ga=0.06; disp(' Do you want assign Connection parameters or get default? '); AssignCon = input(' Get default = press 0, Assign = press (0/1):'); if AssignCon==1 Lt=input([' Chieu dai thep goc ban canh tren > Lt = ']); Tt=input([' Chieu day thep goc ban canh tren > Tt = ']); Kt=input([' K/c tu got thep goc tren den khop deo dau tien -> Kt=']); Gt=input([' K/cach tu got thep goc den tam lo bulong dau tien ->Gt=']); Woc=input([' Chieu rong cua dai oc > Woc = ']); D=input([' Chieu cao dam (m) > D = ']); %========================================================================= = % NEU LIEN KET TSAW (SemiType=2) THI NHAP TIEP DU LIEU CHO THEP GOC BUNG %========================================================================= = if (SemiType==2) Da=input([' Chieu dai thep goc bung (m) > Da = ']); Ta=input([' Chieu day thep goc bung (m) > Ta = ']); Ka=input([' K/c tu got thep goc bung den k/deo dau tien ->Ka = ']); Ga=input([' K/cach tu got thep goc bung den lo bulong ->Ga = ']); end end % Gan kich thuoc thep goc lien ket canh duoi = thep goc lien ket canh tren Ls=Lt; Ts=Tt; Ks=Kt; Gs=Gt; %========================================================================= = % TINH DO CUNG KHOI DAU *** INITIAL STIFFNESS (Rki) %========================================================================= = Ll=(D-Da)/2; EIt=E*Lt*Tt^3/12; HV: Hoàng Quốc Thanh  PL 51  Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu EIa=E*Da*Ta^3/12; D1=D+Ts/2+Tt/2; D3=D/2+Ts/2; G1=Gt-Woc/2-Tt/2; G3=Ga-Woc/2-Ta/2; Rkit=3*EIt*D1^2/(G1*(G1^2+0.78*Tt^2)); if SemiType==1 Rki=Rkit elseif SemiType==2 Rkia=6*EIa*D3^2/(G3*(G3^2+0.78*Ta^2)); Rki=Rkit+Rkia end %========================================================================= = % TINH MOMENT CUC HAN *** ULTIMATE MOMENT _ (Mu) %========================================================================= = if SemiType==1 G2=Gt-Kt-Woc/2-Tt/2; Mos=Fy*Ls*Ts^2/4; Mpt=Fy*Lt*Tt^2/4; Vot=(Fy*Lt*Tt)/2; Vptn=solve('(Vpt/Vot)^4+G2/Tt*(Vpt/Vot)-1=0','Vpt'); Vptn=eval(Vptn); %Loai bo cac nghiem am va nghiem phuc for ii=1:size(Vptn) if isreal(Vptn(ii))==1 if Vptn(ii)>=0 Vpt=Vptn(ii); end end end D2=D+Ts/2+Kt; Mu=Mos+Mpt+Vpt*D2 elseif SemiType==2 G2=Gt-Kt-Woc/2-Tt/2; Mos=Fy*Ls*Ts^2/4; Mpt=Fy*Lt*Tt^2/4; Vot=(Fy*Lt*Tt)/2; Vptn=solve('(Vpt/Vot)^4+G2/Tt*(Vpt/Vot)-1=0','Vpt'); Vptn=eval(Vptn); %Loai bo cac nghiem am va nghiem phuc for ii=1:size(Vptn) if isreal(Vptn(ii))==1 if Vptn(ii)>=0 Vpt=Vptn(ii); end end end D2=D+Ts/2+Kt; Gy=Ga-Ka; Voa=(Fy*Ta)/2; Vpun=solve('(Vpu/Voa)^4+Gy/Ta*(Vpu/Voa)-1=0','Vpu'); Vpun=eval(Vpun); %Loai bo cac nghiem am va nghiem phuc for ii=1:size(Vpun) HV: Hoàng Quốc Thanh  PL 52  Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu if isreal(Vpun(ii))==1 if Vpun(ii)>=0 Vpu=Vpun(ii); end end end Vpa=(Vpu+Fy*Ta/2)/2*Da; D4=(2*Vpu+Fy*Ta/2)/(3*(Vpu+Fy*Ta/2))*Da+Ll+Ts/2; Mu=Mos+Mpt+Vpt*D2+2*Vpa*D4 end %========================================================================= = % TINH THONG SO HINH DANG *** SHAPE PARAMETER _ (n) %========================================================================= = Theta0=Mu/Rki; if SemiType==1 n=max(2.003*log10(Theta0)+6.070, 0.302) elseif SemiType==2 n=max(1.398*log10(Theta0)+4.631, 0.827) end %========================================================================= = % VE DUONG CONG QUAN HE PHI TUYEN MOMENT-GOC XOAY (M-Theta) %========================================================================= = figure(3) hold off; Theta=0:0.00025:0.05; Mc=Rki.*Theta./(1+(Theta./Theta0).^n).^(1/n); plot(Theta,Mc,'-','LineWidth',2); xlabel('Goc xoay \theta (rad)'); ylabel('Moment M (kN-m)'); title('DUONG CONG QUAN HE MOMENT-GOC XOAY (M-\theta)'); hold off; function StiffSemi %************************************************************************** % StiffSemi.m % % XAC DINH MA TRAN DO CUNG KHUNG KHONG XET P-Delta % % Nguoi lap trinh : HOANG QUOC THANH % % Thay huong dan : PGS.PHAN NGOC CHAU % % Ngay hoan : 20/06/2008 % %************************************************************************** global Nelement Njoint C S L E A J Ac Jc Ab Jb Re Ree r0 M_theta Hc Bc tb tc Hcd Bcd tbd tcd global DynaLoad global Ncolumn Nbeam Nstory Nbeamtype Ncolumntype global MperV WperV muy Fy Xuy global Pn LoadFunc global iFrame jFrame global r N global KK MM global K M R Pdelta HV: Hoàng Quốc Thanh  PL 53  Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu global qt Qt ax global P Ux Uxx NjointR global Omega T Factor Phi Mnorm global PP t Giot CC PPi global a0 a1 a2 a3 a4 a5 a6 a7 ut vt at global Xjoint Seb global PnQ Pnq it Connection global Nspan Lspan Hstory Zjoint Reb rb Rec rc rp Rbas Mub nb Uz Ry Q global NSbeam SBeam X q NScolumn SColumn global Rki Mu n global rluu Rluu global SA time groundtype ag T1 T2 T3 q tB tC tD Sd %========================================================================= % % THIET LAP MA TRAN DO CUNG TONG THE TU MA TRAN DO CUNG PHAN TU % %========================================================================= % K=zeros(3*Njoint,3*Njoint); if Connection==1 for i=1:Nelement % Ma tran chuyen he truc toa Ti=[C(i) S(i) 0 0; -S(i) C(i) 0 0; 0 0 0; 0 C(i) S(i) 0; 0 -S(i) C(i) 0; 0 0 1]; % Xac dinh ma tran cung cua phan tu dam-cot nua cung thu i B11=r(i,1)/(4-r(i,1)*r(i,2)); B12=r(i,1)*r(i,2)/(4-r(i,1)*r(i,2)); B22=r(i,2)/(4-r(i,1)*r(i,2)); Ki_1= E/L(i)*[A(i) 0 -A(i) 0; 12*J(i)/L(i)^2*(B11+B12+B22) 6*J(i)/L(i)*(2*B11+B12) 12*J(i)/L(i)^2*(B11+B12+B22) 6*J(i)/L(i)*(B12+2*B22); 6*J(i)/L(i)*(2*B11+B12) 4*J(i)*3*B11 -6*J(i)/L(i)*(2*B11+B12) 2*J(i)*3*B12; -A(i) 0 A(i) 0; -12*J(i)/L(i)^2*(B11+B12+B22) -6*J(i)/L(i)*(2*B11+B12) 12*J(i)/L(i)^2*(B11+B12+B22) 6*J(i)/L(i)*(B12+2*B22); 6*J(i)/L(i)*(B12+2*B22) 2*J(i)*3*B12 -6*J(i)/L(i)*(B12+2*B22) 4*J(i)*3*B22]; %Xac dinh ma tran cung hinh hoc cua phan tu dam-cot thu i Gi_sr=zeros(6,6); if Pdelta == 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 nua cung cho ma tran cung hinh hoc G32=-(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); G33=(r(i,1)/(5*(4-r(i,1)*r(i,2))^2))*(16*r(i,2)^2+25*r(i,1)*r(i,2)^2-96*r(i,1)*r(i,2)+128*r(i,1)28*r(i,2)); HV: Hoàng Quốc Thanh  PL 54  Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu G35=(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); G36=(4*r(i,2)/(5*(4-r(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)); G62=-(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); G63=(4*r(i,1)/(5*(4-r(i,1)*r(i,2))^2))*(16*r(i,2)^2-5*r(i,1)*r(i,2)^2+9*r(i,1)*r(i,2)-28*r(i,2)+8*r(i,1)); G65=(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); G66=(r(i,2)/(5*(4-r(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)); Ci_g=[0 0 0 0; 0 0; G32 G33 G35 G36 ; 0 0 0; 0 0 0; G62 G63 G65 G66]; Gi_sr=Gi*Ci_g; end % Xac dinh MT cung phan tu dam-cot nua cung thu i co xet P-delta Ki=Ti.'*(Ki_1-Gi_sr)*Ti; % Xac dinh MT cung KET CAU tu MT cung phan tu for h=1:3 for m=1:3 K(3*iFrame(i)+m-3,3*iFrame(i)+h-3)=K(3*iFrame(i)+m-3,3*iFrame(i)+h-3)+Ki(m,h); K(3*iFrame(i)+m-3,3*jFrame(i)+h-3)=K(3*iFrame(i)+m-3,3*jFrame(i)+h-3)+Ki(m,h+3); K(3*jFrame(i)+m-3,3*iFrame(i)+h-3)=K(3*jFrame(i)+m-3,3*iFrame(i)+h-3)+Ki(m+3,h); K(3*jFrame(i)+m-3,3*jFrame(i)+h-3)=K(3*jFrame(i)+m-3,3*jFrame(i)+h-3)+Ki(m+3,h+3); end end end else Connection==2 for i=1:3 % Ma tran chuyen he truc toa Ti=[C(i) S(i) 0 0; -S(i) C(i) 0 0; 0 0 0; 0 C(i) S(i) 0; 0 -S(i) C(i) 0; 0 0 1]; D11=12*r(i,1)/(4-r(i,1)); D22=12/(4-r(i,1)); D12=6*r(i,1)/(4-r(i,1)); Ki_1=[A(i)*E/L(i) 0 -A(i)*E/L(i) 0; (D11+2*D12+D22)*E*J(i)/L(i)^3 (D11+D12)*E*J(i)/L(i)^2 -(D11+2*D12+D22)*E*J(i)/L(i)^3 (D12+D22)*E*J(i)/L(i)^2; (D11+D12)*E*J(i)/L(i)^2 D11*E*J(i)/L(i) -(D11+D12)*E*J(i)/L(i)^2 D12*E*J(i)/L(i); -A(i)*E/L(i) 0 A(i)*E/L(i) 0; -(D11+2*D12+D22)*E*J(i)/L(i)^3 -(D11+D12)*E*J(i)/L(i)^2 (D11+2*D12+D22)*E*J(i)/L(i)^3 -(D12+D22)*E*J(i)/L(i)^2; (D12+D22)*E*J(i)/L(i)^2 D12*E*J(i)/L(i) -(D12+D22)*E*J(i)/L(i)^2 D22*E*J(i)/L(i)]; %Xac dinh ma tran cung hinh hoc cua phan tu cot thu i Gi_sr=zeros(6,6); if Pdelta==2 Gi=N(i)/L(i)*[0 0 0 0; HV: Hoàng Quốc Thanh  PL 55  Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu 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; 000000 -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 nua cung cho ma tran cung hinh hoc P32=-(4/(5*L(i)*(4-r(i,1))^2))*(8*r(i,1)^2-13*r(i,1)-32*r(i,1)^2-8+25*r(i,1)+20); P33=(r(i,1)/(5*(4-r(i,1))^2))*(16+25*r(i,1)-96*r(i,1)+128*r(i,1)-28); P35=(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); P36=(4/(5*(4-r(i,1))^2))*(16*r(i,1)^2-5*r(i,1)^2+9*r(i,1)-28*r(i,1)+8); P62=-(4/(5*L(i)*(4-r(i,1))^2))*(8*r(i,1)-13*r(i,1)^2-32-8*r(i,1)^2+25*r(i,1)+20); P63=(4*r(i,1)/(5*(4-r(i,1))^2))*(16-5*r(i,1)+9*r(i,1)-28+8*r(i,1)); P65=(4/(5*L(i)*(4-r(i,1))^2))*(8*r(i,1)-13*r(i,1)^2-32-8*r(i,1)^2+25*r(i,1)+20); P66=(1/(5*(4-r(i,1))^2))*(16*r(i,1)^2+25*r(i,1)^2-96*r(i,1)+128-28*r(i,1)); Ci_g=[0 0 0 0; 0 0; P32 P33 P35 P36; 0 0 0; 0 0 0; P62 P63 P65 P66]; Gi_sr=Gi*Ci_g; end %Xac dinh ma tran cung ptu cot nua cung thu i co xet P-Delta Ki=Ti.'*(Ki_1-Gi_sr)*Ti; %Xac dinh ma tran cung ket cau tu ma tran cung phan tu for h=1:3 for m=1:3 K(3*iFrame(i)+m-3,3*iFrame(i)+h-3)=K(3*iFrame(i)+m-3,3*iFrame(i)+h-3)+Ki(m,h); K(3*iFrame(i)+m-3,3*jFrame(i)+h-3)=K(3*iFrame(i)+m-3,3*jFrame(i)+h-3)+Ki(m,h+3); K(3*jFrame(i)+m-3,3*iFrame(i)+h-3)=K(3*jFrame(i)+m-3,3*iFrame(i)+h-3)+Ki(m+3,h); K(3*jFrame(i)+m-3,3*jFrame(i)+h-3)=K(3*jFrame(i)+m-3,3*jFrame(i)+h-3)+Ki(m+3,h+3); end end end for i=4:Nelement % Ma tran chuyen he truc toa Ti=[C(i) S(i) 0 0; -S(i) C(i) 0 0; 0 0 0; 0 C(i) S(i) 0; 0 -S(i) C(i) 0; 0 0 1]; % Xac dinh ma tran cung cua phan tu dam-cot nua cung thu i B11=r(i,1)/(4-r(i,1)*r(i,2)); B12=r(i,1)*r(i,2)/(4-r(i,1)*r(i,2)); B22=r(i,2)/(4-r(i,1)*r(i,2)); Ki_1= E/L(i)*[A(i) 0 -A(i) 0; 12*J(i)/L(i)^2*(B11+B12+B22) 6*J(i)/L(i)*(2*B11+B12) 12*J(i)/L(i)^2*(B11+B12+B22) 6*J(i)/L(i)*(B12+2*B22); 6*J(i)/L(i)*(2*B11+B12) 4*J(i)*3*B11 -6*J(i)/L(i)*(2*B11+B12) 2*J(i)*3*B12; -A(i) 0 A(i) 0; HV: Hoàng Quốc Thanh  PL 56  Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu -12*J(i)/L(i)^2*(B11+B12+B22) -6*J(i)/L(i)*(2*B11+B12) 12*J(i)/L(i)^2*(B11+B12+B22) 6*J(i)/L(i)*(B12+2*B22); 6*J(i)/L(i)*(B12+2*B22) 2*J(i)*3*B12 -6*J(i)/L(i)*(B12+2*B22) 4*J(i)*3*B22]; %Xac dinh ma tran cung hinh hoc cua phan tu dam-cot thu i Gi_sr=zeros(6,6); if Pdelta == 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 nua cung cho ma tran cung hinh hoc G32=-(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); G33=(r(i,1)/(5*(4-r(i,1)*r(i,2))^2))*(16*r(i,2)^2+25*r(i,1)*r(i,2)^2-96*r(i,1)*r(i,2)+128*r(i,1)28*r(i,2)); G35=(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); G36=(4*r(i,2)/(5*(4-r(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)); G62=-(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); G63=(4*r(i,1)/(5*(4-r(i,1)*r(i,2))^2))*(16*r(i,2)^2-5*r(i,1)*r(i,2)^2+9*r(i,1)*r(i,2)-28*r(i,2)+8*r(i,1)); G65=(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); G66=(r(i,2)/(5*(4-r(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)); Ci_g=[0 0 0 0; 0 0; G32 G33 G35 G36 ; 0 0 0; 0 0 0; G62 G63 G65 G66]; Gi_sr=Gi*Ci_g; end % Xac dinh MT cung phan tu dam-cot nua cung thu i co xet P-delta Ki=Ti.'*(Ki_1-Gi_sr)*Ti; % Xac dinh MT cung KET CAU tu MT cung phan tu for h=1:3 for m=1:3 K(3*iFrame(i)+m-3,3*iFrame(i)+h-3)=K(3*iFrame(i)+m-3,3*iFrame(i)+h-3)+Ki(m,h); K(3*iFrame(i)+m-3,3*jFrame(i)+h-3)=K(3*iFrame(i)+m-3,3*jFrame(i)+h-3)+Ki(m,h+3); K(3*jFrame(i)+m-3,3*iFrame(i)+h-3)=K(3*jFrame(i)+m-3,3*iFrame(i)+h-3)+Ki(m+3,h); K(3*jFrame(i)+m-3,3*jFrame(i)+h-3)=K(3*jFrame(i)+m-3,3*jFrame(i)+h-3)+Ki(m+3,h+3); end end end end HV: Hoàng Quốc Thanh  PL 57  Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu PL.3 PHẦN XUẤT KẾT QUẢ function Output %************************************************************************** % Output.m % % XUAT DU LIEU NHAP VA KET QUA PHAN TICH KHUNG NUT NUA CUNG % % Nguoi lap trinh : HOANG QUOC THANH % % Thay huong dan : PGS.PHAN NGOC CHAU % % Ngay hoan : 20/06/2008 % %************************************************************************** global Nelement Njoint C S L E A J Ac Jc Ab Jb Re Ree r0 M_theta Hc Hb tc tb Hcd Bcd tbd tcd global DynaLoad global Ncolumn Nbeam Nstory Nbeamtype Ncolumntype global MperV WperV muy Fy Xuy global Pn LoadFunc global iFrame jFrame global r N global KK MM global K M R Pdelta global qt Qt ax global P Ux Uxx NjointR global Omega T Factor Phi Mnorm global PP t Giot CC PPi global a0 a1 a2 a3 a4 a5 a6 a7 ut vt at global Xjoint Seb global PnQ Pnq it Connection global Nspan Lspan Hstory Zjoint Reb rb Rec rc Rbas Mub nb Uz Ry Q global NSbeam SBeam X q NScolumn SColumn global Rki Mu n global rluu Rluu global Se time groundtype ag disp('******************************************************************'); disp('* "Output.m" XUAT DU LIEU NHAP VA KET QUA PHAN TICH *'); disp('******************************************************************'); fprintf('\n'); File_Input=input('filename input (*.m) : ','s'); %Nhap ten file du lieu bai toan (input data file) File_Output=input('filename output (*.txt) : ','s'); %Nhap ten file xuat ket qua (output data file) %In tieu de fprintf( '*********************************************************************\n'); fprintf( '*CHUONGTRINH PHANTICH DONG LUC HOC KHUNGTHEP PHANG LIEN KET NUA CUNG*\n'); fprintf( ' Chuong trinh duoc viet bang ngon ngu MATLAB 7.4 \n'); fprintf( '*********************************************************************\n'); fprintf('\n'); fprintf('Tap tin du lieu bai toan : %s\n',File_Input); HV: Hoàng Quốc Thanh  PL 58  Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Chaâu %In Ten file du lieu bai toan file output (ket qua) fprintf('Tap tin ket qua bai toan : %s\n',File_Output); %In Ten file ket qua bai toan file output (ket qua) fprintf('\n'); if Connection==1 if Pdelta==1 fprintf( '*************************************************************\n'); fprintf( '* PHAN TICH KHUNG LIEN KET CUNG KHONG XET PHI TUYEN HINH HOC*\n'); fprintf( '*************************************************************\n'); elseif Pdelta==2 fprintf( '*************************************************************\n'); fprintf( '* PHAN TICH KHUNG LIEN KET CUNG CO XET PHI TUYEN HINH HOC *\n'); fprintf( '*************************************************************\n'); end elseif Connection==2 if Pdelta==1 fprintf( '*************************************************************\n'); fprintf( '* PHAN TICH KHUNG LKET NUA CUNG KHONG XET PHI TUYEN HINH HOC*\n'); fprintf( '*************************************************************\n'); elseif Pdelta==2 fprintf( '*************************************************************\n'); fprintf( '* PHAN TICH KHUNG LKET NUA CUNG CO XET PHI TUYEN HINH HOC *\n'); fprintf( '*************************************************************\n'); end end %************************************************************************** fprintf('\n'); % 1) DU LIEU HINH HOC KET CAU fprintf('1 DU LIEU HINH HOC KET CAU\n'); fprintf('So tang Nstory = %g\n',Nstory); fprintf('So nhip Nspan = %g\n',Nspan); fprintf('Chieu cao tang Hstory = : %g\n',Hstory); fprintf('Chieu rong nhip Lspan = : %g\n',Lspan); fprintf('Tong so nut : %g\n',Njoint); %In Tong so nut file output.m fprintf('Tong so phan tu : %g\n',Nelement); fprintf('Tong so nut lien ket voi dat : %g\n',Njoint-NjointR); if Connection==2 if M_theta==1 HV: Hoàng Quốc Thanh  PL 59  Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu fprintf('Do cung tuyen tinh cua lien ket, R = %g\n', Reb); fprintf('He so ngam cua lien ket, r = %g\n', rb); else M_theta==2 fprintf('Su dung mo hinh thong so (Kishi-Chen)\n'); fprintf('Do cung khoi dau cua lien ket, Re : %g\n', Reb); fprintf('Moment cuc han cua lien ket, Mu : %g\n', Mub); fprintf('Thong so hinh dang cua lien ket, n: %g\n', nb); end end %************************************************************************** fprintf('\n'); % 2) TOA DO NUT fprintf('2 TOA DO NUT \n'); fprintf('Nut T.Do-X T.Do-Z\n'); for i=1:Njoint fprintf('%3.0f %13.3f %13.3f\n',i,Xjoint(i),Zjoint(i)); %Print coordinates of joint i to file output.m end %************************************************************************** fprintf('\n'); % 3) DAC TRUNG VAT LIEU fprintf('3 DAC TRUNG VAT LIEU\n'); fprintf('Khoi luong rieng cua vat lieu MperV = %g\n',MperV); fprintf('Trong luong rieng cua vat lieu WperV = %g\n',WperV); fprintf('Modul dan hoi cua vat lieu E = %g\n',E); fprintf('He so no hong cua vat lieu muy = %g\n',muy); fprintf('Gioi han chay cua vat lieu Fy = %g\n',Fy); fprintf('Ti so can cua vat lieu Xuy = %g\n',Xuy); %In dac trung vat lieu file output.m %************************************************************************** fprintf('\n'); % 4) DAC TRUNG TIET DIEN PHAN TU fprintf('4 DAC TRUNG TIET DIEN PHAN TU\n'); fprintf('PTu N.Dau N.Cuoi D.Tich J \n'); for i=1:Nelement fprintf('%3.0f %6.0f %6.0f %13.3e %13.3e \n', i,iFrame(i),jFrame(i),A(i),J(i)); %In du lieu phan tu thu i file output.m end %************************************************************************** fprintf('\n'); % 5) DIEU KIEN BIEN fprintf('5 DIEU KIEN BIEN\n'); fprintf('(Khoa = ; Tu = 1)\n'); HV: Hoàng Quốc Thanh  PL 60  Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu fprintf('Nut UX UZ RY\n'); for i=1:Njoint fprintf('%3.0f %5.0f %5.0f %5.0f\n', i,Ux(i),Uz(i),Ry(i)); %In du lieu dieu kien bien nut thu i file output.m end %************************************************************************** fprintf('\n'); % 6) TAI TRONG fprintf('6 TAI TRONG\n'); % 6.1) Tai nut fprintf(' 6.1 Tai nut\n'); fprintf('Nut FX FY FZ\n'); for i=1:Njoint if Pn(3*i-2)~=0 | Pn(3*i-1)~=0 | Pn(3*i)~=0 fprintf('%3.0f %12.3e %11.3e %11.3e\n', i,Pn(3*i-2),Pn(3*i-1),Pn(3*i)); %In du lieu tai nut thu i file output.m end end fprintf('\n'); % 6.2) Tai phan tu fprintf(' 6.2 Tai phan tu\n'); fprintf('PTu q Q\n'); for i=1:Nelement if qt(i)~=0 | Qt(i)~=0 fprintf('%3.0f %12.3e %11.3e\n',i,qt(i),Qt(i)); %In du lieu tai nut thu i file output.m end end %************************************************************************** fprintf('\n'); % 7) KET QUA CHUYEN VI fprintf('7 KET QUA CHUYEN VI NUT\n'); fprintf('Nut Ux Uz Ry\n'); for i=1:Njoint fprintf('%3.0f %13.4e %12.4e %12.4e\n',i,q(3*i-2),q(3*i-1),q(3*i)); %In chuyen vi nut thu i file output.m end %************************************************************************** fprintf('\n'); % 8) KET QUA NOI LUC fprintf('8 KET QUA NOI LUC\n'); fprintf('PTu N Q %end HV: Hoaøng Quốc Thanh  M\n'); PL 61  Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu for i=1:Nelement fprintf('%3.0f\n',i); %In so thu tu phan tu file output.m fprintf(' %3.0f %13.4e %12.4e %12.4e\n',iFrame(i),N(i,1),Q(i,1),M(i,1)); %In ket qua noi luc tai dau I cua phan tu thu i file output.m fprintf(' %3.0f %13.4e %12.4e %12.4e\n',jFrame(i),N(i,2),Q(i,2),M(i,2)); %In ket qua noi luc tai dau J cua phan tu thu i file output.m fprintf('\n'); end %************************************************************************** % 9) KET QUA PHAN LUC NUT fprintf('9 KET QUA PHAN LUC NUT\n'); fprintf('Nut Rx Rz My\n'); for i=1:Nelement if (Ux(iFrame(i))==0) fprintf('%3.0f %13.4e %12.4e %12.4e\n',iFrame(i),Q(i),N(i),M(i)); %In ket qua phan luc cua nut lien ket dat thu i file output.m end end HV: Hoàng Quốc Thanh  PL 62  LÝ LỊCH TRÍCH NGANG LÝ LỊCH BẢN THÂN Họ tên : HOÀNG QUỐC THANH Phái : Nam Ngày, tháng, năm sinh : 07/11/1982 Nơi sinh : Bình Dương Địa liên lạc : 24/24 K7, ấp Chánh Lộc 1, xã Chánh Mỹ, TX TDM, Bình Dương Điện thoại liên lạc : 0904.424.643 QUÁ TRÌNH ĐÀO TẠO ĐẠI HỌC : Hình thức đào tạo : đại học qui Thời gian đào tạo : 09/2000 đến 9/2004 Trường đào tạo : TRƯỜNG ĐẠI HỌC BÌNH DƯƠNG Ngành học : Xây dựng dân dụng công nghiệp CAO HỌC : Thời gian học : từ năm 2005 đến Tên luận văn : PHÂN TÍCH KHUNG NHÀ CÔNG NGHIỆP THẤP TẦNG (2 HOẶC TẦNG) CHỊU TẢI TRỌNG ĐỘNG ĐẤT VỚI LIÊN KẾT CHÂN CỘT NỬA CỨNG Cán hướng dẫn : PGS PHAN NGỌC CHÂU QUÁ TRÌNH CÔNG TÁC 11/2004 đến : Khoa Xây dựng trường Đại học Bình Dương Ngày 02 tháng 07 năm 2009 Người khai, HOÀNG QUỐC THANH Luận Văn Tốt Nghiệp GVHD: PGS Phan Ngọc Châu TÓM TẮT Mặc dù tính chất nửa cứng liên kết dầm – cột giới thiệu nghiên cứu hàng chục năm qua, phân tích thiết kế kết cấu thép cho nhà công nghiệp thường giả thiết liên kết dầm – cột liên kết chân cột với móng cứng tuyệt đối khớp lý tưởng Các giả thiết làm đơn giản hóa trình phân tích thiết kế khung thép, nhiên chúng chưa phản ánh ứng xử thật liên kết, nghiên cứu thực nghiệm chứng tỏ ứng xử liên kết nút khung nằm hai trạng thái cứng khớp, liên kết có ứng xử gọi liên kết nửa cứng Trong luận văn này, tác giả nghiên cứu ảnh hưởng liên kết nửa cứng đến ứng xử động lực học khung nhà công nghiệp thấp tầng (khung thép phẳng) chịu tác dụng tải trọng động đất Trong đó, liên kết mô hình hóa lò xo xoay có chiều dài không liên kết dầm- cột nửa cứng có quan hệ moment – góc xoay liên kết theo mô hình tuyến tính mô hình phi tuyến lũy thừa ba thông số Kishi- Chen (1990), liên kết chân cột liên kết nửa cứng có quan hệ moment –góc xoay theo mô hình tuyến tính nh hưởng độ mềm liên kết trình phân tích khung xem xét thông qua hệ số ngàm liên kết Tính phi tuyến hình học được xem xét đến trình phân tích dạng hiệu ứng P-delta Kích thích động lực học động đất xác định thông qua phương pháp phổ phản ứng thiết kế Thủ tục phân tích bước thời gian áp dụng để xác định phản ứng hệ kết cấu Do đặc tính ứng xử phi tuyến liên kết phi tuyến hình học, sau bước phân tích, ma trận độ cứng, ma trận khối lương tương thích, ma trận cản vector tải trọng nút cập nhật theo hệ số ngàm liên kết Hệ phương trình vi phân cân dao động giải phương pháp tích phân số Newmark HV: Hoàng Quốc Thanh ... văn phân tích ảnh hưởng động đất lên nhà công nghiệp thấp tầng Thực chất phân tích nhà công nghiệp thấp tầng phân tích khung thép phẳng không gian chịu tải trọng động đất có liên kết dầm cột liên. .. DỤNG & C.NGHIỆP Phái : Nam Nơi sinh : BÌNH DƯƠNG MSHV : 02105517 I/ ĐỀ TÀI : PHÂN TÍCH KHUNG NHÀ CÔNG NGHIỆP THẤP TẦNG (2 HOẶC TẦNG) CHỊU TẢI TRỌNG ĐỘNG ĐẤT VỚI LIÊN KẾT CHÂN CỘT NỬA CỨNG II/... kết chân cột 35 CHƯƠNG III: MÔ HÌNH PHẦN TỬ LIÊN KẾT NỬA CỨNG 37 III.1 Mô hình phần tử liên kết nửa cứng dầm cột 37 III.2 Hệ số ngàm liên kết 38 III .3 Khảo sát dầm nửa cứng chịu tải trọng phân bố