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

大家好，我做的是一个桩土接触分析问题，想把一个桩从土中拔出，施加两个荷载步，第一步，施加重力荷载，第二步把桩从中拔出，在运行第二步时出了问题，算不下去，请大师不吝指教 ，不胜感激。命令流如下： {~NiGHY  
FINI bM7y}P5`1  
/CLEAR k(1]!c4J0  
-w nlJi1f  
/filname,pile-soil contact v`3q0,,  
/title，pile %^){Z,}M}  
!/config,nres,2000000 ,E,oz{,i(  
/units,si !国际单位制 eh_{-  
-JOtvJIQI  
fe=0.1 !摩擦系数 I;iJa@HWQ  
fkn=10 !法向接触刚度 {n9]ej^  
fkt=3!切向接触刚度 SXX6EIJr|  
ftoln=1 !初始渗透因子 LWgYGXWT"  
fact=2!静动摩擦系数 !K a!f1  
>DX\^86x  
/prep7 2eErvfC[  
0'u2xe  
et,1,solid45 ?K,xx
H  
j8WMGSrrF  
mp,ex,1,7.2e10 !桩的弹性模量，根据管壁厚度折算 9sYN7x  
mp,nuxy,1,0.3 !桩的泊松比 CVfQ  
mp,dens,1,2000 !桩的密度 $1<V'b[E  
!mp,mu,1,fe !桩体性质 MT
{7I"  
oE:9}]N_  
mp,ex,2,2e8 !土的弹性模量，经验取值 bOR1V\Jr$q  
mp,nuxy,2,0.3 !土的泊松比 N&g9z{m7  
mp,dens,2,2000 !土的饱和密度 9x;CJhX
 
tb,dp,2 !Ra.DSL  
tbdata,1,0,30 !粘聚力c为0,摩擦角为30 BR\3ij  
L=Cm0q 3v  
!创建桩模型 UHr0J jQK  
cylind,0,0.1,-0.1,7.1,0,90 H]e%8w))0  
!创建土体模型 vg@kPuOiO  
cylind,0.1,0.6,0,7,0,90 RC(fhqV  
W*A-CkrO  
!桩体划分网格 M$&aNt;  
allsel,all +~~FfIzf#  
vsel,s,,,1 mqL&bmT  
aslv,s !ceT>i90h  
lsla,s r[; .1,(  
lesize,1,,,6 SF$'$6x}  
lesize,2,,,6 #wz1uw[pI!  
lesize,3,,,6 i'Vrx(y3  
lesize,4,,,6 \uxDMKy  
lesize,5,,,6 hAx#5@*5  
lesize,6,,,6 3^p<Wx  
lesize,7,0.05 +$Q.N{LV  
lesize,8,0.05 ,<iJ#$: Sx  
lesize,9,0.05 y\-f{I  
type,1 \-(.cj)?  
mat,1 EHpIbj;n  
mshape,0,3d qMy>:,)Z
 
mshkey,1 AITV+=sN  
vmesh,all W vh3Y,|3  
1=LI))nV  
vsel,s,,,2 ！土层，划分网格 .48Csc-  
vplot 'iy &%?  
aslv,s <>3}<i<[&  
lsla,s ^Mw>'*5^  
lplot E`vCYhf{  
lesize,10,,,6  d+FS  
lesize,11,,,6 ,_HSvs7-  
lesize,12,,,6 E/x2LYH  
lesize,13,,,6 #H9J/k_  
lesize,14,,,6 ;-SFK+)R"  
lesize,15,,,6 2c?-_OCy;  
lesize,16,,,6 jSMvZJX3n  
lesize,17,,,6 MExP'9  
lesize,18,0.05 *n9t~t6GHg  
lesize,19,0.05 !uaV6K  
lesize,20,0.05 6ww4ZH?j  
lesize,21,0.05 aLr\Uq,83  
type,1 &YMz3ugI  
mat,2 3GMRH;/w  
mshape,0,3d s$g"6;_\  
mshkey,1 ;O7CahdF  
vmesh,all o$dnp`E  
Nb.AsIR^  
!创建接触对 Y{c_5YYf  
!用接触向导建立接触对 "4W@p'  
allsel ?UD2

}D[M  
et,2,170 ^kg[n908Nw  
et,3,173 #H]cb#  
r,1 cYBv}ylw}R  
real,1 4 ZD~i e  
keyopt,3,9,0 !消除初始渗透 02g!mJW>}y  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 3SbtN3  
xw?Mc{w  
桩侧，建立目标面 __x2xtrH  
r,1,,,fkn,ftoln, , , C@!C='b,  
rmore,,,,,,fkt Z";~]]$!Y  
RMORE,,,, ,, K9JW&5Q  
RMORE,,,fact,,, X\X*-.]{  
mp,mu,3,fe c
\2+f7o@  
h\+U+?u  
!创建接触对 N}ur0 'J0  
vsel,s,,,1 V6)e Jy  
asel,s,,,3 bWc3a  
aplot Y\e,#y  
cm,target,area Af8&PhyrU  
type,2 K-D{Z7J^l  
real,1 W<Ms0  
mat,3 #vk-zx*v7=  
nsla,s,1 .xXe *dm%  
esurf,all }9n{E-bj*  
esel,s,type,,2 ex_Zw+n  
eplot IyTL|W6  
/psymb,esys,1 ;CbQ}k  
allsel!很重要 @^g/`{j>J  
5DgfrX  
!定义接触面 |&JCf=  
vsel,s,,,2 6=V&3|"  
vplot T /iKz  
asel,s,,,9 jJ^p ?  
aplot 3GEI)!  
cm,contact,area v7rEUS-  
type,3 JffjGf-o  
real,1 lq2Ah=FuN  
mat,3 . J"g.Q  
nsla,s,1 T;6 VI|\  
esurf,all QR'yZ45n4  
esel,s,type,,3 !<!5;f8  
eplot L_fu<W  
/psymb,esys,1 }6#lE,\lM  
allsel 5<o8prtB  
save F Qtlo+3  
bn`1JI@S4  
!侧面施加对称约束 xyGk\= S  
allsel rLMjN#`^  
CSYS,0 H{*~d+:ol  
asel,s,loc,x,0 H,r>@Y  
da,all,symm f.?p"~!  
asel,s,loc,y,0  3*@ sp  
da,all,symm r^3QDoy  
Xg>nb1e  
allsel,all ~T;FO
B%w  
!asel,s,,,2!桩顶面2未约束
sSVgDQ~q  
!da,all,all O*MC"%T  
asel,s,,,6 }UwDHq=  
da,all,ALL l{U-$}  
!约束外侧面 O 0Fw!IQk  
asel,S,,,8 W5a)`
%H  
NSLA,S,1 I[|5 DQ  
CSYS,1 rCGyr}(NC  
NROTAT,ALL HCP'V  
D,ALL,UX #]nH
$Kq  
ALLSEL nSxFz!  
SAVE l7G&[\~  
finish!接触对定义结束 o&2(xI2  
!求解过程 x5q5<-#  
/solu KiXXlaOs  
'J+dTs;0  
!约束结束! B j!{JcM-^  
!定义荷载步 "S6'<~s  
!第一步荷载，重力荷载。 g_z%L?N  
acel,0,0,9.8 5mNd5IM  
nlgeom,on ! 打开大变形选项 &WW|! 6  
autot,on ! 打开自动时间步 '1?
b?nVo  
time,1 ve|:z  
nsub,10 _jmkAmeu  
outres,all,all ?m3,e&pB5  
allsel 8BnI0l=\  
SOLVE yp%7zrU  
!第二部荷载,拉动桩 lp`raNNo  
nsubst,100,500,5 #7S[Ch}O  
outres,all,all 5&5 x[S8  
time,101 VEAf,{)Q  
nlgeom,on ! 打开大变形选项 V:?exJg9  
autot,on ! 打开自动时间步 \iFh-?(  
nsel,s,loc,z,7.1 STMc@MeZU_  
D,ALL,UZ,0.1 ?nOul}y/  
--SlxV/x  
lnsrch,on !自动线性搜索 n6T@A;_g  
pred,on !打开预测求解器 JGH;&UYP  

J!sIxwF  
allsel <u\j4<p  
SOLVE

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