有限元桁架源程序.doc

下载文档

15
0
约3.54千字
约 4页
2017-05-26 发布于河南
举报
版权申诉
保障服务

有限元桁架源程序.doc

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

有限元桁架源程序

C-------TRUSS.FOR---(FOR TRUSS PROBLEMS)-------------------- DIMENSION LOC(3,2),CX(3),CY(3),IFIX(3),AL(3),AREA(3), 1 GK(6,6),F(6),FF(3) COMMON NN,NE,ND,NFIX,E C----------------------------------------------------------- OPEN(5,FILE=T1.TXT,STATUS=OLD) OPEN(6,FILE=OUTT1.) READ(5,*) NN,NE,ND,NFIX,E WRITE(6,105) NN,NE,ND,NFIX,E 105 FORMAT(2X,NN NE ND NFIX E/4I5,E11.4) READ(5,*) (LOC(I,1),LOC(I,2),AREA(I),I=1,NE) WRITE(6,108) 108 FORMAT(/1X,ELEMENT NODE1 NODE2 AREA) WRITE(6,110) (I,LOC(I,1),LOC(I,2),AREA(I),I=1,NE) 110 FORMAT(1X,3I6,E12.4) READ(5,*) (CX(J),CY(J),J=1,NN) WRITE(6,112) 112 FORMAT(/3X,NODE X-COORD Y-COORD) WRITE(6,115) (J,CX(J),CY(J),J=1,NN) 115 FORMAT(3X,I3,3X,2E12.4) READ(5,*) (IFIX(K),K=1,NFIX) WRITE(6,118) (IFIX(K),K=1,NFIX) 118 FORMAT(/1X,IFIX=,5I4) C----------------------------------------------------------- DO 10 I=1,ND 10 F(I)=0.0 F(2)=-10000.0 C----------------------------------------------------------- CALL CST(LOC,CX,CY,IFIX,AL,AREA,GK,F,FF) WRITE(6,120) 120 FORMAT(/3X,NODE,6X,X-DISP,8X,Y-DISP) WRITE(6,125) (I,F(2*I-1),F(2*I),I=1,NN) 125 FORMAT(1X,I5,2E15.4) WRITE(6,130) 130 FORMAT(/1X,ELEMENT,5X,FORCE) WRITE(6,135) (I,FF(I),I=1,NE) 135 FORMAT(1X,I5,E15.4) STOP END C----------------------------------------------------------- SUBROUTINE CST(LOC,CX,CY,IFIX,AL,AREA,GK,F,FF) DIMENSION LOC(NE,2),CX(NN),CY(NN),IFIX(NFIX),AL(NE), 1 AREA(NE),GK(ND,ND),F(ND),FF(NE),EK(4,4),XX(4) COMMON NN,NE,ND,NFIX,E C----------------------------------------------------------- DO 10 I=1,ND DO 10 J=1,ND 10 GK(I,J)=0.0 C----------------------------------------------------------- DO 100 I=1,NE I1=LOC(I,1) I2=LOC(I,2) X12=CX(I2)-CX(I1) Y12=CY(I2)-CY(I1) AL(I)=SQRT(X12**2+Y12**2) CS=X12/AL(I) SN=Y12/AL(I) A0=AREA(I)