刘明超，李 慧，刘俊杰

(1.山东大学土建与水利学院，山东济南 250002;2.北京理工大学宇航学院，北京 100081)

压杆稳定性分析是工程实践中的常见问题。为了保证机构安全可靠地工作，必须使压杆处于直线平衡状态，如果将压杆的工作载荷控制在临界载荷允许的范围内，则压杆就不会失稳。由此可见，临界载荷的确定对于压杆稳定性的分析是非常重要的。

在实际分析压杆稳定性的问题中，临界载荷的确定通常采用手工计算的方法，往往比较繁琐。本文在经典稳定性理论的基础上将计算过程系统化、公式化，绘出了流程图，并编写了压杆稳定性分析的MATLAB 程序，与手工计算相比，更加方便实用。

1 经典稳定性理论

1.1 大柔度压杆临界应力(欧拉)公式

对于满足λ≥λp的压杆，就称之为大柔度压杆，也即细长压杆。其在失稳时处于线弹性范围(即弹性压杆)，临界应力计算可采用欧拉公式:

式中:λ为压杆的柔度。

1.2 中柔度压杆临界应力经验公式

对于满足λs≤λ≤λp的压杆，称为中柔度压杆。其失稳时横截面上的应力超过了材料的比例极限σp，但未超过屈服极限，材料仍满足非线弹性应力应变关系，压杆将发生非线弹性失稳。

对中柔度压杆，临界应力的计算公式是以实验结果为依据归纳出的经验公式，常用如下直线公式:

式中:λ为压杆的柔度;a、b 均为材料常数，可查表得出。

1.3 小柔度压杆临界应力公式

对于满足λ＜λs的压杆，称为小柔度压杆。其失稳时横截面上的应力理论上超出了σs，故其不会发生失稳，它的破坏是由于轴向压缩的强度不够造成的。……