有限元编程算例Word下载.docx
- 文档编号:20949424
- 上传时间:2023-01-26
- 格式:DOCX
- 页数:20
- 大小:146.38KB
有限元编程算例Word下载.docx
《有限元编程算例Word下载.docx》由会员分享,可在线阅读,更多相关《有限元编程算例Word下载.docx(20页珍藏版)》请在冰豆网上搜索。
CJ=CJZ(IE,1)-CJZ(ME,1)
BJ=CJZ(ME,2)-CJZ(IE,2)
AE=(BJ*CM-BM*CJ)/2.0
IF(IASK.LE.1)GOTO50
DO10I=1,3
DO10J=1,6
B(I,J)=0.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.0
D(2,3)=0.0
D(3,1)=0.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.0
DO30K=1,3
S(I,J)=S(I,J)+D(I,K)*B(K,J)
30CONTINUE
IF(IASK.LE.2)GOTO50
DO40I=1,6
DO40J=1,6
EKE(I,J)=0.0
DO40K=1,3
And
**********************************ExchangeB
S***********************************************
EKE(I,J)=EKE(I,J)+B(K,I)*S(K,J)*AE*TE
40CONTINUE
50CONTINUE
SUBROUTINETOTSTI
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)
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.LE.0)GOTO30
TKZ(LDH,LD)=TKZ(LDH,LD)+EKE(LH,L)
SUBROUTINELOAD
DO10I=1,NJ2
P(I)=0.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.0
P(2*IE)=P(2*IE)+PE
P(2*JE)=P(2*JE)+PE
P(2*ME)=P(2*ME)+PE
40CONTINUE50CONTINUE
SUBROUTINESUPPOR
DO60I=1,NZ
MZ=NZC(I)
TKZ(MZ,1)=1.0
DO10J=2,NDD
TKZ(MZ,J)=0.0
IF(MZ-NDD)20,20,30
20JO=MZ
GOTO40
30JO=NDD
40DO50J=2,JO
J1=MZ-J
TKZ(J1+1,J)=0.0
P(MZ)=0.0
60CONTINUE
Page154
SUBROUTINESOLVEQ
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,LD1M=J+I-KTKZ(I,J)=TKZ(I,J)-C*TKZ(K,M)
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,7060JO=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
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)
YL(I)=0.0
YL(I)=YL(I)+S(I,J)*WY(J)
SIGX=YL
(1)
SIGY=YL
(2)
TOXY=YL(3)
PYL=(SIGX+SIGY)/2.0
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)
Page156
30CETA=0.0
40WRITE(6,50)MEO,SIGX,SIGY,TOXY,SIG1,SIG2,CETA
50
FORMAT(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)
FORMAT(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)
输入文件为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
NO
X
Y
1
0.00
6.00
2
1.00
3
2.00
4
3.00
5
5.00
6
7
8
9
4.00
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
JD=U=V=
1-29766.873-1173917.750
-14003.185
-1174018.875
-3753.270
-1179518.125
0.000
-1181719.750
-26382.471
-1072681.500
-10746.993
-1073615.000
-2064.593
-1082360.750
-1085873.250
-13536.995
-964010.125
3372.794
-970055.125
7268.415
-989269.125
-998401.812
7816.581
-835383.438
27176.234
-861713.938
22063.230
-905726.125
-927165.188
29514.479
-665602.875
53419.637
-747340.438
34876.832
-839806.812
-881219.125
29580.273
-416288.719
52944.918
-632601.125
17504.195
-803765.688
-859481.938
-120102.820
-583505.375
-76202.375
-787347.188
-829170.812
E=1
SX=
-1489.530
SY=-101489.383
TAU=
-1489.531
S1=
-1467.348
S2=-101511.562
CET=
179.147
E=
-1475.844
SY=-100654.875
-1790.500
-1443.531
S2=-100687.188
178.966
-7021.670
SY=-101597.672
-3741.688
-6873.875
S2=-101745.469
177.738
-8067.500
SY=-98528.750
-4459.156
-7848.227
S2=-98748.023
177.185
SX=-
13143.328
SY=-99391.750
-1662.500
S1=-13111.293
S2=-99423.781
178.896
SX=-14652.781
SY=-98337.500
-1501.062
S1=-14625.867
S2=-98364.414
178.973
E=7
SX=-2923.122
SY=-109168.297
-5888.469
S1=-2597.762
S2=-109493.656
176.837
E=8
SX=-716.078
SY=-103681.562
-8617.406
S1=0.148
S2=-104397.789
175.249
E=9
SX=-9188.316
SY=-105121.867
-9771.594
S1=-8203.125
S2=-106107.062
174.243
E=10
SX=-12285.000
SY=-95180.250
TAU=-12199.594
S1=-10526.887
S2=-96938.359
171.799
E=11
SX=-14170.516
SY=-95500.750
-5489.531
S1=-13801.664
S2=-95869.602
176.156
E=12
SX=-22797.406
SY=-91347.000
-3902.844
S1=-22575.914
S2=-91568.492
176.752
E=13
SX=-5104.269
SY=-129494.438
-11708.750
S1=-4011.727
S2=-130586.977
174.669
E=14
SX=969.672
SY=-108176.375
TAU=-21424.750
S1=5024.582
S2=-112231.281
169.283
E=15
SX=-14954.572
SY=-110883.469
TAU=-
18383.531
S1=-11552.273
S2=-114285.766
169.515
E=16
SX=-19890.141
SY=-86924.312
TAU=-25131.188
S1=-11514.844
S2=-95299.609
161.569
E=17
SX=-22109.688
SY=-87301.625
10225.406
S1=-20543.453
S2=-88867.859
171.292
E=18
SX=-35190.453
SY=-77219.000
-9162.000
S1=-33280.023
S2=-79129.430
168.222
E=19
SX=-9785.850
SY=-171444.172
-20524.969
S1=-7220.594
S2=-174009.422
172.876
E=20
SX=4594.438
SY=-113592.375
TAU=-46145.688
S1=20477.398
S2=-129475.336
161.007
E=21
SX=-25287.307
SY=-118672.312
-30023.750
S1=-16467.512
S2=-127492.109
163.629
E=22
SX=-30634.422
SY=-71127.188
TAU=-44991.469
S1=-1543.715
S2=-100217.891
147.114
E=23
SX=-34259.609
SY=-71743.438
TAU=-14637.906
S1=-29220.699
S2=-76782.344
161.005
E=24
SX=-43958
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 有限元 编程