パラメータ設定
両端ピン支持・中央集中横荷重を想定。一次モーメント $M_1 = HL/4$、増幅係数 $\psi = 1/(1-P/P_{cr})$、二次モーメント $M_2 = M_1\psi$。$r \to 1$ で発散（座屈）。
柱の変形図（一次曲線と P-Δ 二次曲線）
青●＝柱（両端ピン支持）／赤矢印＝中央集中横荷重 H／橙矢印＝軸荷重 P／細線（青破線）＝一次たわみ曲線（軸力なし）／太線（黄）＝P-Δ 増幅後の二次たわみ曲線／中央の白矢印＝中央点たわみ δ（誇張表示）。
増幅係数 ψ = 1/(1−r) の双曲線
横軸＝軸荷重比 $r = P/P_{cr}$ [0〜1]／縦軸＝増幅係数 $\psi$／青線＝双曲線 $\psi = 1/(1-r)$／黄●＝現在の (r, ψ)／橙破線＝$\psi = 2$ ライン（r = 0.5 で 2 倍）／赤破線＝$r = 1$ 漸近線（座屈点で発散）。
理論・主要公式
両端ピン支持柱のオイラー座屈荷重：
$$P_{cr} = \frac{\pi^2 EI}{L^2}$$
軸荷重比と P-Δ 増幅係数：
$$r = \frac{P}{P_{cr}},\qquad \psi = \frac{1}{1-r}$$
中央集中横荷重 $H$ による一次モーメントと P-Δ 増幅後の二次モーメント：
$$M_1 = \frac{H L}{4},\qquad M_2 = M_1\,\psi = \frac{H L/4}{1 - P/P_{cr}}$$
$E$ はヤング率、$I$ は断面二次モーメント、$EI$ は曲げ剛性、$L$ は柱長、$P$ は軸荷重、$H$ は横荷重。$r \to 1$ で $\psi \to \infty$ となり、座屈直前の発散挙動を表現します。
P-Δ 効果とは
🙋
「P-Δ 効果」って聞き慣れない言葉ですけど、要は何ですか？ 普通の梁の曲げ計算と何が違うんですか？
🎓
ざっくり言うと「軸力 P が横変位 Δ に作用して、追加の曲げモーメント P·Δ が生まれる現象」のことだよ。例えば、定規を立てて上から指で押し下げながら、横からも軽く突くと、横変位がどんどん増幅されるよね。あれが P-Δ の本質。普通の一次解析では「変形前の形状」で釣り合いを取るからこの効果は無視されるけど、軸力が大きいと無視できなくなる。本ツールの既定値（L=4m, P=100kN, H=10kN）でも 8.8% モーメントが増えているのが「計算結果」カードでわかるはず。
🙋
スライダーで P を増やすと、増幅係数 ψ がどんどん大きくなって、最後は無限大って表示されました。これって物理的に何が起きてるんですか？
🎓
それがまさに座屈なんだ。$P = P_{cr}$ になると、$\psi = 1/(1-r)$ の分母がゼロになって発散する。これは「軸力単独でも座屈するから、わずかでも横変形が混入すると無限に増幅される」状態。既定値で $P_{cr} \approx 1234$ kN だから、P を 1200 kN 付近まで上げると ψ が一気に跳ね上がるのが右下のグラフで見られるよ。実務では r > 0.6 になったら線形解析は使えず、必ず幾何剛性を含む二次解析（GNLA）に切り替える。
🙋
P-δ と P-Δ ってよく聞きますけど、違いはなんですか？ ツール名は「P-Δ」になってますけど。
🎓
小文字の P-δ（スモールデルタ）は単一部材の中での局所たわみによる増幅で、本ツールが扱っているのはこちらに近い。大文字の P-Δ（ラージデルタ）は多層建物全体の階高変位 Δ に重力荷重が作用する転倒モーメントのことで、地震応答解析や AISC・Eurocode 8 の耐震設計で重要。慣例的にどちらも「P-Δ 効果」と総称することが多くて、本ツールもその慣例に従っている。物理的な本質は同じ「軸力と変位の連成」なんだ。
🙋
柱長 L を 10 m まで伸ばしたら、P_cr が 200 kN 以下まで下がって、ψ が一気に発散しました。これが「細長い柱は危ない」ってことですか？
🎓
その通り！$P_{cr} = \pi^2 EI / L^2$ で、座屈荷重は長さの2乗に反比例するんだ。だから細長くなると急激に危なくなる。例えば L を 4m から 10m にすると $P_{cr}$ は 1/6.25 倍まで下がる。さらに同じ P=100 kN だと r が 0.5 を超えてしまい、ψ が 2 倍以上に。設計実務では「細長比 $\lambda = L/i$（i は断面回転半径）」で評価し、JIS や AISC では λ > 100 を超える柱には P-Δ 効果を含む二次設計が要求される。本ツールで L を動かすと、この感度を直感的に確認できるよ。
よくある質問
本ツールの $\psi = 1/(1-P/P_{cr})$ は AISC 360 や Eurocode 3 が採用している標準的な「モーメント増幅近似」で、両端ピン支持・等しい端部モーメントもしくは横荷重単独の場合に高精度な近似式です。厳密解は微分方程式 $EI\,y'' + P\,y = M(x)$ を解いた三角関数解（$\sec(kL/2)$ を含む式）で得られますが、$P/P_{cr} < 0.5$ の範囲では両者の差は数 % 以内で、設計実務では本式で十分とされています。$P/P_{cr} > 0.7$ を超える領域では誤差が大きくなり、有限要素法による幾何非線形解析（GNLA）が必要です。
境界条件によって有効長 $L_e = K\cdot L$（K は座屈長さ係数）を導入し、$P_{cr} = \pi^2 EI/(K L)^2$ で計算します。K の値は両端ピン K=1.0、両端固定 K=0.5、固定−ピン K=0.7、片持ち（片端固定・片端自由）K=2.0 です。例えば両端固定の柱は同じ材長でも $P_{cr}$ が 4 倍になり、片持ち柱は 1/4 倍になります。本ツールは最も基本的な両端ピン支持を扱っており、他の境界条件は K で長さを補正することで近似可能です。
主要な設計規準では、安定性指数 $\theta = P\Delta/(V h)$ や軸荷重比 $r = P/P_{cr}$ で判定します。AISC 360 では $B_2$（層レベル増幅係数）が 1.4 を超える場合、Eurocode 3 では $\alpha_{cr} = F_{cr}/F_{Ed} < 10$（弾性）または $< 15$（塑性）の場合に二次解析を義務付けています。本ツールの ψ が 1.4 を超える設定（r > 0.286）が、おおまかな閾値の目安になります。耐震設計では地震時の付加変位が大きくなるため、より厳しい基準（ψ > 1.1 程度から考慮開始）が適用されることもあります。
本ツールは弾性挙動を仮定しているため、塑性化の影響は含まれません。実際の鋼構造柱は降伏応力に近づくと弾性係数が実効的に低下し、座屈荷重も下がります。これを考慮するのが「接線係数理論（Engesser）」や「換算係数理論（Shanley）」で、AISC では「インエラスティック座屈曲線」が標準で組み込まれています。残留応力（圧延・溶接由来）は早期降伏を引き起こし、$P_{cr}$ を理論値の 0.6〜0.8 倍まで下げることもあります。実務では座屈低減係数 $\chi$（Eurocode 3 の曲線 a, b, c, d）でこれらをまとめて補正します。
実世界での応用
高層ビルの耐震設計：30 階建てを超える高層ビルでは、地震時の層間変形 $\Delta$ に上階の自重 $P$（数千トン）が作用し、転倒モーメント $P\Delta$ が無視できないレベルに達します。日本の建築基準法施行令や AISC 341、Eurocode 8 では、安定性指数 $\theta$ を計算して 0.10 を超える階で P-Δ 効果を二次解析に含めることを要求しています。これが現代の高層ビルが「動的非線形時刻歴応答解析」を実施する主要な理由の一つです。
橋梁の橋脚と長大スパン：橋脚は車両・列車の動的荷重と自重を支えつつ、温度変化・収縮・クリープによる長期変位を受けます。スパン 100 m を超える長大橋では橋脚高さも 50 m を超えることがあり、$P_{cr}$ が長さの 2 乗で減少するため P-Δ 効果が支配的に。本州四国連絡橋や明石海峡大橋では、設計時から幾何非線形を含む全橋解析が標準化されています。
原子力プラントの大型機器サポート：原子炉容器や蒸気発生器の支持構造は、機器自重（数百〜数千トン）が高所に集中し、地震時の水平変位で P-Δ モーメントが急増します。耐震審査では SRP（米 NRC の標準審査計画）や JEAG 4601（日本電気協会）に従い、F-Δ 効果を含む弾塑性応答解析が義務付けられています。福島事故以降は基準地震動が大幅に強化され、P-Δ 効果による余力評価がいっそう厳しくなりました。
風力発電タワーと洋上構造物：陸上 100 m 級・洋上 200 m 級の風力発電タワーは、自重と回転トルクに加えて風荷重・波浪荷重を受け、頂部変位が 1 m を超えることもあります。$P_{cr}$ がタワー高さの 2 乗で減るため、現代の大型タワーは細長比 100 を超える領域で運用され、必ず幾何非線形を含む設計が行われます。本ツールで L=10 m、EI=10000 kN·m² と入れると、これに近い細長比領域の感度を確認できます。
よくある誤解と注意点
最も多い誤解は、「軸力が小さいから P-Δ は無視できる」と早合点することです。確かに r が 0.1 程度なら ψ は 1.11 で増幅は 11% に留まり、初期設計段階では無視されがちです。しかし耐震設計では「中地震レベルで r = 0.1」が「大地震レベルで一時的に r = 0.5（ψ = 2.0）まで上昇」する可能性があり、降伏後の塑性ヒンジ形成で剛性が下がると実効 $P_{cr}$ もさらに下がります。設計地震を超える稀有な地震動に対する「2 次設計」では、必ず P-Δ を含むプッシュオーバー解析や時刻歴応答解析が要求されます。
次に多いのが、「線形重ね合わせ原理が成り立つ」と思い込む誤りです。一次解析では「荷重 1.0×P + 1.0×H の応答 = 1.0×P 応答 + 1.0×H 応答」が成り立ちますが、P-Δ を含む二次解析では成り立ちません。$\psi = 1/(1-P/P_{cr})$ が P に対して非線形だからです。本ツールで P=50 kN, H=10 kN と P=100 kN, H=10 kN を比較すると、M₂ が単純な比率では増えないことが確認できます。複数荷重ケースの組み合わせ評価では、必ず各荷重ケースを丸ごと二次解析で評価しなければなりません。
最後に、「P_cr さえ計算すれば十分」と考える誤解です。本ツールで扱う P_cr は理論オイラー値で、断面の局部座屈・横倒れ座屈・初期不整・残留応力・塑性化を考慮していません。実際の鋼構造柱では、Eurocode 3 の座屈曲線（a〜d）や AISC 360 の Q 値補正で、理論値の 60〜90% まで設計座屈荷重が下げられます。さらに、ねじり・横倒れ座屈（LTB）が支配的になる H 形鋼や開断面では、P_cr 単独の評価では不十分で、専用の規準式や FEM 線形座屈解析（LBA）が必要です。本ツールは現象の本質を理解するための入門ツールと位置づけ、設計実務では各国規準の適用と FEM 検証を組み合わせてください。