跪求高手指点迷津：深埋工程初始应力平衡后为什么sxx应力图总是不对？|FLAC/FLAC3D - 领先的岩土技术社区，岩土领域的媒体、社区与应用平台！

离线zwl99jy

技术员

发帖: 48

土币: 84

威望: 1

原创币: 0

只看楼主倒序阅读使用道具楼主发表于: 2012-10-21

1017.zip (1799 K) 下载次数:0 这是ansys建的模型，直接导入了flac~

我的初始平衡命令是这样的

def cal_weiyan
  wk=we/(3*(1-2*wp)),体积模量
  wg=we/(2*(1+wp));剪切模量
end
Set we=6e9 wp=0.35 ;弹性模量
cal_weiyan
model mohr
prop bulk=wk shear=wg fric 35 c 1e10 ten 1e10 range z -150 0;
ini dens 2600
fix z range z -149.9 -150.1
app sxx -2.08e7 grad 0 0 52e3 range x -0.1 0.1
app sxx -2.08e7 grad 0 0 52e3 range x 169.9 170.1
app syy -2.08e7 grad 0 0 52e3 range y -0.1 0.1
app syy -2.08e7 grad 0 0 52e3 range y 99.9 100.1
app szz -1.04e7 grad 0 0 26e3 range z 0 -150
set grav 0 0 -10
set mech ratio 1.0e-5
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
hist unbal                        ;记录计算过程中的体系最大不平衡力
plot cont zdis
plot hist 1                       ;绘制体系最大不平衡力历时曲线
solve
;二阶平衡
model mohr
prop bulk=wk shear=wg fric 35 c 6e5  range z -150 0;
solve
save ph-1021.sav

算出的SXX结果如上图，已经找好长时间的错误了还是没找出，很着急要出结果啊，跪求各位大神们的帮忙啊~~~