力のつり合いと静力学 — 自由体図・モーメントからFEM線形静解析まで
1. 静力学とは — CAEの出発点
静力学（Statics）は「加速度ゼロ＝力のつり合い」の状態にある物体を扱う力学の一分野だ。動いていない物体、あるいは等速直線運動している物体——すなわち $\sum \vec{F} = 0$ が成立する系——が対象となる。
建物の柱、橋梁の主桁、工場クレーンのブーム、航空機の胴体フレームはすべて静力学の問いから設計が始まる。「この荷重に対して反力はいくつか？ 内部応力はどこが最大か？」——これを答えるのが静解析だ。
2. 力のベクトル合成と分解
力はベクトル量で、大きさと向きの両方を持つ。複数の力が作用する場合、合力（合成力）はベクトルの和で求める。
$$\vec{F}_{合} = \sum_{i} \vec{F}_i$$2次元での成分分解:
$$F_x = F\cos\theta, \quad F_y = F\sin\theta$$3次元では方向余弦 $(l, m, n)$ を使って $\vec{F} = F(l\hat{x} + m\hat{y} + n\hat{z})$
よくある成分分解の例
傾斜角 $\theta = 30°$ の斜面上の物体（重量 $W = 1\,\text{kN}$）:
- 斜面平行成分（すべり力）: $W\sin 30° = 0.5\,\text{kN}$
- 斜面垂直成分（垂直抗力）: $W\cos 30° = 0.866\,\text{kN}$
FEMでは荷重ベクトルの成分入力が必要で、この分解を3Dで正確に行うことが境界条件設定の第一歩だ。
3. 力のつり合い条件（静的平衡）
物体が静止しているとき（静的平衡状態）は、次の2条件が同時に満たされている。
$$\sum \vec{F} = 0 \quad \text{（並進平衡：力の和 = 0）}$$ $$\sum \vec{M} = 0 \quad \text{（回転平衡：モーメントの和 = 0）}$$2次元問題では:
$$\sum F_x = 0, \quad \sum F_y = 0, \quad \sum M_O = 0$$この3つの独立な方程式から、最大3つの未知数（反力・内力）を求めることができる。
なぜモーメントのつり合いも必要か
力の和だけがゼロでも、物体が回転してしまうことがある。例えば同じ大きさで反対向きの2力が同一線上にない場合——力の和はゼロだが偶力（Couple）モーメントが生じて回転する。実際の構造物設計では力のモーメントが支配的な設計変数になることが多い（梁の曲げモーメント、ダムの転倒モーメントなど）。
4. モーメント（回転効果）
力 $F$ がモーメントアーム（作用線から支点への垂直距離）$d$ の位置に作用するとき、モーメント $M$ は:
$$M = F \times d \quad \text{[N·m]}$$反時計回りを正（右手法則）とするのが一般的な符号規則。3次元のモーメントベクトルは:
$$\vec{M} = \vec{r} \times \vec{F}$$偶力（Couple）の重要性
偶力は同一大きさで逆向きの2力の組であり、合力はゼロだが純粋な回転効果（モーメント）を生む。偶力モーメントは作用点によらず一定値で「自由ベクトル」の性質を持つ。FEM境界条件でのモーメント荷重入力がこれに相当する。
5. 自由体図（FBD）の描き方
自由体図（Free Body Diagram: FBD）は静力学・動力学問題を解くための最重要ツールだ。物体を周囲から「切り離し」、その物体に作用するすべての外力・反力・モーメントを図示する。
FBDを描くステップ
- 対象を選ぶ：解析したい物体（または部分構造）を系として定義する。
- 切り離す：物体を支持物・隣接物から切り離して孤立させる。
- 外力を描く：重力、荷重、加速度力を正確な作用点・方向で矢印描画。
- 反力を描く：支持点（ピン、ローラー、固定端）の反力を未知数として描く。
- つり合い方程式を立てる：$\sum F_x=0$, $\sum F_y=0$, $\sum M=0$ で未知反力を解く。
FBDとFEM境界条件の対応
FEMモデルの境界条件設定は、実質的にFBDを数値モデルに翻訳する作業だ。
| 支持タイプ | FBD（反力） | FEM境界条件 |
|---|---|---|
| ピン支点（2D） | $F_x$, $F_y$（2成分） | $u_x=0$, $u_y=0$（変位拘束） |
| ローラー支点（2D） | $F_n$（法線方向1成分） | 法線方向 $u_n=0$ のみ |
| 固定端（2D） | $F_x$, $F_y$, $M_z$（3成分） | $u_x=u_y=\theta_z=0$ |
| 固定端（3D） | 6成分（3力+3モーメント） | 全6自由度を拘束 |
6. 静定構造と不静定構造
静定（Statically Determinate）：つり合い方程式だけで全反力・内力が決定できる構造。未知数の数 = 方程式の数。
不静定（Statically Indeterminate）：つり合い方程式だけでは解けない構造。未知数が方程式より多く、変形の適合条件（変位の整合性）を追加しなければならない。
7. トラス構造の節点法・断面法
トラス（Truss）：直線部材を節点でピン結合した構造で、各部材は軸力のみ（引張または圧縮）を受け、曲げモーメントはゼロとなる。橋梁・鉄塔・大型クレームのブーム・宇宙ステーションの骨格構造が典型例だ。
節点法
各節点で力のつり合い方程式（2D: $\sum F_x=0$, $\sum F_y=0$）を立てて部材力を順次求める。節点数を $j$、部材数を $m$、支点反力数を $r$ とすると、静定条件は $m + r = 2j$（2D）。
断面法（リッター法）
特定の部材力だけを素早く求めたいとき、トラスを仮想的に切断し、切断面の部材力を未知数としてつり合い式を立てる。3本までの部材を同時に切断し、モーメントの支点を巧みに選ぶことで連立方程式なしに部材力が求まる。
実例：ハウとプラット形式のトラス橋
| 形式 | 特徴 | 主な用途 |
|---|---|---|
| プラット（Pratt） | 斜材が引張、垂直材が圧縮 | 鉄道橋・道路橋 |
| ハウ（Howe） | 斜材が圧縮、垂直材が引張 | 木製橋梁（木は圧縮に強い） |
| ウォーレン（Warren） | 斜材が引張・圧縮交互 | 軽量橋梁・送電鉄塔 |
| Kトラス | 垂直パネル内にK字ブレース | 長スパン橋梁 |
8. 重心と圧力中心
重心（Center of Gravity, CoG）：物体全体の重力が1点に集中したとみなせる点。
$$\bar{x} = \frac{\sum m_i x_i}{\sum m_i}, \quad \bar{y} = \frac{\sum m_i y_i}{\sum m_i}, \quad \bar{z} = \frac{\sum m_i z_i}{\sum m_i}$$均一密度の場合は形状の幾何学的重心（セントロイド）と一致する。CAEでは質量プロパティ解析（FEMのマスプロパティ出力）で自動計算される。
重心と安定性
物体の重心が支持基盤（台形の底面）の内側にある限り安定だ。重心が外に出ると転倒する。クレーンの転倒事故やフォークリフトの横転はすべて重心移動によるもの。タワークレーンは対面に重い「カウンターウェイト」を配置して常に重心をベース内に保つ。
圧力中心
翼面や水門など分布圧力が作用する面で、合力が作用する等価な1点を圧力中心（Center of Pressure）と呼ぶ。航空機の静安定設計では、圧力中心が重心より後ろにあること（ニュートラルポイント）が重要で、CFD解析で攻撃角ごとの圧力中心位置を求める。
9. FEM線形静解析への完全な架け橋
高校静力学の「$\sum F = 0$」は、FEMの基本方程式 $[K]\{u\} = \{F\}$ の特別な形に相当する。
剛性方程式の組み立て
各要素（ばね、棒、梁、連続体要素）に対して要素剛性行列 $[k]^e$ を求め、全体剛性行列 $[K]$ にアセンブリする。
$$[K]\{u\} = \{F\}$$$[K]$: 全体剛性行列（$N \times N$, $N$ = 総自由度数）
$\{u\}$: 節点変位ベクトル（未知数）
$\{F\}$: 外力ベクトル（既知）
境界条件（固定点の変位ゼロ）を適用して $[K]$ の対応行列・列を消去し（あるいは大値を置いて）連立方程式を解く。残りの節点変位が求まれば、各要素の応力・ひずみを計算できる。
反力の計算
変位 $\{u\}$ が求まった後、支点での反力は:
$$\{R\} = [K]\{u\} - \{F\}$$手計算での反力チェックと一致することを確認するのがFEM検証の基本だ。
線形静解析の仮定と限界
- 変形が小さい（微小変形の仮定）——大変形（つる性の変形など）には幾何学的非線形解析が必要
- 材料は線形弾性——塑性変形には材料非線形解析が必要
- 荷重は時間変化しない——振動・衝撃には動解析が必要
10. 実践：建設重機ブーム構造の設計
クローラークレーンの主ブーム（長さ 30 m、最大吊り荷重 100 t）の強度評価フロー。
荷重条件の設定
最悪ケース（ブーム完全水平、静荷重＋動荷重係数 1.3）:
$$F_{design} = 100 \times 1000 \times 9.81 \times 1.3 \approx 1,275\,\text{kN}$$根元の曲げモーメント: $M = 1275 \times 10^3 \times 30 = 38.25\,\text{MN·m}$
断面の初期設計
箱形断面（外寸 800×800 mm, 板厚 20 mm）の断面係数 $W_Z \approx 0.019\,\text{m}^3$:
$$\sigma_{max} = \frac{M}{W_Z} = \frac{38.25 \times 10^6}{0.019} \approx 2,013\,\text{MPa}$$これは鋼材降伏応力（355 MPa）を大幅に超えるから、より大きな断面またはより高強度鋼（$\sigma_Y = 690\,\text{MPa}$ 以上の高張力鋼）が必要であることがわかる。現実の大型クレーンブームはテーパー箱形鋼+バックステーワイヤーを組み合わせて荷重分散する。
FEM解析による最終確認
- 3Dシェル要素（S4R）でブーム箱形断面をモデル化
- 吊り荷・自重・ワイヤー張力を境界条件として入力
- von Mises応力コンタで最大応力位置を特定
- 安全率 $S = \sigma_Y / \sigma_{max} \geq 1.5$ を確認
- 座屈解析（線形座屈）で圧縮部材のバックリング余裕を確認