1、中南大學(xué)鐵道校區(qū)(土木2012級(jí)研究生結(jié)構(gòu)分析程序設(shè)計(jì)基礎(chǔ)試題(考試方式:開(kāi)卷,交卷時(shí)間:第13周1.如圖示平面應(yīng)力懸臂梁模型2423222120191817161514131211987654321012019181716151413121110987654321圖1 平面應(yīng)力懸臂梁模型試寫(xiě)出其節(jié)點(diǎn)聯(lián)系數(shù)組IT 、單元聯(lián)系數(shù)組LMT 和對(duì)角元地址數(shù)組MAXA 。2、求解線性方程組AX=B 式中的A為對(duì)稱(chēng)正定矩陣,其值為對(duì)稱(chēng)B矩陣為B=3.4,5.5,12.3,13.5,0.6,1.4,5.6T要求寫(xiě)出矩陣分解和回代求解過(guò)程。3.在計(jì)算機(jī)上調(diào)試空間桁架分析程序TRUSS ,并用兩個(gè)算例進(jìn)行考證

2、,要求寫(xiě)出相關(guān)計(jì)算結(jié)果。應(yīng)大家要求,程序已經(jīng)調(diào)試,附錄于下: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=0WRITE(*,*INPUT FILE NAME?;READ(*,*NAME;CALL OPENF(NAMECA

3、LL 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,IPP,ISG,ISMCALL CONKB(NP,NE,NM,NWK,IA(IME+1,A(IX+1,& A(IY+

4、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,NE,NF,NPF,NM,NN,IA(IIT+1,A(IFTOOL+1,& A(IDIST+1,A(IAE+1,IA

5、(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:NUM, STATUS=UNKNOWNOPEN(2,FILE=NAME(1:NUM/.RES,STATUS=UN

6、KNOWNOPEN(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 INTEGER*4 (I-NCOMMON /AT/A(18000000COMMON /IAT/IA(2000000RE

7、AD(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,Number of restrained joints RESTRAINTS=,I5& /5X,Number of con

8、centrative 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(IX+I,A(IY+I,A(IZ+I,I=1,NP714 FORMA T(/5X,GENERA TED JOIN

9、T 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(IME+2*(I-1+1,IA(IME+2*(I-1+2,& IA(INAE+I,I=1,NE606 FORMA

10、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,608(A(IRR+2*(I-1+1, A(IRR+2*(I-1+2,I=1,NR608 FORMA T(/5X,JO

11、INT 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,E12.6RETURNENDCSUBROUTINE MKFORCE(NP,NF,NPF,NCF,NN,IT,PF

12、,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,IFTOOL(LAB=PP(NF*(I-1+JENDIFENDDOENDDORETURNENDCCSUBROUTINE DATAOUT(NP,NE,NPF,DIST,FF,SG,SMIMPLICIT REA

13、L*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(2X,E12.6WRITE(2,716(IE,SG(IE,SM(IE,IE=1,NE716 FORMA T(/5X,SOLVED ELEMENT INTERNAL FORCE

14、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.4RETURNENDCCSUBROUTINE FLMT(NP,NE,NN,NN1,NR,RR,ND,NF,NDF,ME,IT,LMTIMPLICIT REAL*8(A-H,O-ZIMPLICIT I

15、NTEGER*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,KENDDONC=CC=C-NCDO J=1,NFC=C*10.0L=C+0.1C=C-LN=N+1IT(J,I=NELSEIT(J,I=0ENDIFENDDOENDDONN=NNN1=NN+1DO IE=1,NEDO I=1,NDNI=ME(I,IED

16、O 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 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

17、-1IG=IPDO IE=1,NEDO J=1,NDFDO K=1,NDFENDDOEND IFENDDOENDDOMAXA(I=MAXA(I-1+IP-IG+1ENDDONWK=MAXA(NN1-1IDIST=ICKK+NWKIFTOOL=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-NDI

18、MENSION 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,6CKK=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,IEIJ

19、=MAXA(NJ+NJ-NICKK(IJ=CKK(IJ+TAK(I,JENDIF220 CONTINUE10 CONTINUERETURNENDCSUBROUTINE LDLT(A,MAXA,NN,ISH,IOUT,NWK,NNMIMPLICIT REAL*8(A-H,O-ZIMPLICIT INTEGER*4 (I-NDIMENSION A(NWK,MAXA(NNMDO 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+1

20、KI=MAXA(KND=MAXA(K+1-KI-1IF(ND 260,260,270270 KK=MIN0(IC,NDC=0.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(KIWRITE(IOUT,2010 N,CSTOP290 B=B+C*A(KK300 A(KK=CA(KN=A(KN-B304 IF(A(KN 310,310,200GOTO 200320 WRITE(IOUT,2000 N,A(KNSTOP200 CONT

21、INUERETURN2000 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,NN,NWK,NNMIMPLICIT

22、 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(KN=NNDO 500 L=2,NNKL=MAXA(N+1KU=MAXA(N+

23、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,FTOOL(NPF,T(2,6, & TT(6,2,AE(2,NM,ME(2

24、,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,IDIST(NF*(I-1+J=0.0DIST(NF*(I-1+J=FTOOL(LABENDIFENDDOENDDODO IE=1,NEN1=ME(1,IE;N2=ME(2,IEUE=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 FF(I=FF(I-PP(I ENDDO RETURN END C SUBRO