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

我做的是一个桩土分析问题，想计算一下桩土相对滑动时的摩擦力，采用的是1/4三维实体模型，计划分为两个荷载步，第一步利用土体自重挤住桩，第二步把桩从土中拔出。 eXD~L&s[  
第一步可以算，第二步却不行，命令流如下，请大师们不吝指教，不胜感激。 3~7!=s\v  
FINI YXZP-=fB>i  
/CLEAR k6J&4?xZ
  
/filname,pile-soil contact nW^h +  
/title，pile '&hd^9]Lo  
/units,si !国际单位制 &ds+9A  
fe=0.1 !摩擦系数 Mtv{37k~  
fkn=10 !法向接触刚度 XI|k,Ko<  
fkt=3!切向接触刚度 zc=G4F01  
ftoln=1 !初始渗透因子 \NMqlxp2  
fact=2!静动摩擦系数 rM^2yr7H  
/prep7 + kT ]qH  
et,1,solid45 M 87CP=yc  
mp,ex,1,7.2e10         !桩的弹性模量，根据管壁厚度折算 1pg#@h[|t  
mp,nuxy,1,0.3           !桩的泊松比 I]i( B+D  
mp,dens,1,2000         !桩的密度 JV@G9PT  
!mp,mu,1,fe             !桩体性质 \+nV~Pi"A  
mp,ex,2,2e8           !土的弹性模量，经验取值 c-5AI{%bl6  
mp,nuxy,2,0.3          !土的泊松比 9r+'DX?>  
mp,dens,2,2000          !土的饱和密度   |pv$],&&:  
tb,dp,2 71%$&6  
tbdata,1,0,30      !粘聚力c为0,摩擦角为30度
:.-z!  
!创建桩模型 ]{#=WTp]  
cylind,0,0.1,-0.1,7.1,0,90 9) ,|h  
!创建土体模型 9od*N$  
cylind,0.1,0.6,0,7,0,90 u8.Tu7~  
!桩体划分网格 NA/`LaJ  
allsel,all pDYJLh-C  
vsel,s,,,1 {d*OJ/4  
aslv,s 7M4J{}9  
lsla,s |&{S ~^$  
lesize,1,,,6 7l$ u.[  
lesize,2,,,6 {J izCUo_'  
lesize,3,,,6 =cZ24I  
lesize,4,,,6 4ryG_p52l  
lesize,5,,,6 eUZk|be  
lesize,6,,,6 qPD(D{,f$  
lesize,7,0.05 g)^s+Y  
lesize,8,0.05 PuREqa\_[  
lesize,9,0.05 W90!*1  
type,1 i0$*):b  
mat,1 ZLxa|R7  
mshape,0,3d o+}k$i!6  
mshkey,1 &z1U0uk  
vmesh,all >V?0#f45@  
vsel,s,,,2 ！土层，划分网格 
H7.l)'  
vplot Q&vdBO/  
aslv,s /+JnEFf  
lsla,s
/{`"X_.o  
lplot jk&xzJH.  
lesize,10,,,6 ^Z#G_%\Y:  
lesize,11,,,6 \u{4=-C.  
lesize,12,,,6 ]b;m~|9  
lesize,13,,,6 fn,hP_  
lesize,14,,,6 _qf39fM;\  
lesize,15,,,6 {icTfPR4E  
lesize,16,,,6 *|g[Mn  
lesize,17,,,6 (m,H
 5  
lesize,18,0.05 ;veD?|  
lesize,19,0.05 7`;f<QNo  
lesize,20,0.05 o ?vGI=
  
lesize,21,0.05 pXl[I;  
type,1 t4#gW$+^?H  
mat,2 KS93v9|  
mshape,0,3d WYTqQqQk  
mshkey,1 P9'` 2c  
vmesh,all Nz.X$zUmY  
!创建接触对 ur8+k4]\"  
!用接触向导建立接触对 i(;`x  
allsel F/BR#J1  
et,2,170 }R#W<4:  
et,3,173 QHZ",1F  
r,1 xmcZN3 ){+  
real,1 ceE]^X;p  
keyopt,3,9,0    !消除初始渗透 p9?kJKN  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 [Cv./hE
Qi  
桩侧，建立目标面 $ch`.$wx  
r,1,,,fkn,ftoln, , , 9/{zS3h3  
rmore,,,,,,fkt >*h+N? m  
RMORE,,,, ,, $U]T8;5Q  
RMORE,,,fact,,,   CUI+@|]%  
mp,mu,3,fe
wxo  
!创建接触对 z``wqK  
vsel,s,,,1 t kj  
asel,s,,,3 N$]er'`  
aplot o=($'(1  
cm,target,area QcQ%A%VIV  
type,2 zq{UkoME  
real,1 !`I@Rk]`c  
mat,3 i m;6$3  
nsla,s,1 `L`+`B  
esurf,all f; |fS~  
esel,s,type,,2 vHpw?
(]
 
eplot N5=BjXSAg  
/psymb,esys,1 *C5`LgeX  
allsel!很重要 d6QrB"J`  
!定义接触面 h9ScN(|0y  
vsel,s,,,2 e9h@G#  
vplot 0,t%us/q  
asel,s,,,9 &R<aRE:+R  
aplot rvd$4l^  
cm,contact,area h:362&?]  
type,3 50N4J  
real,1 -Y/i h(I^  
mat,3 +n;nvf}(  
nsla,s,1 lfc&#Gi3  
esurf,all {s0%XG1$  
esel,s,type,,3 6Kpq~o   
eplot ,sAAV%">
 
/psymb,esys,1 KN|<yF   
allsel JyMk @Y  
save zg<-%r'$  
!侧面施加对称约束 Q p>b  
allsel f?sm~PwC-  
CSYS,0 v&YeQC>  
asel,s,loc,x,0 #%p44%W  
da,all,symm ^b+>r  
asel,s,loc,y,0 g;63$_
<  
da,all,symm `]eJF|"  
allsel,all M-"j8:en  
!asel,s,,,2!没有约束顶面2 p5fr}#en  
!da,all,all AYA{_^#+3  
asel,s,,,6 _ Ncbo#G  
da,all,ALL {\l  
!约束外侧面
`3rwqcxA  
asel,S,,,8 .Az36wD  
NSLA,S,1 l:V R8g[  
CSYS,1 %f1%9YH  
NROTAT,ALL !RmVb}m  
D,ALL,UX！土侧采用的径向约束。 njy2pDC@  
ALLSEL )Fv.eIBY  
SAVE c 3QgX4vq  
finish!接触对定义结束 +H2Jhgi  
!求解过程 ^U_T<x8{  
/solu :*''ci  
!约束结束! \#9LwC"8;  
!定义荷载步 9b+jT{Tg  
!第一步荷载，重力荷载。 `-L?x2)U  
acel,0,0,9.8 t_3j_`  
nlgeom,on ! 打开大变形选项 | ;tH?E  
autot,on ! 打开自动时间步 2^WJ1: A  
time,1 a+^`+p/5  
nsub,10 2!/Kt O)i^  
outres,all,all C 5.3[  
allsel %F>~2g?$  
SOLVE !2L
X+*;  
!第二部荷载,拉动桩 hC:'L9Y  
nsubst,100，500，5 Kq")\Ha,f  
outres,all,all 6$k#B ~~  
time,101 m+Ye`]  
nlgeom,on ! 打开大变形选项 *DF3juf~
 
autot,on ! 打开自动时间步 {[oNUzcd  
C Bkoky9&  
nsel,s,loc,z,7.1 JY0aE  
D,ALL,UZ,0.1 a$My6Qa#  
lnsrch,on !自动线性搜索 uks75W!}U  
pred,on !打开预测求解器 !/is+ xp  
allsel y7,fFUKl  
SOLVE Pp8S\%z~h  
cfLLFPhv)  

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

这个问题估计你是老大吧

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