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

大家好，我做的是一个桩土接触分析问题，想把一个桩从土中拔出，施加两个荷载步，第一步，施加重力荷载，第二步把桩从中拔出，在运行第二步时出了问题，算不下去，请大师不吝指教 ，不胜感激。命令流如下： 7zOhyl?  
FINI Q-M rH   
/CLEAR & 8ccrw  
u![4=w  
/filname,pile-soil contact
37J\i ]  
/title，pile UENYJ*tnP  
!/config,nres,2000000 #"=%b e3  
/units,si !国际单位制 "1_{c *ck  

"ugX /r$_  
fe=0.1 !摩擦系数 un.G6|S  
fkn=10 !法向接触刚度 >|v=Ba6R0  
fkt=3!切向接触刚度 (r.$%[,.<  
ftoln=1 !初始渗透因子 b$;HI7)/K  
fact=2!静动摩擦系数 2qt=jz\s  
z\X60T  
/prep7 nrxjN(9V%+  
ac2}3$u  
et,1,solid45 l :f9Ih  
Kh>?!`lL  
mp,ex,1,7.2e10 !桩的弹性模量，根据管壁厚度折算 9m%7
dsv  
mp,nuxy,1,0.3 !桩的泊松比 %{N>c:2I$  
mp,dens,1,2000 !桩的密度 $
ijx#a&O  
!mp,mu,1,fe !桩体性质 Wpiv1GZ%c8  
/E=h{|  
mp,ex,2,2e8 !土的弹性模量，经验取值 nwZ[Ygl|  
mp,nuxy,2,0.3 !土的泊松比 P}ehNt*($  
mp,dens,2,2000 !土的饱和密度 OI)&vQ5k  
tb,dp,2 ;m-6.AV  
tbdata,1,0,30 !粘聚力c为0,摩擦角为30 e-nWD  
SIKk|I)  
!创建桩模型 I@8+k&nXS  
cylind,0,0.1,-0.1,7.1,0,90 bx+(.F  
!创建土体模型 gubb .EY  
cylind,0.1,0.6,0,7,0,90 '?&B5C  
8Uj68Jl?  
!桩体划分网格 ^G]KE8  
allsel,all c@0l-R{q  
vsel,s,,,1 v1`bDS?*Q  
aslv,s g_3rEvf"4  
lsla,s jEC'l]l  
lesize,1,,,6 BsoFQw4$9  
lesize,2,,,6 FE$M[^1_  
lesize,3,,,6 eyI-s9#t  
lesize,4,,,6 ?7r
mwy\  
lesize,5,,,6 |NphG|  
lesize,6,,,6 Vvxc8v:  
lesize,7,0.05 `D&#U'wB   
lesize,8,0.05 eukX#0/^  
lesize,9,0.05 nOA,x  
type,1 YHwVj?6W  
mat,1 TMnT#ypf<5  
mshape,0,3d &-EyM*:u!  
mshkey,1 F9>(W#aC  
vmesh,all ;j[q?^ b  
IqR[&T)lj  
vsel,s,,,2 ！土层，划分网格 RW|UQY#  
vplot DzX6U[=  
aslv,s e8WPV  
lsla,s -jcrXskb&N  
lplot hW*^1%1  
lesize,10,,,6 `Gy>tD.#V-  
lesize,11,,,6 rRF+\cP?.  
lesize,12,,,6 \?5[RR  
lesize,13,,,6 +oxqS&$L  
lesize,14,,,6 pn ~/!y  
lesize,15,,,6 8EW`*
+%=  
lesize,16,,,6 ~|$) 1  
lesize,17,,,6 vZ1D3ytfG  
lesize,18,0.05 q\{;_?a  
lesize,19,0.05 'tq\<y
 
lesize,20,0.05 H4K(SGx  
lesize,21,0.05 4nVO.Ud0$X  
type,1 XY$cx~  
mat,2 .{-iq(3  
mshape,0,3d |r<.R>  
mshkey,1 W2X+NacD  
vmesh,all vl#V-UW$4P  
RN cI]oJ  
!创建接触对 ul5::  
!用接触向导建立接触对 |zr)hC  
allsel .q'FSEkMJ  
et,2,170 \k$]GK-  
et,3,173 B8BY3~}]  
r,1 y7}~T!UyfF  
real,1 gvwCoCbb  
keyopt,3,9,0 !消除初始渗透 %oSfL;W7  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 Q xj|lr  
Hsux>+Q  
桩侧，建立目标面  t`&s  
r,1,,,fkn,ftoln, , , !Cv<>_N).  
rmore,,,,,,fkt |eVTxeq  
RMORE,,,, ,, Cbf,X[u  
RMORE,,,fact,,, =:0IHyB#0  
mp,mu,3,fe e2VL/>y`  
U8.0L  
!创建接触对 i >/@]2  
vsel,s,,,1 ,rvZW}=  
asel,s,,,3 ":ws~
Zep  
aplot }gi`?58J6  
cm,target,area 9Rt(G_'  
type,2 G1n>@Y'j''  
real,1 +$pO  
mat,3 l*V72!Mv  
nsla,s,1 CsQ}P)  
esurf,all Xf'  
esel,s,type,,2 @)M.u3{\  
eplot H]]>sE  
/psymb,esys,1 S+E3;' H  
allsel!很重要 .jG.90  
iaq:5||,  
!定义接触面 hF{mm(qyv  
vsel,s,,,2 7zOvoQ}  
vplot n$\6}\k  
asel,s,,,9 E rRMiT  
aplot ;tIIEc  
cm,contact,area
h^3Vd K,  
type,3 DQHGq_unP  
real,1 GQc%
OQc\  
mat,3 W0C$*oe!_i  
nsla,s,1 fe]T9EDA  
esurf,all ^dp[Z,[1z  
esel,s,type,,3 nh_xbo5L[  
eplot O'?lW~CD.>  
/psymb,esys,1 *SpO|*'  
allsel /(}l[jf  
save epm8N /  
rH9|JEz  
!侧面施加对称约束 ?)ZLxLV::  
allsel PeR<FSF ,i  
CSYS,0 ]`+"o[  
asel,s,loc,x,0 Ub-k<]yZ  
da,all,symm lVKF^-i  
asel,s,loc,y,0 boHm1hPKS  
da,all,symm TuT=  
1,sO =p)Yg  
allsel,all r;s3(@[,@  
!asel,s,,,2!桩顶面2未约束 dIe 6:s  
!da,all,all Z_TbM^N  
asel,s,,,6 P|QnZ){  
da,all,ALL k@X As  
!约束外侧面 A 5\"e^>  
asel,S,,,8 os0"haOI9h  
NSLA,S,1 s</ktPtu  
CSYS,1 fTnyCaB  
NROTAT,ALL S&g-  
D,ALL,UX ?_`P;}4#  
ALLSEL 75cr!+  
SAVE 4iB)oR  
finish!接触对定义结束 Vf* B1Zb  
!求解过程 % w/1Uo24  
/solu 9hi(P*%q  
b+fy&rk@-  
!约束结束! Sr-^f
aL  
!定义荷载步 ;jfXU_K  
!第一步荷载，重力荷载。 kO O~%|1CP  
acel,0,0,9.8 '[M2Q"X  
nlgeom,on ! 打开大变形选项 xo@/k  
autot,on ! 打开自动时间步 {`SGB;ho  
time,1 "gVH;<&]  
nsub,10 F_ lj>;}a5  
outres,all,all N=
}Z#  
allsel SOY#, Zu  
SOLVE 7UnO/K7oB.  
!第二部荷载,拉动桩 G7#~=W 2M  
nsubst,100,500,5 Vam8NnZ|r  
outres,all,all t7& GCZ  
time,101 }z%/6`7)|  
nlgeom,on ! 打开大变形选项 MhXm-<4  
autot,on ! 打开自动时间步 5+PBS)pJ]%  
nsel,s,loc,z,7.1 uaMmiR  
D,ALL,UZ,0.1 uTFEI.N  
3xR#,22:}  
lnsrch,on !自动线性搜索 9n2%7dLQ*  
pred,on !打开预测求解器 ,}$x'8v  
j;c^pLUP  
allsel YRX2^v ^[  
SOLVE

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