求助 [复制链接]

new Ug(;\*yg  
gen zone brick size 1 1 2 Z./$}tVUG  
model mohr Kv(2x3("  
prop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0 FyleK+D?  
fix z ran z 0 VRden>vKN  
fix x ran x 0 DmzK* O{  
fix x ran x 1 sZ,xbfZby  
fix y ran y 0 &h^E_]P  
fix y ran y 1 v$
~1{}iI5  
ini dens 2000 xVYy`_|  
ini szz -40e3 grad 0 0 20e3 ran z 0 2 fNR2(8;}  
ini syy -20e3 grad 0 0 10e3 ran z 0 2 z OwKh>]  
ini sxx -20e3 grad 0 0 10e3 ran z 0 2 #o`y<1rN  
set grav 0 0 -10 C:z+8wt  
solve
ybk~m  
plo con sz pY@$N&+W  
Y$JGpeq8w  
既然有dens、grav也就有了重力，干嘛还要初始化szz、syy、sxx ，而且这些应力的值是从哪里得出来的？ Q8-;w{%  
请大家指点迷津 ws}cMX]*  

不初始化也可以的,或者在模型很不规则时,根本没有初始化的必要 5:R$xgc  
ov3FKMG?  
q8Jhs7fv  
初始化的目的,在于加快系统的收敛速度, E5;6ks)  
T'^ Do/  
DtEvt+h  
取值来源, 就是假定 szz 等于重力梯度 , sxx syy 符合 侧压力系数 6DkFIkS  
8NF93tqD6  
p]jkfsCjN  
但是当模型很不规则时,比如高度沿x方向一直在变, 这么处理并无法加快收敛速度 @n|Mr/PAj  
-G'U\EXT  
nj1T
X  
lz 可以搜索我的帖子, 曾玩过这个, 在论坛发了两贴.