[論文レビュー] Construction of analysis-suitable $G^1$ planar multi-patch parameterizations

Mario Kapl, Giancarlo Sangalli+1 more|arXiv (Cornell University)|Jun 10, 2017

Advanced Numerical Analysis Techniques参考文献 24被引用数 3

ひとこと要約

本稿では、与えられた初期C0マルチパッチ幾何から、平面領域の解析に適したG1（AS-G1）マルチパッチパラメータ化を構築する手法を提示する。線形制約を伴う二次計画問題を解くことで、境界、頂点、および頂点における一次微分を保持しつつ、元の幾何にできるだけ近いAS-G1パラメータ化を生成する。数値実験では最適なh収束が確認された。

ABSTRACT

Isogeometric analysis allows to define shape functions of global $C^{1}$ continuity (or of higher continuity) over multi-patch geometries. The construction of such $C^{1}$-smooth isogeometric functions is a non-trivial task and requires particular multi-patch parameterizations, so-called analysis-suitable $G^{1}$ (in short, AS-$G^{1}$) parameterizations, to ensure that the resulting $C^{1}$ isogeometric spaces possess optimal approximation properties, cf. [7]. In this work, we show through examples that it is possible to construct AS-$G^{1}$ multi-patch parameterizations of planar domains, given their boundary. More precisely, given a generic multi-patch geometry, we generate an AS-$G^{1}$ multi-patch parameterization possessing the same boundary, the same vertices and the same first derivatives at the vertices, and which is as close as possible to this initial geometry. Our algorithm is based on a quadratic optimization problem with linear side constraints. Numerical tests also confirm that $C^{1}$ isogeometric spaces over AS-$G^{1}$ multi-patch parameterized domains converge optimally under mesh refinement, while for generic parameterizations the convergence order is severely reduced.

研究の動機と目的

• 平面領域の解析に適したG1（AS-G1）マルチパッチパラメータ化を構築する実用的手法を開発すること。
• 得られるパラメータ化が、初期C0マルチパッチ幾何の境界、頂点、および頂点における一次微分を保持すること。
• AS-G1制約を満たしつつ、初期幾何にできるだけ近いAS-G1パラメータ化を生成すること。
• C1等材性的空間がAS-G1パラメータ化上で最適なh収束率を達成することを数値的に検証すること。
• 摂動に対して感度を示すが、それでもAS-G1幾何の構築の柔軟性と実現可能性を示すこと。

提案手法

• AS-G1パラメータ化の構築を、線形等式制約を伴う二次計画問題として定式化する。
• 初期C0マルチパッチパラメータ化を出発点とし、その境界、頂点、および頂点における一次微分を保持する。
• パッチ界面における幾何的連続性条件から導かれるAS-G1制約を課す。
• 制約の感度が高いため、数値的精度を維持するため、記号計算を用いて最適化問題を数値的に解く。
• パッチごとに順次パラメータ化を構築し、AS-G1条件により界面をC1滑らかに保つ。
• 将来の研究において、条件数の改善と感度の低減を図るため、局所的サポートを持つ基底を採用する。

実験結果

リサーチクエスチョン

• RQ1 境界のみが与えられた任意の平面マルチパッチ領域に対して、AS-G1マルチパッチパラメータ化を構築できるか？
• RQ2 AS-G1制約は、初期幾何に近い品質と近接性をどの程度制限するか？
• RQ3 AS-G1パラメータ化の使用が、等材性的解析における最適なh収束を保証するか？
• RQ4 AS-G1制約は微小な摂動に対してどの程度感度を示すか？また、基底設計によってこれを緩和できるか？
• RQ5 提案手法は体積的または表面ベースのマルチパッチ領域へ拡張可能か？

主な発見

• 提案手法は、視覚的に初期C0幾何と区別がつかないが、すべての必要な連続性制約を満たすAS-G1マルチパッチパラメータ化を効果的に生成した。
• 数値実験により、AS-G1パラメータ化上でのC1等材性的空間が、L2ノルムにおいて最適な収束率O(h⁴)を達成することが確認された。
• これに対して、一般的なC0パラメータ化は収束率が著しく低下し、L2ノルムではO(h¹/²)にまで低下し、L∞ノルムでは収束しないことが示された。
• AS-G1制約は高品質なパラメータ化の構築を妨げず、幾何設計における顕著な柔軟性を示した。
• 初期幾何がAS-G1でなくても、複数のテストケース（自動車、パズルピース、変形円など）を通じて、本手法が頑健で効果的であることが示された。
• AS-G1制約の高い感度を考慮し、記号計算が不可欠であることが明らかになった。