/clear
!/view,1,1,1,1
!/ang,1
/prep7
et,1,plane42 !用来划分平面模型
et,2,SOLID45
!mp,ex,1,2.5e10 !混凝土材料弹性模量
mp,ex,1,2.5e11
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,2,0,0
K,3,2,0.5,0
K,4,1.25,0.5,0
K,5,1.25,3,0
K,6,1.25,3.25,0
K,7,0.75,3.25,0
K,8,0.75,3,0
K,9,0.75,0.5,0
K,10,0,0.5,0
K,11,0,-3.0,0
K,12,4.0,-3.0,0
K,13,4,0,0
K,14,4,3,0
K,15,0,3,0

L,1,2
L,2,3
L,3,4
L,4,9
L,9,10
L,10,1
AL,1,2,3,4,5,6 !生成面1


L,4,5
L,5,8
L,8,9
al,7,8,9,4 !生成面2


L,5,6
L,6,7
L,7,8
al,10,11,12,8 !生成面3


L,1,11
L,11,12
L,12,13
L,13,2
AL,13,14,15,16,1 !生成面4

L,13,14
L,14,5
AL,17,18,7,3,2,16 !生成面5

L,8,15
L,15,10
AL,19,20,5,9 !生成面6


/pnum,area,1
/pnum,line,1
aplot


lsel,s,,,1,6,1 !网格划分
lesize,all,0.25
amesh,1

lsel,s,,,7,9,1
lesize,all,0.25
amesh,2


lsel,s,,,10,12,1
lesize,all,0.25
amesh,3


lsel,s,,,13,16,1
lesize,all,0.25
amesh,4


lsel,s,,,17,18,1
lesize,all,0.25
amesh,5

lsel,s,,,19,20,1
lesize,all,0.25
amesh,6


EXTOPT,ESIZE,5,0, !拉伸的份数为5份
EXTOPT,ACLEAR,1 !清除原始的面网格
type,2 !拉伸后体的单元类型
mat,1 !拉伸后体的材料属性为混凝土材料
VEXT,1,,,,,-1.25 !设置拉伸长度为-1.25并拉伸成体

EXTOPT,ESIZE,1,0, !拉伸的份数为1份
EXTOPT,ACLEAR,1 !清除原始的面网格
type,2 !拉伸后体的单元类型
mat,1 !拉伸后体的材料属性为混凝土材料
VEXT,2,3,1,,,-0.25 !设置拉伸长度为-0.25并拉伸成体

EXTOPT,ESIZE,10,0, !以下是拉伸生成土层的模型
EXTOPT,ACLEAR,1
mat,2
VEXT,4,6,1,,,-2.5

EXTOPT,ESIZE,5,0,
vext,7,,,,,-1.25
mat,2

EXTOPT,ESIZE,9,0,
vext,14,,,,,-2.25
mat,2


allsel
nummrg,elem
nummrg,node
nummrg,kp
numcmp,all
/view,1,1,1,1
/ang,1
/PNUM,MAT,1
eplot
save
finish



/solu
/NERR,99999999,0
antype,static
OUTPR,BASIC,ALL !控制结果的输出
OUTRES,ALL,ALL !
NLGEOM,ON !启动大应变选项

asel,s,loc,x,0
aplot
da,all,symm
asel,s,loc,x,4
da,all,ux
aplot
asel,s,loc,z,-2.5
da,all,UZ
asel,s,loc,z,0
aplot
da,all,symm
asel,s,loc,y,-3.0
da,all,all
allsel


ACEL,0,9.8,0, !设置重力加速度
aSel,S,LOC,y,3.25
aplot
SFA,all,1,pres,-500
AUTOTS,ON !自动时间步长跟踪打开
NSUBST,10,,,1
TIME,1
solve
FINISH


/post1
set,last
plnsol,s,y,0,1 !竖向应力
plnsol,s,eqv,0,1 !等效应力
plnsol,eppl,EQV,0,1 !塑性区