某一剛性基礎，上有方柱，假設在方柱上部有一分布載荷，求基礎及周圍土的受力分布情況。整個模型都埋在土下。
　　基本參數如下：
　　（１）方柱及剛性基礎均為混凝土材料：彈模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