abaqus第二讲：ABAQUS中实体单元课件.ppt

第二讲 ABAQUS中的实体单元;ABAQUS中的单元;ABAQUS中的单元;族 有限元族是一种广泛的分类方法。 同族的单元共享许多基本特征。 在同一族单元中又有许多变异。 ;ABAQUS中的单元;自由度 在有限元分析过程中，单元节点的自由度是基本变量。 自由度的例子： 位移 转动 温度 电势 ;ABAQUS中的单元;积分点 在单元之内，刚度和单元质量在采样点，所谓的“积分点”，进行数值计算。 用于积分这些变量的数值算法将影响单元的行为。 ABAQUS包含“全”积分和“减缩”积分单元。;所谓“完全积分”是指当单元具有规则形状时，所用的Gauss积分点的数目足以对单元刚度矩阵中的多项式进行精确积分。对六面体和四边形单元而言，所谓“规则形状”是指单元的边是直线并且边与边相交成直角，在任何边中的节点都位于边的中点上。完全积分的线性单元在每一个方向上采用两个积分点。因此，三维单元C3D8在单元中采用了2?2?2个积分点。完全积分的二次单元（仅存在于ABAQUS/Standard）在每一个方向上采用3个积分点。对于二维四边形单元，完全积分的积分点位置如图所示。;对于线性完全积分单元，在厚度方向的单元数目并不影响计算结果。误差是由于剪力自锁（shear locking）引起的，这是存在于所有完全积分、一阶实体单元中的问题。;为清楚起见，画出了通过积分点的虚线。显然，上部虚线的长度增加，说明1方向的应力是拉伸的。类似地，下部虚线的长度缩短，说明是压缩的。竖直方向虚线的长度没有改变（假设位移是很小的）；所有这些都与受纯弯曲的小块材料应力的预期状态是一致的。但是，在每一个积分点处，竖直线与水平线之间夹角开始时为90度，变形后却改变了，说明这些点上的剪应力不为零。显然，这是不正确的：在纯弯曲时，这一小块材料中的剪应力应该为零。 产生这种伪剪应力的原因是因为单元的边不能弯曲，它的出现意味着应变能正在产生剪切变形，而不是产生所希望的弯曲变形，因此总的挠度变小：即单元是过于的刚硬。 剪力自锁仅影响受弯曲载荷的完全积分的线性单元的行为。在受轴向或剪切载荷时，这些单元的功能表现很好。而二次单元的边界可以弯曲，故它没有剪力自锁的问题。二次单元预测的自由端位移接近于理论解答。;只有当确信载荷只会在模型中产生很小的弯曲时，才可以采用完全积分的线性单元。 如果对载荷产生的变形类型有所怀疑，则应采用不同类型的单元。 在复杂应力状态下，完全积分的二次单元也有可能发生自锁；因此，如果在模型中应用这类单元，应细心地检查计算结果。 然而，对于模拟局部应力集中的区域，应用这类单元是非常有用的！;;;;一次插值;单元命名约定：例子;比较ABAQUS/Standard和ABAQUS/Explicit单元库 两种程序基本上具有相同的单元族：连续体、壳、梁等等。 除了应力分析，ABAQUS/Standard包括许多可以用于其它分析类型的单元：热传导、土壤固结、声学等等。 在ABAQUS/Explicit中也可以使用声学单元。 对于每个单元族，ABAQUS/Standard包含许多变种。 ABAQUS/Explicit包含几乎所有的一阶单元。 例外：二阶三角形和四面体单元、二阶梁单元 对于两种程序，许多单元选择的准则是一样的。;网格生成;四边形/六面体 vs. 三角形/四面体单元 在生成网格时，选择使用四边形/六面体单元或三角形/四面体单元是非常重要的。 在可能的条件下，尽量使用四边形/六面体单元。 它们以最小的费用给出最好的结果。 在为复杂的几何体建模时，几乎没有任何的选择，必须使用三角形和四面体单元。;一阶三角形/四面体单元(CPE3, CPS3, CAX3, C3D4, C3D6)是质量较差的单元；它们有以下的问题： 较差的收敛率。 如果需要得到较好的结果，通常需要非常细的网格。 即使使用杂交公式，对于不可压材料或几乎不可压材料仍然会产生体积锁闭。 这些单元只能用在远离关心精度的区域，并把它们作为填料。;“常规的”二阶四面体、二阶楔形和六节点壳和薄膜单元(C3D10, C3D15, STRI65, M3D6)不能用于模拟接触问题，除非使用基于罚函数的接触公式。 与“经典”的硬接触相比，在单元角点和中点处，一致压力下面的接触力存在明显的不同。 ;修正的 二阶三角形/四面体单元(C3D10M, 等等)减轻了其它三角形/四面体单元的问题。 好的收敛率—与二阶四边形/六面体单元的收敛率相近。 最小化剪切锁闭和体积锁闭。 利用杂交公式(C3D10MH)，可以用于模拟不可压或几乎不可压材料。 在有限变形问题中，这些单元表现强劲。 一致的接触压力可以使这些单元精确的模拟接触问题。 使用它们!;网格细化和收敛性 使用充分细化的网格，以证明ABAQUS的模拟结果是让人满意的。 粗略的网格可