ティモシェンコ梁理論

最大の違いはせん断変形を考慮することだ。オイラー・ベルヌーイ梁では「断面は中立軸に常に直交する」と仮定したが、ティモシェンコ梁ではこの仮定を外す。断面は中立軸から傾くことが許される。

数学的には：

その通り。EB梁では $\gamma = 0$（せん断変形なし）を強制していた。ティモシェンコ梁ではこの拘束を外すことで、より一般的な梁の挙動を記述する。

2つの連立微分方程式になる：

梁理論では断面のせん断応力を一様と仮定するが、実際のせん断応力分布は放物線状だ。$\kappa$ はこの差を補正する係数。

そう。I形鋼のせん断はほとんどウェブが負担するから、有効せん断面積はウェブ面積に近い。フランジは曲げには寄与するがせん断にはほとんど寄与しない。

できる。全たわみは：

せん断たわみの比率：

計算すると約3%。$L/h = 5$ なら約12%。$L/h = 3$ なら約32%。$L/h < 5$ でせん断変形が無視できなくなるという経験則はこの計算から来ている。

ティモシェンコ梁要素は $w$ と $\theta$ を独立な変数として扱う。EB梁要素では $\theta = dw/dx$ の拘束があったが、ティモシェンコ梁要素ではこの拘束がない。

ただし注意が必要。通常の2節点要素（線形補間）ではせん断ロッキングが起きる。

EB梁のシアロッキングとは逆の現象だ。ティモシェンコ梁要素で純粋な曲げ変形を表現しようとすると、寄生的なせん断ひずみが発生して要素がロックする。結果としてたわみが過小になる。

対策：

そう。FEM要素設計の永遠のテーマだ。EB梁はせん断を無視してロッキングを避け、ティモシェンコ梁はせん断を含めてロッキングと戦う。実用的な解は低減積分だ。

要点：

その通り。ティモシェンコ梁はEB梁を包含する。$L/h$ が大きい細長い梁ではEB梁と同じ結果になるから、ティモシェンコ梁を常に使っても問題ない。だからAbaqusやAnsysのデフォルト梁要素はティモシェンコ梁なんだ。

2節点のティモシェンコ梁要素で、曲げ自由度のみ（$w_1, \theta_1, w_2, \theta_2$）の場合：

完璧な確認だ。$\Phi$ はEB梁からの補正量と考えてよい。$\Phi$ が大きい（太短い梁）ほどせん断変形の影響が大きく、EB梁との差が出る。

2節点の線形補間要素で純粋曲げ（$M$ = 一定）を表現する場合を考えよう。曲げモーメントが一定なら $\theta$ = 一定（$d\theta/dx = M/(EI)$）であるべきだが、$w$ の線形補間では $dw/dx$ も一定。すると $\gamma = dw/dx - \theta = \text{一定} \neq 0$。

そう。このせん断ひずみのエネルギーが余分に蓄えられるため、変形が過小評価される（硬すぎる応答）。

低減積分（1点Gauss）で解決する理由：$\gamma$ を要素中心（$x = L/2$）でのみ評価すると、線形のせん断ひずみの平均値がゼロになる。つまり寄生せん断が消える。

AbaqusのB21/B31は低減積分を内部的に使っており、せん断ロッキングは自動的に回避される。Nastranの CBEAM も同様の対策がされている。ユーザーが意識する必要は通常ない。

$\kappa = 1.0$ はせん断面積 = 全断面積を意味し、せん断剛性を過大評価する。I形断面で $\kappa = 1.0$ のままにすると、せん断たわみが過小になる。PBEAM カードの K1, K2 を正しく設定すること。NastranのPBEAML（断面形状自動計算）を使えばこの問題は回避できる。

シェル要素でH形鋼をモデル化すれば、断面のせん断応力分布は自動的に正確に出る。$\kappa$ の設定は不要。ただしDOF数は梁要素の10〜100倍になる。

使い分け：

要点：

EB梁では不十分な場面を整理しよう。

サンドイッチパネル（表面板＋コア材）はせん断変形が支配的だ。コア材のせん断剛性が低いため、$L/h = 10$ でもせん断変形が全たわみの30〜50%を占めることがある。サンドイッチ梁にEB梁を使うのは完全に間違いだ。

RC造の耐震壁に設けられる連結梁（カップリングビーム）は $L/h = 2 \sim 4$ と非常に太短い。せん断変形が全変形の50%以上を占めるため、ティモシェンコ梁が必須だ。

高次の振動モードでは波長が短くなるため、実効的な $L/h$ が小さくなる。低次モードではEB梁で十分でも、高次モードではせん断変形の影響が顕著になる。動的解析にはティモシェンコ梁を使うのが安全だ。

スタッドコネクタによる不完全合成（partial interaction）ではスタッドのせん断変形が全体たわみに寄与する。これを梁要素で表現するにはティモシェンコ梁の考え方が不可欠だ。

1. $L/h$ を計算 — スパン/梁せい

そう。現代のFEMソルバーのデフォルト梁要素はティモシェンコ梁だ。敢えてEB梁を選ぶ理由はほとんどない。

そう。この足し算でFEMの結果を検証できる。せん断項が全たわみに占める比率を計算して、ティモシェンコ梁の効果を確認する。

NastranのCBEAMは $\kappa = 1.0$ がデフォルト。Abaqusは断面形状から自動計算。Ansysも自動計算。ソルバーによって異なるから、必ずマニュアルを確認すること。

実務で最も多いミスの一つだ。PBEAMLカードで断面形状（BOX, I, T等）を指定すれば自動計算されるが、PBEAM（直接入力）を使う場合はK1, K2を手動で設定しないとデフォルトの1.0になる。

AbaqusとAnsysは断面形状からの自動計算がデフォルトなので、この問題は起きにくい。ただしABaqusのGENERAL SECTIONやAnsysの直接入力では同様の注意が必要。

断面を多数の「ファイバー」に分割し、各ファイバーに独立した応力-ひずみ関係を与える。これにより断面の塑性化の進行を正確に追跡できる。耐震設計でRC柱やSRC柱の非線形解析に広く使われている。

そう。梁要素自体はシンプルだが、断面パラメータの扱いで実用上の差が生じる。ソルバーのマニュアルをよく読んで、デフォルト値を確認する習慣をつけてほしい。

梁理論の拡張と新材料への適用が活発だ。

高次梁理論（Higher-Order Beam Theory, HOBT）がある。ティモシェンコ梁は断面のせん断ひずみを一定と仮定するが、HOBTではせん断ひずみの分布を放物線やそれ以上の多項式で表現する。

代表的な高次梁理論：

ある意味そうだ。イタリアのCarrera教授が開発したCUFは、展開次数を上げるだけでEB梁→ティモシェンコ梁→高次梁→3次元弾性解に連続的に遷移する。断面の変形（局所座屈的な変形も含む）まで1次元のフレームワークで扱える画期的な手法だ。

複合材（CFRP等）のティモシェンコ梁は異方性材料の扱いが鍵だ。繊維角によって曲げ-せん断-ねじりが連成する。

VABS（Variational Asymptotic Beam Section）はジョージア工科大学のYu教授が開発した手法で、3次元の断面を2次元のFEMで解析し、等価な梁の剛性マトリクス（6×6）を導出する。複合材ブレード（風車、ヘリコプター）の設計で広く使われている。

ソフトロボティクスでティモシェンコ梁理論が活用されている。柔軟なロボットアーム（コンティニュアムロボット）の変形をティモシェンコ梁として記述し、リアルタイムでの制御に使う。大変形に対応するためにCosserat梁（3次元の有限回転を含む幾何学的に正確な梁理論）が使われることが多い。

まさにそう。カテーテルナビゲーションのリアルタイムシミュレーションは、Cosserat梁理論＋接触のFEMで実現されている。医療ロボティクスの基盤技術だ。

EB梁のトラブルに加えて、ティモシェンコ梁特有の問題がある。

正常な結果だ。ティモシェンコ梁はせん断変形を含むから、EB梁の理論値より大きいたわみになる。差はせん断たわみ $PL/(GA_s)$ に対応する。

検証方法：

$\kappa$ の設定ミスが最も多い原因だ。

対処法：

以下の症状が出たらせん断ロッキングを疑う：

ただし主要ソルバーのデフォルト梁要素はせん断ロッキング対策済みだから、通常は問題にならない。ユーザーサブルーチンで独自に実装する場合に注意が必要。

確認：

サン・ブナンねじり定数 $J$ は薄肉開断面の場合：

そう。梁要素のトラブルの大半は断面パラメータの設定ミスだ。自動計算を使えばミスの大部分は防げる。それでも結果がおかしいときは、単純な問題（片持ち梁等）で理論値と照合する。基本に忠実であることが最良のデバッグだ。