【JY】轴压稳定两百年:从欧拉公式到现代设计方法（稳定篇）

写在前文

一根细长的杆件在压力作用下会发生什么？这是一个看似简单却蕴含深刻力学原理的问题。当我们逐渐增加压力，杆件并不会立即破坏，而是会在某个临界点突然发生"侧向弯曲"——这就是失稳现象，也是结构工程中最危险、最难以预测的破坏模式之一。为了理解和预测这一现象，欧洲数学家和工程师们花费了近三百年的时间，从欧拉的数学推导到现代的多因素耦合分析，逐步建立起了一套完整的轴压稳定理论体系。本文将系统回顾这段跨越两个多世纪的科学探索历程，揭示理论发展与工程实践之间相互促进的辩证关系。

一、欧拉时代：稳定理论的奠基

1.1 问题的起源

十八世纪中叶，随着工业革命的萌芽，工程师们开始遇到越来越多的细长受压构件设计问题。当时人们已经注意到一个令人困惑的现象：同样截面的柱子，越是细长，越容易被"压弯"，而且这种破坏往往发生得非常突然，事前几乎没有任何征兆。这种破坏模式与传统材料强度理论描述的情况截然不同——杆件并没有因为应力超过材料强度而破坏，而是在应力远低于强度极限的情况下突然丧失承载能力。这一反直觉的现象引起了数学家们的极大兴趣。

在欧拉看来，细长杆件在压力作用下的弯曲行为，本质上是一种平衡状态的转换过程：当压力较小时，直线形态是稳定的平衡状态；但当压力达到某个临界值时，这种平衡状态就会变得不稳定，任何微小的扰动都可能导致杆件偏离直线位置，进入弯曲状态。这一洞见成为后来整个稳定理论的基石。

1.2 欧拉公式的推导

1757年，欧拉在发表的著名论文《关于柱的强度》中，首次推导出了轴压杆临界荷载的基本公式。

将临界压力除以截面面积，可以得到临界应力的表达式：

这就是著名的欧拉双曲线，它描述了临界应力随长细比变化的规律。公式表明，临界应力与长细比的平方成反比，这意味着当长细比增大时，临界应力会急剧下降。这解释了为什么细长杆件会在很低的应力水平下就发生失稳破坏。

1.3 欧拉公式的物理意义

欧拉公式的深刻物理意义在于，它将失稳问题从传统的强度破坏范畴中分离出来，建立起了一套全新的分析框架。失稳不再被视为材料的破坏，而是被理解为一种平衡状态的转换。

这一认识对于结构工程具有深远的意义。它告诉工程师们，在设计中不仅要考虑应力是否超过材料强度，还要考虑结构的稳定性是否能够得到保证。对于细长受压构件，稳定问题往往比强度问题更加关键，这彻底改变了传统的设计思路。

二、理论与现实的博弈：十九世纪的探索

2.1 欧拉公式的局限

十九世纪的工程师们发现，当杆件的长细比较小时，按照欧拉公式计算得到的临界应力往往会超过材料的屈服强度，这显然是不可能的。换言之，欧拉公式只适用于"细长"杆件，而不适用于"短粗"杆件。对于后者，破坏模式实际上已经从失稳转变为材料屈服，欧拉的理论就不再适用了。

2.2 兰金的贡献

1858年，苏格兰工程师兰金提出了一个富有创意的经验公式，他假设杆件的承载能力同时受到失稳和屈服两个因素的控制，并将两者以倒数形式组合起来，得到了兰金公式：

兰金公式的核心思想是，杆件的承载能力不是单纯由失稳或屈服决定的，而是两者的综合效果。当杆件非常细长时，失稳起主导作用，公式退化为欧拉公式；当杆件非常短粗时，屈服起主导作用，公式退化为屈服荷载。这一思路虽然缺乏严格的理论推导，但在工程实践中表现出了良好的适用性。

2.3 泰特迈尔的实验研究

泰特迈尔的研究发现了几个重要现象。首先，他证实了杆件的破坏模式确实随长细比而变化：细长杆发生弹性失稳，短粗杆发生材料屈服，中等长细比的杆件则表现出复杂的弹塑性失稳行为。其次，他发现初始缺陷对承载能力有显著影响。

基于大量试验数据，泰特迈尔整理出了一套实用的稳定系数表格，将杆件按长细比分为不同的类别，每个类别采用不同的经验公式。这种"分档处理"的思想后来被各国规范广泛采纳，发展成为柱子曲线的概念。泰特迈尔的工作标志着稳定研究从纯理论探索走向理论与实验相结合的发展阶段。

三、弹塑性理论的突破：二十世纪的革命性进展

3.1 切线模量理论

1947年，美国学者香利提出了一个革命性的概念——切线模量理论，这一理论彻底改变了人们对弹塑性失稳问题的认识。

香利指出，在弹塑性失稳分析中，应该用切线模量    代替弹性模量，从而得到切线模量理论公式：

这一看似简单的替换实际上蕴含着深刻的物理意义。当杆件发生弹塑性失稳时，截面的一部分已经进入塑性状态，这部分材料的刚度已经降低。此时，杆件整体的等效刚度不再是弹性模量，而应该是反映当前应力状态的切线模量。

香利的理论解决了长期困扰研究者们的理论难题，成为现代稳定理论的基石。

3.2 双模量理论与切线模量理论的争论

实际上，可以证明切线模量理论给出的是弹塑性临界荷载的下限，而双模量理论给出的是上限。真实情况介于两者之间，但由于实际杆件存在的初始缺陷，切线模量理论的预测更加接近实际。这场持续半个多世纪的争论终于有了明确的答案，切线模量理论被广泛接受为弹塑性稳定分析的基础理论。

3.3 残余应力的影响

美国学者比德尔和乌尔等人对此进行了大量深入研究。他们发现，焊接残余应力的峰值可以达到材料屈服强度的30%到50%。当构件承受外部压力时，这些残余应力会与外部荷载产生的应力叠加，使得部分截面提前进入塑性状态。由于塑性材料的刚度远低于弹性材料，这会导致构件整体刚度的降低，从而降低其稳定承载能力。

3.4 柱子曲线的建立

1970年代，欧洲钢结构协会和美国结构稳定研究委员会分别建立了系统的柱子曲线体系。所谓柱子曲线，就是稳定系数随正则化长细比变化的关系曲线。不同的曲线代表不同类型的截面，反映了它们在残余应力、初弯曲等方面的差异。

这种分类方法体现了"区别对待"的设计理念：对于残余应力小、初始缺陷少的构件，可以采用较高的曲线；反之，则需要采用较低的曲线。这种精细化的设计方法，既保证了安全性，又兼顾了经济性，是现代概率极限状态设计思想的具体体现。

四、现代设计方法的形成与发展

4.1 稳定系数的确定方法

现代设计规范中，稳定系数的确定是一个核心问题。稳定系数综合反映了构件抵抗失稳的能力，其确定需要考虑众多因素的影响，包括材料特性、截面几何参数、残余应力分布、初始几何缺陷等。

在弹性失稳阶段，稳定系数可由欧拉公式直接得出：

这表明在双对数坐标系中，弹性失稳区的柱子曲线是一条斜率为-2的直线。

4.3 柱子曲线的计算公式

中国规范中柱子曲线采用分段函数的形式表达。当正则化长细比较小时（    ），采用抛物线公式；当正则化长细比较大时（），采用较为复杂的解析公式。

这套公式是由李开禧教授提出的，也称为"李开禧公式"。公式采用了Perry-Robertson公式的形式，能够很好地拟合试验数据和数值分析结果。

实际设计时，轴心受压构件的稳定计算公式为：

设计人员只需要根据构件的具体情况，选择合适的柱子曲线，然后查表或计算得到相应的稳定系数。这种简便的设计方法，将复杂的理论研究成果转化为工程师可以方便使用的工具，体现了理论与实践的良好结合。

五、前沿研究与新时代的发展

5.1 高强钢的稳定问题

研究表明，高强钢的应力-应变曲线与普通钢材有显著不同。普通钢材具有明显的屈服平台，而高强钢往往没有屈服平台，而是直接进入应变硬化阶段。这种差异导致高强钢构件的稳定行为与普通钢构件有所不同。此外，高强钢的残余应力相对值通常比普通钢材小，这对其稳定性能有有利影响。

目前，针对高强钢的稳定设计方法仍在不断完善之中，是当前稳定研究的热点领域之一。

5.3 可靠度分析与设计

传统的安全系数方法存在一个根本性缺陷：它无法定量评估结构的失效概率。可靠度分析方法 正是为解决这一问题而发展起来的。

在可靠度分析框架下，结构的功能函数可表示为：  。当    时，结构失效。

现代设计规范已经普遍采用了基于可靠度理论的极限状态设计方法。这种方法不仅能够保证结构的安全性，而且可以使不同类型构件的安全储备趋于一致，实现设计的优化。可靠度分析方法的推广应用，标志着结构工程设计从经验走向科学的重大进步。

5.4 新材料与新结构

以铝合金为例，其应力-应变关系通常采用Ramberg-Osgood模型描述：

由于铝合金的非线性应力-应变关系，其稳定分析比钢材更为复杂。

六、结语：一堂生动的历史课

回顾轴压稳定理论两百多年的发展历程，我们可以得到许多深刻的启示。这一历程不仅是科学知识的积累过程，更是理论与实践相互促进、共同发展的生动范例。

从认识论的角度看，稳定理论的发展体现了科学认识从简单到复杂、从片面到全面的一般规律。欧拉公式是理想化条件下的理论产物，它抓住了稳定问题的本质特征，但无法直接应用于实际工程。兰金公式和泰特迈尔的实验研究，是在理论与实际之间架起桥梁的初步尝试。切线模量理论的提出，则标志着人们对弹塑性失稳机理认识的深化。而柱子曲线体系的建立，更是综合考虑多种因素影响的系统化成果。

从方法论的角度看，稳定研究展示了理论分析、实验研究和数值模拟三种方法相互补充、相互验证的重要价值。理论分析能够揭示问题的内在机理，建立清晰的概念框架；实验研究能够发现新现象，验证理论预测的准确性；数值模拟能够解决复杂的实际问题，扩展理论的应用范围。三者的有机结合，是推动稳定理论不断进步的强大动力。

一根柱子看似简单，实则蕴含着深厚的力学原理。从欧拉的数学推导，到现代的多因素耦合分析，稳定理论的发展凝聚了几代学者的心血和智慧。这段历史不仅是结构力学的宝贵财富，更是科学研究方法论的生动教材。当我们仰望高楼大厦中的每一根立柱时，不妨想一想：它承载的不仅是建筑的重量，更是两百年科学探索的历史积淀。