分享一个桩土三维接触的实例

yebaobao0310

桩土相互作用接触模型
/filname,pile-soil contact
/title，pile
/config,nres,2000000   
/units,si           !国际单位制
fe=0.3              !摩擦系数
fk1=10              !法向接触刚度
fk2=1               !法向接触刚度
fd=3                !切向接触刚度  
ftoln=1              !初始渗透因子
A=3.1415926*0.6*0.6   !桩的横截面积              
/prep7
et,1,solid45            
mp,mu,1,0.2      !桩体材料属性
mp,ex,1,2.5e10           
mp,nuxy,1,0.2            
mp,dens,1,2500           
mp,ex,2,2.0e8    !土层2材料属性         
mp,nuxy,2,0.4           
mp,dens,2,2000           
mp,ex,3,2.5e8     !土层3材料属性      
mp,nuxy,3,0.4            
mp,dens,3,2000           
mp,ex,4,3.0e8      !土层4材料属性      
mp,nuxy,4,0.4           
mp,dens,4,2000         
mp,ex,5,1.8e9      !岩层5材料属性     
mp,nuxy,5,0.29         
mp,dens,5,2600           
tb,concr,1           !桩体的参数
tbdata,,0.7,0.9,2.6,-1
tb,dp,2             !土层2的D-P参数     
tbdata,1,19,32,30        
tb,dp,3               !土层3的D-P参数         
tbdata,1,19,32,30
tb,dp,4             !土层4的D-P参数     
tbdata,1,19,32,30      
tb,dp,5            !土层5的D-P参数     
tbdata,1,0,35,0   
cylind,0,0.6,0,21.0,0,90        
cylind,0.6,3.0,0,3.0,0,90        
cylind,0.6,3.0,3.0,8.0,0,90     
cylind,0.6,3.0,8.0,14.0,0,90     
cylind,0.6,3.0,14.0,21.0,0,90         
cylind,0,0.6,21.0,24.0,0,90   !建立桩尖岩层几何体
cylind,0.6,3.0,21.0,24.0,0,90
!布尔运算粘结体所有的除桩之外的部分
vsel,s,,,2,7,1
vglue,2,3,4,5,6,7
allsel,all
numcmp,all
/pnum,volu,1
vplot
vlist
save
/pnum,line,1
allsel,all
vsel,s,,,1      !桩体附属性，划分网格  
aslv,s
lsla,s
lplot
lesize,1,0.1
lesize,2,0.1
lesize,3,,,6
lesize,4,,,6
lesize,5,0.1
lesize,6,0.1
lesize,7,0.5
lesize,8,0.5
lesize,9,0.5
type,1
mat,1
mshape,0,3d
mshkey,1
vmesh,all
vsel,s,,,2      !土层2附属性，划分网格  
vplot
aslv,s
lsla,s
lplot
lesize,10,,,6
lesize,11,,,6
lesize,12,,,6
lesize,13,,,6
lesize,14,,,6
lesize,15,,,6
lesize,16,,,6
lesize,17,,,6
lesize,18,0.5
lesize,19,0.5
lesize,20,0.5
lesize,21,0.5
type,1
mat,2
mshape,0,3d
mshkey,1
vmesh,all
vsel,s,,,4     !土层3附属性，并划分网格  
vplot
aslv,s
lsla,s
lplot
lesize,22,,,6
lesize,23,,,6
lesize,24,,,6
lesize,25,,,6
lesize,43,0.5
lesize,44,0.5
lesize,45,0.5
lesize,46,0.5
type,1
mat,3
mshape,0,3d
mshkey,1
vmesh,all
vsel,s,,,5    !土层4附属性，并划分网格  
vplot
aslv,s
lsla,s
lplot
lesize,26,,,6
lesize,27,,,6
lesize,28,,,6
lesize,29,,,6
lesize,47,0.5
lesize,48,0.5
lesize,49,0.5
lesize,50,0.5
type,1
mat,4
mshape,0,3d
mshkey,1
vmesh,all
vsel,s,,,6    !岩层5附属性，并划分网格  
vplot
aslv,s lsla,s
lplot
lesize,30,,,6
lesize,31,,,6
lesize,32,,,6
lesize,33,,,6
lesize,51,0.5
lesize,52,0.5
lesize,53,0.5
lesize,54,0.5
type,1
mat,5
mshape,0,3d
mshkey,1
vmesh,all
vsel,s,,,3   !桩底岩层6附属性，并划分网格  
vplot
aslv,s lsla,s
lplot
lesize,34,,,6
lesize,35,,,6
lesize,36,,,6
lesize,37,,,6
lesize,41,,,6
lesize,42,,,6
type,1
mat,5
mshape,0,3d
mshkey,1
vmesh,all
vsel,s,,,7   !桩底岩层6附属性，并划分网格  
vplot
aslv,s
lsla,s
lplot
lesize,38,,,6
lesize,57,,,6
lesize,58,,,6
lesize,41,0.5
lesize,42,0.5
lesize,55,0.5
lesize,56,0.5
type,1
mat,5
mshape,0,3d
mshkey,1
vmesh,all
!用接触向导建立接触对
et,2,170
et,3,173
keyopt,3,4,1
keyopt,3,9,0     !消除初始渗透
keyopt,3,12,2    !不分开接触，允许滑动，法向不分开
!桩侧，建立目标面
r,1,,,fk1,ftoln, , ,
rmore,,,,,,5
r,2,,,fk2,ftoln, , ,
rmore,,,,,,5
mp,mu,6,0.2
allsel,all
vsel,s,,,1
aslv,s
asel,s,,,3
/pnum,area,1
aplot
cm,target1,area
type,2
real,1
mat,6
nsla,s,1
esurf,all
! 桩侧，建立接触面  
allsel
vsel,s,,,2
vsel,a,,,4,6,1
aslv,s
asel,s,,,9
asel,a,,,22
asel,a,,,26
asel,a,,,30
aplot
cm,contact1,area
type,3
real,1
mat,6
nsla,s,1
esurf,all
!桩尖，建立目标面
allsel,all
vsel,s,,,3
aslv,s
asel,r,,,15
aplot
nsla,s,1
type,2
real,2
mat,6
cm,target2,elem
esurf,all
esel,s,type,,2  
eplot
!桩尖，建立接触面  
allsel,all
vsel,s,,,1
aslv,s
asel,s,,,2
aplot
nsla,s,1
type,3
real,2
mat,6
cm,contact2,area
esurf,all
esel,s,type,,3  
eplot /psymb,esys,1   
allsel  
gplot
save
finish
!求解过程
/solu
!定义面约束条件asel,s,loc,x,0  !侧面施加对称约束   
da,all,symm
asel,s,loc,y,0   
da,all,symm  
asel,s,loc,z,24  !底面试加约束
da,all,all
allsel,all   !土外层施加约束asel,s,,,8
asel,a,,,20
asel,a,,,24
asel,a,,,28
asel,a,,,32
aplot   
da,all,all
allsel,all
!设置非线性选项
acel,0,0,9.8      !定义重力加速度
neqit,50,         !迭代次数
nropt,Modi        ! 用改变的牛顿-拉普拉斯方程求解
nlgeom,on         ! 打开大变形选项
autot,on          ! 打开自动时间步
lnsrch,on         !自动线性搜索<BR>pred,on           !打开预测求解器
outpr,all,all
OUTRES,all,all
allsel,all
cncheck,detail   !检查接触状态cncheck,summary
/solu
CNCHECK,POST
FINISH
!载荷步,加载阶段
/solu
FirstF=2E+06/A  
LastF=3E+06/A  Increment=2E+05/A
*do,Force,FirstF,LastF,Increment
asel,s,loc,z,0
sfa,all,,pres,Force    allsel,all
nsubst,10,100,5
solve
*enddo
finish
!进入POST1后处理器/post1
set,last
pldisp,3
plnsol,u,x
pinsol,s,x
pinsol,s,y
pinsol,s,z
pinsol,s,xy
pinsol,s,eqvpinsol,epto,eqv
pinsol,cont,pres
pinsol,cont,sfric
finish                  !退出
/exit,all

zhuangdy

好东西，收了

guwu1020

谢啦，留着以后参考使用。

zhuangdy

收下了，谢谢楼主了

AvirusZ

毕设就是做这个，拿回去先看看，感激不尽

rhww

最近在搞毕业论文，或许用得到，谢谢楼主

a347524723

楼主好人啊

rhww

非常感谢楼主的无私奉献

zhangbojie2014

这个东西不错啊，学习下

1061800161

楼主太好了。。

yangheng8855

好东西，先留个言

jack916

感谢！学习了

satana123

zlm1966

好东西，谢谢分享

iamcxd

好东西，学习一下！谢谢了

liujianyuanhapp

学习
111111111111111111111