ミンドリン・ライスナー板理論 — トラブルシューティングガイド
この記事は統合版に移行しました
より充実した内容を plate-mindlin.html でご覧いただけます。
より充実した内容を plate-mindlin.html でご覧いただけます。
ミンドリン・ライスナー板理論 — トラブルシューティングガイド
シェル要素のトラブル
シェル要素（ミンドリン板ベース）でよくあるトラブルを教えてください。
シェル要素は最も広く使われる要素だけに、トラブルも多様だ。
せん断ロッキング
薄板なのにたわみが理論値より小さいです。
完全積分のシェル要素（S4等）で稀にせん断ロッキングが起きることがある。対策：
アワーグラスモード
S4Rで変形がジグザグになります。
低減積分のシェル要素でアワーグラスモードが励起されている。
原因：
- 集中荷重
- メッシュが粗い
- 要素の1辺に荷重が偏る
対策：
- 荷重を分散（RBE3等で）
- メッシュを細かく
- Abaqusのホバーグラス剛性を増加
- S4（完全積分）に切り替え
法線方向の向きが不統一
応力結果が隣接要素でプラスとマイナスが逆転しています。
要素の法線方向が反転している。シェル要素は「上面」「下面」の区別があり、法線方向（正のz方向）が要素ごとに異なると、曲げ応力の符号が逆転する。
対策：
- プリプロセッサで全要素の法線方向を統一
- Abaqusの *NORMAL で法線を明示指定
- Nastranの PARAM,SNORM で法線の自動調整
これは見落とされやすいトラブルですね。
法線の不統一はモーメント図の不連続として現れる。モーメントの符号が要素間で反転していたら、法線方向を確認すべき。
板厚が変化する部位の応力
板厚が変わるところで応力がおかしいです。
シェル要素の板厚は各要素で定義される。板厚が急変する部位では：
- 膜力の不連続（同じ膜応力でも板厚が違うと膜力が不連続）
- 中立面のオフセット（板厚変化で中立面がずれる）
対策：
- 板厚変化部にオフセットを正しく設定
- 板厚変化を滑らかに遷移させる（テーパー）
- 必要に応じてソリッド要素に切り替え
まとめ
シェル要素のトラブル対処、整理します。
法線方向の不統一は盲点になりやすいですね。応力の符号が逆転するなんて…。
シェル要素は「表裏」があるのが2次元平面要素やソリッド要素と本質的に異なる点だ。この「表裏」を常に意識することがシェル要素を使いこなす鍵だ。
Coffee Break よもやま話
ドリリングDOF問題の対処
Mindlin板要素で面内回転（ドリリングDOF）を持たない要素を非平面メッシュに適用すると、節点が面外方向に「ぐらつく」数値的特異点が生じる面外特異性問題が起きる。ANSYS SHELL181のKEYOPT(3)=0はドリリングDOFを自動補剛するが、過剛な補剛で高アスペクト比要素の面内応力が5〜15%誤る事例が2008年のANSYS検証事例集に記載されている。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ミンドリン・ライスナー板理論の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
この記事の評価
ご回答ありがとうございます！
参考に
なった
なった
もっと
詳しく
詳しく
誤りを
報告
報告