A. hexahedron element formulation with a new multi-resolution analysis

本文提出了一种基于相互嵌套位移子空间新颖框架的新型多分辨率六面体单元公式化方法，该方法允许通过可调节的分辨率层级，比传统的单分辨率方法更高效、更合理地调节结构分析精度。

要点

想象一下，你正试图拍摄一张复杂的 3D 物体（例如一座桥梁或一栋建筑）的照片，以观察它在压力下是如何弯曲或振动的。在工程领域，这种方法被称为有限元法 (Finite Element Method, FEM)。

旧方法：像素化的拼图

传统上，工程师将结构分解为由许多小 3D 块组成的网格（就像一个个小骰子）。这些被称为“六面体单元 (hexahedron elements)”。

• 问题所在： 每个方块都是“单分辨率”的，这意味着它拥有固定数量的角点（节点），即进行计算的点。通常只有 8 个角点。
• 局限性： 如果你想观察材料上的微小裂纹或剧烈弯曲，旧方法会迫使你将整个结构切割成更多、更多更小的方块。这就像是为了让照片更清晰，而不得不把整张图片切成数百万个细小的、独立的拼图碎片，然后再把它们粘回去。这种做法既慢又乱，而且每次想要更多细节时都需要重新进行繁琐的操作（重划网格/re-meshing）。

新方法：“智能变焦”镜头

本文的作者（夏一鸣和陈少林）发明了一种新型的 3D 方块，它就像一个具有变焦功能的智能摄像机镜头。他们称之为多分辨率六面体单元 (Multi-Resolution Hexahedron Element)。

以下是其工作原理，使用简单的类比说明：

1. “基础形状”（主模板）
不同于标准的 8 节点方块，他们创建了一个特殊的“基础节点形状函数”。想象一下，取一个标准的 8 角方块，并将其影响力延伸到更大的区域，从而创造出一个“主模板”，这个模板能够感知其周围的环境。这个模板是通过将一个角点的形状扩展到周围空间中来构建的，覆盖了 27 个潜在的点位（类似于一个 3x3x3 的点阵立方体）。

2. “分辨率等级”（变焦旋钮）
这是神奇之处所在。这种新单元有一个被称为分辨率等级 (Resolution Level, RL) 的旋钮。

• 低 RL（缩小视角）： 该单元表现得像传统的 8 节点方块。它简单且快速。
• 高 RL（放大视角）： 你转动旋钮，该单元会在不改变方块本身大小的情况下，自动“长出”更多的内部“节点”（计算点）。
• 类比： 想象一张数字照片。低分辨率的照片是模糊的；高分辨率的照片是清晰的。在这种新方法中，你不需要通过把照片切成更小的碎片来使其变清晰，你只需要调高单个方块的“分辨率等级”，它就能瞬间展现出更多细节（更多的节点）。

3. “嵌套”结构
论文解释说，这些不同细节层级的结构可以完美地嵌套在一起，就像俄罗斯套娃一样。一个高细节版本的方块包含了低细节版本的方块。这种数学上的“嵌套”确保了随着你不断放大，计算过程依然保持稳定和准确。

为什么这种方法更好？（论文的观点）

• 不再需要重划网格： 在旧方法中，为了获得更高的精度，你必须拆解结构并重建网格（re-meshing）。而在这种新方法中，你只需调整分辨率等级。这就像是调整相机的对焦，而不是重建相机本身。
• 简洁性： 这种新形状背后的数学原理非常简单。它保持了一种特殊的属性，即“克罗内克尔 $\delta$ 函数 (Kronecker delta)”（这基本上意味着“我知道我的角点确切的位置”），使得计算边界条件（例如墙壁固定的位置）变得非常容易。
• 高效性： 由于你可以用更少的方块获得更高的精度，计算机的工作量也随之减少。论文声称，这种方法比传统方法或其他使用复杂小波（wavelets）的“多分辨率”方法更快、更合理。

现实世界测试（源自论文）

作者在三种场景下测试了他们的“智能方块”：

1. 悬臂梁： 一根从墙壁伸出的梁。他们展示了其中一个“可缩放”的方块如何达到与数十个传统方块相当的精度。
2. 正方形平板： 一个平坦的板。他们将此方法与一种流行的“小波 (Wavelet)”方法进行了对比，发现他们的方法更容易使用，且精度同样高。
3. 折叠板： 一个复杂的弯曲结构。他们展示了如何能够轻松地在结构的特定部分调整细节水平，而在旧方法中，这需要一个极其庞大且复杂的网格。

核心结论

论文认为，这种全新的多分辨率六面体单元是结构分析的一种优越工具。它将“分辨率等级”视为实现精度的关键，而非增加网格碎片的数量。它声称，与传统的 8 节点单元或其他先进方法相比，该方法更加合理、易于实现且更高效，非常适合解决那些对细节要求极高的复杂工程问题（例如裂纹或应力集中点）。

注：论文提到，未来的工作将侧重于如何连接具有不同分辨率等级的方块（例如，如何将一张放大的照片与一张缩小的照片连接起来），但目前的研究结果侧重于单一单元类型的公式化构建与测试。

技术摘要

技术摘要：一种具有新型多分辨率分析的六面体单元公式化方法

问题陈述
本文探讨了计算力学中传统有限元法（FEM）及现有多分辨率分析（MRA）技术的局限性。传统有限元法被描述为“单分辨率”特征，即单个单元包含固定数量的节点。因此，提高分析精度需要进行人工划分网格和重新划分网格，作者认为这一过程是不合理的、低效的，且不适用于具有局部陡峭梯度（如材料非线性、裂纹或冲击）的问题。相反，现有的多分辨率分析方法，如小波有限元法（WFEM）、无网格法（MFM）和自然单元法（NEM），存在固有缺陷。这些缺陷包括构造节点形状函数的复杂性、缺乏克罗内克尔 δ 性质（导致边界条件处理困难），以及缺乏坚实的数学基础，这往往使得这些方法更趋向于经验主义，难以应用于工程实践。

方法论
作者提出了一种基于新型 MRA 框架的新型多分辨率六面体单元公式化方法。其核心方法论包括以下步骤：

1. 构造基础节点形状函数： 从定义在区域 $[0,1]^3$ 上的传统 8 节点六面体单元的标准等参数形状函数出发，作者将该函数扩展到区域 $[-1,1]^3$。这是通过将等参数单元围绕特定节点（位于原点）在其他七个象限内进行平移来实现的。这种扩展覆盖了 26 个相邻节点，并创建了一个连续且具有克罗内克尔 δ 性质的“基础节点形状函数”（$\varphi$）。
2. 形成相互嵌套的位移子空间序列： MRA 框架通过一系列位移子空间进行构建。这些子空间的基函数是通过对基础节点形状函数 $\varphi$ 进行缩放和平移生成的。
   • 缩放： 函数通过整数参数（$m, n, l$）进行缩放，以控制分辨率等级（RL）。
   • 平移： 缩放后的函数被平移到单元域内的不同节点位置。
   • 这一过程创建了一个子空间序列（$V_{ijk}$），其中 $V_{ijk} \subset V_{(i+1)(j+1)(k+1)}$，在不依赖小波基函数的情况下，建立了一个数学上严谨的、相互嵌套的结构。
3. 单元公式化： 单元内的位移场由分辨率等级（RL）$(m+1) \times (n+1) \times (l+1)$ 定义。总节点数直接由 RL 决定，而非由网格密度决定。单元刚度矩阵、质量矩阵和载荷向量是利用构造的形状函数，根据最小势能原理推导而出的。
4. 实现： 该方法允许自动生成节点周围的单元矩阵。通过变换矩阵将局部单元矩阵转换为全局结构坐标系。

主要贡献

• 新型 MRA 框架： 本文引入了一个用于固体力学的坚实的 MRA 数学基础，该框架不需要小波基函数，而是依赖于扩展等参数形状函数的缩放和平移。
• 可调分辨率等级（RL）： 所提出的单元具有可调的 RL，可调节单元节点数。这使得在单个单元内直接控制分析精度成为可能，与传统的重新划分网格形成对比。
• 克罗内克尔 δ 性质： 与许多现有的 MRA 方法不同，所提出的形状函数保留了克罗内克尔 δ 性质，简化了边界条件的施加。
• 单分辨率与多分辨率的统一： 作者证明了传统的 8 节点六面体单元是所提出的多分辨率单元（具体为 $m=n=l=1$ 时）的一个特例（单分辨率）。

结果
本文通过三个数值算例验证了该方法：

1. 悬臂梁： 在均匀横向载荷下，采用单单元及不同 RL 的所提单元，其端部挠度达到了与显著更高网格密度（例如 $1 \times 1 \times 30$ 网格）的传统单元相当的结果。随着 RL 的增加，结果向解析解收敛。
2. 简支正方形平板： 将所提单元与传统 8 节点单元及二维 B 样条小波插值（BSWI）单元进行了对比。结果表明，增加 RL 提高精度的效果与增加 BSWI 的阶数相似，但实现起来更简单。作者指出，由于自由度（DOF）的非理性增加以及缺乏等参数公式化，3D BSWI 单元通常过于复杂。
3. 折叠正方形板： 该方法被应用于具有自由边界和简支边界的复杂几何形状。分析表明，通过调整特定单元的 RL，可以实现对位移响应的精确建模，且与传统重新划分网格相比，所需的总单元数和变换矩阵乘法次数更少。结果与文献数据一致。

意义与主张
作者声称，所提出的多分辨率六面体单元法比传统的单分辨率 FEM 和其他现有的 MRA 方法更合理、更容易实现且计算效率更高。

• 效率： 通过避免复杂的重新划分网格和单分辨率单元的人工组装，该方法降低了计算开销。围绕节点自动获取矩阵进一步提高了效率。
• 精度控制： 论文断言，结构分析的精度由分辨率等级（RL）而非网格决定。这允许进行便捷的自适应分析，通过调整 RL 而非重建网格来调节精度。
• 适用性： 该方法特别适用于涉及局部陡峭梯度（如裂纹或冲击）的结构问题，在这些问题中需要高局部分辨率，而无需进行全局网格细化。
• 未来工作： 作者谦虚地指出，未来的工作将侧重于处理不同 RL 多分辨率单元之间的界面，可能使用桥接域或扩展的 Serendipity 单元。

总之，本文认为这种新的 MRA 框架弥合了传统 FEM 与先进多分辨率技术之间的鸿沟，为精确的结构计算提供了一个强大的工具，克服了现有方法的局限性。