キルヒホッフ板理論 — トラブルシューティングガイド
より充実した内容を plate-kirchhoff.html でご覧いただけます。
キルヒホッフ板のトラブル
キルヒホッフ板理論に関連するトラブルを教えてください。
直接的なFEM要素のトラブルではなく、理論の適用範囲を超えたときに問題が起きる。
せん断変形の無視が不適切
キルヒホッフ理論のたわみ公式で計算した結果と、FEM（ミンドリン要素）の結果が合いません。
板が厚すぎる（$b/t < 20$）。キルヒホッフ理論はせん断変形を無視するから、厚板ではFEMのほうが正しい（ミンドリン要素がせん断変形を含む）。
確認方法：ミンドリン板のFEMで $t$ を非常に薄くする（$b/t > 100$）。キルヒホッフの理論値に収束するはず。
支持条件の違い
理論値とFEMで支持条件が違うのでは？
キルヒホッフ板の「単純支持」は：
- $w = 0$（面外変位ゼロ）
- $M_n = 0$（法線方向モーメントゼロ）
FEMで $w = 0$ だけ拘束すると単純支持。$w = 0$ かつ $\theta = 0$ を拘束すると固定支持。混同しないこと。
集中荷重による特異性
集中荷重をかけたとき、応力が無限大に発散します。
キルヒホッフ板理論では、点荷重のたわみは $w = P/(8\pi D) \cdot r^2 \ln r$ で有限だが、曲げモーメントは $r \to 0$ で対数的に発散する。FEMではメッシュを細かくするほどモーメントが増加し続ける。
対策：
- 集中荷重を微小面積の分布荷重に変換
- 点荷重近傍の応力は信用しない（Saint-Venantの原理で離れた位置を評価）
- 実構造では完全な点荷重は存在しないから、接触面積を考慮
まとめ
キルヒホッフ板のトラブル対処、整理します。
理論の限界を知っていることが、FEMの結果を正しく解釈する鍵ですね。
まさにそう。キルヒホッフ板理論を知らずに板の曲げをFEMで解析するのは、教科書を読まずに試験を受けるようなものだ。
厚板での過小たわみ問題
キルヒホッフ板要素を板厚比（t/L）が0.1以上の厚板に適用すると、せん断変形を無視するためたわみを過小評価する。t/L=0.2の単純支持正方形板でキルヒホッフ要素を使うと理論値より約10%過少なたわみが出る。この限界を超える場合はMindlin-Reissner要素への切り替えが基本で、Abaqusではt/L>0.1でS4R5（Kirchhoff）からS4R（Mindlin）への自動切替え警告が出る。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——キルヒホッフ板理論の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告