资源描述：
《杆件结构分析程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、c program frameC PROGRAMMAIN INTEGER RR,AA,GG,EP,RP,T REALKK CHARACTER*12FNAME,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),
2、FI3(20,6) COMMON/A5/AA(200),KK(6000),NX,NH COMMON/A6/FF(60),UV(20,6) WRITE(0,101) 101 FORMAT(///5X,'PLEASEINPUTDATAFILENAME:') READ(*,100)FNAME 100 FORMAT(8A) WRITE(0,102) 102 FORMAT(///5X,'PLEASEINPUTOUTPUTDATAFILENAME:') READ(*,103)FNAOT 103
3、 FORMAT(8A) OPEN(5,FILE=FNAME,STATUS='OLD') OPEN(6,FILE=FNAOT,STATUS='new') CALLINPUT CALLFORRR CALLFORAA CALLFORKK CALLDIVKK DO200N9=1,MX CALLFORFF(N9) CALLFORDIS CALLFORFM CALLFORRC 200 CONTINUE STOP END SUBR
4、OUTINEINPUT 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,MX50 FORMAT(2X,'NUMBEROFELEMENTNE=',I2/2X, #'NUMBEROFNODENP=',I2/2X,
5、 #'NUMBEROFSUPPOTNODENR=',I2/2X, #'NUMBEROFMATERIALSETMT=',I2/2X, #'NUMBEROFWORKINGCASEMX=',I2) READ(5,*)((AE(I,J),J=1,3),I=1,MT) WRITE(6,221)221 FORMAT(2X,'N0.OFSET',2X,'MODULUSOFELASTICITY',2X, #'CROSS-AREA',2X,'MOMEMTOFINERTIA') DO15I=1,MT
6、 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,'DATAOFELEMENT') DO17I=1,NE WRITE(6,193)I,(EP(I,J),J=1,3)193 FORMAT(2X,I4,5X,3I8)17 CONTINUE READ(5,*)((RP(I,
7、J),J=1,4),I=1,NR) WRITE(6,224)224 FORMAT(10X,'DATAOFSUPPORT') 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 FORMAT(10X,'NODALCOORDINATE'/5X,'N0.NODE',8
8、X,'X-X' #,8X,'Z-Z') DO227I=1,NP WRITE(6,194)I,(XZ(I,J),J=1,2)194 FORMAT(5X,I4,5X,2F12.6)227 CONTINUE RETURN
