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

大家好，我做的是一个桩土接触分析问题，想把一个桩从土中拔出，施加两个荷载步，第一步，施加重力荷载，第二步把桩从中拔出，在运行第二步时出了问题，算不下去，请大师不吝指教 ，不胜感激。命令流如下： >H]|R }h  
FINI C'joJEo  
/CLEAR 9_%??@^>  
t8~isuiK  
/filname,pile-soil contact p$XvVzW#<  
/title，pile Tnnj8I1v  
!/config,nres,2000000 >#5jO9  
/units,si !国际单位制 R, UYwI  
!<0 `c  
fe=0.1 !摩擦系数 '=ydU+X  
fkn=10 !法向接触刚度 >ZnnGX6$(  
fkt=3!切向接触刚度 R'B-$:u  
ftoln=1 !初始渗透因子 *HQ>tvUh  
fact=2!静动摩擦系数 iz6+jHu'l  
:LTjV"f  
/prep7 AK$i0Rn;pm  
>ti)m >f  
et,1,solid45 4 :M}Vz-  
!=vd:,  
mp,ex,1,7.2e10 !桩的弹性模量，根据管壁厚度折算 xL!@$;J  
mp,nuxy,1,0.3 !桩的泊松比 F%:o6mT  
mp,dens,1,2000 !桩的密度  4{D^ 4G  
!mp,mu,1,fe !桩体性质 Y]"lcr}  
-^$IjK-N  
mp,ex,2,2e8 !土的弹性模量，经验取值 %u*H
No  
mp,nuxy,2,0.3 !土的泊松比 f2i9UZ$=e!  
mp,dens,2,2000 !土的饱和密度 207h$a,  
tb,dp,2 wF{M"$am  
tbdata,1,0,30 !粘聚力c为0,摩擦角为30 U%aDkC+M  
-dH]_  
!创建桩模型 3\a VZx!  
cylind,0,0.1,-0.1,7.1,0,90 TA!6|)BUW  
!创建土体模型 F<Z13]|  
cylind,0.1,0.6,0,7,0,90 ZtOv'nTD  
cIXqnb  
!桩体划分网格 D4U<Rn6N_5  
allsel,all f(*iagEy  
vsel,s,,,1 1<pb=H  
aslv,s X->` ~-aj  
lsla,s C=P}@|K  
lesize,1,,,6 0 TOw4pC  
lesize,2,,,6 NxN~"bfh  
lesize,3,,,6 I\$?'q>  
lesize,4,,,6 Ihx[S!:  
lesize,5,,,6 }ykc AK3U  
lesize,6,,,6 fI-f Gx  
lesize,7,0.05 ePiZHqIsv/  
lesize,8,0.05 W$Z8AZ{E  
lesize,9,0.05 G,,f' >  
type,1 p1 mY!&e(  
mat,1 #c`/ f6z  
mshape,0,3d
y6gao
j  
mshkey,1 UkM#uKr:  
vmesh,all *~8g:;u  
MZdj!(hO  
vsel,s,,,2 ！土层，划分网格 Xrzpn&Y=#  
vplot  lLNI5C  
aslv,s +pbP;zu  
lsla,s Q`z2SYz>  
lplot Z>ztFU  
lesize,10,,,6 FdS'0#$  
lesize,11,,,6 [lE^0_+  
lesize,12,,,6 QTjnXg?Ri  
lesize,13,,,6 (4FZK7Fm  
lesize,14,,,6 %H=^U
8WB  
lesize,15,,,6 F2dwT  
lesize,16,,,6 ;H"OZRQ  
lesize,17,,,6 Bh()?{q  
lesize,18,0.05 *I`Eb7 ^  
lesize,19,0.05 X9ec*x  
lesize,20,0.05 a#+;BH1  
lesize,21,0.05 P&AaD!Qn  
type,1 -;GB Xq  
mat,2 1A N)%  
mshape,0,3d ``Rg0o  
mshkey,1 ^'4I%L"  
vmesh,all X92I==-w  
o>
yo9n%t  
!创建接触对 P7;q^jlB  
!用接触向导建立接触对 s~g]`/h$r  
allsel `k3s
l 0z%  
et,2,170 oX?~  
et,3,173 J&;' gT  
r,1 +@anYtv%7  
real,1 4&H&zST//m  
keyopt,3,9,0 !消除初始渗透 ^*WO*f>y  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 Q57Z~EsF  
$A>]lLo0  
桩侧，建立目标面 gflO0$i  
r,1,,,fkn,ftoln, , , ky[^uQ>0  
rmore,,,,,,fkt /@6T~XY M  
RMORE,,,, ,, ^aY,Wq  
RMORE,,,fact,,, ^fti<Lw5  
mp,mu,3,fe  T_<:  
%wIb@km  
!创建接触对 >S:+&VN`M  
vsel,s,,,1 aY\(R02B  
asel,s,,,3 &M2fcw?  
aplot }Hb_8P  
cm,target,area ,ydn]0SS  
type,2 ?1Uq ud  
real,1 C.+:FY.H  
mat,3 HLl"=m1/>  
nsla,s,1 =${ImMwj  
esurf,all r1.OLn?C  
esel,s,type,,2 8%^W<.Y  
eplot ^fvx2<  
/psymb,esys,1 Ej'N!d.  
allsel!很重要 JiRfLB  
~uEI}z
 
!定义接触面 -7A2@
g  
vsel,s,,,2 &5R-bYGW  
vplot xWKUti i  
asel,s,,,9 %?!TqJT?{  
aplot TU?$yNE  
cm,contact,area -FQS5Zb.!  
type,3 JO\KTWtjO  
real,1 9f$3{ g{m  
mat,3 ,w7ZsI4:[  
nsla,s,1 |}<!O@<|  
esurf,all WRN}>]NgQ  
esel,s,type,,3 z=DK(b;$z  
eplot 2yndna-  
/psymb,esys,1 ux>LciNq  
allsel r+}<]?aT>-  
save XI;F=r}'  
?,
FL"ye  
!侧面施加对称约束 G3de<?K.[V  
allsel dlc'=M  
CSYS,0 !9EbG  
asel,s,loc,x,0 'o-J)+oa  
da,all,symm j4Y] 8  
asel,s,loc,y,0 BmCBC,j<v>  
da,all,symm XC\'8hL:  
=<z~OE'lV  
allsel,all KZwzQ"Hl  
!asel,s,,,2!桩顶面2未约束 qF4DX$$<  
!da,all,all #L xfE<^  
asel,s,,,6 n16TQe"8  
da,all,ALL +!<{80w  
!约束外侧面 <`*v/D7\02  
asel,S,,,8 n|!O .+\b  
NSLA,S,1 t++\&!F  
CSYS,1 luV_  
NROTAT,ALL ^j]"5@f  
D,ALL,UX zl["}I(*n  
ALLSEL .;NoKO7)  
SAVE zjwo"6c>  
finish!接触对定义结束 q@iZo,Yk  
!求解过程 T2i\S9X  
/solu $ 7uxReFZR  
tq5o  
!约束结束! a|@1RH>7H  
!定义荷载步 jp+#N pH  
!第一步荷载，重力荷载。 kl9<l*  
acel,0,0,9.8 pt8#cU
\  
nlgeom,on ! 打开大变形选项 jX 6+~  
autot,on ! 打开自动时间步 f{f|frs  
time,1 }7/Ob)O  
nsub,10 Y# ?M%I%j  
outres,all,all /lN09j  
allsel \DRYqLT`  
SOLVE 4L[-[{2  
!第二部荷载,拉动桩 w0$l3^}z  
nsubst,100,500,5 C~:!WRCz  
outres,all,all wfc+E9E  
time,101 Op:7EdT#  
nlgeom,on ! 打开大变形选项 g;<_GL  
autot,on ! 打开自动时间步 !FJ_\UST0  
nsel,s,loc,z,7.1 h_ J|uu  
D,ALL,UZ,0.1 y{1|@?ii  
b. oA}XP  
lnsrch,on !自动线性搜索 [bJAh ` I  
pred,on !打开预测求解器 8??%H7~  
d
#wK  
allsel SZ"^>}zl=  
SOLVE

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