《有限元法及其应用》，机械工业出版社，2006, 光盘演示算例

11 接触分析算例 (MARC)

(1) 首先定义上下接触表面，在marc main menu->mesh generation->add->points，加入两个关键点(0,0,0)和(1000,0,0)
(2) 进入marc main menu->mesh generation->add->curves，点击关键点1和2，生成直线
(3) 点击marc main menu->mesh generation->curve type，选择curve 类型为circles->center/radius.
(4) 进入marc main menu->mesh generation->add->curves，输入圆心坐标为(300,200,0)，半径为200

(5) 进入marc main menu->mesh generation->duplicate，设置复制的间距为(400,0,0)。点击刚才建立的圆环曲线（编号：2），复制。
(6) 设置复制的间距为(0,400,0)，点击首先生成的直线（编号1），复制。
(7) 下面开始进行单元网格划分。进入marc main menu->mesh generation->automesh->curve divisions，选择划分的方式为Fixed avg length，设定平均长度为30。点击apply curve divisions，选择all: Exist
(8) 进入marc main menu->mesh generation->automesh->2D planar meshing，选择tri mesh，选择刚才建立的两个圆环，完成网格划分。

(9) 下面开始设置材料属性。进入marc main menu->material properties->isotropic，设定弹性模量为200e3，泊松比为0.3，将材料赋予所有单元。点击elements->add all: exist
(10) 接着设置几何属性。进入marc main menu->geometry properties->planar->plane strain，设定厚度为1，将几何属性赋予所有单元。
(11) 下面开始定义接触。进入marc main menu->contact->contact body->deformable。设定摩擦系数为0.5。点击elements->add->all: exist

(12) 下面定义一个table来描述顶部刚体运动。点击tables，设定tables的名称为Load，表格自变量类型为time，加入两个点(0,0)和(1,1)，建立表格。
(13) 点击contact bodies->new，建立新的接触体为刚体(rigid)，输入摩擦系数为0.5，设定刚体运动方式为用位置描述position，设定位置变化为Y方向运动-100，运动方式遵循Load table定义的运动方式。点击curves->add，选择顶部的直线（编号：4）


(14) 点击contact bodies->new，设定摩擦系数为0.5，点击curves->add，点击底部直线（编号：1），点击flip curves，点击底部直线（编号：1）
(15) 进入marc main menu->loadcase，设定分析类型为static，点击solution control，将非正定分析non-positive conditions设置为on。


(16) 进入marc main menu->jobs->mechanical，设定分析的loadcase为lcase1。选择contact control，设定摩擦类型为coulumb，设定摩擦的参考速度为1e-3


(17) 最后选择element types->mechanical->plane strain solid->plane strain full integration->tria3，选择单元类型为6。将该单元类型添加给所有单元。
(18) 进入marc main menu->jobs->run，submit任务，直至出现3004结束标识。
(19) 进入results，选择open default，点击last，点击deformed shape->def only，设置scalar plot为contour band，设定scalar类型为equivalent stress，得到应力分布如图。