Phân tích khung thép phẳng có liên kết nửa cứng bằng phương pháp nâng cao

ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA o0o TRẦN TUẤN KIỆT PHÂN TÍCH KHUNG THÉP PHẲNG CÓ LIÊN KẾT NỬA CỨNG BẰNG PHƯƠNG PHÁP NÂNG CAO Chuyên ngành : Xây dựng dân dụng công nghiệp Mã số ngành: 23.04.10 LUẬN VĂN THẠC SĨ PHẦN THUYẾT MINH TP Hồ Chí Minh, tháng 12 năm 2002 Đại Học Quốc Gia Tp Hồ Chí Minh TRƯỜNG ĐẠI HỌC BÁCH KHOA  CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc Lập - Tự Do - Hạnh Phúc  NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên : TRẦN TUẤN KIỆT Phái : Nam Ngày, tháng, năm sinh : 27 - 05 - 1974 Nơi sinh : Nha Trang - Khánh Hòa Chuyên ngành : Xây Dựng DD-CN Mã số : 23.04.10 TÊN ĐỀ TÀI : Phân tích khung thép phẳng có liên kết nửa cứng pháp phân tích nâng cao phương II NHIỆM VỤ VÀ NỘI DUNG : • Nghiên cứu liên kết nửa cứng khung thép • Nghiên cứu phương pháp phân tích nâng cao dùng để phân tích khung thép phẳng có liên kết nửa cứng • Xây dựng phần mềm để phân tích khung - So sánh kết - So sánh khung cứng (khung có liên kết cứng) khung nửa cứng (khung có liên kết nửa cứng) • Nghiên cứu phân tích khung có liên kết nửa cứng ANSYS So sánh kết III NGÀY GIAO NHIỆM VỤ : IV NGÀY HOÀN THÀNH NHIỆM VỤ : V HỌ VÀ TÊN CÁC BỘ HƯỚNG DẪN : TS BÙI CÔNG THÀNH VI HỌ VÀ TÊN CÁN BỘ NHẬN XÉT : VII HỌ VÀ TÊN CÁN BỘ NHẬN XÉT : CÁN BỘ HƯỚNG DẪN CÁN BỘ NHẬN XÉT CÁN BỘ NHẬN XÉT 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 TRƯỞNG PHÒNG ĐÀO TẠO SĐH tháng năm 2002 CHỦ NHIỆM NGÀNH PGS TS CHU QUỐC THẮNG CÔNG TRÌNH ĐƯC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA TP HỒ CHÍ MINH ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH Người hướng dẫn khoa học : TS BÙI CÔNG THÀNH Phó Trưởng Khoa Kỹ Thuật Xây Dựng Đại học Bách Khoa Tp Hồ Chí Minh Người chấm nhận xét : Người chấm nhận xét : Luận án Cao Học bảo vệ HỘI ĐỒNG BẢO VỆ LUẬN ÁN CAO HỌC TRƯỜNG ĐẠI HỌC BÁCH KHOA TP HỒ CHÍ MINH ngày tháng năm 2002 Có thể tìm hiểu luận án Thư Viện Trường Đại Học Bách Khoa Tp Hồ Chí Minh Lời Cảm Ơn Tôi phải cám ơn trước hết toàn Thầy, Cô Khoa Kỹ thuật Xây dựng trường Đại Học Bách Khoa Tp Hồ Chí Minh, nơi kiến thức nhận thức hình thành suốt năm Đại học năm Cao học Một lòng biết ơn sâu xa xin gởi tới người Thầy trực tiếp hướng dẫn, uốn nắn suy nghó tôi, Ts Bùi Công Thành Với lời khuyên quý báu cách nhận định vấn đề, cách giải vấn đề với kiến thức sâu rộng hướng dẫn tận tình mình, Thầy giúp vượt qua khó khăn để thực luận án Dó nhiên không nhắc tới thân phụ tôi, người Thầy tạo xác lẫn hồn mẹ tôi, người cho học quý giá nhất: biết bước dài đến hiểu biết Lời cảm tạ sau xin dành cho anh chị tôi, bạn bè tôi, đồng nghiệp cổ vũ, khích lệ giúp đỡ suốt trình thực luận án Xin chân thành cảm ơn tất iv TÓM TẮT Trên giới nay, khung thép kết cấu xây dựng thông dụng Để phân tích khung này, theo phương pháp truyền thống, liên kết cột dầm giả thiết ngàm hoàn toàn khớp lý tưởng Tuy nhiên, ứng xử thật lại nằm hai trạng thái Một liên kết gọi liên kết nửa cứng Giải toán khung có liên kết nửa cứng vấn đề thiết cần đặt Vì vậy, luận án tập trung vào hai vấn đề chính: thứ nghiên cứu liên kết nửa cứng thứ hai nghiên cứu phương pháp phân tích khung thép phẳng có liên kết nửa cứng Để giải vấn đề thứ nhất, tác giả nghiên cứu mô hình liên kết, từ xây dựng chương trình máy tính (bằng Matlab) để tính ba thông số liên kết dựa vào mô hình Kishi-Chen, mô hình tương đối đơn giản kết lại phản ánh với thực tế Vấn đề thứ hai giải qua phương pháp phân tích nâng cao (Advanced Analysis Method) Phương pháp phân tích giải đồng thời toán phi tuyến hình học, phi tuyến vật liệu phi tuyến liên kết (liên kết nửa cứng ứng xử theo quy luật phi tuyến) Hàm ổn định dùng để giải toán phi tuyến hình học (ảnh hưởng P-δ P-∆); môđun tiếp tuyến CRC (Column Research Council) với hàm parabol thể chảy dẻo dùng để giải toán phi tuyến vật liệu; mô hình thông số Kishi-Chen dùng để giải toán phi tuyến liên kết Và vậy, việc kiểm tra ổn định cấu kiện riêng lẻ theo quy phạm không cần thiết ổn định phần tử ổn định toàn kết cấu đưa vào xử lý đồng thời xác định độ bền giới hạn kết cấu Đây lý để gọi phương pháp phân tích phương pháp phân tích nâng cao Trên sở lý thuyết đó, tác giả xây dựng chương trình máy tính ngôn ngữ Visual Basic để phân tích khung thép phẳng có liên kết nửa cứng toán ví dụ dùng để kiểm tra tính đắn chương trình toán ví dụ giải hai trường hợp khung cứng (khung có liên kết giả thiết cứng) khung nửa cứng (khung có liên kết nửa cứng) để so sánh khác biệt Phần cuối cùng, tác giả xin giới thiệu cách mô liên kết nửa cứng phần mềm tiếng ANSYS nêu lên số bình luận sử dụng phần mềm để giải toán có liên kết nửa cứng v MỘT SỐ KÝ HIỆU VIẾT TẮT M A , M B , P mômen lực dọc gia tăng đầu phần tử θ A , θ B , e chuyển vị dọc trục góc xoay liên kết gia tăng α hệ số tỷ lệ (scale factor) β hàm kích thước dầm (được coi hệ số dầm) µ thể tính chất vật liệu (được coi hệ số vật liệu) θ = θ r θ : góc xoay tương đối không thứ nguyên θ0 góc xoay tham chiếu (= Mu/Rki) ∆ Chuyển vị ngang ab diện tích bu lông kéo hàng (in2) bf chiều rộng cánh dầm (in) C số khống chế góc nghiêng đường cong Cj số xấp xỉ đường cong từ phân tích hồi quy tuyến tính (Lui, 1985) d chiều cao dầm (in) db đường kính bulông đinh tán EIt, EIa độ cứng uốn tương ứng nhánh thép góc cánh trên, bên bụng gắn với mặt cột fb ứng suất trung bình dầm (ksi) Fby ứng suất chảy dẻo vật liệu làm bu lông (ksi) Fy ứng suất chảy dẻo vật liệu (ksi) g khoảng cách hàng bu lông (gauge distance) (in) g2 khoảng cách khớp dẻo H1 H2 (Hình 2.11) gt khoảng cách từ gót thép góc đến tâm lỗ đinh nhánh thép góc gắn với cột k khoảng cách từ gót thép góc cánh đến đoạn gãy Hình 2.11 li khoảng cách đo từ biên cánh đáy dầm tới biên thép góc bụng dầm ls chiều rộng thép góc cánh lt chiều rộng thép góc cánh m = M M u : mômen liên kết, dạng không thứ nguyên M mômen uốn bậc hai M0 giá trị bắt đầu mômen liên kết vị trí đường cong bắt đầu xấp xỉ với vi đường cong thực nghiệm Mp mômen dẻo Mu mômen cực hạn liên kết n thông số hình dạng đường cong M-θr p khoảng cách cột bu lông (bolt pitch) (in) pf khoảng cách bu lông (in) P lực dọc trục bậc hai Py tải trọng ép chảy dẻo (squash load) Rkf độ cứng củng cố biến dạng liên kết Rki độ cứng liên kết ban đầu Rkp độ cứng liên kết củng cố biến dạng S1, S2 hàm ổn định Sx môđun tiết diện dầm ta chiều dày thép góc cánh tf chiều dày cánh dầm (in) tfc chiều dày cánh cột (in) chiều dày nối đầu dầm (in) ts chiều dày thép góc cánh tt chiều dày thép góc cánh tw chiều dày bụng (in) V0 σylttt/2 V0a giá trị lực cắt biên thép góc bụng V0t lực cắt dẻo nhánh thép góc cánh bỏ qua ảnh hưởng uốn (= σylttt/2) Vpu giá trị lực cắt biên W đường kính lỗ đinh    vii MỤC LỤC TẬP PHẦN THUYẾT MINH CHƯƠNG : TỔNG QUAN 1 LIÊN KẾT NỬA CỨNG PHÂN TÍCH NÂNG CAO 2.1 Ứng xử khung có liên kết nửa cứng 2.2 Các phương pháp thiết kế khung thép truyền thống MỤC TIÊU VÀ PHẠM VI CỦA ĐỀ TÀI CHƯƠNG : ỨNG XỬ VÀ MÔ HÌNH CỦA LIÊN KẾT NỬA CỨNG 2.1 GIỚI THIỆU 2.2 ỨNG XỬ CỦA LIÊN KẾT 2.3 CÁC LOẠI LIÊN KẾT NỬA CỨNG 12 2.4 CÁC MÔ HÌNH LIÊN KẾT 14 2.4.1 Mô hình tuyến tính 15 2.4.2 Mô hình đa thức 15 2.4.3 Mô hình B-spline 16 2.4.4 Mô hình lũy thừa 16 2.4.5 Mô hình số mũ 17 2.4.6 Các mô hình liên kết dựa phương pháp phần tử hữu hạn 19 2.4 MÔ HÌNH LŨY THỪA CỦA KISHI-CHEN 20 2.5.1 Xác định độ cứng liên kết khởi đầu mômen cực hạn 21 2.5.2 Phương trình thực nghiệm xác định thông số hình dạng 27 2.5.3 Chương trình máy tính xác định thông số liên kết 28 2.6 TÓM TẮT CHƯƠNG : PHÂN TÍCH NÂNG CAO 30 31 3.1 GIỚI THIỆU 31 3.2 CÁC PHƯƠNG PHÁP PHÂN TÍCH NÂNG CAO 33 3.2.1 Phương pháp vùng dẻo 35 3.2.2 Phương pháp giả khớp dẻo 36 3.2.3 Phương pháp khớp đàn-dẻo 37 viii 3.2.4 Phương pháp hiệu chỉnh khớp dẻo 37 3.3 TẠI SAO PHẢI PHÂN TÍCH NÂNG CAO ? 38 3.4 TÓM TẮT 39 CHƯƠNG : PHƯƠNG PHÁP KHỚP DẺO HIỆU CHỈNH 4.1 GIỚI THIỆU 40 40 4.2 CÁC NHÂN TỐ THEN CHỐT ẢNH HƯỞNG ĐẾN ỨNG XỬ CỦA KHUNG THÉP 41 4.2.1 Chảy dẻo uốn 41 4.2.2 Sự chảy dẻo kết hợp với ứng suất dư 43 4.2.3 Các ảnh hưởng bậc 44 4.2.4 Sai số hình học 45 4.2.5 Phi tuyến liên kết 45 4.3 PHÂN TÍCH KHỚP DẺO HIỆU CHỈNH BẬC 46 4.3.1 Hàm ổn định xét tới ảnh hưởng bậc 47 4.3.2 Quan hệ chuyển vị - lực gia tăng 49 4.3.3 Độ bền dẻo tiết diện ngang 53 4.3.4 Hiệu chỉnh độ cứng phần tử có khớp dẻo diện 53 4.3.5 Mô hình môđun tiếp tuyến kết hợp với ứng suất dư 55 4.3.6 Mô hình giảm độ cứng hai-mặt kết hợp với uốn 56 4.3.7 Quan hệ lực gia tăng - chuyển vị kể đến liên kết nửa cứng 59 4.3.8 Các phương pháp xét đến sai số hình học 61 4.4 TÓM TẮT 62 CHƯƠNG : CHƯƠNG TRÌNH PHÂN TÍCH KHUNG THÉP PHẲNG CÓ LIÊN KẾT NỬA CỨNG BẰN PHƯƠNG PHÁP NÂNG CAO 63 5.1 GIỚI THIỆU 63 5.2 TỔ CHỨC CHƯƠNG TRÌNH 63 5.3 CÁCH NHẬP LIỆU 65 5.4 CÁCH LƯU DỮ LIỆU 70 5.5 GIẢI BÀI TOÁN 70 5.6 ĐỌC DỮ LIỆU XUẤT RA 70 5.7 TÓM TẮT 72 CHƯƠNG : CÁC VÍ DỤ MINH HỌA 73 6.1 GIỚI THIỆU 73 6.2 BÀI TOÁN VÍ DỤ 73 6.3 BÀI TOÁN VÍ DỤ 82 6.4 BÀI TOÁN VÍ DỤ 86 6.5 BÀI TOÁN VÍ DỤ 89 ix 6.6 BÀI TOÁN VÍ DỤ 93 6.7 BÀI TOÁN VÍ DỤ 97 6.8 BÀI TOÁN VÍ DỤ 101 6.9 SO SÁNH KHUNG CỨNG VÀ KHUNG CÓ LIÊN KẾT NỬA CỨNG 106 6.10 TÓM TẮT 109 CHƯƠNG : MÔ PHỎNG LIÊN KẾT NỬA CỨNG TRONG ANSYS 111 7.1 GIỚI THIỆU 111 7.2 ĐÔI NÉT VỀ ANSYS 111 7.3 VÍ DỤ MINH HỌA CÁCH MÔ PHỎNG LIÊN KẾT NỬA CỨNG 112 7.4 TÓM TẮT 123 CHƯƠNG : KẾT LUẬN VÀ KIẾN NGHỊ 124 8.1 KẾT LUẬN 124 8.2 KIẾN NGHỊ 125 TÀI LIỆU THAM KHẢO 128 TẬP PHỤ LỤC MÃ NGUỒN CÁC CHƯƠNG TRÌNH MÁY TÍNH A CHƯƠNG TRÌNH TÍNH THÔNG SỐ LIÊN KẾT NỬA CỨNG 1 Liên kết kiểu 1 Liên kết kiểu B CHƯƠNG TRÌNH PHÂN TÍCH KHUNG THÉP PHẲNG CÓ LIÊN KẾT NỬA CỨNG I CÀI ĐẶT VÀ SỬ DỤNG CHƯƠNG TRÌNH II MÃ NGUỒN CỦA CHƯƠNG TRÌNH 9 16 formmodau 16 formnhaplieu 16 saveform 60 formketqua 63 formvehinh 64 module 65 openform 91 C CÁCH NHẬP LIỆU CỦA CÁC VÍ DỤ TRONG PHẦN THUYEÁT MINH 94 x IPPY = GoTo L1000 ElseIf (PPY > And IREDUE > 9) Then Print #2, "P/PY cho phan tu"; IDFRM; Print #2, " lon hon 1" Print #2, "Vecto luc nut o trang thai la:" For i = To NDOFS Print #2, RGS(i) Next i MPC = End If ' 'Neu Momen cuc han cua moi noi > momen cuc han cua cau kien ke no ' If (JA = 1) Then F1 = If (JB = 1) Then F2 = If (JA = 2) Then MMPC(1) = JMPC(IA) F1 = Else MMPC(1) = MPC End If If (JB = 2) Then MMPC(2) = JMPC(IB) F2 = Else MMPC(2) = MPC End If ''Neu Momen cuc han cua moi noi < Momen cuc han cua cau kien ke no ''thi thay the MPC boi RCF(IDFRM,IEND) If (IEND 0) Then If (IULTMO = And RCF(IDFRM, IEND) < MPC) Then If (IEND = 1) Then MMPC(1) = Abs(RCF(IDFRM, IEND)) MMPC(2) = MPC ElseIf (IEND = 2) Then MMPC(1) = MPC MMPC(2) = Abs(RCF(IDFRM, IEND)) End If End If End If 'Neu can thiet phai scaling thi RATIO1 = RATIO2 = If (MMPC(1) < 0.00001) Then GoTo L36 If (F1 > (1 + saisoError)) Then PPYT = Abs(RCFTMP(IDFRM, 3)) / PY XM1T = Abs(RCFTMP(IDFRM, 1)) / MP EQN1 = * XM1T / If (PPYT < EQN1) Then F11 = PPYT / + XM1T ElseIf (PPYT > EQN1) Then F11 = PPYT + * XM1T / End If RATIO1 = - (F1 - 1) / (F1 - F11) For i = To JJJH If (IDFRM = JJMEM(i) And JJEND(i) = 1) Then RATIO1 = If (NCODE(3) = 0) Then GoTo L31 If (NCODE(3) = JJHING(i) And ISAME(NCODE(3)) = 2) Then RATIO1=1 L31: Next i End If If (MMPC(2) < 0.00001) Then GoTo L36 83 If (F2 > (1 + saisoError)) Then PPYT = Abs(RCFTMP(IDFRM, 3)) / PY XM2T = Abs(RCFTMP(IDFRM, 2)) / MP EQN2 = * XM2T / If (PPYT < EQN2) Then F22 = PPYT / + XM2T ElseIf (PPYT > EQN2) Then F22 = PPYT + * XM2T / End If RATIO2 = - (F2 - 1) / (F2 - F22) For i = To JJJH If (IDFRM = JJMEM(i) And JJEND(i) = 2) Then RATIO2 = If (NCODE(6) = 0) Then GoTo L32 If (NCODE(6) = JJHING(i) And ISAME(NCODE(6)) = 2) Then RATIO2=1 L32: Next i End If If (Abs(RATIO1 - RATIO2) < 0.00001) Then GoTo L36 If Abs(RATIO1) < Abs(RATIO2) Then RATIO = Abs(RATIO1) Else RATIO = Abs(RATIO2) End If ISCALE = If (RATIO1 < RATIO2) Then Print #2, "Momen tai dau A cua phan tu khung ", IDFRM; Print #2, " lon hon momen deo, vi vay tai gia tang bi giam bang cach"; Print #2, " nhan voi he so ", RATIO Print #2, IDFRM, RATIO ElseIf (RATIO1 > RATIO2) Then Print #2, "Momen tai dau B cua phan tu khung ", IDFRM; Print #2, " lon hon momen deo, vi vay tai gia tang bi giam bang cach"; Print #2, " nhan voi he so ", RATIO Print #2, IDFRM, RATIO End If GoTo L1000 L36: ' 'Kiem tra cac khop deo ' ''Khop deo o ca dau If (Abs(F1 - 1) < saisoError And Abs(F2 - 1) < saisoError Or (N = 2)) Then IFLAGH = ''' RESET RCF va lap DELMPC If (RCF(IDFRM, 1) > 0) Then RCF(IDFRM, 1) = MMPC(1) If (N = 2) Then DELMPC(IDFRM, 1) = MMPC(1) - JMPC(IA) ElseIf (RCF(IDFRM, 1) < 0) Then RCF(IDFRM, 1) = -MMPC(1) If (N = 2) Then DELMPC(IDFRM, 1) = -MMPC(1) + JMPC(IA) End If If (IULTMO = And IEND = 1) Then If (RCF(IDFRM, 2) > 0) Then RCF(IDFRM, 2) = MMPC(2) If (N = 2) Then DELMPC(IDFRM, ElseIf (RCF(IDFRM, 2) < 0) Then RCF(IDFRM, 2) = -MMPC(2) If (N = 2) Then DELMPC(IDFRM, End If If (IULTMO = And IEND = 2) Then ''' UPDATE JMPC cho MMPC If (N = 2) Then JMPC(IA) = MMPC(1) DELMPC(IDFRM, 1) = 2) = MMPC(2) - JMPC(IB) 2) = -MMPC(2) + JMPC(IB) DELMPC(IDFRM, 2) = 84 JMPC(IB) = MMPC(2) End If GoTo L500 End If ''Khop deo chi tai dau A If (Abs(F1 - 1) < saisoError And Abs(F2 - 1) > saisoError Or (N = And NFLAGH = 1)) Then IFLAGH = ''' RESET RCF va hinh DELMPC If (RCF(IDFRM, 1) > 0) Then RCF(IDFRM, 1) = MMPC(1) If (N = And NFLAGH = 1) Then DELMPC(IDFRM, 1) = MMPC(1) JMPC(IA) ElseIf (RCF(IDFRM, 1) < 0) Then RCF(IDFRM, 1) = -MMPC(1) If (N = And NFLAGH = 1) Then DELMPC(IDFRM, 1) = -MMPC(1) + JMPC(IA) End If If (IULTMO = And IEND = 1) Then DELMPC(IDFRM, 1) = ''' RESET JMPC If (N = And NFLAGH = 1) Then JMPC(IA) = MMPC(1) If (N = And NFLAGH = 1) Then DCFI(2) = DCF(2) - DCFTMP(2) If (ISOLVE = 1) Then Call ETSTIF(RCFTMP(IDFRM, 3), EET, E, PY) Call SFUNC(RCFTMP(IDFRM, 3), S1, S2, LO, EET, INERT) XM2TMP = RCFTMP(IDFRM, 2) PP = Abs(RCFTMP(IDFRM, 3)) / PY Call ESTIFF(PP, XM2TMP, MP, PHIB) RCFI(2) = (EET * INERT / LO) * (S1 - S2 ^ / S1) * DCFI(2) * PHIB + PHIB * S2 * DELMPC(IDFRM, 1) / (S1 - (1 - PHIB) * S2 ^ / S1) ElseIf (ISOLVE = 0) Then Call SFUNC(RCFTMP(IDFRM, 3), S1, S2, LO, E, INERT) RCFI(2) = (E * INERT / LO) * (S1 - S2 ^ / S1) * DCFI(2) + S2 * DELMPC(IDFRM, 1) / S1 End If RCF(IDFRM, 2) = RCFTMP(IDFRM, 2) + RCFI(2) End If GoTo L500 End If ''Khop deo chi tai dau B If (Abs(F1 - 1) > saisoError And Abs(F2 - 1) < saisoError Or (N = And NFLAGH = 2)) Then IFLAGH = ''' RESET RCF va hinh DELMPC If (RCF(IDFRM, 2) > 0) Then RCF(IDFRM, 2) = MMPC(2) If (N = And NFLAGH = 2) Then DELMPC(IDFRM, 2) = MMPC(2) JMPC(IB) ElseIf (RCF(IDFRM, 2) < 0) Then RCF(IDFRM, 2) = -MMPC(2) If (N = And NFLAGH = 2) Then DELMPC(IDFRM, 2) = -MMPC(2) + JMPC(IB) End If If (IULTMO = And IEND = 2) Then DELMPC(IDFRM, 2) = 'Print #2, "phan tu ", IDFRM 'Print #2, RCF(IDFRM, 2) ''' RESET JMPC If (N = And NFLAGH = 2) Then JMPC(IB) = MMPC(2) If (N = And NFLAGH = 2) Then DCFI(1) = DCF(1) - DCFTMP(1) If (ISOLVE = 1) Then Call ETSTIF(RCFTMP(IDFRM, 3), EET, E, PY) Call SFUNC(RCFTMP(IDFRM, 3), S1, S2, LO, EET, INERT) XM1TMP = RCFTMP(IDFRM, 1) PP = Abs(RCFTMP(IDFRM, 3)) / PY 85 Call ESTIFF(PP, XM1TMP, MP, PHIA) RCFI(1) = (EET * INERT / LO) * (S1 - S2 ^ / S1) * DCFI(1) * PHIA + PHIA * S2 * DELMPC(IDFRM, 2) / (S1 - (1 - PHIA) * S2 ^ / S1) ElseIf (ISOLVE = 0) Then Call SFUNC(RCFTMP(IDFRM, 3), S1, S2, LO, E, INERT) RCFI(1) = (E * INERT / LO) * (S1 - S2 ^ / S1) * DCFI(1) + S2 * DELMPC(IDFRM, 2) / S1 End If RCF(IDFRM, 1) = RCFTMP(IDFRM, 1) + RCFI(1) End If GoTo L500 End If ''Khong co khop deo o dau nao ca If (Abs(F1 - 1) > saisoError And Abs(F2 - 1) > saisoError Or (N = 0)) Then IFLAGH = GoTo L500 End If L500: Select Case IFLAGH Case 'Khop deo o ca dau For JJ = To JJH If (IDFRM = JMEM(JJ) And JEND(JJ) = 1) Then GoTo L103 If (NCODE(3) = 0) Then GoTo L101 If (NCODE(3) = JHINGE(JJ) And ISAME(NCODE(3)) = 2) Then GoTo L103 L101: Next JJ JJH = JJH + JHINGE(JJH) = NCODE(3) JMEM(JJH) = IDFRM JEND(JJH) = JMPC(JJH) = MMPC(1) L103: For JJ = To JJH ''' Khop deo hinh tai dau B: If (IDFRM = JMEM(JJ) And JEND(JJ) = 2) Then GoTo L106 If (NCODE(6) = 0) Then GoTo L104 If (NCODE(6) = JHINGE(JJ) And ISAME(NCODE(6)) = 2) Then GoTo L106 L104: Next JJ JJH = JJH + JHINGE(JJH) = NCODE(6) JMEM(JJH) = IDFRM JEND(JJH) = JMPC(JJH) = MMPC(2) L106: '' Xac dinh JFLAGH NNN = NNFLAG = For JJ = To JJH If (IDFRM = JMEM(JJ) And JEND(JJ) = 1) Then NNN = NNN + NNFLAG = End If If (IDFRM = JMEM(JJ) And JEND(JJ) = 2) Then NNN = NNN + NNFLAG = End If L107: Next JJ If (NNN = 2) Then JFLAGH(IDFRM) = If (NNN = And NNFLAG = 1) Then JFLAGH(IDFRM) = If (NNN = And NNFLAG = 2) Then JFLAGH(IDFRM) = 86 If (NNN = 0) Then JFLAGH(IDFRM) = GoTo L707 Case 'Khop deo chi tai dau A ''Cap nhat JHINGE, JMEM AND JEND For JJ = To JJH If (IDFRM = JMEM(JJ) And JEND(JJ) = 1) Then GoTo L203 If (NCODE(3) = 0) Then GoTo L201 If (NCODE(3) = JHINGE(JJ) And ISAME(NCODE(3)) = 2) Then GoTo L203 L201: Next JJ JJH = JJH + JHINGE(JJH) JMEM(JJH) = JEND(JJH) = JMPC(JJH) = = NCODE(3) IDFRM MMPC(1) L203: '' Xac dinh JFLAGH NNN = NNFLAG = For JJ = To JJH If (IDFRM = JMEM(JJ) And JEND(JJ) = 1) Then NNN = NNN + NNFLAG = End If L207: Next JJ If (NNN = And NNFLAG = 1) Then JFLAGH(IDFRM) = If (NNN = 0) Then JFLAGH(IDFRM) = GoTo L707 Case 'Khop deo chi tai dau B '' Cap nhat JHINGE, JMEM AND JEND For JJ = To JJH If (IDFRM = JMEM(JJ) And JEND(JJ) = 2) Then GoTo L303 If (NCODE(6) = 0) Then GoTo L301 If (NCODE(6) = JHINGE(JJ) And ISAME(NCODE(6)) = 2) Then GoTo L303 L301: Next JJ JJH = JJH + JHINGE(JJH) JMEM(JJH) = JEND(JJH) = JMPC(JJH) = = NCODE(6) IDFRM MMPC(2) L303: ''Xac dinh JFLAGH NNN = NNFLAG = For JJ = To JJH If (IDFRM = JMEM(JJ) And JEND(JJ) = 2) Then NNN = NNN + NNFLAG = End If L307: Next JJ If (NNN = And NNFLAG = 2) Then JFLAGH(IDFRM) = If (NNN = 0) Then JFLAGH(IDFRM) = GoTo L707 Case 'Khong co khop deo o dau nao ca JFLAGH(IDFRM) = End Select L707: L1000: If (IEND 0) Then THE(IDFRM) = THETAR 87 MOM(IDFRM) = RCF(IDFRM, IEND) End If End Sub Public Sub MULTMV(C, A, B, M, N) 'Chuong trinh tinh tich cua ma tran A(MxN) va vecto cot B(N) For i = To M C(i) = For j = To N C(i) = C(i) + A(i, j) * B(j) Next j Next i End Sub Public Sub ATAN2(GOC, CANHY, CANHX) If CANHX = Then GOC = 3.14159265358979 / ElseIf CANHY >= And CANHX > Then GOC = Atn(CANHY / CANHX) ElseIf CANHY >= And CANHX < Then GOC = 3.14159265358979 + Atn(CANHY / CANHX) ElseIf CANHY 1.5) Then PT = 0.877 * PY / LAMBDA ^ End If ''Xac dinh lai RCT neu dat den hay vuot qua trang thai gioi han If (JJFAIL(IDTRS) = 0) Then '''Luc doc la keo If (RCT(IDTRS) > 0) Then If (RCT(IDTRS) < PY) Then GoTo Thoat ElseIf (RCT(IDTRS) > PY) Then RCT(IDTRS) = PY JFAIL(IDTRS) = End If '''Luc doc la nen ElseIf (RCT(IDTRS) < 0) Then If (Abs(RCT(IDTRS)) < PT) Then GoTo Thoat ElseIf (Abs(RCT(IDTRS)) > PT And Abs(RCT(IDTRS)) < PY) Then RCT(IDTRS) = -PT JFAIL(IDTRS) = ElseIf (Abs(RCT(IDTRS)) > PT And Abs(RCT(IDTRS)) > PY) Then RCT(IDTRS) = -PY JFAIL(IDTRS) = End If End If End If If (JJFAIL(IDTRS) = 1) Then If (RCT(IDTRS) > 0) Then RCT(IDTRS) = PT ElseIf (RCT(IDTRS) < 0) Then RCT(IDTRS) = -PT End If 90 End If If (JJFAIL(IDTRS) = 2) Then If (RCT(IDTRS) > 0) Then RCT(IDTRS) = PY ElseIf (RCT(IDTRS) < 0) Then RCT(IDTRS) = -PY End If End If Thoat: End Sub openform Private Sub cancelbutton_Click() Openform.Hide End Sub Private Sub Dir1_Change() File1.Path = Dir1.Path End Sub Private Sub Drive1_Change() Dir1.Path = Drive1.Drive End Sub Private Sub File1_DblClick() saveroi = 'chi bat dau giai du lieu day du If (Right$(Dir1.Path, 1) = "\") Then Tentaptin$ = Dir1.Path + File1.FileName Else Tentaptin$ = Dir1.Path + "\" + File1.FileName End If Open Tentaptin$ For Input As #1 Input #1, IGEOM Input #1, ILRFD Input #1, NNODE Input #1, NBOUND Input #1, NINCRE Input #1, NCTYPE Input #1, NFTYPE Input #1, NTTYPE Input #1, NUMCNT Input #1, NUMFRM Input #1, NUMTRS formnhaplieu.IGEOMtext.Text = IGEOM formnhaplieu.ILRFDtext.Text = ILRFD formnhaplieu.NNODEtext.Text = NNODE formnhaplieu.NBOUNDtext.Text = NBOUND formnhaplieu.NINCREtext.Text = NINCRE formnhaplieu.NCTYPEtext.Text = NCTYPE formnhaplieu.NFTYPEtext.Text = NFTYPE formnhaplieu.NTTYPEtext.Text = NTTYPE formnhaplieu.NUMCNTtext.Text = NUMCNT formnhaplieu.NUMFRMtext.Text = NUMFRM formnhaplieu.NUMTRStext.Text = NUMTRS ' 'Load du lieu cho grid1 If NCTYPE > Then ReDim daygrid1(1 To NCTYPE, To 3) For i = To NCTYPE For j = To Input #1, daygrid1(i, j) Next j Next i For i = To NCTYPE For j = To formnhaplieu.Grid1.Row = i formnhaplieu.Grid1.Col = j 91 formnhaplieu.Grid1.Text = daygrid1(i, Next j Next i End If ' 'Load du lieu cho grid2 If NFTYPE > Then ReDim daygrid2(1 To NFTYPE, To 6) For i = To NFTYPE For j = To Input #1, daygrid2(i, j) Next j Next i For i = To NFTYPE For j = To formnhaplieu.Grid2.Row = i formnhaplieu.Grid2.Col = j formnhaplieu.Grid2.Text = daygrid2(i, Next j Next i End If ' 'Load du lieu cho grid3 If NTTYPE > Then ReDim daygrid3(1 To NTTYPE, To 5) For i = To NTTYPE For j = To Input #1, daygrid3(i, j) Next j Next i For i = To NTTYPE For j = To formnhaplieu.Grid3.Row = i formnhaplieu.Grid3.Col = j formnhaplieu.Grid3.Text = daygrid3(i, Next j Next i End If ' 'Load du lieu cho grid4 If NUMCNT > Then ReDim daygrid4(1 To NUMCNT, To 3) For i = To NUMCNT For j = To Input #1, daygrid4(i, j) Next j Next i For i = To NUMCNT For j = To formnhaplieu.Grid4.Row = i formnhaplieu.Grid4.Col = j formnhaplieu.Grid4.Text = daygrid4(i, Next j Next i End If ' 'Load du lieu cho grid5 If NUMFRM > Then ReDim daygrid5(1 To NUMFRM, To 5) For i = To NUMFRM For j = To Input #1, daygrid5(i, j) Next j Next i For i = To NUMFRM For j = To formnhaplieu.Grid5.Row = i formnhaplieu.Grid5.Col = j j) j) j) j) 92 formnhaplieu.Grid5.Text = daygrid5(i, Next j Next i End If ' 'Load du lieu cho grid6 If NUMTRS > Then ReDim daygrid6(1 To NUMTRS, To 5) For i = To NUMTRS For j = To Input #1, daygrid6(i, j) Next j Next i For i = To NUMTRS For j = To formnhaplieu.Grid6.Row = i formnhaplieu.Grid6.Col = j formnhaplieu.Grid6.Text = daygrid6(i, Next j Next i End If ' 'Load du lieu cho grid7 If NNODE > Then ReDim daygrid7(1 To NNODE, To 3) For i = To NNODE For j = To Input #1, daygrid7(i, j) Next j Next i For i = To NNODE For j = To formnhaplieu.Grid7.Row = i formnhaplieu.Grid7.Col = j formnhaplieu.Grid7.Text = daygrid7(i, Next j Next i End If ' 'Load du lieu cho grid8 If NBOUND > Then ReDim daygrid8(1 To NBOUND, To 4) For i = To NBOUND For j = To Input #1, daygrid8(i, j) Next j Next i For i = To NBOUND For j = To formnhaplieu.Grid8.Row = i formnhaplieu.Grid8.Col = j formnhaplieu.Grid8.Text = daygrid8(i, Next j Next i End If ' Close #1 Openform.Hide End Sub j) j) j) j) Private Sub OKbutton_Click() File1_DblClick End Sub C CÁCH NHẬP LIỆU CỦA CÁC VÍ DỤ TRONG PHẦN THUYẾT MINH Ví dụ : 93 Nhập liệu cho toán sau : • Nhập IGEOM, ILRFD, tổng số nút, số lần gia tăng tải • Nhập điều kiện biên : • Nhập số loại phần tử khung : • Nhập phần tử khung : • Nhập tải trọng : • Các ô lại, tất nhập hết Bấm nút Giải toán giải Ví dụ : - Số liệu nhập sau (lưu ý nhập điều kiện biên, giải thuật yêu cầu, tất nút phải khóa thành phần chuyển vị thứ 3): 94 Ví dụ : - Cách nhập liệu vào chương trình : Ví dụ : 95 Ví dụ : Ví dụ : 96 Ví dụ : 97 ... LIÊN KẾT NỬA CỨNG CHƯƠNG : PHƯƠNG PHÁP PHÂN TÍCH NÂNG CAO CHƯƠNG : PHÂN TÍCH KHUNG THÉP PHẲNG CÓ LIÊN KẾT NỬA CỨNG BẰNG PHƯƠNG PHÁP NÂNG CAO CHƯƠNG : CHƯƠNG TRÌNH PHÂN TÍCH KHUNG THÉP PHẲNG CÓ... "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" Đề tài tập trung vào vấn đề: 1) Liên kết nửa cứng, 2) Dùng phương pháp phân nâng cao để phân tích khung thép phẳng có liên. .. để phân tích khung thép phẳng có liên kết nửa cứng • Xây dựng phần mềm để phân tích khung - So sánh kết - So sánh khung cứng (khung có liên kết cứng) khung nửa cứng (khung có liên kết nửa cứng)