有限元编程算例fortran范本模板Word文档下载推荐.docx
- 文档编号:20050473
- 上传时间:2023-01-16
- 格式:DOCX
- 页数:18
- 大小:61.13KB
有限元编程算例fortran范本模板Word文档下载推荐.docx
《有限元编程算例fortran范本模板Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《有限元编程算例fortran范本模板Word文档下载推荐.docx(18页珍藏版)》请在冰豆网上搜索。
DO10J=1,6
B(I,J)=0。
10CONTINUE
B(1,1)=-BJ—BM
B(1,3)=BJ
B(1,5)=BM
B(2,2)=—CJ-CM
B(2,4)=CJ
B(2,6)=CM
B(3,1)=B(2,2)
B(3,2)=B(1,1)
B(3,3)=B(2,4)
B(3,4)=B(1,3)
B(3,5)=B(2,6)
Page151
B(3,6)=B(1,5)
DO20I=1,3
DO20J=1,6
B(I,J)=B(I,J)/(2。
0*AE)
20CONTINUE
D(1,1)=EO/(1.0-UN*UN)
D(1,2)=EO*UN/(1。
0-UN*UN)
D(2,1)=D(1,2)
D(2,2)=D(1,1)
D(1,3)=0。
D(2,3)=0。
D(3,1)=0。
D(3,2)=0.0
D(3,3)=EO/(2。
0*(1.0+UN))
DO30I=1,3
DO30J=1,6
S(I,J)=0。
DO30K=1,3
S(I,J)=S(I,J)+D(I,K)*B(K,J)
30CONTINUE
2)GOTO50
DO40I=1,6
DO40J=1,6
EKE(I,J)=0。
DO40K=1,3
**********************************ExchangeBAndS***********************************************
EKE(I,J)=EKE(I,J)+B(K,I)*S(K,J)*AE*TE
40CONTINUE
50CONTINUE
SUBROUTINETOTSTI
COMMON/X1/NJ,NE,NZ,NDD,NPJ,IND,NJ2,EO,UN,GAMA,TE,AE
COMMON/X2/JM(100,3),NZC(50),CJZ(100,2),PJ(100,2),B(3,6),D(3,3),S(3,6),TKZ(200,20),EKE(6,6),P(200)
Page152
DO20I=1,NJ2
DO20J=1,NDD
TKZ(I,J)=0.0
*************NotUnderstanded*****************************
DO30MEO=1,NE
CALLELEST(MEO,3)
DO30I=1,3
DO30II=1,2
LH=2*(I—1)+II
LDH=2*(JM(MEO,I)-1)+II
DO30J=1,3
DO30JJ=1,2
L=2*(J—1)+JJ
LZ=2*(JM(MEO,J)-1)+JJ
LD=LZ-LDH+1
IF(LD。
0)GOTO30
TKZ(LDH,LD)=TKZ(LDH,LD)+EKE(LH,L)
SUBROUTINELOAD
COMMON/X1/NJ,NE,NZ,NDD,NPJ,IND,NJ2,EO,UN,GAMA,TE,AE
COMMON/X2/JM(100,3),NZC(50),CJZ(100,2),PJ(100,2),B(3,6),D(3,3),S(3,6),TKZ(200,20),EKE(6,6),P(200)
DO10I=1,NJ2
P(I)=0。
IF(NPJ.EQ.0)GOTO30
DO20I=1,NPJ
I1=I+1
J=IFIX(PJ(I1,2))
P(J)=PJ(I1,1)
30IF(GAMA.LE.0.0)GOTO50
Page153
DO40MEO=1,NE
CALLELEST(MEO,1)
PE=—GAMA*AE*TE/3。
IE=JM(MEO,1)
JE=JM(MEO,2)
ME=JM(MEO,3)
P(2*IE)=P(2*IE)+PE
P(2*JE)=P(2*JE)+PE
P(2*ME)=P(2*ME)+PE
SUBROUTINESUPPOR
COMMON/X1/NJ,NE,NZ,NDD,NPJ,IND,NJ2,EO,UN,GAMA,TE,AE
COMMON/X2/JM(100,3),NZC(50),CJZ(100,2),PJ(100,2),B(3,6),D(3,3),S(3,6),TKZ(200,20),EKE(6,6),P(200)
DO60I=1,NZ
MZ=NZC(I)
TKZ(MZ,1)=1.0
DO10J=2,NDD
TKZ(MZ,J)=0。
10CONTINUE
IF(MZ-NDD)20,20,30
20JO=MZ
GOTO40
30JO=NDD
40DO50J=2,JO
J1=MZ—J
TKZ(J1+1,J)=0.0
50CONTINUE
P(MZ)=0。
60CONTINUE
Page154
SUBROUTINESOLVEQ
COMMON/X1/NJ,NE,NZ,NDD,NPJ,IND,NJ2,EO,UN,GAMA,TE,AE
COMMON/X2/JM(100,3),NZC(50),CJZ(100,2),PJ(100,2),B(3,6),D(3,3),S(3,6),TKZ(200,20),EKE(6,6),P(200)
NJ1=NJ2—1
DO50K=1,NJ1
IF(NJ2-K—NDD+1)10,10,20
10IM=NJ2
GOTO30
20IM=K+NDD-1
30K1=K+1
DO50I=K1,IM
L=I—K+1
C=TKZ(K,L)/TKZ(K,1)
LD1=NDD—L+1
DO40J=1,LD1
M=J+I-K
TKZ(I,J)=TKZ(I,J)-C*TKZ(K,M)
40CONTINUE
P(I)=P(I)-C*P(K)
P(NJ2)=P(NJ2)/TKZ(NJ2,1)
DO100I1=1,NJ1
I=NJ2-I1
************************************************************************下面一行可能出错
IF(NDD-NJ2+I—1)60,60,70
60JO=NDD
GOTO80
70JO=NJ2—I+1
80DO90J=2,JO
LH=J+I—1
P(I)=P(I)—TKZ(I,J)*P(LH)
90CONTINUE
P(I)=P(I)/TKZ(I,1)
100CONTINUE
Page155
WRITE(6,110)(I,P(2*I-1),P(2*I),I=1,NJ)
************************************************************************************
110FORMAT(2X,3HJD=,3X,2HU=,12X,2HV=/(I4,3X,F16.7,3X,F16。
7))
SUBROUTINESTRESS
COMMON/X1/NJ,NE,NZ,NDD,NPJ,IND,NJ2,EO,UN,GAMA,TE,AE
COMMON/X2/JM(100,3),NZC(50),CJZ(100,2),PJ(100,2),B(3,6),D(3,3),S(3,6),TKZ(200,20),EKE(6,6),P(200)
DIMENSIONWY(6),YL(3)
DO60MEO=1,NE
CALLELEST(MEO,2)
DO10J=1,2
LH=2*(I—1)+J
LDH=2*(JM(MEO,I)—1)+J
WY(LH)=P(LDH)
DO20I=1,3
YL(I)=0。
DO20J=1,6
YL(I)=YL(I)+S(I,J)*WY(J)
SIGX=YL
(1)
SIGY=YL
(2)
TOXY=YL(3)
PYL=(SIGX+SIGY)/2。
SIG=(SIGX—SIGY)**2/4.0+TOXY*TOXY
RYL=SQRT(SIG)
SIG1=PYL+RYL
SIG2=PYL—RYL
IF(SIGY.EQ.SIG2)GOTO30
CETA1=TOXY/(SIGY—SIG2)
CETA=90。
0—57。
29578*ATAN(CETA1)
GOTO40
Page156
30CETA=0.0
40WRITE(6,50)MEO,SIGX,SIGY,TOXY,SIG1,SIG2,CETA
50FORMAT(4X,2HE=,I3/2X,3HSX=,F11。
3,3X,3HSY=,F11.3,3X,4HTAU=,F11。
3/2X,3HS1=,F11。
3,3X,3HS2=,F11.3,3X,4HCET=,F11。
3)
50FORMAT(4X,2HE=,I3/2X,3HSX=,Fll。
3,3X,3HSY=,F11。
3,3X,4HTAU=,F11。
3/2X,3HSl=,Fll。
3,3X,3HS2=,F11.3,3X,4HCET=,F11.3)
60CONTINUE
输入文件为datain
28,36,9,10,4,0
1,0.17,0,1
1,5,2
2,5,6
2,6,3
3,6,7
3,7,4
4,7,8
5,9,6
6,9,10
6,10,7
7,10,11
7,11,8
8,11,12
9,13,10
10,13,14
10,14,11
11,14,15
11,15,12
12,15,16
13,17,14
14,17,18
14,18,15
15,18,19
15,19,16
16,19,20
17,21,18
18,21,22
18,22,19
19,22,23
19,23,20
20,23,24
21,25,22
22,25,26
22,26,23
23,26,27
23,27,24
24,27,28
0,6
1,6
2,6
3,6
0,5
1,5
2,5
3,5
0,4
1,4
2,4
3,4
0,3
1,3
2,3
3,3
0,2
1,2
2,2
3,2
0,1
1,1
2,1
3,1
0,0
1,0
2,0
3,0
7,15,23,31,39,47,49,50,55
-5E4,2
-10E4,4
—10E4,6
-5E4,8
输出结果为:
DATAOUT
NOXY
10.006。
00
21.006.00
32。
006.00
43.006。
50.005。
61。
005。
72。
005.00
83.005。
90。
004.00
101.004.00
112。
004。
123。
130.003。
141。
003。
152.003。
163。
170。
002。
181.002。
192.002.00
203。
002.00
210.001。
221。
001.00
232.001。
243。
001。
250.000。
261.000.00
272.000。
283.000.00
JD=U=V=
1-29766.873-1173917。
750
2-14003。
185-1174018。
875
3—3753。
270—1179518。
125
40.000—1181719.750
5-26382.471—1072681.500
6-10746。
993-1073615.000
7-2064。
593—1082360。
80.000—1085873。
250
9—13536.995—964010.125
103372.794-970055.125
117268。
415—989269。
120。
000-998401。
812
137816。
581-835383.438
1427176.234-861713。
938
1522063。
230—905726。
160。
000—927165.188
1729514。
479—665602。
1853419.637-747340。
438
1934876.832—839806。
200。
000—881219.125
2129580.273-416288。
719
2252944.918-632601。
2317504。
195-803765.688
240.000—859481。
250.0000.000
26—120102。
820-583505。
375
27—76202。
375-787347。
188
280。
000-829170。
E=1
SX=—1489.530SY=-101489。
383TAU=-1489.531
S1=-1467.348S2=-101511。
562CET=179.147
E=2
SX=—1475。
844SY=—100654.875TAU=—1790。
500
S1=—1443.531S2=—100687.188CET=178。
966
E=3
SX=-7021.670SY=—101597.672TAU=—3741。
688
S1=—6873.875S2=-101745.469CET=177.738
E=4
SX=-8067。
500SY=—98528.750TAU=—4459.156
S1=—7848.227S2=—98748。
023CET=177。
185
E=5
SX=—13143。
328SY=—99391.750TAU=-1662.500
S1=—13111。
293S2=—99423。
781CET=178。
896
E=6
SX=—14652。
781SY=-98337.500TAU=-1501.062
S1=-14625。
867S2=—98364。
414CET=178.973
E=7
SX=-2923。
122SY=-109168。
297TAU=—5888。
469
S1=—2597.762S2=—109493。
656CET=176.837
E=8
SX=-716.078SY=—103681。
562TAU=-8617。
406
S1=0。
148S2=—104397.789CET=175。
249
E=9
SX=—9188。
316SY=-105121.867TAU=-9771.594
S1=—8203.125S2=-106107。
062CET=174.243
E=10
SX=—12285。
000SY=-95180.250TAU=—12199.594
S1=-10526.887S2=-96938.359CET=171。
799
E=11
SX=—14170.516SY=—95500。
750TAU=—5489。
531
S1=-13801。
664S2=—95869。
602CET=176。
156
E=12
SX=-22797。
406SY=—91347.000TAU=—3902.844
S1=-22575。
914S2=—91568。
492CET=176.752
E=13
SX=—5104。
269SY=—129494。
438TAU=-11708.750
S1=-4011.727S2=—130586。
977CET=174。
669
E=14
SX=969。
672SY=—108176。
375TAU=-21424.750
S1=5024.582S2=—112231。
281CET=169。
283
E=15
SX=—14954.572SY=-110883。
469TAU=-18383。
S1=-11552.273S2=—114285.766CET=169.515
E=16
SX=-19890.141SY=—86924。
312TAU=-25131。
S1=-11514。
844S2=-95299.609CET=161。
569
E=17
SX=-22109.688SY=—87301。
625TAU=—10225.406
S1=-20543。
453S2=—88867。
859CET=171。
292
E=18
SX=-35190.453SY=—77219.000TAU=-9162。
000
S1=—33280.023S2=—79129。
430CET=168。
222
E=19
SX=-9785。
850SY=—171444。
172TAU=-20524。
969
S1=-7220.594S2=—174009。
422CET=172。
876
E=20
SX=4594.438SY=-113592。
375TAU=—46145。
S1=20477.398S2=—129475.336CET=161.007
E=21
SX=—25287。
307SY=—118672.312TAU=—30023。
S1=-16467.512S2=-127492。
109CET=163。
629
E=22
SX=-30634。
422SY=-71127.188TAU=—44991。
S1=—1543.715S2=—100217.891CET=147。
114
E=23
SX=-34259.609SY=-71743.438TAU=-14637。
906
S1=—29220。
699S2=—76782.344CET=161。
005
E=24
SX=—43958。
047SY=—53418.938TAU=-17697。
562
S1=-30369。
627S2=-67007.359CET=142.482
E=25
SX=-19028。
160SY=-252549。
000TAU=-34958。
S1=-13907.055S2=-257670.094CET=171.666
E=26
SX=3973.812SY=—114063.750TAU=—92238。
344
S1=54459。
047S2=—164548.984CET=151。
307
E=27
SX=-39180。
809SY=-12
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 有限元 编程 fortran 范本 模板