两种计算同一边坡安全系数命令流的区别

阳光明媚

书上第一种是计算边坡安全系数的过程是：建立模型-设置边界条件-生成初始地应力-赋真实材料参数-计算安全系数，命令流如下：
new
;建立模型网格
gen zone brick &
p0 0 0 0 p1 2 0 0 p2 0 0.5 0 p3 0 0 3 size 3 1 3
gen zone brick &
p0 2 0 0 p1 20 0 0 p2 2 0.5 0 p3 2 0 3 &
size 17 1 3 ratio 1.03 1 1
gen zone brick &
p0 2 0 3 p1 20 0 3 p2 2 0.5 3 p3 12 0 13 &
p4 20 0.5 3 p5 12 0.5 13 p6 20 0 13 &
p7 20 0.5 13 size 17 1 17 ratio 1.03 1 1

;设置边界条件
fix x y z ran z -0.1 0.1
fix x ran x 19.9 20.1
fix x ran x -0.1 0.1
fix y

;生成初始地应力
m e
prop den 2000 b 3e9 s 1e9
set grav 0 0 -10
plot con zdis ou on
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0

;安全系数求解
m m
prop den 2000 bulk 1e8 shear 3e7 &  ；材料真实参数
coh 12380 ten 1e6 fric 20 dila 20
solve fos file slope3dfosl.sav associated
plot fos cont ssi ou on vel red

第二种是计算边坡安全系数的过程是：建立模型-设置边界条件-赋真实材料参数-静力平衡计算-计算安全系数，命令流如下：
new
;建立模型网格
gen zone brick &
p0 0 0 0 p1 2 0 0 p2 0 0.5 0 p3 0 0 3 size 3 1 3
gen zone brick &
p0 2 0 0 p1 20 0 0 p2 2 0.5 0 p3 2 0 3 &
size 17 1 3 ratio 1.03 1 1
gen zone brick &
p0 2 0 3 p1 20 0 3 p2 2 0.5 3 p3 12 0 13 &
p4 20 0.5 3 p5 12 0.5 13 p6 20 0 13 &
p7 20 0.5 13 size 17 1 17 ratio 1.03 1 1

;设置边界条件
fix x y z ran z -0.1 0.1
fix x ran x 19.9 20.1
fix x ran x -0.1 0.1
fix y

;静力平衡计算

m m
prop den 2000 bulk 1e8 shear 3e7 &  ;材料真实参数
coh 12380 ten 1e6 fric 20 dila 20

;安全系数求解
solve fos file slope3dfosl.sav associated
plot fos cont ssi ou on vel red
请问这两种方式计算出来的安全系数是一样的吗？先生成初始地应力有什么好处？不能直接赋真实材料计算吗？

江心野鹤

生成初始应力时常赋大值，以免屈服