ラグランジュ乗数法 — CAE用語解説

良い質問だ。最も典型的なのは「接触問題」だ。例えば、2つの部品が押し合っている時、互いにめり込んではいけないという「非貫通条件」が制約となる。この条件を厳密に満たすためにラグランジュ乗数法が使われる。乗数は、接触面に生じる「接触圧力」そのものと解釈できる。

変分原理で考えるとわかりやすい。弾性体のポテンシャルエネルギーΠ(u)が最小になる変位uを求めるが、接触面Γ_cではギャップg(u) ≤ 0という制約が付く。ラグランジュ関数は

その通りで、これが接触問題が非線形で難しい理由の一つだ。ソルバーは反復計算の中で、各接触点が「接触中」か「離れている」かを判定し、条件を満たすλとuを探す。このアルゴリズムとして、Abaqus/Standardで使われる「ラグランジュ乗数法」や、Ansys Mechanicalの「Augmented Lagrange法」が有名だ。

決定的な違いは、制約の厳密さと収束性のトレードオフだ。純粋なラグランジュ乗数法は非貫通条件を（数値誤差の範囲内で）厳密に満たすが、全体の剛性行列にゼロ対角成分が現れ、方程式がやや解きにくくなる。一方、Augmented Lagrange法はペナルティ項を追加した関数

モデルの全体サイズと要素サイズに依存するが、例えば10mm角の部品で要素サイズ0.5mmの場合、許容貫通量は1e-3 mm ~ 1e-5 mm程度に制御されることが多い。多くの実務的な構造解析では、このオーダーの誤差は応力結果に与える影響が無視できる。Ansysのデフォルト設定では、このAugmented Lagrange法が採用されている。

二つのケースが考えられる。第一に、流体を含む密閉容器のシール解析など、幾何学的な隙間が物理的に絶対に許されない場合。第二に、過剰な拘束による剛体モードを正確に扱う必要がある場合だ。例えば、ベアリングのボールとレース間の接触では、Augmented Lagrange法では微小なすべりが生じ、回転摩擦の計算に影響する可能性がある。Abaqus/Standardの「ラグランジュ乗数法（デフォルト）」はこうしたケースで選択される。

いや、いくつか対策がある。まず、接触面の定義を見直せ。要素サーフェスが滑らかで、主従面のメッシュサイズが極端に違わないか（理想は1:1、許容は1:10程度）。次に、接触剛性（ペナルティ剛性）を調整してみろ。Ansysでは「FKN」パラメータ、Abaqusでは「接触剛性スケールファクタ」だ。デフォルトの1.0から0.1や10.0に変えると収束が改善することもある。

その通りトレードオフだ。経験則として、接触面近傍の最も小さな要素寸法の1%以下の貫通量であれば、応力分布への影響は小さい。まずは0.1に下げて実行し、ソフトが出力する貫通量の監視結果（Abaqusの`.msg`ファイル、Ansysのソルバー出力）を確認するんだ。貫通量が許容範囲内で収束するならOK。もし貫通が大きすぎるなら、0.5など中間値で試行錯誤する。

最低3つは確認する。1. **接触圧力（CPRESS）の連続性と特異点**： 圧力分布が物理的か（鋭いピークはメッシュ依存の可能性）。2. **接触状態（CSTATUS）**： 想定通り接触/離脱している領域か。3. **貫通量/ギャップ量（COPEN）**： 許容範囲内か。特に、Abaqusでは「接触制約」の満足度を確認する。これらは、ラグランジュ乗数法が適切に機能しているかの直接的な証拠だ。

用語とデフォルト設定が大きく異なる。**Abaqus/Standard**では、「ラグランジュ乗数法」が面-面接触のデフォルトだ。これは純粋に近いラグランジュ乗数法で、制約を厳密に扱う。一方、**Ansys Mechanical**のデフォルトは「Augmented Lagrange法」であり、ソフト内では単に「ラグランジュ」と表示されることもある。Ansysで純粋なラグランジュ乗数法に近い挙動をさせたい場合は、「純ペナルティ法」を選択し、非常に高いペナルティ剛性を設定するという回り道が必要だ。

COMSOLの接触は「ペナルティ法」が基本だが、「Nitsche法」という選択肢もある。Nitsche法は、ラグランジュ乗数法のように追加変数（乗数）を導入せずに弱形式で制約を課す方法で、安定性と精度のバランスが良い。流体-構造連成（FSI）で移動する界面の「スライディングメッシュ」条件など、マルチフィジックス問題に適応しやすいのが特徴だ。ただし、AnsysやAbaqusほど接触アルゴリズムのチューニングオプションは豊富ではない。

「優れている」ではなく「適している」が正しい。**Abaqus/Standard**は、自動車のボルト締結やゴムシールなど、産業界で要求される厳密な接触解析の実績が豊富で、アルゴリズムが堅牢。**Ansys Mechanical**のAugmented Lagrange法は、収束性が良く、複雑な大変形接触（例えば、ハイパーエラストマーの摺動）で使いやすい。**COMSOL**は、熱接触抵抗や電気接触抵抗など、物理場を組み合わせた特殊な接触条件を柔軟に定義したい研究開発向け、と言えるだろう。

むしろ、ラグランジュ乗数法の「長所」が引き起こす現象だ。この方法は制約を厳密に満たそうとするため、例えば1つの節点を複数の主面が同時に拘束しようとすると、拘束条件が重複（過剰拘束）し、数値的に特異な方程式系が生じる。ペナルティ法ではある程度の「しなり」があるので、このエラーは出にくい。

まず、エラーが出た節点の位置を特定する。よくある原因は二つだ。1. **凹角部分での自己接触**： 鋭い凹角では、一つの従節点が隣接する二つの主面セグメントに同時に接触判定される。対策は、接触面を分割して定義するか、微少な面取りを幾何学的に追加する。2. **対称面と接触面の交線**： 対称境界条件と接触条件が同じ節点を拘束する。対策は、接触面の定義から対称面にごく近い節点を除外する（Abaqusなら`Adjust`オプション）。

その通り、物理的に不適だ。これは数値計算の不安定性のサインで、通常、反復途中の一時的な現象だ。主な原因は、接触剛性が低すぎる、または荷重増分が大きすぎて接触状態の急激な変化が起きていることだ。対策として、荷重ステップをより細かく分割する（Ansysの`Auto Time Stepping`、Abaqusの`Stabilization`を利用）、または接触剛性（FKN等）を1段階上げて再計算してみる。それでも解消しない場合は、モデル自体に（例えば初期ギャップが負＝初期貫通があるなど）問題がないか疑う。

良い指摘だ。乗数は接触面の自由度として追加されるため、接触面の節点数が非常に多い（例えば10万節点）と、無視できないメモリと計算コストが増加する。対策は二つ。第一に、接触面を定義する際、本当に接触する可能性がある領域だけを選択する。第二に、ソルバー設定を見直す。Abaqus/Standardで直接法ソルバーを使っている場合、疎行列直接法（例えば`PARALLEL=ALL`）は乗数によるゼロ対角成分に比較的強い。反復法ソルバー（`ITERATIVE`）を選択する場合は、前処理オプションのチューニングが必要になるかもしれない。