特征值屈曲分析的ANSYS实现方法

1907年，在加拿大劳伦斯河上，跨长五百多米的魁北克大桥，因压杆失稳而导致整座大桥倒塌，两次事故造成88人遇难。

在杆、壳结构的设计中，稳定性分析和强度分析是同样重要的。受轴向压力的细长杆，当压力超过一定数值时，压杆会由原来的直线平衡形式突然变弯，致使结构丧失承载能力，平衡形式的突然变化称为稳定失效，简称失稳或屈曲。工程中的柱、桁架中的压杆、薄壳结构在压力作用下，都可能发生失稳。

屈曲分析是一种用于确定结构开始变得不稳定时的临界载荷和屈曲模态形状的技术。

采用ANSYS进行屈曲分析，在各结构领域已经非常成熟，如桁架设计、船体设计、风电叶片设计、飞机地铁壳体的设计等。

ANSYS提供了非线性屈曲分析和线性屈曲分析两种方法。其中线性屈曲分析又称为特征值屈曲分析，时教科书中传统的弹性屈曲分析方法。非线性屈曲分析方法比线性屈曲分析方法更精确，但是应用起来更复杂。本文介绍特征值屈曲分析基本流程。

首先，需要建立结构的有限元模型，包括定义单元类型、材料属性、建立几何模型，划分网格等。

第二、对结构进行静力学分析，在分析之前需要打开预应力选项，因为该分析需要计算应力刚度矩阵。通常需要为结构施加单位压载荷，由屈曲分析计算出的特征值即为临界载荷。

假设某桁架中某构件，其受力可以简化为上图所示模型，细长杆件承受压力，两端铰支。

材料力学中的欧拉公式，临界载荷为Fcr=，对于两端铰支

该构件的截面形状为正方形，边长为0.03m，杆长度为2m，材料为Q235结构钢，弹性模量200GPa

截面惯性矩    

临界载荷     

下面采用ansys进行特征值屈曲分析，求解其临界应力。

FINISH
/CLEAR
/FILNAME,EIGENVALUE_BUCKLE

/PREP7
ET,1,BEAM188

MP,EX,1,2E11
MP,PRXY,1,0.3

SECTYPE,1,BEAM,RECT
SECOFFSET,CENT
SECDATA,0.03,0.03

K,1,0,0,0
K,2,0,2,0
LSTR,1,2
LESIZE,1,,,10
TYPE,1
MAT,1
SECNUM,1
LMESH,1
FINISH

/SOLU
ANTYPE,STATIC
PSTRES,ON
DK,1,UX
DK,1,UY
DK,1,UZ
DK,1,ROTX
DK,1,ROTY

DK,2,UX
DK,2,UZ
DK,2,ROTX
DK,2,ROTY

FK,2,FY,-1
SOLVE

FINISH

/SOLU
ANTYPE,BUCKLE
BUCOPT,LANB,1
MXPAND,1
OUTPR,ALL,ALL
SOLVE
SAVE
FINISH

/POST1
SET,FIRST
SET,LIST
PLDISP,2
FINISH

ANSYS求解的特征值为33791.4，意思是临界屈曲载荷为33791.4N，与欧拉理论解33276相比，误差1.5% 。