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

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

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

发帖
12
土币
60
威望
9
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2013-03-15
vs j3  
new Urgtg37  
config dyn 5IKL#V `3a  
set dyn off '>(.%@  
;=============================================== j:bgR8 %e  
;建模型 Yh)Isg|0>  
;=============================================== Em7 WDu0  
def setup              ;定义函数setup Xq4|uuS-O  
   yc=1               ;y方向尺寸 F6hmku>\1  
   czz1=-13.16 4m-I5!=O  
   czz11=1.9            ;四分之一车站z方向尺寸   auWXgkwZs/  
   czz12=3            ;四分之一车站z方向尺寸 rbZ[!LA  
   czz13=6.8            ;四分之一车站z方向尺寸 ?U7) XvQ  
   czz14=1.46            ;四分之一车站z方向尺寸     -]:G L>b  
   ttx1=10.05              ;四分之一车站外围土体x方向尺寸 tM j1~ R  
   ttx2=-9.55              ;四分之一车站外围土体-x方向尺寸   0L^u2HZYL  
   ht1=-3.1 {k(eNr,  
   ht2=3.1 8I,QD` xu  
   tt2=68.9 ;78.95             ]8;n{ }X  
   tt1=-68.9 yA[({2%  
   xh=33.74              ；车站以下土体高度 Pj?Dmk~   
   CCDDK L]N:  
end !Ss HAE|  
setup   !q"CV  
;------------------------------------------------------------------------------------------------- &2I*0  
; 一、生成右侧土体 Wn|&cG9  
;------------------------------------------------------------------------------------------------- A4mSJ6K]  
def gen_Sando1 kVS?RHR  
  command W >(vYU  
gen zone brick      p0 0,0,ht1   p1 add ttx1,0,0  p2 add 0,yc,0  p3 add 0,0,ht2 & sn{AwF%  
                    p4 add ttx1,yc,0  p5 add 0,yc,ht2  p6 add ttx1,0,ht2  p7 add ttx1,yc,ht2 & VK3it3FI>3  
                    size 10 1 3                               ;ratio 1 1 1 1.2 P6U%=xaC  
W2wpcc  
        end_command [*m2  
end FnHi(S|A  
gen_Sando1 ]|:uU  
group sand_o1 range x 0 10.05  y 0 1  z -3.1 0 8X&Ya =  
;------------------------------------------------------------------------------------------------- hF&}lPVtv  
def gen_Sando1 <JI& {1  
  command \0FT!} L  
gen zone brick      p0 ttx1,0,ht1   p1 add tt2,0,0  p2 add 0,yc,0  p3 add 0,0,ht2 & `&$B3)Eb  
                    p4 add tt2,yc,0  p5 add 0,yc,ht2  p6 add tt2,0,ht2  p7 add tt2,yc,ht2 & ~=y3Gd B3  
                    size 50 1 3 ;ratio 1.1 1 1 mE#nU(+Ta  
yy(A(}  
        end_command gd^Js 1Z  
end /2NSZO  
gen_Sando1 ZYtiMBJ  
group sand_o1 range x 10.05 78.95 y 0 1 z -3.1 0                   gY9\o#)<  
;------------------------------------------------------------------------------------------------- rU#li0 >  
def gen_Sando2 vi {uy  
  command ?Hy+'sq[  
gen zone brick      p0 ttx1,0,-5   p1 add tt2,0,0  p2 add 0,yc,0  p3 add 0,0,czz11 & $(HjI \%l^  
                    p4 add tt2,yc,0  p5 add 0,yc,czz11  p6 add tt2,0,czz11  p7 add tt2,yc,czz11 &   CP` XUpX`&  
                    size 50 1 2 ;ratio 1.1 1 1 S-Ai3)t6  
5^)_B;.f  
        end_command 3#.\  
end o7xgRSz\  
gen_Sando2 q^bO*bv  
group sand_o2 range x 10.05 78.95 y 0 1 z -5 -3.1                                    ;11.26 13.16 @Z\~  
aU<s<2 O)  
]HvZ$  
;------------------------------------------------------------------------------------------------- !Ua&0s%  
;------------------------------------------------------------------------------------------------- 3x5!a5$Y  
def gen_Sando3 Dl;d33  
  command YL^Z4: p  
gen zone brick      p0 ttx1,0,-8     p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz12 & V_a)jJ  
                    p4 add tt2,yc,0  p5 add 0,yc,czz12  p6 add tt2,0,czz12  p7 add tt2,yc,czz12 &   F!8=FTb  
                    size 50 1 3 ;ratio 1.1 1 1 QB1M3b  
MTmO>V&O  
        end_command uu}-"/<~7  
end Jyu`-=It  
gen_Sando3 6[==BbZ  
group sand_o3 range x 10.05 78.95 y 0 1 z -8 -5                                          ;8.26 11.26 ?q%b*Ek  
;------------------------------------------------------------------------------------------------- am:.NG+  
def gen_Sando4 G :4;y7  
  command (K!4Kp^m  
gen zone brick      p0 ttx1,0,-14.8   p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz13 & Gc;B[/:  
                    p4 add tt2,yc,0   p5 add 0,yc,czz13  p6 add tt2,0,czz13  p7 add tt2,yc,czz13 &   %D^bah f  
                    size 50 1 6 ;ratio 1.1 1 1 En5!"w|j  
Bxv8RB  
        end_command ="B n=>  
end 6An{3 "  
gen_Sando4 }aZr ou3E  
group sand_o4 range x 10.05 78.95 y 0 1 z -14.8 -8                                   ;1.46 8.26 #?m{YT{P  
~n]NyVFP  
;------------------------------------------------------------------------------------------------- 0$2={s4ze  
def gen_Sando5 %c1FwAC  
  command /J;]u3e|  
gen zone brick      p0 ttx1,0,-16.26   p1 add tt2,0,0     p2 add 0,yc,0       p3 add 0,0,czz14 & v>at/ef  
                    p4 add tt2,yc,0    p5 add 0,yc,czz14  p6 add tt2,0,czz14  p7 add tt2,yc,czz14 &   7!- \L7<  
                    size 50 1 2 ;ratio 1.1 1 1 pbdF]>\  
=>YvA>izE  
        end_command (9z|a ,  
end I*c;hfu  
gen_Sando5 "/y|VTV"  
group sand_o5 range x 10.05 78.95 y 0 1 z -16.26 -14.8                            ;0 1.46 E$T(Qu<-  
;------------------------------------------------------------------------------------------------- a4u^f5)@  
;------------------------------------------------------------------------------------------------- A`C-sD >  
yiO31uQt  
def gen_Sando6 B"[{]GP BY  
  command cUC17z2D  
gen zone brick      p0 0,0,-50        p1 add ttx1,0,0    p2 add 0,yc,0     p3 add 0,0,xh & L?.7\a@  
                    p4 add ttx1,yc,0  p5 add 0,yc,xh     p6 add ttx1,0,xh  p7 add ttx1,yc,xh &   h60\ Y 8  
                    size 10 1 20 ;ratio 1.2 1 1 \-G5l+!  
MHNe>C-!q  
        end_command =~$U^IsWA  
end pVz pN8!  
gen_Sando6 e/P4mc)  
group sand_o6 range x 0 10.05  y 0 1 z  -50  -16.26 ^yD"d =z  
;------------------------------------------------------------------------------------------------- ?.j,Bq5At  
def gen_Sando6 :f]!O@.~  
  command -3`Isv  
gen zone brick      p0 ttx1,0,-50    p1 add tt2,0,0  p2 add 0,yc,0    p3 add 0,0,xh & UU"d_~pp  
                    p4 add tt2,yc,0  p5 add 0,yc,xh  p6 add tt2,0,xh  p7 add tt2,yc,xh &   :o<N!*pT  
                    size 40 1 20 ;ratio 1.1 1   V^ Y*xZ  
I1~G$)w#  
        end_command }8 \|1@09  
end @92gb$xT  
gen_Sando6 RSFJu\0}N  
group sand_o6 range  x 10.05 78.95 y 0 1  z -50   -16.26                      ;-33.74 0 x57'Cg \  
;------------------------------------------------------------------------------------------------- gb9[Meg'  
;------------------------------------------------------------------------------------------------- N[]U%9[=2F  
; 二、生成左侧土体 `:R-[>5P8  
;------------------------------------------------------------------------------------------------- ^^'[%ok  
;------------------------------------------------------------------------------------------------- tNYJQ  
def gen_Sand11 &R0OeRToUb  
  command ,?fN#gc :  
gen zone brick      p0 0,0,ht1        p1 add 0,yc,0  p2 add ttx2,0,0      p3 add 0,0,ht2 & |9x%gUm  
                    p4 add ttx2,yc,0  p5 add ttx2,0,ht2  p6 add 0,yc,ht2  p7 add ttx2,yc,ht2 & T[m ~6  
                    size 1 10 3                               ;ratio 1 1 1 1.2 wYa0hNd  
CR.bMF}  
        end_command uH0#rgKt  
end  .?70=8{  
gen_Sand11 | [p68v>  
group sand_o7 range x -9.55 0 y 0 1 z -3.1 0 Io09W^  
;------------------------------------------------------------------------------------------------- v'K % %z  
def gen_Sand2 2h5tBEOX.s  
  command `)LIVi"(D  
gen zone brick      p0 ttx2,0,ht1    p1 add 0,yc,0   p2 add tt1,0,0     p3 add 0,0,ht2 & DqRLx85d1  
                    p4 add tt1,yc,0  p5 add tt1,0,ht2  p6 add 0,yc,ht2  p7 add tt1,yc,ht2 & {$3j/b  
                    size 1 50  3 ;ratio 1 1.1 1 #fO*ROe  
d RHlx QUn  
        end_command `~0P[>|+  
end Ko@zk<~"[  
gen_Sand2 FBYll[8  
group sand_o7 range x -78.45 -9.55 y 0 1 z -3.1 0                   R9bhC9NP  
;------------------------------------------------------------------------------------------------- w< v1 N  
def gen_Sand31 C+T&O  
  command Q +R3H,  
gen zone brick      p0 ttx2,0,-5     p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz11 & #"|"cYi,  
                    p4 add tt1,yc,0  p5 add tt1,0,czz11  p6 add 0,yc,czz11  p7 add tt1,yc,czz11 &   OuMco+C  
                    size 1 50 2 ;ratio 1 1.1 1 10/x'#(  
G OH  
        end_command QCVsVG!sN  
end -*rHB&e  
gen_Sand31 te4F"SEf  
group sand_o8 range x -78.45 -9.55 y 0 1 z -5 -3.1                                    ;11.26 13.16 Nvi Fq  
;------------------------------------------------------------------------------------------------- !T)T_P[  
def gen_Sand32 mu 2 A%"7  
  command <'yf|N!9G  
gen zone brick      p0 ttx2,0,-8     p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz12 & B:B8"ODV  
                    p4 add tt1,yc,0  p5 add tt1,0,czz12  p6 add 0,yc,czz12  p7 add tt1,yc,czz12 &   t1G1(F#&%  
                    size 1 50 3 ;ratio 1 1.1 1 E"G:K`Q  
(zjz]@qJ  
        end_command nXv 7OEpTx  
end M ' a&  
gen_Sand32 20$F$YYuk  
group sand_o9 range x -78.45 -9.55 y 0 1  z -8 -5                                          ;8.26 11.26 ,:Rft  
;------------------------------------------------------------------------------------------------- ;r gH}r  
def gen_Sand33 x%_qJ]o  
  command  )GB3=@  
gen zone brick      p0 ttx2,0,-14.8   p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz13 & dDsjPM;2  
                    p4 add tt1,yc,0   p5 add tt1,0,czz13  p6 add 0,yc,czz13  p7 add tt1,yc,czz13 &   A4~D#V  
                    size 1 50 6 ;ratio 1 1.1 1 Vf"O/o}hq,  
Uzan7A  
        end_command 1M7\:te*  
end c-[Q,c  
gen_Sand33 z\xiACIc  
group sand_o10 range x -78.45 -9.55 y 0 1  z -14.8 -8                                   ;1.46 8.26 _8,vk-,'  
A2}Z *U(;  
;------------------------------------------------------------------------------------------------- l*F!~J3  
def gen_Sand34 fR+Ov8PCq  
  command qf_h b  
gen zone brick      p0 ttx2,0,-16.26   p1 add 0,yc,0     p2 add tt1,0,0       p3 add 0,0,czz14 & Qw3a"k-  
                    p4 add tt1,yc,0    p5 add tt1,0,czz14  p6 add 0,yc,czz14  p7 add tt1,yc,czz14 &   #[gcg]6c  
                    size 1 50 2 ;ratio 1 1.1 1 &&]"Y!r -  
h5LJij J  
        end_command 3g?MEM~  
end ?)Tz'9l  
gen_Sand34 *,O :>Z5I  
group sand_o11 range x -78.45 -9.55 y 0 1  z -16.26 -14.8                            ;0 1.46 Sf+(1_^`t  
;------------------------------------------------------------------------------------------------- hEyX~f  
;------------------------------------------------------------------------------------------------- C ffTv  
Qh)|FQ[s$r  
def gen_Sand4 EpFIKV!  
  command I ybl;u  
gen zone brick      p0 0,0,-50        p1 add 0,yc,0  p2 add ttx2,0,0     p3 add 0,0,xh & 8L]gQ g  
                    p4 add ttx2,yc,0  p5 add ttx2,0,xh   p6 add 0,yc,xh  p7 add ttx2,yc,xh &   {{^Mr)]5K  
                    size 1 10 20 IeR l6r%:  
W /~||s  
        end_command 3Eb nZb  
end XA&tTpfJE  
gen_Sand4 M!xm1-,[  
group sand_o12 range x -9.55 0 y 0 1 z -16.26 -50 gp#bQ  
;------------------------------------------------------------------------------------------------- sYAG,r>h  
def gen_Sand5 "-0pz\a  
  command ?#VP)A  
gen zone brick      p0 ttx2,0,-50    p1 add 0,yc,0  p2 add tt1,0,0    p3 add 0,0,xh & z4%Z6Y  
                    p4 add tt1,yc,0  p5 add tt1,0,xh  p6 add 0,yc,xh  p7 add tt1,yc,xh &   A^T~@AO  
                    size 1 50 20 ;ratio 1 1.1 1 "<cB73tY  
G/LXUhuif  
        end_command 'U|MM;(  
end `@u+u0  
gen_Sand5 z/eU^2V  
group sand_o12 range x -78.45 -9.55 y 0 1 z -16.26 -50                        ;-33.74 0 7|h3.  
;------------------------------------------------------------------------------------------------- oW yN:Qh  
call Flacmodel_haitang.dat C }h<ldlY  
;------------------------------------------------------------------------------------------------- X]+(c_i:hC  
plot add axes Ux%\Y.PPI  
;-------------------------------------------------------------------------------------- ;JPbBwm  
;二、生成接触面的命令 ,p[\fT($]  
;-------------------------------------------------------------------------------------- T!=20!I  
interface 1 face range x -9.55 10.05 y 0 1 z -16.27 -16.25 )$gsU@H -  
interface 2 face range x -9.55 10.05 y 0 1 z -3.11 -3.09 C$AIP\j- )  
interface 3 face range x -9.56 -9.54 y 0 1 z -5 -3.1 @w#gRQCl  
interface 4 face range x -9.56 -9.54 y 0 1 z -8  -5 }\`-G+i{W  
interface 5 face range x -9.56 -9.54 y 0 1 z -14.8 -8 ^|F Vc48{  
interface 6 face range x -9.56 -9.54 y 0 1 z -16.26 -14.8 %U=S6<lbj;  
6*(h9!_T1  
int  1  prop kn  1.570e9  ks  1.570e9  fric  10 coh 10000                   ;prop属性 kn法向刚度 ks剪切刚度 fric为摩擦角 coh内聚力 RAD4q"}k  
int  2  prop kn  1.570e9  ks  1.570e9  fric  15 coh 55000 ery?G-  
int  3  prop kn  1.570e9  ks  1.570e9  fric  26 coh 8000 `aO@N(  
int  4  prop kn  1.570e9  ks  1.570e9  fric  15 coh 55000 j &0fC!k  
int  5  prop kn  1.570e9  ks  1.570e9  fric  14 coh 24000 PtPGi^  
int  6  prop kn  1.570e9  ks  1.570e9  fric  17 coh 56000 Sob $j  
ini y add -2.0 rang group 1 .^)C:XiW  
;=============================================== zP{<0o  
;初始应力计算程序 5ykk11!p$  
;=============================================== ERX|cc  
;定义模型及性质 }(g+:]p-  
model  mohr    rang  group sand_o1                                           {P<BJ52=  
prop bulk 1.875e7   shear   0.967e7 coh 10000   fric 10 rang group sand_o1 Zs3]|bUR  
model  mohr    rang  group sand_o2                                           )bIK0h  
prop bulk 1.875e7   shear   0.967e7 coh 55000   fric 15 rang group sand_o2 :~0^ib<v;  
model  mohr    rang  group sand_o3                                           (Qh7bfd  
prop bulk 2.068e7   shear   0.793e7 coh 8000   fric 26 rang group sand_o3 $3]E8t  
model  mohr    rang  group sand_o4                                           )/+eL RN5G  
prop bulk 1.875e7   shear   0.967e7 coh 55000   fric 15 rang group sand_o4 #8Id:56  
model  mohr    rang  group sand_o5                                           FB =  
prop bulk 1.396e7   shear   0.758e7 coh 24000   fric 14 rang group sand_o5 :,}:c%-^"  
model  mohr    rang  group sand_o6                                           [! dnm1   
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o6 >QZt)<[  
model  mohr    rang  group sand_o7                                           d_W nK{  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o7 mtWx ?x  
model  mohr    rang  group sand_o8                                           f_5R!;  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o8 Ng_rb KXC#  
model  mohr    rang  group sand_o9                                           Q,,fDBN  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o9 Kz3h]/A.  
model  mohr    rang  group sand_o10                                           C5sN[  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o10 53u.p c  
model  mohr    rang  group sand_o11                                           xeP;"J}  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o11     M#lVPXS  
model  mohr    rang  group sand_o12                                           jK C qH$  
prop bulk 2.178e7   shear   1.183e7 coh 56000   fric 17  rang group sand_o12 X!M fJ^)q  
\")YKN=W  
mo elas                 rang group 1 c@o/Cv  
;prop bulk 20.8e9  shear 15.58e9  rang group 1   ;正确 K}5 $;W#  
prop bulk 22.8e9  shear 17.08e9  rang group 1   ;只把柱子等效 c9g\7L,Z  
ini dens 2500 rang group 1 r/q1&*T  
ini dens 1900 rang group sand_o1     ;18000 9J;H.:WH  
ini dens 1990 rang group sand_o2 1TxhEXB  
ini dens 1910 rang group sand_o3 $OFFH[_z  
ini dens 1990 rang group sand_o4 %yuIXOJ  
ini dens 1970 rang group sand_o5 "$:y03V  
ini dens 2010 rang group sand_o6 nxh/&%  
ini dens 2010 rang group sand_o7 &A5[C{x  
ini dens 2010 rang group sand_o8 J*lYH]s  
ini dens 2010 rang group sand_o9 #u=O 5%.  
ini dens 2010 rang group sand_o10 Dgq[g_+l  
ini dens 2010 rang group sand_o11 | C2.Zay  
ini dens 2010 rang group sand_o12 u, SX`6%  
;边界条件 P2:Q+j:PX  
fix x y z rang z -50.1 -49.9 p_40V%y^  
fix x rang x  78.85  79.05 vVN[bD<  
fix x rang x -78.55 -78.35 k2tX$\E  
fix y rang y -.1  .1 OB  i!fLa  
fix y rang y 0.9  1.1 DwrCysIK  
;应力条件 rgZ rE;*;  
ini szz 0  grad 0 0 20e3   rang z -50 0     ;加权平均重度为1993 K?OX  
ini sxx 0  grad 0 0 10e3   rang z -50 0 \FY De  
ini syy 0  grad 0 0 10e3   rang z -50 0 fi4/@tV?$L  
;--------------------------------------------------------------------------------- ]zMBZs  
hist id 1 unbal "$"mWF-  
plot hist 1 <ZvPtW  
; !RUo:b+  
set grav 0 0 -10 \utH*;J|x  
solve ratio 1.0e-4 BiLreZ~"  
save station1.sav
回复
举报
分享到 淘江湖 新浪 QQ微博 QQ空间 开心 人人 豆瓣 网易微博 百度 鲜果 白社会 飞信
离线edelmann
总工

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

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

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

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

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