ABAQUS中实体单元的应用

ABAQUS中实体单元的应用（完整可以看附件）

（此为学习时看到的一个很好的讲解。                  学习交流群：1063594113）

在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－1列出了不同网格情况下自由端位移与梁的理论解3.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方向的应力，s₁₁，是拉伸的。类似地，下部直线的长度缩短，说明s₁₁是压缩的。竖直直线的长度没有改变(假设位移很小)。因此，所有积分点上的s₂₂为零。所有这些结论与受纯弯的小块材料所预计的应力状态是一致的。但是在每一个积分点，竖直线与水平线之间夹角开始时是，变形后改变了。这说明每一点的剪应力 s₁₂不为零。这是不正确的：纯弯时一小块材料中的剪应力应为零。

图4－4  受弯曲材料的变形

 

图4－5  受弯曲的完全积分线性单元的变形

 

 

 

出现这个伪剪应力的原因是因为单元的边不能弯曲。它的存在意味着应变能导致剪切变形，而不是导致弯曲变形，其结果导致总的挠度变小了：即单元太刚硬了。剪力锁闭只影响受弯曲载荷的完全积分线性单元，这些单元的功能在受纵向或剪切荷载时并没有问题。而二次单元的边界可以弯曲(见图4－6)，故它没有剪力锁闭的问题。对表4－1所示的二次单元，计算所得的自由端位移接近于理论解。但是，如果二次单元扭曲或弯曲应力有梯度，则也可能出现某些锁闭现象，而这两种情况在实际问题中是可能发生的。

只有在确认载荷将产生小弯曲时，才可采用完全积分的线性单元。而如果对载荷产生的位移类型有怀疑，则应采用不同的单元类型。在复杂应力状态下，完全积分的二次单元也可能发生锁闭。因此如果在模型中有此类单元，则应细心地检查计算的结果。但是，对于局部应力集中问题，完全积分的线性单元是非常有用的。

 

图4－6 受弯曲的完全积分二次单元的变形

 

4.1.2  减缩积分

 

只有四边形和六面体单元才能采用减缩积分；而所有的楔形体、四面体和三角形实体单元只能采用完全积分，即使它们与减缩积分的六面体或四边形单元用在同一个网格中。

减缩积分单元比完全积分单元在每个方向少用一个积分点。减缩积分的线性单元只在单元中心有一个积分点。（实际上，在ABAQUS中这些一阶单元采用了更精确的均匀应变公式，对此单元计算了其应变分量的平均值。在这里的讨论中此种区别是不重要的）。对减缩积分四边形单元，积分点的位置如图4－7所示：

 

图4－7  采用减缩积分的二维单元的积分点

 

 

利用前叙的四类单元及图4－3所示的四种有限元网格，通过减缩积分来对悬臂梁问题进行计算，其结果列于表4－2。

 

表4－2  减缩积分单元的梁挠度比值

单元	网格尺寸(高度´长度)
1´6	2´12	4´12	8´24
CPS4R	20.3*	1.308	1.051	1.012
CPS8R	1.000	1.000	1.000	1.000
C3D8R	70.1*	1.323	1.063	1.015
C3D20R	1.000	1.000	1.000	1.000

^* 没有刚度抵抗所加载荷

 

线性的减缩积分单元由于存在着所谓沙漏 (hourglassing) 的数值问题而过于柔软。再一次考虑用单个减缩单元模拟受纯弯载荷的小块材料（见图4－8）。

 

图4－8  受弯曲的减缩积分线性单元的位移

 

单元中虚线的长度均没有改变，并且它们的夹角也没有改变，这意味着在单元单个积分点上的所有应力分量都为零。由于单元变形没有产生应变能，所以这种弯曲的变形模式是一个零能量模式。由于单元在此模式下没有刚度，所以不能抵抗此种形式的位移。在粗网格中，这种零能量模式会通过网格扩展出去，从而产生无意义的结果，这就是所谓的沙漏问题。

可在ABAQUS中对减缩积分单元引入少量的人工“沙漏刚度”以限制沙漏模式的扩展。当模型中有更多的单元时，这种刚度在限制沙漏模式方面是更有效的，这意味着只要采用合理的细网格，线性减缩积分单元会给出可接受的结果。对许多应用而言，采用细网格的线性减缩积分单元所产生的误差是在一个可接受的范围内的。这个结果说明当用这类单元来模拟承受弯曲载荷的结构时，在厚度方向上至少应采用四个单元。当在梁的厚度方向只有一个线性减缩积分单元时，所有的积分点都位于中性轴上，从而该模型将不能抵抗弯曲载荷。(这种情况在表4－2中用*标出)。

因为线性减缩积分单元对变形的鲁棒性，因此可在变形很大的模拟中采用剖分较细的此类单元。

二次减缩积分单元也有沙漏模式。然而在正常网格中这种模式几乎不可能扩展出去，并且在网格足够细时基本上不会造成什么问题。由于沙漏问题,C3D20R单元的1´6网格计算发散；若在宽度方向上变为两个单元，即2×6网格，就不会发散，但对于更细的网格，即便在宽度方向上只有一个单元也不会发散。即使在复杂应状态下，二次减缩积分单元对锁闭并不敏感。因此一般来说，除了大应变的大位移问题和一些接触分析问题外，这些单元是应力/位移模拟最佳选择。

 

 

 

 

4.1.3  非协调单元

 

非协调单元是克服完全积分的一阶单元的剪力锁闭问题的一种尝试。既然剪力锁闭是由于单元的位移场不能模拟与弯曲相关的运动学而引起的，那么可以考虑把增强单元变形梯度的附加自由度引入到一阶单元中去。对变形梯度的加强使一阶单元在单元中的变形梯度呈线性变化，如图4－9(a)所示。在标准单元列式中，变形梯度在单元中是常量，见图4－9(b)所示，故标准单元列式必然导致与剪力锁闭相关的非零剪切应力。变形梯度的增强完全是在单元内部的，并且与边节点无关。与直接增强位移场的非协调模式的单元列式不同，在ABAQUS中所采用的列式不会导致图4－10那样的两个单元交界处的重叠或裂隙，进而ABAQUS中的非协调单元列式很容易拓广到非线性有限应变模拟以及某些难以采用增强位移场的场合。

图4－9  位移梯度的变化 (a) 非协调单元（增强位移梯度）和 (b) 采用标准构造的一阶单元

 

 

图4－10 利用增强位移场而不是增强位移梯度所导致的非协调单元的可能运动非协调性。ABAQUS对非协调单元采用了增强位移梯度形式

 

在弯曲问题中，非协调元可得到与二次单元相当的结果，而计算费用却明显降低。但非协调元对单元扭曲很敏感。图4－11表示用有意扭歪的非协调单元来模拟悬臂梁：一种情况是“平行”扭歪，另一种是“交错”扭歪。

图4－12画出了悬臂梁模型的自由端位移相对于单元扭歪水平的曲线。图中比较了三类平面应力单元：完全积分的线性单元、减缩积分的二次单元以及线性非协调单元。象所预见的那样，完全积分的线性单元的结果较差。而减缩积分的二次单元则给出了很好的结果，直到单元扭歪得很严重时其结果才会恶化。

当非协调单元是矩形时，即使在悬臂的厚度方向只有一个单元，也能给出与理论值十分相近的结果。但是即使很小的交错扭歪也使单元过于刚硬。平行扭歪也降低了单元的精度，但程度较小。

 

图 4－11  非协调单元的扭歪网格

 

 

图4－12  平行和交错扭曲对非协调单元的影响

 

非协调单元之所以有用，是因为如果应用得当，则在很低花费时仍可得到较高的精度。但是必须注意保证单元扭歪是非常小的，然而当网格较复杂时这一点是很难保证的；因此，对于具有这种几何形状的模型，应再次考虑应减缩积分的二次单元，因为它们对网格扭歪并不敏感。

 

4.1.4  杂交单元

 

ABAQUS中的每一种实体单元，包括所有的减缩积分单元和非协调单元，都还有杂交单元列式。杂交单元名字前标有字母“H”。

对不可压缩材料（泊松比＝0.5）或非常接近于不可压缩的材料（泊松比>0.495）问题需采用杂交单元。橡胶就是具有不可压缩性质的材料的例子。不能用常规单元来模拟不可压缩材料的响应（除了平面应力情况），这是因为在单元中的压应力是不确定的。现在考虑均匀静水压力作用下的一个        图4－13  在静水压力下的单元单元（图4－13）。                                   

如果材料不可压缩，其体积在载荷作用下并不改变。因此压应力不能由节点位移计算，对于具有不可压缩材料性质的单元，一个纯位移列式是不适定的。

杂交单元包含一个可直接确定单元压应力的附加自由度。其节点位移只用来计算偏（剪）应变和偏应力。

在第8章将给出对橡胶材料的更详细的描述。

 

4.2  选择实体单元

 

对某一具体的模拟计算，如果要想以合理的费用达到精确的结果，则正确地选择单元是非常关键的。在使用ABAQUS的经验日益丰富时，毫无疑问每个用户会建立起自己的单元选择准则来解决具体问题，但若是刚开始使用ABAQUS，则可考虑下面的建议：

l   如果不需要模拟非常大的应变或进行复杂的需改变接触条件的问题，则应采用二次减缩积分单元（CAX8R，CPE8R，CPS8R，C3D20R等）。

l   如果存在应力集中，则应在局部采用二次完全积分单元（CAX8，CPE8，CPS8，C3D20等）。它们可用最低费用提供应力梯度最好的解答。

l   涉及到有非常大的网格扭曲问题（大应变分析），建议采用细网格剖分的线性减缩积分单元（CAX4R，CPE4R，CPS4R，C3D8R等）。

l   对接触问题采用线性减缩积分单元或细分的非协调单元（CAX4I，CPE4I，CPS4II，C3D8I等）。详见第11章。

l   尽可能地减少网格形状的扭歪，形状扭歪的粗网格线性单元会导致非常差的结果。

l   对三维问题应尽可能采用六面体单元。它们以最小费用给出最好的结果。当几何形状复杂时，完全采用六面体单元构造网格往往难以办到；因此可能需要采用楔形和四面体单元。众所周知，这些形状的一阶单元，如C3D6和C3D4，是较差的单元；若要取得较好的精度，需剖分很细的网格，因此，只有在为了完成网格建模而万不得已的情况下才会应用这些单元，即使如此，这些单元也应远离精度要求较高的区域。

l   一些前处理程序包含了自由网格算法，它们可用四面体单元构造任意形状的网格。只要采用二次四面体单元（C3D10），除了接触问题，其结果对小位移问题应该是合理的。C3D10单元的修正单元C3D10M对大变形问题、接触问题有鲁棒性，并表现出最小剪切和体积锁闭性质。但无论采用何种四面体单元，计算所花费的时间都多于采用相应密度的六面体单元。建议不采用只包含线性四面体单元（C3D4）的网格，因为如果不用大量的单元其结果将是不准确的。

附件：里面包括以上所述以及一个实例分析讲解

ABAQUS中实体单元的应用.doc