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

ZQrgYeQl"  
new )hK1W\5  
config dyn xqaw00,s  
set dyn off vX}mwK8  
;=============================================== hin6cac  
;建模型 }i2dXC/  
;=============================================== AmT*{Fz8  
def setup              ;定义函数setup 7=]Y7"XCf  
   yc=1               ;y方向尺寸 |c!lZo/  
   czz1=-13.16 +@K8:}lOW
 
   czz11=1.9            ;四分之一车站z方向尺寸   7.xJ:r|  
   czz12=3            ;四分之一车站z方向尺寸 ,5Tw
5<S  
   czz13=6.8            ;四分之一车站z方向尺寸 `H\NJ,  
   czz14=1.46            ;四分之一车站z方向尺寸     $a+)v#?,  
   ttx1=10.05              ;四分之一车站外围土体x方向尺寸 \fD[Ej  
   ttx2=-9.55              ;四分之一车站外围土体-x方向尺寸   nB86oQ/S  
   ht1=-3.1 r#K"d  
   ht2=3.1 1V1T1  
   tt2=68.9 ;78.95             y1hJVYE2  
   tt1=-68.9 F84?Mi{r2  
   xh=33.74              ；车站以下土体高度 .(zZTyZr  
   0qD.OF)8  
end v7- d+P=  
setup   ^->vUf7PX  
;------------------------------------------------------------------------------------------------- hdee]qLS  
; 一、生成右侧土体 !<
MW*7P=  
;------------------------------------------------------------------------------------------------- vghn+P8  
def gen_Sando1 }@wXm  
  command w^QqYUL${  
gen zone brick      p0 0,0,ht1   p1 add ttx1,0,0  p2 add 0,yc,0  p3 add 0,0,ht2 & DR#[\RzNI  
                    p4 add ttx1,yc,0  p5 add 0,yc,ht2  p6 add ttx1,0,ht2  p7 add ttx1,yc,ht2 & .ZOyZnr Z  
                    size 10 1 3                               ;ratio 1 1 1 1.2 gc{5/U9H*  
6c&OR2HGqO  
        end_command DX#F]8bWl  
end Qmn'G4#@E  
gen_Sando1 %q,^A+
=  
group sand_o1 range x 0 10.05  y 0 1  z -3.1 0 CI,xp  
;------------------------------------------------------------------------------------------------- \{`*`WQF  
def gen_Sando1 Q*AgFF%wn  
  command >eaK@
u-'0  
gen zone brick      p0 ttx1,0,ht1   p1 add tt2,0,0  p2 add 0,yc,0  p3 add 0,0,ht2 & T 9?!.o  
                    p4 add tt2,yc,0  p5 add 0,yc,ht2  p6 add tt2,0,ht2  p7 add tt2,yc,ht2 & JZrUl^8E  
                    size 50 1 3 ;ratio 1.1 1 1 =;HmU.Uek%  
=;A~$[g  
        end_command +v'n[xa1v  
end ~b{j`T  
gen_Sando1 78<QNlKn  
group sand_o1 range x 10.05 78.95 y 0 1 z -3.1 0                   u+uu?.bM  
;------------------------------------------------------------------------------------------------- 9TW  
def gen_Sando2 {1`n^j(>  
  command TVFxEV7Fx  
gen zone brick      p0 ttx1,0,-5   p1 add tt2,0,0  p2 add 0,yc,0  p3 add 0,0,czz11 & .[#bOp*  
                    p4 add tt2,yc,0  p5 add 0,yc,czz11  p6 add tt2,0,czz11  p7 add tt2,yc,czz11 &   <(-4?"1  
                    size 50 1 2 ;ratio 1.1 1 1 'k[gxk|d2  
9 !qVYU42(  
        end_command G6x
2!Ny  
end Q Ph6 p3bg  
gen_Sando2 sOW,hpNW  
group sand_o2 range x 10.05 78.95 y 0 1 z -5 -3.1                                    ;11.26 13.16 MBH/,Yd  
ph=[|P)  
&b&o];a  
;------------------------------------------------------------------------------------------------- ;^:$O6J7T~  
;------------------------------------------------------------------------------------------------- &,@wLy^T  
def gen_Sando3 hk1jxnQh  
  command 5Ai$1'*p  
gen zone brick      p0 ttx1,0,-8     p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz12 & ,@*`2I>`  
                    p4 add tt2,yc,0  p5 add 0,yc,czz12  p6 add tt2,0,czz12  p7 add tt2,yc,czz12 &   J'y*>dW  
                    size 50 1 3 ;ratio 1.1 1 1 WP0{%  
T89VSB~  
        end_command H0i\#)Xs  
end f7QX"p&P  
gen_Sando3 EM.rO/qcW  
group sand_o3 range x 10.05 78.95 y 0 1 z -8 -5                                          ;8.26 11.26 f^X\N/  
;------------------------------------------------------------------------------------------------- uDi#a~m@  
def gen_Sando4 _a,XL<9I  
  command %uLyL4*L(p  
gen zone brick      p0 ttx1,0,-14.8   p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz13 & >~^##bIb  
                    p4 add tt2,yc,0   p5 add 0,yc,czz13  p6 add tt2,0,czz13  p7 add tt2,yc,czz13 &   ZI#Xh5  
                    size 50 1 6 ;ratio 1.1 1 1 Sf0[^"7  
a'ODm6#  
        end_command :7Q, `W9  
end XG}pp`{o  
gen_Sando4 Q CB~x2C  
group sand_o4 range x 10.05 78.95 y 0 1 z -14.8 -8                                   ;1.46 8.26 >zAI#N4  
3$ 1 z  
;------------------------------------------------------------------------------------------------- k|T0Bly3P  
def gen_Sando5 '$n#~/#}  
  command cZXra(AD  
gen zone brick      p0 ttx1,0,-16.26   p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz14 & >jDx-H.N  
                    p4 add tt2,yc,0    p5 add 0,yc,czz14  p6 add tt2,0,czz14  p7 add tt2,yc,czz14 &   !4G<&hvb  
                    size 50 1 2 ;ratio 1.1 1 1 Gk5SG_o  
H=k*;'  
        end_command &g<`i
{_  
end 8RR6f98FF  
gen_Sando5 Jv=G
3=.  
group sand_o5 range x 10.05 78.95 y 0 1 z -16.26 -14.8                            ;0 1.46 ;]^JUmxU[d  
;------------------------------------------------------------------------------------------------- 4\yKd8I  
;------------------------------------------------------------------------------------------------- ^@..\X9  
1)m&6:!b  
def gen_Sando6 da&f0m U  
  command I3V>VLv  
gen zone brick      p0 0,0,-50        p1 add ttx1,0,0    p2 add 0,yc,0     p3 add 0,0,xh & }~28UXb23  
                    p4 add ttx1,yc,0  p5 add 0,yc,xh     p6 add ttx1,0,xh  p7 add ttx1,yc,xh &   rfNt  
                    size 10 1 20 ;ratio 1.2 1 1 f*SAbDE  
gJ>HFid_C  
        end_command
g8_IZ(%:  
end Af"vSL  
gen_Sando6 &vp0zYd+v  
group sand_o6 range x 0 10.05  y 0 1 z  -50  -16.26 L)ry!BuHI  
;------------------------------------------------------------------------------------------------- 3eFBe2  
def gen_Sando6 #FV(a~  
  command $R<eXDW6:  
gen zone brick      p0 ttx1,0,-50    p1 add tt2,0,0  p2 add 0,yc,0    p3 add 0,0,xh & o<-+y\J8K  
                    p4 add tt2,yc,0  p5 add 0,yc,xh  p6 add tt2,0,xh  p7 add tt2,yc,xh &   DweWFipyPi  
                    size 40 1 20 ;ratio 1.1 1   SbW6O_   
\i#0:3s.  
        end_command ba  
end m=<;)  
gen_Sando6 O(E-ox~q  
group sand_o6 range  x 10.05 78.95 y 0 1  z -50   -16.26                      ;-33.74 0 XL7jUi_4:L  
;------------------------------------------------------------------------------------------------- zdUi1 b  
;------------------------------------------------------------------------------------------------- n`hes_{,g  
; 二、生成左侧土体 W=~H_L?/  
;------------------------------------------------------------------------------------------------- 8h2!8'  
;------------------------------------------------------------------------------------------------- 8W_X&X?Q  
def gen_Sand11 I:aG(8Bi)H  
  command XOy#?X/`  
gen zone brick      p0 0,0,ht1        p1 add 0,yc,0  p2 add ttx2,0,0      p3 add 0,0,ht2 & 9jwo f}OU  
                    p4 add ttx2,yc,0  p5 add ttx2,0,ht2  p6 add 0,yc,ht2  p7 add ttx2,yc,ht2 & 4hv'OEl  
                    size 1 10 3                               ;ratio 1 1 1 1.2 ?k:i3$  
]& qmV  
        end_command QYL ';  
end 5&r2a}K  
gen_Sand11 BOp&s>hI  
group sand_o7 range x -9.55 0 y 0 1 z -3.1 0 J;wA  
;------------------------------------------------------------------------------------------------- cAc i2e  
def gen_Sand2 (8(z42  
  command ~L'}!' &.  
gen zone brick      p0 ttx2,0,ht1    p1 add 0,yc,0   p2 add tt1,0,0     p3 add 0,0,ht2 & Zb7:qe<UN  
                    p4 add tt1,yc,0  p5 add tt1,0,ht2  p6 add 0,yc,ht2  p7 add tt1,yc,ht2 & v+*l|!v  
                    size 1 50  3 ;ratio 1 1.1 1 =JnUTc_u  
$'Hg}|53  
        end_command ico(4KSk  
end TGz5t$]I  
gen_Sand2 xQhvs=Zm]  
group sand_o7 range x -78.45 -9.55 y 0 1 z -3.1 0                   ?iBHJ{  
;------------------------------------------------------------------------------------------------- lFY;O !Y5\  
def gen_Sand31 2v<[XNX  
  command f V.(v&  
gen zone brick      p0 ttx2,0,-5     p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz11 & <!vAqqljt  
                    p4 add tt1,yc,0  p5 add tt1,0,czz11  p6 add 0,yc,czz11  p7 add tt1,yc,czz11 &   prY9SQd  
                    size 1 50 2 ;ratio 1 1.1 1 Uq6..<#  
]X)EO49  
        end_command G#4cWn'  
end ^$y_~z3o#7  
gen_Sand31 `&U ['_%  
group sand_o8 range x -78.45 -9.55 y 0 1 z -5 -3.1                                    ;11.26 13.16 VmB/X))  
;------------------------------------------------------------------------------------------------- gU}?Yy  
def gen_Sand32 (IR'~:W  
  command HiG&`:P>q  
gen zone brick      p0 ttx2,0,-8     p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz12 & k|7XC@i]%  
                    p4 add tt1,yc,0  p5 add tt1,0,czz12  p6 add 0,yc,czz12  p7 add tt1,yc,czz12 &   R%Yws2Le2  
                    size 1 50 3 ;ratio 1 1.1 1 z {J1pH_X  
?D6rFUs9;  
        end_command a;
Y9wn  
end Pz"!8b-MN  
gen_Sand32 '4A8\&lQO  
group sand_o9 range x -78.45 -9.55 y 0 1  z -8 -5                                          ;8.26 11.26 _dEf@==  
;------------------------------------------------------------------------------------------------- cZ7b$MZ%9  
def gen_Sand33 E 2DTE  
  command -j9R%+YW<  
gen zone brick      p0 ttx2,0,-14.8   p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz13 & KV0e^c;  
                    p4 add tt1,yc,0   p5 add tt1,0,czz13  p6 add 0,yc,czz13  p7 add tt1,yc,czz13 &   ]eq3cwR[|  
                    size 1 50 6 ;ratio 1 1.1 1 \(LHcvbb  
\0pJ+@\T9  
        end_command ca_8S8lv  
end WiL~b =fT  
gen_Sand33 UmU=3et<Wj  
group sand_o10 range x -78.45 -9.55 y 0 1  z -14.8 -8                                   ;1.46 8.26 VMW?[j  
y*6r&989  
;------------------------------------------------------------------------------------------------- ;.h5; `&  
def gen_Sand34 "t"=9:_t  
  command R@0ELxzA  
gen zone brick      p0 ttx2,0,-16.26   p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz14 & L$x/T3@  
                    p4 add tt1,yc,0    p5 add tt1,0,czz14  p6 add 0,yc,czz14  p7 add tt1,yc,czz14 &   T8FKa4ikn  
                    size 1 50 2 ;ratio 1 1.1 1 .n`MPx'  
'vTD7a^  
        end_command k>Qr
14F  
end OX4+1@$tk  
gen_Sand34 sh?Dxodp9  
group sand_o11 range x -78.45 -9.55 y 0 1  z -16.26 -14.8                            ;0 1.46 EQ>bwEG  
;------------------------------------------------------------------------------------------------- N3H!ptn37  
;------------------------------------------------------------------------------------------------- ?
^&!/,  
>}/"gx  
def gen_Sand4 ls6ywLP{  
  command 0nz k?iP  
gen zone brick      p0 0,0,-50        p1 add 0,yc,0  p2 add ttx2,0,0     p3 add 0,0,xh & s^9N7'  
                    p4 add ttx2,yc,0  p5 add ttx2,0,xh   p6 add 0,yc,xh  p7 add ttx2,yc,xh &   8L 9;VY^Y  
                    size 1 10 20 R#bg{|  
I=^%l7  
        end_command o=_4v^  
end )[)-.{q  
gen_Sand4 p$@=N6)I.k  
group sand_o12 range x -9.55 0 y 0 1 z -16.26 -50 4f"a
/(>*  
;------------------------------------------------------------------------------------------------- GKPqBi[rO  
def gen_Sand5 /BVNJNhz  
  command /kVy#sT|  
gen zone brick      p0 ttx2,0,-50    p1 add 0,yc,0  p2 add tt1,0,0    p3 add 0,0,xh & [:!#F7O-  
                    p4 add tt1,yc,0  p5 add tt1,0,xh  p6 add 0,yc,xh  p7 add tt1,yc,xh &   0r=Lilu{q  
                    size 1 50 20 ;ratio 1 1.1 1 ,9"</\]`  
s/Wg^(&M  
        end_command MnsnW{VGX  
end r/L
3j0  
gen_Sand5 TR@$$RrU  
group sand_o12 range x -78.45 -9.55 y 0 1 z -16.26 -50                        ;-33.74 0 >|7&hj$  
;------------------------------------------------------------------------------------------------- =k`(!r2"#  
call Flacmodel_haitang.dat ,Ql3RO,  
;------------------------------------------------------------------------------------------------- 4#=!VK8ZH  
plot add axes {R,rc!yF  
;-------------------------------------------------------------------------------------- );;
UNO21+  
;二、生成接触面的命令 @Cm"lv.hz  
;-------------------------------------------------------------------------------------- j&Ayk*  
interface 1 face range x -9.55 10.05 y 0 1 z -16.27 -16.25 cuquA ~  
interface 2 face range x -9.55 10.05 y 0 1 z -3.11 -3.09 wy_TFV  
interface 3 face range x -9.56 -9.54 y 0 1 z -5 -3.1 |_~BV&g,N  
interface 4 face range x -9.56 -9.54 y 0 1 z -8  -5 {\ziy4<II  
interface 5 face range x -9.56 -9.54 y 0 1 z -14.8 -8 s:cS 9A8  
interface 6 face range x -9.56 -9.54 y 0 1 z -16.26 -14.8 4!6g[[|&J  
0tB9X9:,  
int  1  prop kn  1.570e9  ks  1.570e9  fric  10 coh 10000                   ;prop属性 kn法向刚度 ks剪切刚度 fric为摩擦角 coh内聚力 ~%Yh`c EP  
int  2  prop kn  1.570e9  ks  1.570e9  fric  15 coh 55000 Zk}e?Grc  
int  3  prop kn  1.570e9  ks  1.570e9  fric  26 coh 8000 Z[`J'}?|  
int  4  prop kn  1.570e9  ks  1.570e9  fric  15 coh 55000 xMhR;lKY  
int  5  prop kn  1.570e9  ks  1.570e9  fric  14 coh 24000 hv .Mf.m  
int  6  prop kn  1.570e9  ks  1.570e9  fric  17 coh 56000 YKl!M/  
ini y add -2.0 rang group 1 $YaL3n  
;=============================================== Hv0sl+  
;初始应力计算程序 4DfTVO"h  
;=============================================== p9_45u`u2  
;定义模型及性质 Qq]UEI `Go  
model  mohr    rang  group sand_o1                                           miHW1h[=  
prop bulk 1.875e7   shear   0.967e7 coh 10000   fric 10 rang group sand_o1 '7'cKp  
model  mohr    rang  group sand_o2                                           VkhK2  
prop bulk 1.875e7   shear   0.967e7 coh 55000   fric 15 rang group sand_o2 \{+7`4g  
model  mohr    rang  group sand_o3                                           3g|O2>*?  
prop bulk 2.068e7   shear   0.793e7 coh 8000   fric 26 rang group sand_o3 foE2rV/Y  
model  mohr    rang  group sand_o4                                           EGGy0ly  
prop bulk 1.875e7   shear   0.967e7 coh 55000   fric 15 rang group sand_o4 4'9h^C&  
model  mohr    rang  group sand_o5                                           XW]|Mv[M  
prop bulk 1.396e7   shear   0.758e7 coh 24000   fric 14 rang group sand_o5 sS(^7GARa  
model  mohr    rang  group sand_o6                                           8yFD2(#  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o6 =GM!M@~,Ab  
model  mohr    rang  group sand_o7                                           Zml9ndzT  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o7 =$Q3!bJ  
model  mohr    rang  group sand_o8                                           !1"~tA!+p=  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o8 ,-DE;l^Q=  
model  mohr    rang  group sand_o9                                           
NM~e  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o9 ~5}b$qL#`  
model  mohr    rang  group sand_o10                                           "Jnq~7]  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o10 nhXa&Nro  
model  mohr    rang  group sand_o11                                           &"C1XM  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o11     rmQGzQnun  
model  mohr    rang  group sand_o12                                           #8|;Q`Or:  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o12 e@^}y4 C  
dGh<R|U3  
mo elas                 rang group 1 uNhAfZ  
;prop bulk 20.8e9  shear 15.58e9  rang group 1   ;正确 5'V'~Q%  
prop bulk 22.8e9  shear 17.08e9  rang group 1   ;只把柱子等效 ieS
5*@^k  
ini dens 2500 rang group 1 Punbw\9!d,  
ini dens 1900 rang group sand_o1     ;18000 .#WF'  
ini dens 1990 rang group sand_o2 PD/JXExK  
ini dens 1910 rang group sand_o3 '}4[m>/  
ini dens 1990 rang group sand_o4 B{1+0k  
ini dens 1970 rang group sand_o5 2#W%--  
ini dens 2010 rang group sand_o6 6x/ X8zu  
ini dens 2010 rang group sand_o7 )vGRfFjw_  
ini dens 2010 rang group sand_o8 pdnL~sv  
ini dens 2010 rang group sand_o9 GJy,)EO6{  
ini dens 2010 rang group sand_o10 N'
m:V  
ini dens 2010 rang group sand_o11 b<.+WkO  
ini dens 2010 rang group sand_o12 PLo.q|%  
;边界条件 ^ad> (W  
fix x y z rang z -50.1 -49.9 Z*]n]eS  
fix x rang x  78.85  79.05 6o A0a\G'  
fix x rang x -78.55 -78.35 /'VCJjzZ
 
fix y rang y -.1  .1 9R;s;2$.  
fix y rang y 0.9  1.1 ocgbBE  
;应力条件 [
EGE|  
ini szz 0  grad 0 0 20e3   rang z -50 0     ;加权平均重度为1993 ~T4=Id  
ini sxx 0  grad 0 0 10e3   rang z -50 0 $X*$,CCIB  
ini syy 0  grad 0 0 10e3   rang z -50 0 nNilTJ   
;--------------------------------------------------------------------------------- 4 <]QMA0
 
hist id 1 unbal rhn*kf{8  
plot hist 1 Cv$TNkP*  
; "v*RY "5#  
set grav 0 0 -10 cS ];?tqrA  
solve ratio 1.0e-4 {/0,lic  
save station1.sav

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

edelmann:土体强度、刚度参数给的太小了，你可以检查塑性区，不正常地大。 <&L;9fr  
我没试。 (2013-03-15 17:38)  SV\x2^Ea0  

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

参数的问题