Augmented Lagrangian法
Augmented Lagrangian法の理論基礎
Augmented Lagrangian法とは
先生、Augmented Lagrangian法はペナルティ法とLagrange乗数法の「いいとこ取り」ですか？
まさにそう。ペナルティ法の簡単さとLagrange乗数法の貫通ゼロ精度を組み合わせた手法。
ペナルティ項 $k_p g_n$ に加えてLagrange乗数 $\lambda_n$ を反復的に更新する。反復を繰り返すと貫通がゼロに収束。
ペナルティ剛性 $k_p$ が小さくても、反復で貫通を減らせるんですね。
$k_p$ への依存性がペナルティ法より小さい。これが最大の利点。AnsysのデフォルトがAugmented Lagrangian。
まとめ
要点：
- ペナルティ＋反復的Lagrange更新 — いいとこ取り
- $k_p$への依存が小さい — ペナルティ法より安定
- 追加DOF不要 — Lagrange乗数法より効率的
- Ansysのデフォルト — 最も広く推奨される手法
Hestenes-Powell法1969年
拡張Lagrangian法（Augmented Lagrangian）は、1969年にM.R. HestenesとM.J.D. Powellが独立に発表したペナルティとLagrange乗数のハイブリッド手法だ。乗数を外部反復ごとに更新することで、純ペナルティ法の条件数悪化を回避しつつ、純乗数法の大規模連立方程式を回避できる。接触FEMへの応用は1980年代後半にSimoとLaursenが体系化した。
Augmented Lagrangian法の数値計算手法
Augmented Lagrangianの実装
アルゴリズム：
1. 初期の$\lambda = 0$でペナルティ法として解く
2. 貫通量 $g_n$ を確認
3. $\lambda$ を更新: $\lambda \leftarrow \lambda + k_p g_n$
4. 更新した$\lambda$で再度解く
5. $g_n$ が十分小さくなるまで反復
外側のループ（$\lambda$の更新）と内側のループ（Newton-Raphson）があるんですね。
二重の反復ループ。内側で平衡、外側で接触拘束を満たす。計算コストはペナルティ法の1.5〜2倍程度。
ソルバー設定
まとめ
乗数更新の反復スキーム
拡張Lagrangian接触では、外部ループでLagrange乗数λを更新し、内部ループで非線形有限要素解を解く2重反復構造を持つ。Laursen & Simho(1993年)の定式化では、乗数更新式 λ_{k+1} = λ_k + ε_N g_N（g_N:浸透量、ε_N:ペナルティ）を採用し、ε_Nを物理的なギャップ許容値から逆算して設定できるため、ペナルティ法よりも収束チューニングが直感的に行える。
Augmented Lagrangian法の実務適用
Augmented Lagrangianの実務
Ansysユーザーは無意識にAugmented Lagrangianを使っていることが多い。デフォルトが最も安定。
実務チェックリスト
航空エンジン翼フレッティング
GE Aviationは2008年頃からANSYS Mechanicalの拡張Lagrangian接触を用いてタービンブレードのダブテール部フレッティング摩耗解析を実施している。サイクル数10⁷回相当の繰り返し接触を累積損傷モデルと組み合わせて評価し、インコネル718製ブレードの摩耗深さを実測値±15μmの精度で再現。試験品製作コストを約30%削減した事例として学会発表された。
Augmented Lagrangian法のソフトウェア比較
Augmented Lagrangianのツール
選定ガイド
ANSYS ALM実装の変遷
ANSYSは1990年代後半にAugmented Lagrange Method（ALM）をCONTA174/TARGE170要素のデフォルト接触アルゴリズムとして採用した。ANSYS 10.0（2005年）ではALMの収束判定基準が接触力ベースに改訂され、厚板の浸透問題が大幅に改善した。現行のANSYS Mechanical 2024では自動的にALMとペナルティ法を切り替える「program controlled」設定がデフォルトとなっている。
Augmented Lagrangian法の先端研究
Augmented Lagrangianの先端
変分的接触の最前線
2020年代に入り、拡張Lagrangian法をisogeometric analysis（IGA）と組み合わせた変分的接触定式化が活発に研究されている。Temizer(2022年)らはNURBSベースの接触面表現とALMを統合し、接触圧力の空間分布をC¹連続で表現することに成功した。従来要素では現れていた接触圧力のギザギザ（oscillation）が消失し、フレッティング疲労寿命予測の精度が実験値との差異20%以内に収まると報告されている。
Augmented Lagrangian法のトラブル対応
Augmented Lagrangianのトラブル
乗数更新発散の対処
拡張Lagrangian法の落とし穴は、ペナルティε_Nを大きく設定しすぎると乗数更新が振動発散することだ。2010年代のNastran SOL 401ユーザーから「乗数が10サイクルで±10⁸まで振れる」という報告が複数あった。対策はε_Nを接触面の局所剛性（E/h）の1/10以下に抑えること、かつ外部ループ収束判定にギャップ残差ノルム基準を加えることで、ほぼ確実に回避できる。
関連トピック
なった
詳しく
報告