[求助]FLAC3D自重应力不平衡，急求指教续！！！ [复制链接]

gen zon brick p0 -30 0 -17.25 p1 -15.3 0 -17.25 p2 -30 79 -17.25 p3 -30 0 -12.75 size 5 16 3 group 第六层土               TB31- ()  
gen zon brick p0 -15.3 0 -17.25 p1 -5 0 -17.25 p2 -15.3 79 -17.25 p3 -15.3 0 -12.75 size 10 16 3 group 第六层土             La[V$+Y  
gen zon brick p0 -5 0 -17.25 p1 -2 0 -17.25 p2 -5 79 -17.25 p3 -5 0 -12.75 size 1 16 3 group 第六层土               YGNP53CU  
gen zon brick p0 -2 0 -17.25 p1 8.3 0 -17.25 p2 -2 79 -17.25 p3 -2 0 -12.75 size 10 16 3 group 第六层土               |`FY1NN   
gen zon brick p0 -30 79 -17.25 p1 -15.3 79 -17.25 p2 -30 85 -17.25 p3 -30 79 -12.75 size 5 2 3 group 第六层土               N8df8=.kw  
gen zon brick p0 -15.3 79 -17.25 p1 -5 79 -17.25 p2 -15.3 85 -17.25 p3 -15.3 79 -12.75 size 10 2 3 group 第六层土           KMax$  
gen zon brick p0 -5 79 -17.25 p1 -2 79 -17.25 p2 -5 85 -17.25 p3 -5 79 -12.75 size 1 2 3 group 第六层土               FGzwhgy  
gen zon brick p0 -2 79 -17.25 p1 8.3 79 -17.25 p2 -2 85 -17.25 p3 -2 79 -12.75 size 10 2 3 group 第六层土               <N~K;n v  
gen zon brick p0 -30 85 -17.25 p1 -15.3 85 -17.25 p2 -30 138 -17.25 p3 -30 85 -12.75 size 5 11 3 group 第六层土             7b+6%fV  
gen zon brick p0 -15.3 85 -17.25 p1 -5 85 -17.25 p2 -15.3 138 -17.25 p3 -15.3 85 -12.75 size 10 11 3 group 第六层土         wUJcmM;  
gen zon brick p0 -5 85 -17.25 p1 -2 85 -17.25 p2 -5 138 -17.25 p3 -5 85 -12.75 size 1 11 3 group 第六层土               hM
!a_'  
gen zon brick p0 -2 85 -17.25 p1 8.3 85 -17.25 p2 -2 138 -17.25 p3 -2 85 -12.75 size 10 11 3 group 第六层土               r5^eNg k  
gen zon brick p0 8.3 0 -17.25 p1 10.8 0 -17.25 p2 8.3 79 -17.25 p3 8.3 0 -12.75 size 1 16 3 group 第六层土               5|)W.*Q  
gen zon brick p0 8.3 79 -17.25 p1 10.8 79 -17.25 p2 8.3 85 -17.25 p3 8.3 79 -12.75 size 1 2 3 group 第六层土               k+*u/neh  
gen zon brick p0 8.3 85 -17.25 p1 10.8 85 -17.25 p2 8.3 138 -17.25 p3 8.3 85 -12.75 size 1 11 3 group 第六层土             =7UsVn#o  
gen zon brick p0 10.8 0 -17.25 p1 49.8 0 -17.25 p2 10.8 79 -17.25 p3 10.8 0 -12.75 size 26 16 3 group 第六层土             "" EQE>d  
gen zon brick p0 10.8 79 -17.25 p1 16.8 79 -17.25 p2 10.8 85 -17.25 p3 10.8 79 -12.75 size 4 4 3 group 第六层土             ^S; -fYW2  
gen zon brick p0 40.8 79 -17.25 p1 49.8 79 -17.25 p2 40.8 85 -17.25 p3 40.8 79 -12.75 size 6 4 3 group 第六层土             Tw<q,O  
gen zon brick p0 22.8 79 -17.25 p1 34.8 79 -17.25 p2 22.8 85 -17.25 p3 22.8 79 -12.75 size 8 4 3 group 第六层土             2GG2jky{/  
gen zon brick p0 10.8 85 -17.25 p1 49.8 85 -17.25 p2 10.8 138 -17.25 p3 10.8 85 -12.75 size 26 11 3 group 第六层土         6_B]MN!(  
gen zon brick p0 -30 0 -35.25 p1 -15.3 0 -35.25 p2 -30 79 -35.25 p3 -30 0 -17.25 size 5 16 9 group 第七层土               xskz)kk  
gen zon brick p0 -15.3 0 -35.25 p1 -5 0 -35.25 p2 -15.3 79 -35.25 p3 -15.3 0 -17.25 size 10 16 9 group 第七层土             R=\IEqqsi  
gen zon brick p0 -5 0 -35.25 p1 -2 0 -35.25 p2 -5 79 -35.25 p3 -5 0 -17.25 size 1 16 9 group 第七层土               n+M<\  
gen zon brick p0 -2 0 -35.25 p1 8.3 0 -35.25 p2 -2 79 -35.25 p3 -2 0 -17.25 size 10 16 9 group 第七层土               ~a2}(]  
gen zon brick p0 -30 79 -35.25 p1 -15.3 79 -35.25 p2 -30 85 -35.25 p3 -30 79 -17.25 size 5 2 9 group 第七层土               ]6j{@z?{  
gen zon brick p0 -15.3 79 -35.25 p1 -5 79 -35.25 p2 -15.3 85 -35.25 p3 -15.3 79 -17.25 size 10 2 9 group 第七层土           N?8!3&TiV  
gen zon brick p0 -5 79 -35.25 p1 -2 79 -35.25 p2 -5 85 -35.25 p3 -5 79 -17.25 size 1 2 9 group 第七层土               '~ 47)fN  
gen zon brick p0 -2 79 -35.25 p1 8.3 79 -35.25 p2 -2 85 -35.25 p3 -2 79 -17.25 size 10 2 9 group 第七层土               f _:A0  
gen zon brick p0 -30 85 -35.25 p1 -15.3 85 -35.25 p2 -30 138 -35.25 p3 -30 85 -17.25 size 5 11 9 group 第七层土             .T`%tJ-Em  
gen zon brick p0 -15.3 85 -35.25 p1 -5 85 -35.25 p2 -15.3 138 -35.25 p3 -15.3 85 -17.25 size 10 11 9 group 第七层土         Zv{'MIv&v  
gen zon brick p0 -5 85 -35.25 p1 -2 85 -35.25 p2 -5 138 -35.25 p3 -5 85 -17.25 size 1 11 9 group 第七层土               iWR)ke  
gen zon brick p0 -2 85 -35.25 p1 8.3 85 -35.25 p2 -2 138 -35.25 p3 -2 85 -17.25 size 10 11 9 group 第七层土               wC'Szni  
gen zon brick p0 8.3 0 -35.25 p1 10.8 0 -35.25 p2 8.3 79 -35.25 p3 8.3 0 -17.25 size 1 16 9 group 第七层土               <F'\lA9  
gen zon brick p0 8.3 79 -35.25 p1 10.8 79 -35.25 p2 8.3 85 -35.25 p3 8.3 79 -17.25 size 1 2 9 group 第七层土               -mh3DhJ,  
gen zon brick p0 8.3 85 -35.25 p1 10.8 85 -35.25 p2 8.3 138 -35.25 p3 8.3 85 -17.25 size 1 11 9 group 第七层土             ~wdGd+ez  
gen zon brick p0 10.8 0 -35.25 p1 49.8 0 -35.25 p2 10.8 79 -35.25 p3 10.8 0 -17.25 size 26 16 9 group 第七层土             r<$y=B  
gen zon brick p0 10.8 79 -35.25 p1 16.8 79 -35.25 p2 10.8 85 -35.25 p3 10.8 79 -17.25 size 4 4 9 group 第七层土             Kc\fu3Q  
gen zon brick p0 40.8 79 -35.25 p1 49.8 79 -35.25 p2 40.8 85 -35.25 p3 40.8 79 -17.25 size 6 4 9 group 第七层土             M"L=L5OH-  
gen zon brick p0 22.8 79 -35.25 p1 34.8 79 -35.25 p2 22.8 85 -35.25 p3 22.8 79 -17.25 size 8 4 9 group 第七层土             {_*yGK48n  
gen zon brick p0 10.8 85 -35.25 p1 49.8 85 -35.25 p2 10.8 138 -35.25 p3 10.8 85 -17.25 size 26 11 9 group 第七层土         {oL>1h,%3?  
save 一直七层土.sav                                       1|=A*T-<M  
;------ 1#墩（分布四根桩） ------                                     {{!-Gr  
generate zone radcyl p0 18.3 80.5 12.15 p1 16.8 80.5 12.15 p2 18.3 80.5 10.15 p3 18.3 79 12.15 size 3 2 2 1 dim 0.75 0.75 0.75 0.75 rat 1 1 1 1.1 group 承台1     Dw"\/p:-3  
generate zone radcyl p0 18.3 80.5 10.15 p1 16.8 80.5 10.15 p2 18.3 80.5 -34.85 p3 18.3 79 10.15 size 3 45 2 1 dim 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩周土1     Q+{n-? :  
generate zone radcyl p0 18.3 80.5 -34.85 p1 16.8 80.5 -34.85 p2 18.3 80.5 -35.25 p3 18.3 79 -34.85 size 3 1 2 1 dim 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩下土1 fill 7zj{wp!  
Nz-&MS  
Q/Rqa5LI:  
group 桩下土1 );YDtGip J  
generate zone radcyl p0 18.3 80.5 12.15 p1 18.3 79 12.15 p2 18.3 80.5 10.15 p3 19.8 80.5 12.15 size 3 2 2 1 dim 0.75 0.75 0.75 0.75 rat 1 1 1 1.2 group 承台1     h{qgEIk&  
generate zone radcyl p0 18.3 80.5 10.15 p1 18.3 79 10.15 p2 18.3 80.5 -34.85 p3 19.8 80.5 10.15 size 3 45 2 1 dim 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩周土1     0> \sQ,T  
generate zone radcyl p0 18.3 80.5 -34.85 p1 18.3 79 -34.85 p2 18.3 80.5 -35.25 p3 19.8 80.5 -34.85 size 3 1 2 1 dim 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩下土1 fill 
:k#HW6p  
eyxW 0}[  
#<x
m.  
group 桩下土1 #O&8A  
generate zone  radcyl p0 18.3 80.5 12.15 p1 19.8 80.5 12.15 p2 18.3 80.5 10.15 p3 18.3 82 12.15 size 3 2 2 1 DIM 0.75 0.75 0.75 0.75 rat 1 1 1 1.1 group 承台1     |w3M7;~eF  
generate zone  radcyl p0 18.3 80.5 10.15 p1 19.8 80.5 10.15 p2 18.3 80.5 -34.85 p3 18.3 82 10.15 size 3 45 2 1 DIM 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩周土1     t_1LL >R  
generate zone  radcyl p0 18.3 80.5 -34.85 p1 19.8 80.5 -34.85 p2 18.3 80.5 -35.25 p3 18.3 82 -34.85 size 3 1 2 1 DIM 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩下土1 fill gRzxLf`K  
/x *3}oI  
(cO:`W6.  
group 桩下土1 3XNCAb2  
generate zone  radcyl p0 18.3 80.5 12.15 p1 18.3 82 12.15 p2 18.3 80.5 10.15 p3 16.8 80.5 12.15 size 3 2 2 1 DIM 0.75 0.75 0.75 0.75 rat 1 1 1 1.1 group 承台1    
B33\?Yj)  
generate zone  radcyl p0 18.3 80.5 10.15 p1 18.3 82 10.15 p2 18.3 80.5 -34.85 p3 16.8 80.5 10.15 size 3 45 2 1 DIM 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩周土1     N2o7%gJw  
generate zone  radcyl p0 18.3 80.5 -34.85 p1 18.3 82 -34.85 p2 18.3 80.5 -35.25 p3 16.8 80.5 -34.85 size 3 1 2 1 DIM 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩下土1 fill K(|}dl:  
*m(=V1"  
C,eu9wOT  
group 桩下土1 4skD(au8  
generate zone cylinder p0 18.3 80.5 62.15 p1 19.05 80.5 62.15 p2 18.3 80.5 15.15 p3 18.3 81.25 62.15 size 3 47 2 group 桩1 nJ;.Td  
generate zon reflect origin 18.3 80.5 12.15 normal 1 0 0 range grou 桩1                           y
f,z$CR  
generate zon reflect origin 18.3 80.5 12.15 normal 0 1 0 range grou 桩1                           m4Zk\,1m.|  
generate zone reflect origin 19.8 82 12.15 normal 0 1 0 range group 承台1 any group 桩周土1 any group 桩下土1 any group 桩1 any                 qxc[M8s  
generate zone reflect origin 19.8 82 12.15 normal 1 0 0 range group 承台1 any group 桩周土1 any group 桩下土1 any group 桩1 any                 -nwypu  
;------ 2#墩（分布四根桩） ------                                     }6ldjCT/,  
generate zone radcyl p0 36.3 80.5 12.15 p1 34.8 80.5 12.15 p2 36.3 80.5 10.15 p3 36.3 79 12.15 size 3 2 2 1 dim 0.75 0.75 0.75 0.75 rat 1 1 1 1.1 group 承台2     qe\5m.k  
generate zone radcyl p0 36.3 80.5 10.15 p1 34.8 80.5 10.15 p2 36.3 80.5 -34.85 p3 36.3 79 10.15 size 3 45 2 1 dim 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩周土2     % ]U  
generate zone radcyl p0 36.3 80.5 -34.85 p1 34.8 80.5 -34.85 p2 36.3 80.5 -35.25 p3 36.3 79 -34.85 size 3 1 2 1 dim 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩下土2 fill $/ ]
,tSm  
n=q76W\  
NHE18_v5  
group 桩下土2 5bb(/YtFy  
generate zone radcyl p0 36.3 80.5 12.15 p1 36.3 79 12.15 p2 36.3 80.5 10.15 p3 37.8 80.5 12.15 size 3 2 2 1 dim 0.75 0.75 0.75 0.75 rat 1 1 1 1.2 group 承台2     ~V6D<  
generate zone radcyl p0 36.3 80.5 10.15 p1 36.3 79 10.15 p2 36.3 80.5 -34.85 p3 37.8 80.5 10.15 size 3 45 2 1 dim 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩周土2     5mR 1@  
generate zone radcyl p0 36.3 80.5 -34.85 p1 36.3 79 -34.85 p2 36.3 80.5 -35.25 p3 37.8 80.5 -34.85 size 3 1 2 1 dim 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩下土2 fill NxILRKwO  
J .<F"r>  
`d(ThP;g  
group 桩下土2 |V(0GB  
generate zone  radcyl p0 36.3 80.5 12.15 p1 37.8 80.5 12.15 p2 36.3 80.5 10.15 p3 36.3 82 12.15 size 3 2 2 1 DIM 0.75 0.75 0.75 0.75 rat 1 1 1 1.1 group 承台2     2fS:- 8N  
generate zone  radcyl p0 36.3 80.5 10.15 p1 37.8 80.5 10.15 p2 36.3 80.5 -34.85 p3 36.3 82 10.15 size 3 45 2 1 DIM 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩周土2     yt2PU_),  
generate zone  radcyl p0 36.3 80.5 -34.85 p1 37.8 80.5 -34.85 p2 36.3 80.5 -35.25 p3 36.3 82 -34.85 size 3 1 2 1 DIM 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩下土2 fill vih9KBT  
RM/ 0A|  
q,%st~  
group 桩下土2 fN2lLn9/u  
generate zone  radcyl p0 36.3 80.5 12.15 p1 36.3 82 12.15 p2 36.3 80.5 10.15 p3 34.8 80.5 12.15 size 3 2 2 1 DIM 0.75 0.75 0.75 0.75 rat 1 1 1 1.1 group 承台2     1Z&(6cDY8M  
generate zone  radcyl p0 36.3 80.5 10.15 p1 36.3 82 10.15 p2 36.3 80.5 -34.85 p3 34.8 80.5 10.15 size 3 45 2 1 DIM 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩周土2     y1#1Ne_  
generate zone  radcyl p0 36.3 80.5 -34.85 p1 36.3 82 -34.85 p2 36.3 80.5 -35.25 p3 34.8 80.5 -34.85 size 3 1 2 1 DIM 0.75 0.75 0.75 0.75 rat 1 1 1 1 group 桩下土2 fill W*Y/l~x}  
 L
"aeG  
$:^td/p J  
group 桩下土2 \{D" !e  
generate zone cylinder p0 36.3 80.5 62.15 p1 37.05 80.5 62.15 p2 36.3 80.5 15.15 p3 36.3 81.25 62.15 size 3 47 2 group 桩2 /mZE/>&~,  
generate zon reflect origin 36.3 80.5 12.15 normal 1 0 0 range grou 桩2                           VxBo1\'  
generate zon reflect origin 36.3 80.5 12.15 normal 0 1 0 range grou 桩2                           :23P!^Y  
generate zone reflect origin 37.8 82 12.15 normal 0 1 0 range group 承台2 any group 桩周土2 any group 桩下土2 any group 桩2 any                 2Khv>#l  
generate zone reflect origin 37.8 82 12.15 normal 1 0 0 range group 承台2 any group 桩周土2 any group 桩下土2 any group 桩2 any                 !5N.B|Nt  
;---------设定桩及承台网格单元接触面---------                                       QZ8IV>  
;----1#墩                                       }-2|XD%]  
interface 1 face range  cylinder end1 18.3 80.5 12.15 end2 18.3 80.5 -34.86 radius 0.76 any &                       -Qe'Y
By:  
    cylinder end1 18.3 80.5 10.16 end2 18.3 80.5 -34.86 radius 0.74 not                         |':{lH6+1  
interface 2 face range  cylinder end1 18.3 80.5 -34.84 end2 18.3 80.5 -34.86 radius 0.76                         Uw:"n]G]D?  
interface 3 face range  cylinder end1 21.3 80.5 10.16 end2 21.3 80.5 -34.86 radius 0.76 any &                       % AgUUn&k  
    cylinder end1 21.3 80.5 10.16 end2 21.3 80.5 -34.86 radius 0.74 not                          0+8e,  
interface 4 face range  cylinder end1 21.3 80.5 -34.84 end2 21.3 80.5 -34.86 radius 0.76                         'N(R_q6MW  
interface 5 face range  cylinder end1 18.3 83.5 10.16 end2 18.3 83.5 -34.86 radius 0.76 any &                       |vC~HJpuv'  
    cylinder end1 18.3 83.5 10.16 end2 18.3 83.5 -34.86 radius 0.74 not                         G+m }MOQP7  
interface 6 face range  cylinder end1 18.3 83.5 -34.84 end2 18.3 83.5 -34.86 radius 0.76                         _u Il  
interface 7 face range  cylinder end1 21.3 83.5 10.16 end2 21.3 83.5 -34.86 radius 0.76 any &                       u%!@(eKM-  
    cylinder end1 21.3 83.5 10.16 end2 21.3 83.5 -34.86 radius 0.74 not                         xYB{;K  
interface 8 face range  cylinder end1 21.3 83.5 -34.84 end2 21.3 83.5 -34.86 radius 0.76                         'c~4+o4co  
;----2#墩                                       nr3==21Om4  
interface 9 face range  cylinder end1 36.3 80.5 12.15 end2 36.3 80.5 -34.86 radius 0.76 any &                       $pz/?>!  
    cylinder end1 36.3 80.5 10.16 end2 36.3 80.5 -34.86 radius 0.74 not                         2&5K.Ui%  
interface 10 face range  cylinder end1 36.3 80.5 -34.84 end2 36.3 80.5 -34.86 radius 0.76                         +cRn%ioVi  
interface 11 face range  cylinder end1 39.3 80.5 10.16 end2 39.3 80.5 -34.86 radius 0.76 any &                       H,NF;QPPC  
    cylinder end1 39.3 80.5 10.16 end2 39.3 80.5 -34.86 radius 0.74 not                         !@"OB~  
interface 12 face range  cylinder end1 39.3 80.5 -34.84 end2 39.3 80.5 -34.86 radius 0.76                         ji="DYtL  
interface 13 face range  cylinder end1 36.3 83.5 10.16 end2 36.3 83.5 -34.86 radius 0.76 any &                       rZpXPI  
    cylinder end1 36.3 83.5 10.16 end2 36.3 83.5 -34.86 radius 0.74 not                         R@2X3s:  
interface 14 face range  cylinder end1 36.3 83.5 -34.84 end2 36.3 83.5 -34.86 radius 0.76                         QsW/X0YBv  
interface 15 face range  cylinder end1 39.3 83.5 10.16 end2 39.3 83.5 -34.86 radius 0.76 any &                       1E$|~   
    cylinder end1 39.3 83.5 10.16 end2 39.3 83.5 -34.86 radius 0.74 not                         1 TXioDs=_  
interface 16 face range  cylinder end1 39.3 83.5 -34.84 end2 39.3 83.5 -34.86 radius 0.76                         wgA_38To  
;-----1#墩-----                                       "Y.y:Vv;  
gen zone cylinder p0 19.8 82 12.15 p1 20.55 82 12.15 p2 19.8 82 15.65 p3 19.8 81.25 12.15 size 3 3 4 group 桥墩1           y)<q/  
gen zone radcyl p0 19.8 82 12.15 p1 22.8 82 12.15 p2 19.8 82 12.75 p3 19.8 79 12.15 size 3 1 4 3 &               OZ&o:/*HM  
             dim 0.75 0.75 0.75 0.75 rat 1 1 1 1.1 group 桥墩1周围土               2A!FDr~cdT  
gen zone reflect origin 19.8 82 12.15 normal 0 1 0 range group 桥墩1 any  group 桥墩1周围土 any                       e}W)LPR!  
gen zone reflect origin 19.8 82 12.15 normal 1 0 0 range group 桥墩1 any  group 桥墩1周围土 any                       ]_$[8#kg  
;-----2#墩-----                                       phz&zlD  
gen zone cylinder p0 37.8 82 12.15 p1 38.55 82 12.15 p2 37.8 82 15.65 p3 37.8 81.25 12.15 size 3 3 4 group 桥墩2           w2'5#`m  
gen zone radcyl p0 37.8 82 12.15 p1 40.8 82 12.15 p2 37.8 82 12.75 p3 37.8 79 12.15 size 3 1 4 3 &               .S4u-  
             dim 0.75 0.75 0.75 0.75 rat 1 1 1 1.1 group 桥墩2周围土               5-A\9UC*@  
gen zone reflect origin 37.8 82 12.15 normal 0 1 0 range group 桥墩2 any  group 桥墩2周围土  any                     oL<St$1  
gen zone reflect origin 37.8 82 12.15 normal 1 0 0 range group 桥墩2 any  group 桥墩2周围土  any                     _VXN#
@y  
save 最终模型.sav                                       KY^
Z  
ini z add -50 range gr 桩1 any gr 桩2 any                             }GIt!PG  
interface 1 prop ks 2e9 kn 2e9 fric 15                               "wc<B4"  
interface 2 prop ks 2e9 kn 2e9 fric 15                               Yr|4Fl~U  
interface 3 prop ks 2e9 kn 2e9 fric 15                               tl>7^hH  
interface 4 prop ks 2e9 kn 2e9 fric 15                               {c0`Um3&>  
interface 5 prop ks 2e9 kn 2e9 fric 15                              
S|}L&A  
interface 6 prop ks 2e9 kn 2e9 fric 15                               4P
o_-4  
interface 7 prop ks 2e9 kn 2e9 fric 15                                AOx[  
interface 8 prop ks 2e9 kn 2e9 fric 15                               C9;kpqNG#u  
interface 9 prop ks 2e9 kn 2e9 fric 15                               S8gs-gL#Og  
interface 10 prop ks 2e9 kn 2e9 fric 15                               6MI8zRX  
interface 11 prop ks 2e9 kn 2e9 fric 15                               d d;T-wa}  
interface 12 prop ks 2e9 kn 2e9 fric 15                               8b=_Y;  
interface 13 prop ks 2e9 kn 2e9 fric 15                               %jM,W}2  
interface 14 prop ks 2e9 kn 2e9 fric 15                               5LMw?P.<  
interface 15 prop ks 2e9 kn 2e9 fric 15                               f *)Z)6E  
interface 16 prop ks 2e9 kn 2e9 fric 15                               LH6vLuf  
;-------------材料参数-------------                                       Q59W#e)  
model mohr                                        =BrRYA  
;第一层土                                       t$ *0{w E  
prop bulk=2.5e6 shear=4e6 cohesion=3.8e4 friction=11.3 density=1870 tension=5.8e4 range z 11.75 12.75                       T^q 0'#/  
;第二层土                                       .G\7cZ  
prop bulk=4.9e6 shear=5e6 cohesion=5.6e4 friction=16 density=1580 tension=7.2e4 range z 8.75 11.75                         jj>]9z  
prop bulk=4.9e6 shear=5e6 cohesion=5.6e4 friction=16 density=1580 tension=7.2e4 range group 燃气管外土体                   :
E?V.  
;第三层土Ⅱ                                       Ir]
\|t  
prop bulk=5.5e6 shear=5.1e6 cohesion=5e4 friction=20 density=1930 tension=7.9e4 range z 0 8.75                             Vw"\
{`  
;第四层土                                       zW nR6*\  
prop bulk=5.2e7 shear=4.4e6 cohesion=6e4 friction=24 density=1860 tension=7.9e4 range z -3.15 0                             p b,. r  
;第五层土                                       ?h2}#wg  
prop bulk=5.9e6 shear=5.6e6 cohesion=5.3e4 friction=21.5 density=1770 tension=7.9e4 range z -12.75 -3.15                   :v 4]D4\o  
prop bulk=4.3e7 shear=4.4e6 cohesion=4e4 friction=24 density=1860 tension=7.9e4 range gr 左线上洞内土体1 any gr 左线上洞内土体2 &                           8;X-)&R  
  gr 左线下洞内土体1 any gr 左线下洞内土体2 any gr 右线上洞内土体1 any &                             paMa+jhQQ  
  gr 右线上洞内土体2 any gr 右线下洞内土体1 any gr 右线下洞内土体2 any                             y
+q5UC|  
;第六层土                                       FgO)DQm  
prop bulk=6.8e6 shear=7.85e6 cohesion=7.5e4 friction=20.5 density=1620 tension=5.6e4 range z -17.25 -12.75                 XX~,>Q}H=  
;第七层土                                       #fM'>$N  
prop bulk=4.5e6 shear=4.56e6 cohesion=5.3e4 friction=24.5 density=2050 tension=7.9e4 range z -35.25 -17.25                 ch]29  
;燃气管                                       ,u!sjx  
prop bulk=1.75e11 shear=8.1e10 density=7800  range gr 燃气管                                 hv+zGID7  
  ;混泥土                                       aQ~s`^D  
prop bulk=1.67e4 shear=7.69e3 density=2400 range &                                   ;wD)hNLAvR  
gr 桩1 any gr 桩2 any &                                 D)Dr__x  
gr 桥墩1 any gr 桥墩2 any &                                 %XTI-B/K  
gr 承台1 any gr 承台2 any                                 !!ya  
apply szz -8e5 range group 桥墩1 any gr 桥墩2 any                               x)VJFuqy  
;施加边界条件                                       .wr>]yN  
set gravity 0 0 -9.8                                   =\d?'dII:  
attach face                                       Q@HV- (A  
;位移边界条件                                       Xm&L BX  
fix z range z -35.26 -35.24                                   Y\tui+?J  
fix x range x -30.01 -29.99                                   }
~q5w{_n  
fix x range x 49.79 49.81                                   !&\INl-Z  
fix y range y -0.99 1.01                                   BerwI 7!=  
fix y range y 137.99 138.01                                   i6N',&jFU  
solve                                       J8D,ZfPN`d  
save 自重应力求解.sav                                   S tyf
B  
  LACC3D

好长的帖子，

不平衡很多时候是由于模型参数设置问题以及网格没有建好（plot sk一下层与层之间有黑面，表示层之间节点没对应，解决办法是改网格size或在两层之间加一个interface）造成的

不平衡很多时候是由于模型参数设置问题以及网格没有建好（plot sk一下层与层之间有黑面，表示层之间节点没对应，解决办法是改网格size或在两层之间加一个interface）造成的

不平衡很多时候是由于模型参数设置问题以及网格没有建好（plot sk一下层与层之间有黑面，表示层之间节点没对应，解决办法是改网格size或在两层之间加一个interface）造成的