原创力文档- 1、本文档共23页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
杆件结构分析程序.
c???? program? \frame\
C???? PROGRAM MAIN
????? INTEGER? RR,AA,GG,EP,RP,T
????? REAL KK
????? CHARACTER*12 FNAME,FNAOT
????? COMMON/A1/NE,NP,NR,MT,MX
????? COMMON/A2/ME(2),RR(20,3),GG(6),T
????? COMMON/A3/EP(20,3),RP(20,4),XZ(30,2),AE(5,3)
????? COMMON/A4/KI2,KI3,MI2(20),MI3(20),FI2(20,3),FI3(20,6)
????? COMMON/A5/AA(200),KK(6000),NX,NH
????? COMMON/A6/FF(60),UV(20,6)
????? WRITE(0,101)
?101? FORMAT(///5X,PLEASE INPUT DATA FILENAME:)
????? READ(*,100) FNAME
?100?? FORMAT(8A)
????? WRITE(0,102)
?102? FORMAT(///5X,PLEASE INPUT OUTPUT DATA FILENAME:)
????? READ(*,103) FNAOT
?103? FORMAT(8A)
????? OPEN(5,FILE=FNAME,STATUS=OLD)
????? OPEN(6,FILE=FNAOT,STATUS=new)
????? CALL INPUT
??? CALL FORRR
????? CALL FORAA
????? CALL FORKK
????? CALL DIVKK
????? DO 200 N9=1,MX
????? CALL FORFF(N9)
????? CALL FORDIS
????? CALL FORFM
????? CALL FORRC
?200? CONTINUE
????? STOP
????? END
????? SUBROUTINE INPUT
????? INTEGER? EP,RP
????? COMMON/A1/NE,NP,NR,MT,MX
????? COMMON/A3/EP(20,3),RP(20,4),XZ(30,2),AE(5,3)
????? READ(5,*)NE,NP,NR,MT,MX
????? WRITE(6,50)NE,NP,NR,MT,MX
????? WRITE(0,50)NE,NP,NR,MT,MX
50??? FORMAT(2X,NUMBER OF ELEMENT NE=,I2/2X,
???? #NUMBER OF NODE NP=,I2/2X,
???? #NUMBER OF SUPPOT NODE NR=,I2/2X,
???? #NUMBER OF MATERIAL SET MT=,I2/2X,
???? #NUMBER OF WORKING CASE MX=,I2)
????? READ(5,*)((AE(I,J),J=1,3),I=1,MT)
????? WRITE(6,221)
221?? FORMAT(2X,N0. OF SET,2X,MODULUS OF ELASTICITY,2X,
???? #CROSS-AREA,2X,MOMEMT OF INERTIA)
????? DO 15 I=1,MT
????? WRITE(6,222) I,(AE(I,J),J=1,3)
222?? FORMAT(2X,I4,10X,3F16.5)
15??? CONTINUE
????? READ(5,*)((EP(I,J),J=1,3),I=1,NE)
????? WRITE(6,229)
229?? FORMAT(10X,DATA OF ELEMENT)
????? DO 17 I=1,NE
????? WRITE (6,193)I,(EP(I,J),J=1,3)
193?? FORMAT(2X,I4,5X,3I8)
17??? CONTINUE
????? READ(5,*)((RP(I,J),J=1,4),I=1,NR)
????? WRITE(6,224)
224?? FORMAT(10X,DATA OF SUPPORT)
????? WRITE(6,96)((RP(I,J),J=1,4),I=1,NR)
????? WRITE(0,96)((RP(I,J),J=1,4),I=1,NR)
96??? FORMAT(2X,4I6)
????? READ(5,*)((XZ(I,J),J=1,2),I=1,NP)
????? WRITE(6,226)
226
文档评论（0）