エネルギーバランスエラー

非線形構造解析、特に弾塑性解析において、外力が物体にする仕事（外部仕事）と、物体内部に蓄えられるひずみエネルギーと塑性散逸エネルギーの合計が一致していないことを示します。例えば、Ansys Mechanicalでは、エネルギーバランス誤差が5%を超えると結果の信頼性に疑問が生じるとされています。式で書くと、

塑性散逸エネルギーは、材料が塑性変形する際に、熱や内部組織の変化として不可逆的に消費されるエネルギーです。一方、弾性ひずみエネルギーは、バネのように蓄えられ、除荷時に回収可能です。例えば、S45C鋼を降伏応力355MPaを超えて引張ると、塑性変形部分ではエネルギーが熱に変わります。JIS G 4051のS45C材料データを用いた解析では、この散逸エネルギーの見積もりが重要です。

数値誤差以外にも、接触問題における摩擦によるエネルギーの散逸や、減衰（レイリー減衰など）のモデル化が不適切な場合、物理的にエネルギーが「抜け落ちて」見えることがあります。例えば、接触面で設定した摩擦係数μ=0.2が実際の挙動と合っていないと、摩擦熱として失われるエネルギーが見積もり誤差となり、バランスエラーに現れます。

密接に関係しています。多くのFEMソルバー（Abaqus/Standardなど）では、各増分ステップの反復計算において、力の残差（Force Residual）とエネルギーの増分（Energy Increment）の両方で収束を判定します。エネルギーバランスエラーは、全ステップを通じた積分値の誤差です。例えば、Abaqusのデフォルト収束基準は、力の残差が平均反力の0.5%、エネルギー増分が全エネルギーの1%以下です。このステップ毎の小さな誤差が累積すると、全体のバランスエラーが大きくなることがあります。

自動増分制御（Automatic Time Stepping）は、解析が発散しそうな急激な変化（例えば、接触の発生や塑性化の開始）の直前に、増分サイズ（時間ステップや荷重ステップ）を自動的に小さくします。これにより、各ステップでのニュートン・ラフソン法の収束性が向上し、力の残差やエネルギーの誤差が小さく抑えられるため、結果として全体のエネルギーバランス誤差の蓄積を減らせます。Ansysでは「Preditor」設定もこれに影響します。

大きく変わります。陰解法（Abaqus/Standard, Ansys Mechanical）では、前述の通り反復計算の誤差累積が主な原因です。一方、陽解法（Abaqus/Explicit, LS-DYNA）では、時間積分に伴う数値的なエネルギー散逸（人工粘性）や発散（Hourglassエネルギー）が主要な誤差源です。LS-DYNAの結果では、Hourglassエネルギーが内部エネルギーの10%を超えないことが一つの目安とされます。陽解法では「エネルギーバランス」そのものが結果の信頼性を判断する最重要指標の一つです。

最低限、以下の5点をチェックリストとして確認します。1) 材料モデル（特に塑性データ）が適切か。降伏応力後の硬化則が定義されているか。2) 接触設定：接触剛性（Penalty Stiffness）が適切か。Abaqusでは「硬い」接触が推奨されますが、過大だと収束不良の原因になります。3) メッシュ品質：特に塑性変形が集中する領域のアスペクト比が10以下か。4) 境界条件：剛体運動が完全に拘束されているか。5) 増分制御：最初から自動ステップを有効にしておく。

まず、エラーが発生した「時系列」を特定します。Abaqusの.staファイルやAnsysのソルバー出力を確認し、エラーが急激に増加した荷重ステップまたは時間ステップを探します。次に、そのステップで何が起きたかを調べます。接触状態の急変、特定要素の大変形・破壊、材料の降伏開始などが典型的です。その領域のメッシュを細かくし、接触定義を見直して再解析します。また、塑性データにバイリニアモデルを使っている場合、多線形モデルに変更することで材料のエネルギー吸収挙動がより正確になり、改善することがあります。

エラーが小さくても、局所的に非物理的な現象が起きていないかの確認が必要です。具体的には、1) 反力の総和が印加荷重と力のつり合い（通常、誤差1%以内）を満たしているか。2) ひずみエネルギー密度の分布が、応力集中部と対応しているか。3) 塑性変形領域が、実験や経験則とおおむね一致しているか。例えば、丸棒のネック部の塑性変形は対称的に発生するはずです。これらの物理的妥当性と、数値的なエネルギー誤差を総合的に判断します。

はい、明確な違いがあります。Abaqus/Standardでは、.msgファイルや.staファイルに「THE ANALYSIS HAS COMPLETED」の後に「ENERGY BALANCE CHECK」という形で、外部仕事、内部エネルギー、バランス誤差率が直接出力されます。デフォルトの警告閾値は明確に定義されていませんが、実務では5%が目安です。一方、Ansys Mechanicalでは、ソルバー出力（Solution Information）の「Force Convergence」や「Energy Error」グラフをモニターし、また「Results」タブで「Energy」を選択して各種エネルギーをプロットし、手動で計算する必要があります。Ansys APDLでは「PRENERGY」コマンドで出力可能です。

COMSOLでは「保存則」の観点から各物理場のバランスを厳密に評価します。「研究」の設定で「計算中に保存則をチェック」を有効にできます。構造力学単体でも、エネルギー残差を評価できますが、より特徴的なのは、熱伝導と構造（熱応力）を連成させた場合などです。この時、熱エネルギーとひずみエネルギーの変換（熱弾性効果）が適切にモデル化されているか、系全体のエネルギー保存が成り立っているかを、各物理インターフェースの「フラックス/ソース項」のバランスから詳細に追跡できます。デフォルトの相対許容誤差は通常0.001に設定されています。

基本原理は同じですが、情報と自動化されたトラブルシューティング機能が少ないため、ユーザーの知識がより求められます。例えば、CalculiX（Abaqusと類似の入力形式）では、.datファイルに「total energy of the system」として内部エネルギーと運動エネルギーが出力されますが、外部仕事は直接出力されないため、反力と変位から自分で計算する必要があります。Code_Asterでは、`.resu`ファイルに「BILAN_ENERGIE」というコマンドで詳細なエネルギー内訳を出力させることができます。収束が悪い場合、ステップサイズや接触剛性を調整するパラメータを手動で細かく設定する必要があり、その最適値を見つける試行錯誤が多くなりがちです。

第一原因は「接触の貫通」です。接触が正しく定義されていないため、部品同士が重なり合い、見かけ上、外力が物体に仕事をしていない（または負の仕事をしている）ように計算され、エネルギー計算が破綻します。Ansysの場合、修正方法は以下の通りです。1) 接触設定を「Frictional」から「Frictionless」に一時的に簡素化。2) 「Normal Stiffness」を「Manual」にし、値（例えば、1.0E+05）を小さくしてから徐々に大きくする。3) 「Interface Treatment」で「Adjust to Touch」を適用。4) 初期接触状態を「プログラム制御」から手動で確認・修正する。

想定内ですが、その大きさに注意が必要です。材料が破壊（要素削除）すると、その要素が担っていたひずみエネルギーが突然システムから消失します。この消失エネルギーが外部仕事から差し引かれず、エラーとして計上されます。Abaqusの延性破壊モデルなどでは、これがある程度発生します。問題は、破壊が非現実的に早く、広範囲で起きていないかです。破壊ひずみの値（例えば0.3）が材料データと合っているか、メッシュ依存性の影響（細かいメッシュほど早く破壊する）を軽減するために「Mesh-Size Dependent」設定を確認してください。エラーが30%を超えるようなら、破壊モデルかメッシュに問題がある可能性が高いです。

最終手段として、以下の2つの方法がありますが、物理的正当性を損なわないよう注意が必要です。1) ソルバー設定の厳密化：Abaqusでは、ステップ設定で「Time Incrementation」の「Maximum number of increments」を大幅に増やし（例えば100,000）、「Increment size」の初期値を小さく（例えば1e-05）する。Ansysでは「Solver Controls」で「Weak Springs」をOffにし、「Large Deflection」をOnにする。2) モデルの物理的簡略化：局所的な大変形部分に剛体要素を適用する、または、問題の領域のみ超弾性材料などより柔軟な材料モデルに置き換えて変形を緩和する。これらは根本解決ではありませんが、エラー要因の切り分けや、やむを得ない場合の暫定対策として使われます。