桩从土中拔出的数值计算 请指教 [复制链接]

我做的是一个桩土分析问题，想计算一下桩土相对滑动时的摩擦力，采用的是1/4三维实体模型，计划分为两个荷载步，第一步利用土体自重挤住桩，第二步把桩从土中拔出。 8|gwH2st~  
第一步可以算，第二步却不行，命令流如下，请大师们不吝指教，不胜感激。 E`BL3+kQ  
FINI Z{|wjZb(  
/CLEAR +as(m  
/filname,pile-soil contact XK>B mq/]  
/title，pile {qK>A?
9  
/units,si !国际单位制 )D Y?Y-n  
fe=0.1 !摩擦系数 %kUIIHV}  
fkn=10 !法向接触刚度 }k$2r3  
fkt=3!切向接触刚度 |?g k%g  
ftoln=1 !初始渗透因子 =98@MX%P  
fact=2!静动摩擦系数 [+UF]m%W  
/prep7 bNi\+=v<Ys  
et,1,solid45 ?FJU>+{">  
mp,ex,1,7.2e10         !桩的弹性模量，根据管壁厚度折算 Ahm*_E2E  
mp,nuxy,1,0.3           !桩的泊松比 d=`hFwD9  
mp,dens,1,2000         !桩的密度 &G:#7HX@-  
!mp,mu,1,fe             !桩体性质 ;>bcI).  
mp,ex,2,2e8           !土的弹性模量，经验取值 YFeF(k!!n  
mp,nuxy,2,0.3          !土的泊松比 /g@!#Dt  
mp,dens,2,2000          !土的饱和密度   i.Yz)Bw  
tb,dp,2 +TL5yuA  
tbdata,1,0,30      !粘聚力c为0,摩擦角为30度 Bg{"{poy  
!创建桩模型 -Z9e}$q$,  
cylind,0,0.1,-0.1,7.1,0,90 NvY%sx,  
!创建土体模型 mGb,oj7l  
cylind,0.1,0.6,0,7,0,90 (V5_q,2  
!桩体划分网格 @uApm~}  
allsel,all 63 F@Ft  
vsel,s,,,1 Eu2@%2}P  
aslv,s ;.+sz(:hm  
lsla,s pKMy:j
 
lesize,1,,,6 f!AcBfaLr  
lesize,2,,,6 @uXF(KDX  
lesize,3,,,6 Yv\>\?865  
lesize,4,,,6 1?\G
6T  
lesize,5,,,6 )cxLpTr  
lesize,6,,,6 K_;'-
B  
lesize,7,0.05 Jsde+G,N  
lesize,8,0.05 R1)v;^B|)  
lesize,9,0.05 :+06M@  
type,1 A&XI1. j6  
mat,1 `ZhDoLpH<  
mshape,0,3d hX;JMQ915  
mshkey,1 e'Njl?>3  
vmesh,all Em?bV
(  
vsel,s,,,2 ！土层，划分网格 `saDeur#X  
vplot >|IUjv2L  
aslv,s >NDI<9<'0}  
lsla,s Gf*|f"O  
lplot sF[7pE  
lesize,10,,,6 &?59{B.mD  
lesize,11,,,6 :(ni/,~Q  
lesize,12,,,6 z$C}V/Ey  
lesize,13,,,6 9\y\{DHd  
lesize,14,,,6 iEjUo, Y[  
lesize,15,,,6 F|nJ3:v  
lesize,16,,,6 F/}(FG<'>I  
lesize,17,,,6 WTK )SKa,.  
lesize,18,0.05 H}vq2|MN  
lesize,19,0.05 SA!P:Q?h  
lesize,20,0.05 Qu]F<H*Y|  
lesize,21,0.05 ;&=c@>!xP#  
type,1 @M=xdZNyJ  
mat,2 B*B}eXUph  
mshape,0,3d xO3-I@  
mshkey,1 f_'#wc6  
vmesh,all X!6oviT|m  
!创建接触对 ,X^I]]  
!用接触向导建立接触对 *7cc4 wGQ  
allsel l<X8Ooan#{  
et,2,170 =zBc@VTp  
et,3,173 c{4Y?SSx  
r,1 Y~,ZBl,  
real,1 xQ~}9Kt\  
keyopt,3,9,0    !消除初始渗透 ,0k3Qi%  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 lfTDpKz3D  
桩侧，建立目标面 [ H|ifi  
r,1,,,fkn,ftoln, , , |%\>+/j$  
rmore,,,,,,fkt /fh[_!qN  
RMORE,,,, ,, cug=k  
RMORE,,,fact,,,   ey!QAEg"X1  
mp,mu,3,fe M4rI]^lJ  
!创建接触对 5=@q!8a*  
vsel,s,,,1 3Q;XvrGA  
asel,s,,,3
:$qa  
aplot KF!?;q0J  
cm,target,area *UxN~?
N|  
type,2 E)ne z  
real,1 u]`ur
#_  
mat,3 >_esLsPWh]  
nsla,s,1 "Zr+>a  
esurf,all V d`}F0WD  
esel,s,type,,2 K-X@3&X}  
eplot Q&\(m[:)  
/psymb,esys,1 ku*H*o~  
allsel!很重要 nI0TvBD  
!定义接触面 Wks?9)Is  
vsel,s,,,2 LKX; ^  
vplot ?xX9o  
asel,s,,,9 nNj<!}HvV  
aplot C]dK/~Z#r  
cm,contact,area L>@:Xo@  
type,3 Fx!NRY_  
real,1 2,T^L(]  
mat,3 ;;f&aujSHD  
nsla,s,1 +0DPhc  
esurf,all @T 5dPmn  
esel,s,type,,3 o%j[]P@4G  
eplot /
U@T#S  
/psymb,esys,1 #I &#x59  
allsel w%'8bH!  
save K (px-jY  
!侧面施加对称约束 LWX,u  
allsel 5oOF|IYi  
CSYS,0 {V
K   
asel,s,loc,x,0 {>r56\!F  
da,all,symm sR;^7(f!m  
asel,s,loc,y,0 Lkf}+aY  
da,all,symm _-6IB>  
allsel,all /l6r4aO2=  
!asel,s,,,2!没有约束顶面2 J n~t>?  
!da,all,all "~+?xke5z  
asel,s,,,6 SN<Dxa8Iy  
da,all,ALL |K(jXZ)  
!约束外侧面 0D==0n  
asel,S,,,8 SJ0IEPk  
NSLA,S,1 G_1`NyI  
CSYS,1 _+=M)lPm  
NROTAT,ALL V(#z{!  
D,ALL,UX！土侧采用的径向约束。 i!KZg74V  
ALLSEL =wcqCW,]  
SAVE **KkPjAO?  
finish!接触对定义结束 G?$0OU  
!求解过程 p3`odmbN  
/solu SSrYFu"  
!约束结束! ca6kqh"  
!定义荷载步 0pW?v:!H  
!第一步荷载，重力荷载。 yH<a;@C  
acel,0,0,9.8 4+1aW BJ2  
nlgeom,on ! 打开大变形选项 X6Wj,a  
autot,on ! 打开自动时间步 0r/
pZ3/  
time,1 U#U'iPy  
nsub,10 ^.?5!9U  
outres,all,all %G43g#pD  
allsel P-Upv6J3  
SOLVE 8n'"RaLQ8  
!第二部荷载,拉动桩 d&G#3}kOb%  
nsubst,100，500，5 @a]O(S>Ub  
outres,all,all }<=4A\LZ  
time,101 !Zi_4 .(4  
nlgeom,on ! 打开大变形选项 Z]^Ooy[pb  
autot,on ! 打开自动时间步 <$+Cd=71\  
Ms61FmA4  
nsel,s,loc,z,7.1 ZvVrbj&  
D,ALL,UZ,0.1 {vf4l4J(  
lnsrch,on !自动线性搜索 ^1 U<,<  
pred,on !打开预测求解器 rEfo)jod  
allsel *f ;">(`o*  
SOLVE aePLP  
|,)=-21&;  

我的qq号 57979269 知道的请指点 不胜感激

这个问题估计你是老大吧

让我试一试，应该是可以做的，估计是你的参数设置有问题。