/FILNAME,biye
/TITLE，guanzhuang
/UNITS,SI
/UNITS,SI !国际单位制
/prep7
ET,1,SOLID45
MP,MU,1,0.2
MP,EX,1,3.8e10
MP,NUXY,1,0.2
MP,DENS,1,2500
MP,EX,2,5.25e6
MP,NUXY,2,0.4
MP,DENS,2,1.93e3
MP,EX,3,2.26e6
MP,NUXY,3,0.4
MP,DENS,3,1.76e3
MP,EX,4,1e7
MP,NUXY,4,0.2
MP,DENS,4,2000
TB,CONCR,1
TBDATA,,0.7,0.9,2.6,-1
TB,DP,2
TB,DP,2,,,
TBMODIF,1,1,51000
TBMODIF,1,2,14
TBMODIF,1,3,0
TBDE,DP,3
TB,DP,3,,,
TBMODIF,1,1,5900
TBMODIF,1,2,2
TBMODIF,1,3,0
TBDE,DP,4
TB,DP,4,,,
TBMODIF,1,1,0
TBMODIF,1,2,28
TBMODIF,1,3,0
/VIEW,1,1,1,1
/ANG,1
/REP,FAST
K,1,1.5,0.866,0
K,2,-0.5,0.866,0
K,3,-1.5,-0.866,0
K,4,0.5,-0.866,0
A,1,2,3,4
VOFFST,1,3.5,,
VOFFST,2,18,,
VOFFST,7,0.5,,
VOFFST,12,5,, ! 土层4
CYLIND,0,0.3,0,3.5,0,360
CYLIND,0,0.3,3.5,21.5,0,360
CYLIND,0,0.3,21.5,22,0,360
VSBV,1,5 ! 土层1
VSBV,2,6 ! 土层2
VSBV,3,7 ! 土层3
CYLIND,0.23,0.3,0,22,0,360 !创建管桩体
ALLSEL,ALL !选择所有实体
NUMCMP,ALL !压缩编号消除空号
/PNUM,VOLU,1 !显示体号码
VPLOT !显示所选体
VGLUE,2,5,1,4 !布尔运算粘结体所有的除桩之外的部分，即土层
ALLSEL,ALL !选择所有实体
NUMCMP,ALL !压缩编号消除空号
/PNUM,VOLU,1 !显示体号码
VPLOT !显示所选体
SAVE !
!采用扫掠方法划分网格
VSEL,S,,,3
AADD,21,30
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,5,0
LSEL,,,,29,32,1
LESIZE,all,,,2
VSEL,,,,3
TYPE,1
MAT,4
VSWEEP,3,15,35 !划分土层4
VSEL,S,,,5
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,1,0
LSEL,,,,21,24,1
LSEL,,,,51,54,1
LESIZE,all,,,2
VSEL,,,,5
TYPE,1
MAT,4
VSWEEP,5,21,20 !划分土层3
VSEL,S,,,1
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,18,0
LSEL,,,,13,16,1
LSEL,,,,45,48,1
LESIZE,all,,,2
VSEL,,,,1
TYPE,1
MAT,3
VSWEEP,1,20,1 !划分土层2
VSEL,S,,,4
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,4,0
LSEL,,,,5,8,1
LSEL,,,,41,44,1
LESIZE,all,,,2
VSEL,,,,4
TYPE,1
MAT,2
VSWEEP,4,1,28 !划分土层1
VSEL,S,,,2
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,300,0
LSEL,,,,63,70,1
LESIZE,all,,,2
VSEL,,,,2
TYPE,1
MAT,1
VSWEEP,2,22,6 !划分管桩
SAVE !
!接触命令流
ET,2,TARGE170 ! 定义目标单元类型
ET,3,CONTA173 ! 定义接触单元类型
KEYOPT,3,4,1
KEYOPT,3,9,0
KEYOPT,3,10,2
KEYOPT,3,12,2
R,1,,,1,1,,,
RMORE,,,,,,5 ! 定义接触单元实常数
! 创建管桩目标面1
ALLSEL,ALL
VSEL,S,,,2
ASLV,S
ASEL,S,,,23
TYPE,2
NSLA,S,1
ESURF,ALL
! 创建土层1、2、3接触面1
ALLSEL,ALL
VSEL,S,,,1
VSEL,S,,,4,5,1
ASLV,S
ASEL,S,,,32
ASEL,S,,,24
ASEL,S,,,33
TYPE,3
NSLA,S,1
ESURF,ALL
ET,4,TARGE170 ! 定义目标单元类型
ET,5,CONTA173 ! 定义接触单元类型
KEYOPT,5,4,1
KEYOPT,5,9,0
KEYOPT,5,10,2
KEYOPT,5,12,2
R,2,,,1,1,,,
RMORE,,,,,,5 ! 定义接触单元实常数
! 创建管桩目标面2
ALLSEL,ALL
VSEL,S,,,2
ASLV,S
ASEL,S,,,26
TYPE,4
NSLA,S,1
ESURF,ALL
! 创建土层1、2、3接触面2
ALLSEL,ALL
VSEL,S,,,1
VSEL,S,,,4,5,1
ASLV,S
ASEL,S,,,31
ASEL,S,,,25
ASEL,S,,,34
TYPE,5
NSLA,S,1
ESURF,ALL
SAVE !
!位移约束
/SOL
ALLSEL,ALL
VSEL,S,,,4
ASEL,S,,,2,5,1
DA,ALL,UX
DA,ALL,UY ! 在土层1四周施加X、Y方向约束
ALLSEL,ALL
VSEL,S,,,1
ASEL,S,,,7,10,1
DA,ALL,UX
DA,ALL,UY ! 在土层2四周施加X、Y方向约束
ALLSEL,ALL
VSEL,S,,,5
ASEL,S,,,11,14,1
DA,ALL,UX
DA,ALL,UY ! 在土层3四周施加X、Y方向约束
ALLSEL,ALL
VSEL,S,,,3
ASEL,S,,,16,19,1
DA,ALL,UX
DA,ALL,UY ! 在土层4四周施加X、Y方向约束
ALLSEL,ALL
VSEL,S,,,3
ASLV,S
ASEL,S,,,15
DA,ALL,ALL ! 在土层4底面施加X、Y、Z方向约束
SAVE !
!施加荷载
/SOLU
ACEL,0,0,9.8 !定义重力加速度
ALLSEL,ALL
ASEL,S,,,6
SFA,ALL,1,PRES,760000 ! 在管桩顶面施加荷载
ALLSEL,ALL
ASEL,S,,,28
SFA,ALL,1,PRES,40000 ! 在土层1顶面施加荷载
/STATUS,SOLU
SOLVE
FINISH
