桩土接触分析请大师指教

wxq9822

我做的是一个桩土分析问题，想计算一下桩土相对滑动时的摩擦力，采用的是1/4三维实体模型，计划分为两个荷载步，第一步利用土体自重挤住桩，第二步把桩从土中拔出。
第一步可以算，第二步却不行，命令流如下，请大师们不吝指教，不胜感激。
FINI
/CLEAR
/filname,pile-soil contact
/title，pile
/units,si !国际单位制
fe=0.1 !摩擦系数
fkn=10 !法向接触刚度
fkt=3!切向接触刚度
ftoln=1 !初始渗透因子
fact=2!静动摩擦系数
/prep7
et,1,solid45
mp,ex,1,7.2e10         !桩的弹性模量，根据管壁厚度折算
mp,nuxy,1,0.3           !桩的泊松比
mp,dens,1,2000         !桩的密度
!mp,mu,1,fe             !桩体性质
mp,ex,2,2e8           !土的弹性模量，经验取值
mp,nuxy,2,0.3          !土的泊松比
mp,dens,2,2000          !土的饱和密度   
tb,dp,2
tbdata,1,0,30      !粘聚力c为0,摩擦角为30度
!创建桩模型
cylind,0,0.1,-0.1,7.1,0,90
!创建土体模型
cylind,0.1,0.6,0,7,0,90
!桩体划分网格
allsel,all
vsel,s,,,1
aslv,s
lsla,s
lesize,1,,,6
lesize,2,,,6
lesize,3,,,6
lesize,4,,,6
lesize,5,,,6
lesize,6,,,6
lesize,7,0.05
lesize,8,0.05
lesize,9,0.05
type,1
mat,1
mshape,0,3d
mshkey,1
vmesh,all
vsel,s,,,2 ！土层，划分网格
vplot
aslv,s
lsla,s
lplot
lesize,10,,,6
lesize,11,,,6
lesize,12,,,6
lesize,13,,,6
lesize,14,,,6
lesize,15,,,6
lesize,16,,,6
lesize,17,,,6
lesize,18,0.05
lesize,19,0.05
lesize,20,0.05
lesize,21,0.05
type,1
mat,2
mshape,0,3d
mshkey,1
vmesh,all
!创建接触对
!用接触向导建立接触对
allsel
et,2,170
et,3,173
r,1
real,1
keyopt,3,9,0    !消除初始渗透
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开
桩侧，建立目标面
r,1,,,fkn,ftoln, , ,
rmore,,,,,,fkt
RMORE,,,, ,,
RMORE,,,fact,,,  
mp,mu,3,fe
!创建接触对
vsel,s,,,1
asel,s,,,3
aplot
cm,target,area
type,2
real,1
mat,3
nsla,s,1
esurf,all
esel,s,type,,2
eplot
/psymb,esys,1
allsel!很重要
!定义接触面
vsel,s,,,2
vplot
asel,s,,,9
aplot
cm,contact,area
type,3
real,1
mat,3
nsla,s,1
esurf,all
esel,s,type,,3
eplot
/psymb,esys,1
allsel
save
!侧面施加对称约束
allsel
CSYS,0
asel,s,loc,x,0
da,all,symm
asel,s,loc,y,0
da,all,symm
allsel,all
!asel,s,,,2!没有约束顶面2
!da,all,all
asel,s,,,6
da,all,ALL
!约束外侧面
asel,S,,,8
NSLA,S,1
CSYS,1
NROTAT,ALL
D,ALL,UX！土侧采用的径向约束。
ALLSEL
SAVE
finish!接触对定义结束
!求解过程
/solu
!约束结束!
!定义荷载步
!第一步荷载，重力荷载。
acel,0,0,9.8
nlgeom,on ! 打开大变形选项
autot,on ! 打开自动时间步
time,1
nsub,10
outres,all,all
allsel
SOLVE
!第二部荷载,拉动桩
nsubst,100，500，5
outres,all,all
time,101
nlgeom,on ! 打开大变形选项
autot,on ! 打开自动时间步

nsel,s,loc,z,7.1
D,ALL,UZ,0.1
lnsrch,on !自动线性搜索
pred,on !打开预测求解器
allsel
SOLVE

253276391

谢谢你的共享！

zsc1984001

谢谢不错的分享

wsf2008

你可以不要自动时间步，自己设定时间步看看，这条建议未必成熟，你可以试一下

zzllww9

谢谢你的共享！

zem

自动线性搜索
能够有效帮助收敛

21people

fe=0.1 !摩擦系数
fkn=10 !法向接触刚度
fkt=3!切向接触刚度
ftoln=1 !初始渗透因子
fact=2!静动摩擦系数
楼主这些参数哪来的？有参考意义么？

秋之风景

刚刚接触，谢谢分享！