キルヒホッフ・ラブ薄肉シェル理論 — トラブルシューティングガイド
この記事は統合版に移行しました
より充実した内容を shell-kirchhoff-love.html でご覧いただけます。
より充実した内容を shell-kirchhoff-love.html でご覧いただけます。
キルヒホッフ・ラブ薄肉シェル理論 — トラブルシューティングガイド
K-Lシェルのトラブル
K-Lシェルに関連するトラブルを教えてください。
直接K-Lシェル要素を使うことは稀なので、ミンドリン系シェルを薄肉で使う際の注意点として整理する。
極薄シェルでのせん断ロッキング
$R/t > 1000$ の極薄シェルで精度が出ません。
ミンドリン系シェル要素のせん断ロッキング対策(MITC法等)が不十分なケースがある。極薄シェルでは:
- 二次要素(S8R, SHELL281)に切り替え — ロッキングに頑健
- STRI3/STRI65(K-L系要素)を検討 — せん断自由度がないのでロッキングの問題がない
- メッシュを細かく — ロッキングの影響を軽減
膜ロッキング
「膜ロッキング」って何ですか?
曲面シェルで曲げ変形に寄生的な膜ひずみが発生する現象。せん断ロッキングとは別の問題で、曲率のあるシェル特有だ。
例えば円筒シェルの純曲げ問題で、低次の要素が曲率効果を正しく表現できず、膜力が人為的に生じる。
対策:
- 二次要素を使う
- メッシュを細かく
- 曲面の曲率に対して十分な要素数
シェル理論の限界を超えた使用
厚いシェルにK-L理論を使ってしまった場合は?
$R/t < 10$ の厚肉シェルではK-L理論(およびミンドリン理論)の仮定が崩れる。板厚方向の応力が無視できなくなるため、ソリッド要素を使うべきだ。
確認方法:
- FEMの $\sigma_z$(板厚方向応力)が $\sigma_x, \sigma_y$ に比べて5%以上あるか
- ソリッド要素の結果と比較してずれがないか
まとめ
K-Lシェルのトラブル対処、整理します。
シェル理論の適用限界を知っていることが、正しい要素選択につながるんですね。
その通り。「どの理論が使えるか」を判断できることが、FEMエンジニアの根本的な能力だ。
Coffee Break よもやま話
KLシェルのC1不適合による誤差
C1連続性を満たさないKLシェル要素(例:通常のQ4シェル)を薄板に使うと「せん断ロッキング」が発生し、変位が実解の1/10以下になることがある。Timoshenko(1959年)の解析解との比較検証でこの問題は古くから知られており、板厚/スパン比が1/100以下になると急激に精度が低下する。診断法はメッシュ収束テスト:要素数を4倍にしても変位がほとんど変わらない場合はロッキングを疑う。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——キルヒホッフ・ラブ薄肉シェル理論の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
この記事の評価
ご回答ありがとうございます!
参考に
なった
なった
もっと
詳しく
詳しく
誤りを
報告
報告