题目：直墙拱形断面洞室。
已知解析：取直墙拱形断面洞室的右半部分用有限元法。用TGEN法划分网格。
模型尺寸：总宽27.2m，总高35.2m，模型原点坐标设在半圆拱圆心，圆心以下22m，以上13.2m。设置21条生成线，每条线的单元生成个数（间隔数）为5，加权（比例）系数仍为1.2。考虑设置4个分区，即5条分区界限生成线，各分区角点的节点号为1-6-36-31、31-36-66-61、61-66-96-91及91-96-126-121。
一．划分网格
用EXCEL表格求21条生成线的起点终点坐标。得：生成线 x1 y1 x2 1 0.00 -4.40 0.00 2 0.68 -4.40 5.44 3 1.36 -4.40 10.88 4 2.04 -4.40 16.32 5 2.72 -4.40 21.76 6 3.40 -4.40 27.20 7 3.40 -3.52 27.20 8 3.40 -2.64 27.20 9 3.40 -1.76 27.20 10 3.40 -0.88 27.20 11 3.40 0.00 27.20 12 9.00 3.36 0.53 27.20 13 18.00 3.23 1.05 27.20 14 27.00 3.03 1.54 27.20 15 36.00 2.75 2.00 27.20 16 45.00 2.40 2.40 27.20 y2 -22.00 -22.00 -22.00 -22.00 -22.00 -22.00 -17.60 -13.20 -8.80 -4.40 0.00 2.64 5.28 7.92 10.56 13.20
17 18 19 20 21
54.00 63.00 72.00 81.00
2.00 1.54 1.05 0.53 3.40
2.75 3.03 3.23 3.36 0.00
21.76 16.32 10.88 5.44 0.00
13.20 13.30 13.20 13.20 13.20
二．输入数据文件TGEN1.DAT为：
21,1.2
5,0.00,-4.40,0.00,-22.00 5,0.68,-4.40,5.44,-22.00 5,1.36,-4.40,10.88,-22.00 5,2.04,-4.40,16.32,-22.00 5,2.72,-4.40,21.76,-22.00 5,3.40,-4.40,27.20,-22.00 5,3.40,-3.52,27.20,-17.60 5,3.40,-2.64,27.20,-13.20 5,3.40,-1.76,27.20,-8.80 5,3.40,-0.88,27.20,-4.40 5,3.40,0.00,27.20,0.00
5,3.36,0.53,27.20,2.64 5,3.23,1.05,27.20,5.28 5,3.03,1.54,27.20,7.92 5,2.75,2.00,27.20,10.56 5,2.40,2.40,27.20,13.20 5,2.00,2.75,21.76,13.20 5,1.54,3.03,16.32,13.30 5,1.05,3.23,10.88,13.20 5,0.53,3.36,5.44,13.20 5,0.00,3.40,0.00,13.20 5,1,6,11,16,21
三．输入文件数据
(1) 控制数据：126,200,16,17 （2）节点坐标
NODAL CO-ORDINATES 1 0.0 -4.4 2 0.0 -6.76508 3 0.0 -9.60318 4 0.0 -13.0089 5 0.0 -17.09576 6 0.0 -22.0 7 0.68 -4.4 8 1.31965 -6.76508 9 2.08722 -9.60318 10 3.00832 -13.0089 11 4.11363 -17.09576 12 5.44 -22.0 13 1.36 -4.4 14 2.63929 -6.76508 15 4.17445 -9.60318 16 6.01663 -13.0089
NODAL 64 15.04158 65 20.56814 66 27.2 67 3.36 68 6.56361 69 10.40795 70 15.02115 71 20.55699 72 27.2 73 3.23 74 6.45108 75 10.31638 76 14.95474 77 20.52077 78 27.2 79 3.03 CO-ORDINATES 0.0 0.0 0.0 0.53 0.81354 1.15379 1.56209 2.05205 2.64 1.05 1.61843 2.30054 3.11907 4.10131 5.28 1.54
17 8.22725 -17.09576 80 6.27796 2.39734 18 10.88 -22.0 81 10.17551 3.42615 19 2.04 -4.4 82 14.85257 4.66073 20 3.95894 -6.76508 83 20.46504 6.14221 21 6.26167 -9.60318 84 27.2 7.92 22 9.02495 -13.0089 85 2.75 2.0
23 12.34088 -17.09576 86 6.03558 3.15029 24 16.32 -22.0 87 9.97828 4.53064 25 2.72 -4.4 88 14.70952 6.18706 26 5.27859 -6.76508 89 20.38701 8.17476 27 8.3489 -9.60318 90 27.2 10.56 28 12.03327 -13.0089 91 2.4 29 16.45451 -17.09576 92 5.73262 30 21.76 -22.0 93 9.73176 31 3.4 -4.4 94 14.53073 32 6.59824 -6.76508 95 20.28949 33 10.43612 -9.60318 96 27.2 34 15.04158 -13.0089 97 2.0 35 20.56814 -17.09576 98 4.65534 36 27.2 -22.0 99 7.84175 37 3.4 -3.52 100 11.66545 38 6.59824 -5.41207 101 16.25388 39 10.43612 -7.68255 102 21.76 40 15.04158 -10.40712 103 1.54 41 20.56814 -13.67661 104 3.52613 42 27.2 -17.6 105 5.90949 43 3.4 -2.64 106 8.76952 44 6.59824 -4.05905 107 12.20156 45 10.43612 -5.76191 108 16.32 46 15.04158 -7.80534 109 1.05 47 20.56814 -10.25746 110 2.37095 48 27.2 -13.2 111 3.9561 49 3.4 -1.76 112 5.85827 50 6.59824 -2.70603 113 8.14087 51 10.43612 -3.84127 114 10.88 52 15.04158 -5.20356 115 0.53 53 20.56814 -6.83831 116 1.1898 54 27.2 -8.8 117 1.98157 55 3.4 -0.88 118 2.93169 56 6.59824 -1.35302 119 4.07183 57 10.43612 -1.92064 120 5.44 58 15.04158 -2.60178 121 0.0 59 20.56814 -3.41915 122 0.0 60 27.2 -4.4 123 0.0 61 3.4 0.0 124 0.0 62 6.59824 0.0 125 0.0 63 10.43612 0.0 126 0.0
（3）单元定义
ELEMENT DEFINITIONS
ELEMENT NODE1 NODE2 NODE3 ELEMENT 1 1 2 7 101 61 2 2 8 7 102 62 2.4 3.8513 5.59286 7.68273 10.19058 13.2 2.75 4.15427 5.83939 7.86154 10.28811 13.2 3.03 4.41008 6.06617 8.05349 10.43827 13.3 3.23 4.56977 6.17748 8.10675 10.42186 13.2 3.36 4.6823 6.26905 8.17316 10.45809 13.2 3.4 4.71692 6.29723 8.19359 10.46923 13.2 NODE1 NODE2 62 67 68 67 NODE3
3 4 5 6 7 8 9 10 11 12 13 2 3 8 103 62 63 68 3 9 8 104 63 69 68 3 4 9 105 63 64 69 4 10 9 106 64 70 69 4 5 10 107 64 65 70 5 11 10 108 65 71 70 5 6 11 109 65 66 71 6 12 11 110 66 72 71 7 8 13 111 67 68 73 8 14 13 112 68 74 73 8 9 14 113 68 69 74 14 9 15 9 16 10 17 10 18 11 19 11 20 12 21 13 22 14 23 14 24 15 25 15 26 16 27 16 28 17 29 17 30 18 31 19 32 20 33 20 34 21 35 21 36 22 37 22 38 23 39 23 40 24 41 25 42 26 43 26 44 27 45 27 46 28 47 28 48 29 49 29 50 30 51 31 52 31 53 32 54 32 55 33 15 14 10 15 16 15 11 16 17 16 12 17 18 17 14 19 20 19 15 20 21 20 16 21 22 21 17 22 23 22 18 23 24 23 20 25 26 25 21 26 27 26 22 27 28 27 23 28 29 28 24 29 30 29 26 31 32 31 27 32 33 32 28 33 34 33 29 34 35 34 30 35 36 35 38 37 32 38 39 38 33 39 40 39 114 69 115 69 116 70 117 70 118 71 119 71 120 72 121 73 122 74 123 74 124 75 125 75 126 76 127 76 128 77 129 77 130 78 131 79 132 80 133 80 134 81 135 81 136 82 137 82 138 83 139 83 140 84 141 85 142 86 143 86 144 87 145 87 146 88 147 88 148 89 149 89 150 90 151 91 152 91 153 92 154 92 155 93 75 74 70 75 76 75 71 76 77 76 72 77 78 77 74 79 80 79 75 80 81 80 76 81 82 81 77 82 83 82 78 83 84 83 80 85 86 85 81 86 87 86 82 87 88 87 83 88 89 88 84 89 90 89 86 91 92 91 87 92 93 92 88 93 94 93 89 94 95 94 90 95 96 95 98 97 92 98 99 98 93 99 100 99
56 33 34 40 156 93 94 100 57 34 41 40 157 94 101 100 58 34 35 41 158 94 95 101 59 35 42 41 159 95 102 101 60 35 36 42 160 95 96 102 61 37 44 43 161 97 104 103 62 37 38 44 162 97 98 104 63 38 45 44 163 98 105 104 64 38 39 45 164 98 99 105 65 39 46 45 165 99 106 105 66 39 40 46 166 99 100 106 67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
(四)节点荷载
NODE
36
42
48
54 40 47 46 167 40 41 47 168 41 48 47 169 41 42 48 170 43 50 49 171 43 44 50 172 44 51 50 173 44 45 51 174 45 52 51 175 45 46 52 176 46 53 52 177 46 47 53 178 47 54 53 179 47 48 54 180 49 56 55 181 49 50 56 182 50 57 56 183 50 51 57 184 51 58 57 185 51 52 58 186 52 59 58 187 52 53 59 188 53 60 59 189 53 54 60 190 55 62 61 191 55 56 62 192 56 63 62 193 56 57 63 194 57 64 63 195 57 58 64 196 58 65 64 197 58 59 65 198 59 66 65 199 59 60 66 200 APPLIED LOADS APPLIED LOADS NODE 7.92 0.0 84 9.504 15.84 0.0 90 9.504 15.84 0.0 96 4.752 15.84 0.0 102 0.0 100 107 106 100 101 107 101 108 107 101 102 108 103 110 109 103 104 110 104 111 110 104 105 111 105 112 111 105 106 112 106 113 112 106 107 113 107 114 113 107 108 114 109 116 115 109 110 116 110 117 116 110 111 117 111 118 117 111 112 118 112 119 118 112 113 119 113 120 119 113 114 120 115 122 121 115 116 122 116 123 122 116 117 123 117 124 123 117 118 124 118 125 124 118 119 125 119 126 125 119 120 126 APPLIED LOADS 0.0 0.0 14.144 28.288
60 15.84 0.0 108 0.0 28.288 66 12.672 0.0 114 0.0 28.288 72 9.504 0.0 120 0.0 28.288 78 9.504 0.0 126 0.0 14.144
（五）节点约束条件
FIXITIES
1 1 0 30 0 1
2 1 0 36 0 1
3 1 0 121 1 0
4 1 0 122 1 0
5 1 0 123 1 0
6 1 1 124 1 0
12 0 1 125 1 0
18 0 1 126 1 0
24 0 1
（六）材料弹性性质参数
100000.0 0.25 1.0
四．有限元计算结果
********************** PLAIN PROBLEM ********************** NODAL DISPLACEMENTS
NODE X-COMP Y-COMP NODE X-COMP 1 6.02955E-28 5.05912E-04 64 3.52661E-04 5.05912E-04 2 1.24586E-27 5.01285E-04 65 4.55571E-04 5.01285E-04 3 1.37341E-27 4.43734E-04 66 5.84015E-04 4.43734E-04 4 1.53831E-27 3.38286E-04 67 2.59432E-04 3.38286E-04 5 1.76901E-27 1.89818E-04 68 2.39239E-04 1.89818E-04 6 9.55371E-28 1.36420E-27 69 2.79573E-04 1.36420E-27 7 3.30830E-05 5.16572E-04 70 3.56334E-04 5.16572E-04 8 5.33626E-05 5.16240E-04 71 4.61602E-04 5.16240E-04 9 6.60307E-05 4.56459E-04 72 5.93462E-04 4.56459E-04 10 7.96198E-05 3.46857E-04 73 2.43422E-04 3.46857E-04 11 9.81679E-05 1.94038E-04 74 2.33451E-04 1.94038E-04 12 1.25034E-04 2.76723E-27 75 2.78571E-04 2.76723E-27 13 6.60023E-05 5.44098E-04 76 3.59052E-04 5.44098E-04 14 9.97756E-05 5.50508E-04 77 4.68030E-04 5.50508E-04 15 1.22488E-04 4.83616E-04 78 6.03864E-04 4.83616E-04 16 1.50180E-04 3.64431E-04 79 2.24518E-04 3.64431E-04 17 1.88707E-04 2.03432E-04 80 2.25365E-04 2.03432E-04 18 2.40359E-04 2.84343E-27 81 2.75878E-04 2.84343E-27 19 9.71639E-05 5.85358E-04 82 3.61176E-04 5.85358E-04 20 1.35805E-04 5.90742E-04 83 4.74918E-04 5.90742E-04 21 1.68233E-04 5.09417E-04 84 6.15368E-04 5.09417E-04 22 2.12403E-04 3.78776E-04 85 2.02200E-04 3.78776E-04 23 2.72106E-04 2.09781E-04 86 2.15103E-04 2.09781E-04 24 3.48558E-04 2.87783E-27 87 2.71689E-04 2.87783E-27 25 1.25351E-04 6.33848E-04 88 3.62854E-04 6.33848E-04 26 1.63401E-04 6.23855E-04 89 4.82260E-04 6.23855E-04 27 2.08513E-04 5.25480E-04 90 6.27788E-04 5.25480E-04 28 2.71011E-04 3.86629E-04 91 1.77505E-04 3.86629E-04 29 3.52373E-04 2.11868E-04 92 2.02937E-04 2.11868E-04 Y-COMP
30 4.54066E-04 2.87433E-27 93 2.66370E-04 2.87433E-27
31 1.50459E-04 6.82968E-04 94 3.64352E-04 6.82968E-04
32 1.86769E-04 6.42982E-04 95 4.90081E-04 6.42982E-04
33 2.47370E-04 5.34424E-04 96 6.40591E-04 5.34424E-04
34 3.28567E-04 3.89515E-04 97 1.51477E-04 3.89515E-04
35 4.31929E-04 2.11424E-04 98 1.79294E-04 2.11424E-04
36 5.59163E-04 1.41704E-27 99 2.29854E-04 1.41704E-27
37 1.97859E-04 7.58630E-04 100 3.07444E-04 7.58630E-04
38 1.90506E-04 7.07456E-04 101 4.11265E-04 7.07456E-04
39 2.51870E-04 6.21941E-04 102 5.33570E-04 6.21941E-04
40 3.31785E-04 5.04287E-04 103 1.19731E-04 5.04287E-04
41 4.33947E-04 3.59185E-04 104 1.48440E-04 3.59185E-04
42 5.60721E-04 1.87143E-04 105 1.88567E-04 1.87143E-04
43 2.34185E-04 8.36693E-04 106 2.46813E-04 8.36693E-04
44 2.07410E-04 7.78965E-04 107 3.29486E-04 7.78965E-04
45 2.56285E-04 7.11317E-04 108 4.31816E-04 7.11317E-04
46 3.35602E-04 6.20372E-04 109 8.36149E-05 6.20372E-04
47 4.37053E-04 5.07480E-04 110 1.09359E-04 5.07480E-04
48 5.63297E-04 3.74327E-04 111 1.39087E-04 3.74327E-04
49 2.58278E-04 9.15523E-04 112 1.78938E-04 9.15523E-04
50 2.26192E-04 8.57987E-04 113 2.37326E-04 8.57987E-04
51 2.64380E-04 8.04200E-04 114 3.23687E-04 8.04200E-04
52 3.40191E-04 7.37672E-04 115 4.32614E-05 7.37672E-04
53 4.41438E-04 6.56271E-04 116 5.93841E-05 6.56271E-04
54 5.67704E-04 5.61307E-04 117 7.62164E-05 5.61307E-04
55 2.70403E-04 9.94582E-04 118 9.76415E-05 9.94582E-04
56 2.39106E-04 9.42038E-04 119 1.29056E-04 9.42038E-04
57 2.73003E-04 9.01000E-04 120 1.83448E-04 9.01000E-04
58 3.46245E-04 8.56571E-04 121 5.55903E-28 8.56571E-04
59 4.47455E-04 8.05374E-04 122 8.51014E-28 8.05374E-04
60 5.74561E-04 7.47870E-04 123 8.87034E-28 7.47870E-04
61 2.72356E-04 1.07436E-03 124 9.96741E-28 1.07436E-03
62 2.42090E-04 1.02730E-03 125 1.17821E-27 1.02730E-03
63 2.78384E-04 9.99660E-04 126 7.18203E-28 1.69179E-03
NODAL REACTIONS
NODE X-COMP Y-COMP NODE X-COMP Y-COMP 1 -6.02954 0.0 30 0.0 -28.74329
2 -12.45863 0.0 36 0.0 -14.17039
3 -13.7341 0.0 121 -5.55903 0.0
4 -15.38307 0.0 122 -8.51014 0.0
5 -17.69007 0.0 123 -8.87034 0.0
6 -9.55371 -13.64197 124 -9.96741 0.0
12 0.0 -27.67233 125 -11.78205 0.0
18 0.0 -28.43428 126 -7.18203 0.0
24 0.0 -28.7783
ELEMENT STRESSES
IEL SIG XX SIG YY TAU XY IEL SIG XX SIG YY TAU XY 1 5.25906 1.93852 0.55736 101 2.29277 9.57014 -1.41552 2 4.42724 1.77965 0.48691 102 1.53587 6.30595 -0.84162 3 5.03427 3.60074 0.40294 103 3.20071 6.9122 -0.36639
4 5 6 7 8 9 10 11 12 13 14 4.18203 3.5475 4.47536 4.42746 3.99474 4.45595 4.11476 4.81602 3.88489 4.86666 3.92167 4.97702 3.85841 5.03752 5.55804 2.90399 4.15461 2.45965 4.7502 4.24644 3.97348 4.22649 0.36227 104 2.91299 5.78695 -0.29153 0.21677 105 3.53017 6.00298 -0.1284 0.21395 106 3.38323 5.48128 -0.13195 0.10131 107 3.6199 5.56307 -0.04804 0.10459 108 3.54747 5.31004 -0.05635 0.03648 109 3.614 5.33298 -0.00357 0.02624 110 3.57996 5.18388 -0.00052 1.5999 111 2.82067 9.45787 -1.9242 1.09197 112 1.8682 6.11165 -1.14345 1.10279 113 3.30258 6.78105 -0.56688 0.69829 114 3.01507 5.72775 -0.40598 15 4.15458 16 3.85852 17 3.88758 18 3.7923 19 3.77642 20 3.73616 21 5.95869 22 3.71081 23 4.33659 24 3.71324 25 3.81982 26 3.70771 27 3.69865 28 3.67732 29 3.6671 30 3.64246 31 6.39765 32 3.34596 33 3.97328 34 3.58693 35 3.64437 36 3.62809 37 3.62363 38 3.62124 39 3.61686 40 3.60482 41 6.86732 42 3.18688 43 3.67882 44 3.5518 45 3.60018 46 3.59972 47 3.59807 48 3.6071 49 3.59478 50 3.59354 51 6.21082 52 3.12316 53 3.5741 54 3.48928 55 3.61893 56 3.54663 4.76981 4.90387 4.99105 5.11638 5.06875 5.17838 4.84112 3.36654 5.24389 4.9056 5.22534 5.23977 5.21261 5.26518 5.23453 5.26991 7.40243 4.08969 5.97164 5.27152 5.44386 5.34271 5.32935 5.31354 5.30039 5.29769 10.09903 4.52367 5.99947 5.36005 5.50519 5.35716 5.35221 5.32023 5.28329 5.28536 10.22242 5.55921 5.70952 5.48352 5.52674 5.36477 0.58085 0.3697 0.27188 0.16931 0.09806 0.03324 2.53101 1.29075 1.32458 0.66552 0.57192 0.30511 0.23216 0.12034 0.07032 0.01951 3.15032 1.12449 1.08916 0.45076 0.36922 0.1897 0.13543 0.06073 0.02666 0.00874 3.41661 0.8209 0.63131 0.28825 0.21596 0.10631 0.06218 0.02679 0.00278 0.00644 3.1548 0.9066 0.30831 0.27569 0.13347 0.08511 115 3.55532 116 3.42152 117 3.62539 3.56146 119 3.62375 120 3.58958 3.33835 122 2.23112 123 3.43896 124 3.13083 125 3.60468 3.46997 127 3.64138 128 3.57986 129 3.62519 130 3.59504 131 3.84561 132 2.61849 133 3.59856 134 3.25213 135 3.67058 136 3.52251 137 3.6612 138 3.5978 139 3.62104 140 3.59684 141 4.32906 142 3.00288 143 3.77489 144 3.37528 145 3.74394 146 3.57397 147 3.67824 148 3.61194 149 3.60886 150 3.59448 151 5.28127 152 3.09713 153 4.0445 154 3.53619 155 3.75192 3.65883 5.94434 5.45295 5.53137 5.28947 5.31283 5.18182 9.03003 5.88707 6.62041 5.65002 5.86777 5.41321 5.48772 5.26789 5.28721 5.18056 8.33433 5.63182 6.45477 5.55872 5.78572 5.36674 5.43874 5.24537 5.25936 5.18319 7.42576 5.37252 6.27767 5.46485 5.70159 5.31962 5.38851 5.22448 5.23183 5.18899 7.01391 4.64656 6.49121 5.33975 5.70552 5.30721 -0.22059 -0.17523 -0.08039 -0.06326 -0.00774 -0.00127 -2.36105 -1.3625 -0.72417 -0.49204 -0.29402 -0.20583 -0.10402 -0.06675 -0.00848 0.00273 -2.68619 -1.50297 -0.83966 -0.55143 -0.34469 -0.22278 -0.11564 -0.06512 -0.00382 0.00719 -2.85535 -1.5645 -0.9262 -0.58384 -0.37467 -0.22614 -0.11544 -0.05909 0.00583 0.00723 -2.88107 -1.73035 -1.12577 -0.72912 -0.44504 -0.29613
118 121 126 156
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
3.61833 5.38867 0.03835 157 3.71432 5.41048 -0.16756 3.57808 5.29017 0.0115 158 3.58823 5.17503 -0.04608 3.61283 5.30175 0.00764 159 3.61496 5.25632 -0.05106 3.58663 5.22826 -0.00261 160 3.58762 5.16701 0.01004 4.21473 9.8158 2.22537 161 6.21493 6.26104 -3.01802 2.42229 5.8185 0.987 162 3.51262 4.20287 -1.8508 3.82883 6.28734 0.65128 163 4.33064 5.95316 -1.33694 3.50511 5.58054 0.26832 164 3.66391 5.19813 -0.93025 3.6004 5.6123 0.11372 165 3.89321 5.64308 -0.65527 3.52987 5.40703 0.08234 166 3.78201 5.41118 -0.47558 3.61391 5.43504 0.02974 167 3.80308 5.3759 -0.29088 3.57124 5.30273 0.01108 168 3.68293 5.26292 -0.15298 3.61562 5.31752 0.00266 169 3.68533 5.24848 -0.08021 3.58857 5.22979 -0.0067 170 3.51327 5.18582 0.01164 2.97876 9.48635 1.27596 171 7.26949 5.26845 -2.80789 1.8406 6.15114 0.77317 172 4.09698 3.58242 -1.84285 3.57596 6.72959 0.60953 173 4.80236 5.29843 -1.40426 3.27734 5.67773 0.28608 174 3.95591 4.84573 -1.07792 3.60831 5.78806 0.14486 175 4.16891 5.36787 -0.81196 3.52354 5.4492 0.07183 176 3.9744 5.31295 -0.6654 3.6128 5.47895 0.01318 177 3.98631 5.31576 -0.44038 3.56604 5.31474 0.00582 178 3.94019 5.36721 -0.35226 3.61859 5.33225 -0.00586 179 3.88265 5.20233 -0.11284 3.58885 5.22547 -0.00787 180 3.62374 5.25586 -0.05292 2.36784 9.30738 0.37817 181 8.1626 4.09137 -2.12847 1.43977 6.36995 0.35306 182 4.77553 2.8491 -1.57997 3.30145 6.99051 0.28584 183 5.42181 4.50229 -1.1538 2.99503 5.77702 0.19141 184 4.39047 4.27724 -1.04166 3.55917 5.96507 0.08166 185 4.6389 4.91955 -0.75433 3.45411 5.48431 0.04974 186 4.265 4.93601 -0.72428 3.615 5.53794 -0.00633 187 4.35463 5.16311 -0.49243 3.56018 5.3214 -0.00248 188 4.26123 5.21763 -0.48194 3.62264 5.34223 -0.0164 189 4.25748 5.1873 -0.21406 3.58761 5.21609 -0.00778 190 4.19994 5.26706 -0.20657 2.21373 9.33302 -0.44417 191 8.57701 2.51314 -0.9816 1.2316 6.38549 -0.17434 192 5.37033 2.0178 -0.87645 3.24933 7.05807 -0.17766 193 5.9796 3.74191 -0.41436 2.81273 5.80802 -0.01046 194 4.84317 3.62353 -0.58488 3.54674 6.05269 -0.07944 195 5.12208 4.42575 -0.27802 3.37375 5.49394 -0.01299 196 4.62749 4.41949 -0.41042 3.62475 5.57761 -0.05261 197 4.77764 4.85933 -0.21924 3.54332 5.32094 -0.01538 198 4.62855 4.87187 -0.27435 3.61832 5.34594 -0.02741 199 4.72209 5.15086 -0.20071 3.58079 5.20034 -0.00968 200 4.90032 5.10888 -0.11927 ********************** THE END **********************
五．分析讨论
1）正确性验证：
A 水平外边界上单元水平应力与水平分布荷载相近（误差不大于0.5MPa）。
单元编号 水平应力 分布荷载 误差值
149 150 159 160 169 170 179 180 189 190 199 200 5.23546 5.19128 5.25927 5.16532 5.24822 5.18474 5.19966 5.25808 5.18549 5.27239 5.14963 5.10778 5.2 5.2 5.2 5.2 5.2 5.2 5.2 5.2 5.2 5.2 5.2 5.2 0.03546 -0.00872 0.05927 -0.03468 0.04822 -0.01526 -0.00034 0.05808 -0.01451 0.07239 -0.05037 -0.09222
B 竖直上边界上单元竖向应力与竖向分布荷载相近（误差不大于0.5MPa）。
单元编号
49 50 59 60 69 70 79 80 89 90 99 100 109 110 119 120 129 130 139 140 149 150 160
竖向应力 3.60973 3.5926 3.61072 3.58976 3.61115 3.58637 3.59949 3.57528 3.62626 3.58121 3.6301 3.59287 3.63032 3.59512 3.62807 3.59936 3.63794 3.59908 3.6286 3.59179 3.6116 3.58933 3.5763
分布荷载 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6 3.6
误差值 0.00973 -0.0074 0.01072 -0.01024 0.01115 -0.01363 -0.00051 -0.02472 0.02626 -0.01879 0.0301 -0.00713 0.03032 -0.00488 0.02807 -0.00064 0.03794 -0.00092 0.0286 -0.00821 0.0116 -0.01067 -0.0237
C 静力平衡：计算模型竖向约束反力之和等于作用于模型上表面竖向荷载之和（误差不大于0.005MN）。
D 静力平衡：计算模型水平约束反力之和等于作用于模型右表面 水平荷载之和（误差不大于0.005MN）。
NODE 1 2 3 4 5 6 12 18 24 30 36 121 122 123 124 125 126
约束反力之和 竖向荷载之和
误差
X-COMP -6.29074 -12.56193 -13.70171 -15.30006 -17.60687 -9.50883
0 0 0 0 0 -5.56073 -8.4903 -8.84377 -9.93748 -11.74792 -7.16952 -126.71986 126.72 0.00014
Y-COMP 0 0 0 0 0
-13.68021 -27.6862 -28.36732 -28.69163 -28.69593 -14.31778
0 0 0 0 0 0
-141.43907 141.44 0.00093
2)分析应力、位移分布：
A 绘出邻近水平轴各单元水平应力与竖向应力对单元横向坐标的关系曲线（Excel）。
邻近水平轴单元有91-100,101-110号单元，共两个图
B 绘出邻近竖向轴各单元水平应力与竖向应力对单元竖向坐标的关系曲线（Excel）。
邻近竖向轴单元有1-10,191-200号单元，共两个图，
C 绘出邻近水平轴各节移与节点横向坐标的关系曲线（Excel）。 邻近水平轴节点共6个61-66，
D 绘出邻近竖直轴各节点竖向位移与节点竖向坐标的关系曲线（Excel）。
竖直轴节点共两部分，1-6,121-126，所以有两个图。
六.分析讨论: