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

大家好，我做的是一个桩土接触分析问题，想把一个桩从土中拔出，施加两个荷载步，第一步，施加重力荷载，第二步把桩从中拔出，在运行第二步时出了问题，算不下去，请大师不吝指教 ，不胜感激。命令流如下： [xh*"wT#g  
FINI MG /,==  
/CLEAR sz--27es  
Z_gC&7+  
/filname,pile-soil contact ?:|-Dq,  
/title，pile BYY RoE[P  
!/config,nres,2000000 ?<Y+peu  
/units,si !国际单位制 *jMk/9oa<N  
c}[+h5  
fe=0.1 !摩擦系数 ;Nr]X  
fkn=10 !法向接触刚度 w'X]M#Q><  
fkt=3!切向接触刚度 @\w,otT  
ftoln=1 !初始渗透因子 \v6M:KR5/  
fact=2!静动摩擦系数 x f4{r+  
kAM1TWbaVQ  
/prep7 M0|'f'  
c9j*n;Q  
et,1,solid45 ](eN@Xi&@  
KOg?FmD  
mp,ex,1,7.2e10 !桩的弹性模量，根据管壁厚度折算 !;6Jng%  
mp,nuxy,1,0.3 !桩的泊松比 rdnRBFt  
mp,dens,1,2000 !桩的密度 r-kMLw/)  
!mp,mu,1,fe !桩体性质 577:u<Yt  
X%bFN  
mp,ex,2,2e8 !土的弹性模量，经验取值 Qp2I[Ioz3  
mp,nuxy,2,0.3 !土的泊松比 ES<{4<Kpx  
mp,dens,2,2000 !土的饱和密度 5efxEt>U  
tb,dp,2 N~pIC2Woo  
tbdata,1,0,30 !粘聚力c为0,摩擦角为30 D;Gq)]O  
~|wos-nM  
!创建桩模型 =e><z9hY  
cylind,0,0.1,-0.1,7.1,0,90 Jdy<w&S  
!创建土体模型 :Vg}V"QR  
cylind,0.1,0.6,0,7,0,90 55t\Bms{  
|3K)$.6~  
!桩体划分网格 .+y>8h3{  
allsel,all xt4)Ya  
vsel,s,,,1 -ng=l;  
aslv,s n:U>Fj>q  
lsla,s 7ou46v|m5  
lesize,1,,,6 hNXBVIL<&  
lesize,2,,,6 h{lDxOH*  
lesize,3,,,6 <Fs-3(V+\  
lesize,4,,,6 JK#vkCkyM  
lesize,5,,,6  ni  
lesize,6,,,6 6zIgQ4Bp24  
lesize,7,0.05 1\dn1Hh  
lesize,8,0.05 <B /5J:o<  
lesize,9,0.05 XE}H3/2  
type,1 }a&mY^  
mat,1 Q36)7=at  
mshape,0,3d >EXb|vw
 
mshkey,1 E)ugLluL  
vmesh,all 'mdMq=VI  
(1pEEq84  
vsel,s,,,2 ！土层，划分网格 +VEU:1Gt  
vplot [ neXFp}S  
aslv,s Uz%Z&K  
lsla,s `-B+JQmen  
lplot :f'&z47  
lesize,10,,,6 {OT:3SS7  
lesize,11,,,6 SW(q$i  
lesize,12,,,6 ,]Zp+>{  
lesize,13,,,6 g"Qh]:  
lesize,14,,,6 e=/&(
Y  
lesize,15,,,6 <Y'>F!?#  
lesize,16,,,6 d0
er^ ~  
lesize,17,,,6 ^T4Ay=~{  
lesize,18,0.05 1'6cGpZY  
lesize,19,0.05 0MIUI<;j  
lesize,20,0.05 %-;bu|  
lesize,21,0.05 i(2y:U3[@  
type,1 FM^9
}*  
mat,2 < /p8r  
mshape,0,3d i}TwOy<4s  
mshkey,1 }bB`(B,m  
vmesh,all f?Bj _z  
qH%")7>  
!创建接触对 `GQ{*_-  
!用接触向导建立接触对  eIPG#A  
allsel @?B6aD|jE  
et,2,170 kpfwqHT  
et,3,173 3$M3Q]z  
r,1 v.Ogf5  
real,1 ]Ql 0v"` F  
keyopt,3,9,0 !消除初始渗透 lx,`hl%  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 O>DNC-m)i{  
#P2;K dDO  
桩侧，建立目标面 Mxz,wfaH>  
r,1,,,fkn,ftoln, , , c"gsB!xh  
rmore,,,,,,fkt sLcY,AH  
RMORE,,,, ,, j kn^Z":  
RMORE,,,fact,,, I#A2)V0P)  
mp,mu,3,fe h.tY 'F  
0x)dnq\  
!创建接触对 M<8ML!N0;t  
vsel,s,,,1 9\O(n>  
asel,s,,,3 ]j.k?P$U}  
aplot oJ;rc{n-  
cm,target,area aSaAC7sFk  
type,2 |ek ak{js  
real,1 32XS`Z  
mat,3  ;\b@)E}  
nsla,s,1 )ctr"&-  
esurf,all @Gjny BJ  
esel,s,type,,2 ?{J!#`tfV  
eplot @YL}km&Fw  
/psymb,esys,1 R"6;NPeo  
allsel!很重要 EZr6oO@Nc  
mVH,HqsXa  
!定义接触面 -$4kBYC l+  
vsel,s,,,2 @Jv# fr  
vplot cqL(^R.  
asel,s,,,9 $Tl<V/  
aplot V7<eQ0;m  
cm,contact,area :\}U9QfCw  
type,3 -@bp4Z=  
real,1 f{Qp  
mat,3 !sIwFv)  
nsla,s,1 ]A:(L9  
esurf,all H7IW"UkBR  
esel,s,type,,3 9Vh
_[^bR  
eplot ubGs/Vzye  
/psymb,esys,1 !L\'Mk/=A  
allsel xks?y.wA  
save Lar r}o=  
+l\<?  
!侧面施加对称约束 F1u)i  
allsel d/j@_3'  
CSYS,0 .LE+/n  
asel,s,loc,x,0 PW[6/7  
da,all,symm |>utWT]S  
asel,s,loc,y,0 H7e /  
da,all,symm jX
cNAl  
h
[U7!aM  
allsel,all eqFvrESN~=  
!asel,s,,,2!桩顶面2未约束 nxRwWj57  
!da,all,all _X%Dw  
asel,s,,,6 YOw?'+8  
da,all,ALL t.=Oj  
!约束外侧面 T\3[F%?  
asel,S,,,8 ~mA7pOHj  
NSLA,S,1 sv=^k(d3  
CSYS,1 x;U|3{Io  
NROTAT,ALL V*qY"[  
D,ALL,UX .Ty,_3+{#p  
ALLSEL ks:{TA27  
SAVE ?`. XK}  
finish!接触对定义结束 f)Y  
!求解过程 BI/y<6#rR  
/solu SfDQ;1?  
R+Lk~X^*l'  
!约束结束! qY}Cg0[@g  
!定义荷载步 |
cu`f{E2]  
!第一步荷载，重力荷载。 AN10U;p/O  
acel,0,0,9.8 #: hVF/  
nlgeom,on ! 打开大变形选项 22vq=RO7Z  
autot,on ! 打开自动时间步 -3k;u  
time,1 TcZN%  
nsub,10 zPBfiK_hV  
outres,all,all K -E`y  
allsel SsRVd^=;x  
SOLVE uGCtLA+sL  
!第二部荷载,拉动桩 E\!n4
9  
nsubst,100,500,5 ;IhPvff  
outres,all,all 3ZN>9`  
time,101 3Uo]>BG  
nlgeom,on ! 打开大变形选项 d$Em\*C  
autot,on ! 打开自动时间步 %BICt @E  
nsel,s,loc,z,7.1 |Mlh;  
D,ALL,UZ,0.1 QK7e|M  
Bm<^rhJ9  
lnsrch,on !自动线性搜索 Y
-~;E3(  
pred,on !打开预测求解器 s!`H  
dJE`9$j
N  
allsel A^vvST%7  
SOLVE

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