桩从土中拔出的数值计算 请指教 [复制链接]

我做的是一个桩土分析问题，想计算一下桩土相对滑动时的摩擦力，采用的是1/4三维实体模型，计划分为两个荷载步，第一步利用土体自重挤住桩，第二步把桩从土中拔出。 RF
U(wek  
第一步可以算，第二步却不行，命令流如下，请大师们不吝指教，不胜感激。 z| Hl*T  
FINI ; =ai]AYW  
/CLEAR F42r]k  
/filname,pile-soil contact ;]D@KxO$dJ  
/title，pile ;cH|9m:Y  
/units,si !国际单位制 lbZ,?wm  
fe=0.1 !摩擦系数 j7K9T  
fkn=10 !法向接触刚度 M}k )Ep9  
fkt=3!切向接触刚度 5?r#6:(yI  
ftoln=1 !初始渗透因子 hPE#l?H@A  
fact=2!静动摩擦系数 AU)"L_ i}  
/prep7 Am'5|  
et,1,solid45 dQVV0)z  
mp,ex,1,7.2e10         !桩的弹性模量，根据管壁厚度折算 mT <4@RrB  
mp,nuxy,1,0.3           !桩的泊松比 BNw^ _j1  
mp,dens,1,2000         !桩的密度 #I|Vyufw  
!mp,mu,1,fe             !桩体性质 (YVl5}V  
mp,ex,2,2e8           !土的弹性模量，经验取值 .]6_  
mp,nuxy,2,0.3          !土的泊松比 .h~M&d!  
mp,dens,2,2000          !土的饱和密度   zmQ V6o=k  
tb,dp,2 UTGR{>=>  
tbdata,1,0,30      !粘聚力c为0,摩擦角为30度 wj!p6D;;S  
!创建桩模型 M
CD]
n  
cylind,0,0.1,-0.1,7.1,0,90  kj~)#KDN  
!创建土体模型 &?IOrHSv!  
cylind,0.1,0.6,0,7,0,90 odMjxWY  
!桩体划分网格 _A|1_^[G(  
allsel,all s-Q-1lKV,  
vsel,s,,,1 *v)JX _  
aslv,s  VT96ph  
lsla,s Z$ Fh4  
lesize,1,,,6 [yM{A<\L  
lesize,2,,,6 c[}h( jkP  
lesize,3,,,6 sesr`,m.,  
lesize,4,,,6 D|/Azy.[  
lesize,5,,,6 .7++wo!,  
lesize,6,,,6 aR6?+`6<  
lesize,7,0.05 )tl=tH/$  
lesize,8,0.05 zNofI$U  
lesize,9,0.05 O+%Y1=S[WQ  
type,1 %H75u
6  
mat,1 ol^V@3[<  
mshape,0,3d ;2q;RT`h  
mshkey,1 #Z;ziM:  
vmesh,all fQi7e5  
vsel,s,,,2 ！土层，划分网格 VQ5
T$,&  
vplot SL" ;\[uI  
aslv,s QMea2q|3$  
lsla,s %_;q<@9)  
lplot IzF7W?k  
lesize,10,,,6 U>7"BpC  
lesize,11,,,6 D-\\L[  
lesize,12,,,6 mVfg+d(  
lesize,13,,,6 rVowH
P  
lesize,14,,,6 x.qn$?3V]  
lesize,15,,,6 EKZ40z`  
lesize,16,,,6 '31pb9@fH  
lesize,17,,,6 1fM=>Z  
lesize,18,0.05 z)xGZ*{=  
lesize,19,0.05 ;Xg6'yxJ  
lesize,20,0.05 Un@B D}@\  
lesize,21,0.05 `Jhu&MWg  
type,1 S!A)kK+  
mat,2 3.U5Each-  
mshape,0,3d 6[,*2a8  
mshkey,1 "(5A5>  
vmesh,all xfCq;?MupW  
!创建接触对 {LYA?w^GT  
!用接触向导建立接触对 3GUO  
allsel u6IEBYG ((  
et,2,170 zk=5uKcPE  
et,3,173 Iz,a Hrq  
r,1 ,u^i0uOg  
real,1 $J"}7+  
keyopt,3,9,0    !消除初始渗透 wr,X@y%(!  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 Pwf2dm$,+  
桩侧，建立目标面 c) Zid1  
r,1,,,fkn,ftoln, , , G 5w:  
rmore,,,,,,fkt _;3xG0+  
RMORE,,,, ,, |UZPn>F~  
RMORE,,,fact,,,   :]iV*zo_
 
mp,mu,3,fe *i|O!h1St  
!创建接触对 z #c)Q  
vsel,s,,,1 @;2,TY>Di  
asel,s,,,3 EzqYHY+_r  
aplot 0M(\xO  
cm,target,area eyBLgJt8P  
type,2 "^VKs_U8o  
real,1 |0:<Z(  
mat,3 S~0 mY} m  
nsla,s,1 /*!K4)$-*2  
esurf,all ; S7 %  
esel,s,type,,2 b/cc\d<  
eplot XSC._)ztEE  
/psymb,esys,1 NF&R}7L  
allsel!很重要 gfG Mu0FjB  
!定义接触面 G/>upnA{w  
vsel,s,,,2 y3@R>@$  
vplot Q4\EI=4P]  
asel,s,,,9 QyQ&xgS  
aplot "f/lm 2<  
cm,contact,area /{ Lo0  
type,3 .vYU4g]  
real,1 \M(0@#-$C  
mat,3 (Em^qN  
nsla,s,1 &Q 7Q1`S  
esurf,all JYA$_T
  
esel,s,type,,3 vggyQf%  
eplot fa9c!xDt  
/psymb,esys,1 GhpVi<FL  
allsel lgFA}p@  
save If>k~aL7I  
!侧面施加对称约束 f{MXH&d 1\  
allsel QxG:NN;jW  
CSYS,0 8sjAr.iT.  
asel,s,loc,x,0 ts/rV#s~  
da,all,symm y<jW7GNt  
asel,s,loc,y,0 
"4"gHs  
da,all,symm f(!:_!m*  
allsel,all vp7J';  
!asel,s,,,2!没有约束顶面2 .<566g}VP  
!da,all,all xU+
c?OLi  
asel,s,,,6 Y ]&D;w  
da,all,ALL Uu ~BErEC  
!约束外侧面 .E}fk,hLB  
asel,S,,,8 _]:wltPv  
NSLA,S,1 k2DT+}u7G  
CSYS,1 1uBnU2E  
NROTAT,ALL hvCX,^LoJ  
D,ALL,UX！土侧采用的径向约束。 E uxD,(  
ALLSEL s"*ZQ0OaD  
SAVE NF+^  
finish!接触对定义结束 TOmq2*,/  
!求解过程 h}r64<Y2{  
/solu %K[_;8  
!约束结束! ;zVtJG`  
!定义荷载步 }qU(G3  
!第一步荷载，重力荷载。 l2Z!;Wm(  
acel,0,0,9.8 MebLY $&8  
nlgeom,on ! 打开大变形选项 .vHSKd{  
autot,on ! 打开自动时间步 &?C% -"|c  
time,1 gz8>uGx&V!  
nsub,10 `Mx&,;x  
outres,all,all '&?cW#J?  
allsel L5$r<t<  
SOLVE :4RD.l  
!第二部荷载,拉动桩 9NXf~-V-  
nsubst,100，500，5 L:UJur%  
outres,all,all t;X  !+  
time,101 AX=$r]
_  
nlgeom,on ! 打开大变形选项 x9S~ns+r  
autot,on ! 打开自动时间步  ;LEO+,6  
yg34b}m{  
nsel,s,loc,z,7.1 uFPJ}m[>5  
D,ALL,UZ,0.1 {jB& e,  
lnsrch,on !自动线性搜索 bV c"'RQ  
pred,on !打开预测求解器 a{*r^m'N  
allsel =_~bSEqyRI  
SOLVE /J{ e_a  
b#\i]2b:  

我的qq号 57979269 知道的请指点 不胜感激

这个问题估计你是老大吧

让我试一试，应该是可以做的，估计是你的参数设置有问题。