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

大家好，我做的是一个桩土接触分析问题，想把一个桩从土中拔出，施加两个荷载步，第一步，施加重力荷载，第二步把桩从中拔出，在运行第二步时出了问题，算不下去，请大师不吝指教 ，不胜感激。命令流如下： Ij#mmj NW  
FINI 2|,$#V=  
/CLEAR >fC&bab  
']nI
a7  
/filname,pile-soil contact ;dzy5o3  
/title，pile t
utk*|S  
!/config,nres,2000000 \tgY2:  
/units,si !国际单位制 +p\E%<uQ
 
Wg}KQ6 6  
fe=0.1 !摩擦系数 9~UR(Ts}l  
fkn=10 !法向接触刚度 j e\!0{  
fkt=3!切向接触刚度 l+Wux$6U  
ftoln=1 !初始渗透因子 $wM..ee  
fact=2!静动摩擦系数 [(n5-#1S  
JO|j?%6YY  
/prep7 k[x-O?$O@  
Mk*4J]PP  
et,1,solid45 %j&vV>2  
p\[!=ZXFr\  
mp,ex,1,7.2e10 !桩的弹性模量，根据管壁厚度折算 FF8jW1  
mp,nuxy,1,0.3 !桩的泊松比 !RKuEg4hQ  
mp,dens,1,2000 !桩的密度 u#ya 8  
!mp,mu,1,fe !桩体性质 gT8(LDJ  
MD[hqshoh  
mp,ex,2,2e8 !土的弹性模量，经验取值 Mq91HmC(@  
mp,nuxy,2,0.3 !土的泊松比 &E`Nu (e  
mp,dens,2,2000 !土的饱和密度 B!&5*f}*  
tb,dp,2 !td!">r46e  
tbdata,1,0,30 !粘聚力c为0,摩擦角为30 n{QyqI  
MXSD8]
je  
!创建桩模型 g(&cq  
cylind,0,0.1,-0.1,7.1,0,90 NO*,}aeG  
!创建土体模型 u$JAjA
 
cylind,0.1,0.6,0,7,0,90 "(mJupI  
;2kQ)Bq"  
!桩体划分网格 2VV>?s  
allsel,all 12^uu)6Xm,  
vsel,s,,,1 0/QDfA?  
aslv,s oRbWqN`F.  
lsla,s NCk-[I?R  
lesize,1,,,6 ,3?=W/Um4  
lesize,2,,,6 8
O^x~[sQ  
lesize,3,,,6 v6)QLp
 
lesize,4,,,6 xsZN@hT  
lesize,5,,,6 wiI@DJ>E  
lesize,6,,,6 I!/EQO|  
lesize,7,0.05 `he# !"  
lesize,8,0.05 j}Tv/O,f  
lesize,9,0.05 [],[LkS  
type,1 'ON/WKJr|W  
mat,1 va@;V+cD  
mshape,0,3d ~|KqG  
mshkey,1 `v?hL
~  
vmesh,all rai'x/Ut}+  
:3M,]W]  
vsel,s,,,2 ！土层，划分网格 ?h`,@~6u  
vplot >9w^C1"  
aslv,s />xEpR3_A  
lsla,s !hHe`  
lplot G4i%/_JU  
lesize,10,,,6 S[L@8z.Sj  
lesize,11,,,6 4<s;xSCL  
lesize,12,,,6 qBk[Afjgz  
lesize,13,,,6 l i<9nMZ<  
lesize,14,,,6 cs2-
jbRn  
lesize,15,,,6
$l,U)  
lesize,16,,,6 _L8&.=4]i  
lesize,17,,,6 oN}\bK  
lesize,18,0.05 ~T}D#}  
lesize,19,0.05 Z*Ffdh>*:&  
lesize,20,0.05 +/&rO,Ql  
lesize,21,0.05 Y5E0n(Z  
type,1 -(57C*#ap  
mat,2 %>K(IRpMW  
mshape,0,3d ^fKKsfIf  
mshkey,1 |e8A)xM]wC  
vmesh,all "0!#De  
UusAsezm:  
!创建接触对 Z(:\Vj"  
!用接触向导建立接触对 jpi,BVTI-X  
allsel JSg=9p$  
et,2,170 :/6:&7s  
et,3,173 bN?*p($/  
r,1 A]"6/Lr9P  
real,1 *effDNE!  
keyopt,3,9,0 !消除初始渗透 Gh_5$@ hF  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 9ZOQNN<ex  
_ (b4|hJ'  
桩侧，建立目标面 kYS#P(1  
r,1,,,fkn,ftoln, , , h6~xz0,u  
rmore,,,,,,fkt N_D=j6B  
RMORE,,,, ,, j&DlI_  
RMORE,,,fact,,, UVXruH  
mp,mu,3,fe ]>)}xfL &,  
u9;3Xn8  
!创建接触对
0uKm)t/  
vsel,s,,,1 LEKE+775  
asel,s,,,3 Rq15AR  
aplot |%4nU#GoB  
cm,target,area em'3 8L|(  
type,2 tDAX pi(  
real,1 `LFT"qnp  
mat,3 T>|+cg  
nsla,s,1 q|YnNk>1  
esurf,all {{A=^rr%C  
esel,s,type,,2 `
mkOjsj &  
eplot '!X`X=  
/psymb,esys,1 pz2E+o  
allsel!很重要 wB8548C}-  
{(-TWh7V  
!定义接触面 K{y`Sb~k  
vsel,s,,,2 "j%Gr:a  
vplot G]l/L\{  
asel,s,,,9 1 =?pL$+G  
aplot d>M0:  
cm,contact,area H"+|n2E^  
type,3 /_<_X 7  
real,1 "% \y$  
mat,3 v'L"sgW6I  
nsla,s,1 !h&h;m/c  
esurf,all jhG6,;1zMI  
esel,s,type,,3 2aivc,m{r  
eplot &}gH!5L m  
/psymb,esys,1 (N}\Wft%  
allsel 2P57C;N8|  
save L%7WHtU*#  
+L?;g pVE&  
!侧面施加对称约束 k;umLyz  
allsel g3n>}\xG>  
CSYS,0 6mZpyt  
asel,s,loc,x,0 x=JZ"|TE  
da,all,symm F[ ^ p~u{  
asel,s,loc,y,0 ,vi6<C\  
da,all,symm (4l M3clF  
bN*zx)f  
allsel,all n\JSt}A  
!asel,s,,,2!桩顶面2未约束 ;).QhHeg>  
!da,all,all `5t~ Vlp  
asel,s,,,6 1%.CtTi  
da,all,ALL .Xta;Py|J  
!约束外侧面 cCtd\/ \  
asel,S,,,8 5k_%%><: q  
NSLA,S,1 i+-Y"vRi  
CSYS,1 Gd&G*x  
NROTAT,ALL I~ SFY>s  
D,ALL,UX +DT tKj  
ALLSEL DKQQZ`PF  
SAVE c1%ki%J#  
finish!接触对定义结束 a;7gy419<p  
!求解过程 mX SLH'  
/solu bxz6 >>  
tG,xG&  
!约束结束! .@(MNq{"6  
!定义荷载步 hEFn>  
!第一步荷载，重力荷载。 D//Ts`}+n  
acel,0,0,9.8 !Je!;mEvI  
nlgeom,on ! 打开大变形选项 M>Ws}Y
 
autot,on ! 打开自动时间步 xs  >Y  
time,1 (B+zh  
nsub,10 9&c *%mm  
outres,all,all P>6wr\9i[  
allsel K0^+2lx  
SOLVE >*w(YB]/$V  
!第二部荷载,拉动桩 z81`Lhg6  
nsubst,100,500,5 Lp||C@h~  
outres,all,all b|Ed@C  
time,101 p t{/|P  
nlgeom,on ! 打开大变形选项 . =A|  
autot,on ! 打开自动时间步 .Wyx#9  
nsel,s,loc,z,7.1 l&Fx< W  
D,ALL,UZ,0.1 .9bP8u2B{  
]s_BOt  
lnsrch,on !自动线性搜索 $/i;UUd  
pred,on !打开预测求解器 /Bnh%6#ab  
& V/t0  
allsel 8-vNXvl  
SOLVE

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