梁-实体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

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

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