負の弾性係数エラー
概要
先生！ 今日は負の弾性係数エラーの話なんですよね？ どんなものなんですか？
理論と物理
弾性係数の物理的意味
「負の弾性係数エラー」と出ました。弾性係数って、材料がどれだけ変形しにくいかを表す値ですよね？ なぜ負の値が物理的にありえないんですか？
その通りです。弾性係数、特にヤング率Eは、引張応力と引張ひずみの比例定数で、
では、ポアソン比も負になることはあるんですか？ 横ひずみと縦ひずみの比ですよね。
良い質問です。通常の等方性材料のポアソン比νは、熱力学的制約から-1 < ν < 0.5の範囲にあります。ほとんどの金属は0.3前後です。しかし、特殊なメタマテリアルでは「負のポアソン比」を示すものも研究されています。それでもヤング率とせん断弾性係数Gは正である必要があり、
体積弾性率Kはどうですか？ 圧力をかけた時に体積が増える材料ってあるんですか？
通常の安定した物質では、圧縮すると体積が減るので、体積弾性率Kも正です。KはEとνから
数値解法と実装
ソルバーによる入力チェック
エラーはいつ検出されるんですか？ 前処理で弾性係数を入力する時？ それともソルバーが計算を始めてから？
ほとんどの商用ソルバーは、入力ファイルの読み込み段階、または剛性マトリクスを組み立てる前の材料モデル初期化段階でチェックします。例えば、Abaqus/Standardは`.inp`ファイルを読み込む際、材料プロパティセクションをスキャンし、Eやνの値が許容範囲外ならジョブ送信前にエラーを出力します。
剛性マトリクスにどう影響するんですか？ 負の値だと具体的に何がまずいのですか？
線形弾性の要素剛性マトリクス
非線形材料や超弾性モデルでは、初期接線剛性が負になることもありますよね？ それとは違うんですか？
全く別の問題です。非線形解析での負の接線剛性は、材料の軟化や座屈などの物理現象を表し、ニュートン・ラフソン法の収束を難しくしますが、入力パラメータ自体は正です。一方、「負の弾性係数エラー」は、線形材料モデルの入力パラメータ（E, ν, Gなど）そのものが負または許容範囲外であるという、根本的な入力ミスを指しています。
実践ガイド
入力ミスの特定と修正
実際にエラーが出た時、まずどこを確認すべきですか？ 入力ファイルが大きいと探すのが大変です。
第一に、エラーメッセージを精読することです。例えば、「Negative or zero elasticity modulus detected for material 'STEEL_1'」とあれば、材料名が特定されています。次に、その材料のプロパティ入力箇所へジャンプします。多くのプリプロセッサでは、材料テーブルの数値に指数表記（2.05e+11）や単位系の不一致（MPaとGPaの混在）による間違いが頻発します。鋼のEを205 MPa（実際は205 GPa）と入力していませんか？
単位系の不一致はよく聞きます。具体的にどう混同するんですか？
典型的なのは、幾何（長さ）をmmでモデリングし、力をNで入力する場合です。この単位系では、応力はN/mm² (=MPa)になります。ここで、鋼のヤング率をカタログ値の210 GPa (=210,000 MPa)のつもりで「210000」と入力すればOKです。しかし、「210」とだけ入力すると、ソフトウェアはそれを210 MPaと解釈し、一見正の値ですが、実際より3桁小さい値になります。一方、「2.1e-4」のような誤入力で負の指数を使うと、非常に小さな正の値になったり、場合によっては「0」と認識されて別のエラーを引き起こすこともあります。
材料プロパティをExcelなどからコピペする時は注意が必要そうですね。他にありがちな入力ミスは？
はい、特にCSVデータの区切り文字や負号の認識ミスです。例えば、温度依存材料で「-20, 2.06e11」という行がある時、ソフトウェアがカンマを小数点と誤認し、「-20.2」を温度、「06e11」を弾性係数と読み、後者が数値として解釈不能になるケースがあります。また、材料定数を計算で定義する場合、
ソフトウェア比較
各ソフトウェアのエラーメッセージと挙動
このエラーは、Ansys、Abaqus、COMSOLなど、ソフトウェアによってメッセージやチェックの厳しさは違うんですか？
はい、異なります。Ansys Mechanical APDLでは、`MP`コマンドで弾性係数を定義する際、負の値を入力しても即時エラーにはならず、求解コマンド`SOLVE`を実行した際に「Small or negative pivot term encountered.」という、より一般的なエラーとして現れることがあります。一方、Abaqus/CAEでは、材料プロパティマネージャーで「200e-9」のような非常に小さい正の値を入力すると、「The modulus may not be less than 1.0E-12」という警告を出し、事実上0とみなします。
COMSOLはどうですか？ マルチフィジックスで材料を多用する印象があります。
COMSOL Multiphysicsでは、材料ノードで「線形弾性材料」を選択した時、ヤング率とポアソン比の入力フィールドが現れます。ここに負の値を入力して「計算」ボタンを押すと、大抵は警告やエラーが表示されます。興味深いのは、COMSOLが「実数」と「複素数」の両方の材料定数を扱える点です。波動伝搬の問題で複素弾性係数（損失を表現）を使うことがありますが、その実数部が負にならないかはユーザーが管理する必要があります。
無償のオープンソースソルバー、例えばCalculiXやCode_Asterではどうでしょう？
CalculiX (CCX) は入力ファイル`.inp`の`*ELASTIC`行を読み、値が正であるかチェックします。負の値があればエラーで停止します。Code_Aster (Salome-Meca環境) では、コマンドファイル`.comm`の`AFFE_MATERIAU`で定義しますが、その前段階のグラフィカルインターフェース（Salome）では基本的な範囲チェックが行われます。ただし、オープンソースはエラーメッセージが技術的で冗長な場合があり、「MATERIAL KEYWORD ELASTIC HAS ILLEGAL VALUE」などと出力され、具体的な材料名や行番号をログから探す必要があります。
トラブルシューティング
エラー解決のステップ
エラーメッセージに材料名が出力されなかった場合、どうやって問題の材料を特定すればいいですか？
系統的なアプローチが必要です。まず、使用している全ての材料のプロパティをリストアップします。多くのソフトウェアは、前処理段階で材料プロパティのサマリーレポートを出力する機能があります（例：Abaqusの`.dat`ファイル、Ansysの`.log`ファイル）。それを確認し、E, ν, G, Kの値をチェックします。もし明らかな負の値が見当たらなければ、温度依存や場の変数依存のプロパティ定義を疑います。全ての温度テーブル、全データポイントを確認してください。
温度依存プロパティで、ある温度点だけ間違えることもあるんですか？
よくあります。例えば、SUS304ステンレス鋼の高温特性を入力する時、参考文献の表から値を転記する際、行を1つずらしてしまい、800°Cのデータを誤って-800°Cの列から取ってくるかもしれません。あるいは、表のフォーマットが「温度(°C), E(GPa)」なのに、「E(GPa), 温度(°C)」の順で入力してしまうと、弾性係数の列に温度値（例えば20）が入り、それは正ですが非現実的に小さな値となり、別の警告を誘発します。
問題の材料を特定して修正した後、再計算する前に他に確認することは？
絶対にやってほしいことが二つあります。第一に、修正した材料が実際に全ての部品/要素に割り当てられているか確認すること。古い材料定義が割り当てられたままのパートが残っている可能性があります。第二に、材料定数間の整合性チェックです。例えば、Eとνを入力した場合、せん断弾性係数Gと体積弾性率Kが計算されますが、これらが常に正になるか、
結局、このエラーは単純な入力ミスが原因なんですね。予防するためのベストプラクティスはありますか？
その通り、ほぼ100%が人的ミスです。予防策としては、(1) 材料ライブラリを作成し、検証済みの値を再利用する。(2) 入力時は単位系を徹底し、数値の後ろにコメントで単位を記入する（例：`2.05e11 ! [Pa]`）。(3) 重要な解析の前には、材料プロパティだけを読み込む簡単なテストモデル（単純な片持ち梁など）を作成し、常識的な変位が得られるか確認する。(4) チームで作業する場合は、材料入力用の標準化されたテンプレートやチェックリストを使用する。これらを徹底すれば、この基本的なエラーに悩まされる時間を大幅に減らせます。
関連トピック
なった
詳しく
報告