有限元法上机实习作业.docx
《有限元法上机实习作业.docx》由会员分享，可在线阅读，更多相关《有限元法上机实习作业.docx（26页珍藏版）》请在冰豆网上搜索。
有限元法上机实习作业
有限元法上机实习作业
（一）
——弹性力学平面问题有限元分析
一、设计目的
掌握用有限元方法分析弹性力学平面问题的步骤和技巧，熟练使用编制的程序解决工程实际问题，并分析计算结果的正确性，分析应力分布和位移分布规律。
二、设计任务
某地下岩石洞室断面为直墙拱形，如图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.57480697.6279110.21226
274.88814-6.262677010.4531013.99670
287.39186-9.48812711.255002.50050
2910.39634-13.35866722.325554.64212
3014.00170-18.00330733.610207.21206
312.80170-2.40330745.1517910.29599
324.88814-4.19168757.0017013.99670
337.39186-6.3377376.665802.64770
3410.39634-8.91299771.194104.76190
3514.00170-12.00330781.828057.29893
362.80170-1.20440792.5888010.34337
374.88814-2.09838803.5017013.99670
387.39186-3.1711681.001702.79670
3910.39634-4.4585082.001704.88314
4014.00170-6.0033083.001707.38686
412.80170-.0033084.0017010.39134
424.88814-.0033085.0017013.99670
437.39186-.00330
NODALNUMBERSOFELEMENT
LENODESLENODESLENODES
11264528293389565762
22764629343390566261
32374729303491575863
43874830353492576362
53484931323693585964
64985032373694586463
74595132333795596065
851095233383796596564
967115333343897616267
10712115434393898616766
1178125534353999626368
128131256354039100626867
13891357363741101636469
149141358374241102636968
159101459373842103646570
1610151460384342104647069
1711121661383943105666772
1812171662394443106667271
1912131763394044107676873
2013181764404544108677372
2113141865414247109686974
2214191866414746110687473
2314151967424348111697075
2415201968424847112697574
2516172169434449113717277
2617222170434948114717776
2717182271444550115727378
2818232272445049116727877
2918192373464752117737479
3019242374465251118737978
3119202475474853119747580
3220252476475352120748079
3321222677484954121767782
3422272678485453122768281
3522232779495055123777883
3623282780495554124778382
3723242881515257125787984
3824292882515756126788483
3924252983525358127798085
4025302984525857128798584
4126273185535459129102030
4227323186535958130000
4327283287545560131000
NODALLOADS
NO.VALUENO.VALUE
499.600011911.2000
5919.20001295.6000
6919.20001309.8000
7919.200014019.6000
8915.200015019.6000
9911.200016019.6000
10911.20001709.8000
NODALDISPLACEMENTS
NOUVNOUV
1.0000E+00.9515E-0344.2233E-03.1328E-02
2.0000E+00.9301E-0345.2565E-03.1302E-02
3.0000E+00.8300E-0346.1601E-03.1481E-02
4.0000E+00.7049E-0347.1558E-03.1456E-02
5.0000E+00.6245E-0348.1779E-03.1449E-02
6.4263E-04.9632E-0349.2129E-03.1451E-02
7.4373E-04.9346E-0350.2570E-03.1472E-02
8.4399E-04.8247E-0351.1442E-03.1538E-02
9.4818E-04.6788E-0352.1423E-03.1520E-02
10.2321E-04.4791E-0353.1679E-03.1534E-02
11.7979E-04.9847E-0354.2098E-03.1567E-02
12.8121E-04.9437E-0355.2629E-03.1625E-02
13.8225E-04.8112E-0356.1289E-03.1581E-02
14.7910E-04.6299E-0357.1306E-03.1574E-02
15-.1016E-04.4245E-0358.1620E-03.1611E-02
16.1169E-03.1019E-0259.2140E-03.1674E-02
17.1128E-03.9513E-0360.2785E-03.1769E-02
18.1155E-03.7823E-0361.1156E-03.1612E-02
19.1176E-03.5289E-0362.1201E-03.1621E-02
20.1591E-04.0000E+0063.1596E-03.1683E-02
21.1509E-03.1067E-0264.2253E-03.1780E-02
22.1379E-03.9471E-0365.3033E-03.1916E-02
23.1415E-03.7414E-0366.9918E-04.1642E-02
24.1179E-03.3989E-0367.1001E-03.1658E-02
25.4028E-04.0000E+0068.1283E-03.1716E-02
26.1756E-03.1141E-0269.1879E-03.1808E-02
27.1633E-03.1037E-0270.2689E-03.1941E-02
28.1740E-03.8776E-0371.7504E-04.1676E-02
29.1746E-03.6094E-0372.7401E-04.1695E-02
30.1486E-03.2285E-0373.9082E-04.1757E-02
31.1924E-03.1220E-0274.1430E-03.1842E-02
32.1817E-03.1140E-0275.2417E-03.1975E-02
33.1956E-03.1030E-0276.4076E-04.1707E-02
34.2088E-03.8584E-0377.4001E-04.1726E-02
35.2067E-03.6220E-0378.4341E-04.1801E-02
36.1967E-03.1310E-0279.8234E-04.1890E-02
37.1872E-03.1259E-0280.2301E-03.2028E-02
38.2043E-03.1192E-0281.0000E+00.1726E-02
39.2284E-03.1098E-0282.0000E+00.1743E-02
40.2486E-03.9788E-0383.0000E+00.1829E-02
41.1721E-03.1418E-0284.0000E+00.1985E-02
42.1708E-03.1386E-0285.2502E-03.2135E-02
43.1921E-03.1358E-02
RFACTIVEFORCESOFNODALCONSTRAINT
NO.VALUENO.VALUE
40-41.941739-4.545105
50-5.066434161-8.261791
1-12.50068163-11.48822
3-19.26372165-12.56956
5-18.17349167-20.91741
7-17.21299
STRESSOFELEMENT
LEXXYYXYS1S2CETA
17.11953.0113.62607.2128
升级会员 