[論文レビュー] A dual lumping procedure for static condensation in mixed NURBS-based isogeometric elements with optimal convergence rates for arbitrary open knot vectors

Lisa Stammen, Wolfgang Dornisch|arXiv (Cornell University)|Feb 25, 2026

Advanced Numerical Analysis Techniques被引用数 0

ひとこと要約

要約: 本研究は、強化された近似デュアル基底関数を用いてせん断パラメータ行列を lump することで混合はさみジオメトリック板の静的凝縮を効率化し、内部連続性が限定された knot ベクトルにも対しても最適収束率を達成する。

ABSTRACT

Locking is a common effect in finite element and isogeometric analysis. In the case of plates, transverse shear locking is most prominent, for shells several other types of locking exist. A common cure are mixed methods that introduce additional fields of unknowns into the variational formulation. These fields reduce constraints and thus alleviate locking significantly. As a drawback, the discretized additional fields increase computational costs significantly. These fields are often eliminated by static condensation, which requires the inverse of a part of the stiffness matrix. In Lagrange-based finite elements, this inverse is computed on element level, due to a discontinuous interpolation of additional fields. Since isogeometric analysis features higher continuity, static condensation must be performed on patch level, which requires a costly matrix inversion on that level. In this contribution, the virtual shear parameters of a mixed isogeometric plate formulation are interpolated by enhanced approximate dual basis functions. This allows to conduct row-sum lumping of the relevant matrix part at a minimal loss of accuracy, since this part becomes diagonal dominant. For a properly chosen integration space, this lumped matrix becomes the identity matrix. Thus, the proposed condensation procedure does not require an inversion anymore. The crucial and novel point is the proposed treatment of knot vectors with limited internal continuity. With the help of several single- and multi-patch examples, both with full and with limited internal continuity, we show that the proposed procedure obtains optimal error convergence rates in all cases, while without these alterations, convergence rates are significantly deteriorated.

研究の動機と目的

• 再構成鎖れた Reissner–Mindlin はおのじ​​混合形式を用いたせん断ロックの解消。
• 大規模パッチレベルの行列を逆行列化せずに凝縮技術を開発。
• 最小限の精度低下で対角優勢な lumping を可能にする拡張近似デュアル基底関数を導入。
• 内部連続性が限定された単一および複数パッチ離散化を調査。
• 実用的な knot ベクトル設定で最適収束率を維持するガイドラインを提供。

提案手法

• 拡張近似デュアル基底関数で仮想せん断パラメータを補間し、関連行列部分の行和 lumping を可能にする。
• 適切な積分空間下で単位行列になる可能性のある対角優勢な凝縮行列を実現。
• NURBS のデュアル基底フレームワークを用い、近似デュアル（1変量および2変量）へ射影する変換マトリクスを構築。
• Xi+ を適応させて knot ベクトルの内部連続性が限定的な場合に最適収束率を回復。
• C0 および C1 連続性を持つ knot を扱い、複数パッチ離散化下の収束を保持するために拡張近似デュアル基底を導入。

実験結果

リサーチクエスチョン

• RQ1 近接デュアル基底モデルに基づく lumping アプローチは、パッチレベル行列の逆行列化なしに混合はさみジオメトリック板の効率的な静的凝縮を実現できるか。
• RQ2 拡張近似デュアル基底関数は、内部連続性が限定的な knot ベクトルに対して単一・複数パッチ設定で最適収束率を回復するか。
• RQ3 knot の連続性が収束に与える影響は何であり、提案フレームワークでそれをどのように緩和できるか。
• RQ4 せん断パラメータの補間における重みを大幅な精度低下なしに削除することは可能か。
• RQ5 このアプローチは他のスプライン基底（例：T-Splines）やパッチ結合にも拡張可能か。

主な発見

• 標準 NURBS 補間を用いた row-sum lumping は大きな誤差と収束低下をもたらす。
• 拡張近似デュアル基底関数は内部連続性が限定的な knot ベクトルに対して最適収束率を回復する。
• Knot の連続性を Xi+ の多重性を増やして減少させると、多くのケースで global-stress の最適収束を回復。
• 拡張デュアルは C0 連続性 knot に対して補正効果をもち、複数パッチ離散化で精度低下を回避。
• せん断パラメータ補間の重みは、効率的な実装のために最小限の accuracy loss で削除可能。
• この手法は他のスプライン基底（例：T-Splines）やパッチ結合シナリオへの適用可能性を示唆。