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

H: ;XU  
new $Yp.BE<}  
config dyn ;JMmr-@  
set dyn off ;Wa&Dg/5`  
;=============================================== u?rX:KkS  
;建模型 x6$3KDQm  
;=============================================== fdHFSnQ g  
def setup              ;定义函数setup 8F'm#0
 
   yc=1               ;y方向尺寸 |WpJen*?Y  
   czz1=-13.16 Jw>na _FJ  
   czz11=1.9            ;四分之一车站z方向尺寸   }:u-l3e  
   czz12=3            ;四分之一车站z方向尺寸 2kk; z0f  
   czz13=6.8            ;四分之一车站z方向尺寸 ?G<?:/CU  
   czz14=1.46            ;四分之一车站z方向尺寸     +md"X@k5*  
   ttx1=10.05              ;四分之一车站外围土体x方向尺寸 fW[RCd
 
   ttx2=-9.55              ;四分之一车站外围土体-x方向尺寸   (Q&O'ng1  
   ht1=-3.1 o\PHs4Ws'7  
   ht2=3.1 @6%7X7m  
   tt2=68.9 ;78.95             u.*}'C>^^v  
   tt1=-68.9 l_8ibLyo  
   xh=33.74              ；车站以下土体高度 4?+jvVq  
   F@#p  
end aL&9.L|1g  
setup   =(Y0wZP|  
;------------------------------------------------------------------------------------------------- IxG7eX!  
; 一、生成右侧土体 gzy|K%K  
;------------------------------------------------------------------------------------------------- W`P>vK@=  
def gen_Sando1 qq_ZkU@xg  
  command P!IXcPKW53  
gen zone brick      p0 0,0,ht1   p1 add ttx1,0,0  p2 add 0,yc,0  p3 add 0,0,ht2 & g_!xD;0  
                    p4 add ttx1,yc,0  p5 add 0,yc,ht2  p6 add ttx1,0,ht2  p7 add ttx1,yc,ht2 & ]g }5p4*&  
                    size 10 1 3                               ;ratio 1 1 1 1.2 ; xQhq*  
G{O{ p  
        end_command nywC]T  
end Kr<UPr  
gen_Sando1 ep0dT3&  
group sand_o1 range x 0 10.05  y 0 1  z -3.1 0 us8HXvvp{  
;------------------------------------------------------------------------------------------------- K`60[bdp  
def gen_Sando1 =6f)sZpPh  
  command ];5Auh0o  
gen zone brick      p0 ttx1,0,ht1   p1 add tt2,0,0  p2 add 0,yc,0  p3 add 0,0,ht2 & 6__HqBQ  
                    p4 add tt2,yc,0  p5 add 0,yc,ht2  p6 add tt2,0,ht2  p7 add tt2,yc,ht2 & | +uc;[`  
                    size 50 1 3 ;ratio 1.1 1 1 1 .[OS  
-3y
  
        end_command B9Wd '  
end V#+F*w?&D  
gen_Sando1 @,}tY ?>a  
group sand_o1 range x 10.05 78.95 y 0 1 z -3.1 0                   ~;YkR'q0_  
;------------------------------------------------------------------------------------------------- Ra6}<
o  
def gen_Sando2 Ed{sC[j=  
  command 9]lyV  
gen zone brick      p0 ttx1,0,-5   p1 add tt2,0,0  p2 add 0,yc,0  p3 add 0,0,czz11 & !qHB?]  
                    p4 add tt2,yc,0  p5 add 0,yc,czz11  p6 add tt2,0,czz11  p7 add tt2,yc,czz11 &   c.5u \I9"  
                    size 50 1 2 ;ratio 1.1 1 1 [Q\(kd*4  
\rO!lvX  
        end_command 3xmPY.  
end a,)/D_{1  
gen_Sando2 "zzb`T[8  
group sand_o2 range x 10.05 78.95 y 0 1 z -5 -3.1                                    ;11.26 13.16 i2;,\FI@t%  
_ b</ ::Tp  
Vg :''!4t2  
;------------------------------------------------------------------------------------------------- XX "3.zW  
;------------------------------------------------------------------------------------------------- .6I'V3:Kg  
def gen_Sando3 YXh!+}  
  command :h/v"2uDN  
gen zone brick      p0 ttx1,0,-8     p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz12 & f"NWv!  
                    p4 add tt2,yc,0  p5 add 0,yc,czz12  p6 add tt2,0,czz12  p7 add tt2,yc,czz12 &   ;j[:tt\k  
                    size 50 1 3 ;ratio 1.1 1 1 'H4?V  
5R%y3::$S  
        end_command B2KBJ4rI[1  
end .<xD'54  
gen_Sando3 &}|`h8JA]K  
group sand_o3 range x 10.05 78.95 y 0 1 z -8 -5                                          ;8.26 11.26 yq<W+b/  
;------------------------------------------------------------------------------------------------- @?;)x&<8?3  
def gen_Sando4 F\GNLi  
  command Y"^.6  
gen zone brick      p0 ttx1,0,-14.8   p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz13 & -N6ek`  
                    p4 add tt2,yc,0   p5 add 0,yc,czz13  p6 add tt2,0,czz13  p7 add tt2,yc,czz13 &   ZR"qrCSw`  
                    size 50 1 6 ;ratio 1.1 1 1 1_]%,  
^<uQ9p^B  
        end_command CJwzjH  
end V]"pM]>3X  
gen_Sando4 o*"Q{Xh#Qd  
group sand_o4 range x 10.05 78.95 y 0 1 z -14.8 -8                                   ;1.46 8.26 (d*||"  
ZG{#CC
=  
;------------------------------------------------------------------------------------------------- *~*"p)`<  
def gen_Sando5 D~|q^Ms,%  
  command |5&7;;$  
gen zone brick      p0 ttx1,0,-16.26   p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz14 & 5*Qzw[[=  
                    p4 add tt2,yc,0    p5 add 0,yc,czz14  p6 add tt2,0,czz14  p7 add tt2,yc,czz14 &   _^ic@h3'X~  
                    size 50 1 2 ;ratio 1.1 1 1 (<@`MPI\@  
rY&#g%B6Fp  
        end_command iel@"E 4  
end %o0H#7'  
gen_Sando5 9'(m"c_  
group sand_o5 range x 10.05 78.95 y 0 1 z -16.26 -14.8                            ;0 1.46 la4%Vqwgu  
;------------------------------------------------------------------------------------------------- feeHXKD|  
;------------------------------------------------------------------------------------------------- ue4{h  
1'iQlnMO@  
def gen_Sando6 #?e
MEws  
  command z]LVq k  
gen zone brick      p0 0,0,-50        p1 add ttx1,0,0    p2 add 0,yc,0     p3 add 0,0,xh & -v|lM8  
                    p4 add ttx1,yc,0  p5 add 0,yc,xh     p6 add ttx1,0,xh  p7 add ttx1,yc,xh &  
0I do_V  
                    size 10 1 20 ;ratio 1.2 1 1 Iz;^D!  
Ts|;5ya5m  
        end_command 
?^^TR/  
end [-81s!#mkw  
gen_Sando6 uq7/G|  
group sand_o6 range x 0 10.05  y 0 1 z  -50  -16.26 VVpJ +  
;------------------------------------------------------------------------------------------------- -j:yEZ4Oy  
def gen_Sando6 M'oZK  
  command GU9p'E  
gen zone brick      p0 ttx1,0,-50    p1 add tt2,0,0  p2 add 0,yc,0    p3 add 0,0,xh & NS TO\36  
                    p4 add tt2,yc,0  p5 add 0,yc,xh  p6 add tt2,0,xh  p7 add tt2,yc,xh &   cS2]?zI  
                    size 40 1 20 ;ratio 1.1 1   }_mMQg2>=  
LyR<cd$W  
        end_command o>T+fBHE  
end SQ]M"&\{y  
gen_Sando6 Fu\#:+5\  
group sand_o6 range  x 10.05 78.95 y 0 1  z -50   -16.26                      ;-33.74 0 i70\`6*;B  
;------------------------------------------------------------------------------------------------- j~>J?w9<O  
;------------------------------------------------------------------------------------------------- L/%{,7l<^?  
; 二、生成左侧土体 R6:m@  
;------------------------------------------------------------------------------------------------- -^;,m=4{3  
;------------------------------------------------------------------------------------------------- h?}S|>9  
def gen_Sand11 Uz[#ye  
  command T&bB8tQk  
gen zone brick      p0 0,0,ht1        p1 add 0,yc,0  p2 add ttx2,0,0      p3 add 0,0,ht2 & )q\6pO@  
                    p4 add ttx2,yc,0  p5 add ttx2,0,ht2  p6 add 0,yc,ht2  p7 add ttx2,yc,ht2 & h|ib*%P_  
                    size 1 10 3                               ;ratio 1 1 1 1.2 KoWG:~>|  
1jAuW~  
        end_command J4z&J SY  
end AAKc8{  
gen_Sand11 t]iKU@3  
group sand_o7 range x -9.55 0 y 0 1 z -3.1 0 )vg@Kc26  
;------------------------------------------------------------------------------------------------- Fpwhy
ls  
def gen_Sand2 PlT_]p  
  command rY1jC\  
gen zone brick      p0 ttx2,0,ht1    p1 add 0,yc,0   p2 add tt1,0,0     p3 add 0,0,ht2 & '<*%<J{(  
                    p4 add tt1,yc,0  p5 add tt1,0,ht2  p6 add 0,yc,ht2  p7 add tt1,yc,ht2 & hJw]hVYa  
                    size 1 50  3 ;ratio 1 1.1 1 _JA)""l%  
&OEBAtc/  
        end_command +_gA"I  
end ;B(16&l=q  
gen_Sand2 gS`Z>+V5!c  
group sand_o7 range x -78.45 -9.55 y 0 1 z -3.1 0                   0]$-}A
YM  
;------------------------------------------------------------------------------------------------- >IA1 \?(  
def gen_Sand31 0>e]i[P.  
  command @+)T"5_Y[  
gen zone brick      p0 ttx2,0,-5     p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz11 & -mo4`F  
                    p4 add tt1,yc,0  p5 add tt1,0,czz11  p6 add 0,yc,czz11  p7 add tt1,yc,czz11 &   ]1|7V|N6  
                    size 1 50 2 ;ratio 1 1.1 1 -7o-d-d F  
yS[HYq  
        end_command ac966<
#  
end IjXxH]2  
gen_Sand31 E
x@}x#3  
group sand_o8 range x -78.45 -9.55 y 0 1 z -5 -3.1                                    ;11.26 13.16 ,_D@ggL-  
;------------------------------------------------------------------------------------------------- qK~]au:C  
def gen_Sand32 )7Qp9Fxo  
  command |z&7KoYK'  
gen zone brick      p0 ttx2,0,-8     p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz12 & (v:8p!QN  
                    p4 add tt1,yc,0  p5 add tt1,0,czz12  p6 add 0,yc,czz12  p7 add tt1,yc,czz12 &   hx/A215L  
                    size 1 50 3 ;ratio 1 1.1 1 C7}iwklcsa  
b^(
)[4M;  
        end_command PI,2b(`h_  
end PL!dkaD^y>  
gen_Sand32 Ml{4)%~Y7f  
group sand_o9 range x -78.45 -9.55 y 0 1  z -8 -5                                          ;8.26 11.26 BQeg-M  
;------------------------------------------------------------------------------------------------- )KkV<$  
def gen_Sand33 T!pZj_ h=  
  command LfK/wSvWw  
gen zone brick      p0 ttx2,0,-14.8   p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz13 & #J8(*!I  
                    p4 add tt1,yc,0   p5 add tt1,0,czz13  p6 add 0,yc,czz13  p7 add tt1,yc,czz13 &   "?W8o[c+  
                    size 1 50 6 ;ratio 1 1.1 1 N=~DSsw  
!x||ObW\H  
        end_command ">S1,rhgS  
end )nK+`{;@!  
gen_Sand33 w\V<6_[vv.  
group sand_o10 range x -78.45 -9.55 y 0 1  z -14.8 -8                                   ;1.46 8.26 i }Zz[b  
7s2*VKr  
;------------------------------------------------------------------------------------------------- r(_Fr#Qn  
def gen_Sand34 0tPw
hJ  
  command * kUb[  
gen zone brick      p0 ttx2,0,-16.26   p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz14 & |lu@rN  
                    p4 add tt1,yc,0    p5 add tt1,0,czz14  p6 add 0,yc,czz14  p7 add tt1,yc,czz14 &   a5d_= :S;  
                    size 1 50 2 ;ratio 1 1.1 1 =}u?1~V  
TV0Y{x*~iH  
        end_command m5HMtoU  
end PGVp1TQ  
gen_Sand34 kGakdLl  
group sand_o11 range x -78.45 -9.55 y 0 1  z -16.26 -14.8                            ;0 1.46 oR7f3';?6  
;------------------------------------------------------------------------------------------------- 8493O x4 O  
;------------------------------------------------------------------------------------------------- /q7$"wP  
*"6A>:rQs  
def gen_Sand4 9}42s+  
  command %.Kr`#lCr  
gen zone brick      p0 0,0,-50        p1 add 0,yc,0  p2 add ttx2,0,0     p3 add 0,0,xh & J~ +p7S  
                    p4 add ttx2,yc,0  p5 add ttx2,0,xh   p6 add 0,yc,xh  p7 add ttx2,yc,xh &   3/(eK%d4Xb  
                    size 1 10 20 jIHY[yDT  
&_j<!3*  
        end_command jZvIqR/  
end nLPd]%78>  
gen_Sand4 se}$/Y}t  
group sand_o12 range x -9.55 0 y 0 1 z -16.26 -50 322-'S3<  
;------------------------------------------------------------------------------------------------- QA)"3g   
def gen_Sand5 w vI v+Q9  
  command nrXKS&6  
gen zone brick      p0 ttx2,0,-50    p1 add 0,yc,0  p2 add tt1,0,0    p3 add 0,0,xh & XaoVv2=G~  
                    p4 add tt1,yc,0  p5 add tt1,0,xh  p6 add 0,yc,xh  p7 add tt1,yc,xh &   "GJ.`Hj  
                    size 1 50 20 ;ratio 1 1.1 1 8,VEuBZ  
0m\( @2E  
        end_command =)N6R  
end HzuG- V  
gen_Sand5 m6 Y0,9
 
group sand_o12 range x -78.45 -9.55 y 0 1 z -16.26 -50                        ;-33.74 0 m`Z.xIA7;  
;------------------------------------------------------------------------------------------------- ^EW6}oj[  
call Flacmodel_haitang.dat ycvgF6Me<  
;------------------------------------------------------------------------------------------------- NqFfz9G)  
plot add axes +JXn   
;-------------------------------------------------------------------------------------- :Dtm+EQ  
;二、生成接触面的命令 J9y}rGO  
;-------------------------------------------------------------------------------------- &NbSG
+t  
interface 1 face range x -9.55 10.05 y 0 1 z -16.27 -16.25 Zw%:mZN   
interface 2 face range x -9.55 10.05 y 0 1 z -3.11 -3.09 CP$,fj  
interface 3 face range x -9.56 -9.54 y 0 1 z -5 -3.1  p1&=D%/  
interface 4 face range x -9.56 -9.54 y 0 1 z -8  -5 >fkV65w{*  
interface 5 face range x -9.56 -9.54 y 0 1 z -14.8 -8 HW'I$ .  
interface 6 face range x -9.56 -9.54 y 0 1 z -16.26 -14.8 %zDi|WZ  
'dv
(  
int  1  prop kn  1.570e9  ks  1.570e9  fric  10 coh 10000                   ;prop属性 kn法向刚度 ks剪切刚度 fric为摩擦角 coh内聚力 Mpk7$=hjc  
int  2  prop kn  1.570e9  ks  1.570e9  fric  15 coh 55000 D`.\c#;cN  
int  3  prop kn  1.570e9  ks  1.570e9  fric  26 coh 8000 02(Ob  
int  4  prop kn  1.570e9  ks  1.570e9  fric  15 coh 55000 qw)O
u]L=  
int  5  prop kn  1.570e9  ks  1.570e9  fric  14 coh 24000 c|(Q[=   
int  6  prop kn  1.570e9  ks  1.570e9  fric  17 coh 56000 U7,.L  
ini y add -2.0 rang group 1 $YJi]:3&  
;=============================================== `bn@;7`X  
;初始应力计算程序 wsc=6/#u  
;=============================================== -*-"kzgd  
;定义模型及性质 iU=:YPE+.  
model  mohr    rang  group sand_o1                                           Ys?0hd<cn  
prop bulk 1.875e7   shear   0.967e7 coh 10000   fric 10 rang group sand_o1 u09D`QPP]  
model  mohr    rang  group sand_o2                                           i1]}Q$  
prop bulk 1.875e7   shear   0.967e7 coh 55000   fric 15 rang group sand_o2 ,W/Y@Sc
C  
model  mohr    rang  group sand_o3                                           62G%.'7  
prop bulk 2.068e7   shear   0.793e7 coh 8000   fric 26 rang group sand_o3 z U*Mk  
model  mohr    rang  group sand_o4                                           AI,E9  
prop bulk 1.875e7   shear   0.967e7 coh 55000   fric 15 rang group sand_o4 73{<;z}i  
model  mohr    rang  group sand_o5                                           d>r]xXB6  
prop bulk 1.396e7   shear   0.758e7 coh 24000   fric 14 rang group sand_o5 AP4s_X+=  
model  mohr    rang  group sand_o6                                           J*ZcZ FbWN  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o6 :`<MlX  
model  mohr    rang  group sand_o7                                           *TCV}=V G  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o7 2Qc_TgWF  
model  mohr    rang  group sand_o8                                           ="Azg8W  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o8 Fta=yH}  
model  mohr    rang  group sand_o9                                           <A`S
C;k\u  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o9 o>m*e7l,  
model  mohr    rang  group sand_o10                                           t+^__~IX  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o10 |V
u`-L'Jz  
model  mohr    rang  group sand_o11                                           @ Yo*h"s  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o11     O
RXH<;^0y  
model  mohr    rang  group sand_o12                                           o4L
VG  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o12 r/0AM}[!*j  
pE0@m-p  
mo elas                 rang group 1 qNMYZ0,  
;prop bulk 20.8e9  shear 15.58e9  rang group 1   ;正确 E>2AG3)  
prop bulk 22.8e9  shear 17.08e9  rang group 1   ;只把柱子等效 }/IP\1bG  
ini dens 2500 rang group 1 ?#nk}=;g8  
ini dens 1900 rang group sand_o1     ;18000 (hRg0Z=  
ini dens 1990 rang group sand_o2 ~*~aFf5  
ini dens 1910 rang group sand_o3 rg]z  
ini dens 1990 rang group sand_o4 j<9^BNl  
ini dens 1970 rang group sand_o5 !.4q{YWcYk  
ini dens 2010 rang group sand_o6 9'|_1Q.b^  
ini dens 2010 rang group sand_o7 oL!C(\ER
h  
ini dens 2010 rang group sand_o8 J%!vhQ  
ini dens 2010 rang group sand_o9 4Yt'I#*  
ini dens 2010 rang group sand_o10 t~qAA\p}o  
ini dens 2010 rang group sand_o11 E5*pD*#  
ini dens 2010 rang group sand_o12 89P7iSV#*  
;边界条件 P86wRq  
fix x y z rang z -50.1 -49.9 0U#m7j  
fix x rang x  78.85  79.05 vAOThj)  
fix x rang x -78.55 -78.35 fl4'dv  
fix y rang y -.1  .1 ygK,t*T20  
fix y rang y 0.9  1.1 R4zOiBi'B  
;应力条件 W&3,XFnI_  
ini szz 0  grad 0 0 20e3   rang z -50 0     ;加权平均重度为1993 55,2eg#{O  
ini sxx 0  grad 0 0 10e3   rang z -50 0 1:u~T@;" `  
ini syy 0  grad 0 0 10e3   rang z -50 0 %/!f^PIwX  
;--------------------------------------------------------------------------------- 60D6UW  
hist id 1 unbal !RjC0,  
plot hist 1 &b-&0rTqz  
; E%Ko[G  
set grav 0 0 -10 !2/o]_K@+  
solve ratio 1.0e-4 fj9&J[
 
save station1.sav

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

edelmann:土体强度、刚度参数给的太小了，你可以检查塑性区，不正常地大。 F9r|EU#;  
我没试。 (2013-03-15 17:38)  v9#F\F/  

N~?(<DyZR  
土体强度指的哪个值？刚度参数是接触面上的刚度还是说体积模量和剪切模量？请赐教，因为我刚学习这个  还不是很懂

参数的问题