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

我做的是一个桩土分析问题，想计算一下桩土相对滑动时的摩擦力，采用的是1/4三维实体模型，计划分为两个荷载步，第一步利用土体自重挤住桩，第二步把桩从土中拔出。 _eE hIQ9  
第一步可以算，第二步却不行，命令流如下，请大师们不吝指教，不胜感激。 O'OVj  
FINI ^Iz(V2  
/CLEAR eVVm"96Q.;  
/filname,pile-soil contact *s"dCc  
/title，pile h)MU^aP  
/units,si !国际单位制 >H*?ktcW  
fe=0.1 !摩擦系数 J=*X%^jX9Z  
fkn=10 !法向接触刚度 :JEzfI1  
fkt=3!切向接触刚度 <DM /"^*  
ftoln=1 !初始渗透因子 0rjH`H]M  
fact=2!静动摩擦系数 u$"dL=s!  
/prep7 xjplJ'jB  
et,1,solid45 ^F*G  
mp,ex,1,7.2e10         !桩的弹性模量，根据管壁厚度折算 $jL{l8x  
mp,nuxy,1,0.3           !桩的泊松比 JS&=V67[  
mp,dens,1,2000         !桩的密度 k
s=l Nz9  
!mp,mu,1,fe             !桩体性质 Q|O! cEW/  
mp,ex,2,2e8           !土的弹性模量，经验取值 QNk\y@yKw  
mp,nuxy,2,0.3          !土的泊松比 xux
j
 
mp,dens,2,2000          !土的饱和密度   ?/SIA9VK  
tb,dp,2 )sh+cfTCb  
tbdata,1,0,30      !粘聚力c为0,摩擦角为30度 "'@iDq%y  
!创建桩模型 O!3MXmaO  
cylind,0,0.1,-0.1,7.1,0,90 5$(qnOi  
!创建土体模型 aQzu[N  
cylind,0.1,0.6,0,7,0,90 }=+J&cR  
!桩体划分网格 ;Tn$c70  
allsel,all I1IuvH6  
vsel,s,,,1 w=]bj0<A=  
aslv,s S:*.,zC  
lsla,s IC1oW)  
lesize,1,,,6 :z^VI M  
lesize,2,,,6 )o:%Zrk  
lesize,3,,,6 aBH!K   
lesize,4,,,6 -UHa;WH  
lesize,5,,,6 )_|;h2I  
lesize,6,,,6 vw;GbQH(  
lesize,7,0.05 7B|ddi7Q>  
lesize,8,0.05 #`1@4,iC  
lesize,9,0.05 }bW"Z2^nB  
type,1 Y6>@zznk  
mat,1 .5i\L OTd  
mshape,0,3d rYp]RX>  
mshkey,1 /!V)2j,  
vmesh,all ^UB<U#8,  
vsel,s,,,2 ！土层，划分网格 $rb #k{  
vplot V*uoGWL]+  
aslv,s y)s+/Teb  
lsla,s aCq ) hR  
lplot .ykCmznf*  
lesize,10,,,6 JDQ7  
lesize,11,,,6 Wsz-#kc\[  
lesize,12,,,6 !,
}F2z?4c  
lesize,13,,,6 8BL]]gT-I  
lesize,14,,,6 }67lL~L  
lesize,15,,,6 'Mfn:n+  
lesize,16,,,6 x*Lm{c5+  
lesize,17,,,6 sS 5aJ}Qs  
lesize,18,0.05 XS L*e  
lesize,19,0.05 u9TiEEof3  
lesize,20,0.05 iiWs]5  
lesize,21,0.05 Znv3h  
type,1 mHm"QBa!  
mat,2 P6n9yJ$,cb  
mshape,0,3d d)_fI*:f  
mshkey,1 vZDM}u  
vmesh,all q;W(;B  
!创建接触对 s/;iZiWK  
!用接触向导建立接触对 q_h=O1W  
allsel ^7qqO%  
et,2,170 /$q9 Kxb  
et,3,173 7S +YQ$_  
r,1 B4hT(;
k  
real,1 [VD)DO5  
keyopt,3,9,0    !消除初始渗透 -^Pn4
y]A)  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 gPo3jwo$  
桩侧，建立目标面 ,)S(SnCF  
r,1,,,fkn,ftoln, , , cgQ4JY/6  
rmore,,,,,,fkt aBBTcN%'  
RMORE,,,, ,, ?)8OC(B8q  
RMORE,,,fact,,,   =B+^-2G8  
mp,mu,3,fe TA2?Ia;@xV  
!创建接触对 :`pgdn  
vsel,s,,,1 u%o2BLx  
asel,s,,,3 3gUGfedi  
aplot p41TSALq  
cm,target,area 5IwQ<V  
type,2 Lrr1) h  
real,1 i`FskEoijq  
mat,3 -f{NVX\<0  
nsla,s,1 ^}F@*A;o  
esurf,all ]TE,N$X  
esel,s,type,,2 qu}&4_`%:V  
eplot U_X/  
/psymb,esys,1 g RX`61  
allsel!很重要 f  _ O  
!定义接触面 ~{8X$xs  
vsel,s,,,2 )L
,.KO  
vplot ?`,Rkg0fe  
asel,s,,,9 y;.U-}e1  
aplot QR.]?t;1  
cm,contact,area
i^i^g5l!  
type,3 a*qf\&Vb|  
real,1 F$|Ec9  
mat,3 +d$l1j  
nsla,s,1 62}rZVJq  
esurf,all 6 4fB$  
esel,s,type,,3 l]bCt b%_  
eplot l2AAEB_C.  
/psymb,esys,1 =GP~h*5es  
allsel b?=>)':f  
save |uT&M`7\{  
!侧面施加对称约束 Po~u-5  
allsel I^CKq?V?:  
CSYS,0 0~)_/yx?S  
asel,s,loc,x,0 qoifzEc`U  
da,all,symm B9Ha6kj  
asel,s,loc,y,0 2l8TX#K  
da,all,symm JdP[ cN  
allsel,all Ah_Ttj  
!asel,s,,,2!没有约束顶面2 KMhEU**  
!da,all,all UP8=V>T02  
asel,s,,,6 5bfd8C  
da,all,ALL YFGQPg  
!约束外侧面 ?HsQ417.H  
asel,S,,,8 _olhCLIR-  
NSLA,S,1 }"sZ)FE  
CSYS,1 "#3p=}]  
NROTAT,ALL l6AG!8H  
D,ALL,UX！土侧采用的径向约束。 o#X=1us  
ALLSEL {m<NPtp910  
SAVE m^9[k,;K  
finish!接触对定义结束 miq"3  
!求解过程 _YT9zG  
/solu !,I530eh7  
!约束结束! o:.6{+|N  
!定义荷载步 HxH.=M8S_  
!第一步荷载，重力荷载。 ssWSY(j]  
acel,0,0,9.8 No'^]r  
nlgeom,on ! 打开大变形选项 ITq$8  
autot,on ! 打开自动时间步 0zL7$Q#c  
time,1 8.g(&F  
nsub,10 Xv5|j/<~p  
outres,all,all 5@:c6(5$  
allsel 1z#0CX}Y/H  
SOLVE /.Fvl;!
J;  
!第二部荷载,拉动桩 ,rO>5$w.  
nsubst,100，500，5 PCx] >&  
outres,all,all y5v}EX`m&  
time,101 #wvmVB.5~  
nlgeom,on ! 打开大变形选项 ^EG@tB $<  
autot,on ! 打开自动时间步 bSKe@4C  
FTA[O.tiG  
nsel,s,loc,z,7.1 \0^ZNa?  
D,ALL,UZ,0.1 =s\RK   
lnsrch,on !自动线性搜索 Py3Y*YP  
pred,on !打开预测求解器 0VA$ Ige  
allsel ?puZqVu5  
SOLVE Hf ]w  
_i-
(`5  

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

这个问题估计你是老大吧

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