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

大家好，我做的是一个桩土接触分析问题，想把一个桩从土中拔出，施加两个荷载步，第一步，施加重力荷载，第二步把桩从中拔出，在运行第二步时出了问题，算不下去，请大师不吝指教 ，不胜感激。命令流如下： [qYr~:`-[  
FINI 
`S {&gl  
/CLEAR 8`q7Yss6F  
hS]w A"\87  
/filname,pile-soil contact ?E([Nc0T  
/title，pile |YJ83nSO~  
!/config,nres,2000000 3A#Tn7  
/units,si !国际单位制 iM+`7L'  
Y %JQ  
fe=0.1 !摩擦系数 qA)OkR'm  
fkn=10 !法向接触刚度 Mbi]EZ  
fkt=3!切向接触刚度 Gt`7i(  
ftoln=1 !初始渗透因子 *M.xVUPr  
fact=2!静动摩擦系数 ( ayAP  
y?$DDD  
/prep7 1+x" 5<(W  
tb3VqFx  
et,1,solid45 4qQ,1&!]S  
X]0>0=^  
mp,ex,1,7.2e10 !桩的弹性模量，根据管壁厚度折算 @+u>rS|IB  
mp,nuxy,1,0.3 !桩的泊松比 mayJwBfU  
mp,dens,1,2000 !桩的密度 TQa}Ps  
!mp,mu,1,fe !桩体性质 I.<c{4K5  

@R[{  
mp,ex,2,2e8 !土的弹性模量，经验取值 .G0 N+)  
mp,nuxy,2,0.3 !土的泊松比 >Fel) a  
mp,dens,2,2000 !土的饱和密度 pCIzpEsRs  
tb,dp,2 ^jdU4  
tbdata,1,0,30 !粘聚力c为0,摩擦角为30 H];B?G';C  
mDB  
!创建桩模型 5(1:^:LGK  
cylind,0,0.1,-0.1,7.1,0,90 mB?x_6#d9  
!创建土体模型 o5 L^  
cylind,0.1,0.6,0,7,0,90 c`;\sW-_W  
rO1!h%&o"  
!桩体划分网格 }Mcb\+[  
allsel,all RZ".?  
vsel,s,,,1 T<AT
&4  
aslv,s ZO0 Ee1/  
lsla,s %,$n^{v  
lesize,1,,,6 zG!nqSDG  
lesize,2,,,6 U.%Kt,qB  
lesize,3,,,6 UX 1 )((  
lesize,4,,,6 #/)t]&n  
lesize,5,,,6 *XVwTW[a  
lesize,6,,,6 :=*deZ<  
lesize,7,0.05 ~wc:/UM|  
lesize,8,0.05 E~'mxx~i  
lesize,9,0.05 &p0e)o~Ux  
type,1 / HTY>b  
mat,1 ()
3\(d5e  
mshape,0,3d 7s9h:/Lu  
mshkey,1 B/? L$m  
vmesh,all pSZ2>^";  
[K!9x
M6  
vsel,s,,,2 ！土层，划分网格 okNo-\Dh!  
vplot Tb/TP3N  
aslv,s J=4>zQ
LW  
lsla,s M6Fo.eeK3  
lplot }aCa2%  
lesize,10,,,6 /7Q|D sa  
lesize,11,,,6 yV30x9i!2  
lesize,12,,,6 8)i""OD@I  
lesize,13,,,6 x&}]8S)  
lesize,14,,,6 AOf4y&B>q  
lesize,15,,,6 ~E<2gMKjO  
lesize,16,,,6 4w4B\Na>l  
lesize,17,,,6 ,Kit@`P%  
lesize,18,0.05 NW3c_]`=  
lesize,19,0.05  >@ t  
lesize,20,0.05 9>""x
t
 
lesize,21,0.05 :!ya&o  
type,1 e|~MJu+1  
mat,2 ]S|FK>U[  
mshape,0,3d s='+[*&&  
mshkey,1 nQ%HtXt;  
vmesh,all v9$!v^U"D  
IgT`on3Y  
!创建接触对 3r\8v`^>  
!用接触向导建立接触对 sE1cvAw9l  
allsel Z6jEj9?O  
et,2,170 Ic&h8vSU  
et,3,173 G=dzP}B'WA  
r,1 a QH6akH  
real,1 1v o)]ff  
keyopt,3,9,0 !消除初始渗透 '6
KvB  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 WWT1=#"  
"L~(%Nx3  
桩侧，建立目标面 IRdt:B|@  
r,1,,,fkn,ftoln, , , !_S>ER  
rmore,,,,,,fkt %|Ps|iV  
RMORE,,,, ,, jYsAL=oh,*  
RMORE,,,fact,,, 5pO|^Gj1  
mp,mu,3,fe :5d>^6eoB?  
S63Zk0(25  
!创建接触对 kc,"w\ ai  
vsel,s,,,1 ~j0rORy]  
asel,s,,,3 ngY%T5-  
aplot !>UlvT-  
cm,target,area +)]YvZ6%[,  
type,2 UZZJtQt  
real,1 n|`3d~9$&  
mat,3 xad`-vw  
nsla,s,1 V"FQVtTx7  
esurf,all gcLz}84  
esel,s,type,,2 "fFSZ@,r  
eplot |WS)KR !  
/psymb,esys,1
R1jl<=  
allsel!很重要 {]vD@)k  
}
d$-:l,w  
!定义接触面 AeJ

;g  
vsel,s,,,2 rIb+c=|F  
vplot ,d G.67  
asel,s,,,9 Zg;$vIhn  
aplot v=`yfCX-qX  
cm,contact,area >rG>Bz^Pu  
type,3 50UdY9E_v}  
real,1 %36x'Dn?  
mat,3 OvdT* g=8*  
nsla,s,1 P rt} 01$  
esurf,all N2r zHK  
esel,s,type,,3 .nV2n@SR  
eplot %>_[
b,  
/psymb,esys,1 r8k(L{W  
allsel |rmelQ-  
save )U^=`* 7  
F_-Lu]*  
!侧面施加对称约束 qTO6I5u  
allsel )7q$PcY  
CSYS,0 e3~MU6  
asel,s,loc,x,0 -@TY8#O#-  
da,all,symm SnCwoxK  
asel,s,loc,y,0 |diI(2w  
da,all,symm P<<$o-a"  
J0G@]H  
allsel,all p|UL<M9{a]  
!asel,s,,,2!桩顶面2未约束 jnH\}IB  
!da,all,all lFBdiIw  
asel,s,,,6 QV4|f[Ki%  
da,all,ALL .6?"<zdPU  
!约束外侧面 gLiJ&H  
asel,S,,,8 9>}&dQ8  
NSLA,S,1 %&ejO=r  
CSYS,1 Z
9{~t  
NROTAT,ALL 8-$t7bV5  
D,ALL,UX _0H oJ  
ALLSEL {I
#_0Q,i  
SAVE ,GbmL8P7Y  
finish!接触对定义结束  tm1=  
!求解过程 ^@fD{]I  
/solu %L|fTndKH  
4;`z6\u9-  
!约束结束! T:}Q3  
!定义荷载步 fv+]iK<{  
!第一步荷载，重力荷载。 ^BsT>VSH6  
acel,0,0,9.8 X1o",,N^M  
nlgeom,on ! 打开大变形选项 oOFTQB_6  
autot,on ! 打开自动时间步 AGhenDNV  
time,1 5Rqdo\vE  
nsub,10 =?
sG~  
outres,all,all 5#Er& 6s  
allsel 6jE.X  
SOLVE M!&Hn,22  
!第二部荷载,拉动桩 ;?h[WIy  
nsubst,100,500,5 ^z,3#gK  
outres,all,all W"{Ggk`  
time,101 8
NN+Z<  
nlgeom,on ! 打开大变形选项 U S^% $Z:  
autot,on ! 打开自动时间步 ]HKt7 %,  
nsel,s,loc,z,7.1 x6ghO-s  
D,ALL,UZ,0.1 !V|{(>+<  
LIg1U  
lnsrch,on !自动线性搜索 /V
iY:-8s  
pred,on !打开预测求解器 GT{4L]C  
/A_:`MAZ  
allsel r"H::A  
SOLVE

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