Q3T Prisms: A Linear-Quadratic Solid Shell Element for Elastoplastic Surfaces

一句话总结

提出一种沿厚度方向采用二次插值的三棱柱实体壳单元（Q3T），在面内保持线性、沿厚度保持二次的”线性-二次”混合运动学下，用统一的三维实体单元框架高效且准确地模拟可产生永久形变的弹塑性曲面。

研究背景

对可变形曲面（布料、纸张、金属薄板等）的力学建模是图形学几十年来的核心课题。这类物体的共同特点是一个维度（厚度）远小于另外两个维度，因此人们通常不直接用三维实体单元离散，而是退化到二维的薄壳（thin-shell）模型：只保留中面自由度，弯曲刚度由中面曲率的变化推导而来。

薄壳模型高效，但也有内在局限。它把厚度方向的力学”积分掉”，用中面上的曲率量近似弯曲响应。对于纯弹性、小厚度、以弯曲为主的场景这很好用；然而当曲面发生塑性形变——例如纸张被揉皱后留下的折痕、金属薄板被压出的凹坑——材料在厚度方向上呈现明显不均匀的应力/应变分布（内外侧受拉受压不同），这种沿厚度的梯度正是塑性折痕得以形成和保留的关键。纯中面的薄壳运动学难以自然表达这种厚度方向的分层行为。

另一条路线是直接用三维实体单元堆叠厚度，但为了刻画弯曲需要沿厚度放置多层单元，自由度和病态程度都急剧上升，且薄结构下实体单元容易出现各种”锁死”（locking）现象。本文的目标是取两者之长：保留实体单元能自然表达厚度方向应力梯度、直接套用三维本构（含塑性）的优点，同时避免多层堆叠的高昂代价。

方法

整体思路是设计一种单层三棱柱实体壳单元：面内用三角形的线性形函数，沿厚度方向引入二次插值（quadratic through-the-thickness，即 Q3T），使得单层单元就能表达厚度方向上非线性的应变分布，从而捕捉弯曲与塑性折痕所需的分层力学，而无需在厚度上堆叠多层网格。

flowchart LR
    A[三角网格曲面 + 厚度] --> B[三棱柱实体壳单元]
    B --> C[面内线性 / 厚度二次<br/>混合插值 Q3T]
    C --> D[三维变形梯度 F]
    D --> E[弹塑性乘法分解<br/>F = Fe · Fp]
    E --> F[弹性能 + 塑性流动/硬化]
    F --> G[增量隐式时间积分 / 能量极小化]

关键设计一：线性-二次混合运动学。单元在面内（沿三角形两条重心方向）采用线性位移场，在厚度方向采用二次位移场。二次项让单元内部沿厚度的应变可以非线性变化，这正是表达弯曲主导变形与塑性梯度所必需的自由度；而面内保持线性则控制住了单元的自由度规模与刚度矩阵稠密度。这种”低维方向用低阶、薄方向用高阶”的非对称设计，是单层单元同时兼顾精度与效率的核心。变形通过标准的变形梯度描述：

\[\boldsymbol{F} = \frac{\partial \boldsymbol{x}}{\partial \boldsymbol{X}},\]

其中 \(\boldsymbol{X}\) 为参考构型坐标、\(\boldsymbol{x}\) 为当前构型坐标；由于厚度方向插值为二次，\(\boldsymbol{F}\) 在单元内沿厚度是坐标的函数而非常量。

关键设计二：三维弹塑性本构。因为单元本质上是一个真实的三维实体，塑性可以直接套用连续介质力学中成熟的弹塑性框架，对变形梯度做乘法分解：

\[\boldsymbol{F} = \boldsymbol{F}_e\,\boldsymbol{F}_p,\]

其中 \(\boldsymbol{F}_e\) 为弹性部分、\(\boldsymbol{F}_p\) 为塑性部分。弹性能只由弹性部分决定：

\[\Psi = \Psi(\boldsymbol{F}_e) = \Psi(\boldsymbol{F}\,\boldsymbol{F}_p^{-1}),\]

塑性部分则按屈服准则与流动法则演化（当应力越过屈服面时发生塑性流动，并可带硬化）。相比薄壳需要专门推导”弯曲塑性/薄膜塑性”的降维本构，实体壳单元天然继承三维塑性模型，建模更直接、物理更透明。

关键设计三：抗锁死与单层离散。薄结构的实体单元通常受剪切锁死、体积锁死等困扰。Q3T 通过沿厚度的二次运动学缓解弯曲相关的锁死，使得单层三棱柱即可给出可接受的弯曲响应，避免了传统实体离散沿厚度堆叠多层带来的自由度膨胀。系统整体通过对总势能（弹性能 + 外力势 + 惯性项）做增量隐式积分/能量极小化求解，可与图形学中常见的隐式时间积分与优化式求解器对接。

实验结果

说明：本篇为在无法获取论文全文正文的情况下、基于摘要与公开信息整理的概要，未能提取到具体的定量实验数据，因此下表为方法定位层面的定性对照，不含论文中的具体数值。

模型	厚度方向运动学	塑性建模方式	折痕/凹陷等塑性梯度
传统薄壳 (thin-shell)	仅中面（无厚度自由度）	需专门推导降维塑性本构	难以自然表达
多层实体堆叠	多层线性	直接三维本构	可表达但自由度高、易锁死
本文 Q3T 棱柱	单层、厚度二次	直接三维弹塑性乘法分解	单层即可捕捉

论文的核心主张是：在与薄壳相近的离散规模下，Q3T 单层棱柱单元能够更忠实地再现弹塑性曲面在弯曲与永久变形（如揉皱、压凹、折痕保留）中的行为，同时规避多层实体堆叠的开销与锁死问题。

亮点与局限

亮点：

• “面内线性 + 厚度二次”的非对称插值设计，用单层单元就获得表达厚度方向应变梯度的能力，是效率与精度之间的巧妙折中。
• 作为真正的三维实体单元，可直接复用连续介质力学中成熟的弹塑性本构（乘法分解 + 屈服/流动/硬化），无需为壳专门推导降维塑性模型。
• 目标聚焦于弹塑性曲面这一薄壳模型的传统弱项，弥补了纯中面模型在塑性折痕、凹陷等厚度分层现象上的表达短板。

局限（部分为基于方法定位的推断）：

• 二次厚度运动学相比纯薄壳仍引入额外自由度与积分点，单元级成本高于最简薄壳。
• 单层棱柱对极端厚度比、复杂多层材料或强各向异性材料的适用范围有待更多验证。
• 本概要未能获取原文的具体基准、对比方法与定量指标，量化结论以论文原文为准。

延伸思考

这项工作体现了曲面仿真中一个反复出现的张力：把厚度”积分掉”换取效率（薄壳），还是保留厚度换取物理保真（实体）。Q3T 的答案是”只在需要的方向上保留高阶自由度”——面内维度用低阶、厚度这一关键方向用二次，从而在两条路线之间找到一个甜点。这种按物理各向异性来分配离散精度的思路，对其它降维仿真（如杆、梁、织物的多尺度建模）同样具有启发意义。

另一个值得关注的点是塑性建模的”接口”价值：因为单元回归到三维实体，它天然与连续介质塑性、损伤、断裂等成熟本构对接。这意味着后续可以较低成本地把更丰富的材料行为（率相关塑性、各向异性硬化、乃至与接触/自碰撞耦合的褶皱演化）叠加进来，使其成为模拟金属薄板成形、纸张揉皱、包装材料等工程与图形场景的统一底座。