最近在开始用flac 5.0做动三轴实验，下面是完整的代码，有几个问题请教各位：
------------------------------------------------------------------------------
config gw dyn
grid 5,10
gen 0 0 0 0.10 0.05 0.10 0.05 0
model e
;
prop dens=1600 bulk=2E7 shear=1E7 cohesion=0.0
prop poros 0.5 fric=34.0 dilation=0.0 tension=0.0
;
fix x y j 1
def sigma3
sigma3=-50e3
end
;
apply nstress sigma3 i=1
apply nstress sigma3 i=6
apply nstress sigma3 j=11
;
ini sxx -50000
ini syy -50000
;ini szz -50000
set gravity=9.81
set dyn=off
set flow=off
history 999 unbalanced
solve
save state1.sav
;equilibrium
;
;dynamic analysis
rest state1.sav
model finn
;
set dyn=on
set flow=on
set large
;
water dens 1000 bulk 2e9 tens 1e10
;prop dens=2000 bulk=3E7 shear=2E7 cohesion=0.0
;prop fric=35.0 dilation=0.0 tension=0.0
prop perm 3e-8 poros 0.4
prop ff_latency=50
;
ini xvel=0 yvel=0 xdisp=0 ydisp=0
;
def sin_wave
while_stepping
sin_wave=ampl*sin(2*pi*freq*dytime)
end
;
def zs_load
zss_load=sigma3-sin_wave
command
apply nstress zss_load j=11
endcommand
end
;
ini dy_damp hyst default -3.325 0.823
;
;parameters for Byrne formular
prop ff_switch=1
def _setCoeff_Byrne
ff_c1_=8.7*exp(-1.25*ln(n_60))
ff_c2_=0.4/ff_c1_
ff_c3_ = 0.0000
end
;
set ncwrite 10
;set dydt=1e-5
set freq=1.0 ampl=25e3
set n_60=7
_setCoeff_Byrne
prop ff_c1=ff_c1_
prop ff_c2=ff_c2_
prop ff_c3=ff_c3_
;
;monitor some variables
his reset
his dytime
his pp i 3 j 5
his nstep 20
set step 1000000
zs_load
solve dytime 50
-------------------------------------------------------------
总体错误，无论怎么调参数，都是bad geometry错误。问题：
1）第一步重力及等压固结平衡时，需要初始化z方向的应力么？因为是2d，感觉不用管z方向的，但是通过查看手册中的振动台试验，发现那个例子中在初始化应力的时候初始化了z方向，不知道究竟是否需要初始化z方向？
通过学习，发现2d里面也要设置z方向，所以这点现在确定了。
2）在第一步静力平衡时已经设置了材料参数，如果材料参数不变，第二步动力荷载的时候还需要设置材料参数么？在我的尝试中，如果去掉第二步中的材料参数设置，程序在第一个网格处马上就会报错，如果加上材料参数设置，程序不会在第一个网格处报错，但是还是会出现bad geometry错误，请问这是为什么？如果第一步是等压固结，第二步中的材料参数需要适当改变么？
3）在自定义sin_wave动力荷载的时候，需要whilestepping么？看到例子中有的地方用了，有的地方没用，还没有深刻理解这个命令。
这个命令现在也明白了，所以也不是个问题了。
4）如果总是bad geometry错误的话，大概调试方向是什么？
多谢各位大侠！！！！
