[論文レビュー] $H^2$-conformal approximation of Miura surfaces

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

Advanced Materials and Mechanics被引用数 4

ひとこと要約

本稿では、均質化極限におけるミウラテッセレーション表面をモデル化する非線形楕円型PDEに対する解の存在および一意性を確立する。$H^2$-適合有限要素法とニュートン反復を組み合わせた手法を提案し、$H^2$-ノルムにおける1次収束を証明するとともに、非軸対称な非自明なミウラ表面、例えば変形された放物面や円錐型の形状のロバストな計算を実現する。

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.

研究の動機と目的

• 均質化極限におけるミウラテッセレーション形状を記述する非制約的$H^2$-適合PDEに対する解の存在および一意性を確立すること。
• 一般のディリクレ境界条件を有するこの非線形楕円型PDEの解を近似するロバストな数値法を開発すること。
• 提案された$H^2$-適合有限要素法の収束速度が、解析的および非解析的テストケースにおいて正しく評価されることを検証すること。
• 従来の手法では到達できなかった非自明で非軸対称なミウラ表面幾何形状（例：変形された放物面、円錐型の形状）を本手法がロバストに計算できることを示すこと。

提案手法

• 領域および境界データに$\mathcal{C}^{2,\alpha}$正則性を仮定したもとで、$\mathbb{R}^3$における強形式PDE $p(\varphi_x)\varphi_{xx} + q(\varphi_y)\varphi_{yy} = 0$ の解の存在および一意性を証明する。
• 解の近似に$H^2$-適合有限要素離散化を導入し、離散解の$\mathcal{C}^1$連続性を保証する。
• 得られた非線形離散系を解くためにニュートン法を採用し、高精度な近似解への収束を実現する。
• 正則性条件の下で、有限要素近似の$H^2$-ノルムにおける1次収束率が厳密に証明される。
• 構造的三角形メッシュと適応的メッシュ細分化を用いて、周期的境界条件を含む複雑な幾何形状に対応する。
• 解析解を用いた数値的収束率の検証および、境界条件が異なる非解析的表面におけるロバスト性の確認を実施する。

実験結果

リサーチクエスチョン

• RQ1 ミウラテッセレーション表面をモデル化する非制約的$H^2$-適合PDEに対して、どのような条件下で一意な解が存在するか。
• RQ2 ニュートン反復を用いた$H^2$-適合有限要素法が、この非線形楕円型PDEに対して最適収束を達成できるか。
• RQ3 本手法を用いてロバストに計算可能な非自明で非軸対称なミウラ表面形状はどのようなものか。
• RQ4 制約$|\varphi_y|^2 > 1$が自動的に満たされない場合、本手法の性能はどのように変化するか。物理的実現可能性にどのような影響を与えるか。
• RQ5 従来のアルゴリズムでは不可能であった、変形された放物面や平らに押しつぶされた円錐型の形状を本手法が計算可能か。

主な発見

• 領域およびディリクレ境界データに$\mathcal{C}^{2,\alpha}$正則性を仮定したもとで、非制約的PDEに対する$\mathcal{C}^{2,\alpha}$解の存在および一意性が証明された。
• $H^2$-適合有限要素離散化の$H^2$-ノルムにおける1次収束率が、正則性条件の下で厳密に確立された。
• 数値実験により、解析解上での理論的収束率が確認され、手法の精度が検証された。
• 本手法は、領域全体で$|\varphi_{h,y}| > 1$を満たす変形された放物面を正しく計算でき、PDEが全領域で満たされていることを確認した。
• 本手法は、$\mathcal{C}^1$ではあるが$\mathcal{C}^2$でない非軸対称な円錐型表面を生成した。$|\varphi_{h,y}| < 1$であるため、これは異なる方程式の解であることが示され、従来の分類結果と整合的であった。
• 表面が平坦化する際、$y$方向に領域が拡大すると、本手法は収束しなくなることが確認され、$|\varphi_y|^2 > 1$制約に関連する理論的限界が裏付けられた。