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

大家好，我做的是一个桩土接触分析问题，想把一个桩从土中拔出，施加两个荷载步，第一步，施加重力荷载，第二步把桩从中拔出，在运行第二步时出了问题，算不下去，请大师不吝指教 ，不胜感激。命令流如下： J0eJRs  
FINI A&:i$`m,  
/CLEAR ief~*:5  
8 FqhSzw  
/filname,pile-soil contact M#PutrH  
/title，pile fZ*LxL  
!/config,nres,2000000 [4_JK  
/units,si !国际单位制 Vif0z*\e{  
_FJ,, /~  
fe=0.1 !摩擦系数 aj71oki)  
fkn=10 !法向接触刚度 GkT:7`|C  
fkt=3!切向接触刚度 cxNb!G  
ftoln=1 !初始渗透因子  d Xiv8B1  
fact=2!静动摩擦系数 <e s>FD  
UY!N"[&  
/prep7 OZ/!=;  
lhk[U!>#  
et,1,solid45 Y8Bc &q}  
I,:R~^qJ8v  
mp,ex,1,7.2e10 !桩的弹性模量，根据管壁厚度折算 (np %urx!  
mp,nuxy,1,0.3 !桩的泊松比 h5:>o  
mp,dens,1,2000 !桩的密度 mF?GQls`  
!mp,mu,1,fe !桩体性质 k[6xuyY]  
z DP  
mp,ex,2,2e8 !土的弹性模量，经验取值 >brf7h  
mp,nuxy,2,0.3 !土的泊松比 i"h '^6M1  
mp,dens,2,2000 !土的饱和密度 QX'EMyK$  
tb,dp,2 huu v`$~y  
tbdata,1,0,30 !粘聚力c为0,摩擦角为30 \GR M,c  
yJ
heni  
!创建桩模型 (L{Kg U&{$  
cylind,0,0.1,-0.1,7.1,0,90 \1AtBc&  
!创建土体模型 t~luBUF  
cylind,0.1,0.6,0,7,0,90 RO([R=.`/  
d=TZaVL$$  
!桩体划分网格 Nf#8V|  
allsel,all PjH[8:,  
vsel,s,,,1 "w)Y0Qq*z  
aslv,s (|6Y1``  
lsla,s ,+u.FQv~  
lesize,1,,,6 yU-^w^4  
lesize,2,,,6 sg.8Sd"]7  
lesize,3,,,6 8$|<`:~J  
lesize,4,,,6 Y7vTs
eq  
lesize,5,,,6 "k8Yc<`u  
lesize,6,,,6 kHO2&"6  
lesize,7,0.05 wIrjWU2  
lesize,8,0.05 |M*jo<C  
lesize,9,0.05 GwOn&EpY!  
type,1 D^?jLfW8  
mat,1 23lLoyN  
mshape,0,3d *Cdw"n  
mshkey,1 "Ih>>|r  
vmesh,all NF}QQwG3  
LqH<HGMFD  
vsel,s,,,2 ！土层，划分网格 fZ9EE3  
vplot ;g3z?Uz)  
aslv,s Kk}|[\fW  
lsla,s AoaRlk-#  
lplot u,zA^%  
lesize,10,,,6 uhyw?#f  
lesize,11,,,6 P87!+pB(  
lesize,12,,,6 F:PaVr3q  
lesize,13,,,6 pjC2jlwm*  
lesize,14,,,6 Q}KNtNCpx  
lesize,15,,,6 ^w0V{qF{  
lesize,16,,,6 |H:JwxH  
lesize,17,,,6
4%8}vCs  
lesize,18,0.05 Ywj=6 +;  
lesize,19,0.05 ")/TbTVu  
lesize,20,0.05 E<[_L!2  
lesize,21,0.05 y_^w|  
type,1 0_}^IiG  
mat,2 vybQ}dscn  
mshape,0,3d `EV"
/&`  
mshkey,1 ? m$uqi  
vmesh,all G}o?lo\#h  
M18H1e@Al  
!创建接触对 \9U4V>p  
!用接触向导建立接触对 g/)$-Z)Nu  
allsel ;(a\F  
et,2,170 c0J=gZiP  
et,3,173 2?",2x09  
r,1 2v!ucd}  
real,1 }pE8G#O&  
keyopt,3,9,0 !消除初始渗透 'Zq$W]i  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 s m42  
'~dE0ohWb  
桩侧，建立目标面 {K{&__Nk  
r,1,,,fkn,ftoln, , , 8q
)wT0A~  
rmore,,,,,,fkt )+^1QL  
RMORE,,,, ,, u9
WQ0.  
RMORE,,,fact,,, `w_?9^7mH  
mp,mu,3,fe F1V[8I.0  
5e!YYt>  
!创建接触对 ;Nw.  
vsel,s,,,1
dx=\Pq  
asel,s,,,3 aR }|^ex  
aplot 2b+0}u>a  
cm,target,area #z|\AmZ\  
type,2 _1z|QC  
real,1 e=11EmN9  
mat,3 ;WQ@dC  
nsla,s,1 oB9t&yM  
esurf,all fIrl?X']  
esel,s,type,,2 _Qy3A T~  
eplot `O-LM
e  
/psymb,esys,1 | v? pS  
allsel!很重要 V\ud4  
l@+WGh  
!定义接触面 x+V;UD=mH  
vsel,s,,,2 zNE"5  
vplot B,sv! p+q5  
asel,s,,,9 MC.,n$O}6  
aplot \V@Hf"=j  
cm,contact,area s*R\!L  
type,3 zE{@
'  
real,1 .xtjB8gc  
mat,3 |aiP7C  
nsla,s,1 -.A8kJ  
esurf,all S- Mh0o"  
esel,s,type,,3 =e7,d$i  
eplot q% Eze  
/psymb,esys,1 (fgX!G[W  
allsel &"dT/5}6  
save tuA
,t  
gf]biE"k  
!侧面施加对称约束 ;!<WL@C~  
allsel xCH,d:n=  
CSYS,0 m178S3  
asel,s,loc,x,0 B;'Dh<J1
 
da,all,symm Jtj_R
l !  
asel,s,loc,y,0 fH.:#O:  
da,all,symm (7RxCo=X  
-GkNA"2M[  
allsel,all :$&%Pxm  
!asel,s,,,2!桩顶面2未约束 V/Hjd`n)`i  
!da,all,all DPvM|n`TW  
asel,s,,,6 73/kyu-0%  
da,all,ALL %H:uE*WZ  
!约束外侧面 $'3`$   
asel,S,,,8 `nccRy<l  
NSLA,S,1 2Mc/ah  
CSYS,1 [>"bL$tlo*  
NROTAT,ALL 4[za|t  
D,ALL,UX I__|+%oC  
ALLSEL .BGM1ph}~  
SAVE v*%#Fp,g8  
finish!接触对定义结束 Q
Rnkj]b  
!求解过程 jsSxjf;O  
/solu YV*s1t/  
RQI?\?o  
!约束结束! NLFSw  
!定义荷载步 -ea>}S  
!第一步荷载，重力荷载。 kLVf}J~?  
acel,0,0,9.8 jw$3cwddH  
nlgeom,on ! 打开大变形选项 E6
n3[Z  
autot,on ! 打开自动时间步 O=}4?Xv  
time,1 |pJ)w  
nsub,10 )/H=m7}1h  
outres,all,all ud K)F$7  
allsel ]gYnw;W$  
SOLVE >N>WOLbb7(  
!第二部荷载,拉动桩 0es[!  
nsubst,100,500,5 M:|Z3p K  
outres,all,all 2"lDKjj  
time,101 TEl:;4  
nlgeom,on ! 打开大变形选项 Z%~j
)  
autot,on ! 打开自动时间步 t>|Y-i3cb  
nsel,s,loc,z,7.1 vCT5do"C&  
D,ALL,UZ,0.1 <}-[9fW  
|ofegO}W7  
lnsrch,on !自动线性搜索 `k65&]&d  
pred,on !打开预测求解器 W HO;;j  
T9]|*~ ,T  
allsel ([zt}uf  
SOLVE

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