中南大学铁道校区(土木2012级研究生结构分析程序设计基础试题(考试方式:开卷,交卷时间:第13周1.如图示平面应力悬臂梁模型2423222120191817161514131211987654321012019181716151413121110987654321图1平面应力悬臂梁模型试写出其节点联系数组IT、单元联系数组LMT和对角元地址数组MAXA。2、求解线性方程组[A]{X}={B}式中的[A]为对称正定矩阵,其值为4.50.21.300005.30000010.25.101.70[]8.40000.6003.105.6A-????????-??=????????????对称{B}矩阵为{B}=[3.4,5.5,12.3,13.5,0.6,1.4,5.6]T要求写出矩阵分解和回代求解过程。3.在计算机上调试空间桁架分析程序TRUSS,并用两个算例进行考证,要求写出相关计算结果。应大家要求,程序已经调试,附录于下:CCCTRUSS.FORCC*******************************************************************************SPACIALTRUSSSTRUCTUREANALYSIS**DebuggedbyWYL2013.05*******************************************************************************PROGRAMTRUSSIMPLICITREAL*8(A-H,O-ZIMPLICITINTEGER*4(I-NCHARACTERNAME*40COMMON/AT/A(18000000COMMON/IAT/IA(2000000A=0;IA=0WRITE(*,*'INPUTFILENAME?';READ(*,*NAME;CALLOPENF(NAMECALLDATAIN(NP,NE,NF,ND,NDF,NPF,NM,NR,NCF,&IME,INAE,IIT,ILMT,IMAXA,&IX,IY,IZ,IRR,IAE,IPF,ICKKCALLFLMT(NP,NE,NN,NN1,NR,A(IRR+1,ND,NF,NDF,&IA(IME+1,IA(IIT+1,IA(ILMT+1CALLFMAXA(NN1,NE,IA(ILMT+1,IA(IMAXA+1,NWK,NPF,NDF,&ICKK,IDIST,IFTOOL,IFF,IPP,ISG,ISMCALLCONKB(NP,NE,NM,NWK,IA(IME+1,A(IX+1,&A(IY+1,A(IZ+1,A(IAE+1,IA(INAE+1,&IA(ILMT+1,IA(IMAXA+1,A(ICKK+1,NN1CALLMKFORCE(NP,NF,NPF,NCF,NN,&IA(IIT+1,A(IPF+1,A(IPP+1,A(IFTOOL+1CALLLDLT(A(ICKK+1,IA(IMAXA+1,NN,1,3,NWK,NN1CALLRESOLVE(A(ICKK+1,A(IFTOOL+1,IA(IMAXA+1,NN,NWK,NN1CALLDISPLS(NP,NE,NF,NPF,NM,NN,IA(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+1CALLDATAOUT(NP,NE,NPF,A(IDIST+1,A(IFF+1,A(ISG+1,A(ISM+1CALLCLOSEFENDCSUBROUTINEOPENF(NAMECHARACTERNAME*40NUM=0DOI=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(1CLOSE(2CLOSE(3RETURNENDCSUBROUTINEDATAIN(NP,NE,NF,ND,NDF,NPF,NM,NR,NCF,&IME,INAE,IIT,ILMT,IMAXA,&IX,IY,IZ,IRR,IAE,IPF,ICKKIMPLICITREAL*8(A-H,O-ZIMPLICITINTEGER*4(I-NCOMMON/AT/A(18000000COMMON/IAT/IA(2000000READ(1,*NP,NE,NM,NR,NCFWRITE(2,701NP,NE,NM,NR,NCF701FORMAT(//1X,'###OUTPUTOFORIGINALINPUTINFORMATION###'&//5X,'NumberofjointsJOINTS=',I5&/5X,'NumberofelementsELEMENTS=',I5&/5X,'NumberofmaterialpropertygroupsPROPERTYTYPES=',I5&/5X,'NumberofrestrainedjointsRESTRAINTS=',I5&/5X,'NumberofconcentrativeforcedjointsNCF=',I5C--------------FORMPOINTER----------------------------------NF=3ND=2NDF=ND*NFNPF=NP*NFIME=0INAE=IME+2*NEIIT=INAE+NEILMT=IIT+NF*NPIMAXA=ILMT+NDF*NEIX=0IZ=IY+NPIRR=IZ+NPIAE=IRR+2*NRIPF=IAE+2*NMICKK=IPF+4*NCFREAD(1,*(A(IX+I,A(IY+I,A(IZ+I,I=1,NPWRITE(2,714(I,A(IX+I,A(IY+I,A(IZ+I,I=1,NP714FORMAT(//5X,'GENERATEDJOINTCOORDINATESDATA'&/1X,'JOINT',15X,'X',13X,'Y',13X,'Z'&/(4X,I5,3X,3(2X,E12.6READ(1,*(A(IAE+2*(I-1+1,A(IAE+2*(I-1+2,I=1,NMREAD(1,*(IA(IME+2*(I-1+1,IA(IME+2*(I-1+2,IA(INAE+I,I=1,NEWRITE(2,606(I,A(IAE+2*(I-1+1,A(IAE+2*(I-1+2,I=1,NMWRITE(2,607(I,IA(IME+2*(I-1+1,IA(IME+2*(I-1+2,&IA(INAE+I,I=1,NE606FORMAT(/5X,'ELEMENTMATERAILPROPERTIESDATA'&/2X,'NO.',10X,'E',10X,'Ax'&/(2X,I3,2(1X,E11.5607FORMAT(/5X,'TRUSSELEMENTDEFINITIONDATA'&/2X,'NO.',10X,'JOINT_1',10X,'JOINT_2',10X,'NAE'&/(2X,I3,3(10X,I5READ(1,*(A(IRR+2*(I-1+1,A(IRR+2*(I-1+2,I=1,NR...
