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

大家好，我做的是一个桩土接触分析问题，想把一个桩从土中拔出，施加两个荷载步，第一步，施加重力荷载，第二步把桩从中拔出，在运行第二步时出了问题，算不下去，请大师不吝指教 ，不胜感激。命令流如下： LF
{8hC[  
FINI !F Zg' 9  
/CLEAR E8lq2r=  
&-
My[t  
/filname,pile-soil contact U9\w)D|+eE  
/title，pile +sq_fd ;'D  
!/config,nres,2000000 eT2Tg5Etc  
/units,si !国际单位制 bq8h?Q  
n-GoG(s..b  
fe=0.1 !摩擦系数 IO2@^jup  
fkn=10 !法向接触刚度 T;,
,!  
fkt=3!切向接触刚度 CLuQ=-[
|  
ftoln=1 !初始渗透因子 r0pwKRE~t  
fact=2!静动摩擦系数 k2#|^N  
p%#'`*<a_  
/prep7 Nn?$}
g  
D mky!Cp
 
et,1,solid45 ]1[:fQF7/L  
IF6-VFY:6  
mp,ex,1,7.2e10 !桩的弹性模量，根据管壁厚度折算
\IqCC h  
mp,nuxy,1,0.3 !桩的泊松比 q(2ZJn13f  
mp,dens,1,2000 !桩的密度 ^nFP#J)_5  
!mp,mu,1,fe !桩体性质 5Z:qU{[  
HHs!6`R$0c  
mp,ex,2,2e8 !土的弹性模量，经验取值 Ue>{n{H"y  
mp,nuxy,2,0.3 !土的泊松比 B)P]C5KRD  
mp,dens,2,2000 !土的饱和密度 Q/h-Khm
z  
tb,dp,2 lPtML<a  
tbdata,1,0,30 !粘聚力c为0,摩擦角为30 eq<!  
(Y^tky$9  
!创建桩模型 ]9W7]$  
cylind,0,0.1,-0.1,7.1,0,90 Pj!f^MN  
!创建土体模型 0fqycGSmU  
cylind,0.1,0.6,0,7,0,90 /cvMp#<]  
`z7,HJ.0c  
!桩体划分网格 VZ 7(6?W  
allsel,all wAL}c(EHO  
vsel,s,,,1 *!dA/sid  
aslv,s (I+e@UUiL  
lsla,s *=T(ncR['  
lesize,1,,,6 hR.vJ2oa  
lesize,2,,,6 h 'Hnq m  
lesize,3,,,6 M23r/eg]  
lesize,4,,,6 9}4~3_gv;M  
lesize,5,,,6 kZi/2UA5Z  
lesize,6,,,6 (I35i!F+tY  
lesize,7,0.05 MhB>bnWXR  
lesize,8,0.05 v?%vB#A^  
lesize,9,0.05 |WOc0M[U  
type,1 Cbm^: _LR  
mat,1 Imq-5To#  
mshape,0,3d Ji6.
-[:  
mshkey,1 $3%+N|L  
vmesh,all lJT"
aXt'M  
| :i
d/  
vsel,s,,,2 ！土层，划分网格 L1y71+iqU  
vplot 33&\E- Q>  
aslv,s j5]ul!ji  
lsla,s - leYR`P  
lplot YhC|hDC  
lesize,10,,,6 6GqC]rd*:  
lesize,11,,,6 f/iMI)J  
lesize,12,,,6 3=*ur( Qy  
lesize,13,,,6 )8S
WU)/  
lesize,14,,,6 ,(f W0d#  
lesize,15,,,6 j0(jXAc;UB  
lesize,16,,,6 ;P/ 4.|<  
lesize,17,,,6 k]`-Y E  
lesize,18,0.05 @7V~CNB+  
lesize,19,0.05 j['B9
vG  
lesize,20,0.05 GQQp(%T  
lesize,21,0.05 *JQ*$$5  
type,1 =q`T|9v  
mat,2 evNe6J3  
mshape,0,3d )}t't"  
mshkey,1 m4Ue)  
vmesh,all ;K%/sIIke  
 _+(@?  
!创建接触对  ts=:r  
!用接触向导建立接触对 Chx+p&!  
allsel z0#2?o  
et,2,170 6<R[hIWpZ}  
et,3,173 )feZ&G]  
r,1 K>b4(^lf  
real,1 D!.1R!(Z  
keyopt,3,9,0 !消除初始渗透 %H2ios[UO  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 6]kBG?m0  
a60rJ#GD  
桩侧，建立目标面 ]^>:)q  
r,1,,,fkn,ftoln, , , ?:-:m'jdU  
rmore,,,,,,fkt E-FR w  
RMORE,,,, ,, '3WtpsKA  
RMORE,,,fact,,, 7mM
MVz2  
mp,mu,3,fe cDE5/!  
qMA-#  
!创建接触对 I~@8SSO,vH  
vsel,s,,,1 7')W+`o8eL  
asel,s,,,3 X{OWDy  
aplot slUi)@b  
cm,target,area SgehOu  
type,2 UxzF5V5  
real,1 di0@E<@1:  
mat,3 w`V6vYd@  
nsla,s,1 js@L%1r#L  
esurf,all `[U.BVP'  
esel,s,type,,2 w D r/T3  
eplot ]
sk=V.GGQ  
/psymb,esys,1 ?+.C@_QZQ  
allsel!很重要 GZNN2 '  
.9PT)^2  
!定义接触面 0<NS1y  
vsel,s,,,2 '?L^Fa_H  
vplot %A=/(%T>  
asel,s,,,9 O=(F46 M  
aplot H"6x/&s.=k  
cm,contact,area 7bihP@I!  
type,3 f:<BUqa  
real,1 V|$PO Qa3  
mat,3 31alQ\TH  
nsla,s,1 o5NmNOXm  
esurf,all e5P9P%1w  
esel,s,type,,3 hr6j+p:  
eplot ^MF 2Q+  
/psymb,esys,1 X.k8w\~  
allsel iaeNY;T  
save %oTBh*K'o  
kwc Cf2  
!侧面施加对称约束 F}ATY!  
allsel f!w/zC .  
CSYS,0 >19s:+  
asel,s,loc,x,0 A@)Q-V8*9s  
da,all,symm
to
</  
asel,s,loc,y,0 "g&f:[a/  
da,all,symm ]h(Iun  
?{J1Uw<  
allsel,all ,.cR@5qI  
!asel,s,,,2!桩顶面2未约束 c]aU}[s1  
!da,all,all m{ !$_z8:  
asel,s,,,6 pF-_yyQ  
da,all,ALL 4=Ru{ewRV  
!约束外侧面 fJc(  
asel,S,,,8 az0=jou<Zl  
NSLA,S,1 EOXkMr  
CSYS,1 BoYY^ih
 
NROTAT,ALL 0't)-Pj+,  
D,ALL,UX Yr_B(n  
ALLSEL D5@=#/
?*  
SAVE nsU7cLf"^V  
finish!接触对定义结束 s)r!3HS  
!求解过程 `D~oY=  
/solu _Fa\y ZX  
G&{yM2:E  
!约束结束! r*vh3.Agl  
!定义荷载步 @M4c/k}  
!第一步荷载，重力荷载。 &'W7-Z\j-  
acel,0,0,9.8 } u;{38~  
nlgeom,on ! 打开大变形选项 |ffHOef  
autot,on ! 打开自动时间步 "gM!/<~  
time,1 zMbFh_dcq  
nsub,10 qm!oJL  
outres,all,all
Q
jY}$  
allsel +?U[362>  
SOLVE cLf90|YFp  
!第二部荷载,拉动桩 p;}`PW  
nsubst,100,500,5 `J,>#Y6(J  
outres,all,all 0<!9D):Bb  
time,101 )Y)_T&O  
nlgeom,on ! 打开大变形选项 #RR;?`,L}  
autot,on ! 打开自动时间步 it\$Pih]  
nsel,s,loc,z,7.1 |JIlp"[  
D,ALL,UZ,0.1 KMIe%2:b5  
SED52$zA  
lnsrch,on !自动线性搜索 ) ~=pt&+  
pred,on !打开预测求解器 zU&Iy_Ke.  
+ m-88  
allsel J\*uW|=F  
SOLVE

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