初始应力例子(附命令流)

zs814811599

n
;***************************************************************
;导入模型
impgrid W8-43.flac3d
;***************************************************************
plot bl gr

plot add axes red
;通过旋转坐标

ini x add -157.35
ini z add  -10
ini y add  152.56
;通过旋转坐标
def xiugai_11  
p_gp=gp_head
loop while p_gp # null   
    w_temp1=gp_xpos(p_gp)   
   w_temp2=gp_ypos(p_gp)     
  w_temp3=gp_zpos(p_gp)     
   gp_xpos(p_gp)=w_temp1     
   gp_ypos(p_gp)=-w_temp3   
   gp_zpos(p_gp)=w_temp2        
   p_gp = gp_next(p_gp)
endloop xiugai_11
end
xiugai_11
fix  x y z ran z 0
fix x ran x 0
fix x ran x 690
fix y


model mohr ran group Jiyan
model mohr ran group Fgc
model null ran group Tl

prop bulk 7.8e6 shear 3.0e6 coh 10e10 tension 1e10 ran group Fgc
ini dens 1500 ran group Fgc
prop bulk 3.91e6 shear 1.5e6 coh 10e10 tension 1e10  ran group Jiyan
ini dens 1800 ran group Jiyan
set grav 0 0 -9.8
hist id=1 unbal
solve
prop bulk 7.8e6 shear 3.0e6 coh 10e3 fric 15 ran group Fgc
prop bulk 3.91e6 shear 1.5e6  coh 10e3 fric 20 ran group Jiyan

solve
save W8-4.sav

大家看看这个有事么问题呢？ 谢谢大家哈！

zs814811599

大家帮忙看看吧  很急的