我这个路基初始应力程序对吗？

ketaokeqi

我做的东西不难，我是新手，所以不会，希望大家帮帮忙！
  我做的模型是一个路基，我取了一个截面，土质是粉土，地下水位在路基的上表面，我只有“没有地下水时土质的参数”，不知道能不能直接带入我这个有地下水的模型？我现在只想计算初始应力，我编的程序不知道对不对，最重要的是参数问题，我快愁死了！
new                                       ；建立模型
config fluid
set fluid off
gen zone brick size    40  10   30     &         
                    p0 0   0    0      &  
                    p1 18  0    0      &
                    p2 0   2.25 0      &
                    p3 0   0    15   group base
gen zone brick size    20  10   30  ratio 1.1 &
                    p0 18  0    0      &
                    p1 60  0    0      &
                    p2 18  2.25 0      &
                    p3 18  0    15   group base1   
     
gen zone reflect dip 90 dd 90 origin 0 0 0  
mod m range group base                ； 设置模型
mod m range group base1
mod null range group ludi
pro bulk=24.664e6
shear=11.383e6
co=0.15e5
fri=33.1
ten 0 range group base
pro bulk=24.664e6
shear=11.383e6
co=0.15e5
fri=33.1
ten 0 range group base1
：问题1.红色的数值我直接取的没有地下水时的参数，不知道对不对
      
ini dens 2000 ran z 0 15  
；问题2.  2000也是我自己想的   
ini pp 1.5e5 grad 0 0 -10e3 ran z 0 15   
；问题3.我如果是按9.8作为重力加速度是不是要改成9.8？前面就变成9.8*15e3
model fl_iso                                                                                                                               prop por 0.4 perm 1e-10
ini fden 1000 ften -1e10
ini sat 1 ran z  0 15
；问题4.上面这句需要吗？是不是多余呀？
fix x y z range z -0.1 0.1
fix y range y -0.1 0.1
fix y range y 2.15 2.35
fix x range x 59.9 60.1
fix x range x -59.9 -60.1
set gravity=(0.0, 0.0, -10)
solve
save 初始应力.sav
我怎么觉得全是错的啊？救命呀！

ketaokeqi

高手们快来帮帮我吧

ketaokeqi

我都等了一天了，怎么没人来啊？