我现在做的东西里面要求地震作用下的边坡安全系数。我准备在动力计算完之后在进行安全系数求解。
下面是我的命令流，请各位前辈看一下问题出在什么地方？
new
config dyn
set dyn off
gen zone brick p0 0 0 0 p1 15 0 0 p2 0 0.5 0 p3 0 0 6 size 20 1 6
gen zone brick p0 15 0 0 p1 48 0 0 p2 15 0.5 0 p3 15 0 6 size 40 1 6 ratio 1.03 1 1
gen zone brick p0 15 0 6 p1 48 0 6 p2 15 0.5 6 p3 33 0 15 p4 48 0.5 6 p5 33 0.5 15 p6 48 0 15 p7 48 0.5 15 size 40 1 9 ratio 1.03 1 1
fix x y z range z -0.1 0.1
fix x range x 47.9 48.1
fix x range x -0.1 0.1
fix y
model elas
prop density 1905 bulk 1e9 shear 8e5
set gravity 0 0 -10.0
hist unbal
solve
set dyn on
free x y z ran z -.1 .1 ;去除静力边界
apply nquiet squiet dquiet ran z -.1 .1 ;在底部施加静态边界
;定义动荷载（以水平加速度形式）
table 1 read acc-16s.txt
apply xaccel 10 hist table 1 ran z -0.1 0.1 ;加动力荷载
apply ff;施加自由场边界
;给自由网格分组，方便显示主体网格，进行后处理
group ff_corner
group ff_side range x 0 48
group ff_side range y 0 0.5
group main_grid range x 0 48 y 0 0.5
set dyn time=0 ;设置动力计算从0s开始
hist reset
hist unbal
hist dytime
set dyn damp local .157 ;局部阻尼
solve age 9
plot cont ssi range gr main_grid
plot cont disp range gr main_grid
set dyn off
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
del xacc
mod null range x=-1,0
mod null range x=48,49
mod null range y=-1,0
mod null range y=0.5,1.5
fix x y z range z -0.1 0.1
fix x range x 47.9 48.1
fix x range x -0.1 0.1
fix y
model mohr
prop density 1905 bulk 1.9608e7 shear 7.5188e6 coh 22830 fric 10.27 dila 10.27
solve fos file slope3dfos1.sav associated
