材料モデルの不安定性（Drucker安定性違反）

良い質問だ。Drucker安定性は、材料が安定した塑性挙動を示すための公準で、増分塑性理論の根幹をなす。簡単に言えば、外部から材料に加えられる「仕事」が常に正でなければならないという条件だ。数学的には、応力増分

物理的に非現実的な解を導き、数値計算が破綻するからだ。例えば、降伏面が凹んだ形状（非凸性）を持つ材料モデルを使うと、この条件に違反しやすい。具体的には、ある応力状態から別の応力状態に移る際に、材料が外部に対して正の仕事を「してしまう」、つまりエネルギーを生み出すような挙動を許してしまう。FEMのソルバーはこれを検知すると、剛性マトリクスが正定値でなくなり、Newton-Raphson反復が発散して「解が収束しない」エラーを引き起こす。Abaqus/Explicitであれば、安定時間増分が極端に小さくなり、事実上計算が進まなくなる。

古典的な例は、Mohr-CoulombやDrucker-Pragerといった摩擦角を考慮したモデルだ。特に、三軸引張り状態や低い平均応力下で、降伏面に「くぼみ」が生じる。もう一つの例は、異方性塑性モデル、例えばHill'48モデルでR値（異方性パラメータ）の組み合わせによっては非凸領域が生じ得る。実務では、地盤のせん断破壊解析でMohr-Coulombを使う際、パラメータ設定（特に摩擦角φと粘着力cの比）に注意が必要で、φが大きすぎるとDrucker安定性違反のリスクが高まる。

商用ソルバーは様々だ。Abaqus/Standardの多くの塑性モデルでは、材料接線剛性マトリクス（Jacobian）を構成する際に、安定性を暗黙的に仮定している。違反が起こると、接線剛性マトリクスが正定値でなくなり、直接法ソルバー（例えば、デフォルトの疎行列ソルバー）でピボットエラーが発生したり、反復法ソルバーで収束しなくなる。一方、Abaqus/Explicitでは、材料応答から算出される音速が虚数になり、安定時間増分Δtが計算不能になることで顕在化する。

その通りだ。ユーザーサブルーチンでは、応力更新アルゴリズムと一貫性接線剛性（Consistent Tangent Modulus）の導出が全ての鍵を握る。例えば、Drucker-Pragerモデルを実装する際、非関連流れ則（塑性ポテンシャル≠降伏関数）を使うと、接線剛性マトリクスが非対称になり、Drucker安定性を満たさない可能性が高まる。これを回避するには、関連流れ則を使うか、数値的に「正則化」を行う必要がある。具体的には、接線剛性マトリクスの負の固有値をクリッピングして小さな正の値に置き換えるなどの手法が論文で紹介されている。

一つのアプローチは、凹んだ部分を微小な円弧で置き換える「角の平滑化」だ。Mohr-Coulombの六角形の降伏面を、Drucker-Pragerの円錐面で近似するのもその一例だが、精度が落ちる。より洗練された方法では、降伏関数Fを修正せず、塑性変形の際の接線剛性

まず、摩擦角を持つモデル（Mohr-Coulomb, Drucker-Prager）を使う場合だ。地盤材料では摩擦角φが30度を超えることが多いが、CAEでは安定的に計算するために、φを実際の値より低く設定することがある。例えば、実際の砂がφ=40度だとしても、初期検討ではφ=30度で解析し、挙動を確認する。また、Drucker-Pragerモデルでは、圧子応力pとせん断応力qの空間で、降伏面の傾き（摩擦角に相当）βと、三軸引張り強さと三軸圧縮強さの比Kを適切に設定する。Abaqusマニュアルでは、安定性のためには

その懸念はもっともだ。だからこそ、パラメータを安易にいじるのではなく、以下のステップを踏む。1) 単純な単要素試験（Uniaxial, Pure Shear）で材料モデル自体の挙動を検証する。2) 実際の構造解析で不安定性エラーが出た場合、そのエラーが発生する「応力パス」を特定する。多くの場合、低い平均応力下でのせん断負荷が原因だ。3) その問題の応力領域だけをカバーするために、複数の材料モデルを領域ごとに割り当てるか、より高度なキャップモデル（Cambridge Clayモデルなど）への切り替えを検討する。ANSYSのDPモデルでは「Tensile Strength」パラメータを0より大きく設定することで、引張り側の角を丸め、安定化を図れる。

ある。代表的なのは「軸対称」または「平面ひずみ」状態の単一要素に、様々な平均応力とせん断応力の組み合わせ（異なる応力パス）を強制変位で与えるテストだ。例えば、一定の圧子応力を保ちながらせん断変形を増加させる試験（排水せん断試験のシミュレーション）を行う。この時、反力として計算される応力-ひずみ曲線が単調増加しているか、ソルバーの収束履歴に振動がないかを確認する。COMSOL Multiphysicsの「材料試験」機能や、Abaqusの「Material Calibration」モジュールを使えば、この種の検証を比較的容易に行える。

大きな違いがある。まずAbaqusは、組み込みのDrucker-Pragerモデルにおいて、ユーザーが「Smooth」オプションを選択できる。これは降伏面の角を自動的に平滑化し、非凸性を低減する機能だ。また、Abaqus/Explicitでは、材料から計算される体積音速とせん断音速を監視し、虚数値になるとエラーメッセージを出力する。一方、Ansys Mechanical APDLのDPモデルは、より古典的な実装に近く、引張り切断面（Tension Cut-off）を明示的に定義しないと、低圧力域で不安定になりやすい。その代わり、ユーザーがパラメータの制約を理解して設定する必要がある。

大きく異なる。Code_Aster（Salome-Mecaプラットフォーム）はドキュメントが非常に詳細で、理論マニュアル（U）に材料モデルの安定性条件について言及されていることがある。しかし、計算エンジン自体は安定性違反に対してより「無防備」で、計算が突然停止したり、非物理的な変形が発生したりする可能性が高い。CalculiXのDrucker-Prager実装は比較的シンプルで、安定化のための特別なオプションは少ない。これらのソフトを使う場合は、ユーザー自身が材料パラメータの選択により一層注意を払い、単純なケースから慎重に検証を積み重ねる必要がある。商用ソフトの「ブラックボックス的な安定化処理」に頼れないからだ。

決定的な証拠は、エラーが発生する「応力状態」を調べることだ。Abaqusでは.msgファイルや.datファイルを確認し、収束に失敗した増分ステップでの、問題の要素の応力成分を出力させる。特に、平均応力（Pressure）が引張り側（正）またはゼロに近く、せん断応力が働いている場合、疑わしい。また、「負の固有値」は全体剛性マトリクスの話だが、その原因が特定の要素の材料剛性（接線剛性）の異常にある場合が多い。対策の第一歩は、その要素の材料パラメータ、特に摩擦角や粘着力を見直し、より保守的な値に変更して再計算してみることだ。

いくつかの工学的対策がある。1) メッシュを粗くする：不安定性はひずみの局所化に伴って発生することが多い。メッシュを粗くすることでひずみ集中が緩和され、数値的に安定することがある。ただし、結果の精度は落ちる。2) 非関連流れ則から関連流れ則に変更する：塑性ポテンシャルを降伏関数と一致させると、多くの場合で安定性が向上する。ただし、体積膨張率が理論値から変わってしまう。3) ソルバー設定を変更する：Abaqus/Standardで「非対称マトリクスソルバー」を使用する。非関連流れ則では接線剛性が非対称になるため、これがデフォルトの対称ソルバーと相容れず、不安定性を助長している可能性がある。

強く関連する。Explicit法の安定条件は