关于软土地基沉降计算问题！

shulichen

我做的是海域的堤坝沉降计算，模型很简单，就只是堤坝与软土地基。海域的软土基本全是淤泥，并且很厚，计算到堤底以下100m。堤坝为堆石坝，模型如下图所示。
        我的问题是：1.淤泥应取什么模型？我取的弹性模型，计算结果达60多cm，我知道是不对的，但是用摩尔库仑的话，堤坝自身沉降总是很大，总沉降也很大（原意是不考虑大坝自身的沉降的）；2.要考虑大变形吗?我计算了很久了，才想起来要考虑大变形的问题
       下面是我的命令流，请大侠指点阿！着急

;---------------------------------------------------
group dibas range group 2 z 0 -8.5
group dibax range group 2 z -8.5 -122.5
group soil range group dibas not group dibax not
;---------------------------------------------------
config fluid
;---------------------------------------------------
;-------------------------------------------------------------------------------------------------
def constants
                 
    c_E_db=2.5e8    ;堤坝的变形模量
    c_poss1=0.35       ;柏松比
    c_poss_db=0.25

    c_fric1=15          ;内摩擦角
    c_fricdb=45
                    
    c_coh_db=0.0        ;堤坝的粘聚力
                                    
    c_bulk_db=c_E_db/(3*(1.0-2*c_poss_db))  ;堤坝的体积模量
                                          
    c_shear_db=c_E_db/(2*(1.0+c_poss_db))   ;堤坝的剪贴模量

    c_density1=1.7e3                 ;密度
    c_density2=2.0e3

    c_perm1=1.22e-11                ;渗透系数
    c_perm2=2.25e-8
   
    c_poro1=0.5                     ;孔隙率
    c_poro2=0.3
   

    c_w_density=1e3                 ;水的密度
                           
end
constants
;-------------------------------------------------------------
def constants1
    zz=-8.5-6.12
    c_coh1=[3.7305+0.535*(-zz)]*1000
   
    c_E1=3.21*1000000                      ;第一层土的变形模量
    c_bulk1=c_E1/(3*(1.0-2*0.35))
    c_shear1=c_E1/(2*(1.0+0.35))
end
def constants2
    zz=-8.5-6.12-12.24
    c_coh2=[3.7305+0.535*(-zz)]*1000
   
    c_E2=3.43*1000000                     ；第二层土的变形模量
    c_bulk2=c_E2/(3*(1.0-2*0.35))
    c_shear2=c_E2/(2*(1.0+0.35))
end
def constants3
    zz=-8.5-6.12-12.24*2
    c_coh3=[3.7305+0.535*(-zz)]*1000
  
    c_E3=5.61*1000000
    c_bulk3=c_E3/(3*(1.0-2*0.35))
    c_shear3=c_E3/(2*(1.0+0.35))
end
def constants4
    zz=-8.5-6.12-12.24*3
    c_coh4=[3.7305+0.535*(-zz)]*1000
   
    c_E4=7.76*1000000
    c_bulk4=c_E4/(3*(1.0-2*0.35))
    c_shear4=c_E4/(2*(1.0+0.35))
end
def constants5
    zz=-8.5-6.12-12.24*4
    c_coh5=[3.7305+0.535*(-zz)]*1000
   
    c_E5=9.87*1000000
    c_bulk5=c_E5/(3*(1.0-2*0.35))
    c_shear5=c_E5/(2*(1.0+0.35))
end
def constants6
    zz=-8.5-6.12-12.24*5
    c_coh6=[3.7305+0.535*(-zz)]*1000
   
    c_E6=11.97*1000000
    c_bulk6=c_E6/(3*(1.0-2*0.35))
    c_shear6=c_E6/(2*(1.0+0.35))
end
def constants7
    zz=-8.5-6.12-12.24*6
    c_coh7=[3.7305+0.535*(-zz)]*1000
   
    c_E7=14.06*1000000
    c_bulk7=c_E7/(3*(1.0-2*0.35))
    c_shear7=c_E7/(2*(1.0+0.35))
end
def constants8
    zz=-8.5-6.12-12.24*7
    c_coh8=[3.7305+0.535*(-zz)]*1000
   
    c_E8=16.13*1000000
    c_bulk8=c_E8/(3*(1.0-2*0.35))
    c_shear8=c_E8/(2*(1.0+0.35))
end
def constants9
    zz=-8.5-6.12-12.24*8
    c_coh9=[3.7305+0.535*(-zz)]*1000
   
    c_E9=18.19*1000000
    c_bulk9=c_E9/(3*(1.0-2*0.35))
    c_shear9=c_E9/(2*(1.0+0.35))
end
def constants10
    zz=-8.5-6.12-12.24*9
    c_coh10=[3.7305+0.535*(-zz)]*1000
   
    c_E10=20.24*1000000
    c_bulk10=c_E10/(3*(1.0-2*0.35))
    c_shear10=c_E10/(2*(1.0+0.35))
end
constants1
constants2
constants3
constants4
constants5
constants6
constants7
constants8
constants9
constants10
;---------------------------------------------------
model null range group dibas
;---------------------------------------------------
model ela  range z -8.5 -130.9
prop bulk c_bulk1 shear c_shear1   rang z  -8.5   -20.74
prop bulk c_bulk2 shear c_shear2   rang z -20.74  -32.98
prop bulk c_bulk3 shear c_shear3   rang z -32.98  -45.22
prop bulk c_bulk4 shear c_shear4   rang z -45.22  -57.46
prop bulk c_bulk5 shear c_shear5   rang z -57.46  -69.7
prop bulk c_bulk6 shear c_shear6   rang z -69.7   -81.94
prop bulk c_bulk7 shear c_shear7   rang z -81.94  -94.18
prop bulk c_bulk8 shear c_shear8   rang z -94.18  -106.42
prop bulk c_bulk9 shear c_shear9   rang z -106.42 -118.66
prop bulk c_bulk10 shear c_shear10   rang z -118.66 -130.9
ini density c_density1
set grav 0 0 -10
fix x range x 107.78 107.76
fix x range x -107.78 -107.76
fix y
fix x y z range z -130.7 -130.9
;---------------------------------------------------
model fl_iso
prop perm c_perm1 poro c_poro1
ini fdensity 1000

ini pp -8.5e4 grad 0 0 -1e4 range z -8.5 -130.9
fix pp range z -8.5
fix pp range x 107.78 107.76
fix pp range x -107.78 -107.76
fix pp range y 0.1 -0.1
fix pp range y 0.9 1.1
;--------------------------------------------------
;----------------ini balance-------------------------------------------------
hist n=5
hist unbal
hist gp zdis 0 0 0
plot create hist12
plot add hist 1 2
plot show
set fluid off
ini fmod 0
sol
sav 初始平衡624

;----------------------------------------------------------------------------
ini xdis 0 ydis 0 zdis 0
;----------------------------------------------------------------------------
model null range group dibax
;----------------------------------------------------------------------------
model ela   range group dibas any group dibax any
prop bulk c_bulk_db  shear c_shear_db      range group dibas any group dibax any
ini density c_density2   range group dibas any group dibax any
model fl_iso  range group dibas any group dibax any
ini fmod 2e9
ini pp 0 range group dibas
ini fdensity 1000              range group dibas any group dibax any
prop perm c_perm2 poro c_poro2 range group dibas any group dibax any

set fluid pcut on
set fluid implicit on
set fluid dt 10
set fluid on mech on
set mech force 10
sol age 5.9e7 ;1.8年

hello_zy

在海里的防波堤沉降1、2米是很正常的，你可以用分层总和法算算比较一下

xjiang01

plaxis maybe simple than flac.

shulichen

谢谢关注，这两种方法我都算过了，分层综合法计算为3.75m ,plaxis计算结果是17cm，用的是软土模型，计算过程也比较简单。
      但我还是想问一下，这个用flac计算要用大变形计算吗？软土用MC 行吗？至少用弹性模型模拟是不合理的吧？

shulichen

急啊，高手过来给个建议啥！

faye100

海堤的沉降ms本来就很大啊。。。
我以前用adina试过mc材料，收敛不了，也是算类似的。。。
求讨论啊
~~~ 4# shulichen lz

liuai

我现在做软土地区的路基沉降，实际最终沉降也要在1-2米吧，模型我现在用M-C但还没有算出来，我和你一样也在考虑是否考虑大变形问题将nlgeom on 了。至于模型，我同事说最好 用修正剑桥模型，但是这个模型的参数很难取，只能用别人的数据，所以最近一直在用M-c. 也有可能像前面朋友说的那样，M-c根本不收敛。不过那个朋友说其它模型结果不太正常没有都收敛不了。