菜鸟求助！！急

kappa119

小弟现在在做薄壁管桩的桩土分析，时间比较紧，由于ansys学的不精，请各位大侠多多指点。我写的命令流总是有问题，我感觉是接触单元的问题，但又不知道是什么问题，请大家帮我看看我的命令流，看哪里需要改进的，谢谢大家！


/filname,zhuangmo4
/title，121212
/config,nres,2000000
/config,fsplit,750
/view,1,1,1,1
/angle,1,65  
/units,si   
!接触系数参数确定
fe1=0.4            !摩擦系数0.4
fk1=-6.19e9              !法向接触刚度0.01~10,其下材料采用TB塑性定义时法向刚度将缩减100倍
                      !法向接触刚度，值越低越易收敛
fd=3                !切向接触刚度  
ftoln=0.1              !初始渗透因子，通常取值小于0.2，适用于K2取值0、2、3时，值越高越易收敛
fkt=-0.65e9          !由摩擦系数和法向刚度确定，K2取值3或K10等于1、2时fkt=MU*PRES/SLTO
/PREP7
!**有限元单元类型**
ET,1,SOLID65
ET,2,TARGE170
ET,3,CONTA173
KEYOPT,3,12,2
!**地层材料模型**
!***粉质粘土***
MP,DENS,1,1950
MP,EX,1,1E7
MP,NUXY,1,0.3
MP,MU,1,0.3 !设置摩擦系数为0.3
TB,DP,1
TBDATA,1,3.0E4,21
!***砼桩***
MP,DENS,2,2500
MP,EX,4,3e10
MP,NUXY,4,0.2
!**地层几何模型并划分网格**
cylind,50,200,0,-1400,0,90
cylind,,38,0,-1400,0,90
cylind,50,200,-1400,-1800,0,90
cylind,,38,-1400,-1800,0,90
cylind,38,50,-1400,-1800,0,90
ALLSEL,ALL !合并重复节点并重新排序
VGLUE,ALL
NUMCMP,ALL
!***桩模型***
cylind,38,50,0,-1400,0,90
!**有限元网格划分**
!***土层***
vsel,s,,,1
lesize,1,,,10
lesize,2,,,10
lesize,3,,,10
lesize,4,,,10
lesize,5,,,10
lesize,6,,,10
lesize,7,,,10
lesize,8,,,10
lesize,9,,,70
lesize,10,,,70
lesize,11,,,70
lesize,12,,,70
type,1
mat,1
vmesh,all
vsel,s,,,2
lesize,13,,,10
lesize,14,,,10
lesize,15,,,10
lesize,16,,,10
lesize,17,,,10
lesize,18,,,10
lesize,19,,,70
lesize,20,,,70
lesize,21,,,70
type,1
mat,1
vmesh,all
vsel,s,,,3
lesize,26,,,10
lesize,27,,,10
lesize,28,,,10
lesize,29,,,20
lesize,30,,,20
lesize,31,,,20
type,1
mat,1
vmesh,all
vsel,s,,,4
lesize,22,,,10
lesize,23,,,10
lesize,24,,,10
lesize,25,,,10
lesize,32,,,20
lesize,33,,,20
lesize,34,,,20
lesize,35,,,20
type,1
mat,1
vmesh,all
vsel,s,,,5
lesize,36,,,3
lesize,37,,,3
lesize,38,,,3
lesize,39,,,3
type,1
mat,1
vmesh,all
!***桩***
vsel,s,,,6
lesize,41,,,10
lesize,43,,,10
lesize,44,,,10
lesize,46,,,10
lesize,40,,,3
lesize,42,,,3
lesize,45,,,3
lesize,47,,,3
lesize,48,,,70
lesize,49,,,70
lesize,50,,,70
lesize,51,,,70
type,2
mat,1
vmesh,all
!***接触部分***
r,2,,,fk1,ftoln, 0,
rmore,,,1.0e20,0.0,1.0,fkt           !fkt=MU*PRES/SLTO,SLTO为接触对平均接触长度的1%
RMORE,10e3,0,1.0,,1.0,1e-8
RMORE,0,1e-8,1.0,0.0,,1.0

r,3,,,fk1,ftoln, 0,
rmore,,,1.0e20,0.0,1.0,fkt           !fkt=MU*PRES/SLTO,SLTO为接触对平均接触长度的1%
RMORE,10e3,0,1.0,,1.0,1e-8
RMORE,0,1e-8,1.0,0.0,,1.0
!******建立目标面*******
!****桩部分接触面****
ASEL,S,,,27
CM,_TARGET,AREA
NSLA,S,1 !选出需要加接触的节点
ESLN,S,0
TYPE,2 !桩为刚性面TARGET170单元
REAL,2
ESURF,ALL
!****土体部分接触面****
ASEL,S,,,4
CM,_TARGET,AREA
NSLA,S,1 !选出需要加接触的节点
ESLN,S,0
TYPE,3 !土为柔性面CONTAC173单元
REAL,2
ESURF,ALL
!****桩部分接触面****
ASEL,S,,,28
CM,_TARGET,AREA
NSLA,S,1 !选出需要加接触的节点
ESLN,S,0
TYPE,2 !桩为刚性面TARGET170单元
REAL,3
ESURF,ALL
!****土体部分接触面****
ASEL,S,,,9
CM,_TARGET,AREA
NSLA,S,1 !选出需要加接触的节点
ESLN,S,0
TYPE,3 !土为柔性面CONTAC173单元
REAL,3
ESURF,ALL
ALLSEL,ALL !全部选定，避免前面操作的漏选
FINISH
!*求解*
/SOLU
!**施加约束**
ASEL,S,LOC,X,0
DA,ALL,all

ASEL,S,LOC,Y,0
DA,ALL,ALL
ASEL,S,,,3
ASEL,A,,,20
DA,ALL,ALL
ASEL,S,LOC,Z,-1800
DA,ALL,ALL
csys,1
nsel,s,loc,x,500
d,all,all
ALLSEL,ALL !全部选定，避免前面操作的漏选
!**求解参数设置**
ANTYPE,STATIC !分析类型为静力
ACEL,,,9.8 !施加重力加速度9.8
NLGEOM,ON !打开大变形选项
NROPT,FULL,,ON !打开牛顿-辛普生选项，加快收敛速度
CNCHECK,ADJUST
SOLVE

rodge

如果觉得接触命令有问题，就用定义接触的向导模式，一般不会出错

kappa119

哦，多谢！！！