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

我做的是一个桩土分析问题，想计算一下桩土相对滑动时的摩擦力，采用的是1/4三维实体模型，计划分为两个荷载步，第一步利用土体自重挤住桩，第二步把桩从土中拔出。 _xo;[rEw8  
第一步可以算，第二步却不行，命令流如下，请大师们不吝指教，不胜感激。 #\3X;{  
FINI )=#zMdK&  
/CLEAR 9%?a\#C  
/filname,pile-soil contact {_jbFJ  
/title，pile 9u_D@A"aC`  
/units,si !国际单位制 7<93n`byM  
fe=0.1 !摩擦系数 :-T[)Q+-3  
fkn=10 !法向接触刚度 D.w6/DxaXa  
fkt=3!切向接触刚度 )JR&  
ftoln=1 !初始渗透因子 sN/8OLc  
fact=2!静动摩擦系数 .JWN\\  
/prep7 ?T>)7Y)  
et,1,solid45 ~HUO$*U4<  
mp,ex,1,7.2e10         !桩的弹性模量，根据管壁厚度折算 D[K!xq  
mp,nuxy,1,0.3           !桩的泊松比 "Q1oSpF  
mp,dens,1,2000         !桩的密度 \8Yv}wQ  
!mp,mu,1,fe             !桩体性质 B5#>ieM*  
mp,ex,2,2e8           !土的弹性模量，经验取值 iGIry^D  
mp,nuxy,2,0.3          !土的泊松比 f"G-',O<  
mp,dens,2,2000          !土的饱和密度   j<d,7  
tb,dp,2 AZbFj-^4  
tbdata,1,0,30      !粘聚力c为0,摩擦角为30度 *V@t]d$=#  
!创建桩模型 )p 8P\Rl  
cylind,0,0.1,-0.1,7.1,0,90 34aSRFsk*  
!创建土体模型 4{ED~w|  
cylind,0.1,0.6,0,7,0,90 .oe\wJS6  
!桩体划分网格 zIc_'Z,b  
allsel,all WWVQJ{,}  
vsel,s,,,1 9G{#a#Z.  
aslv,s sbq:8P#  
lsla,s S2~im?^21  
lesize,1,,,6 v}u
zUY  
lesize,2,,,6 Vx=tP.BO]  
lesize,3,,,6 XlxM
.;i0H  
lesize,4,,,6 a)r
T3gl  
lesize,5,,,6 81S0:=   
lesize,6,,,6 VP6_}9:9   
lesize,7,0.05 ,egbU(:l  
lesize,8,0.05 Ak@y"!wnM  
lesize,9,0.05 eY'RDQa  
type,1 b'(Hwc\ t  
mat,1 J L2g!n= K  
mshape,0,3d idY Xv)R  
mshkey,1 nMH:7[x3  
vmesh,all cIXqnb  
vsel,s,,,2 ！土层，划分网格 cbv%1DT3  
vplot d'e\
tO  
aslv,s f(*iagEy  
lsla,s ~Z$Ro/;l  
lplot d@+u&xrd  
lesize,10,,,6 X;yThb`iI  
lesize,11,,,6 7m;<b$  
lesize,12,,,6 Y]+KsiOL  
lesize,13,,,6 E{IY7Xz^>  
lesize,14,,,6 }_9yemP  
lesize,15,,,6 {:`XhPS<B  
lesize,16,,,6 ~{{@m]P  
lesize,17,,,6 l#3jJn  
lesize,18,0.05 x8RiYi+  
lesize,19,0.05 (Cbm*VL  
lesize,20,0.05 fI-f Gx  
lesize,21,0.05 H:,Hr_;nC  
type,1 p!+bn,?G  
mat,2 (6X{ &  
mshape,0,3d JYa3xeC;  
mshkey,1 _*u$U  
vmesh,all @BI;H V%k  
!创建接触对 p)*x7~3e  
!用接触向导建立接触对 @eN,m {b  
allsel @6{F4  
et,2,170 GBGna3  
et,3,173 -}"nb-RR\  
r,1 l+<AM%U\ V  
real,1 N=c{@h  
keyopt,3,9,0    !消除初始渗透 wo\O0?d3{  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 } v3w-  
桩侧，建立目标面 #++D|oE  
r,1,,,fkn,ftoln, , , \qB.>f"%p|  
rmore,,,,,,fkt QzV%m0  
RMORE,,,, ,, r/RX|M  
RMORE,,,fact,,,   g ~>nT>6  
mp,mu,3,fe NdM \RD_R  
!创建接触对 Rh{`#dI~=  
vsel,s,,,1 Gn 1  
asel,s,,,3 5ih5=qX  
aplot Twpk@2=l  
cm,target,area 
/O[Z  
type,2 rh6 e  
real,1 ?fNUmk^A<  
mat,3 M8f[ck  
nsla,s,1 m%)S<L7 l  
esurf,all LBnlaH.  
esel,s,type,,2 4gn|zSe>^  
eplot ]N6
UY  
/psymb,esys,1 >\!k~Zi  
allsel!很重要 G(

#EW+  
!定义接触面 GCp90  
vsel,s,,,2 e}2[g  
vplot MN_1^T5  
asel,s,,,9 R{{?wr6b$  
aplot lYy:A%yDT  
cm,contact,area fTd":F  
type,3 e J:#vX86  
real,1 <E7y:%L[Go  
mat,3 ,A>cL#Oe  
nsla,s,1 @g1T??h   
esurf,all 
dXDuO  
esel,s,type,,3 ^'4I%L"  
eplot 6\6g-1B`  
/psymb,esys,1 cp4~`X  
allsel o>
yo9n%t  
save &IkHP/  
!侧面施加对称约束 q`HK4~i,  
allsel s~g]`/h$r  
CSYS,0 rz5@E  
asel,s,loc,x,0 f!JSb?#3  
da,all,symm pm:#@sl  
asel,s,loc,y,0 gTg[!}_;\N  
da,all,symm >,nK  
allsel,all 2Kw i4R  
!asel,s,,,2!没有约束顶面2 X'usd$[.  
!da,all,all 2o~UA\:+=  
asel,s,,,6 Q8NrbMrl  
da,all,ALL fP{IW`t}]  
!约束外侧面 bl4I4RB  
asel,S,,,8 ]v0=jm5A  
NSLA,S,1 ^-s`$lTp  
CSYS,1 ?n
D]p!  
NROTAT,ALL zu``F]B  
D,ALL,UX！土侧采用的径向约束。 CZ,2Rq  
ALLSEL 6x@4gPy[  
SAVE *.+N?%sAP)  
finish!接触对定义结束 %`]fZr A]#  
!求解过程 (_eM:H=e>  
/solu %wIb@km  
!约束结束! Uk4G9}I  
!定义荷载步 U (A#}  
!第一步荷载，重力荷载。 aY\(R02B  
acel,0,0,9.8 `&SBp }W}  
nlgeom,on ! 打开大变形选项 cd)}a_9  
autot,on ! 打开自动时间步 ?P Mi#H  
time,1 q(}
#{OO  
nsub,10 ;et(Yi;9  
outres,all,all FK;2u$
:  
allsel !FeNx*31i  
SOLVE q=+wQ[a<  
!第二部荷载,拉动桩 yH.Z%*=xQa  
nsubst,100，500，5 g3\13<  
outres,all,all @0iXqM#jH  
time,101 J%8hf%!ud  
nlgeom,on ! 打开大变形选项 'PdUSv|lH  
autot,on ! 打开自动时间步 T@Th?  
-p f9Wk  
nsel,s,loc,z,7.1 \`8?=_ST  
D,ALL,UZ,0.1 [.>g.p,;  
lnsrch,on !自动线性搜索 +$B#] ,  
pred,on !打开预测求解器 $H1igYc  
allsel [ejl #'*5  
SOLVE S5wkBdr{  
42wZy|oqp  

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

这个问题估计你是老大吧

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