[論文レビュー] Nitsche's method method for mixed dimensional analysis: conforming and non-conforming continuum-beam and continuum-plate coupling

Vinh Phu Nguyen, Pierre Kerfriden+2 more|Open Repository and Bibliography (University of Luxembourg)|Aug 13, 2013

Advanced Numerical Analysis Techniques参考文献 68被引用数 15

ひとこと要約

本稿では、等化的および非等化的定式化を用いて、混合次元解析における固体連続体とビーム・プレートを結合するNitsche法に基づく手法を提示する。この手法は、安定的かつ一貫性のある界面条件の弱い強制と安定化パラメータを用いることで、等化的および非等化的結合を可能にし、特に複合構造物の適応的モデリングにおいて高い精度と柔軟性を示す。

ABSTRACT

A Nitche's method is presented to couple different mechanical models. They include coupling of a solid and a beam and of a solid and a plate. Both conforming and non-conforming formulations are presented. In a non-conforming for- mulation, the structure domain is overlapped by a refined solid model. Applications can be found in multi-dimensional analyses in which parts of a structure are modeled with solid elements and others are modeled using a coarser model with beam and/or plate elements. Discretisations are performed using both standard Lagrange elements and high order NURBS (Non Uniform Rational Bsplines) based isogeometric elements. We present various examples to demonstrate the performance of the method.

研究の動機と目的

• 3次元連続体と1次元ビームまたは2次元プレートとの混合次元解析（MDA）における、強固で弱い強制による結合手法の開発。
• ラグランジュ乗数、ペナルティ法、マルチポイント制約などの従来手法の限界を克服し、より安定的で柔軟な代替手法を提供すること。
• 特に重複メッシュによる離散化を支援する、等化的および非等化的結合を可能にする。
• 高精度な固体要素と次元削減されたビーム／プレートモデルの効率的結合を可能にすることで、等化的解析（IGA）の適用範囲を大規模複合構造物に拡大すること。
• 剛性行列の正定値性を保証するための安定化パラメータの理論的基礎を提供すること。

提案手法

• Nitsche法の定式化を用いて、連続体とビーム／プレートモデル間の界面で変位の連続性と応力の釣り合いを弱く強制する。
• ラグランジュ乗数の代わりに物理的フラックスを用い、双線形形式の安定化と強制性の確保のための一貫したペナルティ項を追加する。
• より高い精度と滑らかさを実現するため、標準ラグランジュ有限要素と高次NURBSに基づく等化的要素の両方を用いて実装する。
• 非等化的結合のため、固体領域をビームまたはプレート領域と重ね合わせ、レベルセット法を用いて要素の交差を検出し、界面項を正確に積分する。
• 双線形形式の数値解析を実施し、グローバル剛性行列の正定値性を保証するための安定化パラメータの理論的下界を導出する。
• 古典的および一次剪切変形理論に基づくビーム／プレート理論の両方をサポートし、厚さのある構造物のモデリングを可能にする。

実験結果

リサーチクエスチョン

• RQ1 混合次元解析において、3次元固体要素と1次元ビームまたは2次元プレートモデルとの間で、どのように安定的かつ一貫性のある結合を達成できるか？
• RQ2 混合次元結合における剛性行列の正定値性を保証するための、Nitsche法における安定化パラメータの最適選択は何か？
• RQ3 モルタル法、ペナルティ法、MPCなどの従来手法と比較して、Nitsche法に基づく結合の精度、安定性、実装の複雑さの観点での性能はいかがなものか？
• RQ4 重複メッシュを用いた非等化的結合は、複合構造物の適応的モデリングにおいて、効率的に実装可能か？
• RQ5 高次NURBS要素の使用は、混合次元問題の精度と収束性をどの程度向上させるか？

主な発見

• 提案手法は、線形ラグランジュ要素および高次NURBSに基づく等化的離散化において、混合次元問題で最適収束率を達成する。
• 数値例では、結合された固体プレートおよび固体ビームモデルと参照解との間に優れた一致が得られ、ベンチマークケースでは変位および応力の相対誤差が1%未満であった。
• 安定化パラメータは理論的に下限値が導出され、実験的数値（例：$1 imes 10^6$）が安定かつ正確な結果をもたらすことが確認された。
• 重複メッシュを用いた非等化的結合により、メッシュ再構築なしに固体領域をシフトさせても正確な結果を維持できる柔軟なモデル適応性が実証された。
• ラグランジュ乗数空間の必要性やペナルティパラメータのチューニングを回避する点で、従来手法に比べ実装の簡便性に優れている。
• 回転自由のキルホフNURBSプレート要素とNitsche結合を組み合わせることで、自由度を最小限に抑えつつ、薄板の高精度なモデリングが可能になった。