基于PERA SIM的压杆特征值屈曲分析

工程结构中有很多构件可以看作杆件，即一个方向的尺寸远大于另外两个方向的尺寸。此类构件受压时即为压杆。对于大柔度杆，通常临界压应力小于弹性极限，即发生弹性失稳。

通常，我们会利用有限元法进行压杆的特征值屈曲分析。压杆的特征值屈曲分析属于线性分析。特征值屈曲分析得到的是屈曲载荷和相应的失稳模态，分析简单，计算速度快，在实际工程中应用大。

PERA SIM Mechanical是安世亚太自主开发的一款机械仿真分析软件，基于PERA SIM通用仿真软件架构，可以实现与流体、电磁、声学等物理场的耦合分析计算。PERA SIM Mechanical求解器可以完成如下分析功能：结构静力学分析、模态分析、瞬态动力学分析、谐响应分析、反应谱分析、屈曲分析、随机振动分析、热分析(稳态+瞬态)以及并行计算。

本文借助结构有限元软件中的屈曲分析模块完成了压杆特征值屈曲分析，展现了软件丰富的操作功能，并且与国际成熟软件的计算结果对比，验证了计算的准确性，为学者和工程师提供了特征值屈曲分析的一种新方法。

2.1 压杆失稳现象

让我们回忆一下几类关于平衡的概念：稳定平衡、不稳定平衡以及条件稳定平衡。

稳定平衡：凡能在被移动离开它的平衡位置后，仍试图恢复其原来位置，从而恢复到原来的平衡状态的物体，它原来的平衡状态叫“稳定平衡”。

不稳定平衡：处于平衡状态的物体，由于受到某种外界微小的作用，如果物体稍有偏离就不能恢复到原来的平衡状态，这种情况叫“不稳定平衡”。

条件稳定平衡：在一定条件下，可以实现稳定平衡；否则，则处于不稳定平衡。

构件的稳定性是指构件在外力作用下保持其原有平衡状态的能力。平衡形式的突然变化称为稳定失效，简称失稳或屈曲。因此，屈曲问题本质上是平衡性问题。一般的，如拉杆或内压圆筒，呈现稳定平衡；如压杆或外压圆筒，呈现条件平衡。

此处，特别讨论受压杆件的稳定性问题。受压应力作用的结构，当受到外界扰动时，随着压应力的增大，会突然丧失原有几何形状而失稳。对应的压应力为临界压应力，对应的载荷为临界载荷。扰动是普遍存在的，如材料不均匀缺陷，成形制造几何缺陷等。因此失稳(又称屈曲)对于受压应力的构件而言是内在的特性。为了防止屈曲失效，需确定结构的临界载荷或临界压应力，使实际载荷或压应力小于许用值。

当构件失稳时，一种内在的本质是原本沿截面厚度均匀分布的压应力，随着变形的增大跃变为弯曲应力，截面因抗弯能力不够而不能维持原有几何形状。

当结构的抗弯截面尺寸较小时，如细长的大柔度杆，薄壁圆筒等，失稳时截面的压应力往往低于材料的弹性极限，这种失稳称为弹性失稳。但当结构的抗弯截面尺寸较大时如大柔度杆，壁厚较厚的圆筒，失稳时截面的压应力往往高于材料的弹性极限出现塑性变形，这种失稳称为弹塑性失稳。

临界载荷或临界压应力的大小首先与抗弯刚度有关，对于弹性失稳，弹性模量越大，抗弯刚度越大，抗弹性失稳能力越强。

2.2 特征值屈曲分析

特征值屈曲分析预测了理想弹性结构的理论屈曲强度。该方法与教科书中的弹性屈曲分析方法相一致。欧拉柱的特征值屈曲分析与经典欧拉解相匹配。缺陷和非线性行为阻碍了大多数现实世界结构实现其理论弹性屈曲强度。因此，特征值屈曲通常产生非保守结果，因为它没有考虑到这些影响。

虽然特征值屈曲分析是非保守的，但与非线性屈曲解决方案相比，它具有计算成本低的优点，并且可以提供近似的(尽管是非保守的)屈曲条件预测。下图为屈曲的加载曲线(线性和非线性)。

图1 屈曲加载曲线(线性和非线性)

在特征值屈曲分析中，求解以下特征值问题，得到屈曲载荷乘子λi和屈曲模态ψi。在该解中，假设[K]和[S]矩阵为常数，即解是线性的。

其中，{ψi}为特征向量—屈曲模态；[S]为预应力刚度矩阵；λi为载荷乘子；[K]为结构刚度矩阵。

注意事项：

1）需要施加载荷，进行预应力计算，得到预应力刚度矩阵。

2）求矩阵的特征值，获取使行列式值为零的载荷乘子，即可得到临界载荷。

2.3 压杆稳定理论解

对于大柔度压杆，通常临界压应力小于弹性极限，发生弹性失稳。屈曲载荷可用欧拉公式求解。

   ，

μ—支撑系数，两端简支时取1.0；两端固支时取0.5；一端固支一端简支时取0.707。