Caveats on formulating finite elasto-plasticity in curvilinear coordinates

この論文は、有限ひずみ・有限変形を伴う軸対称弾性塑性問題の曲線座標系における数値実装において、変形勾配やヤコビアン、シフターなどの取り扱いに生じる追加項を明確化し、微分幾何学的手法ではなく標準的な直交座標系からの基底変換を用いた実用的な有限要素法への定式化手順を提示するものである。

要約

この論文は、**「丸い形（円筒形）のものを、巨大な力や変形を伴ってシミュレーションする際、数学の落とし穴にハマらないための注意書き」**です。

専門用語を捨てて、日常の例え話を使って解説しますね。

1. 背景：なぜ「丸い形」は特別なのか？

私たちが普段、建物の強度や土の動きをコンピューターで計算する時（有限要素法といいます）、世界を「箱（直方体）」の集まりとして考えるのが一番簡単です。これは**「直交座標（カーテンの格子）」**のようなものです。

しかし、現実には**「トンネル」「冷却塔」「土の円柱」など、丸い形（円筒形）のものがたくさんあります。これらを「箱」で無理やり表現すると、計算が非効率になります。そこで、「円筒座標（半径、角度、高さ）」**という、丸い形に合わせたルールを使おうとします。

ここが問題の始まりです。
「箱」のルール（直交座標）をそのまま「丸い形」に当てはめようとすると、**「見えない追加のルール」**が突然現れて、計算が破綻してしまいます。

2. 核心：「変形勾配」という魔法の鏡

この論文が扱っている最大のテーマは**「変形勾配（Deformation Gradient）」という概念です。これを「物体を拡大・縮小・回転させる魔法の鏡」**だと想像してください。

• 直交座標（箱）の場合：
  この鏡はシンプルです。「横に 2 倍、縦に 1 倍」のように、単純な数字の掛け算で済みます。
• 円筒座標（丸）の場合：
  ここで**「しわ寄せ（Shifter）」という見えない役者が現れます。
  丸い形では、半径が広がると「円周の長さ」も自然に変わります。この「半径が変わることで生じる、見えない長さの変化」を補正する役者が「しわ寄せ」**です。

論文の重要な発見：
多くの研究者は、この「しわ寄せ」の存在を忘れがちです。
「直交座標の計算式を丸い形にそのままコピー＆ペーストしただけ」でプログラムを作ると、**「変形したはずの物体の体積が、計算上は消えたり、増えたりする」**というバグが起きます。
例えば、土の円柱を圧縮したのに、計算上は体積が増えたり、逆に消えたりしてしまうのです。

3. プラスチック変形：「元に戻らない」泥沼

さらに、この論文は**「塑性（Plasticity）」**という現象にも触れています。

• 弾性（Elastic）： ゴムのように、力を抜くと元に戻る状態。
• 塑性（Plastic）： 粘土のように、力を抜いても元に戻らない状態。

この「元に戻らない」状態を計算する時、「中間の仮の姿（ストレスフリーな状態）」という、目に見えない架空の空間を想像する必要があります。
円筒座標では、この「架空の空間」のルール（計量テンソル）も、通常の箱とは異なります。論文は、この「見えない空間」と「現在の空間」の関係を正しく結びつけるための、ステップバイステップのレシピを提供しています。

4. 具体的な解決策：論文が提案していること

著者たちは、難しい微分幾何学（曲がった空間の数学）を使うのではなく、**「あえて直交座標（箱）の考え方に戻り、そこから円筒座標（丸）への『翻訳』を丁寧に行う」**という実用的なアプローチを取りました。

彼らが提案する手順は以下の通りです：

1. 変形を計算する時： 単純な微分ではなく、「クリストッフェル記号」という**「曲がった道のりの補正係数」**を必ず足す。
2. 体積を計算する時： 単純な行列の式ではなく、「計量係数（半径の長さなど）」を掛け合わせて補正する。
3. 応力（力）を計算する時： 「ヤコビアン（体積変化率）」を、弾性部分と塑性部分に分けて正しく計算する。

5. 検証：実際に試してみた結果

論文の最後には、**「厚い円筒の壁」**をシミュレーションする実験があります。

• 方法 A： 3 次元全体を「箱」で計算する（非常に時間がかかる）。
• 方法 B： この論文の新しいルールを使って、2 次元の「円筒断面」だけで計算する（非常に高速）。

結果：
両者の計算結果は**「ほぼ完全に一致」しました。
つまり、「この論文のルールさえ守れば、3 次元の重たい計算を、2 次元の軽い計算で、同じ精度で済ませられる」**ことが証明されたのです。

まとめ：この論文が私たちに教えてくれること

この論文は、**「丸い形を計算する時は、直線の感覚（直交座標）をそのまま使うと痛い目を見る」**と警告しています。

• アナロジー：
  地球儀（丸い世界）の地図を引く時、平面の紙（直交座標）のルールをそのまま使うと、グリーンランドがアフリカより大きく見えてしまいます（メルカトル図法の歪み）。
  これと同じで、**「丸い物体の変形を計算する時、平面のルールをそのまま使うと、体積や力が歪んで計算されてしまう」**のです。

著者たちは、この歪みを正しく補正するための**「地球儀用の正しい地図の描き方（数式）」**を、エンジニアが実際にプログラムに落とし込めるように、具体的かつ丁寧に解説してくれました。

これにより、トンネルやダム、生体組織のシミュレーションなどが、より正確に、かつ低コストで計算できるようになるのです。

技術概要

以下は、提示された論文「Caveats on formulating finite elasto-plasticity in curvilinear coordinates（曲線座標系における有限弾塑性の定式化に関する留意点）」の技術的サマリーです。

1. 問題の背景と課題

固体力学における数値解析（特に有限要素法：FEM）では、計算コストを削減するために、平面ひずみ、平面応力、軸対称の 3 つの 2 次元近似が頻繁に用いられています。軸対称問題は、地盤工学（杭基礎、トンネル掘削）、土木構造物（冷却塔、サイロ）、金属加工、生体力学など広範な分野で重要です。

しかし、有限ひずみ（大変形）を伴う弾塑性解析を軸対称条件下で実施する際、以下の課題が存在します。

• 座標系の違いによる誤解: 連続体力学はテンソル解析に基づき座標不変性を持っていますが、FEM の実装はユーザーが選択する基底（通常は直交座標）における行列表現に依存します。
• 曲線座標系の複雑さ: 軸対称解析では円筒座標（曲線座標）が用いられます。直交座標系では自明または暗黙的に扱われる量（変形勾配、ヤコビアン、シフターなど）が、非直交・非正規化された基底を持つ曲線座標系では、追加の項や補正を必要とします。
• 弾塑性の難しさ: 等積（isochoric）な塑性変形を考慮する場合、応力状態は現在の弾性変形だけでなく、中間構成（中間状態）の幾何学的変化にも依存します。これを曲線座標系で正しく定式化し、数値的に安定した線形化（consistent linearisation）を行うことは容易ではありません。

既存の多くの実装では、直交座標系の定式化を単純に曲線座標系に拡張しようとしており、これにより変形勾配や体積変化の定義に誤りが生じ、FEM プログラムの根本的な不整合を引き起こすリスクがあります。

2. 手法とアプローチ

本論文は、幾何学的な多様体（Riemannian manifold）の枠組みを採用するのではなく、標準的な直交座標表現を基盤とし、明示的な基底変換（change of basis）を用いて曲線座標形式を導出するという実用的なアプローチをとっています。

• モデル問題: 空洞膨張問題（cavity expansion）を初期モデルとして、変形勾配、ヤコビアン、シフター（shifter）の役割を明確化しました。
• 運動学の定式化:
  • 変形勾配（Deformation Gradient）: 曲線座標系における成分の定義を慎重に行い、直交座標系との関係を明確にしました。特に、変形勾配が単なるベクトルの微分（勾配）ではなく、点から点への写像の微分であることを強調し、その転置や逆行列の計算において計量テンソル（metric coefficients）の重要性を指摘しました。
  • シフター（Shifter）: 異なる構成（初期状態と現在状態）の基底間を結びつけるテンソルとして定義し、軸対称条件下でのその挙動を解析しました。
  • ヤコビアン（Jacobian）: 体積変化を表すヤコビアンは、変形勾配の行列式だけでなく、計量係数の行列式を含む補正項が必要であることを示しました。
• 弾塑性の定式化:
  • 変形勾配の乗法的分解（$F = F^e \cdot F^p$）に基づき、中間構成のメトリックを定義しました。
  • 弾性ヤコビアン（$J_e$）と塑性ヤコビアン（$J_p$）を区別し、特に塑性変形が等積である場合の Cauchy 応力の依存関係を扱いました。
  • 対数ひずみ（logarithmic strain）と弹性左 Cauchy-Green 変形テンソルを用いることで、ヤコビアンの線形化を効率的かつ正確に行う手法を提案しました。
• 平衡方程式と弱形式:
  • 一般曲線座標系における応力測度（Cauchy 応力、第 1 ピオラ - キルヒホフ応力）と運動量保存則の強形式・弱形式を導出しました。
  • これらを軸対称条件に特化し、2 次元の弱形式（積分に $2\pi r$ の項が現れる）を導きました。
• 数値実装:
  • 更新ラグランジュ（UL）法と全ラグランジュ（TL）法の両方のアルゴリズムフローを提示し、ニュートン - ラプソン法（NR 法）による反復解法における線形化の詳細（共変微分の扱いなど）を記述しました。

3. 主要な貢献

1. 曲線座標系における運動量の厳密な定義: 変形勾配、その転置・逆行列、ヤコビアン、シフターが曲線座標系でどのように定義・計算されるべきかを、空洞膨張問題を通じて詳細に解説しました。特に、直交座標系からの単純な拡張では生じる誤差（例えば、変形勾配の非対角成分やヤコビアンの値の誤り）を指摘しました。
2. 有限ひずみ弾塑性の整合的線形化: 塑性変形を含む場合のヤコビアンの分解と、それに基づく応力 - ひずみ関係の線形化手法を提案しました。これにより、ニュートン - ラプソン法の収束性を保ちつつ、大変形弾塑性解析を軸対称条件下で安定的に実行可能にしました。
3. 実用的な FEM 実装ガイド: 幾何学的な抽象論に頼らず、実際の FEM コード実装に直結するステップバイステップの手順とアルゴリズム（更新ラグランジュ法と全ラグランジュ法の両方）を提供しました。
4. 3 次元解析との検証: 厚肉円筒の空洞膨張問題において、提案された軸対称 2 次元定式化と完全な 3 次元解析の結果を比較しました。

4. 結果

• 精度の検証: 厚肉円筒の数値例において、提案された軸対称 2 次元定式化と 3 次元解析の結果（Cauchy 応力の不変量、ヤコビアン、収束挙動）を比較しました。その結果、両者の間に実質的な差異は見られず、提案手法が 3 次元解析と同等の精度を持つことが確認されました。
• 計算効率: 軸対称定式化は 2 次元メッシュのみを使用するため、3 次元解析に比べて計算コストが大幅に低減されます。
• 収束性: 提案された線形化手法を用いたニュートン - ラプソン法は、非線形性の高い大変形問題においても、3 次元解析と同様に安定して収束しました（最大 7 回、平均約 3.5 回の反復で収束）。

5. 意義と結論

本論文は、軸対称有限ひずみ弾塑性解析を曲線座標系で実施する際に、**「直交座標系の定式化を安易に拡張することの危険性」を明確に示し、「基底変換と計量テンソルを正しく扱うことの重要性」**を強調しています。

• 理論的意義: 変形勾配が単なる微分ではなく、点間の写像の微分であるという本質的な理解に基づき、曲線座標系における運動量と応力の関係を再構築しました。
• 実用的意義: 地盤工学や構造工学など、軸対称大変形問題に直面するエンジニアや研究者に対し、誤りのない FEM 実装のための具体的な指針とアルゴリズムを提供しました。
• 結論: 曲線座標系、特に軸対称条件における有限ひずみ解析では、変形勾配、ヤコビアン、シフターなどの幾何学的量の定義を慎重かつ正確に行うことが不可欠です。これらの定義の誤りは、その後の微分や応力測度の計算に誤差を蓄積させ、全ラグランジュ法と更新ラグランジュ法の間の不整合や、数値的不安定性を引き起こす可能性があります。本論文で提示された手法は、これらの課題を解決し、ロバストな有限要素解析を可能にします。