首页 岩土在线 评论员 岩土论坛 注册岩土 工程市场 岩土问答 岩土人才网 岩土领域 岩土企业
论坛首页    职业区    学术与写作    工程技术区    软件区    资料区    商务合作区    社区办公室
 您好！ 欢迎 登录 或 注册 最新帖子 邀请注册 活动聚焦 统计排行 社区服务 帮助
 
  • 帖子
  • 日志
  • 用户
  • 版块
  • 群组
帖子
热搜： 注册岩土规范勘察基坑
发帖 回复
返回列表
  • 2150阅读
  • 3回复

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

 上一主题 下一主题
离线chuhachian
实习生
 

发帖
12
土币
60
威望
9
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2013-03-15
uzb|yV'B  
new } PL{i  
config dyn 8LF=l1=~  
set dyn off #lik: ?  
;=============================================== ,Y| ;V  
;建模型 ' Hj([N  
;=============================================== G,+3(C  
def setup              ;定义函数setup fg ,vTpBk  
   yc=1               ;y方向尺寸 Ya~Th)'>q  
   czz1=-13.16 *'?V>q,  
   czz11=1.9            ;四分之一车站z方向尺寸   Y_C6*T%  
   czz12=3            ;四分之一车站z方向尺寸 1}Guhayy  
   czz13=6.8            ;四分之一车站z方向尺寸 uMm`j?Y23q  
   czz14=1.46            ;四分之一车站z方向尺寸     9]T61Z{OW1  
   ttx1=10.05              ;四分之一车站外围土体x方向尺寸 {|0YcL  
   ttx2=-9.55              ;四分之一车站外围土体-x方向尺寸   -|u yJh  
   ht1=-3.1 a; a1>1  
   ht2=3.1 nm_taER  
   tt2=68.9 ;78.95             }s"].Xm^2  
   tt1=-68.9 U:@tdH+A7  
   xh=33.74              ；车站以下土体高度 C \5yo  
   jT]R"U/Q  
end ~("bpS#ZgD  
setup   ?N9Z;_&^.  
;------------------------------------------------------------------------------------------------- -ert42fN  
; 一、生成右侧土体 mQt0?c _  
;------------------------------------------------------------------------------------------------- CX2qtI8N?  
def gen_Sando1 PB*G#2W  
  command FQ 0 ;%Z  
gen zone brick      p0 0,0,ht1   p1 add ttx1,0,0  p2 add 0,yc,0  p3 add 0,0,ht2 & ./zzuKO8XK  
                    p4 add ttx1,yc,0  p5 add 0,yc,ht2  p6 add ttx1,0,ht2  p7 add ttx1,yc,ht2 & [j TU nP  
                    size 10 1 3                               ;ratio 1 1 1 1.2 InRRcn(  
w`< {   
        end_command W@z xGH$z>  
end @+ T33X)h%  
gen_Sando1 2^=.f?_YR  
group sand_o1 range x 0 10.05  y 0 1  z -3.1 0 I d8MXdV  
;-------------------------------------------------------------------------------------------------  6shN%  
def gen_Sando1 Jw;G_dQ[  
  command ;P}007;  
gen zone brick      p0 ttx1,0,ht1   p1 add tt2,0,0  p2 add 0,yc,0  p3 add 0,0,ht2 & eC<?g  
                    p4 add tt2,yc,0  p5 add 0,yc,ht2  p6 add tt2,0,ht2  p7 add tt2,yc,ht2 & zRgGSxn  
                    size 50 1 3 ;ratio 1.1 1 1 S&&Q U #  
Jq"3xj   
        end_command rmC7!^/  
end vV=rBO0a?  
gen_Sando1 }4piZ ch  
group sand_o1 range x 10.05 78.95 y 0 1 z -3.1 0                   [5!{>L`  
;------------------------------------------------------------------------------------------------- (i\{hq/  
def gen_Sando2 pKLNBR|  
  command OrL4G `O  
gen zone brick      p0 ttx1,0,-5   p1 add tt2,0,0  p2 add 0,yc,0  p3 add 0,0,czz11 & xY.?OHgG/  
                    p4 add tt2,yc,0  p5 add 0,yc,czz11  p6 add tt2,0,czz11  p7 add tt2,yc,czz11 &   `|&0j4(Pg  
                    size 50 1 2 ;ratio 1.1 1 1 *>:<  
@o1#J` rv  
        end_command yK"HHdYTV  
end (]?M=?0\  
gen_Sando2 "9X!Ewm"P  
group sand_o2 range x 10.05 78.95 y 0 1 z -5 -3.1                                    ;11.26 13.16  6cjCn  
Pd;8<UMk  
*q\>DE=7  
;------------------------------------------------------------------------------------------------- x1Z'_Qw  
;------------------------------------------------------------------------------------------------- (+xT5 2  
def gen_Sando3 4M0p:Ey '  
  command mBB"e"o  
gen zone brick      p0 ttx1,0,-8     p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz12 & RkTYvAk|kY  
                    p4 add tt2,yc,0  p5 add 0,yc,czz12  p6 add tt2,0,czz12  p7 add tt2,yc,czz12 &   hW~UJ/$  
                    size 50 1 3 ;ratio 1.1 1 1 y"8,jm  
<e S+3,  
        end_command Xwu&K8q21  
end tCRsaDK>  
gen_Sando3 A|4 3W =  
group sand_o3 range x 10.05 78.95 y 0 1 z -8 -5                                          ;8.26 11.26 A"qDc  
;------------------------------------------------------------------------------------------------- aMT=pGU  
def gen_Sando4 #RyTa /L  
  command C]3:&dx9  
gen zone brick      p0 ttx1,0,-14.8   p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz13 & )Pc>+} D  
                    p4 add tt2,yc,0   p5 add 0,yc,czz13  p6 add tt2,0,czz13  p7 add tt2,yc,czz13 &   {md5G$* %  
                    size 50 1 6 ;ratio 1.1 1 1 0k_3]Li=(  
}Q\+w,pJgN  
        end_command `PeC,bp  
end YUTh*`1k<  
gen_Sando4 g-u4E^,*|  
group sand_o4 range x 10.05 78.95 y 0 1 z -14.8 -8                                   ;1.46 8.26 pVzr]WFx  
`SZ-o{  
;------------------------------------------------------------------------------------------------- +~:OUR*>  
def gen_Sando5 r? }|W2^%  
  command CRiqY_gBf  
gen zone brick      p0 ttx1,0,-16.26   p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz14 & !?J- Y  
                    p4 add tt2,yc,0    p5 add 0,yc,czz14  p6 add tt2,0,czz14  p7 add tt2,yc,czz14 &   e\-,e+  
                    size 50 1 2 ;ratio 1.1 1 1 5-H"{29  
lqO>Q1_{K  
        end_command 8$RiFD ,  
end A@Zqh<,Ud  
gen_Sando5 0"GLgj:9  
group sand_o5 range x 10.05 78.95 y 0 1 z -16.26 -14.8                            ;0 1.46 M+j*5wNy  
;------------------------------------------------------------------------------------------------- $Fi1Bv)  
;------------------------------------------------------------------------------------------------- r}(mjC"o  
b?!S$Sxz  
def gen_Sando6 e%)MIAS0  
  command ;;C2t&(  
gen zone brick      p0 0,0,-50        p1 add ttx1,0,0    p2 add 0,yc,0     p3 add 0,0,xh & 6#qt%t%?D  
                    p4 add ttx1,yc,0  p5 add 0,yc,xh     p6 add ttx1,0,xh  p7 add ttx1,yc,xh &   uvR l`"Y  
                    size 10 1 20 ;ratio 1.2 1 1 *)?'!  
*c%{b3T_  
        end_command "~zLG"  
end oLq N  
gen_Sando6 UxF9Ko( ]d  
group sand_o6 range x 0 10.05  y 0 1 z  -50  -16.26 '6g-]rE[  
;------------------------------------------------------------------------------------------------- sV0NDM0  
def gen_Sando6 M$!-B,1BX  
  command GjN6Af~}  
gen zone brick      p0 ttx1,0,-50    p1 add tt2,0,0  p2 add 0,yc,0    p3 add 0,0,xh & {KK/mAp{  
                    p4 add tt2,yc,0  p5 add 0,yc,xh  p6 add tt2,0,xh  p7 add tt2,yc,xh &   92C; a5s  
                    size 40 1 20 ;ratio 1.1 1   ZP@NV|B  
7hLh}  
        end_command De{ZQg)  
end >o3R~ [  
gen_Sando6 .!+7|us8l\  
group sand_o6 range  x 10.05 78.95 y 0 1  z -50   -16.26                      ;-33.74 0 2qVoe}F  
;------------------------------------------------------------------------------------------------- 9pr.`w  
;------------------------------------------------------------------------------------------------- @KHY8y7  
; 二、生成左侧土体 f;OB"p  
;------------------------------------------------------------------------------------------------- o!&+ _BKw  
;------------------------------------------------------------------------------------------------- 4hfq7kq7(  
def gen_Sand11 Vo.~1^  
  command O~?d;.b  
gen zone brick      p0 0,0,ht1        p1 add 0,yc,0  p2 add ttx2,0,0      p3 add 0,0,ht2 & fo~*Bp()-E  
                    p4 add ttx2,yc,0  p5 add ttx2,0,ht2  p6 add 0,yc,ht2  p7 add ttx2,yc,ht2 & %h,&ND  
                    size 1 10 3                               ;ratio 1 1 1 1.2 M^ e}w!U  
eR`Q7]j] -  
        end_command 5yj#9H  
end 48 0M|^  
gen_Sand11 OTAe#]#  
group sand_o7 range x -9.55 0 y 0 1 z -3.1 0 amX1idHo^  
;------------------------------------------------------------------------------------------------- O:~J_Wwl!  
def gen_Sand2 1D!MXYgm1b  
  command Nq6; z)$  
gen zone brick      p0 ttx2,0,ht1    p1 add 0,yc,0   p2 add tt1,0,0     p3 add 0,0,ht2 & WjSu4   
                    p4 add tt1,yc,0  p5 add tt1,0,ht2  p6 add 0,yc,ht2  p7 add tt1,yc,ht2 & !&.-{ _$  
                    size 1 50  3 ;ratio 1 1.1 1 |TQ4:P1T  
i6P$>8jBQ-  
        end_command =\MAz[IDj  
end H?ug-7k/  
gen_Sand2 mQSn*;9\T3  
group sand_o7 range x -78.45 -9.55 y 0 1 z -3.1 0                   YRv96|c,  
;------------------------------------------------------------------------------------------------- )%kiM<})  
def gen_Sand31 W|E %  
  command @Jqo'\~&  
gen zone brick      p0 ttx2,0,-5     p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz11 & AOwmPHEL  
                    p4 add tt1,yc,0  p5 add tt1,0,czz11  p6 add 0,yc,czz11  p7 add tt1,yc,czz11 &   M0?%r`  
                    size 1 50 2 ;ratio 1 1.1 1 IAN={";p  
iweT @P`  
        end_command ([^f1;ncm  
end XWNo)#_3  
gen_Sand31 [}l 90lP  
group sand_o8 range x -78.45 -9.55 y 0 1 z -5 -3.1                                    ;11.26 13.16 2AMb-&po&f  
;------------------------------------------------------------------------------------------------- E.0J94>iM  
def gen_Sand32 QctzIC#;k  
  command `|v/qk7 ^?  
gen zone brick      p0 ttx2,0,-8     p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz12 & Jk7 Am-.0  
                    p4 add tt1,yc,0  p5 add tt1,0,czz12  p6 add 0,yc,czz12  p7 add tt1,yc,czz12 &   D6fd(=t1Z  
                    size 1 50 3 ;ratio 1 1.1 1 rz`"$g+#  
p&4n"hC  
        end_command D#`>p  
end x&+&)d  
gen_Sand32 0%q H=do6  
group sand_o9 range x -78.45 -9.55 y 0 1  z -8 -5                                          ;8.26 11.26 D dCcsYm,  
;------------------------------------------------------------------------------------------------- }L%2K"8?}  
def gen_Sand33 *XYp~b  
  command 4b, +;  
gen zone brick      p0 ttx2,0,-14.8   p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz13 & darbL_1  
                    p4 add tt1,yc,0   p5 add tt1,0,czz13  p6 add 0,yc,czz13  p7 add tt1,yc,czz13 &   oIj -Y`92!  
                    size 1 50 6 ;ratio 1 1.1 1 5}! 36SO\  
oSjYp(h:  
        end_command , )TnIByM  
end 0ZLLbEfnPB  
gen_Sand33 %]4=D)Om  
group sand_o10 range x -78.45 -9.55 y 0 1  z -14.8 -8                                   ;1.46 8.26 QoxQ"r9Wh  
jY=M{?h''  
;------------------------------------------------------------------------------------------------- MR5[|kHJT  
def gen_Sand34 Jt<J#M<}7  
  command %BT]h3dcSS  
gen zone brick      p0 ttx2,0,-16.26   p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz14 & 5')]Y1J  
                    p4 add tt1,yc,0    p5 add tt1,0,czz14  p6 add 0,yc,czz14  p7 add tt1,yc,czz14 &   u~JR]T  
                    size 1 50 2 ;ratio 1 1.1 1 ]Wjcr2Wq  
a({N}ZDo  
        end_command ;R<V-gab  
end Ro `Xs.X  
gen_Sand34 ,!PV0(F(  
group sand_o11 range x -78.45 -9.55 y 0 1  z -16.26 -14.8                            ;0 1.46 Nz*,m'-1e  
;------------------------------------------------------------------------------------------------- B&1E&Cv_8  
;------------------------------------------------------------------------------------------------- -II03 S1  
@[f$MRp\  
def gen_Sand4 *WFd[cKE  
  command hM>.xr  
gen zone brick      p0 0,0,-50        p1 add 0,yc,0  p2 add ttx2,0,0     p3 add 0,0,xh & >Y+m54EE  
                    p4 add ttx2,yc,0  p5 add ttx2,0,xh   p6 add 0,yc,xh  p7 add ttx2,yc,xh &   8TU(5:xJo  
                    size 1 10 20 gNDMJ^`  
zuFPG{^\#  
        end_command DWCf+4  
end qzO5p=}  
gen_Sand4 >M##q?.  
group sand_o12 range x -9.55 0 y 0 1 z -16.26 -50 z[O*f#t  
;------------------------------------------------------------------------------------------------- F ~7TE91C  
def gen_Sand5 vCK+v r!  
  command 5DkEJk7a  
gen zone brick      p0 ttx2,0,-50    p1 add 0,yc,0  p2 add tt1,0,0    p3 add 0,0,xh & o Q*LP{M  
                    p4 add tt1,yc,0  p5 add tt1,0,xh  p6 add 0,yc,xh  p7 add tt1,yc,xh &   "3a}~J<g  
                    size 1 50 20 ;ratio 1 1.1 1 tGbx/$Y   
ufw[Ei$I:  
        end_command 6V@_?a-K  
end s5Wb iOF  
gen_Sand5 @6aJh< c  
group sand_o12 range x -78.45 -9.55 y 0 1 z -16.26 -50                        ;-33.74 0 K_:2sDCaN  
;------------------------------------------------------------------------------------------------- T5I#7LN#  
call Flacmodel_haitang.dat rce._w }  
;------------------------------------------------------------------------------------------------- ?4aW^l6/  
plot add axes a"t~ K  
;-------------------------------------------------------------------------------------- %q9"2] cR  
;二、生成接触面的命令 1A#/70Mo  
;-------------------------------------------------------------------------------------- T2tvU*[=  
interface 1 face range x -9.55 10.05 y 0 1 z -16.27 -16.25 OQKc_z'"  
interface 2 face range x -9.55 10.05 y 0 1 z -3.11 -3.09 X8R:9q_  
interface 3 face range x -9.56 -9.54 y 0 1 z -5 -3.1 ,q7FK z{  
interface 4 face range x -9.56 -9.54 y 0 1 z -8  -5 59"tHb6E  
interface 5 face range x -9.56 -9.54 y 0 1 z -14.8 -8 EQw7(r|v:  
interface 6 face range x -9.56 -9.54 y 0 1 z -16.26 -14.8 >LH}A6dUC  
Di}M\!-[  
int  1  prop kn  1.570e9  ks  1.570e9  fric  10 coh 10000                   ;prop属性 kn法向刚度 ks剪切刚度 fric为摩擦角 coh内聚力 c6h+8QS  
int  2  prop kn  1.570e9  ks  1.570e9  fric  15 coh 55000 F?cwIE\J  
int  3  prop kn  1.570e9  ks  1.570e9  fric  26 coh 8000 ;+#Nb/M  
int  4  prop kn  1.570e9  ks  1.570e9  fric  15 coh 55000 rO{?.#~  
int  5  prop kn  1.570e9  ks  1.570e9  fric  14 coh 24000 % -+7=x  
int  6  prop kn  1.570e9  ks  1.570e9  fric  17 coh 56000 8Z "f"  
ini y add -2.0 rang group 1 3)2{c  
;=============================================== J6>tGKa+e  
;初始应力计算程序 wf\7sz  
;=============================================== _%\%  
;定义模型及性质 BgDWl{pm  
model  mohr    rang  group sand_o1                                           EAxdF u  
prop bulk 1.875e7   shear   0.967e7 coh 10000   fric 10 rang group sand_o1 :V@)A/}uk  
model  mohr    rang  group sand_o2                                           T{-gbo`Yji  
prop bulk 1.875e7   shear   0.967e7 coh 55000   fric 15 rang group sand_o2 + 660/ e8N  
model  mohr    rang  group sand_o3                                           lk R^2P  
prop bulk 2.068e7   shear   0.793e7 coh 8000   fric 26 rang group sand_o3 (ov&iNx  
model  mohr    rang  group sand_o4                                           Of$R+n.  
prop bulk 1.875e7   shear   0.967e7 coh 55000   fric 15 rang group sand_o4 )~ &gBX  
model  mohr    rang  group sand_o5                                           XJy.xI>;  
prop bulk 1.396e7   shear   0.758e7 coh 24000   fric 14 rang group sand_o5 ab.B?bx  
model  mohr    rang  group sand_o6                                           0_Elxc  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o6 nG{o$v_|  
model  mohr    rang  group sand_o7                                           fBz|-I:k +  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o7 5~im.XfiVx  
model  mohr    rang  group sand_o8                                           @0C[o9  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o8 dV}]\ 8N  
model  mohr    rang  group sand_o9                                           CPeu="[  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o9 \1n (Jr.<  
model  mohr    rang  group sand_o10                                           &@BAVc z  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o10 9Nx%Sdu  
model  mohr    rang  group sand_o11                                           Ai^0{kF6  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o11     l%?4L/J)#  
model  mohr    rang  group sand_o12                                           Zg>]!^X8  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o12  ylS6D  
,w9| ?%S  
mo elas                 rang group 1 guf*>qNr  
;prop bulk 20.8e9  shear 15.58e9  rang group 1   ;正确 DO+~    
prop bulk 22.8e9  shear 17.08e9  rang group 1   ;只把柱子等效 .)p%|A#^  
ini dens 2500 rang group 1 ]:']  
ini dens 1900 rang group sand_o1     ;18000 -AolW+Y  
ini dens 1990 rang group sand_o2 kCoE;)y$  
ini dens 1910 rang group sand_o3 y9LO;{(  
ini dens 1990 rang group sand_o4 ]%FP*YU4O  
ini dens 1970 rang group sand_o5 M&gi$Qs[E  
ini dens 2010 rang group sand_o6 @,c` #,F/  
ini dens 2010 rang group sand_o7 ^Pu:&:ki  
ini dens 2010 rang group sand_o8 KK6z3"tk5  
ini dens 2010 rang group sand_o9 $d4&H/u^  
ini dens 2010 rang group sand_o10 >msQ@Ch  
ini dens 2010 rang group sand_o11 ^K_FGE0ec  
ini dens 2010 rang group sand_o12 )54a' Hp  
;边界条件 h;y}g/HZ  
fix x y z rang z -50.1 -49.9 qK2jJ3)>  
fix x rang x  78.85  79.05 Qe4 % A  
fix x rang x -78.55 -78.35 Hi/[  
fix y rang y -.1  .1 > LN*3&W  
fix y rang y 0.9  1.1 V\e1NS  
;应力条件 ._<, Eodv  
ini szz 0  grad 0 0 20e3   rang z -50 0     ;加权平均重度为1993 ^,5%fl  
ini sxx 0  grad 0 0 10e3   rang z -50 0 +uTl Lu;MT  
ini syy 0  grad 0 0 10e3   rang z -50 0 #`K{vj  
;--------------------------------------------------------------------------------- )l! `k  
hist id 1 unbal ;_|4c7  
plot hist 1 >Bdh`Ot-!  
; 6U$e;cr6  
set grav 0 0 -10 HD2C^V2@M  
solve ratio 1.0e-4 \Y8 sIs  
save station1.sav
回复
举报
分享到 淘江湖 新浪 QQ微博 QQ空间 开心 人人 豆瓣 网易微博 百度 鲜果 白社会 飞信
离线edelmann
总工

发帖
2978
土币
53316
威望
7629
原创币
0
只看该作者 1楼 发表于: 2013-03-15
土体强度、刚度参数给的太小了，你可以检查塑性区，不正常地大。 q4G$I?4  
我没试。
回复
举报
离线chuhachian
实习生

发帖
12
土币
60
威望
9
原创币
0
只看该作者 2楼 发表于: 2013-03-24
回 edelmann 的帖子
edelmann:土体强度、刚度参数给的太小了，你可以检查塑性区，不正常地大。 gy~2LY!}  
我没试。 (2013-03-15 17:38)  =1(7T.t  
i_l{#*t  
土体强度指的哪个值？刚度参数是接触面上的刚度还是说体积模量和剪切模量？请赐教，因为我刚学习这个  还不是很懂
回复
举报
离线wangdan
助理工程师

发帖
87
土币
368
威望
47
原创币
0
只看该作者 3楼 发表于: 2014-05-20
参数的问题
回复
举报
发帖 回复
返回列表
快速回复
限100 字节
温馨提示：欢迎交流讨论，请勿纯表情、纯引用!
 
上一个 下一个

      https://beian.mps.gov.cn/ 粤公网安备 44010602012919号 广州半山岩土网络科技有限公司 粤ICP备2024274469号

      工业和信息化部备案管理系统网站