计算结构力学第五章.ppt

下载文档 降价啦

12
0
约5.61千字
约 57页
2017-06-24 发布于湖北
举报
版权申诉
保障服务

计算结构力学第五章.ppt

1、本文档共57页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

计算结构力学第五章

第五章结构刚度矩阵与荷载向量主程序,数据输入等计算各单元的L,C,S 计算XSA(NE,7) 形成单元刚度矩阵形成结构刚度矩阵结束形成JW数组形成MW数组程序设计框图 4、形成结构刚度[K]的源程序设计 C THE PROGRAM OF KJEXAM DIMENSION JW ( 3,20 ),J TX (4,20), JH ( 2,20 ), * MW (6), JMH ( 20 ) REA L * 8 CX (20), SY ( 20 ), SL (20),EA ( 5 ), * X (20),Y(20),XSA ( 20,7 ),ZK (50,50), El (5) WRITE(*,*) FINDING THE MW OF * EI.EMENT OPEN( 1, FILE = KJE. DA T ) 数据文件名为KJE.DAT READ( 1,*) NE, NJ,NJT,NAI READ( 1,* ) ( (JH( I,J ),I = 1.2),J =1,NE ) READ( 1,* ) ((JTX( IJ),i=1,4),J=1,NJT) READ( 1,* ) (JMH(I),I= 1,NT) READ( l,* ) (EA(I),I= 1,NAI) READ( l,* ) (EI(I),I=I,NAI) READ(1,* ) (X(I),I=I,NJ) READ( l,* ) ( Y(l),I=1, NJ) CALL QMW ( NJ, NJT,JTX,JW,N) DO 10 M= I,NE CALL QMW ( M,NE,NJ,JH, JW,MW) WRITE( *, 100)M, ( MW(I), I = 1,6) 10 CALL DCH ( NE, NJ,JH,X, Y, SL,CX, S Y ) CALL QXS(NE,NAI,JH,JMH,SL, CX,SY,EA, * EI,XSA) CALL KJX 1 ( NE, N, XSA,NJ, JH,JW,ZK ) WRITE( *,* )THIS IS STRUCTURAL * MATRAIX ZK’ WRITE(*,50) ((ZK(I,J),I= I,N),J= I,N) 100 FORMAT ( 1X,’ ELEMT, NO = 15,/6 X, ’MW * =’,6I5) 50 FORMAT(9 F9.1 ) STOP 44 END SUBROUTINE KJX 1 (NE, N, XSA, NJ, JH, * JW,ZK ) DIMENSION JH(2,NE),JW(3,NJ ),MW(6) REAL * 8 DK ( 6 6), ZK (50, 50),XSA( NE,7), * XS (7) CALL ZERO2 ( ZK, 50,50) DO 20 M=1,NE DO 25 I= 1,7 XS ( I ) =XSA ( M,I) 25 CONTINUE CALL DKX( XS,DK ) WRITE( *,* )THIS IS ELEMENT * MATRIX WRHE( *,11) ((DK(I,J),l=l,6),J=l,6) 11 FORMAT(1X,6FI0.2) CALL QMW ( M, NE, NJ,JH,JW,MW ) DO 30 I=1,6 L= MW(I) IF ( L. LE. 0 ) GO TO 30 DO 40 J=1,6 K= MW(J) IF (K. LE.0) GO TO 40 ZK(L, K) = ZK (L,K)+DK (I, J) 40 CONTINUE 30