桩基础模型

finish
/clear,start

/prep7
ET,1,PLANE42
ET,2,SOLID45
MP,EX,1,1.0e9 ! 桩体材料属性
MP,PRXY,1,0.167
mp,dens,1,2437

MP,EX,2,1e8 ! 土体材料属性
MP,PRXY,2,0.1
mp,dens,2,2000
tb,dp,2
tbdata,1,12000,27,27

tb,creep,2
tbdata,1,1/1.55e6,1,0,0,10.88/600,2

!几何模型建立
PCIRC,0.035,0,0,90,
PCIRC,0.035,0.135,0,90,
LESIZE,2,,,8
LESIZE,3,,,8
LESIZE,1,,,8
mshape,0
mshkey,1
AMESH,1
LESIZE,5,,,20
LESIZE,7,,,20
LESIZE,4,,,8
LESIZE,6,,,8
mshape,0
mshkey,1
AMESH,2
ALLSEL
WPSTYLE
WPOFF,0,0,-0.25
PCIRC,0.135,0.035,0,90,
PCIRC,0,0.035,0,90,
ASEL,S,,,3,4,1
APLOT
LESIZE,14,,,8
LESIZE,13,,,8
LESIZE,12,,,8
mshape,0
mshkey,1
AMESH,4
LESIZE,11,,,20
LESIZE,9,,,20
LESIZE,10,,,8
LESIZE,8,,,8
mshape,0
mshkey,1
AMESH,3
ALLSEL
!进行拖拉
EXTOPT,ESIZE,20,1 !拖拉桩(地面以上）
EXTOPT,ACLEAR,0
TYPE,2
MAT,1
VEXT,1,,,0,0,0.1,,,,
EXTOPT,ESIZE,50,1 !拖拉桩(地面以下）
EXTOPT,ACLEAR,1
TYPE,2
MAT,1
VEXT,1,,,0,0,-0.25,,,,
EXTOPT,ESIZE,50,1 !拖拉土体
EXTOPT,ACLEAR,1
TYPE,2
MAT,2
VEXT,2,,,0,0,-0.25,,,,
EXTOPT,ESIZE,20,1 !拖拉土(端部1外）
EXTOPT,ACLEAR,1
TYPE,2
MAT,2
VEXT,3,,,0,0,-0.1,,,,
EXTOPT,ESIZE,20,1 !拖拉土(端部2内）
EXTOPT,ACLEAR,1
TYPE,2
MAT,2
VEXT,4,,,0,0,-0.1,,,,
ALLSEL
VSEL,S,,,3,5,1 !合并体3-5的重复单元及节点
VPLOT
ESLV,S !选择当前有效体上的单元
NSLV,S,1 !选择当前有单元上的节点
EPLOT
NUMMRG,ALL
NUMCMP,ALL
ALLSEL
VSEL,S,,,1,2,1 !合并体1-2的重复单元及节点
VPLOT
ESLV,S !选择当前有效体上的单元
NSLV,S,1 !选择当前有单元上的节点
EPLOT
NUMMRG,ALL
NUMCMP,ALL
ALLSEL
/PNUM,AREA,1 !显示实体单元上的面编号及颜色
/VIEW,1,1,1,1
/ANG,1
/PNUM,MAT,1,2
EPLOT
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ET,3,170 !刚性目标面的单元类型
ET,4,173 !柔性接触面的单元类型
r,10,,,0.05, !修改刚度
REAL,10
MP,MU,2,0.4 !库仑摩擦系数为0.25
KEYOPT,4,9,0
KEYOPT,4,12,1
VSEL,S,,,2
VPLOT
ASEL,S,,,9,10,1
APLOT
CM,TARGET,AREA
TYPE,3
NSLA,S,1
ESURF,ALL
esel,s,type,,3
eplot
/psymb,esys,1
allsel
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
VSEL,S,,,3,5,1
VPLOT
ASEL,S,,,15
ASEL,a,,,4
APLOT
CM,contact,AREA
TYPE,4
NSLA,S,1
ESURF,ALL
ESEL,S,TYPE,,4
EPLOT
/PSYMB,ESYS,1
allsel
GPLOT
SAVE
FINISH
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!


/SOLU !求解
rate,off
time,1.0e-6
!施加约束

nSEL,S,LOC,Z,-0.35
D,ALL,ALL
ALLSEL
ASEL,S,,,11
ASEL,a,,,14
ASEL,a,,,19
ASEL,a,,,23
ASEL,a,,,7
APLOT
DA,ALL,SYMM
ALLSEL
ASEL,S,,,12
ASEL,a,,,16
ASEL,a,,,21
ASEL,a,,,24
ASEL,a,,,8
APLOT
DA,ALL,SYMM
ALLSEL
ASEL,S,,,13
ASEL,A,,,18
APLOT
DA,ALL,ux
DA,ALL,uy
!nsla,s,0
!d,all,all



antype,static
deltim,0.1,0.05,0.2
nlgeom,on
autots,on
pred,on
lnsrc

h,on
nropt,full
cnvtol,f,,0.02,2,0.5
outres,all,9

!!!!!!!!!!!!!!!!!!!!!!!
!施加重力加速度及桩顶分布荷载
ACEL,0,0,9.8
nSEL,S,LOC,Z,0.1
esln,s
eplot
SFe,ALL,6,PRES,,10400
/PSF,PRES,NORM,2,0,1
!设置分析选项并求解

allsel
solve
finish

/solu
rate,on
bfunif,temp,40

kbc,0
nsub,10
allsel

*do,t,1,60,1
time,t
solve
*enddo