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

我做的是一个桩土分析问题，想计算一下桩土相对滑动时的摩擦力，采用的是1/4三维实体模型，计划分为两个荷载步，第一步利用土体自重挤住桩，第二步把桩从土中拔出。 {: /}NpA$  
第一步可以算，第二步却不行，命令流如下，请大师们不吝指教，不胜感激。 X'ag)|5ot  
FINI ?<!|  
/CLEAR IGl9g_18  
/filname,pile-soil contact kiEa<-]  
/title，pile M`_0C38  
/units,si !国际单位制 w)f#V s  
fe=0.1 !摩擦系数 J.a]K[ci  
fkn=10 !法向接触刚度 :#Wd~~d  
fkt=3!切向接触刚度 x2xRBkRg=  
ftoln=1 !初始渗透因子 )=+|i3]U  
fact=2!静动摩擦系数 sJZiI}Xc  
/prep7 F9PxSk_\9  
et,1,solid45 G|Ti4_w  
mp,ex,1,7.2e10         !桩的弹性模量，根据管壁厚度折算 V~GDPJ+  
mp,nuxy,1,0.3           !桩的泊松比 i-1op> Y  
mp,dens,1,2000         !桩的密度 3";q[&F9y  
!mp,mu,1,fe             !桩体性质 `5*}p#G  
mp,ex,2,2e8           !土的弹性模量，经验取值 MgZ/(X E  
mp,nuxy,2,0.3          !土的泊松比 =_CzH(=f#  
mp,dens,2,2000          !土的饱和密度   4#D,?eA7  
tb,dp,2 rq{$,/6.  
tbdata,1,0,30      !粘聚力c为0,摩擦角为30度 dtDFoETz  
!创建桩模型 }BEB1Q}L  
cylind,0,0.1,-0.1,7.1,0,90 /ZX}Nc g  
!创建土体模型 Wtnfa{gP%  
cylind,0.1,0.6,0,7,0,90 6ujWNf  
!桩体划分网格 F?0Ykjh3  
allsel,all I9^x,F"E]  
vsel,s,,,1 OUnA;_  
aslv,s &oNAv-m^GD  
lsla,s vx =&QavL  
lesize,1,,,6 Z,gk|M3.  
lesize,2,,,6 #!=tDc &  
lesize,3,,,6 -"x$ZnHU  
lesize,4,,,6 VbYdZCC  
lesize,5,,,6 ]Wup/o  
lesize,6,,,6 )%TmAaj9d  
lesize,7,0.05 W/N7vAx X  
lesize,8,0.05 F,kZU$  
lesize,9,0.05 5x
i
EPh  
type,1 F59 TZI  
mat,1 ).O)p9  
mshape,0,3d W9&=xs6  
mshkey,1 ! nx{ X  
vmesh,all }e1ZbmW  
vsel,s,,,2 ！土层，划分网格 0GLM(JmK  
vplot NEs:},)o  
aslv,s ?ub35NLa  
lsla,s xT8?&Bx  
lplot P \I|,  
lesize,10,,,6 iZmcI;?u  
lesize,11,,,6 5P bW[  
lesize,12,,,6 >P(.:_^p  
lesize,13,,,6 %d@z39-;  
lesize,14,,,6 kh<2BOV  
lesize,15,,,6 [),ige  
lesize,16,,,6 F4QVAOM]U  
lesize,17,,,6 y_)FA"IkE  
lesize,18,0.05 :jf3HG  
lesize,19,0.05
Ry&6p>-  
lesize,20,0.05 4euO1=  
lesize,21,0.05 Wwo0%<2y  
type,1 %#+Hl0,Tt  
mat,2 e-;}3
66}  
mshape,0,3d 8ag!K*\V<  
mshkey,1 JF]JOI6.e  
vmesh,all [E_9V%^  
!创建接触对 sOY:e/_F  
!用接触向导建立接触对 lE;!TQj:X  
allsel +@UV?"d  
et,2,170 bA 2pbjg=  
et,3,173 42{~Lhxt  
r,1 9r9NxKuAO  
real,1 gYj'(jB  
keyopt,3,9,0    !消除初始渗透 Z+SRXKQ  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 xdPx{"C 3  
桩侧，建立目标面 / {%%"j  
r,1,,,fkn,ftoln, , , DU^loB+  
rmore,,,,,,fkt y =@N|f!  
RMORE,,,, ,, P?<y%c<  
RMORE,,,fact,,,   4H/OBR  
mp,mu,3,fe 7<4qQ.deE  
!创建接触对 SbZ6t$"  
vsel,s,,,1 XW/o<[91  
asel,s,,,3 [g,}gyeS(  
aplot f);FoVa6  
cm,target,area \V:^h[ad  
type,2
vO=fP_  
real,1 z:O8Ls^\T  
mat,3 cQ|NJ_F{1  
nsla,s,1 )7@0[>  
esurf,all XppOU  
esel,s,type,,2 )oZ dj`  
eplot P>T"cv  
/psymb,esys,1 "@kaHIf[  
allsel!很重要 NK+o1  
!定义接触面 f$( e\++  
vsel,s,,,2 KvSG;  
vplot 6!o1XQr=Z  
asel,s,,,9 4i bc  
aplot hTkyz la  
cm,contact,area xw%0>K[  
type,3 'w/hw'
F6  
real,1 7)m9"InDI  
mat,3 ]9-\~Mwh  
nsla,s,1 1C.VnzRnJ  
esurf,all 2oW"'43X  
esel,s,type,,3 M|-)GvR$J  
eplot -j(6;9"7]|  
/psymb,esys,1 ICCc./l|  
allsel A&{Nh` q
 
save M5B# TAybC  
!侧面施加对称约束 #ob/p#k  
allsel zs;JJk^  
CSYS,0 G}*hM$F  
asel,s,loc,x,0 a*;b^Ze`v  
da,all,symm :k"]5>(^  
asel,s,loc,y,0 ?2a$*(  
da,all,symm
Dq xs+
 
allsel,all k)u[0}  
!asel,s,,,2!没有约束顶面2 *j=% #  
!da,all,all =Qq+4F)MD  
asel,s,,,6 ;S{(]K7i  
da,all,ALL Xj*Wu_  
!约束外侧面 '-6~tWC~7  
asel,S,,,8 hZ3bVi)L\  
NSLA,S,1 :Tc^y%b0  
CSYS,1 E`q_bn  
NROTAT,ALL }u|q0>^8  
D,ALL,UX！土侧采用的径向约束。 YIE<pX4Q7)  
ALLSEL $]1=\I  
SAVE 9uY'E'm*  
finish!接触对定义结束 AW .F3hN)  
!求解过程 <3iMRe  
/solu 0:+E-^X  
!约束结束! 0(Ij%Wi,  
!定义荷载步 DIvHvFss  
!第一步荷载，重力荷载。 )jj0^f1!j  
acel,0,0,9.8 ?%86/N>  
nlgeom,on ! 打开大变形选项 J,G lIv.A  
autot,on ! 打开自动时间步 w!CNRtM:~  
time,1 QJNFA}*>  
nsub,10 6zkaOA46V  
outres,all,all mOSv9w#,  
allsel z#N@ 0R  
SOLVE 4Hg9N}  
!第二部荷载,拉动桩 3T 9j@N77  
nsubst,100，500，5 NA*&#X#~  
outres,all,all ^8tEach  
time,101 l6B@qYLZ  
nlgeom,on ! 打开大变形选项 C~[,z.FvO  
autot,on ! 打开自动时间步 <{pz<io)  
)"LJ hLg  
nsel,s,loc,z,7.1 ex|F|0k4}  
D,ALL,UZ,0.1 wr4:Go`  
lnsrch,on !自动线性搜索 ijcm2FJcG  
pred,on !打开预测求解器 NI5``BwpO  
allsel N [@?gFtT  
SOLVE fM}#ON>Z  
)[ ,A_3E  

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

这个问题估计你是老大吧

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