梁-实体MPC接触小例

首先简单说一下用于梁-实体MPC接触的conta175单元。这是一个只有单一节点的单元,所以创建的时候只需要用E,而不是ESURF!!我这一个习惯性的ESURF废了我将近一个多小时,我说怎么用GUI都搞出来的东西,命令流怎么就总提示 No surface elements were generated呢,只有一个节点,自然没有surface elements啦,囧啊。

再有,虽然ansys帮助文档中,对于175的keyopt(3)和keyopt(4)有推荐值,但对于MPC算法而言,接触刚度什么的都是浮云,不用理会,只需按照如下小例中的方法设定即可。

finish
/clear
/prep7
et,1,solid95 !实体单元
mp,ex,1,3.5e4
mp,nuxy,1,0.2
mp,dens,1,2600

et,2,beam44 !梁单元
mp,ex,2,3.2e4
mp,dens,2,2500
mp,prxy,2,0.1667

et,3,175 !contact175 接触单元,准备做梁-实体的MPC
KEYOPT,3,2,2 !使用MPC多点约束算法
KEYOPT,3,12,5 !设置接触行为是;绑定粘结

et,4,170 !Target170

r,1 !空实常数,用于赋给实体
r,2 !空实常数,用于识别接触对

block,-3,0,0,2,-4,4 !创建实体
type,1,
mat,1
real,1
esize,0.5
mshape,0
mshkey,1
vmesh,all

*get,nodemax,node,0,num,max !创建梁
n,nodemax+1,0.1,1,0
ngen,6,1,nodemax+1,,,1
type,2
mat,2
sectype,1,beam,csolid,pier
secdata,0.3
secnum,1
*do,i,1,5
e,nodemax+i,nodemax+i+1
*enddo
allsel,all
/eshape,1
/VIEW,1,1,2,3
eplot

!开始创建MPC
nsel,s,loc,x,0
nsel,u,,,nodemax+1
nplot
type,4
real,2
esurf !生成实体接触点上的Target单元

nsel,s,,,nodemax+1
nplot
type,3
real,2
e,nodemax+1 !生成梁端部接触点上的contact单元

/solu
allsel,all
asel,s,loc,y,0
da,all,all
allsel,all
f,nodemax+6,fy,-2
solve

/post1
plnsol,s,1

梁端节点和实体可以有公共节点,也可以没有

结果显示,梁顺利地把荷载传递到实体上

Mark

这个人很懒,还没想好自我介绍

相关推荐

暂无评论

微信扫一扫,分享到朋友圈

梁-实体MPC接触小例