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

我做的是一个桩土分析问题，想计算一下桩土相对滑动时的摩擦力，采用的是1/4三维实体模型，计划分为两个荷载步，第一步利用土体自重挤住桩，第二步把桩从土中拔出。 zuk"  
第一步可以算，第二步却不行，命令流如下，请大师们不吝指教，不胜感激。 6zv;lx0<D&  
FINI _O"L1Let  
/CLEAR uC- A43utv  
/filname,pile-soil contact qw5&Y$((  
/title，pile IKGTsA;  
/units,si !国际单位制  oHOW5  
fe=0.1 !摩擦系数 {K<uM'ww>  
fkn=10 !法向接触刚度 3mk=ZWwv  
fkt=3!切向接触刚度 5dqQws-,?1  
ftoln=1 !初始渗透因子 d3Y(SPO  
fact=2!静动摩擦系数 -MDOZz\  
/prep7 _2]O^$L  
et,1,solid45 O+p]3u  
mp,ex,1,7.2e10         !桩的弹性模量，根据管壁厚度折算 pF|8OB%  
mp,nuxy,1,0.3           !桩的泊松比 qZXyi'(d  
mp,dens,1,2000         !桩的密度 @+iC/
  
!mp,mu,1,fe             !桩体性质 Cy<T Vk8  
mp,ex,2,2e8           !土的弹性模量，经验取值 : #3OcD4  
mp,nuxy,2,0.3          !土的泊松比 m dC`W&r  
mp,dens,2,2000          !土的饱和密度   y!SF/i?Py  
tb,dp,2 /%62X{=>
;  
tbdata,1,0,30      !粘聚力c为0,摩擦角为30度 D!X>O}
  
!创建桩模型 D*CIE\+  
cylind,0,0.1,-0.1,7.1,0,90 S|~i>  
!创建土体模型 nrF%wH/5  
cylind,0.1,0.6,0,7,0,90 c/b%T  
!桩体划分网格 UHr{  
allsel,all tp#Z@5=  
vsel,s,,,1 8oa)qaG1  
aslv,s
Tku/OG'  
lsla,s GXarUjs  
lesize,1,,,6 ~l(G6/R  
lesize,2,,,6 'CS^2Z  
lesize,3,,,6 qo61O\qm  
lesize,4,,,6 xR*5q1j  
lesize,5,,,6 ?hxK/%)  
lesize,6,,,6 4sj9Z:  
lesize,7,0.05 +Y^-e.UO  
lesize,8,0.05 ?68$3;  
lesize,9,0.05 ?F)_T  
type,1 Uc\|X;nkRk  
mat,1 A *a{  
mshape,0,3d A(q~{  
mshkey,1 fB_4f{E  
vmesh,all ;:AG2z
E!  
vsel,s,,,2 ！土层，划分网格 |)28=Z|Z  
vplot ^Yf3"D?&  
aslv,s 9@(O\xr  
lsla,s h~,x7]
w6  
lplot -"2<h:#  
lesize,10,,,6 d|>9rX+f  
lesize,11,,,6 tU9rCL:P  
lesize,12,,,6 L.SDMz  
lesize,13,,,6 lO551Y^  
lesize,14,,,6 T
{hyt  
lesize,15,,,6 |}M~kJ)  
lesize,16,,,6 <dq,y>  
lesize,17,,,6 GA'*5
8  
lesize,18,0.05 U#8\#jo  
lesize,19,0.05 )W[KD,0+j  
lesize,20,0.05 "CIpo/ebL  
lesize,21,0.05 ;X+G6F'  
type,1 bq c;.4$  
mat,2 FI3sLA  
mshape,0,3d K.:6YXVs<  
mshkey,1 d8xk&za  
vmesh,all A28ZSL  
!创建接触对 j'G"ZPw1  
!用接触向导建立接触对 '}wG"0  
allsel O4lxeiRgC  
et,2,170 " 7l jc  
et,3,173 r0~7v1rG  
r,1 ~z\pI|DQ  
real,1 1joc<EI  
keyopt,3,9,0    !消除初始渗透 RH<@c^ S  
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 nvU+XCx  
桩侧，建立目标面 Q'FX:[@x-S  
r,1,,,fkn,ftoln, , , 3m-edpH  
rmore,,,,,,fkt ?whRlh  
RMORE,,,, ,, 3c1o,2  
RMORE,,,fact,,,   ({*.!ty  
mp,mu,3,fe 0|],d?-h  
!创建接触对 ,$hQ(yF  
vsel,s,,,1 .a 'ETNY:>  
asel,s,,,3 ~*66 3pA  
aplot B ]|5?QP-  
cm,target,area OU0xZ=G  
type,2 I>N-95  
real,1 pKOT Qf  
mat,3 )OARO  
nsla,s,1 T+&x{+gZ  
esurf,all t:?<0yfp&  
esel,s,type,,2 F!z! :yp  
eplot V/QTYy1  
/psymb,esys,1 /ChJ~g
"  
allsel!很重要 Y}V)4j  
!定义接触面 [+%p!T  
vsel,s,,,2 5*G8W\ $  
vplot sRkz WMl  
asel,s,,,9 9b{g+lMZo  
aplot ')5jllxv  
cm,contact,area y>)mSl@1y  
type,3 9`nP
(~  
real,1 #Z2>TN  
mat,3 DI$mD{  
nsla,s,1 /u" cl2|  
esurf,all b?+Yo>yF8  
esel,s,type,,3 9*GL@_c  
eplot $5q{vy  
/psymb,esys,1 B`RW-14g  
allsel M]8eW  
save [#H
8=  
!侧面施加对称约束 )w}*PL  
allsel : L6-{9$  
CSYS,0
SdEb[  
asel,s,loc,x,0 30gZ_8C>}  
da,all,symm N7j]yvE  
asel,s,loc,y,0 2uI`$A:  
da,all,symm (&N$W&  
allsel,all %k1q4qOG]^  
!asel,s,,,2!没有约束顶面2 O"2wV +9  
!da,all,all yK @X^jf  
asel,s,,,6 Ggxrj'r  
da,all,ALL &9jUf:gJ0  
!约束外侧面 nU">> 1!U  
asel,S,,,8 `9G$p|6  
NSLA,S,1 i?P]}JENM  
CSYS,1 `#IT24!  
NROTAT,ALL *;\ K5  
D,ALL,UX！土侧采用的径向约束。 ?(`nBlWQ5  
ALLSEL I?:+~q}lZr  
SAVE hg86#jq%  
finish!接触对定义结束 #!<+:y'S?  
!求解过程 ]0/~6f  
/solu 34O+#0<y~  
!约束结束! FZ}C;yUPD  
!定义荷载步 eog,EP"a8Y  
!第一步荷载，重力荷载。 I5|S8d<  
acel,0,0,9.8 aaqjE  
nlgeom,on ! 打开大变形选项 OFPd6,(E  
autot,on ! 打开自动时间步 $[^ KCNB  
time,1 iz,q8}/(  
nsub,10 `OF;>u*:  
outres,all,all Qbe{/  
allsel 0ghwFo  
SOLVE ^5zS2nm  
!第二部荷载,拉动桩 tM?I()Y&P  
nsubst,100，500，5 Z?G3d(YT  
outres,all,all v#AO\zYKd  
time,101 `P)64So-1  
nlgeom,on ! 打开大变形选项 .q;RNCUt  
autot,on ! 打开自动时间步 @Ig,_i\UY:  
K\%"RgF@&  
nsel,s,loc,z,7.1 *FG4!~<e  
D,ALL,UZ,0.1 ,_X,V!  
lnsrch,on !自动线性搜索 L>R!A3G1  
pred,on !打开预测求解器 hTNYjXj  
allsel }O-%kl  
SOLVE brhJ&|QDE  
hzf}_
1  

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

这个问题估计你是老大吧

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