有限元法上机实习作业.docx
《有限元法上机实习作业.docx》由会员分享,可在线阅读,更多相关《有限元法上机实习作业.docx(23页珍藏版)》请在冰豆网上搜索。
有限元法上机实习作业
有限元法上机实习作业
(一)
——弹性力学平面问题有限元分析
一、设计目的
掌握用有限元方法分析弹性力学平面问题的步骤和技巧,熟练使用编制的程序解决工程实际问题,并分析计算结果的正确性,分析应力分布和位移分布规律。
二、设计任务
某地下岩石洞室断面为直墙拱形,如图SII.1所示。
岩石的弹性参数E=1.0105MPa,μ=0.25,在远处作用有水平均布荷载qx和竖向均布荷载qy,不计岩石重力。
取5rX5(r+h)计算范围,按平面应变问题取二分之一洞室计算岩石的位移、应变和应力,并对计算结果进行必要的讨论。
三、设计条件
直墙拱形断面洞室。
洞室断面如图SII.1所示。
图SII.1直墙拱形断面洞室示意图
序号
学号
r(m)
h(m)
qx(MPa)
qy(MPa)
15
086907
2.8
4.8
3.2
5.6
四、设计结果
1、单元划分
1)节点坐标
.00170,-4.80330
.00174,-8.38033
.00178,-12.67277
.00184,-17.82369
.00190,-24.00480
.74530,-4.80330
1.25879,-8.38033
1.87497,-12.67277
2.61439,-17.82369
3.50170,-24.00480
1.40170,-4.80330
2.44494,-8.38033
3.69682,-12.67277
5.19909,-17.82369
7.00180,-24.00480
2.09110,-4.80330
3.64885,-8.38033
5.51815,-12.67277
7.76131,-17.82369
10.45310,-24.00480
2.80170,-4.80330
4.88814,-8.38033
7.39186,-12.67277
10.39634,-17.82369
14.00170,-24.00480
2.80170,-3.57480
4.88814,-6.26267
7.39186,-9.48812
10.39634,-13.35866
14.00170,-18.00330
2.80170,-2.40330
4.88814,-4.19168
7.39186,-6.33773
10.39634,-8.91299
14.00170,-12.00330
2.80170,-1.20440
4.88814,-2.09838
7.39186,-3.17116
10.39634,-4.45850
14.00170,-6.00330
2.80170,-.00330
4.88814,-.00330
7.39186,-.00330
10.39634,-.00330
14.00170,-.00330
2.71740,.67870
4.81954,1.20366
7.34211,1.83362
10.36920,2.58956
14.00170,3.49670
2.50470,1.25170
4.64647,2.32193
7.21659,3.60621
10.30073,5.14734
14.00170,6.99670
2.24420,1.67330
4.43449,3.30793
7.06285,5.26949
10.21687,7.62336
14.00170,10.44800
1.98230,1.97590
4.22138,4.21524
6.90828,6.90246
10.13256,10.12711
14.00170,13.99670
1.67670,2.24040
3.31165,4.43047
5.27359,7.05856
7.62791,10.21226
10.45310,13.99670
1.25500,2.50050
2.32555,4.64212
3.61020,7.21206
5.15179,10.29599
7.00170,13.99670
.66580,2.64770
1.19410,4.76190
1.82805,7.29893
2.58880,10.34337
3.50170,13.99670
.00170,2.79670
.00170,4.88314
.00170,7.38686
.00170,10.39134
.00170,13.99670
2)单元定义
1,2,6
2,7,6
2,3,7
3,8,7
3,4,8
4,9,8
4,5,9
5,10,9
6,7,11
7,12,11
7,8,12
8,13,12
8,9,13
9,14,13
9,10,14
10,15,14
11,12,16
12,17,16
12,13,17
13,18,17
13,14,18
14,19,18
14,15,19
15,20,19
16,17,21
17,22,21
17,18,22
18,23,22
18,19,23
19,24,23
19,20,24
20,25,24
21,22,26
22,27,26
22,23,27
23,28,27
23,24,28
24,29,28
24,25,29
25,30,29
26,27,31
27,32,31
27,28,32
28,33,32
28,29,33
29,34,33
29,30,34
30,35,34
31,32,36
32,37,36
32,33,37
33,38,37
33,34,38
34,39,38
34,35,39
35,40,39
36,37,41
37,42,41
37,38,42
38,43,42
38,39,43
39,44,43
39,40,44
40,45,44
41,42,47
41,47,46
42,43,48
42,48,47
43,44,49
43,49,48
44,45,50
44,50,49
46,47,52
46,52,51
47,48,53
47,53,52
48,49,54
48,54,53
49,50,55
49,55,54
51,52,57
51,57,56
52,53,58
52,58,57
53,54,59
53,59,58
54,55,60
54,60,59
56,57,62
56,62,61
57,58,63
57,63,62
58,59,64
58,64,63
59,60,65
59,65,64
61,62,67
61,67,66
62,63,68
62,68,67
63,64,69
63,69,68
64,65,70
64,70,69
66,67,72
66,72,71
67,68,73
67,73,72
68,69,74
68,74,73
69,70,75
69,75,74
71,72,77
71,77,76
72,73,78
72,78,77
73,74,79
73,79,78
74,75,80
74,80,79
76,77,82
76,82,81
77,78,83
77,83,82
78,79,84
78,84,83
79,80,85
79,85,84
2、应力分布和位移分布规律
PLANEPROBLEM
LXM=1NJ=85NE=129NZ=11NPJ=14
LDD=14E0=.10E+06EMU=.25OU=.00TE=1.00
NODALCO-ORDINATES
NOXYNOXY
1.00170-4.803304410.39634-.00330
2.00174-8.380334514.00170-.00330
3.00178-12.67277462.71740.67870
4.00184-17.82369474.819541.20366
5.00190-24.00480487.342111.83362
6.74530-4.803304910.369202.58956
71.25879-8.380335014.001703.49670
81.87497-12.67277512.504701.25170
92.61439-17.82369524.646472.32193
103.50170-24.00480537.216593.60621
111.40170-4.803305410.300735.14734
122.44494-8.380335514.001706.99670
133.69682-12.67277562.244201.67330
145.19909-17.82369574.434493.30793
157.00180-24.00480587.062855.26949
162.09110-4.803305910.216877.62336
173.64885-8.380336014.0017010.44800
185.51815-12.67277611.982301.97590
197.76131-17.82369624.221384.21524
2010.45310-24.00480636.908286.90246
212.80170-4.803306410.1325610.12711
224.88814-8.380336514.0017013.99670
237.39186-12.67277661.676702.24040
2410.39634-17.82369673.311654.43047
2514.00170-24.00480685.273597.05856
262.80170-3.57