ティモシェンコ梁理論 — トラブルシューティングガイド
より充実した内容を beam-timoshenko.html でご覧いただけます。
ティモシェンコ梁のトラブル
ティモシェンコ梁要素でよくあるトラブルを教えてください。
EB梁のトラブルに加えて、ティモシェンコ梁特有の問題がある。
EB梁の理論解と一致しない
理論値 $\delta = PL^3/(3EI)$ と比較してFEMのたわみが大きいです。
正常な結果だ。ティモシェンコ梁はせん断変形を含むから、EB梁の理論値より大きいたわみになる。差はせん断たわみ $PL/(GA_s)$ に対応する。
検証方法：
- ティモシェンコ梁の理論値（EB＋せん断項）と比較して一致するか確認
- $L/h$ を十分大きくして（100程度）解析し、EB梁の理論値に収束するか確認
せん断補正係数のミス
せん断変形が過大または過小に出ます。
$\kappa$ の設定ミスが最も多い原因だ。
| 状況 | 考えられる原因 |
|---|---|
| せん断変形が過大 | $\kappa$ が小さすぎる（またはゼロ） |
| せん断変形が過小 | $\kappa = 1.0$（デフォルト）のまま |
| I形断面で結果がおかしい | $\kappa$ をウェブ面積ではなく全断面積で計算 |
対処法：
- 断面形状自動計算を使う — PBEAML(Nastran), *BEAM SECTION(Abaqus)
- 手動設定する場合は理論値と照合 — 矩形: 5/6, 円形: 6/7, I形: $A_w/A$ 程度
せん断ロッキング
せん断ロッキングが起きているかどうかは、どう判別しますか？
以下の症状が出たらせん断ロッキングを疑う：
- たわみがEB梁の理論値より小さい（硬すぎる応答）
- メッシュを細かくすると結果が大きく変わる
- $L/h$ が大きい梁でもEB梁より硬い結果が出る
ただし主要ソルバーのデフォルト梁要素はせん断ロッキング対策済みだから、通常は問題にならない。ユーザーサブルーチンで独自に実装する場合に注意が必要。
梁のねじりが出ない
ティモシェンコ梁でねじりを解析したいのですが、ねじり変形が出ません。
確認：
- ねじり定数 $J$ が設定されているか — $J = 0$ だとねじり剛性がゼロ
- ねじりの自由度が拘束されていないか — 端部で全自由度を拘束すると回転も止まる
- 偏心荷重が正しく与えられているか — せん断中心からオフセットした荷重がねじりを発生
H形鋼のねじり定数はどう計算しますか？
サン・ブナンねじり定数 $J$ は薄肉開断面の場合：
フランジとウェブの板幅 $b_i$ と板厚 $t_i$ の3乗の和。フィレットがあれば10〜20%増加する。断面形状自動計算を使えば正確に出る。
まとめ
ティモシェンコ梁のトラブル対処、整理します。
結局「断面形状自動計算を使え」に帰着しますね。
そう。梁要素のトラブルの大半は断面パラメータの設定ミスだ。自動計算を使えばミスの大部分は防げる。それでも結果がおかしいときは、単純な問題（片持ち梁等）で理論値と照合する。基本に忠実であることが最良のデバッグだ。
剪断ロッキングの診断法
ティモシェンコ梁要素で細長い梁を解析すると、フルインテグレーション要素では剪断ロッキングで変位が著しく過小評価される。診断法として「参照解の1/100以下なら要疑い」が実務の目安だ。Abaqus B21要素は1点減少積分でロッキングを自動回避するが、NASTRANのCBEAMはユーザーが積分点数を手動で設定する必要がある。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ティモシェンコ梁理論の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告