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

1、有限元法上机实习作业有限元法上机实习作业（一）弹性力学平面问题有限元分析一、设计目的掌握用有限元方法分析弹性力学平面问题的步骤和技巧，熟练使用编制的程序解决工程实际问题，并分析计算结果的正确性，分析应力分布和位移分布规律。二、设计任务某地下岩石洞室断面为直墙拱形，如图SII.1所示。岩石的弹性参数E 1.0105MPa，= 0.25，在远处作用有水平均布荷载qx和竖向均布荷载qy，不计岩石重力。取5rX5（r+h）计算范围，按平面应变问题取二分之一洞室计算岩石的位移、应变和应力，并对计算结果进行必要的讨论。三、设计条件直墙拱形断面洞室。洞室断面如图SII.1所示。图SII.1 直墙拱形断面洞室

2、示意图序号学号r (m)h (m)q x (MPa)qy (MPa)150869072.84.83.25.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, -1

3、2.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, -

4、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

5、 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.90

6、828, 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.5017

7、0, 13.99670 .00170, 2.79670 .00170, 4.88314 .00170, 7.38686 .00170, 10.39134 .00170, 13.996702）单元定义 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

8、, 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, 3

9、2, 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,

10、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

11、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, 842、应力分布和位移分布规律 PLANE PROBLE

12、M LXM= 1 NJ= 85 NE=129 NZ= 11 NPJ=14 LDD= 14 E0= .10E+06 EMU= .25 OU= .00 TE= 1.00 NODAL CO-ORDINATES NO X Y NO X Y 1 .00170 -4.80330 44 10.39634 -.00330 2 .00174 -8.38033 45 14.00170 -.00330 3 .00178 -12.67277 46 2.71740 .67870 4 .00184 -17.82369 47 4.81954 1.20366 5 .00190 -24.00480 48 7.34211 1.8

13、3362 6 .74530 -4.80330 49 10.36920 2.58956 7 1.25879 -8.38033 50 14.00170 3.49670 8 1.87497 -12.67277 51 2.50470 1.25170 9 2.61439 -17.82369 52 4.64647 2.32193 10 3.50170 -24.00480 53 7.21659 3.60621 11 1.40170 -4.80330 54 10.30073 5.14734 12 2.44494 -8.38033 55 14.00170 6.99670 13 3.69682 -12.67277

14、 56 2.24420 1.67330 14 5.19909 -17.82369 57 4.43449 3.30793 15 7.00180 -24.00480 58 7.06285 5.26949 16 2.09110 -4.80330 59 10.21687 7.62336 17 3.64885 -8.38033 60 14.00170 10.44800 18 5.51815 -12.67277 61 1.98230 1.97590 19 7.76131 -17.82369 62 4.22138 4.21524 20 10.45310 -24.00480 63 6.90828 6.9024

15、6 21 2.80170 -4.80330 64 10.13256 10.12711 22 4.88814 -8.38033 65 14.00170 13.99670 23 7.39186 -12.67277 66 1.67670 2.24040 24 10.39634 -17.82369 67 3.31165 4.43047 25 14.00170 -24.00480 68 5.27359 7.05856 26 2.80170 -3.57480 69 7.62791 10.21226 27 4.88814 -6.26267 70 10.45310 13.99670 28 7.39186 -9

16、.48812 71 1.25500 2.50050 29 10.39634 -13.35866 72 2.32555 4.64212 30 14.00170 -18.00330 73 3.61020 7.21206 31 2.80170 -2.40330 74 5.15179 10.29599 32 4.88814 -4.19168 75 7.00170 13.99670 33 7.39186 -6.33773 76 .66580 2.64770 34 10.39634 -8.91299 77 1.19410 4.76190 35 14.00170 -12.00330 78 1.82805 7

17、.29893 36 2.80170 -1.20440 79 2.58880 10.34337 37 4.88814 -2.09838 80 3.50170 13.99670 38 7.39186 -3.17116 81 .00170 2.79670 39 10.39634 -4.45850 82 .00170 4.88314 40 14.00170 -6.00330 83 .00170 7.38686 41 2.80170 -.00330 84 .00170 10.39134 42 4.88814 -.00330 85 .00170 13.99670 43 7.39186 -.00330 NO

18、DAL NUMBERS OF ELEMENT LE NODES LE NODES LE NODES 1 1 2 6 45 28 29 33 89 56 57 62 2 2 7 6 46 29 34 33 90 56 62 61 3 2 3 7 47 29 30 34 91 57 58 63 4 3 8 7 48 30 35 34 92 57 63 62 5 3 4 8 49 31 32 36 93 58 59 64 6 4 9 8 50 32 37 36 94 58 64 63 7 4 5 9 51 32 33 37 95 59 60 65 8 5 10 9 52 33 38 37 96 59

19、 65 64 9 6 7 11 53 33 34 38 97 61 62 67 10 7 12 11 54 34 39 38 98 61 67 66 11 7 8 12 55 34 35 39 99 62 63 68 12 8 13 12 56 35 40 39 100 62 68 67 13 8 9 13 57 36 37 41 101 63 64 69 14 9 14 13 58 37 42 41 102 63 69 68 15 9 10 14 59 37 38 42 103 64 65 70 16 10 15 14 60 38 43 42 104 64 70 69 17 11 12 16

20、 61 38 39 43 105 66 67 72 18 12 17 16 62 39 44 43 106 66 72 71 19 12 13 17 63 39 40 44 107 67 68 73 20 13 18 17 64 40 45 44 108 67 73 72 21 13 14 18 65 41 42 47 109 68 69 74 22 14 19 18 66 41 47 46 110 68 74 73 23 14 15 19 67 42 43 48 111 69 70 75 24 15 20 19 68 42 48 47 112 69 75 74 25 16 17 21 69

21、43 44 49 113 71 72 77 26 17 22 21 70 43 49 48 114 71 77 76 27 17 18 22 71 44 45 50 115 72 73 78 28 18 23 22 72 44 50 49 116 72 78 77 29 18 19 23 73 46 47 52 117 73 74 79 30 19 24 23 74 46 52 51 118 73 79 78 31 19 20 24 75 47 48 53 119 74 75 80 32 20 25 24 76 47 53 52 120 74 80 79 33 21 22 26 77 48 4

22、9 54 121 76 77 82 34 22 27 26 78 48 54 53 122 76 82 81 35 22 23 27 79 49 50 55 123 77 78 83 36 23 28 27 80 49 55 54 124 77 83 82 37 23 24 28 81 51 52 57 125 78 79 84 38 24 29 28 82 51 57 56 126 78 84 83 39 24 25 29 83 52 53 58 127 79 80 85 40 25 30 29 84 52 58 57 128 79 85 84 41 26 27 31 85 53 54 59

23、 129 10 20 30 42 27 32 31 86 53 59 58 130 0 0 0 43 27 28 32 87 54 55 60 131 0 0 0 NODAL LOADS NO. VALUE NO. VALUE 49 9.6000 119 11.2000 59 19.2000 129 5.6000 69 19.2000 130 9.8000 79 19.2000 140 19.6000 89 15.2000 150 19.6000 99 11.2000 160 19.6000 109 11.2000 170 9.8000 NODAL DISPLACEMENTS NO U V N

24、O U V 1 .0000E+00 .9515E-03 44 .2233E-03 .1328E-02 2 .0000E+00 .9301E-03 45 .2565E-03 .1302E-02 3 .0000E+00 .8300E-03 46 .1601E-03 .1481E-02 4 .0000E+00 .7049E-03 47 .1558E-03 .1456E-02 5 .0000E+00 .6245E-03 48 .1779E-03 .1449E-02 6 .4263E-04 .9632E-03 49 .2129E-03 .1451E-02 7 .4373E-04 .9346E-03 50

25、 .2570E-03 .1472E-02 8 .4399E-04 .8247E-03 51 .1442E-03 .1538E-02 9 .4818E-04 .6788E-03 52 .1423E-03 .1520E-02 10 .2321E-04 .4791E-03 53 .1679E-03 .1534E-02 11 .7979E-04 .9847E-03 54 .2098E-03 .1567E-02 12 .8121E-04 .9437E-03 55 .2629E-03 .1625E-02 13 .8225E-04 .8112E-03 56 .1289E-03 .1581E-02 14 .7

26、910E-04 .6299E-03 57 .1306E-03 .1574E-02 15 -.1016E-04 .4245E-03 58 .1620E-03 .1611E-02 16 .1169E-03 .1019E-02 59 .2140E-03 .1674E-02 17 .1128E-03 .9513E-03 60 .2785E-03 .1769E-02 18 .1155E-03 .7823E-03 61 .1156E-03 .1612E-02 19 .1176E-03 .5289E-03 62 .1201E-03 .1621E-02 20 .1591E-04 .0000E+00 63 .1

27、596E-03 .1683E-02 21 .1509E-03 .1067E-02 64 .2253E-03 .1780E-02 22 .1379E-03 .9471E-03 65 .3033E-03 .1916E-02 23 .1415E-03 .7414E-03 66 .9918E-04 .1642E-02 24 .1179E-03 .3989E-03 67 .1001E-03 .1658E-02 25 .4028E-04 .0000E+00 68 .1283E-03 .1716E-02 26 .1756E-03 .1141E-02 69 .1879E-03 .1808E-02 27 .16

28、33E-03 .1037E-02 70 .2689E-03 .1941E-02 28 .1740E-03 .8776E-03 71 .7504E-04 .1676E-02 29 .1746E-03 .6094E-03 72 .7401E-04 .1695E-02 30 .1486E-03 .2285E-03 73 .9082E-04 .1757E-02 31 .1924E-03 .1220E-02 74 .1430E-03 .1842E-02 32 .1817E-03 .1140E-02 75 .2417E-03 .1975E-02 33 .1956E-03 .1030E-02 76 .407

29、6E-04 .1707E-02 34 .2088E-03 .8584E-03 77 .4001E-04 .1726E-02 35 .2067E-03 .6220E-03 78 .4341E-04 .1801E-02 36 .1967E-03 .1310E-02 79 .8234E-04 .1890E-02 37 .1872E-03 .1259E-02 80 .2301E-03 .2028E-02 38 .2043E-03 .1192E-02 81 .0000E+00 .1726E-02 39 .2284E-03 .1098E-02 82 .0000E+00 .1743E-02 40 .2486

30、E-03 .9788E-03 83 .0000E+00 .1829E-02 41 .1721E-03 .1418E-02 84 .0000E+00 .1985E-02 42 .1708E-03 .1386E-02 85 .2502E-03 .2135E-02 43 .1921E-03 .1358E-02 RFACTIVE FORCES OF NODAL CONSTRAINT NO. VALUE NO. VALUE 40 -41.94173 9 -4.545105 50 -5.066434 161 -8.261791 1 -12.50068 163 -11.48822 3 -19.26372 165 -12.56956 5 -18.17349 167 -20.91741 7 -17.21299 STRESS OF ELEMENT LE XX YY XY S1 S2 CETA 1 7.1195 3.0113 .6260 7.2128 2.