有限元法上机实习作业.docx

有限元法上机实习作业.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.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.21282.