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

大家好，我做的是一个桩土接触分析问题，想把一个桩从土中拔出，施加两个荷载步，第一步，施加重力荷载，第二步把桩从中拔出，在运行第二步时出了问题，算不下去，请大师不吝指教 ，不胜感激。命令流如下： (x9d7$2  
FINI #ej^K |Qx  
/CLEAR s)%RmsdL  
GDC`\
cy  
/filname,pile-soil contact Xr~6_N{J  
/title，pile X ^>o/U  
!/config,nres,2000000 ,J?Hdy:R  
/units,si !国际单位制 -}k'a{sj=  
lpkg(J#&  
fe=0.1 !摩擦系数 a#W:SgE?Y  
fkn=10 !法向接触刚度  G~T]
m .  
fkt=3!切向接触刚度 9&5\L  
ftoln=1 !初始渗透因子 06 gE;iT  
fact=2!静动摩擦系数 2X2,(D!  
%7BVJJp2  
/prep7 e!yUA!x`u  
N m-{$U  
et,1,solid45 Nd/iMV6V;  
j.4oYxK!s/  
mp,ex,1,7.2e10 !桩的弹性模量，根据管壁厚度折算 +]Of f^s  
mp,nuxy,1,0.3 !桩的泊松比 V%&t'H{  
mp,dens,1,2000 !桩的密度 th8f  
!mp,mu,1,fe !桩体性质 |E-0P=h  
:qy`!QPUm  
mp,ex,2,2e8 !土的弹性模量，经验取值 w[PWJ! <  
mp,nuxy,2,0.3 !土的泊松比 y|b|_eE?{  
mp,dens,2,2000 !土的饱和密度 qOz,iR?}  
tb,dp,2 $DC*&hqpt  
tbdata,1,0,30 !粘聚力c为0,摩擦角为30 W;*rSK|(Sc  
L'4o
b4r{L  
!创建桩模型 PPT"?lt*&  
cylind,0,0.1,-0.1,7.1,0,90 &NV[)6!  
!创建土体模型 )NZ6!3[@  
cylind,0.1,0.6,0,7,0,90 
E!'H,#"P  
%>'2E!%  
!桩体划分网格 J) v~  
allsel,all $enh>!mU  
vsel,s,,,1 _#9:cH*  
aslv,s u4B,|_MK  
lsla,s f&t]O$  
lesize,1,,,6 *!UY;InanX  
lesize,2,,,6 ,-A8;DW]^J  
lesize,3,,,6 5=Mm=HyI2  
lesize,4,,,6 phSF.WC  
lesize,5,,,6 xI'<4lo7Z  
lesize,6,,,6 !mK
[kXo  
lesize,7,0.05 \/4ipU.  
lesize,8,0.05 +NPk9jn  
lesize,9,0.05 &|P@$O>  
type,1 dC@aQi6{6  
mat,1 vOsd>3"  
mshape,0,3d 9Qp39(l:  
mshkey,1 cs`/^2Vf"#  
vmesh,all O z%K*  
Y."ujo#bB  
vsel,s,,,2 ！土层，划分网格 c+AZ(6O?\  
vplot %a+X\\v2  
aslv,s 1(M0C[P  
lsla,s UC^&& 2maI  
lplot )'\Jp 7*3  
lesize,10,,,6 [.B)W);  
lesize,11,,,6 L7mN&Xr  
lesize,12,,,6 _lb ^  
lesize,13,,,6 1Nx.
aji  
lesize,14,,,6 ME~ga,|K  
lesize,15,,,6 vTjgW?9  
lesize,16,,,6 boo,KhW'Y  
lesize,17,,,6 R|H9AM ~E  
lesize,18,0.05
eA&hiAP/  
lesize,19,0.05 ~c[}%Ir>  
lesize,20,0.05 a&)0_i:r  
lesize,21,0.05 _Jj/"?  
type,1 Pgg6(O9}B^  
mat,2 qie7iE`o  
mshape,0,3d i^sK+v  
mshkey,1 YE&"IH]lF  
vmesh,all zvL&V .>  
La?q>  
!创建接触对 ~\/>b}^uf'  
!用接触向导建立接触对 $*K
5  
allsel 0CI?[R\  
et,2,170 vP&dvAUF  
et,3,173 @lAOi1m,,  
r,1 Z$0r+phQk=  
real,1 b
].:2  
keyopt,3,9,0 !消除初始渗透 @~fg[)7M  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 ws tI8">  
MK[l*=\s  
桩侧，建立目标面 I#@iA!  
r,1,,,fkn,ftoln, , , hC<X\yxe  
rmore,,,,,,fkt Vn? %w~0!  
RMORE,,,, ,, t$r^'
ZN  
RMORE,,,fact,,, ?UQVmE&  
mp,mu,3,fe RH ow%2D  
bv0B  
!创建接触对 ,5'LbO-  
vsel,s,,,1 :SG9ygq'  
asel,s,,,3 >b,o yM  
aplot XEV-D9n  
cm,target,area dN;kYWRK  
type,2 k(o[T),_%0  
real,1 xBU\$ToC  
mat,3 )gV+BHK  
nsla,s,1 ;OmmXygl  
esurf,all \(.&E`r  
esel,s,type,,2 Jl&bWp^3  
eplot Y5=~>*e  
/psymb,esys,1 j11\t  
allsel!很重要 !U}A1)  
,Ihuo5>/z  
!定义接触面 @B ~![l  
vsel,s,,,2 :gY$/1SYD  
vplot +GI[ Kq  
asel,s,,,9 C<fWDLwYqV  
aplot WKX5Dl  
cm,contact,area ;_K+b,  
type,3 cO<]%L0  
real,1 KgVit+4u/  
mat,3 ];6c/#2x  
nsla,s,1 "e g`3v  
esurf,all rwFR5  
esel,s,type,,3 %@$h?HP  
eplot [y}/QPR  
/psymb,esys,1 8+ eZU<\B(  
allsel ^G=wRtS  
save i9k7rEW^  
@8Drhx  
!侧面施加对称约束 y#HD1SZ  
allsel (p
`'Okw  
CSYS,0 !
^!<Xz;  
asel,s,loc,x,0 C=@BkneQ  
da,all,symm MLY19;e  
asel,s,loc,y,0 zy4AFW  
da,all,symm >1a-}>r  
E!Hq%L!/  
allsel,all Vj4 if@Z  
!asel,s,,,2!桩顶面2未约束 xq=+M!V  
!da,all,all >joGGT  
asel,s,,,6 F/ 2@%,2n  
da,all,ALL O;f^'N  
!约束外侧面 [e3|yE6  
asel,S,,,8 4C[,S|J  
NSLA,S,1 -'JTVfm.  
CSYS,1 fOJk+? c  
NROTAT,ALL ;|w &n  
D,ALL,UX Rp A76ug  
ALLSEL z=!$3E ecr  
SAVE # j*$ `W;  
finish!接触对定义结束 C!XI0d  
!求解过程 !$AVlMnJ  
/solu rfYu8-  
J"|)?$d]z  
!约束结束! vX)6N#D!  
!定义荷载步 <qZXpQ#  
!第一步荷载，重力荷载。 t*<vc]D  
acel,0,0,9.8 w>; :mf  
nlgeom,on ! 打开大变形选项 COFs?L.`  
autot,on ! 打开自动时间步 +@]1!|@(  
time,1 ]l+Bg;F#V  
nsub,10 /V~L:0%  
outres,all,all \l{*1lQ`  
allsel P~_CDh.N  
SOLVE )(*A1C[  
!第二部荷载,拉动桩 *x`z5_yfO  
nsubst,100,500,5 i2.y)K)  
outres,all,all FFbMG:>:  
time,101 2iI"|k9M  
nlgeom,on ! 打开大变形选项 8_>\A= E  
autot,on ! 打开自动时间步 51.F,uY  
nsel,s,loc,z,7.1 :84ja>`c  
D,ALL,UZ,0.1 a\vf{2  
l}335;(  
lnsrch,on !自动线性搜索 riZFcV
sB  
pred,on !打开预测求解器 W)^:*z  
G6JyAC
9j  
allsel
'15j$q  
SOLVE

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