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

大家好，我做的是一个桩土接触分析问题，想把一个桩从土中拔出，施加两个荷载步，第一步，施加重力荷载，第二步把桩从中拔出，在运行第二步时出了问题，算不下去，请大师不吝指教 ，不胜感激。命令流如下： o}mhy`}  
FINI S& ,Ju%  
/CLEAR =0s`4Y"+  
2pxWv )0  
/filname,pile-soil contact AF*ni~  
/title，pile Bt*&L[&57  
!/config,nres,2000000 Sr ztTfY  
/units,si !国际单位制 s]F?=yEp  
,);= (r9  
fe=0.1 !摩擦系数 Ypn%[sSOp  
fkn=10 !法向接触刚度 q#Yg0w~  
fkt=3!切向接触刚度 8"?Vcw&  
ftoln=1 !初始渗透因子 2rR@2Vsw2  
fact=2!静动摩擦系数 x%IXwP0  
SmT+L,:D  
/prep7 %S >xSqX  
=9DhO7I'  
et,1,solid45 |z!q r}i  
!~lVv&YO  
mp,ex,1,7.2e10 !桩的弹性模量，根据管壁厚度折算 tf$PaA  
mp,nuxy,1,0.3 !桩的泊松比 Ad}-I%Ie  
mp,dens,1,2000 !桩的密度 fH#F"^A  
!mp,mu,1,fe !桩体性质 a2cx  
IL+#ynC  
mp,ex,2,2e8 !土的弹性模量，经验取值 IKMsY5i  
mp,nuxy,2,0.3 !土的泊松比 @MN}^umx`  
mp,dens,2,2000 !土的饱和密度 a gmeiJT  
tb,dp,2 YiIddQ  
tbdata,1,0,30 !粘聚力c为0,摩擦角为30 $n=W2WJ6f  
EFa{O`_@U  
!创建桩模型 dAYI DE  
cylind,0,0.1,-0.1,7.1,0,90 >~`Y
 
!创建土体模型 g^
s+C Z  
cylind,0.1,0.6,0,7,0,90 /32Ta  
l]a^"4L4`o  
!桩体划分网格 8FO1`%8Oe  
allsel,all &57qjA,8<  
vsel,s,,,1 Ek"YM[  
aslv,s -vv
 
lsla,s nu4Pc  
lesize,1,,,6 ax-=n(   
lesize,2,,,6 lEO?kn.:z  
lesize,3,,,6 SjZd0H0  
lesize,4,,,6 T$gkq>!j<E  
lesize,5,,,6 J>f /u:.  
lesize,6,,,6 + $k07mb\  
lesize,7,0.05 pj!k|F9  
lesize,8,0.05 ,3wI~j=  
lesize,9,0.05 QT`fix{  
type,1 [WcS[](ob  
mat,1 3mpP|b"  
mshape,0,3d '4af ],  
mshkey,1 c_~XL^B@  
vmesh,all
EG F:xl  
er(8}]X8Q  
vsel,s,,,2 ！土层，划分网格 ,$HHa
oog  
vplot -(6eVI  
aslv,s X tJswxw`K  
lsla,s !Bag}|#  
lplot n (OjjRm  
lesize,10,,,6 jFgZ}Xp  
lesize,11,,,6 9U!J
K3d  
lesize,12,,,6 ,)Yao;Cvd  
lesize,13,,,6 y' 2<qj  
lesize,14,,,6 6o$Z0mG  
lesize,15,,,6 ^it4z gx@  
lesize,16,,,6 `ffj8U  
lesize,17,,,6 x*:n4FZ7b  
lesize,18,0.05 +nU',E  
lesize,19,0.05 -s{R/6:  
lesize,20,0.05 wG_4$kyj  
lesize,21,0.05 cB{%u '  
type,1 Y3D3.T6Q  
mat,2 ym
HKcQ  
mshape,0,3d a a4$'8s  
mshkey,1 #1MKEfv(~  
vmesh,all ;F|jG}M
"  
BeAk21xb  
!创建接触对 3QH(4N  
!用接触向导建立接触对 i"/r)>"b  
allsel 2K Pqu:lv  
et,2,170 o\d |CE;>  
et,3,173 {>S4#^@}  
r,1  1U  
real,1 QI.t&sCh5  
keyopt,3,9,0 !消除初始渗透 VZ;@S3TS  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 x 96}#0'  
Lg8]dBXu  
桩侧，建立目标面 Ec+22X  
r,1,,,fkn,ftoln, , , #sL
/y  
rmore,,,,,,fkt -H4PRCDH  
RMORE,,,, ,, n>n"{!  
RMORE,,,fact,,,
^Wf S\M`  
mp,mu,3,fe vXc!Zg~  
av$  
!创建接触对 O(o
dNQy~  
vsel,s,,,1 qv.n99?]  
asel,s,,,3 ZH% we  
aplot hUBF/4s\  
cm,target,area 8*vFdoE_oO  
type,2 :|=- (z  
real,1 L+7*NaPY*  
mat,3
m.!wsw  
nsla,s,1 @$Yb#$/  
esurf,all !f@XDW&R  
esel,s,type,,2 m!;G/s*  
eplot 'w:bs!  
/psymb,esys,1 c+H)ed>  
allsel!很重要 !h?=Wv ==]  
&h\7^=s.  
!定义接触面 2av*o~|J*:  
vsel,s,,,2 `6UtxJSx  
vplot L -YNz0A  
asel,s,,,9 5 1"8Py  
aplot E_]k>bf\  
cm,contact,area h+5@I%WX  
type,3 / xs9.w8-  
real,1 %iX/y  
mat,3 %TDXF_.[  
nsla,s,1 A
=0@UqM  
esurf,all 7r['  
esel,s,type,,3 \x(.d.l/  
eplot 9Eu #lV  
/psymb,esys,1 /-><k,mL?  
allsel zPHy2H$28  
save /e .D/;]  
tGC2 ^a#~  
!侧面施加对称约束 d[S#Duz<&  
allsel YQ+^  
CSYS,0 y)tYSTJK  
asel,s,loc,x,0 V'C-'Ythwf  
da,all,symm Ajg\aof0{  
asel,s,loc,y,0 (U_wp's  
da,all,symm LfHzT<)|  
?|lIXz  
allsel,all %2}C'MqS  
!asel,s,,,2!桩顶面2未约束 bb}?h]a  
!da,all,all 8ZW?|-i  
asel,s,,,6 $HxS:3D%D  
da,all,ALL i+1Qf  
!约束外侧面 &>g~-s  
asel,S,,,8 zA,/@/'(  
NSLA,S,1 S*xhX1yUi  
CSYS,1 )@_ugW-j
 
NROTAT,ALL O2% `2h  
D,ALL,UX A=y"x$%-_  
ALLSEL dqxd3,Z  
SAVE
nc&Jmo7  
finish!接触对定义结束 MUB37  
!求解过程 DFc [z"[  
/solu NocFvF7\  
nPOO3!<{  
!约束结束! |AC6sfA+  
!定义荷载步 ^aHh{
BQ%  
!第一步荷载，重力荷载。 uOre,AQR  
acel,0,0,9.8 Y1WHy*s?  
nlgeom,on ! 打开大变形选项 pqH4w(;  
autot,on ! 打开自动时间步 6rT4iC3Q{  
time,1 wFb@1ae\  
nsub,10 GG+5/hU  
outres,all,all ;NP[_2|-,  
allsel Fg4@On[,i  
SOLVE vd(dNu&,<  
!第二部荷载,拉动桩 Zih ?Bm  
nsubst,100,500,5 :3 y_mf>  
outres,all,all ,Bf(r  
time,101 "wINBya'M  
nlgeom,on ! 打开大变形选项 LipxAE?O  
autot,on ! 打开自动时间步 jtCZfFD?
 
nsel,s,loc,z,7.1 `|?]CkP  
D,ALL,UZ,0.1 8YKQItK  
) $#(ZL^m  
lnsrch,on !自动线性搜索 r0{]5JZt/  
pred,on !打开预测求解器 Z/7dg-$?'0  
yd*3)6=  
allsel 4.'JLArw  
SOLVE

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