平面応力問題 — トラブルシューティングガイド
より充実した内容を plane-stress.html でご覧いただけます。
平面応力解析のトラブル
平面応力解析でよくあるトラブルを教えてください。
平面応力は基本的な解析だが、初心者がつまずきやすいポイントがある。
平面応力と平面ひずみを間違える
これが最も多いミスですか?
間違いなく最多だ。結果の影響:
- 応力 — ポアソン比 $\nu = 0.3$ で約10%の差
- 非圧縮に近い材料($\nu \to 0.5$、ゴムなど)では致命的な差 — 平面ひずみで $\nu \to 0.5$ にすると剛性が無限大に発散する
確認方法:
- Abaqus: 要素名の頭文字 CPS = Plane Stress, CPE = Plane Strain
- Ansys: KEYOPT(3) = 0: Plane Stress, 1: Axisymmetric, 2: Plane Strain
- Nastran: PSHELLのMID参照で平面応力/ひずみが暗黙に決まる
ゴムのような非圧縮材料で平面ひずみ要素を使うと発散するんですか?
$\nu = 0.5$ では体積変化がゼロ(非圧縮)の拘束が入る。平面ひずみの完全積分要素では拘束条件が過剰になり、体積ロッキングが起きる。低減積分要素やハイブリッド要素(Abaqusの CPE4H)を使う必要がある。
応力集中値が理論と合わない
Kirsch問題(円孔のある板の一軸引張)で $K_t = 3.0$ になりません。
原因チェック:
1. メッシュが粗い — 孔の周囲に最低16要素(二次要素)。1次要素なら32以上
2. 板が有限サイズ — Kirschの理論解は無限板が前提。$d/W < 0.1$ でないと一致しない
3. 応力の読み取り位置 — 最大応力は孔の赤道上($\theta = 90°$)。要素の積分点から外挿した値か、節点平均値かで結果が変わる
4. 要素タイプ — CSTでは応力勾配を全く表現できない
応力の読み取り方法で結果が変わるんですか。
変わる。FEMの応力は積分点で最も正確だ。節点平均値(各要素の値を平均したもの)は滑らかになるが、応力集中のピーク値を過小評価する。非平均化(unaveraged)の応力コンターを見て、要素間の不連続がないかチェックすること。
変位が理論と合わない
一様引張の板で変位が $\delta = PL/AE$ と合いません。
確認項目:
- 板厚の設定 — デフォルトの1.0のままになっていないか。断面積 $A = t \times W$ が正しくないと変位が変わる
- 単位系 — $E$ の単位と荷重・寸法の単位が整合しているか
- ポアソン効果 — 一軸引張でも横収縮($\varepsilon_y = -\nu \varepsilon_x$)が生じるため、$y$ 方向の拘束状態が影響
要素が歪んで精度が落ちる
自動メッシュで歪んだ要素が出ました。どの程度まで許容できますか?
2次元要素の品質指標:
| 指標 | 理想値 | 許容範囲 | 影響 |
|---|---|---|---|
| アスペクト比 | 1.0 | < 5.0 | 精度低下 |
| スキューネス | 0° | < 45° | 精度低下 |
| ヤコビアン比 | 1.0 | > 0.3 | 負で要素退化 |
| ワーピング | 0 | < 0.1 | 精度低下(3Dで問題) |
ヤコビアンが負になると何が起きますか?
要素が裏返っている(節点の順序が逆)ことを意味する。計算は走るかもしれないが結果は完全に間違い。自動メッシュ後に必ず品質チェックを行うこと。
まとめ
平面応力のトラブル対処、整理します。
平面応力は簡単に見えて、基本的なミスが結果に直結する分野ですね。
FEMの全ての解析は基本の上に成り立っている。2次元平面応力で基本を徹底することが、複雑な3次元解析で間違いを防ぐ最良のトレーニングだ。
平面応力での厚み縮小の見落とし
平面応力ではポアソン効果によりεz=−ν(εx+εy)/(1−ν)が生じ、板は面外に変形する。FEMソルバーはこれを内部で自動計算するが、後処理で厚み変化を忘れて接触ギャップを評価すると過小なギャップ量を見込む誤りが起きる。プレス成形シミュレーションで2010年代にニデック社の解析部門が板厚変化の後処理マクロを整備し、設計部門との情報連携を標準化した。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——平面応力問題の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告