资源描述:
《结构程序设计基础作业》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、结构程序设计基础3、结构分析程序设计实践题:上机调试教材中的TRUSS程序,并要求给出具体算例。解答:CCCTRUSS.FORC*******************************************************************************SPACIALTRUSSSTRUCTUREANALYSIS*******************************************************************************P
2、ROGRAMTRUSSIMPLICITREAL*8(A-H,O-Z)IMPLICITINTEGER*4(I-N)CHARACTERNAME*40COMMON/AT/A(18000000)COMMON/IAT/IA(2000000)A=0;IA=0WRITE(*,*)'INPUTFILENAME?';READ(*,*)NAME;CALLOPENF(NAME)CALLDATAIN(NP,NE,NF,ND,NDF,NPF,NM,NR,NCF,&IME,INAE,IIT,ILMT,IMAXA,&IX,IY,
3、IZ,IRR,IAE,IPF,ICKK)CALLFLMT(NP,NE,NN,NN1,NR,A(IRR+1),ND,NF,NDF,&IA(IME+1),IA(IIT+1),IA(ILMT+1))CALLFMAXA(NN1,NE,IA(ILMT+1),IA(IMAXA+1),NWK,NPF,NDF,&ICKK,IDIST,IFTOOL,IFF,IPP,ISG,ISM)CALLCONKB(NP,NE,NM,NWK,IA(IME+1),A(IX+1),&A(IY+1),A(IZ+1),A(IAE+1),IA
4、(INAE+1),&IA(ILMT+1),IA(IMAXA+1),A(ICKK+1),NN1)CALLMKFORCE(NP,NF,NPF,NCF,NN,&IA(IIT+1),A(IPF+1),A(IPP+1),A(IFTOOL+1))CALLLDLT(A(ICKK+1),IA(IMAXA+1),NN,1,3,NWK,NN1)CALLRESOLVE(A(ICKK+1),A(IFTOOL+1),IA(IMAXA+1),NN,NWK,NN1)CALLDISPLS(NP,NE,NF,NPF,NM,NN,IA
5、(IIT+1),A(IFTOOL+1),&A(IDIST+1),A(IAE+1),IA(IME+1),IA(INAE+1),A(IX+1),A(IY+1),&A(IZ+1),A(IPP+1),A(IFF+1),A(ISG+1),A(ISM+1))CALLDATAOUT(NP,NE,NPF,A(IDIST+1),A(IFF+1),A(ISG+1),A(ISM+1))CALLCLOSEFENDCSUBROUTINEOPENF(NAME)CHARACTERNAME*40NUM=018结构程序设计基础DOI
6、=1,40IF(NAME(I:I).NE.'')NUM=NUM+1ENDDOOPEN(1,FILE=NAME(1:NUM),STATUS='UNKNOWN')OPEN(2,FILE=NAME(1:NUM)//'.RES',STATUS='UNKNOWN')OPEN(3,FILE=NAME(1:NUM)//'.ERO',STATUS='UNKNOWN')RETURNENDCCSUBROUTINECLOSEFCLOSE(1)CLOSE(2)CLOSE(3)RETURNENDCSUBROUTINEDATA
7、IN(NP,NE,NF,ND,NDF,NPF,NM,NR,NCF,&IME,INAE,IIT,ILMT,IMAXA,&IX,IY,IZ,IRR,IAE,IPF,ICKK)IMPLICITREAL*8(A-H,O-Z)IMPLICITINTEGER*4(I-N)COMMON/AT/A(18000000)COMMON/IAT/IA(2000000)READ(1,*)NP,NE,NM,NR,NCFWRITE(2,701)NP,NE,NM,NR,NCF701FORMAT(//1X,'###OUTPUTOFO
8、RIGINALINPUTINFORMATION###'&//5X,'NumberofjointsJOINTS=',I5&/5X,'NumberofelementsELEMENTS=',I5&/5X,'NumberofmaterialpropertygroupsPROPERTYTYPES=',I5&/5X,'NumberofrestrainedjointsRESTRAINTS=',I5&/5X,'NumberofconcentrativeforcedjointsNCF=