圆柱壳水下振动的截断模态规律性研究

2018-12-21李永清张焱冰

舰船科学技术 2018年12期

仝博，李永清，朱锡，张焱冰

(海军工程大学舰船工程系，湖北武汉 430033)

0 引言

圆柱壳结构是船舶领域被广泛应用的一种工程结构形式，许多研究者对其水下振动已做了大量研究。线性范围内，圆柱壳在外力作用下的振动响应可看作是其振动模态的叠加[1]，大多数值方法均是基于这一思想，比如瑞利-里兹法[2]、传递矩阵法[3]、基于模态叠加的有限元法等。数值计算追求最小的计算成本达到较高的计算精度。在运用模态叠加法研究结构振动响应时，模态阶数的选取决定计算的经济性和准确性。选取的模态阶数越高，计算结果越准确，但消耗的计算时长却越高。当模态数量大于一定值时计算结果变化不大，趋于稳定，把该值定义为截断模态数。对截

断模态数的研究通常是在理论方法基础上对数值计算的相关参数的研究。比如王宇等[4]采用Love壳体理论分析了薄壁短圆柱壳在3种边界条件下的振动位移响应，指出截取前八阶模态就可达到相当精度。王献忠等[5]结合精细积分法和传递矩阵法，提出精细传递矩阵法法，在对方法有效性验证的同时开展了模态收敛性分析。李榆银等[6]基于NASTRAN软件对薄壁圆柱壳进行强迫振动分析时对截断模态数进行了灵敏度研究，得到了有效的截断波数。

也有学者专门研究了结构的模态截断对结果的影响。李兴泉等[7]基于有限元法对频率截断和有效模态质量截断这2种截断方法进行对比研究。张淼等[8]针对重频阻尼系统提出了高精度截模态算法，并对该算法进行了模态灵敏度和误差分析。金国光等[9]对高速凸轮结构进行动力学分析的同时进行了模态截断研究，基于有限元法研究了不同截取阶数对计算精度和计算速度的影响。也有学者未采用截断模态数确定结构湿表面有限元网格尺度，而是以主模态分量波长作为参考尺度确定的网格尺度划分原则[10]。但对于宽频的振动响应计算，主模态变化较大，仍要以最高频率确定网格尺度。以上对截断模态的研究评判标准大多为单个节点的位移响应或声压响应，不具有代表性，且研究对象为单一圆柱壳，不具有普适性。

本文基于瑞利-里兹法求解圆柱壳振动方程，研究了35个不同尺度的圆柱壳在环频率以下的水下振动响应，研究分析了截断模态数随长径比和环频率的变化关系，并且对其规律性进行总结，通过有限元仿真验证了其正确性，结论对于圆柱壳水下振动问题的研究具有重要的参考意义。

1 圆柱壳水下振动方程推导

研究对象为水下有限长薄壁圆柱壳，长度为L，半径R，厚度为h，h/R≤0.01，材料密度为 ρs，泊松比为μ，模量为E(1+η)， η为损耗因子，本文取0.01。壳体完全浸没于水中，两端简支在无限长刚性圆柱障板上，如图1所示。本文基于Donnell壳体理论研究圆柱壳径向激励下的振动问题，其运动方程为：

其中：位移系数矩阵各项表达式具体形式可查询文献[12]；u，v，w分别为轴向，环向和径向位移；Fr为径向载荷力，作用点位于壳体中部；pr为流体载荷。

图1 有限长圆柱壳模型Fig.1 Finite cylindrical shell model

为了分离圆柱壳的轴向和环向振型，采用分离变量法对轴向和环向位移进行分解，且基于里兹法的思想，位移可假定为一系列模态振型的线性叠加，其公式如下（忽略时间项eiωt）：

其中：m，n分别为轴向模态半波数和环向模态波数，km与边界条件有关，对于本文的对称激励情况，圆柱壳的振动响应仅仅是对称模态的叠加，即轴向半波数m为奇数的情况，因此将k取为 (2m+1)，m=0，1，m2……。

由模态展开法可对激励力和圆柱壳周围的流体载荷进行分解：

将式（2）和式（3）代入式（1）中，可得

式中，为机械阻抗，表达式为

ZMmn

求解表面振速，还需要事先求得声辐射阻抗，求解声辐射阻抗的方法主要有流固耦合方法和声固耦合法[13]2种。本文基于声固耦合法计算Helmholtz方程[14]求解圆柱壳表面声压载荷，表达式为：

将式（3）代入上述方程，依据贝塞尔方程形式求解得声辐射阻抗为：

此时将式（5）和式（7）代入式（4）便可求得圆柱壳表面振速表达式。

圆柱壳表面均方振速公式为：

其中：

将式（9）代入式（8），依据余弦函数的正交性推导可得均方振速为：

表面均方振速级公式：

2 圆柱壳的模态截断

在基于模态叠加法对圆柱壳进行水下振动特性问题研究时发现，不同尺度的圆柱壳，其截断模态数量不同。因此针对不同尺寸的圆柱壳，根据上述圆柱壳振动理论方程，基于Mathematica软件进行编程计算，得到圆柱壳表面均方振速。以表面均方振速级为评判标准，进行轴向半波数m和环向波数n的无关性研究。

2.1 圆柱壳模型

本文以35个不同尺度的圆柱壳为对象，计算它们在环频率以下的水下振动响应，研究轴向截断模态数m、环向截断模态数n随长径比L/R和环频率fr的变化关系。圆柱壳模型参数见表1。其中，环频率公式为：

表1 圆柱壳参数Tab.1 Geometric parameters of cylindrical shells

2.2 截断模态研究

在结构的动态力学特性计算中，为提高计算效率，在不影响计算精度的情况下，往往选取有限个模态数目作为主模态截断，进而通过模态叠加得到结构表面的振动响应。本文目标是研究不同尺度的圆柱壳截断模态数的选取规律，由于计算样本较多，为提高计算效率，将环频率以下的1/3倍频程中心频率作为计算频率点，1/3倍频程中心频率计算式为：

以均方振速级为精度判断指标，具体要求是在环频率以下，模态截断时的计算结果与准确解误差小于1%。理论上，m和n取无穷大时可得到准确解，但在实际计算中难以实现，因此可取一个相对较大的值进行计算作为准确解。以L=1 m，R=1 m圆柱壳为例，经计算分析，当轴向和环向截断模态数达到100时，就能得到准确解。图2和图3分别反映了该圆柱壳轴向和环向截断模态数取不同值时均方振速曲线对比。图中可看出，随着截断模态数的增加，均方振速级曲线与准确解的吻合频率区间不断扩大，当轴向截断模态数m达到8时，当环向截断模态数n达到56时，就与m=n=100时的均方振速级曲线完全吻合，最大误差低于1%，因此可将m=8，n=56作为该圆柱壳的截断模态数。

图2 不同轴向截断模态数时均方振速曲线对比（L=R=1 m）Fig.2 Comparison of mean square vibration velocity curves at different axial truncated modes (L=R=1 m)

图3 不同轴向截断模态数时均方振速曲线对比（L=R=1 m）Fig.3 Comparison of mean square vibration velocity curves at different circumferential truncated modes (L=R=1 m)

对于表1中每个尺寸的圆柱壳都采用上述方法选取截断模态，进行大量计算后将结果进行整理分析，其中专门研究了截断模态数随环频率的变化关系，如图4和图5所示。图中可看出，长径比相同的情况下，随着环频率的增大，轴向和环向截断模态数均呈下降趋势，且下降趋势逐渐平缓。对于轴向截断模态m，长径比越大，截断模态数越高；但对于环向截断模态n，不同长径比圆柱壳其截断模态数随环频率变化曲线趋于重合，说明环向截断模态数对长径比变化不敏感，它仅与环频率有关，也就是与圆柱壳的半径相关。

图4 不同长径比时圆柱壳轴向截断模态数m随环频率变化趋势对比Fig.4 Change trend of the axial truncated modal number with the ring frequency at different ratios of length to diameter

图5 不同长径比圆柱壳环向截断模态数n随环频率变化趋势对比Fig.5 Change trend of the circumferential truncated modal number with the ring frequency at different ratios of length to diameter

为研究圆柱壳振动截断模态在轴向和环向的相对关系，建立关系式研究 ζ与环频率的关系。简单推导可得：，即 ζ代表了轴向波长和r环向波长的比值。图6反映了不同尺度的圆柱壳轴向和环向截断模态波长之比。由图可知，每个圆柱壳截断模态对应的轴向波长与环向波长之比大体相同，均在2～2.5之间，平均值为2.14，即在保证圆柱壳环频率以下的表面均方振速级最大误差不超过1%的情况下，截断模态满足轴向波长大约为环向的2倍。

图6 轴向和环向截断模态波长比值Fig.6 Wavelength ratios of axial and circumferential truncated modes

3 有限元仿真对比分析

上述对截断模态的研究是基于瑞利-里兹法求解圆柱壳振动方程，为验证理论模型的可靠性和截断模态选取规律的正确性，基于有限元方法进行圆柱壳水下振动特性计算，与理论算法进行对比。

算例对象：圆柱壳材料为铝，长L为2 m，半径R为0.8 m，厚度为2 mm，则长径比L/R为2.5，环频率为1 000 Hz。有限元计算采用结构有限元和流体有限元相结合的方法，基于直接稳态动力学法，对模型的原始方程直接积分计算。其中，圆柱壳周围流场域截断方案已在文献[15]进行了探讨，此处直接取柱形流场域半径为5R，如图7所示。

与理论模型不同的是，采用流体有限元法计算圆柱壳水下振动，壳体两端必须进行封闭处理，为了消除两端面对壳体振动的影响，端面除了边缘一圈节点简支，其余所有节点进行固支约束，如图8所示。