第三章有限元和刚性体.docVIP

第三章 有限元和刚性体 有限元和刚性体是ABAQUS模型的基本元素。有限元是可变形的。刚性体不可变形，只作刚体运动。用户对有限元多少有一些了解，而对有限元软件中刚性体概念可能感到有些陌生。 ABAQUS软件引入刚性体是为了提高计算效率。在ABAQUS中任何物体或物体的局部都可以定义为刚性体，大多数单元类型可为刚性体定义所用（参考ABAQUS/Stardard用户手册24－1节所列的表）刚性体与变形体结合的优点是对刚性体运动的完全描述可只用一个参考点，所以不超过六个自由度即可定位一个刚性体。与此相反，变形单元需要许多自由度，并要求大量的单元计算来确定变形。当某一局部变形可忽略或对它的变形不感兴趣时，则可在建模时构造一个刚性体部件，这样做就会极大地节省计算时间而并不影响整体结果。 3.1 有限元 ABAQUS有各种各样的单元，其庞大的单元库提供了一套强大的工具来解决许多不同类型的问题，本节介绍影响单元特性的五个方面问题。 3.1.1 单元的表征 每一个单元都由下面几个特性来表征： ·单元族 ·自由度（和单元族直接相关） ·节点数 ·数学描述，即单元列式 ·积分 ABAQUS中每一种单元都有自己特有的名字，例如T2D2，S4R和C3D81。单元的名字标志着一种单元的五个特性。单元命名的规则将在本章里说明。 单元族 图3－1给出了应力分析中最常用的单元族。单元族之间一个明显的区别是每一个单元族所假定的几何类型不同。 图3－1 常用单元族 在这本指南里将用到的单元族有实体单元、壳单元、梁单元、桁架和刚性体单元，这些将在以后的各章里详细讨论。其它的单元族在这本指南中没有讲到；如果对应用它们感兴趣，请查阅ABAQUS/Standard用户手册的Part V。 单元名字里开始的字母标志着这种单元属于哪一个单元族。例如，S4R中的S表示它是壳单元，C3D81中的C表示它是实体单元。 自由度 自由度（dof）是分析中计算的基本变量。对于壳和梁单元的应力/位移模拟分析，自由度是每一节点处的平动和转动。对于热传导模拟分析，自由度为每一节点处的温度；因此，热传导分析要求应用与应力分析不同的单元，因为它们的自由度不同。 ABAQUS中自由度的排序规则如下： 1 1方向的平动 2 2方向的平动 3 3方向的平动 4 绕1轴的转动 5 绕2轴的转动 6 绕3轴的转动 7 开口截面梁单元的翘曲 8 声压或孔隙压力 9 电势 11 温度（或物质扩散分析中归一化浓度），对梁和壳，指厚度方向第一点温度。 12 梁和壳厚度上其它点的温度 方向1，2，3分别对应于整体坐标的1－，2－，和3－方向，除非已经在节点处定义了局部坐标系。 轴对称单元是一个例外，其位移和转动自由度指的是： 1 r-方向的平动 2 z-方向的平动 3 r-z平面内的转动 方向r和z分别对应于整体坐标的1-和2-方向，除非已经在节点处定义了局部坐标系。如何在节点处定义局部坐标系，请参考第5章的讨论。 在本指南中，限于结构方面的应用，所以只讨论单元的平动和转动自由度。其它类型的单元（如热传导单元）可参考ABAQUS/Standard用户手册。 节点数——插值的阶数 ABAQUS仅在单元的节点处计算位移或任何其它的自由度。在单元内的任何其它点处，位移是节点位移插值获得的。通常插值的阶数由单元采用的节点数决定。仅在角点处在节点的单元，例如图3－2（a）所示的8节点实体单元，在每一方向上采用线性插值，因此常常称这类单元为线性单元或一阶单元。具有边中点节点的单元，如图3－2（b）所示的20节点实体单元，采用二次插值，因此常常被称为二次单元或二阶单元。一般情况下单元的节点数在其名字中清楚地标记着。8节点实体单元，如你前面已经看到过的那样，叫做C3D8；8节点一般壳单元叫做S8R。梁单元族的记法稍有不同：插值的阶数在单元的名字中标记着。这样，一阶三维梁单元叫做B31，而二阶三维梁单元叫做B32。 图3－2线性实体单元和二次实体单元 单元列式 单元列式是指用来定义单元行为的数学理论。ABAQUS中所有的应力/位移单元行为都是基于拉格朗日或物质描述的：在整个分析过程中和一个单元相关的物质保持和这个单元相关，而且物质不能穿越单元边界。在欧拉或空间描述中，单元在空间固定，而物质在单元之间流动。欧拉方法通常用于流体力学分析。ABAQUS运用欧拉方法来模拟对流换热，这本指南中不讨论这一题目。 为了适用于不同类型的物理行为，ABAQUS中的某些单元族包含具有几种不同列式的单元。例如，壳单元族有三个类别：一类具有一般壳体理论的列式，一类具有薄壳理论的列式，另外一类具有厚壳理论的列式（这些单元的列式在第5章解释）。 某些单元族除了有标准的