求助 [复制链接]

new ]DcFySyv  
gen zone brick size 1 1 2 ;TYBx24vD'  
model mohr uFE)17E  
prop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0 6Oq7#3]  
fix z ran z 0 )e{aN+  
fix x ran x 0 Da|z"I x  
fix x ran x 1 )8AXm  
fix y ran y 0 I,tud
!p`  
fix y ran y 1 w:0E(z  
ini dens 2000 iTwm3V P  
ini szz -40e3 grad 0 0 20e3 ran z 0 2 ;pAK_>
 
ini syy -20e3 grad 0 0 10e3 ran z 0 2 V88p;K$+  
ini sxx -20e3 grad 0 0 10e3 ran z 0 2 eFgA 8kY)  
set grav 0 0 -10 7
dWS  
solve K0~rN.C!0  
plo con sz It(_v  
j%kncGS  
既然有dens、grav也就有了重力，干嘛还要初始化szz、syy、sxx ，而且这些应力的值是从哪里得出来的？ %EH)&k  
请大家指点迷津 h{Y",7]!  

不初始化也可以的,或者在模型很不规则时,根本没有初始化的必要 -).C  
=X}J6|>X  
vM={V$D&  
初始化的目的,在于加快系统的收敛速度, vx =&QavL  
2?C)&  
wYea\^co  
取值来源, 就是假定 szz 等于重力梯度 , sxx syy 符合 侧压力系数  mh%VrAq  
5x
i
EPh  
CIWO7bS
 
但是当模型很不规则时,比如高度沿x方向一直在变, 这么处理并无法加快收敛速度 }MySaL>  
w0. u\  
xT8?&Bx  
lz 可以搜索我的帖子, 曾玩过这个, 在论坛发了两贴.