材料モデルと要素タイプの不適合

その通りです。これはソフトウェアの単なる「嫌がらせ」ではなく、数値解析の理論的基盤に起因します。例えば、平面応力要素（CPS4など）に、体積ひずみに依存する塑性モデル（例えば、Gurson-Tvergaard-Needlemanモデル）を割り当てようとすると、理論的に矛盾が生じます。平面応力状態では厚さ方向の応力はゼロと仮定しますが、GTNモデルは空洞の体積変化を計算するために3次元の応力・ひずみ状態を必要とするからです。

よくある組み合わせの不適合例を挙げましょう。Abaqusの「超弾性材料モデル」（Mooney-RivlinやOgdenモデル）を、シェル要素（S4Rなど）に直接割り当てることはできません。超弾性モデルは非圧縮性を仮定することが多く、その定式化には3次元の変形勾配テンソルが必要です。シェル要素は厚さ方向の積分点が少なく、この情報を完全に表現できません。代わりに、連続体シェル要素（SC8Rなど）やソリッド要素を使う必要があります。

材料の応力-ひずみ関係の計算は、各積分点で独立に行われます。例えば、完全積分要素（AbaqusのC3D8）は要素内に8つの積分点を持ち、減積分要素（C3D8R）は1つしか持ちません。クリープや塑性のような履歴依存性の高い材料モデルでは、積分点ごとに内部状態変数（塑性ひずみ、ダメージパラメータなど）を保存・更新します。減積分要素ではこの情報が1点に集約されるため、局所的なひずみ集中やバンディング現象を捉えきれず、結果が不正確になったり、発散したりするリスクがあります。Ansysのマニュアルでも、高次の材料非線形解析には完全積分要素の使用を推奨しています。

ほとんどの商用CAEソフト（Abaqus/Standard, Ansys Mechanical, MSC Marc）は、入力ファイル（.inp, .dat）の読み込みと前処理の段階でチェックします。要素ライブラリと材料モデルライブラリのマトリクスを照合し、許可されていない組み合わせを早期に検出してエラーメッセージを出力します。例えば、Abaqusでは要素タイプが「CPS4」と定義されているセクションに「超弾性」材料を割り当てると、.datファイル生成時に「THE MATERIAL *MATERIAL NAME* IS NOT AVAILABLE FOR PLANE STRESS ELEMENTS」といったエラーで停止します。

非常に重要な点です。非圧縮性材料（ゴム、生体組織、金属塑性）を扱う場合、体積ひずみが非常に小さくなる（またはゼロになる）ため、変位ベースの標準要素では体積ロッキングという数値的不安定性が発生します。これを回避するために、圧力自由度を追加した混合ベースの定式化（u-p定式化）を用いた要素が開発されています。例えば、AbaqusのC3D8H（HはHybrid）、COMSOLの「混合要素」、Ansysの「Pure Displacement」ではなく「Mixed u-P」オプションがこれに該当します。これらの要素は、材料モデル側で非圧縮性パラメータを適切に定義する必要があります。

直接的には「不適合エラー」にはなりませんが、間接的に重大な影響を与えます。最も典型的なのは、ポアソン比νが0.5に近い（非圧縮性）材料です。線形弾性の範囲では、体積弾性率Kとせん断弾性率Gの関係は

はい、以下の5項目を確認することを習慣づけてください。

エラーで止まらない場合、よりたちの悪い「間違った結果」が出力されることがあります。具体例を挙げます。金属板のプレス成形シミュレーションで、シェル要素に等方性硬化則のみの材料モデルを使ったとします。実際の材料にはバウシュンガー効果（包応力の影響）があるため、反転曲げ時の応力予測が大きく外れます。正しくは、キネマティック硬化または混合硬化モデルを用いるべきです。この場合、ソフトはエラーを出しませんが、金型のスプリングバック量や成形限界の予測が不正確になり、実物試作で初めて問題が発覚するという最悪のシナリオになります。

はい、可能です。例えば、Abaqusの複合シェル要素（例えばS4）では、「Composite Layup」を定義し、各積層（Ply）ごとに異なる材料を割り当てられます。ただし、各材料モデルが「シェル要素で使用可能」である必要があります。ここで注意点が2つ。第一に、面内特性（繊維方向の弾性率など）は定義できても、厚さ方向のせん断特性（Interlaminar Shear）を正確に表現するには「連続体シェル要素」や「ソリッド要素」を用いたモデリングを検討する必要があります。第二に、各層の材料が破壊基準（Hashin, Puckなど）を持つ場合、その破壊モデルが複合シェル要素と互換性があるか、マニュアルで確認が必要です。

大きな違いがあります。まず、**Abaqus**は非常に厳格で、理論的に疑わしい組み合わせはほとんどエラーまたは警告として事前に排除します。要素名に機能がエンコードされています（例: C3D8Hの「H」はHybrid=混合要素）。非圧縮性材料には「H」要素を使え、と明示されます。一方、**Ansys Mechanical**はやや融通的で、例えば超弾性材料をシェル要素に割り当てようとすると、「この材料モデルはソリッド要素を推奨します」という警告は出ますが、強制的に実行させるオプションがある場合もあります。ただし、結果の信頼性は保証されません。Ansysでは混合要素は「Element Technology」の設定で「Mixed u-P」を選択することで有効になります。

根本的な哲学が異なります。**COMSOL**は「物理場」を第一に考えます。「固体力学」物理場インターフェースを選択した時点で、使用可能な材料モデル（線形弾性、超弾性、塑性など）が決まり、さらにその中で「平面応力」「平面ひずみ」「3D」などの「空間次元」を選択します。その後、メッシュを切ると、COMSOLが背後でその物理場と空間次元に最適な要素タイプ（多くの場合、ラグランジュ要素）を自動的に割り当てます。ユーザーが要素次数（2次など）を選択できますが、材料モデルと要素の根本的な不適合が起きるリスクは低く設計されています。ただし、高度な非線形材料を扱う際の要素技術（減積分 vs 完全積分）の選択肢は、Abaqusほど細かくはないかもしれません。

商用ソフトよりユーザーの責任が重くなります。例えば、**CalculiX**（Abaqusと互換性の高い入力形式）では、非対応の組み合わせを指定しても、実行時エラーやセグメンテーションフォルトで異常終了する可能性が高く、親切なエラーメッセージは期待できません。**Code_Aster**はドキュメントが非常に詳細で、各要素（例えばMECA3D8）と各材料モデル（例えばELAS_HYPER）の互換性表が提供されています。しかし、それを自分で調べ、入力ファイル（.comm）に正しく記述する必要があります。オープンソースでは「ソフトが守ってくれる」という安全網が薄いため、理論的な理解がより重要になります。

それはおそらく、材料「Steel」の塑性定義に問題があります。Abaqusで塑性を定義する際、「Plastic」プロパティで応力-塑性ひずみデータを入力しますが、このデータが「単軸引張試験」に基づく真応力-真ひずみデータであることが前提です。平面ひずみ状態では、変形が多軸的になるため、この単軸データをそのまま用いることが理論的に適切でないとソフトが判断した可能性があります。対策としては、まず材料が「Elastic」プロパティのみで定義されていないか確認し、塑性データを入力している場合は、それが適切な形式かマニュアルで再確認してください。また、平面ひずみ用の特殊な材料モデルが必要なわけではなく、標準的な弾塑性モデルはCPE4で通常使用可能です。入力データのフォーマット誤りが疑われます。

ここがAnsysの複雑な点です。SOLID185は3次元ソリッド要素ですが、複合材料の層別破壊（Layer-wise Failure）をシミュレーションする機能は、**SOLID185自体には組み込まれていません**。Ansysでは、複合材料の積層と破壊解析には、専用の「層別複合材料要素」である**SOLID186 LAYERED** または **SHELL181** を使用する必要があります。Chang-Chang破壊基準は、これらの要素と連携する形で実装されています。対策：要素タイプをSOLID186 LAYEREDに変更し、層定義と材料方向を設定し直してください。要素タイプの変更は、メッシュのクリアと再生成が必要になる場合が多いです。

直接の原因は要素の形状歪みですが、背後に材料と要素の不適合が潜んでいることがあります。典型的なのは、**大変形・非圧縮性材料（ゴム）に標準の変位ベースソリッド要素（例：AbaqusのC3D8R）を使った場合**です。要素がせん断変形に耐えられず、アワーグラッシングモードが発生して異常に歪み、このエラーを引き起こします。対策は以下の順で試します。