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

大家好，我做的是一个桩土接触分析问题，想把一个桩从土中拔出，施加两个荷载步，第一步，施加重力荷载，第二步把桩从中拔出，在运行第二步时出了问题，算不下去，请大师不吝指教 ，不胜感激。命令流如下： WOeLn[  
FINI vv5rA 6+  
/CLEAR ;l < amB  
IEkbVIA(  
/filname,pile-soil contact [;:ocy  
/title，pile E4.A$/s8[  
!/config,nres,2000000 YZ<5-C  
/units,si !国际单位制 ^![7X'!;pt  
#P!M"_z  
fe=0.1 !摩擦系数 /,$V/q+  
fkn=10 !法向接触刚度 )}_}D+2  
fkt=3!切向接触刚度 muFWFq&yP  
ftoln=1 !初始渗透因子 4TiH
h  
fact=2!静动摩擦系数 t9lf=+%s  
V<PH5'^$j  
/prep7 iVmy|ewd  
;.R) uCd{=  
et,1,solid45 0|=y#`;,Z  
'b
g'^PN>z  
mp,ex,1,7.2e10 !桩的弹性模量，根据管壁厚度折算 s8';4z  
mp,nuxy,1,0.3 !桩的泊松比 K;w]sN+I  
mp,dens,1,2000 !桩的密度 `v2Xp3o4f  
!mp,mu,1,fe !桩体性质 Jv]$@>#  
ztp2j%'  
mp,ex,2,2e8 !土的弹性模量，经验取值 QWC C  
mp,nuxy,2,0.3 !土的泊松比 ;!3: 3;  
mp,dens,2,2000 !土的饱和密度 l=$?#^^ /  
tb,dp,2 [ThAvQ_$  
tbdata,1,0,30 !粘聚力c为0,摩擦角为30 'coY`B; 8  
iYlkc  
!创建桩模型 gieso
f  
cylind,0,0.1,-0.1,7.1,0,90 l$R9c+L=  
!创建土体模型 lWP]}Uy=5~  
cylind,0.1,0.6,0,7,0,90 g{2~G
6%;0  
u1.0-Y?  
!桩体划分网格 %8YUK/(|n  
allsel,all ?papk4w  
vsel,s,,,1 Q`-Xx  
aslv,s S&J5QZjC  
lsla,s ~JS@$#  
lesize,1,,,6 S{{wcH$n'i  
lesize,2,,,6 xjrL@LO#  
lesize,3,,,6 |PGTP#O
<  
lesize,4,,,6 3`NSSS  
lesize,5,,,6 n+2>jY  
lesize,6,,,6 9}a&:QTHR  
lesize,7,0.05 8-@HzS%  
lesize,8,0.05 ]U}B~Y  
lesize,9,0.05 g\iSc~%?  
type,1 |nqN95'u+]  
mat,1 =gZA9@]W2  
mshape,0,3d ,TYFPulYcp  
mshkey,1 r!|h3*YA  
vmesh,all nmyDGuzk  
38:5g_  
vsel,s,,,2 ！土层，划分网格 vRDs~'f  
vplot -uhVw_qq#  
aslv,s \H1t<B,  
lsla,s $)6M@S  
lplot \pP1k.~UnC  
lesize,10,,,6 nu2m5RYx  
lesize,11,,,6 Q6N?cQtOT  
lesize,12,,,6 *,=8x\Shp  
lesize,13,,,6 x]4>f[>*>  
lesize,14,,,6 /R>YDout}  
lesize,15,,,6 ]a[2QQ+g  
lesize,16,,,6 Ua~8DdW  
lesize,17,,,6 9eR";Wm])  
lesize,18,0.05 v0T?c53?  
lesize,19,0.05 DwWm(8&6;}  
lesize,20,0.05 /neY2D6  
lesize,21,0.05 =CjWPZShV  
type,1 h*3{IHAQ  
mat,2 ^E~1%Md.  
mshape,0,3d Deq@T {  
mshkey,1 Kp;a(D  
vmesh,all 9XUk.Nek  
`f<w+u  
!创建接触对 V`1x!
[\  
!用接触向导建立接触对 $P2*q
pqy  
allsel NH*"AE;  
et,2,170 2~vvE  
et,3,173 vjA!+_I6  
r,1 BQs\!~Ux2  
real,1 /fU-0a8  
keyopt,3,9,0 !消除初始渗透 [lu+"V,<LJ  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 uWP0(6 %  
cxx8I  
桩侧，建立目标面 @3~Wukc  
r,1,,,fkn,ftoln, , , Y3-f68*(  
rmore,,,,,,fkt "_=t1UE  
RMORE,,,, ,, 1{1mL-I;  
RMORE,,,fact,,, KUUZN  
mp,mu,3,fe ^ q3H  
B$7Cjv  
!创建接触对 ~aXJ5sY"f&  
vsel,s,,,1 05 .EI)7  
asel,s,,,3 j9-.bGtm?.  
aplot ]3Jb$Q@  
cm,target,area Tr;&bX5]H  
type,2 ,|j\x  
real,1 E!>MJlA:k6  
mat,3 Jid_&\  
nsla,s,1 +"|TPKas  
esurf,all 2Zg%4/u,Zp  
esel,s,type,,2 _x$\E  
eplot ]LMtZUz  
/psymb,esys,1 HYZp=*eb  
allsel!很重要 ofsLx6Po  
o|_9%o52'  
!定义接触面 (nYGN$qC9  
vsel,s,,,2 Xo6zeLHO  
vplot by$mD_sr  
asel,s,,,9 y I[kaH"J  
aplot n]>L"D,  
cm,contact,area 6_:KFqc W  
type,3 mcs!A/]<  
real,1 Digx#'#jf  
mat,3 Ri*mu*r\}  
nsla,s,1 p4$4;)  
esurf,all _:Jma  
esel,s,type,,3 Z2-"NB  
eplot j%':M  
/psymb,esys,1 }C$D-fH8sW  
allsel NFf`V  
save Q@"!uB.e  
xhB-gG=  
!侧面施加对称约束 Pjq()\/[Z  
allsel \KKE&3=  
CSYS,0 q2;CvoF  
asel,s,loc,x,0 A-h[vP!v|  
da,all,symm o)'06FF\$  
asel,s,loc,y,0 )q[Wzx_ j<  
da,all,symm <L2z|%`  
]H<}6}
Gd  
allsel,all *q@3yB}  
!asel,s,,,2!桩顶面2未约束 >x ]{cb/m  
!da,all,all cqcH1aSv  
asel,s,,,6 O>9-iqP>`d  
da,all,ALL :w7?]y6~S  
!约束外侧面 i1GQ=@  
asel,S,,,8 <)"2rxX&5  
NSLA,S,1 MVEh<_  
CSYS,1 P!{J28dj  
NROTAT,ALL Bp4#"y2  
D,ALL,UX }{[JS=A^  
ALLSEL |lyspD  
SAVE
" ^!=e72  
finish!接触对定义结束 %UG|R:  
!求解过程 =>O{hT^F  
/solu Dw,LB>Eq,  
:z\||f  
!约束结束! hb>uHUb&  
!定义荷载步 *g4Uo{  
!第一步荷载，重力荷载。 {J/+KK  
acel,0,0,9.8 ~RV9'v4  
nlgeom,on ! 打开大变形选项 yvnvIy  
autot,on ! 打开自动时间步 2<d'!cm  
time,1 q3|SZoN  
nsub,10 \)2'+R  
outres,all,all  l .m #  
allsel kc2 8Q2  
SOLVE 7>EMr}f C  
!第二部荷载,拉动桩 Wxi;Tq9C@_  
nsubst,100,500,5 8[a N5M]  
outres,all,all jPg[LZQ'  
time,101 _j~y;R)  
nlgeom,on ! 打开大变形选项 c57`mOe/b  
autot,on ! 打开自动时间步 B|=S-5pv*  
nsel,s,loc,z,7.1 ]#rNz"  
D,ALL,UZ,0.1 7(H?3)%0  
V h5\'Sn  
lnsrch,on !自动线性搜索 4%6@MQ[  
pred,on !打开预测求解器 _6]tbni?v  
\JP9lJ3<  
allsel 8rNf4]5@X(  
SOLVE

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