求助 [复制链接]

new 8SKDL[rN  
gen zone brick size 1 1 2 [x)BQX'  
model mohr !TG"A
W  
prop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0 o.j;dsZ  
fix z ran z 0 0e/~H^,SQ  
fix x ran x 0 ExnszFX*  
fix x ran x 1 Kk??}  
fix y ran y 0 eL-92]]e  
fix y ran y 1 NGb`f-:jw  
ini dens 2000 [vIO  
ini szz -40e3 grad 0 0 20e3 ran z 0 2 Cf%)W:Q9  
ini syy -20e3 grad 0 0 10e3 ran z 0 2 Y3-15:-  
ini sxx -20e3 grad 0 0 10e3 ran z 0 2 Zl* HT%-5  
set grav 0 0 -10  a_Xh(d$  
solve {~d4;ht1Y  
plo con sz $v>- @  
7G
PBn}{W  
既然有dens、grav也就有了重力，干嘛还要初始化szz、syy、sxx ，而且这些应力的值是从哪里得出来的？ \W$
>EH  
请大家指点迷津 sb8SG_c.  

不初始化也可以的,或者在模型很不规则时,根本没有初始化的必要 zJtYy4jI)  
&o)j@5Y?  
N)&3(A@  
初始化的目的,在于加快系统的收敛速度, t0hg!_$bq  
!9V_U  
Soq 'B?>  
取值来源, 就是假定 szz 等于重力梯度 , sxx syy 符合 侧压力系数 G?YKm1:w  
<51(q_f  
%hOe `2#$  
但是当模型很不规则时,比如高度沿x方向一直在变, 这么处理并无法加快收敛速度 O(wt[AE
A  
4.R >mN[  
-_A0<A.  
lz 可以搜索我的帖子, 曾玩过这个, 在论坛发了两贴.