命令流：隧道开挖后为什么最大沉降却在边界附近？

zdks_2000

已经自重平衡了

;--------------------------------前处理--------------------------------------
restore tunnel&cp_grave.sav

plot block group
plot add axes  red
plot show

model mohr

def derive2      ;计算体积模量,剪切模量
        b_mod1=e_mod1/(3.0*(1.0-2.0*p_ratio1))         
        s_mod1=e_mod1/(2.0*(1.0+p_ratio1))
                b_mod2=....
        s_mod2=        ....
             b_mod3=....
        s_mod3=        ....   
              b_mod4=....
        s_mod4=        ....

                 

end
set e_mod1=3.15e6 p_ratio1=0.3243 e_mod2=6.51e6 p_ratio2=0.2537 e_mod3=2.07e6 p_ratio3=0.3976  &
    e_mod4=4.55e6 p_ratio4=0.2857  &
  
   e_mod10=34.5e9   p_ratio10=0.25   
derive2



;1土体   
prop bulk  b_mod1  shear s_mod1   fric 17 coh 11e3 ten 30e3   range group 1soil

;2土体   
prop bulk  b_mod2  shear s_mod2  fric 27 coh 10e3 ten 30e3    range group 2soil

;3土体   
prop  bulk  b_mod3  shear s_mod3  fric 10 coh 10e3 ten 38e3    range group 3soil

;4土体  
prop  bulk  b_mod4  shear s_mod4  fric 20.5 coh 20e3 ten 72e3  range group 4soil


;assign boundary conditions 边界条件
fix x range x -0.1 0.1
fix x range x -75.1 -74.9
fix z range z -0.1 0.1
fix z range z -10.1 -9.9
fix x y z range y -57.5 -57.3

;1土体  
ini density 1770   range group 1soil

;2土体  
ini density 1840  range group 2soil


;3土体  
ini density 1680 range group 3soil


;4土体   
ini density 1810  range group 4soil


;开挖主隧道
model null range group 11soil


def derive3     ;计算体积模量,剪切模量
        
        b_mod10=e_mod10/(3.0*(1.0-2.0*p_ratio10))         
        s_mod10=e_mod10/(2.0*(1.0+p_ratio10))
      
        

end
set     e_mod10=34.5e9   p_ratio10=0.25        
derive3







; assign initial stress state  初始应力赋值
;水土分算，取浮容重

;水压力
apply syy -1.43e4 rang y -0.1 0.1  ;上覆水位14.3m
;1-2土侧压力系数0.5                                                              
ini syy -1.43e4 grad 0 0 1.77e4 range y -4.8 0      
ini sxx -0.715e4 grad 0 0 1.77e4 range y -4.8 0        
ini szz -0.7151e4 grad 0 0 1.77e4 range y -4.8 0        







;监测
hist unbal




;y位移等值线
plot con ydisp

solve

save tunnelandcp_digtunnel.sav

[ 本帖最后由 zdks_2000 于 2008-9-9 13:31 编辑 ]

zdks_2000

请高手指点

zdks_2000

自己顶起来。是不是网格问题呀

splend

有没有平衡地应力？

zdks_2000

第一步已经完成自重平衡了，望高手指点

[ 本帖最后由 zdks_2000 于 2008-9-9 09:51 编辑 ]

strive

估计是边界条件造成的！

建议仔细检查边界条件，这句fix x y z range y -57.5 -57.3, 取掉x z 再试试！

haiou983

fix x y range x -0.1 0.1        ；x起点面沿x轴，y轴方向不产生位移
fix x y range x -75.1 -74.9    ；x终点面沿x轴，y轴方向不产生位移
fix x y z range z -10.1 -9.9   ；z轴底面沿x轴，y轴z轴方向不产生位移
fix y range y (前边界)           ；在y的前端面沿y轴方向产生位移
fix y range y （后边界）    ；在y的后端面沿y轴方向产生位移

不知道对不对，我模拟隧道是这么定义的

zdks_2000

呵呵 谢谢各位跟帖，我的grav方向错了 ，因为坐标转换，grav没认真看。

zhaozhenlin

坐标怎么回事？竖向不是z方向？？？

bcoolior

应该是坐标轴定义有点问题

abblcn

估计是边界条件造成的！