求助 [复制链接]

new <*Kh=v  
gen zone brick size 1 1 2 Fxth>O`$  
model mohr 6`baQ!xc.  
prop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0 v&b.Q:h*'  
fix z ran z 0 hi D7tb=g~  
fix x ran x 0 jvc?hUcLKT  
fix x ran x 1 UtBlP+bE?y  
fix y ran y 0 0ij~
e<  
fix y ran y 1 3 [)s;e  
ini dens 2000 P[GX}~_k  
ini szz -40e3 grad 0 0 20e3 ran z 0 2 iVi3 :7*  
ini syy -20e3 grad 0 0 10e3 ran z 0 2 0Q- Mxcj  
ini sxx -20e3 grad 0 0 10e3 ran z 0 2 kbp( a+5  
set grav 0 0 -10 Xm=^\K3  
solve 2]aZe4H.  
plo con sz %X,B-h^  
~{,vg4L  
既然有dens、grav也就有了重力，干嘛还要初始化szz、syy、sxx ，而且这些应力的值是从哪里得出来的？ $:BK{,\  
请大家指点迷津 H@2"ove-uC  

不初始化也可以的,或者在模型很不规则时,根本没有初始化的必要 <{NYD.  
`zzKD2y  
@"{'j  
初始化的目的,在于加快系统的收敛速度, NrJ_6sjF0g  
*L;pcg8{  
1ntkM?  
取值来源, 就是假定 szz 等于重力梯度 , sxx syy 符合 侧压力系数 YD='M.n\  
!V]MLA`  
k$-~_^4m  
但是当模型很不规则时,比如高度沿x方向一直在变, 这么处理并无法加快收敛速度 $D/bU lFx  
wsJ%* eYf  
-q&7J' N  
lz 可以搜索我的帖子, 曾玩过这个, 在论坛发了两贴.