yebaobao0310 发表于 2011-6-8 15:31:02

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

桩土相互作用接触模型
/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/AIncrement=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 发表于 2012-2-21 21:04:03

好东西，收了

guwu1020 发表于 2012-3-8 20:51:49

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

zhuangdy 发表于 2012-3-19 14:41:56

收下了，谢谢楼主了

AvirusZ 发表于 2012-4-13 09:24:07

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

rhww 发表于 2012-4-14 20:51:17

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

a347524723 发表于 2014-1-13 15:01:05

楼主好人啊

rhww 发表于 2014-2-20 19:36:16

非常感谢楼主的无私奉献

zhangbojie2014 发表于 2014-3-4 10:35:09

这个东西不错啊，学习下

1061800161 发表于 2014-4-20 22:01:09

楼主太好了。。

yangheng8855 发表于 2015-3-29 00:47:20

好东西，先留个言

jack916 发表于 2015-4-16 10:21:33

感谢！学习了

satana123 发表于 2015-5-30 17:52:51

:):(:D:P:o

zlm1966 发表于 2015-7-15 08:43:46

好东西，谢谢分享

iamcxd 发表于 2017-8-10 23:19:40

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

liujianyuanhapp 发表于 2018-8-9 21:49:26

学习
111111111111111111111
页: [1]
查看完整版本: 分享一个桩土三维接触的实例