[软件应用交流]请教桩土接触分析，桩从土中拔出模拟 [复制链接]

大家好，我做的是一个桩土接触分析问题，想把一个桩从土中拔出，施加两个荷载步，第一步，施加重力荷载，第二步把桩从中拔出，在运行第二步时出了问题，算不下去，请大师不吝指教 ，不胜感激。命令流如下： ^uV=|1<%  
FINI 'vP"&lrn  
/CLEAR Pg''>6w>  
OsrHA  
/filname,pile-soil contact }!`_Bz:  
/title，pile ]B2%\}c  
!/config,nres,2000000 "hbCP4  
/units,si !国际单位制 k#oe:u
`<  
#n_gry!5  
fe=0.1 !摩擦系数 >WcOY7  
fkn=10 !法向接触刚度 Jyz*W!kI  
fkt=3!切向接触刚度 "9^OT  
ftoln=1 !初始渗透因子 q*^m8  
fact=2!静动摩擦系数 x+Ws lN2a  
T4JG5  
/prep7 CVAX?c{
 
P9W!xvV`w  
et,1,solid45 N 4!18{/2  
A)5;ae  
mp,ex,1,7.2e10 !桩的弹性模量，根据管壁厚度折算 Q?g#?z&Pu\  
mp,nuxy,1,0.3 !桩的泊松比 +I<^w)  
mp,dens,1,2000 !桩的密度 _;!$1lM[  
!mp,mu,1,fe !桩体性质 "Dt: 8Nf^  
['%$vnS5S  
mp,ex,2,2e8 !土的弹性模量，经验取值 Q"Pl)Q\  
mp,nuxy,2,0.3 !土的泊松比 pXhN?joe  
mp,dens,2,2000 !土的饱和密度 )w_hbU_Pb&  
tb,dp,2 )8Jf
Bz
R  
tbdata,1,0,30 !粘聚力c为0,摩擦角为30 9OS~;9YR  
qlNB\~HCe  
!创建桩模型 zMg(\8  
cylind,0,0.1,-0.1,7.1,0,90 k9*6`w  
!创建土体模型 K_Q-9j  
cylind,0.1,0.6,0,7,0,90 )Y](Mj!D  
L!c.1Rf_  
!桩体划分网格 EK%J%NY  
allsel,all \z8j6 h  
vsel,s,,,1 /S #Z.T~~  
aslv,s :'$V7LZ5  
lsla,s Gf->N `N  
lesize,1,,,6 M669G;w(K  
lesize,2,,,6 8 U<$u,WS  
lesize,3,,,6 `'vNHY  
lesize,4,,,6 \dHdL\f  
lesize,5,,,6 kM;}$*?  
lesize,6,,,6 sJ>JHv  
lesize,7,0.05 a !yBEpMo  
lesize,8,0.05 .gJv})Vi  
lesize,9,0.05 hU~up a<dD  
type,1 Xt%y>'.  
mat,1 ^&z3zFTp  
mshape,0,3d !^x;4@Ejm  
mshkey,1 N0V`xrS  
vmesh,all d(_;@%p1X  
$ dR@Q?_{  
vsel,s,,,2 ！土层，划分网格 j9d^8)O,  
vplot INRP@Cp1  
aslv,s 8o  SL3  
lsla,s PiVp(; rtQ  
lplot '3uN]-A>D  
lesize,10,,,6 O:{N5+HVG  
lesize,11,,,6 =j!nt8]8  
lesize,12,,,6 _, r6t  
lesize,13,,,6 cD5N'3  
lesize,14,,,6 tt03gU`  
lesize,15,,,6 O4g2s8k  
lesize,16,,,6 KlGmO;k  
lesize,17,,,6 ww5UQs2sn  
lesize,18,0.05  84g8$~M  
lesize,19,0.05 sDZ<XA  
lesize,20,0.05 B
GrV,h^  
lesize,21,0.05 "Q.KBX v/  
type,1 UfNcI[xr  
mat,2 n|'}W+  
mshape,0,3d Njmb{L]Cps  
mshkey,1 q&nEodv>+  
vmesh,all ?Z2_y-  
Ywo
=w:'  
!创建接触对 cl{kCSZo.z  
!用接触向导建立接触对 MFtC2*
 
allsel IQ $/|b/  
et,2,170 n_/;j$h  
et,3,173 }? :T*CJ  
r,1 5{|tE!  
real,1 g@Z7f y7  
keyopt,3,9,0 !消除初始渗透 ,GYK3+}Z  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 T!2gOe  
[!S%nYs&8L  
桩侧，建立目标面 i4dy0j
fN  
r,1,,,fkn,ftoln, , , ($X2SIZh  
rmore,,,,,,fkt [KW9J}]  
RMORE,,,, ,, ;[q>  
RMORE,,,fact,,,
nkO
4~p  
mp,mu,3,fe +'"NKZ.>TT  

?58*#'r  
!创建接触对 Ge=|RAw3  
vsel,s,,,1 iGw\A!}w\  
asel,s,,,3 )~{8C:  
aplot ,opS)C$  
cm,target,area *?x[pqGq  
type,2 nj<nW5[  
real,1 VD90JU]X<  
mat,3 G Tz>}
@W  
nsla,s,1 Srom@c  
esurf,all mcb|N_#n/  
esel,s,type,,2 \B Uno6  
eplot cR6Rb[9 N  
/psymb,esys,1 !F
08F>@D  
allsel!很重要 qir8RPW  
_x]q`[Dih  
!定义接触面 ,c&%/
"i:w  
vsel,s,,,2 Yc-gJI*1  
vplot O|mWQp^?q  
asel,s,,,9 C>*1f|<  
aplot 87-z=>IU  
cm,contact,area Blox~=cW  
type,3 w gkY\Q  
real,1 WxJV zHtR  
mat,3 ~(-df>  
nsla,s,1 El^V[s'3  
esurf,all mum4Uj  
esel,s,type,,3 EG J/r  
eplot cq4sgQ?sW  
/psymb,esys,1 AkEt=v
I  
allsel b~C^cM  
save Ewa/6=]LA  
YfUo=ku  
!侧面施加对称约束 &`2$,zX#  
allsel ZPlY]e  
CSYS,0 c9ea%7o{0a  
asel,s,loc,x,0 ,CP&o  
da,all,symm Vif)e4{Pn  
asel,s,loc,y,0 rebWXz7  
da,all,symm ~93#L_V_O  
!a7YM4D  
allsel,all q!
as~{!  
!asel,s,,,2!桩顶面2未约束 _ YcIGOL  
!da,all,all C,) e7  
asel,s,,,6 M=sGPPj  
da,all,ALL e8U6D+jY  
!约束外侧面  (2dkmn  
asel,S,,,8 ^5Ob(FvU  
NSLA,S,1 o9+fAH`D  
CSYS,1 4vMjVbr  
NROTAT,ALL We@wN:  
D,ALL,UX /_V4gwb}|-  
ALLSEL Jl fIYf~  
SAVE `OHdo$Y9  
finish!接触对定义结束 GDwijZw  
!求解过程 )5ev4Qf  
/solu h%ba!  
<y<   
!约束结束! :OD-L)Or  
!定义荷载步 KSR'X0'  
!第一步荷载，重力荷载。 h/NI5  
acel,0,0,9.8 axM(3k.n  
nlgeom,on ! 打开大变形选项 %g7B*AX]  
autot,on ! 打开自动时间步 b" kL)DL1L  
time,1 |o#pd\  
nsub,10 IQyw>_~]  
outres,all,all -uhg7N[3  
allsel m/"}Y]
n!  
SOLVE [u;>b?[{  
!第二部荷载,拉动桩 LrhQG  
nsubst,100,500,5 o(@^V!}V  
outres,all,all /&o<kY  
time,101 ~kOXMLRg  
nlgeom,on ! 打开大变形选项
,D93A  
autot,on ! 打开自动时间步 T1A/>\Ns  
nsel,s,loc,z,7.1 0O|l7mCr%I  
D,ALL,UZ,0.1 t $u.  
F @uOXNz)  
lnsrch,on !自动线性搜索 Io4Ss1="  
pred,on !打开预测求解器 NI2-*G_M  
Y.#:l<  
allsel uX8G<7O^  
SOLVE

很详细啊！看来作者坐了不少工作。