ansys土石坝静力分析，结果不合规律，请大家群主给点建议吧

杜晓帆

土石坝静力分析，运用ansys自带的DP准则，单元生死模拟施工过程，在此过程中，对死单元的节点进行了约束。但对于同一个工程，分别进行了二维和三维模型进行了计算，情况如下：（1）二维模型计算结果：

                               
登录/注册后可看大图

                               
登录/注册后可看大图

个人认为二维计算的结果（沉降和水平位移）还是很符合实际情况的。但是运用相同的求解程序，三维计算结果如下：

                               
登录/注册后可看大图

                               
登录/注册后可看大图

三维计算结果存在的问题：
（1）土石坝整体沉降值很小，明显错误；
（2）土石坝中间部位的最大沉降范围也太大了，不知道是什么地方的错误；
（3）三维模型的水平位移相对二维模型的值也偏小不少。

小弟最近在做这方面的研究，卡了很长一段时间，希望大家给点意见。求解程序如下：
/solu
antype,static   !设定为静力求解
nsubst,10   
pred,on   !打开时间步长预测器
nropt,full   !设定牛顿-拉普森选项
nlgeom,on   !打开大位移效果
lnsrch,on   !打开线性搜索
outres,all,all   !输出所有项
acel,0,10,0

!=================================地基===============================
time,1
allsel
esel,s,mat,,1,3
ekill,all
esel,all
esel,s,live
nsle,s
nsel,inve
d,all,all
nsel,all
esel,all
solve
save,moxing60,db

!=================================BT-1===============================
time,2
allsel
cmsel,,BT-1
ealive,all
nsle,all
ddele,all,all
allsel
/replot
solve
save,moxing60,db

!=================================BT-2===============================
time,3
allsel
cmsel,,BT-2
ealive,all
nsle,all
ddele,all,all
allsel
/replot
solve
save,moxing60,db

!=================================BT-3===============================
time,4
allsel
cmsel,,BT-3
ealive,all
nsle,all
ddele,all,all
allsel
/replot
solve
save,moxing60,db

!=================================BT-4===============================
time,5
allsel
cmsel,,BT-4
ealive,all
nsle,all
ddele,all,all
allsel
/replot
solve
save,moxing60,db

!=================================BT-5===============================
time,6
allsel
cmsel,,BT-5
ealive,all
nsle,all
ddele,all,all
allsel
/replot
solve
save,moxing60,db

!==========================设定水压力（1190.0m）=======================
time,7
allsel
cmsel,,BT-6
ealive,all
nsle,all
ddele,all,all
allsel
/replot
solve
save,moxing60,db

!=================================BT-7===============================
time,8
allsel
cmsel,,BT-7
ealive,all
nsle,all
ddele,all,all
allsel
/replot
solve
save,moxing60,db

!=================================BT-8===============================
time,9
allsel
cmsel,,BT-8
ealive,all
nsle,all
ddele,all,all
allsel
/replot
solve
save,moxing60,db

!=================================BT-9===============================
time,10
allsel
cmsel,,BT-9
ealive,all
nsle,all
ddele,all,all
allsel
/replot
solve
save,moxing60,db

!==========================设定水压力（1241.0m）=====================
time,11
allsel
cmsel,,XQ
eplot   
cmsel,,SYSH-4
nplot
/PSF,PRES,,2,0,1               
sfgrad,pres,0,y,0,-9800   !设定水压力（1241.0m）               
sf,all,pres,SYSH-4
allsel
/replot
solve
/go
save,moxing60,db