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

大家好，我做的是一个桩土接触分析问题，想把一个桩从土中拔出，施加两个荷载步，第一步，施加重力荷载，第二步把桩从中拔出，在运行第二步时出了问题，算不下去，请大师不吝指教 ，不胜感激。命令流如下： !Pfr,a  
FINI Pw`8Wj  
/CLEAR w;:*P  
j[J-f@F \Y  
/filname,pile-soil contact #r~#
I}U  
/title，pile (m(JK^  
!/config,nres,2000000 ">,|V-H  
/units,si !国际单位制 A&Usddcp  
jZkcBIK2  
fe=0.1 !摩擦系数 H;k~oIsk  
fkn=10 !法向接触刚度 Ww+IWW@  
fkt=3!切向接触刚度 2*l/3VW  
ftoln=1 !初始渗透因子 6Vnsi%{  
fact=2!静动摩擦系数 fW1CFRHH  
:vQrOn18p  
/prep7 }?_?V&K|  
qvKG-|j  
et,1,solid45 _FU_Ubkr  
A>
;bHf@  
mp,ex,1,7.2e10 !桩的弹性模量，根据管壁厚度折算 (Y?gn)*t  
mp,nuxy,1,0.3 !桩的泊松比 .|>3k'<l  
mp,dens,1,2000 !桩的密度 goOCu  
!mp,mu,1,fe !桩体性质 bSi%2Onj  
x,@B(9No  
mp,ex,2,2e8 !土的弹性模量，经验取值 U-(01-  
mp,nuxy,2,0.3 !土的泊松比 S3*`jF>q  
mp,dens,2,2000 !土的饱和密度 pG^  
tb,dp,2 vDhh>x(  
tbdata,1,0,30 !粘聚力c为0,摩擦角为30 rCbDu&k]  
jB Z&Ad@e  
!创建桩模型 4=.89T#<  
cylind,0,0.1,-0.1,7.1,0,90 ^Hnb}L  
!创建土体模型 1N#| }ad  
cylind,0.1,0.6,0,7,0,90 H0gbSd+  
x$.^"l-vX  
!桩体划分网格 )9'K($  
allsel,all :tB1D@Cb6  
vsel,s,,,1 c&
?m>2^6  
aslv,s %XoiVl
T@:  
lsla,s (!7sE9rP  
lesize,1,,,6 r|fL&dtr  
lesize,2,,,6 S\YTX%Xm}  
lesize,3,,,6 U}e!Wjrc  
lesize,4,,,6 PI:4m%[  
lesize,5,,,6 .*?wF  
lesize,6,,,6 FQ7T'G![  
lesize,7,0.05 s_OF(o  
lesize,8,0.05 5z8d} I  
lesize,9,0.05 Fg5kX  
type,1
kYqU9cB~  
mat,1 IS{wtuA.  
mshape,0,3d i$:*Pb3mV  
mshkey,1 'qb E=  
vmesh,all F^t DL:  
Vvn
2 Ep  
vsel,s,,,2 ！土层，划分网格  gmO!  
vplot JPc+rfF  
aslv,s R'bTN|Cq  
lsla,s ;4^Rx  
lplot \378rQU  
lesize,10,,,6 +`7i'ff  
lesize,11,,,6
:uq\+(9  
lesize,12,,,6 o WrKM  
lesize,13,,,6 `iAF3:  
lesize,14,,,6 5T_n %vz  
lesize,15,,,6 nwB_8mN|  
lesize,16,,,6 Pw7]r<Q  
lesize,17,,,6 nQX:T;WL@  
lesize,18,0.05 q77;ZPfs8  
lesize,19,0.05 "3)C'WlEy/  
lesize,20,0.05 x=hiQ>BIO0  
lesize,21,0.05 U0 Yll4E  
type,1 b8`)y<7  
mat,2 M=.n7RY-  
mshape,0,3d xd q?/^E  
mshkey,1 g:'xae/]S  
vmesh,all G~^r)fm_  
< Mn ;  
!创建接触对 SIF/-{i(X  
!用接触向导建立接触对 fNli  
allsel 6y%qVx#!  
et,2,170 (lBCO?`fx  
et,3,173 dUeN*Nq&(,  
r,1 BOb">6C  
real,1 %Q__!D[  
keyopt,3,9,0 !消除初始渗透 |"X*@s\'  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 U3ADsdn  
uB]7G0g:  
桩侧，建立目标面 ~:rl=o}  
r,1,,,fkn,ftoln, , , Z58X5"  
rmore,,,,,,fkt <y2U3;t  
RMORE,,,, ,, Ga-k  
RMORE,,,fact,,, F 5bj=mI  
mp,mu,3,fe b@gc{R}7  
}l(&}#dY  
!创建接触对 ?mwt~_s9  
vsel,s,,,1 w=0(<s2  
asel,s,,,3 a?I= !js  
aplot j\[dx^\=  
cm,target,area {NmWQyEv  
type,2 ajbA\/\G;  
real,1 !=*g@mgF  
mat,3 o8V5w!+#  
nsla,s,1 P%V'4pc  
esurf,all zsEc(  
esel,s,type,,2 *)$Uvw E  
eplot <.x{|p  
/psymb,esys,1 Z, zWuE3  
allsel!很重要 u:6Ic)7'  
|sJ[0z  
!定义接触面 vjbASFF0=  
vsel,s,,,2 ,8S/t+H  
vplot .KB^3pOpx  
asel,s,,,9 [N-Di"  
aplot }Sm(]y  
cm,contact,area lK?uXr7^  
type,3 :T^a&)aL%  
real,1 !fV+z%:  
mat,3 \#2Z)Kz  
nsla,s,1 0#7>o^2  
esurf,all |{z:IQLv  
esel,s,type,,3 a5dLQxb  
eplot *<$*"p  
/psymb,esys,1 Q,,e+exbb5  
allsel G)AqbY  
save MD}w Y><C  
!j
8FIY'[  
!侧面施加对称约束
@+&LYy72  
allsel 7cuE7"  
CSYS,0 >2y':fO  
asel,s,loc,x,0 a{L%7  
da,all,symm G*?8MTP8![  
asel,s,loc,y,0 \0^Kram>
 
da,all,symm b<tNk]7  
N~nziY*C,*  
allsel,all paA(C|%{  
!asel,s,,,2!桩顶面2未约束 +C^nO=[E  
!da,all,all q\9JgD)  
asel,s,,,6 z 4e7PW|  
da,all,ALL =Pyj%4Rs  
!约束外侧面 w49t9~  
asel,S,,,8 Yj<a" Gr4[  
NSLA,S,1 k90YV(  
CSYS,1 I {SjlN}d  
NROTAT,ALL pE3?"YO  
D,ALL,UX Ij7p'a  
ALLSEL i7CX65&b  
SAVE H9Gh>u]}  
finish!接触对定义结束 U"~>jZKk  
!求解过程 ^Y?k0z  
/solu G~]Uk*M q  
B5,N7z34F  
!约束结束! H-!,yte  
!定义荷载步 9sM!`Lz{  
!第一步荷载，重力荷载。 .y'>[  
acel,0,0,9.8 dUD[e,?  
nlgeom,on ! 打开大变形选项 IY1//9  
autot,on ! 打开自动时间步 zF@/K`  
time,1 O"+gQXe  
nsub,10 5p,RI&nlN  
outres,all,all &.F4b~A7  
allsel SjK  
SOLVE h<h%*av|  
!第二部荷载,拉动桩 %6t
:(z  
nsubst,100,500,5 av(6wht8  
outres,all,all j\ZXG=j  
time,101 f'F?MINJP  
nlgeom,on ! 打开大变形选项 p0]=QH  
autot,on ! 打开自动时间步 H
kg2P,2  
nsel,s,loc,z,7.1 qRu~$K  
D,ALL,UZ,0.1 qfX6TV5J}!  
mupT<_Y  
lnsrch,on !自动线性搜索 ynp8rf  
pred,on !打开预测求解器 ,10=  
0RzEY!9g+  
allsel 6u?>M9  
SOLVE

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