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

V,<3uQD9a  
new D.zEE-cGyb  
config dyn 5w)tsGX\  
set dyn off tCAh?nR  
;=============================================== n0_q-
8r  
;建模型 Ah,X?0+  
;=============================================== $<.\,wW*'w  
def setup              ;定义函数setup \kU0D  
   yc=1               ;y方向尺寸 6]yYiz2Xn  
   czz1=-13.16 n=#AH;42  
   czz11=1.9            ;四分之一车站z方向尺寸   =T'N6x5@  
   czz12=3            ;四分之一车站z方向尺寸 XBE+O7  
   czz13=6.8            ;四分之一车站z方向尺寸 yrOWC  
   czz14=1.46            ;四分之一车站z方向尺寸     `LU[+F8<  
   ttx1=10.05              ;四分之一车站外围土体x方向尺寸 !63
p?Q=  
   ttx2=-9.55              ;四分之一车站外围土体-x方向尺寸   ES[H^}|Gi  
   ht1=-3.1 p]/HZS.-b  
   ht2=3.1 +G';no\h  
   tt2=68.9 ;78.95             @:\Iw"P  
   tt1=-68.9 F.2<G.9  
   xh=33.74              ；车站以下土体高度 j+He8w-4  
   p R=FH#  
end fvn`$  
setup   wQ /IT}-  
;------------------------------------------------------------------------------------------------- qjQR0MC  
; 一、生成右侧土体 P(BjXMd  
;------------------------------------------------------------------------------------------------- v#1}( hb  
def gen_Sando1 jnfktDV'  
  command u`I&&  
gen zone brick      p0 0,0,ht1   p1 add ttx1,0,0  p2 add 0,yc,0  p3 add 0,0,ht2 & @>n7  
                    p4 add ttx1,yc,0  p5 add 0,yc,ht2  p6 add ttx1,0,ht2  p7 add ttx1,yc,ht2 & b(JQ>,hX  
                    size 10 1 3                               ;ratio 1 1 1 1.2 `)Z"||8K  
EkotVzR5  
        end_command [IT*>;b+?  
end fYrC;&n  
gen_Sando1 c5x2FM z  
group sand_o1 range x 0 10.05  y 0 1  z -3.1 0 F!DDlYUz.  
;------------------------------------------------------------------------------------------------- p4QQ5O$;  
def gen_Sando1 0$)uOUVJ  
  command Vmq:As^a  
gen zone brick      p0 ttx1,0,ht1   p1 add tt2,0,0  p2 add 0,yc,0  p3 add 0,0,ht2 & "mlVs/nsyG  
                    p4 add tt2,yc,0  p5 add 0,yc,ht2  p6 add tt2,0,ht2  p7 add tt2,yc,ht2 & bBwQ1,c$  
                    size 50 1 3 ;ratio 1.1 1 1 hTZ&  
f]Jn\7j4  
        end_command h3-y}.VjG  
end M0Y#=u.  
gen_Sando1 qxKW%{6o  
group sand_o1 range x 10.05 78.95 y 0 1 z -3.1 0                   28,Hd!{  
;------------------------------------------------------------------------------------------------- t5: 1' N9P  
def gen_Sando2 G?`{OW3:_  
  command xNlxi  
gen zone brick      p0 ttx1,0,-5   p1 add tt2,0,0  p2 add 0,yc,0  p3 add 0,0,czz11 & oasp/Y.p  
                    p4 add tt2,yc,0  p5 add 0,yc,czz11  p6 add tt2,0,czz11  p7 add tt2,yc,czz11 &   DFDlp  
                    size 50 1 2 ;ratio 1.1 1 1 FXMrD,qVg  
Vwxb6,}Z  
        end_command e ]-fb{oVH  
end c>I(6$  
gen_Sando2 6Ih8~Hu  
group sand_o2 range x 10.05 78.95 y 0 1 z -5 -3.1                                    ;11.26 13.16 D d['e  
*AK{GfP_  
AA=zDB<N  
;------------------------------------------------------------------------------------------------- []^fb,5a  
;------------------------------------------------------------------------------------------------- r<< ]41  
def gen_Sando3 sz;B-1^6  
  command {A<pb{<u  
gen zone brick      p0 ttx1,0,-8     p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz12 & zt&"K0X|  
                    p4 add tt2,yc,0  p5 add 0,yc,czz12  p6 add tt2,0,czz12  p7 add tt2,yc,czz12 &   'W'['TV  
                    size 50 1 3 ;ratio 1.1 1 1  zciL'9  
+H"[WZ5  
        end_command dsJMhB_41U  
end -|F(qf  
gen_Sando3 >\f'QQ  
group sand_o3 range x 10.05 78.95 y 0 1 z -8 -5                                          ;8.26 11.26 B^|^hZZ>  
;------------------------------------------------------------------------------------------------- i'bviD  
def gen_Sando4 ;iA6[uz  
  command ]#$l"ss,  
gen zone brick      p0 ttx1,0,-14.8   p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz13 & fkKk/M>1  
                    p4 add tt2,yc,0   p5 add 0,yc,czz13  p6 add tt2,0,czz13  p7 add tt2,yc,czz13 &   73C7g< Mx  
                    size 50 1 6 ;ratio 1.1 1 1 }4$k-,1S  
K=Q<G:+&V  
        end_command |
q+3X)Y  
end 3m;*gOLk6  
gen_Sando4 dWKjVf  
group sand_o4 range x 10.05 78.95 y 0 1 z -14.8 -8                                   ;1.46 8.26 HN
XMM  
%?2:1o  
;------------------------------------------------------------------------------------------------- Dx1f<A1  
def gen_Sando5 2!@ER i  
  command 0c{-$K}  
gen zone brick      p0 ttx1,0,-16.26   p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz14 & 6KH&-ffd  
                    p4 add tt2,yc,0    p5 add 0,yc,czz14  p6 add tt2,0,czz14  p7 add tt2,yc,czz14 &   W20qn>{z  
                    size 50 1 2 ;ratio 1.1 1 1 O@9<7@h+Nl  
_8QHx;}  
        end_command "W71#n+[  
end le7!:4/8  
gen_Sando5 M\m6|P  
group sand_o5 range x 10.05 78.95 y 0 1 z -16.26 -14.8                            ;0 1.46 $3yzB9\a"  
;------------------------------------------------------------------------------------------------- d ,"L8  
;------------------------------------------------------------------------------------------------- \d:AV(u  
5~ip N/)E  
def gen_Sando6 VbU*&{j  
  command 3'*}ZDC  
gen zone brick      p0 0,0,-50        p1 add ttx1,0,0    p2 add 0,yc,0     p3 add 0,0,xh & ,ma4bqRMc  
                    p4 add ttx1,yc,0  p5 add 0,yc,xh     p6 add ttx1,0,xh  p7 add ttx1,yc,xh &   QV%,s!_b  
                    size 10 1 20 ;ratio 1.2 1 1 vV$6fvS  
?z{Z!Bt?=)  
        end_command ta  
end S+*%u/;l  
gen_Sando6 U
ngK9uB~  
group sand_o6 range x 0 10.05  y 0 1 z  -50  -16.26 i3dV2^O  
;------------------------------------------------------------------------------------------------- :qAF}|6  
def gen_Sando6 )@\Eibt2oH  
  command |'+ [ '  
gen zone brick      p0 ttx1,0,-50    p1 add tt2,0,0  p2 add 0,yc,0    p3 add 0,0,xh & sjW;Nsp  
                    p4 add tt2,yc,0  p5 add 0,yc,xh  p6 add tt2,0,xh  p7 add tt2,yc,xh &   jhx
@6[  
                    size 40 1 20 ;ratio 1.1 1   W{!Slf  
bnZ`Wc*5b  
        end_command -|cB7P  
end <&CzM"\Em  
gen_Sando6  26[.te9  
group sand_o6 range  x 10.05 78.95 y 0 1  z -50   -16.26                      ;-33.74 0 =@\Li)Y  
;------------------------------------------------------------------------------------------------- 6'a1]K  
;------------------------------------------------------------------------------------------------- B[IqLD'6  
; 二、生成左侧土体 L"x9O'U  
;------------------------------------------------------------------------------------------------- yN/
Uyhq  
;------------------------------------------------------------------------------------------------- 0}4FwcCr\  
def gen_Sand11 :c/](M  
  command Ld4Jp`Zg  
gen zone brick      p0 0,0,ht1        p1 add 0,yc,0  p2 add ttx2,0,0      p3 add 0,0,ht2 & :k&5Z`>)  
                    p4 add ttx2,yc,0  p5 add ttx2,0,ht2  p6 add 0,yc,ht2  p7 add ttx2,yc,ht2 & k62KZ5| D  
                    size 1 10 3                               ;ratio 1 1 1 1.2 w~|z0;hC  
#J w\pOn  
        end_command a? <Ar#)j  
end XEpwk,8*g  
gen_Sand11 <&m `)FJ  
group sand_o7 range x -9.55 0 y 0 1 z -3.1 0 y%Wbm&h  
;------------------------------------------------------------------------------------------------- K/A1g.$  
def gen_Sand2 *Q`y'6S  
  command wK*b2r}0/  
gen zone brick      p0 ttx2,0,ht1    p1 add 0,yc,0   p2 add tt1,0,0     p3 add 0,0,ht2 & /4{.J=R}  
                    p4 add tt1,yc,0  p5 add tt1,0,ht2  p6 add 0,yc,ht2  p7 add tt1,yc,ht2 & WoSJp5By$  
                    size 1 50  3 ;ratio 1 1.1 1 VGL!)1
b  
dk@j!-q^  
        end_command 6$'0^Ftm'  
end ,/1[(^e  
gen_Sand2 41s[p56+@  
group sand_o7 range x -78.45 -9.55 y 0 1 z -3.1 0                   sqjv3=
}  
;------------------------------------------------------------------------------------------------- 'kE^oX_  
def gen_Sand31 !fT3mI6u\  
  command #- z(]Y,y  
gen zone brick      p0 ttx2,0,-5     p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz11 &  Z5[f  
                    p4 add tt1,yc,0  p5 add tt1,0,czz11  p6 add 0,yc,czz11  p7 add tt1,yc,czz11 &   wf&1,t3Bgn  
                    size 1 50 2 ;ratio 1 1.1 1  gU%R9  
-MTO=#5z  
        end_command <pp<%~_Z  
end wPRs.(]_  
gen_Sand31 R1/87eB
 
group sand_o8 range x -78.45 -9.55 y 0 1 z -5 -3.1                                    ;11.26 13.16 $?Yw{%W  
;------------------------------------------------------------------------------------------------- a"pejW`m  
def gen_Sand32 M].D27  
  command ?YA5g' l  
gen zone brick      p0 ttx2,0,-8     p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz12 & 67b[T~92o  
                    p4 add tt1,yc,0  p5 add tt1,0,czz12  p6 add 0,yc,czz12  p7 add tt1,yc,czz12 &   SHt#%3EU  
                    size 1 50 3 ;ratio 1 1.1 1 f<K7m  
5;yVA  
        end_command n%%u0a%  
end c ;@k\6  
gen_Sand32 L\/YS;Y  
group sand_o9 range x -78.45 -9.55 y 0 1  z -8 -5                                          ;8.26 11.26 wJb"X=i*  
;------------------------------------------------------------------------------------------------- $ 8WJ$73  
def gen_Sand33 (Gp|K6
  
  command $/XR
/  
gen zone brick      p0 ttx2,0,-14.8   p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz13 & Xw<5VIAHm;  
                    p4 add tt1,yc,0   p5 add tt1,0,czz13  p6 add 0,yc,czz13  p7 add tt1,yc,czz13 &   +`$[h2Z=:  
                    size 1 50 6 ;ratio 1 1.1 1 "~=\AB=+Z  
i#1~<U  
        end_command }rbsarG@  
end B_uAa5'  
gen_Sand33 ]Yd7  
group sand_o10 range x -78.45 -9.55 y 0 1  z -14.8 -8                                   ;1.46 8.26 x}*Y =Xh  
z;KUIWg  
;------------------------------------------------------------------------------------------------- j[q$;uSD  
def gen_Sand34 /`YbHYNF[  
  command Z+J;nl  
gen zone brick      p0 ttx2,0,-16.26   p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz14 & .|Bmg6g*  
                    p4 add tt1,yc,0    p5 add tt1,0,czz14  p6 add 0,yc,czz14  p7 add tt1,yc,czz14 &   VjhwafYC  
                    size 1 50 2 ;ratio 1 1.1 1 |&xjuBC  
{I~[a#^  
        end_command 8~s-@3J  
end J`@#yHL  
gen_Sand34 cc|"^-j-7  
group sand_o11 range x -78.45 -9.55 y 0 1  z -16.26 -14.8                            ;0 1.46 \y*,N^wu  
;------------------------------------------------------------------------------------------------- YTo^Q&  
;------------------------------------------------------------------------------------------------- b/Q\ .!  
D
|BP]j}6  
def gen_Sand4 Qx%]u8s  
  command 4t;m^I
v  
gen zone brick      p0 0,0,-50        p1 add 0,yc,0  p2 add ttx2,0,0     p3 add 0,0,xh & dtT2h>h9  
                    p4 add ttx2,yc,0  p5 add ttx2,0,xh   p6 add 0,yc,xh  p7 add ttx2,yc,xh &   my(yN|  
                    size 1 10 20 h1uD>heGl  
{Vj25Gt  
        end_command [AHZOA  
end zMI0W&P M  
gen_Sand4 }^/;8cfLY  
group sand_o12 range x -9.55 0 y 0 1 z -16.26 -50 H%:u9DlEK/  
;------------------------------------------------------------------------------------------------- |}y6U< I  
def gen_Sand5 }bxx]rDl  
  command $}/Q%
r  
gen zone brick      p0 ttx2,0,-50    p1 add 0,yc,0  p2 add tt1,0,0    p3 add 0,0,xh & :.,3Zw{l  
                    p4 add tt1,yc,0  p5 add tt1,0,xh  p6 add 0,yc,xh  p7 add tt1,yc,xh &   NA`E
G,2  
                    size 1 50 20 ;ratio 1 1.1 1 68v59)0U  
_-.~>C  
        end_command *;I F^u1  
end iTq~^9G  
gen_Sand5 .
Iret:  
group sand_o12 range x -78.45 -9.55 y 0 1 z -16.26 -50                        ;-33.74 0 mlmXFEC  
;------------------------------------------------------------------------------------------------- F< #!83*%  
call Flacmodel_haitang.dat D;l)&"|r?  
;------------------------------------------------------------------------------------------------- l8_TeO  
plot add axes ?fXlrJ  
;-------------------------------------------------------------------------------------- @`X-=GCl  
;二、生成接触面的命令 $CYB&|d  
;-------------------------------------------------------------------------------------- Nu%JI6&R  
interface 1 face range x -9.55 10.05 y 0 1 z -16.27 -16.25 /`Wd+  
interface 2 face range x -9.55 10.05 y 0 1 z -3.11 -3.09 &3;yho8v@  
interface 3 face range x -9.56 -9.54 y 0 1 z -5 -3.1 c*"TmDY  
interface 4 face range x -9.56 -9.54 y 0 1 z -8  -5 b@&ydgmaQ  
interface 5 face range x -9.56 -9.54 y 0 1 z -14.8 -8 :8<\]}J  
interface 6 face range x -9.56 -9.54 y 0 1 z -16.26 -14.8 D|Tv`47ntu  
fDa$TbhjI  
int  1  prop kn  1.570e9  ks  1.570e9  fric  10 coh 10000                   ;prop属性 kn法向刚度 ks剪切刚度 fric为摩擦角 coh内聚力 Cq0S8Or0  
int  2  prop kn  1.570e9  ks  1.570e9  fric  15 coh 55000 +yYv"J  
int  3  prop kn  1.570e9  ks  1.570e9  fric  26 coh 8000 # Y*cLN`Y7  
int  4  prop kn  1.570e9  ks  1.570e9  fric  15 coh 55000 bof{R{3q  
int  5  prop kn  1.570e9  ks  1.570e9  fric  14 coh 24000 G)E#wh_S^  
int  6  prop kn  1.570e9  ks  1.570e9  fric  17 coh 56000 Cl+TjmOV\`  
ini y add -2.0 rang group 1 = +uUWJ&1G  
;=============================================== 2Rp'ju~O)/  
;初始应力计算程序 pSq3\#Twr  
;=============================================== r|WoM39bp  
;定义模型及性质 V %cU@  
model  mohr    rang  group sand_o1                                           5^\f[}  
prop bulk 1.875e7   shear   0.967e7 coh 10000   fric 10 rang group sand_o1 `[+nz rLkO  
model  mohr    rang  group sand_o2                                           :_qgpE<  
prop bulk 1.875e7   shear   0.967e7 coh 55000   fric 15 rang group sand_o2 Hkv4t5F  
model  mohr    rang  group sand_o3                                           FB0y  
prop bulk 2.068e7   shear   0.793e7 coh 8000   fric 26 rang group sand_o3 ]bq<vI%  
model  mohr    rang  group sand_o4                                           5[Ryc[  
prop bulk 1.875e7   shear   0.967e7 coh 55000   fric 15 rang group sand_o4 ~!,Q<?  
model  mohr    rang  group sand_o5                                           |x AwiF_  
prop bulk 1.396e7   shear   0.758e7 coh 24000   fric 14 rang group sand_o5 ly d[GfJ  
model  mohr    rang  group sand_o6                                           Q1Jkt  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o6 4^KoHeM6  
model  mohr    rang  group sand_o7                                           V+()`>44  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o7 1MV\ ^l_  
model  mohr    rang  group sand_o8                                           $UZ4,S?V  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o8 G
e|& H]W  
model  mohr    rang  group sand_o9                                           hr<E%J1k%  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o9 KJwkkCE/=  
model  mohr    rang  group sand_o10                                           &B^vHH  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o10 ^;2dZgJ4^  
model  mohr    rang  group sand_o11                                           lxRzyx  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o11     I 
[J0r  
model  mohr    rang  group sand_o12                                           *mM+(]8US  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o12 m4@y58n=  
*XmOWV2Y_  
mo elas                 rang group 1 !Cj1:P  
;prop bulk 20.8e9  shear 15.58e9  rang group 1   ;正确 "mAVkq~  
prop bulk 22.8e9  shear 17.08e9  rang group 1   ;只把柱子等效 UC^Bn1  
ini dens 2500 rang group 1 #lax0IYY=  
ini dens 1900 rang group sand_o1     ;18000 6/9h=-w&  
ini dens 1990 rang group sand_o2 bl>MD8bzLE  
ini dens 1910 rang group sand_o3 ,~Y5vnaOQ  
ini dens 1990 rang group sand_o4 >NN|vj  
ini dens 1970 rang group sand_o5 N b(f  
ini dens 2010 rang group sand_o6 v UAYYe  
ini dens 2010 rang group sand_o7 $a6&OH/  
ini dens 2010 rang group sand_o8 C61KY7iyR  
ini dens 2010 rang group sand_o9 !}()mrIlP  
ini dens 2010 rang group sand_o10 'nNw  
ini dens 2010 rang group sand_o11 5hE8b {V  
ini dens 2010 rang group sand_o12 U/M(4H3>H  
;边界条件 =L$};ko  
fix x y z rang z -50.1 -49.9 3:a}<^DuCS  
fix x rang x  78.85  79.05 FeS6>/  
fix x rang x -78.55 -78.35 i&^?p|eKa  
fix y rang y -.1  .1 VwoCRq*  
fix y rang y 0.9  1.1 W|T"'M_  
;应力条件 AAsl)  
ini szz 0  grad 0 0 20e3   rang z -50 0     ;加权平均重度为1993 :&dY1.<N+  
ini sxx 0  grad 0 0 10e3   rang z -50 0 /|y3M/;F  
ini syy 0  grad 0 0 10e3   rang z -50 0 2I:vie  
;--------------------------------------------------------------------------------- ^sq3@*hCw  
hist id 1 unbal Kbc-$oneR  
plot hist 1 Y.yM1 z  
; Dd*T5A?  
set grav 0 0 -10 fp9ksxb@m  
solve ratio 1.0e-4 )\uy 0+b  
save station1.sav

土体强度、刚度参数给的太小了，你可以检查塑性区，不正常地大。 ls,gQ]B:P  
我没试。

edelmann:土体强度、刚度参数给的太小了，你可以检查塑性区，不正常地大。 gMI%z2]'-  
我没试。 (2013-03-15 17:38)  ,{8~TVO  

{%S1x{U}W-  
土体强度指的哪个值？刚度参数是接触面上的刚度还是说体积模量和剪切模量？请赐教，因为我刚学习这个  还不是很懂

参数的问题