fortran語(yǔ)言編寫的有限元結(jié)構(gòu)程序

fortran語(yǔ)言編寫的有限元結(jié)構(gòu)程序

ID:8977078

大小:106.50 KB

頁(yè)數(shù):11頁(yè)

時(shí)間:2018-04-13

fortran語(yǔ)言編寫的有限元結(jié)構(gòu)程序_第1頁(yè)
fortran語(yǔ)言編寫的有限元結(jié)構(gòu)程序_第2頁(yè)
fortran語(yǔ)言編寫的有限元結(jié)構(gòu)程序_第3頁(yè)
fortran語(yǔ)言編寫的有限元結(jié)構(gòu)程序_第4頁(yè)
fortran語(yǔ)言編寫的有限元結(jié)構(gòu)程序_第5頁(yè)
資源描述:

《fortran語(yǔ)言編寫的有限元結(jié)構(gòu)程序》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫(kù)

1、算例一計(jì)算簡(jiǎn)圖及結(jié)果輸出用平面剛架靜力計(jì)算程序下圖結(jié)構(gòu)的內(nèi)力。各桿EA,EI相同。已知:計(jì)算簡(jiǎn)圖如下:(1)輸入原始數(shù)據(jù)控制參數(shù)3,5,8,7,1,2(NE,NJ,N,NW,NPJ,NPF)結(jié)點(diǎn)坐標(biāo)集結(jié)點(diǎn)未知量編號(hào)單元桿端結(jié)點(diǎn)編號(hào)及單元EA、EI結(jié)點(diǎn)荷載非結(jié)點(diǎn)荷載(2)輸出結(jié)果NE=3NJ=5N=8NW=7NPJ=1NPF=2NODEXYXXYYZZ10.00000.000000020.00004.000012330.00004.000012444.00004.000056754.00000.0000008ELEMENTNODE-INODE-JEAEI11

2、20.400000E+070.160000E+052340.400000E+070.160000E+053540.400000E+070.160000E+05CODEPX-PY-PM7.-15.0000ELEMENTINDAQ1.2.2.0000-18.00002.1.4.0000-25.0000NODEUVCETA10.000000E+000.000000E+000.000000E+002-0.221743E-02-0.464619E-04-0.139404E-023-0.221743E-02-0.464619E-040.357876E-024-0.22

3、2472E-02-0.535381E-04-0.298554E-0250.000000E+000.000000E+000.658499E-03ELEMENTNQM1N1=46.4619Q1=10.7119M1=-6.8477N2=-46.4619Q2=7.2881M2=0.00002N1=7.2881Q1=46.4619M1=0.0000N2=-7.2881Q2=53.5381M2=14.15233N1=53.5381Q1=7.2881M1=0.0000N2=-53.5381Q2=-7.2881M2=-29.1523算例二計(jì)算簡(jiǎn)圖及結(jié)果輸出用平面剛架靜力計(jì)

4、算程序下圖結(jié)構(gòu)的內(nèi)力。已知:桁架單元的抗拉剛度為,平面剛架單元的抗拉剛度為已知:,抗彎剛度為。計(jì)算簡(jiǎn)圖如下:(1)輸入原始數(shù)據(jù)控制參數(shù)(NE,NJ,N,NW,NPJ,NPF)結(jié)點(diǎn)坐標(biāo)集結(jié)點(diǎn)未知量編號(hào)單元桿端結(jié)點(diǎn)編號(hào)及單元EA、EI非結(jié)點(diǎn)荷載(2)輸出結(jié)果NE=5NJ=4N=8NW=7NPJ=0NPF=1NODEXYXXYYZZ10.00000.000000124.00000.000023434.0000-3.000056048.00000.0000708ELEMENTNODE-INODE-JEAEI1120.600000E+070.184000E+0622

5、40.600000E+070.184000E+063310.200000E+070.000000E+004320.200000E+070.000000E+005340.200000E+070.000000E+00ELEMENTINDAQ1.1.4.0000-20.0000NODEUVCETA10.000000E+000.000000E+000.312593E-032-0.202759E-04-0.253871E-03-0.144928E-033-0.202759E-04-0.185440E-030.000000E+004-0.405518E-040.000

6、000E+00-0.227378E-04ELEMENTNQM1N1=30.4138Q1=37.1896M1=0.0000N2=-30.4138Q2=42.8104M2=11.24152N1=30.4138Q1=2.8104M1=-11.2415N2=-30.4138Q2=-2.8104M2=0.00003N1=-38.0173Q1=0.0000M1=0.0000N2=38.0173Q2=0.0000M2=0.00004N1=45.6207Q1=0.0000M1=0.0000N2=-45.6207Q2=0.0000M2=0.00005N1=-38.0173Q

7、1=0.0000M1=0.0000N2=38.0173Q2=0.0000M2=0.0000C主程序C(一)輸入原始數(shù)據(jù)DIMENSIONJE(2,100),JN(3,100),JC(6),EA(100),EI(100),X(100),$Y(100),PJ(2,50),PF(4,100)REAL*8KE(6,6),KD(6,6),T(6,6),P(300),KB(200,20),F(6),FO(6),$D(6),BL,SI,CO,S,COPEN(5,FILE='RPF1.TXT')open(6,file='jieguo1.dat',status='new')

8、READ(5,*)NE,NJ,N,NW,NPJ,NPFREAD(5

當(dāng)前文檔最多預(yù)覽五頁(yè),下載文檔查看全文

此文檔下載收益歸作者所有

當(dāng)前文檔最多預(yù)覽五頁(yè),下載文檔查看全文
溫馨提示:
1. 部分包含數(shù)學(xué)公式或PPT動(dòng)畫的文件,查看預(yù)覽時(shí)可能會(huì)顯示錯(cuò)亂或異常,文件下載后無(wú)此問(wèn)題,請(qǐng)放心下載。
2. 本文檔由用戶上傳,版權(quán)歸屬用戶,天天文庫(kù)負(fù)責(zé)整理代發(fā)布。如果您對(duì)本文檔版權(quán)有爭(zhēng)議請(qǐng)及時(shí)聯(lián)系客服。
3. 下載前請(qǐng)仔細(xì)閱讀文檔內(nèi)容,確認(rèn)文檔內(nèi)容符合您的需求后進(jìn)行下載,若出現(xiàn)內(nèi)容與標(biāo)題不符可向本站投訴處理。
4. 下載文檔時(shí)可能由于網(wǎng)絡(luò)波動(dòng)等原因無(wú)法下載或下載錯(cuò)誤,付費(fèi)完成后未能成功下載的用戶請(qǐng)聯(lián)系客服處理。