学术报告：《International Journal of Plasticity》2007年关于RESS单元应用于板料成形仿真的研究

一、 主要作者、研究机构及发表情况

本研究的主要作者是来自葡萄牙阿威罗大学（University of Aveiro）机械工程系的R.J. Alves de Sousa、R.P.R. Cardoso、R.A. Fontes Valente、J.J. Grácio，以及来自美国美铝技术中心（Alcoa Technical Center）材料科学部的J.W. Yoon。该研究成果于2007年发表在著名固体力学与材料本构行为期刊《International Journal of Plasticity》上。

二、 学术背景与研究目的

该研究属于计算固体力学和板料成形工艺数值仿真领域。其核心背景在于解决大规模工业级板料成形（如冲压、液压成形）数值模拟中长期存在的效率与精度平衡问题。

传统上，学者们主要采用两类单元进行薄壁结构模拟：传统的三维实体（3D solid）单元和壳（shell）单元。传统实体单元（如八节点六面体单元）虽可直接应用三维本构模型、方便处理双面接触问题，但其在模拟弯曲主导的薄壁结构时，极易产生多种“闭锁”（Locking）现象（如体积闭锁、横向剪切闭锁、厚度闭锁），导致结果刚度过大、严重失真。为避免闭锁，常需采用完全积分（2x2x2个高斯点），但计算成本高昂。而壳单元采用平面应力假设，专门针对薄壁结构优化，计算效率高，但在处理三维应力状态、双面接触以及厚度方向应变变化时存在天然的理论简化限制。

为了结合二者的优点，“实体-壳”（Solid-Shell）单元概念应运而生。这类单元保持了实体单元的八节点拓扑结构和三维属性，但在数学上进行了特殊处理以提高弯曲性能。然而，开发一种既能避免各种闭锁、又能高效计算、还能精确描述厚度方向应力应变（对回弹预测至关重要）的实体-壳单元，仍是一个挑战。特别是，要准确模拟板料成形中的回弹（Springback），通常需要在厚度方向布置多个积分点，传统做法是堆叠多层单元，但这显著增加了计算量。

因此，本研究的核心目的是评估并展示一种新近提出的“缩减增强实体-壳”（Reduced Enhanced Solid-Shell， RESS）单元在大型弹塑性板料成形模拟中的性能。具体目标包括：验证该单元在接触、各向异性塑性以及回弹效应等复杂工况下的计算精度与效率；展示其单层单元内即可沿厚度方向设置任意数量积分点的独特能力，从而在保证精度的同时显著提升计算效率。

三、 详细工作流程

本研究并非传统意义上的实验研究，而是一项计算方法和数值验证研究。其工作流程围绕RESS单元的理论阐述、算法实现及其在一系列标准板料成形基准案例中的应用验证展开。核心流程如下：

1. RESS单元理论与算法开发（基于前期工作） 本研究的基础是作者团队在2005和2006年发表于《International Journal for Numerical Methods in Engineering》的两篇论文中提出的RESS单元。该单元的核心创新在于巧妙结合了两种现有技术： * 增强假定应变法（Enhanced Assumed Strain, EAS）：通过引入内部增强变量来扩展应变场，以克服闭锁问题。RESS单元的独特之处在于，它仅采用一个增强参数（通常EAS单元需要10个甚至更多），即可有效解决体积闭锁和厚度闭锁，这大大简化了计算。 * 缩减积分（Reduced Integration, RI）方案：在单元面内（即ξ-η平面）仅使用一个高斯积分点（位于单元中心），以此消除横向剪切闭锁，并大幅提高计算效率。

最关键的特色在于，RESS单元的积分方案是“半缩减”的：面内一点积分，但沿厚度方向（ζ坐标）可以设置任意数量的高斯积分点。这意味着仅需单层网格，就能以可控的成本精确捕捉厚度方向的应力梯度，这对于预测回弹至关重要。 为了应对面内一点积分带来的“沙漏”（Hourglass）或“零能”模态（即虚假的变形模式）问题，研究采用了Cardoso等人（2002）发展的物理稳定化程序。该方法基于对单元应变位移算子进行分解，推导出抵抗沙漏模式的附加刚度和力向量，无需经验参数，计算成本低。 本构模型方面，研究集成了适用于大变形弹塑性的通用模型，特别强调了塑性各向异性的处理。研究采用了先进的各向异性屈服函数，如Barlat等人（1991）提出的Yld91和（2005）提出的Yld2004-18p，这些模型能够更精确地描述铝合金等材料的微观结构导致的宏观各向异性行为。应力更新算法采用了基于增量变形理论的多阶段返回映射方法，能够处理非二次屈服函数和大应变增量。

2. 数值验证与性能评估 研究通过四个经典的板料成形模拟案例，系统地评估了RESS单元的性能。所有模拟均使用作者团队开发的隐式求解有限元软件CEREBRO进行。 * 案例一：Al2008-T4铝合金液压胀形试验 * 研究目的：验证单元在双轴应力状态下预测厚度分布和各向异性材料响应的能力。 * 建模：利用对称性，分析四分之一模型。采用Yld91屈服函数描述材料各向异性，沿厚度方向使用7个积分点。 * 对比：将模拟得到的极点位移 与液压压力的关系曲线、以及厚度分布，与实验结果以及文献中基于连续体单元和壳单元的结果进行对比。同时，与各向同性的Von Mises模型结果对比，以凸显考虑各向异性的重要性。 * 案例二：AA2090-T3铝合金圆柱杯拉深及凸耳预测 * 研究目的：测试单元在涉及双面接触的复杂成形工艺中，预测由材料各向异性引起的“凸耳”（Earing）现象的能力。 * 建模：由于恒定压边力，属于典型双面接触问题。利用正交各向异性，分析四分之一模型，共1125个RESS单元，沿厚度5个积分点。 * 对比：分别使用Yld91和Yld2004-18p屈服函数进行模拟，将得到的杯体高度轮廓（即凸耳轮廓）与实验数据对比。重点观察不同屈服函数对凸耳数量（4个还是6个）的预测准确性。 * 案例三：方形模具管材液压成形 * 研究目的：考察单元在预测成形压力历程和最终厚度分布方面的准确性。 * 建模：分析四分之一模型。使用40个RESS单元，但为了精确描述塑性行为和厚度变化，在厚度方向设置了多达15个积分点。材料采用各向同性的Mises屈服准则和幂律硬化模型。 * 对比：将模拟得到的管材圆角半径随内压变化的关系，以及成形后的厚度分布，与文献（Hwang and Chen， 2005）中的实验数据及有限元结果进行对比。 * 案例四：无约束圆柱弯曲与回弹 * 研究目的：验证单元在弯曲主导且回弹严重的工况下，仅用单层单元即可准确预测回弹量的能力。 * 建模：模型几何参考Numisheet 2002基准题。使用RESS单元，沿厚度7个积分点。 * 对比：将模拟得到的回弹前后的参考角度与实验结果（Yoon et al., 2002）进行直接对比。

四、 主要结果

研究在所有四个验证案例中都获得了与实验数据或公认的参考文献高度一致的结果，充分证明了RESS单元的精度和鲁棒性。

1. 在液压胀形试验中，使用Yld91屈服函数的RESS单元模拟结果，其压力 -极点位移曲线与实验数据吻合极佳。相比之下，采用Von Mises模型的模拟结果则与实验存在明显偏差，这证明了在板料成形模拟中考虑材料平面各向异性的必要性。同时，RESS单元的结果与早期使用传统实体单元和壳单元的高质量模拟结果一致，表明该单元在捕获双轴应力状态和厚度应变方面是准确的。

2. 在圆柱杯拉深案例中，RESS单元成功模拟了双面接触过程。更重要的是，凸耳预测结果清晰地展示了不同屈服函数的影响：使用Yld91屈服函数预测出4个凸耳，而使用更先进的Yld2004-18p屈服函数则成功预测出与实验相符的6个凸耳。这一结果与文献中使用壳单元（Yld91）和实体单元（Yld2004-18p）的结论一致，表明RESS单元能够无缝集成复杂的三维各向异性本构模型，并有效应用于实际成形问题。

3. 在管材液压成形案例中，RESS单元的模拟结果（无摩擦假设）与参考文献中有限元模拟结果（低摩擦系数）高度吻合。模拟准确预测了随着内压增加，管材圆角半径逐渐减小并趋于稳定的趋势。在厚度分布上，RESS单元与参考文献的有限元结果均显示最小厚度出现在接触区与自由膨胀区的过渡部位（约25°处），而实验数据显示最小厚度在拐角处（30°）。这种差异可能源于实际实验中的摩擦条件与模拟假设的差异。总体而言，该案例验证了RESS单元处理复杂几何和载荷历史的有效性。

4. 在回弹预测案例中，RESS单元取得了令人印象深刻的结果。模拟得到的回弹前和回弹后的参考角度分别为37°和67°，与实验测量的35°和68°非常接近。这一结果的关键意义在于，它仅通过单层RESS单元（7个厚度方向积分点）就实现了高精度的回弹预测，避免了传统实体单元需要堆叠多层网格带来的计算负担。

五、 研究结论与价值

本研究系统评估了新型缩减增强实体-壳（RESS）单元在大变形弹塑性板料成形模拟中的应用性能。研究得出以下核心结论： 1. 成功验证：RESS单元在几何非线性、材料非线性（包括复杂塑性各向异性）、接触以及回弹等一系列板料成形关键问题上，均展现出与实验结果和既有可靠模拟方法相媲美的精度和鲁棒性。 2. 独特优势得以体现：该单元最具特色的能力——在单层单元内沿厚度方向设置任意数量积分点——在模拟中得到充分发挥。这使得在需要高精度厚度方向分辨率（如回弹分析）时，无需采用计算代价高昂的多层网格，从而在保证精度的前提下显著提升了计算效率。 3. 综合性能优越：RESS单元有效结合了实体单元和壳单元的优点：它像实体单元一样直接使用三维本构、方便处理双面接触、自然计算厚度应变；同时，它通过引入最小数量（仅一个）的EAS增强参数和面内一点积分方案，高效克服了低阶实体单元的各类闭锁问题，并采用了物理稳定化方法控制沙漏模式，获得了与壳单元相当的弯曲计算效率。 4. 应用前景明确：研究表明，RESS单元为工业级板料成形过程的数值模拟提供了一种高效且准确的解决方案，特别适合于涉及严重弯曲、回弹、各向异性材料以及双面接触的复杂工况。

六、 研究亮点

1. 方法学的创新性结合：将“单参数EAS增强”与“面内缩减积分+厚度方向多积分点”方案相结合，是RESS单元设计的核心创新。这种设计以最小的额外计算开销（仅一个内部变量），系统性地解决了闭锁问题，并赋予了单元在积分点布置上的灵活性。
2. 显著的计算效率优势：相比需要多层网格的传统实体单元方法，以及相比难以处理双面接触和全三维应力的传统壳单元，RESS单元在目标应用场景中找到了一个更优的平衡点，为实现高精度、高效率的工业成形仿真提供了新工具。
3. 先进的材料模型集成：研究没有停留在简单的各向同性塑性模型，而是成功集成了当时最先进的各向异性屈服准则（Yld91， Yld2004-18p），并验证了其在预测实际材料行为（如凸耳）中的关键作用，提升了研究的工程实用价值。
4. 全面的验证体系：通过从简单胀形到复杂拉深、从压力预测到回弹计算的一系列经典基准案例，多角度、全方位地验证了单元性能，论证充分有力。

七、 其他有价值内容

本研究还详细阐述了RESS单元理论框架中的一些关键技术细节，例如： * 协同旋转变换：采用基于极分解的协同旋转坐标系来满足大变形下的客观性要求，并使非线性算法的推导更为直接。 * 物理稳定化程序的细节：附录中给出了应变-位移算子分解的完整数学表达式，为其他研究者理解和实现该方法提供了便利。 * 本构积分算法：简要介绍了处理大变形和非二次屈服函数的“增量变形理论”和“多阶段返回映射”策略，强调了其对较大应变增量的适应性。

这篇论文不仅是一份详尽的数值验证报告，更是一份关于高效高精度实体-壳单元设计理念与实现技术的出色阐述，对计算力学和先进制造领域的科研人员与工程师具有重要的参考价值。