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

我做的是一个桩土分析问题，想计算一下桩土相对滑动时的摩擦力，采用的是1/4三维实体模型，计划分为两个荷载步，第一步利用土体自重挤住桩，第二步把桩从土中拔出。 f61]`@Bk  
第一步可以算，第二步却不行，命令流如下，请大师们不吝指教，不胜感激。 ug>]U ~0  
FINI P}b Dn;  
/CLEAR !vi4* @:  
/filname,pile-soil contact @#J H=-06  
/title，pile M|aQ)ivh3  
/units,si !国际单位制 Y-?51g[u  
fe=0.1 !摩擦系数 LIc*tsl  
fkn=10 !法向接触刚度 ;2 \<M6  
fkt=3!切向接触刚度 e1Dj0s?i~K  
ftoln=1 !初始渗透因子 gS0,')w  
fact=2!静动摩擦系数 ZlUFJ*pk  
/prep7 o<Ke3?J\  
et,1,solid45 VdYOm  
mp,ex,1,7.2e10         !桩的弹性模量，根据管壁厚度折算 8~rT  
mp,nuxy,1,0.3           !桩的泊松比 :K5V/-[|V1  
mp,dens,1,2000         !桩的密度 g8B&u u #  
!mp,mu,1,fe             !桩体性质 qRWJ-T:!F  
mp,ex,2,2e8           !土的弹性模量，经验取值 i$2MjFC-  
mp,nuxy,2,0.3          !土的泊松比 047*gn.b  
mp,dens,2,2000          !土的饱和密度   }:tAKO=+  
tb,dp,2 -1dIZy  
tbdata,1,0,30      !粘聚力c为0,摩擦角为30度 1Z=;Uy\  
!创建桩模型 yzODF>KJ  
cylind,0,0.1,-0.1,7.1,0,90 )E_!rR  
!创建土体模型 I%C]>ZZh  
cylind,0.1,0.6,0,7,0,90 _p?I{1O  
!桩体划分网格 y;*My#  
allsel,all $4&8U~Zs  
vsel,s,,,1 AZ]Z,s6  
aslv,s J#_\+G i  
lsla,s 6NP`P jR  
lesize,1,,,6 &7JEb]1C  
lesize,2,,,6 Gf!t< =T   
lesize,3,,,6 ~p0e=u  
lesize,4,,,6 %Gnd"SGs  
lesize,5,,,6 E%KC'TN^D  
lesize,6,,,6 nT(!HD
H  
lesize,7,0.05 1"N/ZKF-x  
lesize,8,0.05 UX+?0K  
lesize,9,0.05 }qTv&Z3$  
type,1 ,(zcl$A[  
mat,1 B&to&|jf  
mshape,0,3d 4h[2C6 \+`  
mshkey,1 sT@u3^>  
vmesh,all 9Vh_XBgP  
vsel,s,,,2 ！土层，划分网格 (gv=P>:  
vplot !;6W!%t.|  
aslv,s i]V F'tG  
lsla,s DWHOSXA4  
lplot 4d,qXSKty  
lesize,10,,,6 S;G"L$&\  
lesize,11,,,6 h:eN>yW  
lesize,12,,,6 75' Ua$  
lesize,13,,,6 w`2_6[,9  
lesize,14,,,6 9iiU,}M`j  
lesize,15,,,6 g5?r9e  
lesize,16,,,6 w?*'vF_2:#  
lesize,17,,,6  ||bA  
lesize,18,0.05 P DNt4=C  
lesize,19,0.05 1j_aH#Fz:  
lesize,20,0.05 V/+H_=|  
lesize,21,0.05 4"`=huQ  
type,1 Sf#\6X<B  
mat,2 \Or]5ogT'  
mshape,0,3d t>fA!K%{  
mshkey,1 6uv'r;U]  
vmesh,all aA!@;rR<yU  
!创建接触对 X:iG[iU*  
!用接触向导建立接触对 8JFnB(3xU  
allsel 1ZGQhjcx  
et,2,170 t;bZc s  
et,3,173 mJU>f-l  
r,1 ajg7xF{l)  
real,1 U!3uaz'  
keyopt,3,9,0    !消除初始渗透 !^rITiy  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 &^"s=g.  
桩侧，建立目标面 gt(X!iN]  
r,1,,,fkn,ftoln, , , jKe$&.q@  
rmore,,,,,,fkt B`t)rBy  
RMORE,,,, ,, >:(6{}b  
RMORE,,,fact,,,   0EF,uRb  
mp,mu,3,fe 'lSnyW{  
!创建接触对 S8rW'}XJ=H  
vsel,s,,,1 %>oT7|x  
asel,s,,,3 89?3,k  
aplot U<#$w{d:  
cm,target,area pRb+'v&_k  
type,2 hA$c.jJr.Z  
real,1 YLr%vnO*NS  
mat,3 )L&n)w  
nsla,s,1 >&4I.nA  
esurf,all y?rK5Yos  
esel,s,type,,2 $NH`Iu9t  
eplot PkZf(=-X  
/psymb,esys,1 0YgFjd 5  
allsel!很重要 6T5A31 Q  
!定义接触面 @8T Vr2uy  
vsel,s,,,2 %`8KG(F^  
vplot qhv4R|)  
asel,s,,,9 fwz5{>ON
]  
aplot il 8A&`%  
cm,contact,area D"1vw<Ak  
type,3 P W0q71  
real,1 j X^&4f  
mat,3 w0F:%:
/  
nsla,s,1 C"V?yDy2~  
esurf,all m7bn%j-{$f  
esel,s,type,,3 X}ey0
)g%  
eplot w%!k?t,*]  
/psymb,esys,1 )4VLm  
allsel |~1rKzZwF  
save [U_Q 2<H  
!侧面施加对称约束
}Etd#">  
allsel mRix0XBI~  
CSYS,0 1`7]C+Pv  
asel,s,loc,x,0 l[ZQ7$kL  
da,all,symm +"*l2
E]5  
asel,s,loc,y,0 ?et0W|^k  
da,all,symm IDL^0:eg<.  
allsel,all OdtbVF~  
!asel,s,,,2!没有约束顶面2 y'i:%n}I  
!da,all,all ?ZD{e|:u  
asel,s,,,6 @9 qzn&A  
da,all,ALL rVc zO+E  
!约束外侧面 Q7OnhGA  
asel,S,,,8 t/EMBfLc  
NSLA,S,1 S:"z<O  
CSYS,1 o)$Q]N##  
NROTAT,ALL Vb"T],N1m  
D,ALL,UX！土侧采用的径向约束。 ^loF#d=s  
ALLSEL Hj-<{#,  
SAVE H-PW(  
finish!接触对定义结束 ;RTrRh0v  
!求解过程 3tx0y
 
/solu s0`]!7D<  
!约束结束! !kjr>
:)x  
!定义荷载步 Q*oA{eZY  
!第一步荷载，重力荷载。 v>yGsJnV'  
acel,0,0,9.8 g6k&c"%IQ(  
nlgeom,on ! 打开大变形选项 v{\n^|=])  
autot,on ! 打开自动时间步 '=@H2T6=  
time,1 Es ZnGuY  
nsub,10 H@Or
X  
outres,all,all iLI.e rm  
allsel 8=u+BDG  
SOLVE 1GyAQHx,  
!第二部荷载,拉动桩 Oa3=+_C~$1  
nsubst,100，500，5 ),W(TL  
outres,all,all I*`=[nR  
time,101 .jrR4@  
nlgeom,on ! 打开大变形选项 a`GN@ 8  
autot,on ! 打开自动时间步 9, sCJ5bb"  
/I`!iK  
nsel,s,loc,z,7.1 RLeSA\di  
D,ALL,UZ,0.1 -hJ>wGI  
lnsrch,on !自动线性搜索 %<bG%V(  
pred,on !打开预测求解器 9Le/'ovq  
allsel JXD?a.vy^q  
SOLVE v\r7.l:hf  
$TH'"XK  

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

这个问题估计你是老大吧

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