空間桁架結(jié)構(gòu)程序設(shè)計(jì)(Fortran).doc

《空間桁架結(jié)構(gòu)程序設(shè)計(jì)(Fortran).doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、空間桁架靜力分析程序及算例1、變量及數(shù)組說明輸入數(shù)據(jù)控制數(shù)據(jù)NF單個(gè)節(jié)點(diǎn)的自由度數(shù)NP結(jié)構(gòu)離散節(jié)點(diǎn)的總數(shù)NE結(jié)構(gòu)離散單元的總數(shù)NM結(jié)構(gòu)中單元不同的特征數(shù)類的總數(shù)NR結(jié)構(gòu)受約束節(jié)點(diǎn)的總數(shù)NCF結(jié)構(gòu)受外荷載作用的節(jié)點(diǎn)總數(shù)ND一個(gè)單元的節(jié)點(diǎn)總數(shù)幾何數(shù)據(jù)X(NP)節(jié)點(diǎn)X坐標(biāo)數(shù)組Y(NP)節(jié)點(diǎn)Y坐標(biāo)數(shù)組Z(NP)節(jié)點(diǎn)Z坐標(biāo)數(shù)組ME(ND,NE)單元節(jié)點(diǎn)信息存儲(chǔ)矩陣ME(1,NE)存儲(chǔ)桿件始端節(jié)點(diǎn)號(hào)ME(2,NE)儲(chǔ)存桿件末端節(jié)點(diǎn)號(hào)RR(2,NR)結(jié)構(gòu)約束信息矩陣RR(1,NR)存放受有約束的節(jié)點(diǎn)號(hào)RR(2,NR)存放節(jié)點(diǎn)位移約束情況單元特征數(shù)據(jù)AE(2,IN)單元特征數(shù)類數(shù)組AE(

2、1,IN)單元的彈性模量AE(2,IN)單元的橫截面面積NAE(NE)單元特征類信息存儲(chǔ)數(shù)組荷載數(shù)據(jù)PF(4,NCF)外荷載信息數(shù)組PF(1,NCF)存放外荷載作用的節(jié)點(diǎn)號(hào)PF(2,NCF)存放X方向的外荷載PF(3,NCF)存放Y方向的外荷載PF(4,NCF)存放Z方向的外荷載輸出數(shù)據(jù)位移DIST(NPF)節(jié)點(diǎn)位移數(shù)組DIST(NF*I-2)存放I節(jié)點(diǎn)X方向的位移DIST(NF*I-1)存放I節(jié)點(diǎn)Y方向的位移DIST(NF*I)存放I節(jié)點(diǎn)Z方向的位移力SG(NE)單元力數(shù)組SM(NE)單元截面應(yīng)力數(shù)組FL(NF*NR)支座反力數(shù)組FL(NF*I-2)存放受約束的I節(jié)點(diǎn)X

3、方向的反力FL(NF*I-1)存放受約束的I節(jié)點(diǎn)Y方向的反力FL(NF*I)存放受約束的I節(jié)點(diǎn)Z方向的反力中間變量NPF=NF*NP二維總剛度矩陣的最大行數(shù)NDF=ND*NF一個(gè)單元的自由度總數(shù)（2*3=6）IN單元特征類總數(shù)AKE(2,2)單元在局部坐標(biāo)系中的剛度局矩陣BL桿件單元長度T(2,6)坐標(biāo)轉(zhuǎn)換矩陣TAK(6,6)單元在總體坐標(biāo)系中的剛度矩陣IT(NF,NP)節(jié)點(diǎn)聯(lián)系數(shù)組LMT(NDF,NE)單元聯(lián)系數(shù)組MAXA(NPF)結(jié)構(gòu)二維總剛度矩陣主對(duì)角元地址數(shù)組NWK結(jié)構(gòu)一維總剛度矩陣的總?cè)萘緾KK(NWK)結(jié)構(gòu)一維總剛度矩陣NN結(jié)構(gòu)矩陣方程的方程總數(shù)（去掉約束）N

4、NMNNM=NN+1V(NN)已知節(jié)點(diǎn)荷載列陣數(shù)組，回代完成后為存放結(jié)構(gòu)位移PP(NPF)所有節(jié)點(diǎn)荷載列陣數(shù)組2、空間桁架結(jié)構(gòu)有限元分析程序源代碼!主程序（讀入文件，調(diào)用總計(jì)算程序，輸出結(jié)果）CHARACTERIDFUT*20,OUTFUT*20WRITE(*,*)'InputDataFilename:'READ(*,*)IDFUTOPEN(11,FILE=IDFUT,STATUS='OLD')WRITE(*,*)'OutputFilename:'READ(*,*)OUTFUTOPEN(12,FILE=OUTFUT,STATUS='UNKNOWN')WRITE(12,*)

5、'*****************************************'WRITE(12,*)'*ProgramforAnalysisofSpaceTrusses*'WRITE(12,*)'*SchoolofCivilEngineeringCSU*'WRITE(12,*)'*2012.6.25DesignedByMuZhaoxiang*'WRITE(12,*)'*****************************************'WRITE(12,*)''WRITE(12,*)'*********************************

6、********'WRITE(12,*)'*************TheInputData****************'WRITE(12,*)'*****************************************'WRITE(12,100)READ(11,*)NF,NP,NE,NM,NR,NCF,NDWRITE(12,110)NF,NP,NE,NM,NR,NCF,ND100FORMAT(6X,'TheGeneralInformation'/2X,'NF',5X,'NP',5X,'NE',5X,'NM',5X,'NR',&5X,'NCF',5X,'ND'

7、)110FORMAT(2X,I2,6I7)NPF=NF*NPNDF=ND*NFCALLANALYSE(NF,NP,NE,NM,NR,NCF,ND,NPF,NDF)END!********************************************************************!總計(jì)算程序SUBROUTINEANALYSE(NF,NP,NE,NM,NR,NCF,ND,NPF,NDF)DIMENSIONX(NP),Y(NP),Z(NP),MM(NE),ME(ND,NE),IT(NF,NP),RR(ND