[求助]求助  我这老是不收敛 请问是哪里的问题？ [复制链接]

Pt+_0
OsR  
new kn.z8%^(  
config dyn edQ><lz  
set dyn off f WXzK<  
;=============================================== =
g:\R$lQ  
;建模型 p\K5B,  
;=============================================== VbYapPu4b!  
def setup              ;定义函数setup >smaR^m  
   yc=1               ;y方向尺寸 _?"J.
i  
   czz1=-13.16 JiR|+6"7  
   czz11=1.9            ;四分之一车站z方向尺寸   yrX]w3kr%  
   czz12=3            ;四分之一车站z方向尺寸 l?;S>s*\?  
   czz13=6.8            ;四分之一车站z方向尺寸 K]l)z* I  
   czz14=1.46            ;四分之一车站z方向尺寸     nDvWOt  
   ttx1=10.05              ;四分之一车站外围土体x方向尺寸 plq\D.C  
   ttx2=-9.55              ;四分之一车站外围土体-x方向尺寸   u[DV{o  
   ht1=-3.1 xT/&'$@{)  
   ht2=3.1 j3><J  
   tt2=68.9 ;78.95             b+>godTi_  
   tt1=-68.9 LmE-&  
   xh=33.74              ；车站以下土体高度 a=R-F!P)  
   mVLGQlvVK  
end @[Th{HTc.G  
setup   BJ5#!I%h  
;------------------------------------------------------------------------------------------------- <PxEl4  
; 一、生成右侧土体 `g~-5Z~J  
;------------------------------------------------------------------------------------------------- QZfnoKz  
def gen_Sando1 AXCJFqk;  
  command mqL+W  
gen zone brick      p0 0,0,ht1   p1 add ttx1,0,0  p2 add 0,yc,0  p3 add 0,0,ht2 & hGeRM4zVZZ  
                    p4 add ttx1,yc,0  p5 add 0,yc,ht2  p6 add ttx1,0,ht2  p7 add ttx1,yc,ht2 & gNeCnf#Xa  
                    size 10 1 3                               ;ratio 1 1 1 1.2 |Th{*IJ<,  
rgCId@R  
        end_command g
nGw7V  
end P$bo8*  
gen_Sando1 ra1hdf0"  
group sand_o1 range x 0 10.05  y 0 1  z -3.1 0 EbQ}w"{  
;------------------------------------------------------------------------------------------------- W=*\4B]  
def gen_Sando1 NO1PGen  
  command ^BZdR<;  
gen zone brick      p0 ttx1,0,ht1   p1 add tt2,0,0  p2 add 0,yc,0  p3 add 0,0,ht2 & s5HbuyR^  
                    p4 add tt2,yc,0  p5 add 0,yc,ht2  p6 add tt2,0,ht2  p7 add tt2,yc,ht2 & ;1nd~0
o  
                    size 50 1 3 ;ratio 1.1 1 1 ?kSs7e>  
q,GL#L  
        end_command 21qhlkdc  
end (hoqLL\}k  
gen_Sando1 92i#It}-/  
group sand_o1 range x 10.05 78.95 y 0 1 z -3.1 0                   xjYFTb}!  
;------------------------------------------------------------------------------------------------- tj3p71%  
def gen_Sando2 ;z68`P-  
  command BG"6jQh  
gen zone brick      p0 ttx1,0,-5   p1 add tt2,0,0  p2 add 0,yc,0  p3 add 0,0,czz11 & *))
|ZE6jI  
                    p4 add tt2,yc,0  p5 add 0,yc,czz11  p6 add tt2,0,czz11  p7 add tt2,yc,czz11 &   EA\~m*k  
                    size 50 1 2 ;ratio 1.1 1 1 M<nn+vy`  
809-p_)B  
        end_command Z6p>R;9n  
end kAoai|m@R  
gen_Sando2 yT4|eHl  
group sand_o2 range x 10.05 78.95 y 0 1 z -5 -3.1                                    ;11.26 13.16 b%BwGS(z  
Evc 9k  
;g9+*$Gw  
;------------------------------------------------------------------------------------------------- sXWMXQ3  
;------------------------------------------------------------------------------------------------- ^RyrUb  
def gen_Sando3 qA30G~S  
  command ,x/j&S9!  
gen zone brick      p0 ttx1,0,-8     p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz12 & O_ cK4  
                    p4 add tt2,yc,0  p5 add 0,yc,czz12  p6 add tt2,0,czz12  p7 add tt2,yc,czz12 &   "'Q:%_;  
                    size 50 1 3 ;ratio 1.1 1 1 *K;~V  
I zT%Kq  
        end_command 2+.m44>Ti  
end k8TMdWW  
gen_Sando3 z!%}0  
group sand_o3 range x 10.05 78.95 y 0 1 z -8 -5                                          ;8.26 11.26 >&R|t_ypw  
;------------------------------------------------------------------------------------------------- *sTQ9 Kr  
def gen_Sando4 .JqIAC~  
  command ]:;gk&P  
gen zone brick      p0 ttx1,0,-14.8   p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz13 & ]%."  
                    p4 add tt2,yc,0   p5 add 0,yc,czz13  p6 add tt2,0,czz13  p7 add tt2,yc,czz13 &   ":Q^/;D}U  
                    size 50 1 6 ;ratio 1.1 1 1 &Lw| t_y  
<+a\'Xc  
        end_command [o~w>,a  
end e/6oC~#]  
gen_Sando4 9SPu 4i  
group sand_o4 range x 10.05 78.95 y 0 1 z -14.8 -8                                   ;1.46 8.26 3-05y!vbcE  
|Bid(`t.  
;------------------------------------------------------------------------------------------------- -[=`bHo  
def gen_Sando5 0czy:d,M%  
  command X:A\{^~  
gen zone brick      p0 ttx1,0,-16.26   p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz14 & :MVD83?4  
                    p4 add tt2,yc,0    p5 add 0,yc,czz14  p6 add tt2,0,czz14  p7 add tt2,yc,czz14 &   ?2M15Q  
                    size 50 1 2 ;ratio 1.1 1 1 g5.Z B@j  
c
1`o3gb  
        end_command ]WG\+1x9  
end TsQMwV_h  
gen_Sando5 <Wd$6  
group sand_o5 range x 10.05 78.95 y 0 1 z -16.26 -14.8                            ;0 1.46 2+enRR~  
;------------------------------------------------------------------------------------------------- E=y
#

~W  
;------------------------------------------------------------------------------------------------- L`\ILJz  
M@8(h=  
def gen_Sando6 6T-(GHzfHJ  
  command R?W8l5CIk  
gen zone brick      p0 0,0,-50        p1 add ttx1,0,0    p2 add 0,yc,0     p3 add 0,0,xh & #L"h>,b  
                    p4 add ttx1,yc,0  p5 add 0,yc,xh     p6 add ttx1,0,xh  p7 add ttx1,yc,xh &   j{vzCRa>8  
                    size 10 1 20 ;ratio 1.2 1 1 tua+R_"  
sYz:(hZS  
        end_command Ii)TCSt9U?  
end xASjw?  
gen_Sando6 S%4K-I  
group sand_o6 range x 0 10.05  y 0 1 z  -50  -16.26 D|xSO~M5  
;------------------------------------------------------------------------------------------------- nT>?}/S  
def gen_Sando6 &zZSWNW  
  command Oj:`r*z43  
gen zone brick      p0 ttx1,0,-50    p1 add tt2,0,0  p2 add 0,yc,0    p3 add 0,0,xh & ^%L$$V nG  
                    p4 add tt2,yc,0  p5 add 0,yc,xh  p6 add tt2,0,xh  p7 add tt2,yc,xh &   AXyuXB  
                    size 40 1 20 ;ratio 1.1 1   3eB2=_V`  
SG~R!kN}Q  
        end_command 46vz=# ,6L  
end X7G6y|4;w  
gen_Sando6 0ode&dB  
group sand_o6 range  x 10.05 78.95 y 0 1  z -50   -16.26                      ;-33.74 0 {XVSHUtw  
;------------------------------------------------------------------------------------------------- Q_dFZ  
;------------------------------------------------------------------------------------------------- eg3{sDv,  
; 二、生成左侧土体 P|\,kw>l  
;------------------------------------------------------------------------------------------------- 7G/"!ePW6`  
;------------------------------------------------------------------------------------------------- Y4_i=}\*vf  
def gen_Sand11 pO^ 6p%  
  command 5XhV+t g.  
gen zone brick      p0 0,0,ht1        p1 add 0,yc,0  p2 add ttx2,0,0      p3 add 0,0,ht2 & (<ejJPWT  
                    p4 add ttx2,yc,0  p5 add ttx2,0,ht2  p6 add 0,yc,ht2  p7 add ttx2,yc,ht2 & 4&\m!s  
                    size 1 10 3                               ;ratio 1 1 1 1.2 &"BKue~q@p  
@*oi1
_q  
        end_command ,FTF@h-Cs  
end loZfzN&6A  
gen_Sand11 */1z=  
group sand_o7 range x -9.55 0 y 0 1 z -3.1 0 Na=q(OKN  
;------------------------------------------------------------------------------------------------- JwO+Dd  
def gen_Sand2 ukw'$Yt2  
  command m*'#`vIbb  
gen zone brick      p0 ttx2,0,ht1    p1 add 0,yc,0   p2 add tt1,0,0     p3 add 0,0,ht2 & dL"v*3Fy  
                    p4 add tt1,yc,0  p5 add tt1,0,ht2  p6 add 0,yc,ht2  p7 add tt1,yc,ht2 & %63<Iz"  
                    size 1 50  3 ;ratio 1 1.1 1 ?$ 3=m)s  
43eGfp'  
        end_command b7$?'neH/.  
end gnv4.f:  
gen_Sand2 CB~&!MdMr  
group sand_o7 range x -78.45 -9.55 y 0 1 z -3.1 0                   [L8gG.wy  
;------------------------------------------------------------------------------------------------- cX*^PSM  
def gen_Sand31 3laSPih[.  
  command u^ T2  
gen zone brick      p0 ttx2,0,-5     p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz11 & G -;Yua2\  
                    p4 add tt1,yc,0  p5 add tt1,0,czz11  p6 add 0,yc,czz11  p7 add tt1,yc,czz11 &   T:si?7CR  
                    size 1 50 2 ;ratio 1 1.1 1 ]?kf;A@  
0<Y)yNsV  
        end_command ':Te#S  
end 0iYe>u  
gen_Sand31 Cc^t&Eg  
group sand_o8 range x -78.45 -9.55 y 0 1 z -5 -3.1                                    ;11.26 13.16 xZkLN5I{  
;------------------------------------------------------------------------------------------------- ~"-wS
Am  
def gen_Sand32 csC3Wm{v  
  command sB6UlX;b:  
gen zone brick      p0 ttx2,0,-8     p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz12 & Z5+0?X0i  
                    p4 add tt1,yc,0  p5 add tt1,0,czz12  p6 add 0,yc,czz12  p7 add tt1,yc,czz12 &   ''Hq-Ng  
                    size 1 50 3 ;ratio 1 1.1 1 ISl'g'o  
6ul34\;  
        end_command a^2?W
 
end pY2nv/  
gen_Sand32 \
^+sgg{  
group sand_o9 range x -78.45 -9.55 y 0 1  z -8 -5                                          ;8.26 11.26
6} 9A0  
;------------------------------------------------------------------------------------------------- xkX, l{6  
def gen_Sand33 O:#to  
  command htjJ0>
&  
gen zone brick      p0 ttx2,0,-14.8   p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz13 & 8Hhe&B  
                    p4 add tt1,yc,0   p5 add tt1,0,czz13  p6 add 0,yc,czz13  p7 add tt1,yc,czz13 &   kxKBI{L  
                    size 1 50 6 ;ratio 1 1.1 1 e0D;]  
'K0Y@y  
        end_command NmeTp?)m  
end 4U((dx*m  
gen_Sand33 QKVZ![Y!s  
group sand_o10 range x -78.45 -9.55 y 0 1  z -14.8 -8                                   ;1.46 8.26 ?.T=(-  
M4QMD;Ez  
;------------------------------------------------------------------------------------------------- )P$
IXA\  
def gen_Sand34 C}Khh`8@5.  
  command Nk7Q  
gen zone brick      p0 ttx2,0,-16.26   p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz14 & 1:,aFp>qr  
                    p4 add tt1,yc,0    p5 add tt1,0,czz14  p6 add 0,yc,czz14  p7 add tt1,yc,czz14 &   P"- ,^?6  
                    size 1 50 2 ;ratio 1 1.1 1 wj/r)rv E  
X\h]N  
        end_command tDi<n}  
end #hai3>9|B  
gen_Sand34 ?Z;knX\?J  
group sand_o11 range x -78.45 -9.55 y 0 1  z -16.26 -14.8                            ;0 1.46 Hi?],5,/  
;------------------------------------------------------------------------------------------------- 4[6A~iC_  
;------------------------------------------------------------------------------------------------- E_h9y  
'\9A78NV{;  
def gen_Sand4 Cc=`:ED+  
  command $rdA0%
;  
gen zone brick      p0 0,0,-50        p1 add 0,yc,0  p2 add ttx2,0,0     p3 add 0,0,xh & 9 Hm!B )Y  
                    p4 add ttx2,yc,0  p5 add ttx2,0,xh   p6 add 0,yc,xh  p7 add ttx2,yc,xh &   f[.]JC+,  
                    size 1 10 20 bC&_OU:  
UZ<!(g.  
        end_command _+UD>u{  
end _uRgKoiy  
gen_Sand4  ~d }-  
group sand_o12 range x -9.55 0 y 0 1 z -16.26 -50 W4Eo1 E  
;------------------------------------------------------------------------------------------------- L<E`~\C'  
def gen_Sand5 F Hv|6zUX  
  command bNqjjg  
gen zone brick      p0 ttx2,0,-50    p1 add 0,yc,0  p2 add tt1,0,0    p3 add 0,0,xh & `T-(g1:9  
                    p4 add tt1,yc,0  p5 add tt1,0,xh  p6 add 0,yc,xh  p7 add tt1,yc,xh &   `+<5QtD  
                    size 1 50 20 ;ratio 1 1.1 1 @A)gsDt9A  
pdE=9l'  
        end_command /_LUys/0  
end kJ~^ }o  
gen_Sand5 ~2pct
qMA  
group sand_o12 range x -78.45 -9.55 y 0 1 z -16.26 -50                        ;-33.74 0 MOj 0"x)  
;------------------------------------------------------------------------------------------------- >iq^Ts  
call Flacmodel_haitang.dat uy B ?-Y+  
;------------------------------------------------------------------------------------------------- RY*6TYX!  
plot add axes Tj.;\a|d  
;-------------------------------------------------------------------------------------- I3SLR  
;二、生成接触面的命令 BqR8%F  
;-------------------------------------------------------------------------------------- *u-
TNg  
interface 1 face range x -9.55 10.05 y 0 1 z -16.27 -16.25 D3I;5m`_  
interface 2 face range x -9.55 10.05 y 0 1 z -3.11 -3.09  yXDf;`J  
interface 3 face range x -9.56 -9.54 y 0 1 z -5 -3.1 nGRF<2!  
interface 4 face range x -9.56 -9.54 y 0 1 z -8  -5 GE"#.J4z  
interface 5 face range x -9.56 -9.54 y 0 1 z -14.8 -8 7OT}V}iP  
interface 6 face range x -9.56 -9.54 y 0 1 z -16.26 -14.8 tnp]wZ  
p<$z!|7m  
int  1  prop kn  1.570e9  ks  1.570e9  fric  10 coh 10000                   ;prop属性 kn法向刚度 ks剪切刚度 fric为摩擦角 coh内聚力 PPohpdd)  
int  2  prop kn  1.570e9  ks  1.570e9  fric  15 coh 55000 8(BLS{-"<  
int  3  prop kn  1.570e9  ks  1.570e9  fric  26 coh 8000 bzZEwMc6  
int  4  prop kn  1.570e9  ks  1.570e9  fric  15 coh 55000 bJ9>,,D  
int  5  prop kn  1.570e9  ks  1.570e9  fric  14 coh 24000 #fa~^]EM]  
int  6  prop kn  1.570e9  ks  1.570e9  fric  17 coh 56000 GwpJxiFgk  
ini y add -2.0 rang group 1 gP<
l  
;=============================================== aeSXHd?+(  
;初始应力计算程序 QtRKmry{
 
;=============================================== 4Jw0m#UN1  
;定义模型及性质 _$1W:!f4  
model  mohr    rang  group sand_o1                                           t.]oLG22r  
prop bulk 1.875e7   shear   0.967e7 coh 10000   fric 10 rang group sand_o1 ><$hFrR!  
model  mohr    rang  group sand_o2                                           I@[.W!w  
prop bulk 1.875e7   shear   0.967e7 coh 55000   fric 15 rang group sand_o2 f~E'0f_  
model  mohr    rang  group sand_o3                                           -0>@jfP^D  
prop bulk 2.068e7   shear   0.793e7 coh 8000   fric 26 rang group sand_o3 5:O"T  
model  mohr    rang  group sand_o4                                           N 4Dyec\  
prop bulk 1.875e7   shear   0.967e7 coh 55000   fric 15 rang group sand_o4 gllXJM^ -  
model  mohr    rang  group sand_o5                                           u%&zY97/  
prop bulk 1.396e7   shear   0.758e7 coh 24000   fric 14 rang group sand_o5 ecr pv+  
model  mohr    rang  group sand_o6                                           w;X-i.%`  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o6 qgu.c`GmW  
model  mohr    rang  group sand_o7                                           T9u/|OP  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o7 .>&kAf.  
model  mohr    rang  group sand_o8                                           B=9|g1e  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o8 "5vFa7y  
model  mohr    rang  group sand_o9                                           sB /*gO  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o9 #w#B'  
model  mohr    rang  group sand_o10                                           Fm*O&6W\@A  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o10 w
KwireOs  
model  mohr    rang  group sand_o11                                           s7=]!7QGS!  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o11     '*22j ]  
model  mohr    rang  group sand_o12                                           p#3P`I>ZrT  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o12 n2$*Z6.G  
lGs fs(  
mo elas                 rang group 1 *
F&C`]  
;prop bulk 20.8e9  shear 15.58e9  rang group 1   ;正确 /VJ[1o^  
prop bulk 22.8e9  shear 17.08e9  rang group 1   ;只把柱子等效 O10h(Wg  
ini dens 2500 rang group 1 \5J/?  
ini dens 1900 rang group sand_o1     ;18000 R,tR{| 8  
ini dens 1990 rang group sand_o2 aG,N>0k8  
ini dens 1910 rang group sand_o3 wWwY.}j  
ini dens 1990 rang group sand_o4 Y(<>[8S m  
ini dens 1970 rang group sand_o5 KaOS!e'  
ini dens 2010 rang group sand_o6 u+S
*D\p<`  
ini dens 2010 rang group sand_o7 HmQuRW  
ini dens 2010 rang group sand_o8 W[+E5I  
ini dens 2010 rang group sand_o9 5?QR  
ini dens 2010 rang group sand_o10 oZ!rK/qoA  
ini dens 2010 rang group sand_o11 ]` 3;8,  
ini dens 2010 rang group sand_o12 4j/8Otn  
;边界条件 c,e 0+  
fix x y z rang z -50.1 -49.9 7:;P>sF@  
fix x rang x  78.85  79.05 C#(4>'  
fix x rang x -78.55 -78.35 Dtelr=/s  
fix y rang y -.1  .1 V" I+E  
fix y rang y 0.9  1.1 p=eSHs{>A  
;应力条件 QarA.Ne~  
ini szz 0  grad 0 0 20e3   rang z -50 0     ;加权平均重度为1993 M,6m*  
ini sxx 0  grad 0 0 10e3   rang z -50 0 Ww@Rewo  
ini syy 0  grad 0 0 10e3   rang z -50 0 (/c9v8Pr(7  
;--------------------------------------------------------------------------------- IX-ir  
hist id 1 unbal JgEpqA12  
plot hist 1 VTD'D+t  
; qdzc"-gH`  
set grav 0 0 -10 cp[k[7XGD  
solve ratio 1.0e-4 E_-CsL%  
save station1.sav

土体强度、刚度参数给的太小了，你可以检查塑性区，不正常地大。 q S2#=  
我没试。

edelmann:土体强度、刚度参数给的太小了，你可以检查塑性区，不正常地大。  >Xh9{/o  
我没试。 (2013-03-15 17:38)  H 2\KI(  

3-srt^>w*  
土体强度指的哪个值？刚度参数是接触面上的刚度还是说体积模量和剪切模量？请赐教，因为我刚学习这个  还不是很懂

参数的问题