貌似网格划分有问题，计算不收敛，请高手指教，先谢过了

可圈可点

finish
/clear
/prep7
Et,1,solid45    !承台-桩
et,2,solid186    !土体
MP,EX,1,3.5e10    !承台
MP,nuxy,1,0.17
Mp,dens,1,2500
MP,EX,2,2.10e11   !桩
MP,nuxy,2,0.3
Mp,dens,2,7850

MP,EX,3,2.6e7       !淤泥质亚粘土
MP,nuxy,3,0.27
Mp,dens,3,790
tb,dp,2
Tbdata,1,11200,21,30    !定义DP材料及材料参数
MP,EX,4,7.4e6       !基岩，参数待定
MP,nuxy,4,0.35
Mp,dens,4,1990
MP,EX,5,7.4e6       !桩基混凝土，参数待定
MP,nuxy,5,0.35
Mp,dens,5,1990

K,1,0,0
k,2,40,0
k,3,40,20
k,4,0,20
*do,i,1,3,1
  l,i,i+1
*enddo
l,4,1
a,1,2,3,4
aplot
!在GUI中拉伸成体
VEXT,1,,,,,100

!建立第一排四根桩
wpoff,11,57.5,25
wprot,0,90,0
cyl4,0,0,1,,,,49.5           !第一排一桩
wpoff,6,0,0
cyl4,0,0,1,,,,49.5       !第一排二桩
wpoff,6,0,0
cyl4,0,0,1,,,,49.5   !第一排三桩
wpoff,6,0,0
cyl4,0,0,1,,,,49.5              !第一排四桩
wpcsys,-1
allsel,all,all
vplot
!复制排桩基
allsel
vgen,6,2,5,1,,,10
!减体
vsel,all
vplot
vsel,s,,,2,25,1
CM,zhuang,VOLU
allsel
VSBV,1,zhuang
vsel,all
!划分土体
allsel
wpcsys,-1
wpoff,0,5,0
wprota,,90
allsel
vsbw,26
allsel
vplot
!建立桩
!建立第一排四根桩
CSYS,0
wpcsys,-1
wpoff,11,57.5,25
wprot,0,90,0
cyl4,0,0,1,,,,49.5           !第一排一桩
wpoff,6,0,0
cyl4,0,0,1,,,,49.5       !第一排二桩
wpoff,6,0,0
cyl4,0,0,1,,,,49.5   !第一排三桩
wpoff,6,0,0
cyl4,0,0,1,,,,49.5              !第一排四桩
wpcsys,-1
allsel,all,all
vplot
!复制排桩基
allsel
vgen,6,3,6,1,,,10
allsel
!建立承台
CSYS,0
wpcsys,-1
wpoff,9,55.5,20
block,0,22,0,3.5,0,60
allsel
!承台与体进行减运算
vsel,all
vsel,s,,,3,26,1
CM,zhuang,VOLU
allsel
VSBV,27,zhuang,,DELETE,KEEP
nummrg,all
numcmp,all
allsel
!划分网格
!选择承台
vsel,all
allsel
VSEL,S,,,27
type,1
mat,1
esize,1
mshape,1,3d
mshkey,0
vmesh,all
!选择钢管桩
vsel,all
allsel
vsel,s,,,3,26
type,1
mat,2
esize,1
mshape,1,3d
mshkey,0
vmesh,all
!划分地基

allsel
vsel,s,,,1
vplot
type,2
mat,4
esize,1.2
mshape,1,3d
mshkey,0
vmesh,all

allsel
vplot
vsel,s,,,2
vplot
type,2
mat,3
esize,2
mshape,1,3d
mshkey,0
vmesh,all

!建立摩擦单元
/prep7
Et,7,170   !刚性面
Et,8,174     !柔性面
KEYOPT,8,9,0
KEYOPT,8,12,0
r,20        !定义单元实常数
real,20
Mp,mu,2,0.32

!定义刚性接触面
allsel
vplot
VSEL,s,MAT,,2
VSEL,R,LOC,Y,8,57.5   ! 选择体
ASLV,S                     !选中依附体上的面
aplot                      !显示所选中的面
Cm,target,are             ! 由所选择的（体面线点等等）组成一个元件
Type,7                    !给随后生成的单元激活一个单元类型号
mat,2                     !指定材料编号
real,20                    ! 单元实常数属性指示器
Nsla,s,1                   !指定仅面内节点背选中
esurf,all                   !在已存在的选中单元的自由面上覆盖产生单元
ESEL,S,TYPE,,7             !选择单元类型为7的单元
EPLOT                       !显示单元网格
/PSYMB,ESYS,1               ! 显示接触单元
asel,s,loc,y,57.5
nsla,s,1                   !选择面上的节点
esln,s,1                   !选择节点上的单元
eplot                      !显示单元
esel,r,type,,7              !选择单元类型为7  单元
eplot                  !显示单元
EDELE,all                    !删除桩顶面的接触单元
allsel
vplot
!!!!!!!!!!!!!!!!!!!!定义柔性接触面!!!!!!!!!!!!!!!!!!!!
VSEL,S,mat,,3,4
ASLV,S
ASEL,R,LOC,Z,24,76
ASEL,u,LOC,y,0,8
ASEL,u,LOC,x,0,10
ASEL,u,LOC,x,30,40
ASEL,u,LOC,y,20
APLOT
Cm,contact,area
type,8
real,20
nsla,s,1
nplot
esurf,all
ESEL,S,TYPE,,8
EPLOT
/PSYMB,ESYS,1
allsel
vplot
esel,s,type,,7
eplot
allsel
esel,s,type,,8
eplot
allsel
save
finish
!加载约束
/solu
vsel,s,mat,,3,4
vplot
aslv,s,1
aplot
asel,r,loc,z,0,0
aplot
da,all,all
allsel
vsel,s,mat,,3,4
vplot
aslv,s,1
aplot
asel,r,loc,z,100
aplot
da,all,all
vsel,s,mat,,3,4
vplot
aslv,s,1
aplot
asel,r,loc,x,0,0
aplot
da,all,all
allsel
vsel,s,mat,,3,4
vplot
aslv,s,1
aplot
asel,r,loc,x,40
aplot
da,all,all

allsel
vsel,s,mat,,3,4
vplot
aslv,s,1
aplot
asel,r,loc,y,0
aplot
da,all,all
allsel
vsel,all
allsel
vplot
!!!!!!!风暴潮浪力 601kN
f=601000
wpcsys,-1
nsel,s,loc,y,40.5
nplot
*GET,N_NUM,NODE,,COUNT
f,all,fx,f*24/N_NUM
allsel,all,all
vplot

!设置分析类型控制选项求解
antype,static      !静态分析
nlgeom,on         !计入大变形
autot,on           !打开自动时间步
acel,,9.8
time,1
nsub,20
OUTRES,ALL,ALL
solve
finish
!读入最后荷载步结果
/post1
set,last
esel,s,type,1
plnsol,u,x
plnsol,s,x      !查看桩-承台结果项
esel,s,type,2   
plnsol,u,z
plnsol,s,eqv      !查看土体结果项
esel,s,type,8       !查看接触状态
etable,st1,nmisc
pletab,st1
etable,cz1,nmisc,45  !查看接触压力
pletab,cz1