BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ XÂY DỰNG TRƯỜNG ĐẠI HỌC KIẾN TRÚC HÀ NỘI - ĐINH TRỌNG QUANG TÍNH TOÁN KẾT CẤU KHUNG THÉP NHÀ CÔNG NGHIỆP BẰNG PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN LUẬN VĂN THẠC SỸ KỸ THUẬT XÂY DỰNG CÔNG TRÌNH DÂN DỤNG & CÔNG NGHIỆP Hà Nội - 2016 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ XÂY DỰNG TRƯỜNG ĐẠI HỌC KIẾN TRÚC HÀ NỘI - ĐINH TRỌNG QUANG KHÓA: 2014 – 2016 TÍNH TOÁN KẾT CẤU KHUNG THÉP NHÀ CÔNG NGHIỆP BẰNG PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN Chuyên ngành: Kỹ thuật Xây dựng Công trình DD&CN Mã số: 60.58.02.08 LUẬN VĂN THẠC SỸ KỸ THUẬT XÂY DỰNG CÔNG TRÌNH DÂN DỤNG & CÔNG NGHIỆP NGƯỜI HƯỚNG DẪN KHOA HỌC: TS PHẠM VĂN TRUNG TS LÊ HỮU THANH Hà Nội - 2016 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ XÂY DỰNG TRƯỜNG ĐẠI HỌC KIẾN TRÚC HÀ NỘI ĐINH TRỌNG QUANG KHÓA: 2014 – 2016 TÍNH TOÁN KẾT CẤU KHUNG THÉP NHÀ CÔNG NGHIỆP BẰNG PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN Chuyên ngành: Kỹ thuật Xây dựng Công trình DD&CN Mã số: 60.58.02.08 LUẬN VĂN THẠC SỸ KỸ THUẬT XÂY DỰNG CÔNG TRÌNH DÂN DỤNG & CÔNG NGHIỆP NGƯỜI HƯỚNG DẪN KHOA HỌC: TS PHẠM VĂN TRUNG TS LÊ HỮU THANH XÁC NHẬN CỦA CHỦ TỊCH HỘI ĐỒNG CHẤM LUẬN VĂN Hà Nội – 2016 LỜI CÁM ƠN Trước hết xin chân thành gửi lời cảm ơn sâu sắc giúp đỡ dẫn hữu ích thường xuyên thầy cô giáo Trường Đại học Kiến Trúc Hà Nội, giúp đỡ tạo điều kiện Khoa Đào tạo sau đại học trình học tập hoàn thành luận văn Tôi xin bày tỏ lòng biết ơn sâu sắc tới TS Phạm Văn Trung TS Lê Hữu Thanh bảo tận tình cho nhiều góp ý đắn tạo điều kiện thuận lợi, cung cấp tài liệu động viên để hoàn thiện luận văn cách hoàn chỉnh tốt Tôi xin trân thành cảm ơn thầy, cô Hội đồng chấm luận văn tham gia đóng góp ý kiến quý báu cho luận văn Tôi xin chân thành cảm ơn bạn bè, đồng nghiệp gia đình động viên, khích lệ, tạo điều kiện giúp đỡ suốt trình thực hoàn thành luận văn Tôi xin chân thành cảm ơn! Hà Nội, ngày 10 tháng năm 2016 Học viên Đinh Trọng Quang LỜI CAM ĐOAN Tôi xin cam đoan Luận văn thạc sỹ công trình nghiên cứu khoa học độc lập Các số liệu khoa học, kết nghiên cứu Luận văn trung thực có nguồn gốc rõ ràng TÁC GIẢ LUẬN VĂN Đinh Trọng Quang MỤC LỤC Lời cảm ơn Lời cam đoan Mục lục Danh mục bảng, biểu Danh mục hình vẽ, đồ thị MỞ ĐẦU  Tính cấp thiết đề tài  Mục đích nghiên cứu đề tài  Phạm vi nghiên cứu  Phương pháp nghiên cứu  Ý nghĩa khoa học thực tiễn đề tài  Cấu trúc luận văn đề tài NỘI DUNG CHƯƠNG 1: KẾT CẤU KHUNG NHÀ CÔNG NGHIỆP 1.1 Tổng quan khung nhà công nghiệp 1.1.1 Khái niệm 1.1.2 Phân loại 1.1.3 Chi tiết cấu tạo phận khung nhà công nghiệp 11 1.2 Tải trọng tác dụng 14 1.2.1 Tĩnh tải (TT) 14 1.2.2 Hoạt tải (HT) 15 1.2.3 Tải trọng gió 15 1.2.4 Các loại tải trọng khác 16 1.2.5 Tổ hợp tải trọng 16 1.3 Tính toán nội lực khung nhà công nghiệp 17 1.3.1 Tính toán nội lực khung phương pháp giải tích 17 1.4.1 Tính toán nội lực khung phương pháp số 21 CHƯƠNG 2: TÍNH TOÁN KẾT CẤU KHUNG NHÀ CÔNG NGHIỆP BẰNG PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN 24 2.1 Phương pháp PTHH đặc điểm tính toán 24 2.1.1 Các khái niệm 24 2.1.2 Các phương trình PTHH 24 2.1.3 Kết nối phần tử hệ kết cấu 28 2.1.4 Áp đặt điều kiện biên 29 2.2 Tính toán nội lực khung nhà công nghiệp PP PTHH 30 2.2.1 Phần tử có tiết diện không thay đổi 30 2.2.2 Phần tử có tiết diện thay đổi theo chiều dài 35 2.3 Chương trình tính toán kết cấu khung thép nhà công nghiệp 40 2.3.1 Khái niệm chung 40 2.3.2 Ngôn ngữ lập trình Fortran 42 2.3.3 Lập trình giao diện GID 44 2.3.4 Sơ đồ khối chương trình 44 2.4 Thiết lập chương trình tính 46 2.4.1 Xây dựng giao diện đầu vào chương trình GID 46 2.4.2 Thiết lập chương trình tính ngôn ngữ Fortan 52 CHƯƠNG 3: VÍ DỤ TÍNH TOÁN 63 3.1 Khái niệm chung 63 3.2 Ví dụ tính toán 63 3.2.1 Số liệu tính toán 63 3.2.2 Trình tự tính toán 65 3.3 Ví dụ tính toán 72 3.3.1 Khung có tiết diện không thay đổi 73 3.3.2 Khung có tiết diện thay đổi 75 3.4 Phân tích đánh giá kết tính toán 77 KẾT LUẬN VÀ KIẾN NGHỊ 79 TÀI LIỆU THAM KHẢO PHỤ LỤC DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT STT Chữ viết tắt Tên đầy đủ DOF Degree of freedom (Bậc tự do) Dir Direction (Phương tác dụng) Fortran GUI GID Phần mềm xử lý dao diện đồ họa HT Hoạt tải PTHH Phần tử hữu hạn TT Tĩnh tải THTT Tổ hợp tải trọng Formula Translating System – Ngôn ngữ lập trình sử dụng cho phương pháp số Graphical User Interface – Giao diện người dùng đồ họa DANH MỤC CÁC BẢNG, BIỂU Số liệu bảng, biểu Tên bảng, biểu Trang Bảng 3.1 Đặc trưng hình học mặt cắt ngang khung 63 Bảng 3.2 Thông số đầu vào vật liệu 63 Bảng 3.3 So sánh kết chuyển vị So sánh kết nội lực 69 70 Bảng 3.4 Bảng 3.5 So sánh phản lực liên kết chân cột 70 Bảng 3.6 Đặc trưng hình học mặt cắt ngang kèo cột 71 Bảng 3.7 Thông số đầu vào vật liệu 72 Bảng 3.8 Đặc trưng hình học mặt cắt ngang kèo cột 72 Bảng 3.9 So sánh kết chuyển vị đỉnh kèo 73 Bảng 3.10 So sánh kết nội lực chân vai cột 74 Bảng 3.11 So sánh phản lực liên kết chân cột 74 Bảng 3.12 So sánh kết chuyển vị đỉnh kèo 75 Bảng 3.13 So sánh kết nội lực chân vai cột 75 Bảng 3.14 So sánh phản lực liên kết chân cột 76 READ(5,1000)(DL(I),I=1,8) IF (DL(8).NE.0.0D0) READ(5,1000)DL(9) WRITE(6,2000)(DL(I),I=1,8) IF (DL(8).NE.0.0D0 )WRITE(6,2001)DL(9) RETURN C COMPUTE ELEMENT STIFFNESS MATRIX [K] CALL STIF06 (S,DL,ST,SA,VV,EL,NST,ISW) RETURN C COMPUTE AND PRINT ELEMENT FORCES, i.e {R}=[k][V]{U} C DETERMINE [k][V] CALL STIF06 (S,DL,ST,SA,VV,EL,NST,ISW) DO 402 I=1,12 402 UT(I) = 0.0D0 DO 405 I=1,12 DO 403 J=1,6 403 UT(I) = UT(I) + SA(I,J)*UL(J,1) DO 404 J=7,12 404 UT(I) = UT(I) + SA(I,J)*UL(J-6,2) 405 CONTINUE C MIDSPAN MOMENTS (+VE; HOGGING MOMENT, -VE; SAGGING MOMENT) FLEX2 = (UT(11) - UT(5))*0.5D0 FLEX3 =-(UT(12) - UT(6))*0.5D0 C DUE TO LOCAL DISTRIBUTED FORCES IF (DL(8).NE.0.0D0) THEN XM = DL(8)*EL*EL/12.0D0 XV = DL(8)*EL/2.0D0 XN = XV IN = IDINT(DL(9)) GOTO (408,406,407),IN C LOCAL FORCES IN 2-DIRECTION 406 UT(2) = UT(2) - XV UT(6) = UT(6) - XM UT(8) = UT(8) - XV UT(12) = UT(12) + XM FLEX3 =-(UT(12) - UT(6))*0.5D0 + DL(8)*EL*EL/8.0D0 GOTO 409 C LOCAL FORCES IN 3-DIRECTION 407 UT(3) = UT(3) - XV UT(5) = UT(5) + XM UT(9) = UT(9) - XV UT(11) = UT(11) - XM FLEX2 = (UT(11) - UT(5))*0.5D0 + DL(8)*EL*EL/8.0D0 GOTO 409 C LOCAL FORCES IN 1-DIRECTION (AXIAL DIRECTION) !!! CHECKING LATER 408 UT(1) = UT(1) - XN UT(7) = UT(7) + XN GOTO 409 ENDIF C PRINT FORCES 409 MCT = MCT - IF (MCT.GT.0) GOTO 410 WRITE(6,3000) MCT = 50 410 WRITE(6,4000)N,MA,IX(1),(UT(I),I=1,6),FLEX2,FLEX3,IX(2), (UT(I),I=7,12) IF (PSOFT.NE.0.0D0) WRITE(7,4001)IX(1),(UT(I),I=1,6), 3IX(2),(UT(I),I=7,12) RETURN C COMPUTE CONSISTENT AND LUMPED MASS MATRICES CALL MASS06 (NST,DL,VV,EL,ST,SA,S,P) RETURN C COMPUTE INTERNAL FORCE VECTOR (RHS) IF (DL(8).NE.0.0D0) THEN C DUE TO LOCAL DISTRIBUTED FORCES XM = DL(8)*EL*EL/12.0D0 XV = DL(8)*EL/2.0D0 XN = XV IN = IDINT(DL(9)) GOTO (611,612,613),IN C LOCAL FORCES IN 1-DIRECTION (AXIAL DIRECTION) 611 GOTO 616 C LOCAL FORCES IN 2-DIRECTION 612 DO 605 I=1,3 P(I) = XV*VV(I,2) P(I+3) = XM*VV(I,3) P(I+6) = P(I) 605 P(I+9) =-P(I+3) GOTO 616 C LOCAL FORCES IN 3-DIRECTION 613 DO 615 I=1,3 P(I) = XV*VV(I,3) P(I+3) =-XM*VV(I,2) P(I+6) = P(I) 615 P(I+9) =-P(I+3) GOTO 616 ENDIF 616 CONTINUE C DUE TO SUPPORT MOVEMENTS C CHECKING SUPPORT MOVEMENT NENS = NEN + NENE = NEN + NEL DO 601 J=1,NDF DO 601 I=NENS,NENE IF (UL(J,I).NE.0.0D0) GOTO 66 601 CONTINUE RETURN C FORM -[K]*{Usupport} 66 CALL STIF06 (S,DL,ST,SA,VV,EL,NST,ISW) CALL EINTF (P,S,UL,NST,NDF,NEL,NENS,NENE) IF (ISW.EQ.9) GOTO RETURN C COMPUTE GEOMENTRIC STIFFNESS MATRIX CALL GEOM06 (S,DL,UL,VV,ST,SA,EL,NST,NDF) RETURN C RETURN C RETURN 1000 FORMAT(8F10.0) 2000 FORMAT(10X,'AXIAL AREA OF MEMBER =',E14.6/ * 10X,'TORSIONAL MOMENT OF INERTIA =',E14.6/ * 10X,'MOMENT OF INERTIA ABOUT 2-AXIS =',E14.6/ * 10X,'MOMENT OF INERTIA ABOUT 3-AXIS =',E14.6/ * 10X,'MODULUS OF ELASTICITY * 10X,'SHEAR MODULUS * 10X,'DENSITY * 10X,'UNIFORM LOAD =',E14.6/ =',E14.6/ =',E14.6/ =',E14.6) 2001 FORMAT(10X,'DIRECTION OF UNIFORM LOAD (LOCAL)=',F14.0/) 3000 FORMAT(//5X, *'3D BEAM ELEMENT (RESULTANTS + IN + LOCAL DIRECTIONS)'// *' ELM MAT NODE AXIAL 2-SHEAR 3-SHEAR TORSION',2X, *'2-MOMENT 3-MOMENT'/) 4000 FORMAT(2I5,1X,'NODE-',I1,1X,6E10.4/11X,'MIDDEL', * 41X,2E10.4/11X,'NODE-',I1,1X,6E10.4) 4001 FORMAT(1X,I5,1X,6E10.4/ * 1X,I5,1X,6E10.4) END C ============================================================ SUBROUTINE UNIVEC (X0,I,J,C,XL,NDM) C C COMPUTE UNIT VECTOR ALONG NODE I,J BY SUNIL C IMPLICIT REAL*8 (A-H,O-Z) DIMENSION C(1),X0(NDM,1) X = X0(1,J) - X0(1,I) Y = X0(2,J) - X0(2,I) Z = X0(3,J) - X0(3,I) XL = DSQRT(X*X+Y*Y+Z*Z) C(1) = X/XL C(2) = Y/XL C(3) = Z/XL C -RETURN END C ====================================================================== SUBROUTINE CROSPR (A,B,C) C C COMPUTE UNIT VECTOR NORMAL TO PLANE A-B C (FROM CROSS PRODUCT OF VECTORS A AND B) BY SUNIL C IMPLICIT REAL*8 (A-H,O-Z) DIMENSION A(3),B(3),C(3) X = A(2)*B(3) - A(3)*B(2) Y = A(3)*B(1) - A(1)*B(3) Z = A(1)*B(2) - A(2)*B(1) XL = DSQRT(X*X+Y*Y+Z*Z) C(1) = X/XL C(2) = Y/XL C(3) = Z/XL C -RETURN END C ====================================================================== SUBROUTINE STIF06 (S,DL,ST,SA,VV,EL,NST,ISW) C C FORMULATION OF ELEMENT STIFFNESS MATRIX C IMPLICIT REAL*8 (A-H,O-Z) DIMENSION S(NST,1),DL(1) DIMENSION VV(3,3),ST(6,6),SA(12,12) C -DO 300 I=1,12 DO 300 J=1,12 300 SA(I,J) = 0.0D0 DO 301 I=1,6 DO 301 J=1,6 301 ST(I,J) = 0.0D0 C FORM LOCAL ELEMENT STIFFNESS MATRIX [k] DO 313 IS = 1,4 GOTO (302,303,304,305),IS C SUNIL 302 ST(1,1) = DL(1)*DL(5)/EL ST(4,4) = DL(2)*DL(6)/EL TEMP = 4.0D0*DL(5)/EL ST(5,5) = TEMP*DL(3) ST(6,6) = TEMP*DL(4) TEMP = 3.0D0/EL/EL ST(3,3) = TEMP*ST(5,5) ST(2,2) = TEMP*ST(6,6) TEMP = 1.5D0/EL ST(2,6) = TEMP*ST(6,6) ST(6,2) = ST(2,6) ST(3,5) =-TEMP*ST(5,5) ST(5,3) = ST(3,5) II = IJ = GOTO 309 303 DO 306 I=1,4 306 ST(I,I) =-ST(I,I) ST(5,5) = ST(5,5)/2.0D0 ST(6,6) = ST(6,6)/2.0D0 ST(5,3) =-ST(5,3) ST(6,2) =-ST(6,2) II = IJ = GOTO 309 304 ST(2,6) =-ST(2,6) ST(3,5) =-ST(3,5) ST(5,3) =-ST(5,3) ST(6,2) =-ST(6,2) II = IJ = GOTO 309 305 DO 308 I=1,4 308 ST(I,I) =-ST(I,I) ST(5,5) = ST(5,5)*2.0D0 ST(6,6) = ST(6,6)*2.0D0 ST(5,3) =-ST(5,3) ST(6,2) =-ST(6,2) II = IJ = C COMPUTE [k][V]; ([k]-LOCAL STIFF, [V]-COOR TRANSFORMATION MATRIX) 309 IA = DO 312 IC=1,2 IB = DO 311 ID=1,2 DO 310 I=1,3 DO 310 J=1,3 DO 310 K=1,3 310 SA(I+IA+II,J+IB+IJ) = SA(I+IA+II,J+IB+IJ) + +ST(I+IA,K+IB)*VV(J,K) 311 312 IB = IB + IA = IA + 313 CONTINUE IF (ISW.EQ.4 OR ISW.EQ.7) RETURN C COMPUTE [K]=[V]T[k][V] IA = DO 316 IC=1,4 IB = DO 315 ID=1,4 DO 314 I=1,3 DO 314 J=1,3 DO 314 K=1,3 314 S(I+IA,J+IB) = S(I+IA,J+IB)+VV(I,K)*SA(K+IA,J+IB) 315 IB = IB + 316 IA = IA + C -RETURN END C ====================================================================== SUBROUTINE MASS06 (NST,D,V,XL,ST,SA,S,P) C C FORMED CONSISTENT AND LUMPED MASS MATRICES BY SUNIL 00/11/14 C IMPLICIT REAL*8 (A-H,O-Z) DIMENSION S(NST,1),V(3,3),ST(6,6),SA(12,12),D(1),P(1) C INITIALIZE ARRAYS DO 300 I=1,12 DO 300 J=1,12 SA(I,J) = 0.0D0 300 CONTINUE DO 301 I=1,6 DO 301 J=1,6 ST(I,J) = 0.0D0 301 CONTINUE C FORM LOCAL CONSISTENT MASS MATRIX [m] DO 313 IS=1,4 GOTO (1,2,3,4),IS TMASS = D(7)*D(1)*XL/420.0D0 ST(1,1) = TMASS*140.0D0 ST(2,2) = TMASS*156.0D0 ST(3,3) = ST(2,2) ST(4,4) = TMASS*140.0D0*D(2)/D(1) ST(5,5) = TMASS*4.0D0*XL*XL ST(6,6) = ST(5,5) ST(2,6) = TMASS*22.0D0*XL ST(6,2) = ST(2,6) ST(3,5) =-ST(2,6) ST(5,3) =-ST(2,6) II = IJ = GOTO 309 ST(2,6) =-ST(2,6) ST(6,2) =-ST(6,2) ST(3,5) =-ST(3,5) ST(5,3) =-ST(5,3) II = IJ = GOTO 309 ST(1,1) = ST(1,1)/2.0D0 ST(2,2) = TMASS*54.0D0 ST(3,3) = ST(2,2) ST(4,4) = ST(4,4)/2.0D0 ST(5,5) =-TMASS*3.0D0*XL*XL ST(6,6) = ST(5,5) ST(2,6) = TMASS*13.0D0*XL ST(6,2) =-ST(2,6) ST(3,5) =-ST(2,6) ST(5,3) = ST(2,6) II = IJ = GOTO 309 ST(2,6) =-ST(2,6) ST(6,2) =-ST(6,2) ST(3,5) =-ST(3,5) ST(5,3) =-ST(5,3) II = IJ = C COMPUTE [m][V]; ([m]-LOCAL MASS, [V]-COOR TRANSFORM MATRIX) 309 IA = DO 312 IC=1,2 IB = DO 311 ID=1,2 DO 310 I=1,3 DO 310 J=1,3 DO 310 K=1,3 SA(I+IA+II,J+IB+IJ) = SA(I+IA+II,J+IB+IJ) + +ST(I+IA,K+IB)*V(J,K) 310 CONTINUE IB = IB + 311 CONTINUE IA = IA + 312 CONTINUE 313 CONTINUE C COMPUTE [M]=[V]T[m][V] IA = DO 316 IC=1,4 IB = DO 315 ID=1,4 DO 314 I=1,3 DO 314 J=1,3 DO 314 K=1,3 S(I+IA,J+IB) = S(I+IA,J+IB)+V(I,K)*SA(K+IA,J+IB) 314 CONTINUE IB = IB + 315 CONTINUE IA = IA + 316 CONTINUE C FORM LUMPED MASS MATRIX (INCLUDING ROTORY AND TORSIONAL INERTIA) XMASS = D(7)*D(1)*XL/2.0D0 P(4) = XMASS*D(2)/D(1) P(5) = XMASS*XL*XL/12.0D0 P(6) = XMASS*XL*XL/12.0D0 DO 502 I=1,3 P(I) = XMASS P(I+6)= XMASS 502 CONTINUE DO 504 I=1,3 P(I+9) = P(I+3) 504 CONTINUE C -RETURN END C ====================================================================== SUBROUTINE GEOM06 (S,D,UL,V,ST,SA,XL,NST,NDF) C C FORM GEOMETRIC STIFFNESS MATRIX BY SUNIL 00/11/23 C IMPLICIT REAL*8 (A-H,O-Z) DIMENSION S(NST,1),D(1),UL(NDF,1) DIMENSION V(3,3),ST(6,6),SA(12,12) C DETERMINE AXIAL FORCE: FORCE CALL STIF06 (S,D,ST,SA,V,XL,NST,7) UT = 0.0D0 AXIAL = 0.0D0 DO 705 J=1,6 705 UT = UT + SA(1,J)*UL(J,1) DO 715 J=7,12 715 UT = UT + SA(1,J)*UL(J-6,2) IF (D(8).NE.0.0D0 AND D(9).EQ.1.0D0) AXIAL = 0.5D0*D(8)*XL FORCE = UT - AXIAL !!! CHECKING LATER C REINITIALIZE ARRAYS DO 725 I=1,12 DO 725 J=1,12 SA(I,J) = 0.0D0 725 CONTINUE DO 735 I=1,6 DO 735 J=1,6 ST(I,J) = 0.0D0 735 CONTINUE DO 745 I=1,NST DO 745 J=1,NST S(I,J) = 0.0D0 745 CONTINUE C FORM LOCAL GEOMENTRIC STIFFNESS MATRIX [kg] DO 313 IS=1,4 GOTO (1,2,3,4),IS CON = FORCE/XL/30.0D0 TMP1 = CON*30.0D0 TMP2 = TMP1*D(2)*D(6)/D(1)/D(5)/XL ST(1,1) = TMP1 ST(4,4) = TMP2 ST(2,2) = CON*36.0D0 ST(3,3) = ST(2,2) ST(5,5) = CON*4.0D0*XL*XL ST(6,6) = ST(5,5) ST(2,6) = CON*3.0D0*XL ST(6,2) = ST(2,6) ST(3,5) =-ST(2,6) ST(5,3) =-ST(2,6) II = IJ = GOTO 309 ST(2,6) =-ST(2,6) ST(6,2) =-ST(6,2) ST(3,5) =-ST(3,5) ST(5,3) =-ST(5,3) II = IJ = GOTO 309 ST(2,2) =-ST(2,2) ST(3,3) =-ST(3,3) ST(5,5) =-ST(5,5)/4.0D0 ST(6,6) = ST(5,5) ST(6,2) =-ST(2,6) ST(5,3) =-ST(3,5) ST(1,1) =-ST(1,1) ST(4,4) =-ST(4,4) II = IJ = GOTO 309 ST(2,6) =-ST(2,6) ST(6,2) =-ST(6,2) ST(3,5) =-ST(3,5) ST(5,3) =-ST(5,3) II = IJ = C COMPUTE [kg][V]; ([kg]-LOCAL GEO STIFF, [V]-COOR TRANSFORM MAT.) 309 IA = DO 312 IC=1,2 IB = DO 311 ID=1,2 DO 310 I=1,3 DO 310 J=1,3 DO 310 K=1,3 SA(I+IA+II,J+IB+IJ) = SA(I+IA+II,J+IB+IJ) + 310 +ST(I+IA,K+IB)*V(J,K) CONTINUE IB = IB + 311 CONTINUE IA = IA + 312 CONTINUE 313 CONTINUE C COMPUTE [KG]=[V]T[kg][V] IA = DO 316 IC=1,4 IB = DO 315 ID=1,4 DO 314 I=1,3 DO 314 J=1,3 DO 314 K=1,3 S(I+IA,J+IB) = S(I+IA,J+IB)+V(I,K)*SA(K+IA,J+IB) 314 CONTINUE IB = IB + 315 CONTINUE IA = IA + 316 CONTINUE C -RETURN END ... khung nhà công nghiệp 17 1.3.1 Tính toán nội lực khung phương pháp giải tích 17 1.4.1 Tính toán nội lực khung phương pháp số 21 CHƯƠNG 2: TÍNH TOÁN KẾT CẤU KHUNG NHÀ CÔNG NGHIỆP BẰNG PHƯƠNG... hình khung nhà công nghiệp Kết cấu khung mái dốc đối xứng Khung nhà công nghiệp có tầng lửng Khung công nghiệp có cột Nhà công nghiệp có cột chịu moment đỉnh kèo Khung công nghiệp có mái cong Khung. .. toán tính toán nội lực kết khung thép nhà công nghiệp bằng phương pháp phần tử hữu hạn  Sử dụng kết hợp phần mềm GID Fortran để thiết lập, xây dựng chương trình tính nội lực khung thép nhà công