参数设置
假定两端铰接、跨中集中横向荷载。一阶弯矩 $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 时翻倍）；红色虚线 = $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=4 m, P=100 kN, H=10 kN）下，「计算结果」卡片已经显示弯矩比一阶值大了 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$，屈曲荷载与长度平方成反比，所以柱一变细长就急剧危险。L 从 4 m 增大到 10 m，$P_{cr}$ 约缩到 1/6.25。同样的 P = 100 kN，r 就超过 0.5，ψ 翻倍以上。设计规范用「长细比 $\lambda = L/i$」（i 为截面回转半径）评估；中国 GB、AISC、Eurocode 均规定 λ > 100 的柱必须做二阶分析。在本工具拖动 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/(KL)^2$。典型值为两端铰接 K=1.0、两端固定 K=0.5、一端铰一端固定 K=0.7、悬臂 K=2.0。同样长度的两端固定柱，$P_{cr}$ 增大到 4 倍；悬臂柱仅为 1/4。本工具针对最基本的两端铰接情形，其他支承条件可通过 K 等效修正柱长后使用。
主要规范用稳定指数 $\theta = P\Delta/(Vh)$ 或轴力比 $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 起就要求考虑 P-Δ。
本工具假设弹性行为，未计入塑性化的影响。实际钢柱在应力接近屈服时有效刚度下降，使屈曲荷载降低，这就是 Engesser 切线模量理论与 Shanley 折减模量理论描述的非弹性屈曲；AISC 标准已经把它纳入设计公式。轧制或焊接产生的残余应力会引起早期屈服，把 $P_{cr}$ 降至理论值的 60〜80%。各国规范一般通过屈曲折减系数 $\chi$（Eurocode 3 的 a, b, c, d 曲线）综合考虑这些因素。
实际工程应用
高层建筑的抗震设计：30 层以上的高层建筑，地震时楼层层间位移 $\Delta$ 与上部楼层重力荷载（数千吨）相乘形成的倾覆弯矩不能忽略。中国 GB 50011、AISC 341 与 Eurocode 8 都要求当稳定指数 $\theta$ 超过 0.10 时必须把 P-Δ 效应纳入二阶分析。这也是现代高层建筑普遍采用动力非线性时程分析的主要原因之一。
桥墩与大跨桥梁：桥墩承受车辆与列车动荷载、自重，同时还经历温度、收缩、徐变等长期位移。跨径超过 100 m 的大跨桥梁桥墩高度往往超过 50 m，$P_{cr}$ 随长度平方下降，P-Δ 效应支配设计。本州四国连络桥与明石海峡大桥从设计阶段就标准化采用全桥几何非线性分析。
核电厂大型设备支承结构：反应堆压力容器、蒸汽发生器等设备自重达数百至数千吨，重心常处于较高位置。地震时水平位移使 P-Δ 弯矩显著增大。抗震审查依据美国 NRC 的 SRP 与日本电气协会的 JEAG 4601，要求做含二阶效应的弹塑性反应分析。福岛事故后基准地震动大幅强化，P-Δ 效应裕度评估的要求更加严格。
风力发电塔与海洋结构物：陆上 100 m 级、海上 200 m 级风电塔承受自重、旋转扭矩、风载与波浪载，塔顶位移可超过 1 m。$P_{cr}$ 随塔高平方下降，现代大型塔身长细比已超过 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}$。超过设防地震的稀有地震作用下，必须进行包含 P-Δ 的 pushover 或时程分析。
第二个误区是「线性叠加原理仍然成立」。在一阶分析中「$1.0P + 1.0H$ 的响应 = $1.0P$ 响应 + $1.0H$ 响应」成立；但含 P-Δ 的二阶分析中不再成立，因为 $\psi = 1/(1-P/P_{cr})$ 对 P 非线性。在本工具中对比 P=50 kN, H=10 kN 与 P=100 kN, H=10 kN，可以看出 $M_2$ 不是简单按比例增大。多荷载组合评估必须对每个组合分别做二阶分析。
最后一个误区是「只算 $P_{cr}$ 就够了」。本工具的 $P_{cr}$ 是理论欧拉值，未计入局部屈曲、侧向扭转屈曲、初始几何缺陷、残余应力与塑性化。实际钢柱按 Eurocode 3 屈曲曲线（a〜d）或 AISC 360 的 Q 系数修正后，设计屈曲荷载会降至理论值的 60〜90%。H 形钢与开口截面以侧向扭转屈曲（LTB）为主，单凭 $P_{cr}$ 是不够的，需要专门的规范公式或有限元线性屈曲分析（LBA）。请把本工具视为理解物理本质的入门，工程实务中务必结合各国规范与有限元验证。