[论文解读] $H^2$-conformal approximation of Miura surfaces

Frédéric Marazzato|arXiv (Cornell University)|Feb 17, 2022

Advanced Materials and Mechanics被引用 4

一句话总结

本文在均匀化极限下建立了描述Miura折纸表面的无约束非线性椭圆型PDE解的存在性与唯一性。提出了一种与牛顿迭代相结合的$H^2$-符合有限元方法，证明了在$H^2$-范数下的一阶收敛性，并成功实现了对复杂非轴对称Miura表面（包括变形双曲面和锥形类似形状）的鲁棒计算。

ABSTRACT

The Miura ori is a very classical origami pattern used in numerous applications in Engineering. A study of the shapes that surfaces using this pattern can assume is still lacking. A constrained nonlinear partial differential equation (PDE) that models the possible shapes that a periodic Miura tessellation can take in the homogenization limit has been established recently and solved only in specific cases. In this paper, the existence and uniqueness of a solution to the unconstrained PDE is proved for general Dirichlet boundary conditions. Then a $H^2$-conforming discretization is introduced to approximate the solution of the PDE coupled to a Newton method to solve the associated discrete problem. A convergence proof for the method is given as well as a convergence rate. Finally, numerical experiments show the robustness of the method and that non trivial shapes can be achieved using periodic Miura tessellations.

研究动机与目标

• 在均匀化极限下，建立描述Miura折纸构型的无约束$H^2$-符合PDE解的存在性与唯一性。
• 开发一种鲁棒的数值方法，用于近似求解具有通用狄利克雷边界条件的该非线性椭圆型PDE。
• 在解析与非解析测试案例上，验证所提$H^2$-符合有限元方法的收敛速率。
• 展示该方法在计算此前难以实现的复杂非轴对称Miura表面几何形状（如变形双曲面和锥形类似形状）方面的能力。

提出的方法

• 在域与边界数据具有$\mathcal{C}^{2,\alpha}$正则性假设下，证明了强形式PDE $p(\varphi_x)\varphi_{xx} + q(\varphi_y)\varphi_{yy} = 0$在$\mathbb{R}^3$中解的存在性与唯一性。
• 提出一种$H^2$-符合有限元离散化方法，用于近似求解PDE，确保离散解具有$\mathcal{C}^1$-连续性。
• 采用牛顿法求解所得非线性离散系统，实现对高精度近似解的收敛。
• 在正则性条件下，严格证明了有限元近似在$H^2$-范数下的一阶收敛速率。
• 使用结构化三角形网格与自适应网格加密技术处理复杂几何形状，包括周期性边界条件。
• 通过解析解数值验证收敛速率，并在具有不同边界条件的非解析曲面上展示方法的鲁棒性。

实验结果

研究问题

• RQ1 在何种条件下，描述Miura折纸表面的无约束$H^2$-符合PDE存在唯一解？
• RQ2 结合牛顿迭代的$H^2$-符合有限元方法能否在该非线性椭圆型PDE上实现最优收敛？
• RQ3 该方法能够鲁棒计算哪些类型的复杂非轴对称Miura表面几何形状？
• RQ4 当约束$|\varphi_y|^2 > 1$未被自动满足时，该方法表现如何？其对物理可实现性有何影响？
• RQ5 该方法能否计算先前算法无法实现的曲面，例如变形双曲面或扁平化的锥形类似形状？

主要发现

• 在域与狄利克雷边界数据具有正则性假设的条件下，证明了无约束PDE存在唯一$\mathcal{C}^{2,\alpha}$解。
• 严格建立了$H^2$-符合有限元离散化在$H^2$-范数下的一阶收敛速率。
• 数值实验验证了在解析解上的理论收敛速率，证实了方法的精度。
• 该方法成功计算出一个变形双曲面，其在整个定义域内满足$|\varphi_{h,y}| > 1$，确认其在全定义域内满足PDE。
• 该方法生成了一个非轴对称的锥形类似曲面，其为$\mathcal{C}^1$但非$\mathcal{C}^2$，且满足$|\varphi_{h,y}| < 1$，表明其求解的是另一类方程，与先前分类结果一致。
• 当曲面趋于平坦且$y$方向域较大时，该方法无法收敛，证实了与$|\varphi_y|^2 > 1$约束相关的理论限制。