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

1、中南大学铁道校区(土木2012级研究生结构分析程序设计基础试题(考试方式:开卷,交卷时间:第13周1.如图示平面应力悬臂梁模型2423222120191817161514131211987654321012019181716151413121110987654321图1 平面应力悬臂梁模型试写出其节点联系数组IT 、单元联系数组LMT 和对角元地址数组MAXA 。2、求解线性方程组AX=B 式中的A为对称正定矩阵,其值为4.50.21.300005.30000010.2 5.10 1.708.40000.6003.105.6A -=对称B矩阵为B=3.4,5.5,12.3,13.5,0.6,1

2、.4,5.6T要求写出矩阵分解和回代求解过程。3.在计算机上调试空间桁架分析程序TRUSS ,并用两个算例进行考证,要求写出相关计算结果。应大家要求,程序已经调试,附录于下:CCC TRUSS.FORCC* * SPACIAL TRUSS STRUCTURE ANAL YSIS * * Debugged by WYL 2013.05 * * PROGRAM TRUSSIMPLICIT REAL*8 (A-H,O-ZIMPLICIT INTEGER*4 (I-NCHARACTER NAME*40COMMON /A T/A(18000000COMMON /IAT/IA(2000000A=0;IA=

3、0WRITE(*,*INPUT FILE NAME?;READ(*,*NAME;CALL OPENF(NAMECALL DATAIN(NP,NE,NF,ND,NDF,NPF,NM,NR,NCF,& IME,INAE,IIT,ILMT,IMAXA,& IX,IY,IZ,IRR,IAE,IPF,ICKKCALL FLMT(NP,NE,NN,NN1,NR,A(IRR+1,ND,NF,NDF,& IA(IME+1,IA(IIT+1,IA(ILMT+1CALL FMAXA(NN1,NE,IA(ILMT+1,IA(IMAXA+1,NWK,NPF,NDF, & ICKK,IDIST,IFTOOL,IFF,I

4、PP,ISG,ISMCALL CONKB(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,NN1CALL MKFORCE(NP,NF,NPF,NCF,NN,& IA(IIT+1,A(IPF+1,A(IPP+1,A(IFTOOL+1CALL LDLT(A(ICKK+1,IA(IMAXA+1,NN,1,3,NWK,NN1CALL RESOLVE(A(ICKK+1,A(IFTOOL+1,IA(IMAXA+1,NN,NWK,NN1CALL DISPLS(NP,N

5、E,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+1CALL DATAOUT(NP,NE,NPF,A(IDIST+1,A(IFF+1,A(ISG+1,A(ISM+1CALL CLOSEFENDCSUBROUTINE OPENF(NAMECHARACTER NAME*40NUM=0DO I=1,40IF(NAME(I:I.NE. NUM=NUM+1ENDDOOPEN(1,FILE=NAME(1

6、:NUM, STATUS=UNKNOWNOPEN(2,FILE=NAME(1:NUM/.RES,STATUS=UNKNOWNOPEN(3,FILE=NAME(1:NUM/.ERO,STATUS=UNKNOWNRETURNENDCCSUBROUTINE CLOSEFCLOSE(1CLOSE(2CLOSE(3RETURNENDCSUBROUTINE DA TAIN(NP,NE,NF,ND,NDF,NPF,NM,NR,NCF,& IME,INAE,IIT,ILMT,IMAXA,& IX,IY,IZ,IRR,IAE,IPF,ICKKIMPLICIT REAL*8 (A-H,O-ZIMPLICIT IN

7、TEGER*4 (I-NCOMMON /AT/A(18000000COMMON /IAT/IA(2000000READ(1,*NP,NE,NM,NR,NCFWRITE(2,701NP,NE,NM,NR,NCF701 FORMA T(/1X,#OUTPUT OF ORIGINAL INPUT INFORMA TION# & /5X,Number of joints JOINTS=,I5& /5X,Number of elements ELEMENTS=,I5 & /5X,Number of material property groups PROPERTY TYPES=,I5& /5X,Numb

8、er of restrained joints RESTRAINTS=,I5& /5X,Number of concentrative forced joints NCF=,I5C-FORM POINTER-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

9、(IX+I,A(IY+I,A(IZ+I,I=1,NP714 FORMA T(/5X,GENERA TED JOINT COORDINA TES DATA & /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(

10、IME+2*(I-1+1,IA(IME+2*(I-1+2,& IA(INAE+I,I=1,NE606 FORMA T(/5X,ELEMENT MA TERAIL PROPERTIES DATA & /2X,NO.,10X, E,10X,Ax& /(2X,I3,2(1X,E11.5607 FORMA T(/5X,TRUSS ELEMENT DEFINITION DATA& /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,60

11、8(A(IRR+2*(I-1+1, A(IRR+2*(I-1+2,I=1,NR608 FORMA T(/5X,JOINT RESTRAINTS DATA& /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,NCF609 FORMA T(/5X,CONCENTRA TIVE FORCED JOINTS DATA & /2X, JOINT,10X,Fx, 10X,Fy,10X,Fz& /(2X,F7.0,3(1X

12、,E12.6RETURNENDCSUBROUTINE MKFORCE(NP,NF,NPF,NCF,NN,IT,PF,PP,FTOOLIMPLICIT REAL*8 (A-H,O-ZIMPLICIT INTEGER*4 (I-NDIMENSION IT(NF,NP,PF(4,NCF,PP(NPF, FTOOL(NPFPP=0;FTOOL=0DO I=1,NCFNOD=PF(1,IDO J=1,NFPP(NF*(NOD-1+J=PF(J+1,IENDDODO I=1,NPDO J=1,NFLAB=IT(J,IIF(LAB.GT.0.AND.LAB.LE.NN THENFTOOL(LAB=PP(NF

13、*(I-1+JENDIFENDDOENDDORETURNENDCCSUBROUTINE DATAOUT(NP,NE,NPF,DIST,FF,SG,SMIMPLICIT REAL*8 (A-H,O-ZIMPLICIT INTEGER*4 (I-NDIMENSION DIST(NPF,FF(NPF,SG(NE,SM(NEWRITE(2,715(I,(DIST(3*(I-1+J,J=1,3,I=1,NP715 FORMA T(/5X,SOLVED JOINT DISPLACEMENTS DATA & /1X, JOINT ,3X,8X,Dx,12X,Dy,12X,Dz& /(4X,I5,3X,3(2

14、X,E12.6WRITE(2,716(IE,SG(IE,SM(IE,IE=1,NE716 FORMA T(/5X,SOLVED ELEMENT INTERNAL FORCE DATA & /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,NP717 FORMA T(/5X,SOLVED JOINT REACTION DATA& /1X, JOINT ,3X,8X,Rx,12X,Ry,12X,Rz& /(4X,I5,3X,3(2X,f12.4RETURNENDCC

15、SUBROUTINE FLMT(NP,NE,NN,NN1,NR,RR,ND,NF,NDF,ME,IT,LMTIMPLICIT REAL*8(A-H,O-ZIMPLICIT INTEGER*4(I-NC This program forms the joint&element numbering matrix IT&LMTDIMENSION IT(NF,NP,LMT(NDF,NE,ME(ND,NE,RR(2,NRNN=0;NN1=0;IT=0;LMT=0N=0DO I=1,NPC=0DO K=1,NRKR=RR(1,KIF(KR.EQ.I C=RR(2,KENDDONC=CC=C-NCDO J=

16、1,NFC=C*10.0L=C+0.1C=C-LIF(L.EQ.0THENN=N+1IT(J,I=NELSEIT(J,I=0ENDIFENDDOENDDONN=NNN1=NN+1DO IE=1,NEDO I=1,NDNI=ME(I,IEDO J=1,NFLMT(I-1*NF+J,IE=IT(J,NIENDDOENDDOENDDORETURNENDCSUBROUTINE FMAXA(NN1,NE,LMT,MAXA,NWK,NPF,NDF, & ICKK,IDIST,IFTOOL,IFF,IPP,ISG,ISMC This program forms the MDE address matrix

17、MAXA of KIMPLICIT REAL*8 (A-H,O-ZIMPLICIT INTEGER*4 (I-NDIMENSION MAXA(NPF,LMT(NDF,NEMAXA=0;NWK=0MAXA(1=1DO I=2,NN1IP=I-1IG=IPDO IE=1,NEDO J=1,NDFIF(LMT(J,IE.EQ.IP THENDO K=1,NDFIF(LMT(K,IE.GT.0.AND.LMT(K,IE.LE.IG IG=LMT(K,IEENDDOEND IFENDDOENDDOMAXA(I=MAXA(I-1+IP-IG+1ENDDONWK=MAXA(NN1-1IDIST=ICKK+N

18、WKIFTOOL=IDIST+NPFIFF=IFTOOL+NPFIPP=IFF+NPFISG=IPP+NPFISM=ISG+NERETURNENDCSUBROUTINE CONKB(NP,NE,NM,NWK,ME,X,Y,Z,AE,NAE, & LMT,MAXA,CKK,NN1IMPLICIT REAL*8 (A-H,O-ZIMPLICIT INTEGER*4 (I-NDIMENSION CKK(NWK,X(NP,Y(NP,Z(NP,AE(2,NM, & NAE(NE,LMT(6,NE,ME(2,NE,MAXA(NN1,& AKE(2,2,T(2,6,TT(6,2,AK(6,2,TAK(6,6

19、CKK=0DO 10 IE=1,NETAK=0CALL FKE(NP,NE,NM,IE,X,Y,Z,ME,NAE,AE,AKECALL FT(IE,NP,NE,X,Y,Z,ME,TCALL MAT(2,6,T,TTAK=MA TMUL(TT,AKETAK=MA TMUL(AK,TDO 220 I=1,6DO 220 J=1,6NI=LMT(I,IENJ=LMT(J,IEIF(NJ-NI.GE.0.AND.NI*NJ.GT.0 THENIJ=MAXA(NJ+NJ-NICKK(IJ=CKK(IJ+TAK(I,JENDIF220 CONTINUE10 CONTINUERETURNENDCSUBROU

20、TINE LDLT(A,MAXA,NN,ISH,IOUT,NWK,NNMIMPLICIT REAL*8(A-H,O-ZIMPLICIT INTEGER*4 (I-NDIMENSION A(NWK,MAXA(NNMIF(NN.EQ.1 RETURNDO 200 N=1,NNKN=MAXA(NKL=KN+1KU=MAXA(N+1-1KH=KU-KLIF(KH304,240,210210 K=N-KHIC=0KLT=KUDO 260 J=1,KHKLT=KLT-1IC=IC+1KI=MAXA(KND=MAXA(K+1-KI-1IF(ND 260,260,270270 KK=MIN0(IC,NDC=0

21、.0DO 280 L=1,KK280 C=C+A(KI+L*A(KLT+LA(KLT=A(KLT-C260 K=K+1240 K=NB=0.0DO 300 KK=KL,KUK=K-1KI=MAXA(KC=A(KK/A(KIIF(ABS(C.LT.1.0E+07 GOTO 290WRITE(IOUT,2010 N,CSTOP290 B=B+C*A(KK300 A(KK=CA(KN=A(KN-B304 IF(A(KN 310,310,200310 IF(ISH.EQ.0 GOTO 320IF(A(KN.EQ.0.0 A(KN=-1.0E-16GOTO 200320 WRITE(IOUT,2000

22、N,A(KNSTOP200 CONTINUERETURN2000 FORMA T(/ Stop-stiffness matrix not positive + definite,/,nonpositive pivot for equation,+ I4,/, pivot =,E20.102010 FORMAT(/, Stop-sturm sequence check failed + because of multiplier growth for column+ number,I4,/, Multiplier = ,E20.8ENDCCSUBROUTINE RESOLVE(A,V,MAXA,

23、NN,NWK,NNMIMPLICIT REAL*8(A-H,O-ZIMPLICIT INTEGER*4 (I-NDIMENSION A(NWK,V(NN,1,MAXA(NNMNIP=1DO IP=1,NIPDO 400 N=1,NNKL=MAXA(N+1KU=MAXA(N+1-1IF(KU-KL 400,410,410410 K=NC=0.0DO 420 KK=KL,KUK=K-1420 C=C+A(KK*V(K,IPV(N,IP=V(N,IP-C400 CONTINUEDO 480 N=1,NNK=MAXA(N480 V(N,IP=V(N,IP/A(KIF(NN.EQ.1RETURNN=NN

24、DO 500 L=2,NNKL=MAXA(N+1KU=MAXA(N+1-1IF(KU-KL 500,510,510510 K=NDO 520 KK=KL,KUK=K-1520 V(K,IP=V(K,IP-A(KK*V(N,IP500 N=N-1ENDDORETURNENDCCccSUBROUTINE DISPLS(NP,NE,NF,NPF,NM,NN,IT,FTOOL, & DIST,AE,ME,NAE,X,Y,Z,PP,FF,SG,SMIMPLICIT REAL*8 (A-H,O-ZIMPLICIT INTEGER*4 (I-NDIMENSION IT(NF,NP,DIST(NPF,FTOO

25、L(NPF,T(2,6, & TT(6,2,AE(2,NM,ME(2,NE,NAE(NE,UE(6,U(2,& AKE(2,2,FE1(2,FE(6,FF(NPF,X(NP,Y(NP,Z(NP,& PP(NPF,SG(NE,SM(NESG=0;SM=0;FF=0DO I=1,NPDO J=1,NFLAB=IT(J,IIF(LAB.EQ.0 THENDIST(NF*(I-1+J=0.0ELSEIF(LAB.GT.0.AND.LAB.LE.NN THENDIST(NF*(I-1+J=FTOOL(LABENDIFENDDOENDDODO IE=1,NEN1=ME(1,IE;N2=ME(2,IEUE=

26、0DO J=1,NFUE(J=DIST(NF*(N1-1+JUE(NF+J=DIST(NF*(N2-1+JENDDOCALL FT(IE,NP,NE,X,Y,Z,ME,TCALL FKE(NP,NE,NM,IE,X,Y,Z,ME,NAE,AE,AKEU=MA TMUL(T,UEFE1=MATMUL(AKE,UCALL MAT(2,6,T,TTFE=MA TMUL(TT,FE1DO J=1,NFFF(NF*(N1-1+J=FF(NF*(N1-1+J+FE(JFF(NF*(N2-1+J=FF(NF*(N2-1+J+FE(NF+JENDDOISW=NAE(IEAO=AE(2,ISWSG(IE=FE1(2SM(IE=FE1(2/AO ENDDO DO I=1,NPF F