刚性基础三维仿真分析(土体为DP材料)

某一刚性基础，上有方柱，假设在方柱上部有一分布载荷，求基础及周围土的受力分布情况。整个模型都埋在土下。
　　基本参数如下：
　　（１）方柱及刚性基础均为混凝土材料：弹模E＝2.5ｅ10，泊松比为0.16，密度为2700。
　　（２）土为DP材料：弹模E＝1.0ｅ7，泊松比为0.45，密度为2000.
粘聚力C＝10，内摩擦角30，膨胀角30
　　（３）外载：桩顶分布力Q=50,重力加速度为9.8
命令流程序为：
/cle
/prep7
/title,zhuang,h-method
et,1,plane42
et,2,solid45
mp,ex,1,2.5e10 !混凝土材料属性
mp,nuxy,1,0.16
mp,dens,1,2700
mp,ex,2,10e6 !土体材料属性
mp,nuxy,2,0.45
mp,dens,2,2000
tb,dp,2 ! 土体本构模型
tbdata,1,10,30,30
!建立关键点
k,1
k,2,1.5
k,3,1.5,0.4
k,4,1.1,0.4
k,5,1.1,0.8
k,6,0.7,0.8
k,7,0.7,1.2
k,8,0.3,1.2
k,9,0.3,2
k,10,0,2
k,11,4,2
k,12,4,0
k,13,4,-2
k,14,0,-2
k,15,0.3,3
k,16,0,3
!建立各条线，面元素
*do,i,1,9,1
l,i,i+1
*enddo
l,10,1
al,1,2,3,4,5,6,7,8,9,10
l,9,11
*do,i,11,13,1
l,i,i+1
*enddo
l,14,1
l,2,12
al,2,3,4,5,6,7,8,11,12,16
al,1,16,13,14,15
l,9,15
l,15,16
l,16,10
al,9,17,18,19
/pnum,area,1
/pnum,line,1
aplot
!划分网格
lsel,s,,,2,8,1
lesize,all,0.2
lsel,s,,,9
lesize,all,0.15
lsel,s,,,10
lesize,all,0.2
lsel,s,,,1
lesize,all,,,8
amesh,1
lsel,s,,,11
lesize,all,0.2
lsel,s,,,12
lesize,all,0.2
lsel,s,,,16
lesize,all,0.2
amesh,2
lsel,s,,,15
lesize,all,0.2
lsel,s,,,14
lesize,all,,,21
lsel,s,,,13
lesize,all,0.2
amesh,3
lsel,s,,,17,19,2
lesize,all,0.2
lsel,s,,,18
lesize,all,,,2
amesh,4
!生成体模型
extopt,esize,3,0
extopt,aclear,1
type,2
mat,1
vext,1,4,3,,,-0.6
extopt,esize,6,0
extopt,aclear,1
mat,2
vext,2,3,1,,,-1.2
extopt,esize,3,0
vext,5,,,,,-0.6
allsel
nummrg,all
numcmp,all
/view,1,1,1,1
/ang,1
/pnum,mat,1
eplot
save
fini
!进入求解器，施加边界条件，外部载荷
/solu
asel,s,loc,x,0
aplot
da,all,symm
asel,s,loc,x,4
da,all,ux
aplot
asel,s,loc,z,0
da,all,symm
asel,s,loc,z,-1.2
da,all,uz
asel,s,loc,y,-2
da,all,all
allsel
!施加载荷
acel,0,9.8,0
nsel,s,loc,y,3
esln,s
eplot
sfe,all,4,pres,,50
antype,static
time,1
nsub,100
outres,all,all
alls
gplot
solv
fini