中南大学研究生结构程序设计基础考试试题概要

中南大学铁道校区(土木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,NRWRITE(2,608(A(IRR+2*(I-1+1,A(IRR+2*(I-1+2,I=1,NR608FORMAT(/5X,'JOINTRESTRAINTSDATA'&/2X,'JOINT',10X,'RESTRAINT',&/(2X,F7.0,10X,F9.3READ(1,*((A(IPF+4*(I-1+J,J=1,4,I=1,NCFWRITE(2,609((A(IPF+4*(I-1+J,J=1,4,I=1,NCF609FORMAT(/5X,'CONCENTRATIVEFORCEDJOINTSDATA'&/2X,'JOINT',10X,'Fx',10X,'Fy',10X,'Fz'&/(2X,F7.0,3(1X,E12.6RETURNENDCSUBROUTINEMKFORCE(NP,NF,NPF,NCF,NN,IT,PF,PP,FTOOLIMPLICITREAL*8(A-H,O-ZIMPLICITINTEGER*4(I-NDIMENSIONIT(NF,NP,PF(4,NCF,PP(NPF,FTOOL(NPFPP=0;FTOOL=0DOI=1,NCFNOD=PF(1,IDOJ=1,NFPP(NF*(NOD-1+J=PF(J+1,IENDDODOI=1,NPDOJ=1,NFLAB=IT(J,IIF(LAB.GT.0.AND.LAB.LE.NNTHENFTOOL(LAB=PP(NF*(I-1+JENDIFENDDOENDDORETURNENDCCSUBROUTINEDATAOUT(NP,NE,NPF,DIST,FF,SG,SMIMPLICITREAL*8(A-H,O-ZIMPLICITINTEGER*4(I-NDIMENSIONDIST(NPF,FF(NPF,SG(NE,SM(NEWRITE(2,715(I,(DIST(3*(I-1+J,J=1,3,I=1,NP715FORMAT(//5X,'SOLVEDJOINTDISPLACEMENTSDATA'&/1X,'JOINT',3X,8X,'Dx',12X,'Dy',12X,'Dz'&/(4X,I5,3X,3(2X,E12.6WRITE(2,716(IE,SG(IE,SM(IE,IE=1,NE716FORMAT(//5X,'SOLVEDELEMENTINTERNALFORCEDATA'&/1X,'ELEMENT',3X,8X,'Nx',8X,'STRESS'&/(4X,I5,3X,2(2X,F12.6WRITE(2,717(I,(FF(3*(I-1+J,J=1,3,I=1,NP717FORMAT(//5X,'SOLVEDJOINTREACTIONDATA'&/1X,'JOINT',3X,8X,'Rx',12X,'Ry',12X,'Rz'&/(4X,I5,3X,3(2X,f12.4RETURNENDCCSUBROUTINEFLMT(NP,NE,NN,NN1,NR,RR,ND,NF,NDF,ME,IT,LMTIMPLICITREAL*8(A-H,O-ZIMPLICITINTEGER*4(I-NCThisprogramformsthejoint&elementnumberingmatrixIT&LMTDIMENSIONIT(NF,NP,LMT(NDF,NE,ME(ND,NE,RR(2,NRNN=0;NN1=0;IT=0;LMT=0N=0DOI=1,NPC=0DOK=1,NRKR=RR(1,KIF(KR.EQ.IC=RR(2,KENDDONC=CC=C-NCDOJ=1,NFC=C*10.0L=C+0.1C=C-LIF(L.EQ.0THENN=N+1IT(J,I=NELSEIT(J,I=0ENDIFENDDOENDDONN=NNN1=NN+1DOIE=1,NEDOI=1,NDNI=ME(I,IEDOJ=1,NFLMT((I-1*NF+J,IE=IT(J,NIENDDOENDDOENDDORETURNENDCSUBROUTINEFMAXA(NN1,NE,LMT,MAXA,NWK,NPF,NDF,&ICKK,IDIST,IFTOOL,IFF,IPP,ISG,ISMCThisprogramformstheMDEaddressmatrixMAXAof[K]IMPLICITREAL*8(A-H,O-ZIMPLICITINTEGER*4(I-NDIMENSIONMAXA(NPF,LMT(NDF,NEMAXA=0