求助 [复制链接]

new 952l1c!
  
gen zone brick size 1 1 2 /i)>
|U 4  
model mohr N;S1s0FN  
prop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0 s8.SEk|pB  
fix z ran z 0 `6xr:s  
fix x ran x 0 @hw
e  
fix x ran x 1 7m4*dBTr  
fix y ran y 0 Yfr4<;%  
fix y ran y 1 n?:%>Os$  
ini dens 2000 g[Q+DT  
ini szz -40e3 grad 0 0 20e3 ran z 0 2 .]YTS
 
ini syy -20e3 grad 0 0 10e3 ran z 0 2 7g|EqJ7  
ini sxx -20e3 grad 0 0 10e3 ran z 0 2 UvF5u(o  
set grav 0 0 -10 CB^.N>'  
solve YY4q99^K  
plo con sz i]15g@  
7S2"e[-x  
既然有dens、grav也就有了重力，干嘛还要初始化szz、syy、sxx ，而且这些应力的值是从哪里得出来的？ 5tl($j  
请大家指点迷津 K$KVm^`  

不初始化也可以的,或者在模型很不规则时,根本没有初始化的必要 U!(es0rX  
qw 03]a  
J0Gjo9L  
初始化的目的,在于加快系统的收敛速度, Zo,066'+[.  
Y?$  
S_TD o  
取值来源, 就是假定 szz 等于重力梯度 , sxx syy 符合 侧压力系数 '&nQ~=3  
fejC,H4I  
=*R6O,  
但是当模型很不规则时,比如高度沿x方向一直在变, 这么处理并无法加快收敛速度 O7p=N8V  
MdN0 Y@Ll  
p|?FA@ 3  
lz 可以搜索我的帖子, 曾玩过这个, 在论坛发了两贴.