応力テンソル — CAE用語解説

良い質問だ。決定的な違いは「座標変換に対する振る舞い」にある。単なる行列は座標を変えれば成分がバラバラに変わるが、応力テンソルの成分は特定のルールに従って変換される。例えば、引張試験片を45度傾けて応力を測ると、垂直応力だけでなくせん断応力も現れる。この関係は、

失われているのではなく、特定の破壊基準に基づいて「要約」している。ミーゼス応力（相当応力）は、

応力テンソルの固有値・固有ベクトルを計算している。具体的には、特性方程式

本質的な理由は、変位の近似関数（形状関数）と応力の関係にある。変位uは形状関数Nで補間される：

ごまかしではなく、数値結果をより現実に近づけるための「後処理」だ。代表的な手法は「節点平均化」。各積分点で求めた要素ごとの応力を、共有する節点に持ち寄り、平均値を計算する。Abaqusでは「Averaging Threshold」を設定でき、デフォルトは75%だ。これは、隣接要素間の応力差が75%以内なら平均化し、超えるなら平均化せずに表示する。鋭い応力集中部では平均化をオフにしないと、ピーク応力が見逃される危険がある。

ソルバーが変位uを求めた「後」に、ポスト処理段階で計算される。具体的な流れは：1. ソルバーが全体剛性行列Kと荷重ベクトルFから変位uを求解。2. 各要素ループ内で、要素節点変位u_eから要素ひずみε_eを計算（Bマトリックス使用：ε_e = B u_e）。3. 構成則マトリックスDを使って要素応力σ_e = D ε_eを計算。この計算は「要素定式化」を担当するモジュールが行い、Ansysでは「Element Library」、Abaqusでは「Element Routine」が該当する。結果は積分点（ガウス点）レベルでまず得られる。

材料の破壊モードと設計規格に従う。大原則は：
1. 延性材料（鋼、アルミ）の静強度：ミーゼス相当応力。自動車車体（SPCC）や機械架台（SS400）など。許容値は降伏強さ（例えば245MPa）に安全率（1.5〜2）をかける。
2. 脆性材料（鋳鉄、ガラス）：最大主応力（第一主応力）。鋳鉄製ブラケット（FC250）など。引張強さを基準にする。
3. 疲労強度：せん断応力や主応力振幅が重要。JIS B 8265（圧力容器構造）などでは修正グッドマン線図と組み合わせる。
まずは使用材料のデータシートと関連規格（JIS, ISO, ASME）を確認せよ。

理論的には「応力特異点」がある限り、メッシュを細かくすれば応力は発散する。正解は「設計上有意味な値」を抽出することだ。実務的な手順：
1. 幾何学的な切り欠き根部の曲率半径Rを測定（例：R=1mm）。
2. そのRに対して、要素サイズを1/4, 1/8, 1/16と細かくし、最大応力の変化を追う。
3. 変化が例えば5%以内に収まったメッシュサイズを採用する。あるいは、構造応力法を用いる。溶接部評価のIIW（国際溶接学会）勧告や、Autodesk Nastranの「Structural Stress」機能は、特異点近傍の応力を線形外挿して、メッシュ依存性を低減する手法を採用している。

すぐに「信じる」対象を決めるのではなく、不一致の原因を特定する。チェックリスト：
1. 境界条件：理論式は「単純支持」だが、FEMでは「面拘束」で剛体運動を抑制していないか？面拘束は局部剛性を上げ、応力を過小評価する。
2. 荷重の与え方：集中荷重と分布荷重の違い。
3. 断面二次モーメント：理論値で使った値と、FEMの3Dモデルの実際の値が一致しているか？
4. 応力評価位置：理論値は「中立軸から最も遠い点」だが、FEMで見ているのは節点平均化された値かもしれない。
まずは、梁や板など単純形状で、メッシュと境界条件を厳密に合わせた検証モデル（ベンチマーク）を作成し、ソルバー自体の信頼性を確認する作業が不可欠だ。

出力される成分自体は同じだが、デフォルトの座標系と出力制御に違いがある。
・Ansys：デフォルトは「全体直交座標系（Global Cartesian）」。X, Y, Z方向の垂直応力とせん断応力を出力。ユーザーが「局部座標系」を定義して、その方向の応力成分を出力できる機能が強力。例えばボルトの軸方向応力を出す時など。
・Abaqus：デフォルトも全体座標系だが、「出力変数リクエスト」で「S」と指定すると、全ての応力成分が出力される。「出力変数」のフィルタリングが細かく設定可能。
また、Ansys APDLでは「RSYS」コマンドで結果座標系を切り替えられるのが古参ユーザーに好まれる。

COMSOLの特徴は、異なる物理場で生じる応力を自然に足し合わせられる「多重物理量結合」にある。例えば、
1. 熱応力：温度場から計算された熱ひずみによる応力テンソルと、機械荷重による応力テンソルを自動的に線形重ね合わせる。
2. 圧電効果：電気場と応力場の結合。生成される応力テンソルは電気変数に依存する。
3. 独自の物理場インターフェース：ユーザーがPDEを定義すれば、それに対応する「応力」のような変数も定義できる。
他のソフトでは別々の解析結果を手動で合成する必要がある場合でも、COMSOLでは一つの連成問題として定式化・求解されるため、整合性が取りやすい。

計算アルゴリズムのコア部分は確立されているので、単純な線形静解析であれば正しく計算できる。CalculiX（Abaqusと互換性のある入力形式）は実績が豊富だ。ただし、注意点は：
1. 要素ライブラリ：商用ソフト（Ansys：200種以上）に比べて対応要素が限られる。複雑な複合材シェル要素や超弾性要素はないか貧弱。
2. 後処理と検証機能：応力の平滑化オプションや、局部座標系での出力、自動収束判定などの機能が商用ソフトより劣る。
3. 材料モデル：非線形、クリープ、疲労などの高度な材料モデルが不足している。
学習や単純部品の検討には十分だが、製品の設計判断を下す際には、商用ソフトでの検証か、実験との対照が必須となる。

主な原因は3つ：
1. 剛体運動（リジッドボディモーション）：拘束不足。最も多い。変位結果で物体全体が明らかに動いていないか確認。Ansysでは「弱いばね」オプションをオンにすると検出されやすくなる。
2. 材料定数の単位系不一致：ヤング率を200GPa (=200,000 MPa) とすべきところを200と入力し、長さ単位mmと組み合わせると、応力が1000倍違う。
3. 不適切な要素変形：例えば、ソリッド要素でアスペクト比が1000を超えるような極端に細長い要素があると、数値誤差で異常な応力が計算される。メッシュ品質をチェックせよ。

必ずしも粗さだけが原因ではない。これは「要素ごとの非連続な応力」が可視化されている状態だ。先述の通り、FEMで計算された応力は要素内部で連続でも、要素境界では隣接要素と値が一致しない。対策は：
1. 節点平均化の適用：ほとんどのソフトでデフォルトでON。Abaqusでは「Plot Contours on Deformed Shape」のオプションで「Averaged」を選択。
2. 平均化してもジャギーが残る場合は、メッシュ自体が不適切な可能性。応力勾配が急峻な領域（切り欠き根部）のメッシュが粗すぎる。その領域を特定し、メッシュを局所細分化する。
3. 異種要素の混在：ソリッド要素とシェル要素が接合されている場合、接合面で応力が不連続になるのは自然。その場合は別々に評価する。

即座に破壊を意味しない。 線形解析は「応力-ひずみ関係が常に線形」と仮定しているので、降伏を超えた後の応力再配分や局部塑性化を考慮できない。降伏強さを超える応力が生じた場合の正しい手順は：
1. 応力集中の要因を確認：鋭い角や微小な切り欠きなど、現実には存在しない幾何学的特異点による過大評価ではないか。
2. 弾塑性解析への移行：材料に塑性域のデータ（例えば二直線硬化則）を定義し、非線形静解析を実行する。Ansysなら「Multilinear Kinematic Hardening」などを使用。
3. 許容応力の再考：一部の規格（例えば建築）では、限界状態設計法を用い、一時的に降伏を超えても全体崩壊に至らないケースを許容する場合がある。
線形解析で降伏強度を超える結果は、「この領域はより精緻な評価が必要」という警告サインと捉えるべきだ。