白 旭，李金华

(1.江苏科技大学 船舶与海洋工程学院，江苏 镇江 212003；2.齐鲁工业大学（山东省科学院），山东省科学院海洋仪器仪表研究所，山东 青岛 266061)

0 引 言

圆柱壳结构具有特殊的几何构型，对于轴压、径压、扭转以及组合载荷均具有很好的承载作用，在船舶与海洋工程等各类领域中都有十分广泛的应用。圆柱壳的主要失效模式为失稳，针对其失稳稳定性问题，特别是轴压作用下的圆柱壳失稳载荷和模态，一直是学者研究最为活跃的课题之一[1]。

早期针对轴压作用下薄壁圆柱壳的稳定性试验表明结构实际的失稳承载能力和理论上理想的完善壳体相比较小，原因就在于初始缺陷的存在会引起轴压薄壁圆柱壳临界失稳载荷很大程度的折减。祝恩淳等[2]针对轴压作用下薄壁圆柱壳的试验数据，采取统计分析的方法给出了结构失稳应力的经验计算公式。Lauren Kougias[3]采用有限元数值方法研究了壳体不圆度对轴压圆柱壳稳定性的影响，认为圆柱壳不圆度与壳体厚度相当时，临界承载能力降低可达75%。Muggeridge等[4]研究认为局部轴对称缺陷对圆柱壳失稳临界载荷的影响很大。龚友根等[5]研究凹坑形式的轴对称缺陷对轴压圆柱壳失稳承载能力的影响。王小明[6]研究了纵骨对环肋圆柱壳壳板稳定性的影响。

一系列的研究表明，研究分析方法对于问题的适用性以及结果准确度存在突出的影响。针对圆柱薄壳结构的失稳问题，Croll 等[7-8]提出的刚度折减方法（Reduced Stiffness Method，RSM）在一定程度上理论值可以与实验值很好的吻合，随即引发了基于RSM 开展不同类型以及不同承载的薄壁壳体结构失稳载荷的研究。Rossana C.Jaca 等[9]应用RSM 法研究了风载荷作用下圆柱壳结构的失稳失稳问题，认为采用RSM 方法得到的失稳下界与非线性分析结果较为一致。Hongtao Wang 等[10]采用RSM 方法对圆柱壳体结构的失稳承载能力进行了优化研究，认为在给定材料和几何参数的范围内，改善了失稳临界载荷。Eduardo M.Sosa 等[11]基于RSM 方法分析了风载荷作用下含有初始缺陷的圆锥壳结构失稳临界载荷。Xu Bai 等[12]分析了刚度折减对环肋圆柱壳承载能力的影响，研究结果表明周向薄膜刚度对圆柱壳的稳定性影响最大，轴向的弯曲刚度对圆柱壳稳定性的影响较小。

在现有基础上，本文采取RSM 针对具有轴对称缺陷的轴压圆柱壳的稳定性进行研究。

1 轴压下完善圆柱壳的失稳载荷

1.1 失稳载荷

轴压作用下完善薄壁圆柱壳模型如图1 所示。圆柱壳两端自由支持在刚性周界，设定长度为l，壁厚为t，半径为 r，均匀轴向压力为 p，材料的弹性模量为E，位移为(U,V,W)。

依据薄壳理论，圆柱壳失稳前的应力状态下：

失稳后的平衡状态为：

图1 轴压完善圆柱壳模型Fig.1 Perfect cylindrical shell model under axial compression

失稳后，总应力、应变、弯矩、转角可以写为：

其中公式右边第1项为失稳前，第2项为失稳后相对增量。

将应变、应力分成线性项部分(ε′x,ε′θ,ε′xθ)，(n′x,n′θ,n′xθ)和二次项部分(ε′′x,ε′′θ,ε′′xθ)，(n′′x,n′′θ,n′′xθ)，忽略高阶项，则圆柱壳总势能为：

其中：

圆柱壳位移选取为：

其中：i 为圆周方向的波数；j 为轴线方向的半波数；uij，vij，wij为任意常数。

依据文献[13]中圆柱壳几何方程和胡克定律，可以得到：

其中：λ=jπ/L，L=l/r。

其中：α=12R2，R=r/t。

由于总势能二次项部分在任意条件下的偏导为0，因而可转化为特征值方程问题:

其中：Λ=p/E。

式（16）存在任意解的条件是系数行列式为0，因而

1.2 结果分析

图2 给出了应用式（18）得到的轴压作用下完善圆柱壳结构的失稳载荷（L=2.88，R=250，t=5，以下计算均以此参数为基础），可以看出轴压圆柱壳存在多种失稳模态。

图2 轴压圆柱壳失稳载荷图谱Fig.2 Instability load spectrum of axial compression cylindrical shell

从图3 和图4 中可以看出，承受轴向压力圆柱壳的失稳载荷仅与圆柱壳半径与厚度的比R 有关，与长度和半径的比L 无关。因此在式（18）中，求其最小值，得

与经典弹性临界应力公式相同。

2 轴压圆柱壳的刚度折减分析

依据式（11），轴压完善圆柱壳失稳总势能中的V2项可以表述为：

图3 不同L 下的失稳载荷Fig.3 Instability load at different L

图4 不同R 下的失稳载荷Fig.4 Instability load at different R

通过壳体的稳定条件，将轴向压力下圆柱壳稳定性问题转化为求解特征方程的特征值问题：

式中：Km和 Kb分别是薄膜刚度矩阵和弯曲刚度矩阵；KG几何刚度矩阵；参数 Λc是特征值而向量Φ是特征向量，对应于圆柱壳的失稳载荷和失稳模态。

当圆柱壳各刚度折减时，式（22）应写为：

其中：Λ′和Φ′为折减刚度后的新的特征值和特征向量。β为薄膜刚度折减系数（0＜β≤1），ξ 为弯曲刚度折减系数（0＜ξ≤1）。

因此式（18）中的系数将变为：

可得含有折减系数的临界载荷方程：

通过分析，临界载荷方程分子的两项分别代表薄膜刚度和弯曲刚度对临界载荷的贡献，在公式（24）中，求其最小值得：

图5 和图6 给出了不同刚度折减系数对轴向受压薄壁圆柱壳失稳载荷的影响。

可以看出，无论是薄膜刚度还是弯曲刚度的降低都会影响轴压圆柱壳的失稳模态。轴压圆柱壳失稳时周向波数随着薄膜刚度的减少而减少，随着弯曲刚度的降低而增加。

图5 弯曲刚度折减失稳载荷图谱Fig.5Instability load spectrum of bending stiffness reduction

图6 薄膜刚度折减失稳载荷图谱Fig.6 Instability load spectrum of film stiffness reduction

3 含轴对称缺陷圆柱壳失稳分析

3.1 基于RSM 的失稳载荷计算

含轴对称初始缺陷圆柱壳模型如图7 所示。设定长度为l，半径为 r，初始缺陷的幅值为 w0，缺陷的曲率为ρ，静水压力为 p，壳板的厚度为t。

图7 含轴对称缺陷轴压圆柱壳模型Fig.7 Axially compressed cylindrical shell model with axisymmetric defects

缺陷幅值相对于半径为小量，忽略对圆柱壳应力分布的影响，壳体失稳前应力和应变仍然为式（1）和式（2），而几何方程如下：

因此相对于完善壳体，式（13）～式（15）分别表述为：

其中：α=12R2，R=r/t。

根据壳体稳定条件得到系数行列式的系数为：

得含轴对称缺陷的轴压圆柱壳失稳临界载荷：

3.2 轴对称缺陷对失稳载荷的影响

为便于分析初始缺陷对薄膜刚度和弯曲刚度的影响，令

式（31）可写为：

图8 缺陷在轴压圆柱壳失稳的影响Fig.8 Effect of defects on the instability of cylindrical shells under axial compression

由图8 中可以看出，a2对圆柱壳的影响很小，可以忽略不计，因此将式（34）简化为：

从式（34）可以看出，轴对称缺陷是经过圆柱壳失稳时轴向波数二次方放大，所以很小的初始缺陷将会对圆柱壳的稳定承载能力产生较大影响。

由图9 可以看出，与无缺陷时轴向受压圆柱壳的失稳模式不同，当纵向失稳波数是1 时，含轴对称初始缺陷的轴压受压圆柱壳的失稳载荷有最小值。比较式（34）和式（25）可得含轴对称缺陷圆柱壳承受轴向均匀压力时，薄膜刚度折减系数与轴对称初始缺陷的关系为：

图9 含轴对称缺陷轴压圆柱壳的失稳临界载荷图谱Fig.9 Instability critical load spectrum of axially compressed cylindrical shells with axisymmetric defects

4 结 语

采用RSM 方法针对轴压载荷作用下完善圆柱壳以及具有轴对称缺陷圆柱壳结构的稳定性进行了分析，主要结论如下：

1）轴压作用下完善圆柱壳结构的失稳载荷与半径、厚度的比有关，与长度、半径比无关。

2）轴压圆柱壳结构的薄膜刚度、弯曲刚度均影响其失稳模态，失稳周向波数随着薄膜刚度的减少而减少，随着弯曲刚度的降低而增加。

3）由于轴对称缺陷通过圆柱壳失稳的轴向波数二次方放大来影响薄膜刚度，故对其临界承载能力影响很大。