安徽建筑工程学院计算结构力学5.pptVIP

下载本文档

3
0
约5.64千字
约 56页
2017-05-19 发布于北京
举报
版权申诉

安徽建筑工程学院计算结构力学5.ppt

1、本文档共56页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

第五章结构刚度矩阵与荷载向量与未知量i在同一单元的未知量叫做未知量i的相关未知量。若i，j相关，则Kij≠0 若i，j不相关，则Kij=0 主程序,数据输入等计算各单元的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.EMENTS 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( I,J),I=1,4),J=1,NJT) READ( 1,* ) (JMH(I),I= 1,NE) READ( l,* ) (EA(I),I= 1,NAI) READ( l,* ) (EI(I),I=1,NAI) READ(1,* ) (X(I),I=1,NJ) READ( l,* ) ( Y(l),I=1, NJ) CALL QJW ( NJ, NJT,JTX,JW,N) DO 10 M= 1,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 = I5,/6 X, ’MW * =’,6I5) 50 FORMAT(9 F9.1 ) STOP 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),ZK(50,50) 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 * MATRAIX‘ WRHE( *,15) M 15 FORMAT(1X,’ ELEMENT NO.=‘,I5) 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