ABAQUS中实体单元的应用.doc

ABAQUS中实体单元的应用 在ABAQUS的单元库中，应用最广泛的是应力/位移实体单元族。对三维单元，可以选择六面体、四面体和楔形体；对二维单元则可在三角形与四边形之间进行选择。这些基本的单元形状，每一种都有线性和二次的两类选择。对六面体和四边形，还可选择完全积分或减缩积分。最后，还可选用标准元或杂交元列式。另外对线性六面体或四边形单元，还有个附加的功能，可选择非协调模式，而对二次的三角形或四面体单元可以应用修正列式。 若列出所有种类的单元，所面临的实体单元的总数目是相当大的，仅三维单元而言就超过20种。模拟的精度将强烈地依赖于所采用的单元类型。特别是在初次使用时，在这些单元中选择哪一个最为合适很可能是一件令人苦恼的事情。然而，用户会逐渐把这个工作看作是从一个20多件的工具组中，有能力选择最恰当的工具或单元来完成的一个有价值的工作。 这一章讨论了不同的单元列式和积分水平对一个特定分析的精度的影响。同时也讨论了一些选择实体单元的一般性原则。这些讨论提供了获得更多应用ABAQUS经验和知识的基础。在本节末的例子将允许用户应用这些知识建立和分析一个连接柄构件的模型。 4.1 单元列式和积分 通过图4-1所示的悬臂梁，可阐明单元阶数（线性或二次），单元列式及积分水平等因素对结构模拟精度的影响。这是评估一个给定单元的性能的经典测试。因为该构件相对是细长的，我们通常用梁单元来对它建立模型。但在这里我们用这个测试来帮助评估各种实体单元的效率。 梁长150mm，宽2.5mm，高5mm；一端固定；自由端承受5N的荷载。材料的杨氏模量E为70GPa，泊松比为0.0。采用梁的理论，在载荷P作用下，梁自由端的挠度为 其中，是长度，b是宽度，d是梁的高度。 P = 5N时自由端挠度是3.09mm。 图4-1 自由端受集中载荷的悬臂梁 4.1.1 完全积分 所谓“完全积分”是指当单元具有规则形状时，所用的Gauss积分点的数目足以对单元刚度矩阵中的多项式进行精确积分。对六面体和四边形单元而言，所谓“规则形状”是指单元的边相交成直角，而任何的节点位于边的中点。线性单元如要完全积分，则在每一方向需要两个积分点。因此，三维单元C3D8在单元中排列了2(2(2个积分点。而二次单元如要完全积分则在每一方向需要3个积分点。在完全积分的二维四边形单元中积分点的位置如图4-2所示。 图4－2 如图4－3(一阶与二阶)和网格密度对结果精度的影响。 表4－13.09mm的比值。 用线性单元CPS4和C3D8所得的挠度值是如此之差以至于其结果是不可用的。网格越粗，结果的精度越差，但即使网格划分得相当细(8(24)，得到的位移仍只是理论值的56%。注意到对线性完全积分单元而言，在厚度方向单元的剖分数并不会造成什么差异。这是由剪力锁闭引起的，它是对所有完全积分的一阶实体单元都存在的问题。 图4－3 表4－1 单元 网格尺寸(高度(长度) 1(6 2(12 4(12 8(24 CPS4 0.074 0.242 0.242 0.561 CPS8 0.994 1.000 1.000 1.000 C3D8 0.077 0.248 0.243 0.563 C3D20 0.994 1.000 1.000 1.000 正如我们已经看到的，剪力锁闭使单元在弯曲时过于刚硬。对之可作如下解释：考虑一个受纯弯的结构中的一小块材料，材料将产生的弯曲如图4－4。 因为线性单元的边不能弯曲，所以，如果用单个单元来模拟小块材料，则其变形后的形状如图4－5 为清楚起见，画出了通过积分点的虚线。很明显，上部直线的长度增加，这说明1方向的应力，(11，是拉伸的。类似地，下部直线的长度缩短，说明(11是压缩的。竖直直线的长度没有改变(假设位移很小)。因此，所有积分点上的(22为零。所有这些结论与受纯弯的小块材料所预计的应力状态是一致的。但是在每一个积分点，竖直线与水平线之间夹角开始时是，变形后改变了。这说明每一点的剪应力 (12不为零。这是不正确的：纯弯时一小块材料中的剪应力应为零。 图4－4 图4－5 (见图4－6)4－1 只有在确认载荷将产生小弯曲时，才可采用完全积分的线性单元。而如果对载荷产生的位移类型有怀疑，则应采用不同的单元类型。在复杂应力状态下，完全积分的二次单元也可能发生锁闭。因此如果在模型中有此类单元，则应细心地检查计算的结果。但是，对于局部应力集中问题，完全积分的线性单元是非常有用的。 图4－6 4.1.2 减缩积分 只有四边形和六面体单元才能采用减缩积分；而所有的楔形体、四面体和三角形实体单元只能采用完全积分，即使它们与减缩积分的六面体或四边形单元用在同一个网格中。 减缩积分单元比完全积分单元在每个方向少用一个积分点。减缩积分的线性单元只在单元中心