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

我做的是一个桩土分析问题，想计算一下桩土相对滑动时的摩擦力，采用的是1/4三维实体模型，计划分为两个荷载步，第一步利用土体自重挤住桩，第二步把桩从土中拔出。 %jJG>T  
第一步可以算，第二步却不行，命令流如下，请大师们不吝指教，不胜感激。
s3N'0
2G  
FINI 4IK(7  
/CLEAR _{ue8kGt  
/filname,pile-soil contact z9f-.72"X  
/title，pile ,O5NLg-  
/units,si !国际单位制 1}+3dB_s  
fe=0.1 !摩擦系数 E*&vy  
fkn=10 !法向接触刚度 (le9q5Qr.  
fkt=3!切向接触刚度 Ha#=(9.  
ftoln=1 !初始渗透因子 810|Tj*U%  
fact=2!静动摩擦系数 Ng&%o  
/prep7 c?Y*Y   
et,1,solid45 - nm"of\o  
mp,ex,1,7.2e10         !桩的弹性模量，根据管壁厚度折算 .N;=\C
*  
mp,nuxy,1,0.3           !桩的泊松比 2YL?,uLS  
mp,dens,1,2000         !桩的密度 ;._ l0Jw  
!mp,mu,1,fe             !桩体性质 U)TUOwF  
mp,ex,2,2e8           !土的弹性模量，经验取值 >-?f0K  
mp,nuxy,2,0.3          !土的泊松比 299H$$WS,Z  
mp,dens,2,2000          !土的饱和密度   =>S]q71  
tb,dp,2 !vi>U|rh  
tbdata,1,0,30      !粘聚力c为0,摩擦角为30度 5PCqYN(:B  
!创建桩模型 D_2:k'4  
cylind,0,0.1,-0.1,7.1,0,90 J6"9v;V  
!创建土体模型 ]|pe>:gf'  
cylind,0.1,0.6,0,7,0,90 -]Bq|qTH[(  
!桩体划分网格 _oL?*ks  
allsel,all >tS'Q`R
 
vsel,s,,,1 umBICC]CU  
aslv,s d7^}tM  
lsla,s k?yoQL*  
lesize,1,,,6 b#c:u2  
lesize,2,,,6 r wL`Czs  
lesize,3,,,6 &N9 a<w8+  
lesize,4,,,6 1dY}\Sp  
lesize,5,,,6 zC:ASt  
lesize,6,,,6 PN%zIkbo  
lesize,7,0.05 b)#hSjWO#  
lesize,8,0.05 ^
S<Y>Nm]  
lesize,9,0.05 -:^U_FL8un  
type,1 Y>z>11yEB0  
mat,1 NSMyliM1Y  
mshape,0,3d W.jGGt\<\  
mshkey,1 BU)U/A8iS  
vmesh,all o)|flI'vT  
vsel,s,,,2 ！土层，划分网格 wVXS%4|v  
vplot ')Zvp7>$  
aslv,s &<g|gsG`  
lsla,s 7O2/z:$f  
lplot f^ZRT@`O  
lesize,10,,,6 8LJ8 }%*  
lesize,11,,,6 >~rTqtKd  
lesize,12,,,6 ,oe <  
lesize,13,,,6 ?
5__oT  
lesize,14,,,6 u]wZQl#-  
lesize,15,,,6 3d8L6GJ  
lesize,16,,,6 .8g)av+  
lesize,17,,,6 eu|YCYj)g  
lesize,18,0.05 Eh`7X=Z7E  
lesize,19,0.05 y8Ir@qp5  
lesize,20,0.05 Ufj`euY  
lesize,21,0.05 ?PxP% $hS  
type,1 m,28u3@r  
mat,2 ~hH REI&  
mshape,0,3d cU (D{~  
mshkey,1 w_c"@CjkE  
vmesh,all Y|m+dT6  
!创建接触对 X56q-|  
!用接触向导建立接触对 jwe*(k]z  
allsel wo}H'Q}Hj  
et,2,170 lgAoJ[  
et,3,173 *U-4Sy  
r,1 5<k"K^0QS  
real,1 ~Gp[_ %K  
keyopt,3,9,0    !消除初始渗透 h f)?1z4  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 .<?GS{6 N  
桩侧，建立目标面 mM~qBrwL  
r,1,,,fkn,ftoln, , , yF:1( 4  
rmore,,,,,,fkt @n/\L<]t  
RMORE,,,, ,, 0JS?;fk  
RMORE,,,fact,,,   iozt&~o  
mp,mu,3,fe t,Lrfv])  
!创建接触对 S>+|OCl";  
vsel,s,,,1 udH7}K v  
asel,s,,,3 234p9A@  
aplot @s>Czm5  
cm,target,area M.D1XX1/  
type,2 N];NAMp  
real,1 1nM  #kJ"  
mat,3 eeg)N1\  
nsla,s,1 ldcqe$7,  
esurf,all iXkF1r]i  
esel,s,type,,2 68|E9^`l  
eplot qbr$>xH  
/psymb,esys,1 iU918!!N  
allsel!很重要 urc| D0n  
!定义接触面 LP^$AAy  
vsel,s,,,2 +QavYqPF  
vplot H'5)UX@LP  
asel,s,,,9 A QU+mo  
aplot eIF5ZPSZi  
cm,contact,area SGRp3,1\4%  
type,3 ?,Xw[pR  
real,1 Jrf=@m\dk  
mat,3 ;O5zUl-`  
nsla,s,1 KkyVSoD\  
esurf,all Ty\R=y}}  
esel,s,type,,3 BZ#(   
eplot ;C#F>SG\S  
/psymb,esys,1 Y Uc+0  
allsel (#c*M?g3  
save pad*oPH,  
!侧面施加对称约束 f`(UQJ  
allsel gaxsv[W>^  
CSYS,0 V-P#1Kkh  
asel,s,loc,x,0 F;EwQjTF  
da,all,symm ;;Y!^^g  
asel,s,loc,y,0 P:S.~Jq  
da,all,symm NYUL:Tp  
allsel,all uc{Ihw  
!asel,s,,,2!没有约束顶面2 v"$L702d$\  
!da,all,all g/_5unI}u  
asel,s,,,6 YqD=>P[O  
da,all,ALL !TH) +zi  
!约束外侧面 ^e5=hH-%  
asel,S,,,8 Kn{4;Xk\  
NSLA,S,1 m 0C@G5  
CSYS,1 _ye |Y  
NROTAT,ALL X05/uX{  
D,ALL,UX！土侧采用的径向约束。 /N+dQe  
ALLSEL MKCsv+   
SAVE o[4}h:> dq  
finish!接触对定义结束 /HE
w-M9z  
!求解过程 l4YbKnp]  
/solu s[*rzoA  
!约束结束! (/YHk`v2  
!定义荷载步 #zy:a%  
!第一步荷载，重力荷载。 <nf@U>wlw  
acel,0,0,9.8 Es`Px_k  
nlgeom,on ! 打开大变形选项 !,uE]gwLw  
autot,on ! 打开自动时间步 s)t@ol  
time,1 e]aDP1n3t  
nsub,10 M?49TOQA  
outres,all,all wm@@$  
allsel ;d$rdFA_  
SOLVE j_[tu!~  
!第二部荷载,拉动桩 qq`4<0I>  
nsubst,100，500，5 +E+p"7  
outres,all,all
//MUeTxR  
time,101 C&f= ywi0  
nlgeom,on ! 打开大变形选项  dFc':|  
autot,on ! 打开自动时间步 l30EKoul)  
sdrfsrNvB-  
nsel,s,loc,z,7.1
?81c 4w  
D,ALL,UZ,0.1 ]cvwIc">  
lnsrch,on !自动线性搜索 @{e}4s?7od  
pred,on !打开预测求解器 0auYG><=  
allsel 3%|&I:tI  
SOLVE FUzzB94a  
i"FtcP^  

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

这个问题估计你是老大吧

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