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

大家好，我做的是一个桩土接触分析问题，想把一个桩从土中拔出，施加两个荷载步，第一步，施加重力荷载，第二步把桩从中拔出，在运行第二步时出了问题，算不下去，请大师不吝指教 ，不胜感激。命令流如下： ![jP)WgF  
FINI nJC}wh2d#  
/CLEAR =Ybu_>  
xbA2R4|  
/filname,pile-soil contact `rQl{$9IC  
/title，pile 3|3lUU\I  
!/config,nres,2000000 ? GW3E  
/units,si !国际单位制 FgaBwd^W  
m!(K  
fe=0.1 !摩擦系数 jX@9849@  
fkn=10 !法向接触刚度 1t'\!  
fkt=3!切向接触刚度 CB)#; |aDB  
ftoln=1 !初始渗透因子 "rJL ^ \r  
fact=2!静动摩擦系数 YdaJ&  
mj~CCokF{?  
/prep7 Vtri"G8 aB  
Y [S^&pF  
et,1,solid45 !I&Sy]G  
inrL'z
 
mp,ex,1,7.2e10 !桩的弹性模量，根据管壁厚度折算 (^\i(cfu6Q  
mp,nuxy,1,0.3 !桩的泊松比 8^T$6A[b  
mp,dens,1,2000 !桩的密度 fsu"Lc  
!mp,mu,1,fe !桩体性质 {eV_+@dT  
w+H=Xh4t  
mp,ex,2,2e8 !土的弹性模量，经验取值 u1<kdTxA N  
mp,nuxy,2,0.3 !土的泊松比  f;a6ux#  
mp,dens,2,2000 !土的饱和密度 ;_*F [ }w  
tb,dp,2 U5=J;[w}N  
tbdata,1,0,30 !粘聚力c为0,摩擦角为30 K)OlCpHc  
cKAl 0_[f"  
!创建桩模型 %Kp}Wo6  
cylind,0,0.1,-0.1,7.1,0,90 na)ceN2h  
!创建土体模型 (FHh,y~v  
cylind,0.1,0.6,0,7,0,90 T94$}- 5/)  
dI|/Xm>  
!桩体划分网格 1qF.0  
allsel,all d0az#Yg!  
vsel,s,,,1 XwMC/]lK<  
aslv,s AQZ\
Kcr  
lsla,s I*"]!z1  
lesize,1,,,6 } q(0uzaG  
lesize,2,,,6 ;'}xD5]  
lesize,3,,,6 =QRZ(2Wq  
lesize,4,,,6 B;Vl+}R  
lesize,5,,,6 ~qH@Kz\%  
lesize,6,,,6 )=@ XF0  
lesize,7,0.05 ^\%%9jY  
lesize,8,0.05 +a;j>hh  
lesize,9,0.05 ^bGi_YC  
type,1 i|Wn*~yFOO  
mat,1 e#^by(1@}  
mshape,0,3d RJM(+5xQ|  
mshkey,1 plAt +*&  
vmesh,all /
2
N%Z  
cPSu!u}D  
vsel,s,,,2 ！土层，划分网格 eKOTxv{  
vplot EbH
eP  
aslv,s 
hRu%> =7  
lsla,s 2$=HDwv  
lplot L_|Y_=r."  
lesize,10,,,6 3WS %H17  
lesize,11,,,6 +/
tD$  
lesize,12,,,6 ~JiA  
lesize,13,,,6 GS%Dn^l  
lesize,14,,,6 $'Z!Y;Ue  
lesize,15,,,6 k@'.d)y0`  
lesize,16,,,6 0M p>X  
lesize,17,,,6 MiRB*eA  
lesize,18,0.05 ]g
ZjV  
lesize,19,0.05 lvlH5Fc  
lesize,20,0.05 KNhH4K2iP8  
lesize,21,0.05 %iv'/B8  
type,1 DGnswN%n1  
mat,2 EzaOg|  
mshape,0,3d lLv0lf  
mshkey,1 uP
Pe"$  
vmesh,all RZjR d   
3}Pa,uN  
!创建接触对 =%p{"<  
!用接触向导建立接触对 N8TO"`wdbs  
allsel ?~Des"F6)1  
et,2,170 8:bNFgJD  
et,3,173 sEa:p:!  
r,1 +FR"Gt$g  
real,1 T}*'9TB  
keyopt,3,9,0 !消除初始渗透 <[bDNe["?  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 0a5P@;"a  
I\_R& v  
桩侧，建立目标面 'I2)-=ZL6  
r,1,,,fkn,ftoln, , , 19_F\32  
rmore,,,,,,fkt [A47OR  
RMORE,,,, ,, $oEDyC  
RMORE,,,fact,,, C Qebb:y  
mp,mu,3,fe FKP^f\!M  
A 3l1$t#w  
!创建接触对 {U9jA
_XX  
vsel,s,,,1 E7D DMU  
asel,s,,,3 RWE%?`  
aplot &knnWm"  
cm,target,area FW4 hqgE@  
type,2 {ig@Iy~DT  
real,1 PZ
s  
mat,3 ?W(f%/B#  
nsla,s,1 x$E l7=.  
esurf,all tznT*EQr  
esel,s,type,,2 Zl)|x%z  
eplot moGbBkO  
/psymb,esys,1 {)M4h?.2  
allsel!很重要 NKRXY~zHh  
5V0=-K  
!定义接触面 ;&!l2UB%  
vsel,s,,,2 ~oI49Q&{  
vplot Dcf`+?3  
asel,s,,,9 uH#NJoRO  
aplot KME #5=~  
cm,contact,area v3p..A~XZ.  
type,3 k)J7) L  
real,1 t65!2G"<  
mat,3 Smp+}-3O  
nsla,s,1 a5iMCmL
+  
esurf,all b#M<b.R)  
esel,s,type,,3 m`|Z1CT  
eplot 1NTe@r!y  
/psymb,esys,1 U7W ct %  
allsel y.Py>GJJ1S  
save Vswi /(  
P|?z1JUd  
!侧面施加对称约束 >Et?7@   
allsel H_Yy.yi  
CSYS,0 _F,OS<>  
asel,s,loc,x,0 TsD >m  
da,all,symm v7-'H/d.  
asel,s,loc,y,0 <j>@Fg#q  
da,all,symm d3\8BKp  
gDv]n^&  
allsel,all ;WhB2/5v  
!asel,s,,,2!桩顶面2未约束 d7&PbITN  
!da,all,all 4Y]`> ;
w  
asel,s,,,6 =P!Vi6[gF~  
da,all,ALL ^V:YNUqp#  
!约束外侧面 `'>>[*06:a  
asel,S,,,8 La!PGZ{  
NSLA,S,1 #df43_u  
CSYS,1 \=@}(<4  
NROTAT,ALL /X?Nv^Hy  
D,ALL,UX Wi[Y@  
ALLSEL Z`W.(gua  
SAVE ;KhYh S(q  
finish!接触对定义结束 buoz La  
!求解过程 .q=X58tHu  
/solu b7
n~z1$  
`XnFc*L 1  
!约束结束! }8svd#S+  
!定义荷载步 17GyE=Uu  
!第一步荷载，重力荷载。 oTL "]3`'  
acel,0,0,9.8 ,uw&)A  
nlgeom,on ! 打开大变形选项 y|aWUX/a  
autot,on ! 打开自动时间步 yDKX,  
time,1 ]:njP3r  
nsub,10 0MOAd!N  
outres,all,all ; ^$RG
  
allsel )sQbDA|p  
SOLVE Ub"\LUu  
!第二部荷载,拉动桩 "n\!y~:  
nsubst,100,500,5 &.}zZ/  
outres,all,all n5b
N/  
time,101 H\S,^)drJ?  
nlgeom,on ! 打开大变形选项 &`9lIVB,K  
autot,on ! 打开自动时间步 fVkl-<?x  
nsel,s,loc,z,7.1 $$4% .J26Z  
D,ALL,UZ,0.1 kO4C^pl"v  
X>CYKRtb  
lnsrch,on !自动线性搜索 DFiexOb  
pred,on !打开预测求解器 (X8N?t
J  
L]VK9qB  
allsel T&c[m!}X|t  
SOLVE

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