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

大家好，我做的是一个桩土接触分析问题，想把一个桩从土中拔出，施加两个荷载步，第一步，施加重力荷载，第二步把桩从中拔出，在运行第二步时出了问题，算不下去，请大师不吝指教 ，不胜感激。命令流如下： `wLmGv+V  
FINI _1$+S0G;  
/CLEAR !$o9:[B  
,Qe`(vU*s  
/filname,pile-soil contact tE=$#  
/title，pile 6XVJ/q
Z  
!/config,nres,2000000 "rQ?2?  
/units,si !国际单位制 z5[Qh<M  
,?GAFgK:  
fe=0.1 !摩擦系数 _8z  
fkn=10 !法向接触刚度 etP`q:6^c  
fkt=3!切向接触刚度 AlgVsE%Va  
ftoln=1 !初始渗透因子 Qm?o^%a  
fact=2!静动摩擦系数 _j_c&  
uC`)?f*I  
/prep7 %o0.8qVJi  
UZc{ Av  
et,1,solid45 rf&nTDaWI  
1 7hXg"B  
mp,ex,1,7.2e10 !桩的弹性模量，根据管壁厚度折算 r3Ih]|FK#  
mp,nuxy,1,0.3 !桩的泊松比 e=Kr>~q=  
mp,dens,1,2000 !桩的密度 lCK:5$ z0  
!mp,mu,1,fe !桩体性质 )#cGePA  
?[O Sy.6  
mp,ex,2,2e8 !土的弹性模量，经验取值 WOGMtT%  
mp,nuxy,2,0.3 !土的泊松比 N%?8Bm~dP  
mp,dens,2,2000 !土的饱和密度 iJu$&u  
tb,dp,2 j7Lw(AJ  
tbdata,1,0,30 !粘聚力c为0,摩擦角为30 f^XfIH_#  
eJ$ {`&J  
!创建桩模型 v#KE"m  
cylind,0,0.1,-0.1,7.1,0,90 3+%L[fW`/  
!创建土体模型 ds QGj&  
cylind,0.1,0.6,0,7,0,90 ';x5 $5k'  
HmEU;UbO-  
!桩体划分网格 ?DJ,YY9P  
allsel,all Hn(L0#Oqy  
vsel,s,,,1 Dr"F5Wbg  
aslv,s en gh3TZC  
lsla,s !1?Nc}T0Q&  
lesize,1,,,6  qZP>h4  
lesize,2,,,6 }
EHmVPe  
lesize,3,,,6 B3Esfk  
lesize,4,,,6 tZY(r {  
lesize,5,,,6 L_7-y92<W  
lesize,6,,,6 o"A%dC_  
lesize,7,0.05 N+ R/ti  
lesize,8,0.05 YJrZ  
lesize,9,0.05 o#&;,9  
type,1 ~ +z'pK~c  
mat,1 nbz?D_  
mshape,0,3d NvlG@^&S  
mshkey,1 kv3jbSK
CT  
vmesh,all !ly]{DTmm  
#9Dixsl*Q  
vsel,s,,,2 ！土层，划分网格 s"]LQM1|  
vplot %!1:BQ,p,i  
aslv,s l4Y}<j\;  
lsla,s .W%{j()op  
lplot t%<d}QuHW  
lesize,10,,,6 ;H8A"$%n~  
lesize,11,,,6 u7<s_M3%N  
lesize,12,,,6 Z$5@r2d)  
lesize,13,,,6 r&ex<(I{  
lesize,14,,,6 xpVYNS{c+|  
lesize,15,,,6 @7';bfsix  
lesize,16,,,6 "ealYve
u  
lesize,17,,,6 #S4{,  
lesize,18,0.05 jW+L0RkX  
lesize,19,0.05 8[;U|SR"  
lesize,20,0.05 I6FglVQ6  
lesize,21,0.05 C"ZCX6p+$  
type,1 v ~%6!Tr  
mat,2 {*%'vVv+  
mshape,0,3d t1 .6+  
mshkey,1 E8Wgm 8  
vmesh,all dKyX70Zy9  
QJ s/0iw  
!创建接触对 X"{s"Mc0G  
!用接触向导建立接触对 asI:J/%+2  
allsel EkqsE$52  
et,2,170 b
-
+iL  
et,3,173 e&u HU8k*  
r,1 ]E)gMf   
real,1 '&}B"1  
keyopt,3,9,0 !消除初始渗透 _A\c 6#  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开
\
^3\_T&6  
9oY%v7  
桩侧，建立目标面 *Ta*0Fr=9|  
r,1,,,fkn,ftoln, , , h^bbU.  
rmore,,,,,,fkt ]baO{pJi  
RMORE,,,, ,, O.*,e  
RMORE,,,fact,,, h
C8'6h  
mp,mu,3,fe [ua{qJ9  
OY6lt.t  
!创建接触对 pFD L5  
vsel,s,,,1 70m}+R(`  
asel,s,,,3 y)!K@  
aplot nlK"2/W  
cm,target,area VHB5
 
type,2 )fR'1_  
real,1 ?+b )=Z  
mat,3 w.Ft-RXA W  
nsla,s,1 *A O/$K@Ma  
esurf,all (Y!@,rKd  
esel,s,type,,2 jO$3>q  
eplot r/fLm8+  
/psymb,esys,1 '8wA+N6Zr7  
allsel!很重要 FbHk6(/
)  
^'C1VQ%  
!定义接触面 3m?3I2k  
vsel,s,,,2 >t O(S  
vplot p}}o#a~V),  
asel,s,,,9 4S@^ym  
aplot SUfl`\O  
cm,contact,area Jpws1~  
type,3 h{kAsd8 G  
real,1 hn-9l1~!h  
mat,3 5+'1 :Sa(i  
nsla,s,1 4?,N;Q  
esurf,all b4s.`%U  
esel,s,type,,3 X5527`?e  
eplot 6<R!`N 6  
/psymb,esys,1 7\K=8G  
allsel $AUC#<*C  
save X)xQKkL0  
,)nO  
!侧面施加对称约束 \:\rkc9LI  
allsel ](@Tbm8  
CSYS,0 %X O97  
asel,s,loc,x,0 a}+|2k_  
da,all,symm +EJIYvkFm  
asel,s,loc,y,0 r+;op_
 
da,all,symm Qhsh{muw(  
Edc< 8-  
allsel,all HkD6aJ:kA!  
!asel,s,,,2!桩顶面2未约束 z'JtH^^Z  
!da,all,all !i
A0u  
asel,s,,,6 EV z>#GC  
da,all,ALL pog*}@OS  
!约束外侧面 7g cr$&+e  
asel,S,,,8 cWP34;NNM  
NSLA,S,1 Qq(/TA0$-  
CSYS,1 P|_>M SO1'  
NROTAT,ALL E=3<F_3W  
D,ALL,UX )VID ;l;4  
ALLSEL {xp/1?Mo*  
SAVE l_y:IY$"  
finish!接触对定义结束 #7+oM8b  
!求解过程 3@*J=LGhKc  
/solu Oj6-  
a1_7plg  
!约束结束! mTT1,|  
!定义荷载步 { ,qm=Xjq  
!第一步荷载，重力荷载。 X<m#:0iD  
acel,0,0,9.8 ?`,<l#sj  
nlgeom,on ! 打开大变形选项 0{j]p^'<  
autot,on ! 打开自动时间步 \m~?mg"#  
time,1 r*fZS$e  
nsub,10 UYn5Pix  
outres,all,all Wqy|Y*$qT  
allsel _$]3&P  
SOLVE #HpF\{{v  
!第二部荷载,拉动桩 $+!/=8R)  
nsubst,100,500,5 @-Tt<pl'L  
outres,all,all ,&rlt+wE  
time,101 D.9qxM"Z>  
nlgeom,on ! 打开大变形选项 *Fb]lM7D  
autot,on ! 打开自动时间步 qI"Xh" c?  
nsel,s,loc,z,7.1 +] >o@  
D,ALL,UZ,0.1 JwL}|o6  
\3LP@;Phn  
lnsrch,on !自动线性搜索 9q5jqFQ  
pred,on !打开预测求解器 s4/4o_[W  
f^b.~jXSR}  
allsel ov*?[Y7|~  
SOLVE

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