目录
0 引言
工程中经常遇到两个部件或者多个部件之间相互接触,当他们相互接触的时候,物体之间在法向存在相互作用力,如果接触面之间存在摩擦,比如放在斜坡中静止不动的滑块,沿着接触面的切线方向也会产生剪力用来抵抗物体的运动。通常接触分析是用来确定接触面的面积和接触应力。
在有限元分析中,接触是一种不连续的约束,他可以使得力从一个模型的一部分传递到另一部分,只有当两个表面接触的时候才附加接触条件,所以这种约束是不连续的。但是当两个接触的面分开的时候,这种约束也就不存在了。 所以在分析的过程中必须清楚的知道什么时候两个表面接触产生接触约束以及什么时候两个表面分开解除接触约束。
1 接触面直接的相互作用
接触面之间存在两种相互作用:一种是接触面之间的法向作用,另一种是接触面之间的切向作用。其中切向作用包括接触面之间的相对运动和存在的摩擦剪切应力。
两个接触面分开产生的距离叫做间隙。当接触变为0的时候,接触就起作用。当接触压力变为0或者负值的时候,接触面就分离,约束撤离,这个行为叫做“硬”接触。
除了以上要确定某一点是否发生接触之外,还需要计算两个表面之间的相互滑动关系。在ABAQUS中有小滑动量和有限滑动量,虽然也没有明确的界限什么叫做小滑动,具体数值是多少并没定义,但是通常遵循的原则是两个接触面(面——面接触)或一点和一面接触时,他们之间的滑动量不超过一个典型的单元尺度的很小部分,就可以认为是小滑动。
当两个表面接触的时候,接触面之间要传递切向力,所以在分析时需要考虑阻止相对滑动的摩擦力。通常我们用库伦摩擦来表示接触面的相互作用的模型,计算公式是,其中u表示摩擦因子,P表示正压力。当两个接触面之间没有发生相对滑动的时候,摩擦力小于库伦摩擦力,相对滑动量为0。
2 在ABAQUS中定义接触
第一步:创建两个构件之间的接触面。
第二步:创建相互接触面之间的相互作用。
第三步:每一个相互作用调用一个接触属性。
2.1 定义接触面
实体单元中的接触面
对于二维和三维实体单元,可以指定部件的区域形成接触面或由ABAQUS自动确定部件的自由面。对于前者可选择部件副本的面形成接触面,对于后者在定义接触面时只需简单地选择整个部件副本,ABAQUS 将略去实体内单元表,只保留与表面有关的单元。
壳、膜和刚性单元上的接触面
对于壳、膜和刚性单元,必须指明单元的哪个面来形成接触面。单元正法向方向的面称为SPOS,而单元负法向方向的面则称为SNEG。
刚性接触面
刚性接触面是刚性体的表面。刚性接触面可以定义为一个解析面或者基于刚性体的单元表面定义。
解析刚性接触面有三种基本形式。在二维模型中给出的解析刚性接触面是一个二维的分段刚性面。接触面的横截面轮廓线可在二维平面上用直线、圆弧,和曲线定义。三维的刚性接触面的横截面可用相同的方式在用户指定的平面上定义。这样这个横截面可以绕-个轴扫掠成回旋面或沿一个矢量拖拉成三维的面。
解析型的刚性接触面的优点在于只用少量的几个点便可定义,并且计算效率高。然而在三维情况下,创建的形状受到限制。
离散形式的刚性面是基于构成刚性体的单元的,这样它可以创建比解析刚性面更为复杂的刚性接触面。离散的刚性面创建的方法与可变形体的面的创建方法相同。
2.2 接触相互作用
默认设置是有限滑动公式。如果两个表面相对滑动量比单元面特征尺度小很多,使用小滑动公式计算效率更高。每一个相互作用必须调用接触属性。
从面和主面
ABAQUS在使用单纯的主从算法:从面上的节点不能侵入到主面的任何部分。该算法对主面没有限制;主面可以侵入从面。
通常主面、从面的选择规则为:
- 从面的网格应该划分的更加精细;
- 如果主面和从面的网格精度差不多,应该选择较为软的材料作为从面,较硬的材料作为主面。