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

我做的是一个桩土分析问题，想计算一下桩土相对滑动时的摩擦力，采用的是1/4三维实体模型，计划分为两个荷载步，第一步利用土体自重挤住桩，第二步把桩从土中拔出。 'C2X9/!,  
第一步可以算，第二步却不行，命令流如下，请大师们不吝指教，不胜感激。 # zbAA<f  
FINI ?;0w1  
/CLEAR z?DI4O#Up  
/filname,pile-soil contact 7a_tT;f;  
/title，pile n,_q6/!  
/units,si !国际单位制 0Kytg\p}  
fe=0.1 !摩擦系数 <Cbi5DtR  
fkn=10 !法向接触刚度 lIUaGz|  
fkt=3!切向接触刚度 IGQcQ/M  
ftoln=1 !初始渗透因子 2]}4)_&d<e  
fact=2!静动摩擦系数 -/'_XR@1  
/prep7 [Ep%9(SgA'  
et,1,solid45 <(c_[o/  
mp,ex,1,7.2e10         !桩的弹性模量，根据管壁厚度折算 L]kd.JJvy  
mp,nuxy,1,0.3           !桩的泊松比 5mYX#//:  
mp,dens,1,2000         !桩的密度 r&/M')}?Lw  
!mp,mu,1,fe             !桩体性质
Qu;cl/&  
mp,ex,2,2e8           !土的弹性模量，经验取值 9{KL^O?g  
mp,nuxy,2,0.3          !土的泊松比 'OTQiI^t=  
mp,dens,2,2000          !土的饱和密度   ,ut-Di=6  
tb,dp,2 * ",/7(  
tbdata,1,0,30      !粘聚力c为0,摩擦角为30度 CVt:tV  
!创建桩模型 rd:WF(]  
cylind,0,0.1,-0.1,7.1,0,90 f
)x(sk  
!创建土体模型 G=ly .  
cylind,0.1,0.6,0,7,0,90 {gkzo3  
!桩体划分网格 =G,wR'M  
allsel,all EQTJ=\WFF  
vsel,s,,,1 V*<`!w  
aslv,s Oe~x,=X)  
lsla,s ("r\3Mvs  
lesize,1,,,6 l{o,"P"  
lesize,2,,,6 .V   
lesize,3,,,6 LpYG!Kl  
lesize,4,,,6 Rl/5eE8  
lesize,5,,,6 !ay:
h Iv  
lesize,6,,,6 5w+KIHhN|  
lesize,7,0.05 p.^qB]%  
lesize,8,0.05 $&~moAl  
lesize,9,0.05  B8~JUGD  
type,1 2t,N9@u=UN  
mat,1 X;&
Iu{&=  
mshape,0,3d J{!U;r!6  
mshkey,1 /f}!G  
vmesh,all |Fi{]9(G2  
vsel,s,,,2 ！土层，划分网格 je`Ysben  
vplot =f/CBYNw@V  
aslv,s JJZu%9~[  
lsla,s 0;Oe&Y  
lplot >2t.7UhDI  
lesize,10,,,6 yCvP-?2  
lesize,11,,,6 4Y5lP00!}  
lesize,12,,,6 srCpgs]h  
lesize,13,,,6 |8q:sr_  
lesize,14,,,6 -Is;cbfLj/  
lesize,15,,,6 .US=fWyrb  
lesize,16,,,6 j"F?^0aR,Q  
lesize,17,,,6 ~~\C.6c#  
lesize,18,0.05 L}a-c(G+8  
lesize,19,0.05 ?`lD|~  
lesize,20,0.05 &pzf*|}  
lesize,21,0.05 \5iMr[s  
type,1
}NJKkj?  
mat,2 RH}i=  
mshape,0,3d qeUT]* w  
mshkey,1 {U'\2Ge<m  
vmesh,all QJ,[K_  
!创建接触对 tP][o494\&  
!用接触向导建立接触对 5(=5GkE)>  
allsel B%^W$7 q  
et,2,170 9,wD  
et,3,173 bt{b%r  
r,1 uo8
[,'  
real,1 Ls`[7w  
keyopt,3,9,0    !消除初始渗透 omMOA  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 MO>9A,&f  
桩侧，建立目标面 Cvp!(<<gK  
r,1,,,fkn,ftoln, , , 9$?Sts}6&  
rmore,,,,,,fkt CBr(a'3{Z  
RMORE,,,, ,, , UsY0YC  
RMORE,,,fact,,,   3%[;nhbA7  
mp,mu,3,fe i$5<>\g  
!创建接触对 (hJ&`Tt  
vsel,s,,,1 2z9s$tp  
asel,s,,,3 4OaU1Y[  
aplot "P9(k>  
cm,target,area J6I:UML  
type,2 PS}'LhZ  
real,1 [} zzG@g,J  
mat,3 KcvstC`  
nsla,s,1 kz\Ss|jl  
esurf,all )! eJW(  
esel,s,type,,2 \47djm
G-  
eplot AxtmG\o>  
/psymb,esys,1 r@a]fTf  
allsel!很重要 X`6"^ xme  
!定义接触面 YO'a
X  
vsel,s,,,2 7 'q *(v  
vplot bEKhU\@=J  
asel,s,,,9 QdrZi.qKH  
aplot 2O 2HmL  
cm,contact,area smUSR4VK  
type,3 21$E.x 6  
real,1 @9tzk [  
mat,3 nSv@FT'~z  
nsla,s,1 <I#nwoHN  
esurf,all = ;cTm5d;T  
esel,s,type,,3 w7@TM%nS  
eplot s(Bcw`'#  
/psymb,esys,1 |z7V1xF  
allsel @PT([1C  
save k5%W8dI  
!侧面施加对称约束 ZuFcJ?8i  
allsel >t_h/:JZ)  
CSYS,0 Vak\N)=u  
asel,s,loc,x,0 "2~L  
da,all,symm L</k+a?H!  
asel,s,loc,y,0 oo\7\b#Jx  
da,all,symm GoLK 95"]  
allsel,all $<QrV,T  
!asel,s,,,2!没有约束顶面2 @jxP3:s  
!da,all,all V,h}l"  
asel,s,,,6 QR2S6
7-  
da,all,ALL (^NYC$ZxM=  
!约束外侧面 ~].?8C.>*  
asel,S,,,8 \$4 [qG=  
NSLA,S,1 CkV5PU  
CSYS,1 )_YB8jUR-X  
NROTAT,ALL 5%%e$o+  
D,ALL,UX！土侧采用的径向约束。 o(k{Ed  
ALLSEL 4`B3Kt`o  
SAVE 'Qh1$X)R7a  
finish!接触对定义结束 _a#k3r  
!求解过程 T-LX>*  
/solu ,v%'2[}  
!约束结束! kV+%(Gl8  
!定义荷载步 vSu dT  
!第一步荷载，重力荷载。 c'.XC}  
acel,0,0,9.8 KdBpfPny@  
nlgeom,on ! 打开大变形选项 jF}-dfe  
autot,on ! 打开自动时间步 >qz#&  
time,1 Kl,NL]]4*5  
nsub,10 "Ccyj/  
outres,all,all U`aB&[=$  
allsel %s! |,Cu  
SOLVE k2@]nW"S  
!第二部荷载,拉动桩 H76iBJ66  
nsubst,100，500，5 'u:
-~nSX)  
outres,all,all s IFE:/1,  
time,101 BsU}HuQZQ  
nlgeom,on ! 打开大变形选项 g<N;31:c\  
autot,on ! 打开自动时间步 ,v<7O_A/e  
^)(-7H  
nsel,s,loc,z,7.1 ]rG/?1'^i  
D,ALL,UZ,0.1 .* )e24`  
lnsrch,on !自动线性搜索 /9e?uC6  
pred,on !打开预测求解器 .
P <3+  
allsel +CM>]Ze  
SOLVE byFO^pce  
4*ZY#7h  

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

这个问题估计你是老大吧

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