梁-实体MPC接触小例

这两天算是感受到MPC的好处了,用完了实体-实体的MPC又来搞梁-实体的,这两天先弄完手头上积压的工作,然后把各种图素之间的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

 2011-3-20 10-44-29

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

2011-3-20 10-47-04

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

长河

本名孙勇,目前在哈工大交通学院桥梁系读博士,正在毕业中,10kN.com博主。一家之言,仅供参考,若有不妥,欢迎探讨。

相关推荐

6 条评论

  1. J-yumin

    长河大哥 我想问你要是体单元包裹着梁单元该怎么用MPC处理啊~纠结了好久了

  2. 专业!

  3. 作为一个建筑学的学生,基本看不懂你发的结构模型。

  4. 现在学习中的站,请指点

  5. 这个是什么东西的模型吗

  6. 沙发,围观了

微信扫一扫

微信扫一扫

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

梁-实体MPC接触小例