资源描述:
《杆件结构分析程序》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、programframeCPROGRAMMAININTEGERRR,AA,GG,EP,RP,TREALKKCHARACTER*12FNAME,FNAOTCOMMON/A1/NE,NP,NR,MT,MXCOMMON/A2/ME(2),RR(20,3),GG(6),TC0MM0N/A3/EP(20,3),RP(20,4),XZ(30,2),AE(5,3)C0MM0N/A4/KI2,KI3,MI2(20),MI3(20),FI2(20,3),FI3(20,6)C0MM0N/A5/AA(200),KK(6000),NX,NHC0MM0N/A6/FF(60),UV(20,6)WRITE(0
2、,101)101FORMAT(///5X,JPLEASEINPUTDATAFTLENAME/)READ(*,100)FNAME100FORMAT(8A)WRITE(0,102)102FORMAT(///5X,?PLEASEINPUTOUTPUTDATAFILENAME/)READ(*,103)FNAOT103FORMAT(8A)OPEN(5,FILE二FNAME,STATUS二'OLD')OPEN(6,F1LE=FNAOT,STATUS二'new')CALLINPUTCALLFORRRCALLFORKKCALLDIVKKDO200N9=1,MXCALLFORFF(N9)CALLFOR
3、DISCALLFORFMCALLFORRC200CONTINUESTOPENDSUBROUTINEINPUTINTEGEREP,RPC0MM0N/A1/NE,NP,NR,MT,MXC0MM0N/A3/EP(20,3),RP(20,4),XZ(30,2),AR(5,3)READ(5,*)NE,NP,NR,MT,MXWRITE(6,50)NE,NP,NR,MT,MXWRITE(0,50)NE,NP,NR,MT,MX50FORMAT(2X,'NUMBEROFELEMENTNE=,,T2/2X,#'NUMBEROFNODENP二',I2/2X,#'NUMBEROFSUPPOTNODENR二'
4、,I2/2X,#'NUMBEROFMATERIALSETMT二',12/2X,#'NUMBEROFWORKINGCASEMX二',12)READ(5,*)((AE(I,J),J二1,3),I二1,MT)WRITE(6,221)221FORMAT(2X,'NO.OFSET',2X,'MODULUSOFELASTICITY',2X,#'CROSS-AREA,,2X,'MOMEMTOFINERTIA,)DO15[二1,MTWRITE(6,222)I,(AE(I,J),J=l,3)222FORMAT(2X,14,1OX,3F16.5)15CONTINUEREAD(5,*)((EP(I,J),
5、J=l,3),I=l,NE)WRITE(6,229)229FORMAT(1OX,'DATAOFELEMENT,)DO171=1,NEWRITE(6,193)1,(EP(I,J),J=l,3)193FORMAT(2X,14,5X,318)17CONTINUEREAD(5,*)((RP(I,J),J二1,4),I二1,NR)WRITE(6,224)224FORMAT(1OX,'DATAOFSUPPORT')WRITE(6,96)((RP(I,J),J二1,4),1=1,NR)WRITE(0,96)((RP(I,J),J二1,4),I二1,NR)96FORMAT(2X,416)READ(5
6、,*)((XZ(I,J),J=l,2),1=1,NP)WRITE(6,226)226FORMAT(1OX,'NODALCOORDINATE,/5X,'NO.NODE',8X,'X-X'#,8X,'Z-Z')DO2271=1,NPWRITE(6,194)1,(XZ(I,J),J=l,2)194FORMAT(5X,14,5X,2F12.6)227CONTINUERETURNENDSUBROUTINEFORRRINTEGERRR,AA,GG,EP,RP,TREALKKCOMMON/A1/NE,NP,NR,MT,MXCOMMON/A2/ME(2),RR(20,3),GG(6),TCOMMON
7、/A3/EP(20,3),RP(20,4),XZ(30,2),AE(5,3)COMMON/A5/AA(200),KK(6000),NX,NHDO483T=1,NPDO483J二1,3RR(I,J)=l483CONTINUEDO4701=1,NRTZ1=RP(T,1)IW1=RP(I,2)IW2=RP(I,3)1W3=RP(1,4)IF(IWl.EQ.0)RR(IZl,1)=0IF(IW2.EQ.0)RR(IZl,2)=0IF(IW3.EQ.0)RR(IZ1