复杂边界条件下中厚层合扇形板振动特性分析

2020-08-24庞福振霍瑞东李海超叶开富王雪仁

哈尔滨工程大学学报 2020年7期

庞福振，霍瑞东，李海超，叶开富,3，王雪仁,3

(1.哈尔滨工程大学船舶工程学院，黑龙江哈尔滨 150001； 2.北京航天计量测试技术研究所，北京 100076； 3.中国人民解放军92578部队，北京 100161)

层合扇形板结构在航空、航天、机械制造等领域应用广泛。因此开展层合扇形板的振动特性分析研究，获取层合扇形板结构任意边界条件下的典型特征频率，对指导相应结构在工程中的应用具有重要意义。本文基于一阶剪切变形理论建立数学模型，设板厚方向的剪切应变为常数，通过选取适当的剪切修正系数，即可进行中厚板振动特性分析。

当前已有众多研究复合扇形板振动特性的数值方法，如有限元法[1]、微分求积法[2]、改进三角级数法[3]、广义微分求积法[4]、动力松弛法[5]，微分求积法[6]等，上述方法虽然可以有效求解复合扇形板振动特性问题，但大多只能应用于刚固、简支、自由等几种特定类型的经典边界，对于各种经典边界、弹性边界及其任意组合的复杂边界下的振动特性研究较少。Li[7]提出改进傅里叶级数法被广泛应用于一般边界条件下的板和壳的振动分析中，本文在该方法的基础上进行了拓展。采用弹簧刚度法模拟一般边界条件，采用基于一阶剪切变形理论的Hamilton方程推导运动方程和边界方程。板结构的位移和旋转分量均展开为三角级数，其由标准三角级数和辅助多项式构成。引入辅助多项式函数以消除传统三角级数的位移及其导数在边界处的不连续性并保证其快速收敛。通过将本方法的计算结果与相关文献进行对比，验证了本方法的准确性和可靠性，并给出了任意边界条件下层合扇形板振动特性的结果。

1 层合扇形板自由振动理论推导

1.1 模型描述

层合扇形板由厚度为h、内径为a、外径为b、径向宽度为R、扇形角为φ的多层扇形中厚板层叠而成，其几何尺寸在正交柱坐标系(r,θ,z)中定义，为方便分析还建立了局部坐标系(s,θ,z)，其几何模型由图1(a)给出。本文在各层扇形板边缘处的中性面处施加3组线性弹簧(ku,kv,kw)和2组旋转弹簧(KR,Kθ)，并通过改变边界弹簧的刚度值来模拟不同的边界条件[8]，见图1(b)。其中下标deg_0、deg_1、r_0和r_1分别表示边界弹簧在θ=0、θ=φ、s=0和s=R条件下的弹簧刚度。记第k层材料正交各向异性方向与r轴间的夹角为层压角α，见图1(c)。

图1 中厚层合扇形板几何模型示意Fig.1 The geometric model of laminated sector plates

1.2 位移关系和应力-应变关系

基于一阶剪切变形理论，层合扇形板的位移函数可表示为：

(1)

式中：u0、v0和w0分别表示参考平面上相应点在r、θ和z方向的位移；Ψr和Ψθ分别表示绕r和θ方向且垂直于参考面的旋转分量；t表示时间。层合扇形板的应变分量可表示为：

(2)

根据胡克定律，层合扇形板的应力-应变关系可表示为：

(3)

式中：σr、σθ表示r、θ方向上的法向应力；τrθ、τrz和τθz表示柱坐标系中r、θ和z方向上的剪切应力；Qij(i,j=1,2,4,5,6)表示第k层扇形板的刚度系数。

沿z方向对应力积分，可得力和力矩的积分式：

(4)

式中：Nr、Nθ和Nrθ表示各方向的合力；Mr、Mθ和Mrθ表示各方向的合力矩；Qr、Qθ表示各方向的剪切力；κ为剪切校正因子，本文取κ=5/6[9]。

1.3 扇形板中的能量关系

层合扇形板内的应变能可定义为：

(5)

将式(2)～(4)代入式(5)，可得到用中性面处位移和旋转分量表示的板内应变能。

层合扇形板中的动能可表示为：

(6)

将式(1)中的u、v和w代入式(6),并对z积分可得层合扇形板内部的动能：

(7)

其中：

(8)

存储在边界弹簧中的应变能可表示为：

(9)

1.4 控制方程和边界条件

根据Hamilton原理可得层合扇形板的边界条件和控制方程，其拉格朗日能量泛函L表示为：

L=T-Us-Ub-Ubs-Usp

(10)

由：

(11)

可得扇形板的控制方程：

(12)

层合扇形板的边界条件表示为：

(13)

1.5 位移形函数与方程求解

本文层合扇形板位移和旋转分量均使用三角级数来表示，但传统三角级数在端点处存在不连续现象，因此引入辅助多项式改善其收敛性，辅助函数的引入使位移容许函数能满足求解域内更高阶导数的存在，从而有效地克服边界上可能存在的不连续问题。

扇形板位移函数的三角级数展开式表示为：

(14)

从式(14)中可以看出，位移函数分量中最高存在二阶导数，这意味着至少需要允许函数的二阶导数连续才能满足板上任意点的连续性。为此本文选取以下简单的辅助函数，该函数满足求解域内二阶导数连续且三阶导数存在：

(15)

易证得：

(16)

(K-ω2M)E=0

(17)

式中:K、M和E分别表示扇形板的刚度矩阵、质量矩阵和未知系数矩阵，对上式求解可得到层合扇形板的固有频率，进而可求得相应的振型。

2 数值结果和讨论

2.1 收敛性研究

2.1.1 边界弹簧刚度的影响

本文对扇形板的边界使用弹簧来模拟边界约束，设定恰当的弹簧刚度值可模拟扇形板任意的边界条件。通过讨论弹性边界约束参数对层合扇形板无量纲频率参数Ω=ωb2/h(ρm/E2)1/2的影响规律，探讨了不同边界下合适的弹簧刚度取值，进而验证了弹簧刚度法在此类问题应用中的可行性。本算例中模型的几何和材料属性参数为：a/b=0.5，h/b=0.1，E1/E2=15，E2=1×1010Pa，G12=G13=0.6E2，G23=0.5E2，μ12=0.25，ρ=1 500 kg/m3。为简化分析过程，取弹性约束参数Гu、Гv、Гw、ГR和Гθ表示相应弹簧刚度与抗弯刚度D=E2h3/12(1-μ122)的比值，即Гu=ku/D、Гv=kv/D、Гw=kw/D、ГR=KR/D、Гθ=Kθ/D。此外，定义参数ΔΩ=ΩΓλ-ΩΓλ=10-8，用于表示无量纲频率Ω与参考值的差值，参考值取ΩΓλ=10-8。

图2给出了非对称型双层层合扇形板[0°/90°]前三阶频率参数ΔΩ随弹性约束参数Гλ(λ=u,v,w,R,θ)的变化关系。本算例所取边界条件为：s=0、s=R处自由，θ=0处刚固，θ=φ处弹性约束。从图2可以看出，弹性约束参数Гλ小于10-1时，ΔΩ几乎为零，此时可视自由边界条件；弹性约束参数Гλ在101～106范围内时，ΔΩ与Гλ呈正相关关系，此时可视为弹性边界条件；弹性约束参数Гλ在107以上后ΔΩ基本趋于稳定，此时扇形板的无量纲频率求解已达到收敛，该范围可视为刚固边界条件。

图2 双层非对称层合扇形板[0°/90°]的频率参数ΔΩ随弹性约束参数Γλ的变化关系Fig.2 The frequency parameters ΔΩ versus the elastic restraint parameters Γλ for laminated sector plate[0°/90°]

通过大量分析计算，给出了本算例所给工况下几种经典边界条件对应各弹簧刚度的合理取值(如表1)。

表1 几种边界条件下各弹簧刚度值的合理取值Table 1 The suitable values of spring stiffness value under several boundary conditions

为了检验以上弹簧刚度取值的合理性，将以弹簧模拟不同边界来进行振动特性分析，并与相关文献进行数值对比。为简洁起见，使用字符串来表示边界条件，例如FCSC表示在扇形板s=0、θ=0、s=R和θ=φ处具有F、C、S和C边界条件。除非另有说明，本文层合扇形板无量纲频率参数表示为：Ω=ωb2/h(ρm/E2)1/2，扇形板的材料和几何参数如下：E2=1×1010Pa，E1/E2=15，G12=G13=0.6E2，G23=0.5E2，μ12=0.25，ρ=1 500 kg/m3。

2.1.2 截断数的影响

本节探讨三角级数截断数对层合扇形板无量纲频率收敛性的影响。扇形板的位移函数由三角级数来表示，显然其位移和旋转分量上三角级数的展开阶数M和N对计算收敛性有直接关系，理论上展开阶数越多计算越精确，但相应的计算效率会大幅降低，因此选取合适的展开阶数具有重要的探讨价值。

表2给出了CCCC边界条件下双层层合扇形板[0°/90°]前六阶固有频率随截断数影响的结果。从表中可以看出本文方法收敛性良好，M=N=8和M=N=16的结果之间的最大差异小于0.077%。因此如不做特殊说明，本文方法在后续计算中三角级数的截断数统一取M=N=12。

表2 CCCC边界下双层[0°/90°]层合扇形板前六阶频率参数收敛性(φ=2π/3, a/b=0.5, h/b=0.2, E1/E2=40)Table 2 Convergence of the first six frequency parameters for a laminated sector plate [0°/90°]with CCCC boundary conditions

2.2 有效性验证及特性分析

表3给出了2种方法对不同经典边界条件下多层层合扇形板[0°/90°]5的频率参数计算结果的比较(φ=π/3,a/b=0.1,h/b=0.2,E1/E2=40)。从表3可知本方法与参考文献相比最大相对误差不超过0.634%，由此可证得本方法的有效性。

表3 不同边界下层合扇形板[0°/90°]5Ω值与文献[10]对比Table 3 Comparison of the Ω for a laminated sector plate [0°/90°]5 with different boundary conditions

表4～6给出了刚固、简支、自由等边界条件下层合扇形板无量纲频率与厚径比、正交各向异性比、层压角等结构参数间的关系。表4表明无量纲频率随扇形板厚径比的增加而降低(φ=π/2,a/b=0.5,E1/E2=15)。表5表明频率参数随正交各向异性比的增加而增加(φ=π/2,a/b=0.5,h/b=0.1)。表6表明扇形板的频率参数随层压角的增大而减小(E1/E2=15,a/b=0.5,h/b=0.1)。图3给出了CCCC和CFCF边界条件下层合板前三阶阵型图(E1/E2=15,a/b=0.5,h/b=0.1,φ=2π/3)。由图3可知，同阶阵型下，边界条件的不同会导致其模态特性不尽相同。

表4 不同边界下层合扇形板[0°/90°]2Ω值与厚径比的关系Table 4 The Ω for a laminated sector plate[0°/90°]2 with thickness-to-radius ratios in different boundary conditions

表5 不同边界下层合扇形板[0°/90°]2 Ω值与各向异性比关系Table 5 The Ω for a laminated sector plate[0°/90°]2 with orthotropy ratios in different boundary conditions

表6 不同边界下扇形板[0°/90°/90°/0°] Ω值与层压角的关系Table 6 The Ω for a laminated sector plate[0°/90°/90°/0°] with included angles in different boundary conditions

层合扇形板在实际工程应用中通常会有更为复杂的一般边界条件约束，表7和表8给出了在一般弹性边界约束下层合扇形板的前四阶固有频率。本算例中扇形板的几何尺寸除φ=π外均与图3相同。从表7和表8可以看出，不同的层合方案对层合扇形板的频率也有明显影响。图4给出了CE3CE3和E3E3E3E3边界条件下层合扇形板[0°/90°/90°/0°]前三阶模态的阵型图，由图4可知，边界条件的不同使得同阶阵型层合板振动模态呈现明显差异。