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

我做的是一个桩土分析问题，想计算一下桩土相对滑动时的摩擦力，采用的是1/4三维实体模型，计划分为两个荷载步，第一步利用土体自重挤住桩，第二步把桩从土中拔出。 gW>uR3Ca4  
第一步可以算，第二步却不行，命令流如下，请大师们不吝指教，不胜感激。 @k,z:~[C=  
FINI $8UUzk  
/CLEAR &O9 |#YUq  
/filname,pile-soil contact CE#gfP  
/title，pile /fZeWU0W  
/units,si !国际单位制 o4m\~as)Y  
fe=0.1 !摩擦系数 pGJ>O/%  
fkn=10 !法向接触刚度 ^@3sT,M,S  
fkt=3!切向接触刚度 =e9>FWf>  
ftoln=1 !初始渗透因子 LhOa{1SY  
fact=2!静动摩擦系数 +jS|2d  
/prep7 _jLL_GD  
et,1,solid45 fZ0M%f  
mp,ex,1,7.2e10         !桩的弹性模量，根据管壁厚度折算 (.D~0a JU  
mp,nuxy,1,0.3           !桩的泊松比 5ymk\Lw  
mp,dens,1,2000         !桩的密度 e+~\+:[?  
!mp,mu,1,fe             !桩体性质 [fAV5U  
mp,ex,2,2e8           !土的弹性模量，经验取值 v|uAzM{73  
mp,nuxy,2,0.3          !土的泊松比 xF;kTBRi  
mp,dens,2,2000          !土的饱和密度   +S>j0m<*  
tb,dp,2 #$18*?tLv|  
tbdata,1,0,30      !粘聚力c为0,摩擦角为30度 }4 )H   
!创建桩模型 L^}_~PO N5  
cylind,0,0.1,-0.1,7.1,0,90 F5P[dp-`1  
!创建土体模型 06O2:5zF  
cylind,0.1,0.6,0,7,0,90 j(m.$:  
!桩体划分网格 Fv~20G(O  
allsel,all yE:+Lo`>  
vsel,s,,,1 ZCmgs4W!  
aslv,s ^hc!FD  
lsla,s c6h.iBJ
'  
lesize,1,,,6 ,K9*%rW)  
lesize,2,,,6 llcb~  
lesize,3,,,6 ,{M^-3C  
lesize,4,,,6 i,>khc  
lesize,5,,,6 w,#>G07D  
lesize,6,,,6 o /j*d3  
lesize,7,0.05 w>pq+og&  
lesize,8,0.05 NQBpX  
lesize,9,0.05 qTdheX/  
type,1 vB.E3r=  
mat,1 1S{D6#bE  
mshape,0,3d &"yx<&c}  
mshkey,1 vdAr|4^qB  
vmesh,all L2\#w<d  
vsel,s,,,2 ！土层，划分网格 3EJj9}#x"'  
vplot 0?s|i :  
aslv,s "W6uV!  
lsla,s >V-A;S:  
lplot @.$|w>>T  
lesize,10,,,6 rDm~h~u5  
lesize,11,,,6 #.rdQ,)<  
lesize,12,,,6 ojaws+(& y  
lesize,13,,,6 K"|l@Q[  
lesize,14,,,6 dP3CG8w5  
lesize,15,,,6 ZP@ $Q%up  
lesize,16,,,6 EK{Eo9l  
lesize,17,,,6 mk.9OhYY  
lesize,18,0.05 |2O]R s  
lesize,19,0.05 ) i;1*jK  
lesize,20,0.05 zr-*$1eu  
lesize,21,0.05 .nrMfl_  
type,1 6 Q%jA7  
mat,2 8IlunJ  
mshape,0,3d v- 2:(IV  
mshkey,1 `=4r+  
vmesh,all BmbyH{4  
!创建接触对 W0dSsjNio  
!用接触向导建立接触对 o_XflzC  
allsel .c8g:WB<  
et,2,170 C'9Cr}cZ.  
et,3,173 arIf'
CG6  
r,1 a=J^  
real,1 uxXBEq;  
keyopt,3,9,0    !消除初始渗透 @5N]ZQ9  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 smlpD3?va  
桩侧，建立目标面 BF\XEm?!  
r,1,,,fkn,ftoln, , , Jl(&!?j  
rmore,,,,,,fkt LInz<bc<(  
RMORE,,,, ,, YWe{juXSw  
RMORE,,,fact,,,   &5\iM^
 
mp,mu,3,fe dG@%jD)  
!创建接触对 %RTBV9LIXr  
vsel,s,,,1 Lt u'W22  
asel,s,,,3 ?9!6%]2D  
aplot CyJEY-  
cm,target,area 95ZyP!  
type,2 T$>WE= Y  
real,1 9]k @Q_  
mat,3 }JF13beU  
nsla,s,1 U;YC}r  
esurf,all [$mHv,~  
esel,s,type,,2 {#ZlM  
eplot SWH2  
/psymb,esys,1 RSfQNc9Z  
allsel!很重要 2GP=&K/A  
!定义接触面 [)H&
'5 +F  
vsel,s,,,2 ,|3MG",@@h  
vplot
D9<!mH  
asel,s,,,9 N4v~;;@(  
aplot NSxoF3  
cm,contact,area n`#tKwWHYx  
type,3 H=<S 9M  
real,1 ,vr? 2k  
mat,3 BL0 {HV!  
nsla,s,1 caIL&G,  
esurf,all m4**~xfC  
esel,s,type,,3 bp* ^z,w  
eplot Zq^At+8+  
/psymb,esys,1 +[M6X} TQ  
allsel .!Oo|m`V@  
save R cAwrsd  
!侧面施加对称约束 CuFSeRe  
allsel DVG(Vw  
CSYS,0 {&cJDqz5=  
asel,s,loc,x,0 ^NRl//  
da,all,symm &q3"g*q  
asel,s,loc,y,0 FEW14U'O  
da,all,symm '9laa=H%8  
allsel,all fa-IhB1!K  
!asel,s,,,2!没有约束顶面2 N@2dA*T,  
!da,all,all \z>f
b%YW  
asel,s,,,6 ohRjvJ'v|  
da,all,ALL q3mJ782p]  
!约束外侧面 D[4u+g?[}>  
asel,S,,,8 r)lEofX,g+  
NSLA,S,1 Bn^0^J-  
CSYS,1 b+%f+zz*h  
NROTAT,ALL 3_ r*y9l  
D,ALL,UX！土侧采用的径向约束。 r A`V}>Xj  
ALLSEL CnU*Jb  
SAVE P2NQHX  
finish!接触对定义结束 eX?OYDDC0j  
!求解过程 Tl%`P_J)-S  
/solu 02f~En}>6  
!约束结束! lNy.g{2f<m  
!定义荷载步 ;!=G   
!第一步荷载，重力荷载。 UxbjA- U[  
acel,0,0,9.8 Ok|*!
!T  
nlgeom,on ! 打开大变形选项 8hu<E4]L  
autot,on ! 打开自动时间步 sQ=]NF)\  
time,1 hB"fhX  
nsub,10 {Bk[rCl  
outres,all,all P60~V"/P  
allsel >W%Em
nLK  
SOLVE p?*Q- f  
!第二部荷载,拉动桩 iIvc43YV%  
nsubst,100，500，5 9%k2'iV7  
outres,all,all zpzK>DH(  
time,101 z
kt+7,vI  
nlgeom,on ! 打开大变形选项 <->{  
autot,on ! 打开自动时间步 a+[RS]le  
J28
M@cn  
nsel,s,loc,z,7.1 Tre]"2l  
D,ALL,UZ,0.1 SbND Y{5RO  
lnsrch,on !自动线性搜索 !F*5M1Kjd  
pred,on !打开预测求解器 7TgOK   
allsel \MsTB|Z  
SOLVE GD
&uQ`Y5  
.!Qki@  

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

这个问题估计你是老大吧

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