圆柱壳方程的奇性解.docxVIP

圆柱壳方程的奇性解 图1中所示的圆形壳体广泛应用于各种设计结构。对于t-r较低的情况，可以通过薄壳理论获得分析公式。圆柱薄壳的定解方程是用中面法向位移un或者切面内力的Airy应力函数ue001φ表示的八阶偏微分方程, 由于壳中几何量与力学量之间满足静力-几何比拟关系,un与ue001φ满足的方程形式相同, 因此这两个方程可以化为复位移-应力函数χ(ξ?φ)=un+i4μ2EΤRφ表示的四阶偏微分方程 (?4+?2-4μ2iο2οξ2)χ=16μ4[R2qnEΤ-?α?2(11-νRΘ0+iR24μ2ΤΘ1)].(1) 式中:ξ=x/R,φ是圆柱壳的量纲1主坐标系, 4μ2=[12 (1-ν2) ]1/2R/T,qn为法向分布面力,?α为热膨胀系数。假设圆柱壳的温度Θ(ξ,φ,ζ) 沿厚度线性分布, 即表示为Θ0(ξ,φ) +ζΘ1(ξ,φ) 。 方程 (1) 称为Morley方程, 它具有薄壳理论本身的精度O(T/R) 。 通常将其一般解分解为一个特解叠加对应的齐次方程通解, 由于方程左端项的四阶偏微分算子无法分解, 在量纲1极坐标系(α=√ξ2+φ2?β=arctan(φ/ξ))中很难求得通解。当α(或φ) 很小时, 方程 (1) 可简化为Donnell扁壳方程, 其齐次形式为 (?4-4μ2iο2οξ2)χ=(?2+2μ√iοοξ)(?2-2μ√iοοξ)χ=0.(2) 式