Prediction on hull vibration and underwater noise radiation based on modified LMM
-
摘要:目的 大质量法(LMM)可将设备机脚的激励加速度转化为力载荷,是结构振动响应与声辐射仿真分析的常用方法。为提高LMM在船体声振响应计算方面的应用精度,提出对LMM的修正。方法 首先,分析LMM的适用原理,得出合理的大质量倍数区间; 然后,基于船舶的阻尼分析,引入瑞利阻尼质量比例系数对LMM加以修正; 最后,将修正前、后的LMM计算结果分别与某油船测试的数据进行对比,并基于修正大质量法(M-LMM)计算所得船体谐响应结果,进一步计算船舶水下辐射声功率与垂直声场指向。结果 结果显示,在20~200 Hz分析频段内,应用2种LMM计算所得的振动响应频谱曲线趋势相近,但M-LMM的计算结果更接近于实测值,平均误差仅2.70%。结论 M-LMM可提高船体响应仿真精度,能为船舶前期设计阶段提供更准确的设计依据。Abstract:Objectives Large Mass Method (LMM) is capable to approximately convert excitation acceleration of equipments' seat into force loads, which is a common method for structural vibration response and acoustic radiation simulation analysis. In order to improve the accuracy of LMM applied in ship's noise & vibration response computation, a modified LMM(M-LMM)is presented.Methods Its application principles in hull structure response were analyzed, based on which the appropriate large mass ratio ensuring the application precision of LMM was obtained. Rayleigh damping was then taken into consideration, and LMM was modified with ratio coefficient. The experiment comparing results from original LMM as well as M-LMM with tested ones, together with the computations of underwater noise radiation power and vertical directivity pattern of sound field, were followed.Results For the frequency range of 20 Hz-200 Hz, the analysis results got from M-LMM were closer to measured values compared to that from LMM although their computed vibration response spectrums were roughly similar.Conclusions It is crystal clear that M-LMM can guarantee a simulation of main hull response with higher accuracy(average deviation 2.7%), providing the pre-design stage of vessels with more necessary information.
-
0. 引言
船体声振响应预测关系到船舶的舒适性及隐身性, 是舰船研究的一个重要领域。通过提高声振预测精度, 设计人员可在前期有针对性地对结构设计方案进行优化, 从而提升船舶建造的整体效率。一些主流的商业仿真软件, 例如ANSYS, Nastran等可快速、直观地预报相应船体及舱段结构的声振响应, 因此得到了广泛应用。但运用这些仿真软件进行声振预测存在一个问题, 即是在进行主激励设备振动的台架试验时, 一般只能获得机脚传至基础表面的振动加速度信号, 无法直接获得机脚激振力, 而主流仿真软件所能施加的载荷则为激振力。针对此问题, 目前解决的方法主要有惯性载荷法、相对运动法、拉格朗日乘子法[1]、大质量法(LMM)等。其中, LMM适用于结构多点激励分析, 且能克服相对运动法求解大型复杂结构的不便, 故备受青睐。LMM最初主要用于模拟地震激励对地面物体运动特性的影响[2], 而后在汽车、船舶等结构分析中得到了应用[3-4]。
不过, 虽然LMM操作简单、使用方便, 但由于缺少严密的数学推理和清晰的物理概念, 其求解精度和适用性在很长一段时间备受争议[5]。为此, Paraskevopoulos等[6]研究了LMM在特定情况下的适用性; Kim和Jhung[7]针对单自由度和多自由度模型, 对LMM进行了数理分析并给出了详尽的数学推导过程。在LMM求解精度方面, 周国良等[8]研究指出, 在考虑系统阻尼时, 应用LMM会产生附加阻尼力并导致误差; 许玲玲等[9]在此基础上, 在基于LMM解决离散元法的多点激励问题时引入了质量比例阻尼系数, 减小了计算误差; 黄炳臣等[10]将此改进思路运用到核电反应堆冷却系统的时程计算中, 补偿了大质量点因瑞利阻尼产生的运动阻力, 从而提升了计算精度。在船舶领域, 王迪等[11]基于LMM, 计算了某船舶舱段的水下辐射噪声; Zhang等[12]则以某油船为研究对象, 应用LMM计算了船体结构振动响应, 结果显示仿真结果与实测值较为吻合。
以上针对船舶声振响应计算的文献中, 均未提出对于LMM的修正。为此, 本文将在参考文献[8-10]的思路的基础上, 首先将其应用到船舶声振计算领域, 引入质量比例阻尼系数, 对LMM的等效加速度予以修正, 以进一步提升仿真精度; 然后分别应用修正及未修正的LMM进行全船振动响应计算, 通过与实测振动响应数值的对比, 证明修正大质量法(M-LMM)的精准度更优; 最后, 以全船谐响应计算结果作为边界条件计算外场辐射噪声, 进一步讨论2种LMM计算得到的辐射噪声结果。
1. 基本原理
1.1 船体结构响应的大质量模型
如图 1所示, 在设备机脚激励引起的船体结构响应分析中, 激励由设备结构传向船体, 因此将设备设为大质量。
建立大质量单自由度简化模型, 设备大质量M与质量块m建立刚性连接, 对于整个由M-m构成的振动系统, 将大质量结构与基底的振动微分方程联立表示为矩阵形式, 有
\left[ {\begin{array}{*{20}{c}} M&0\\ 0&m \end{array}} \right]\left\{ {\begin{array}{*{20}{l}} {{{\mathit{\boldsymbol{\ddot x}}}_{\rm{s}}}}\\ {{{\mathit{\boldsymbol{\ddot x}}}_{\rm{b}}}} \end{array}} \right\} + \left[ {\begin{array}{*{20}{c}} k&{ - k}\\ { - k}&k \end{array}} \right]\left\{ {\begin{array}{*{20}{l}} {{\mathit{\boldsymbol{x}}_{\rm{s}}}}\\ {{\mathit{\boldsymbol{x}}_{\rm{b}}}} \end{array}} \right\} = \left\{ {\begin{array}{*{20}{c}} 0\\ {M\ddot z} \end{array}} \right\} (1) 式中:k为质量块m与大质量M间的连接刚度; xs为大质量M的位移; xb为质量块m的位移; \ddot z为大质量M的加速度; M\ddot z=F, 为加诸于基底的激振力。
参考式(1)的矩阵形式, 多自由度大质量模型的总振动微分方程可以写为
\begin{array}{*{20}{c}} {\left[ {\begin{array}{*{20}{c}} M&0&0& \cdots &0&0\\ 0&{{m_1}}&0& \cdots &0&0\\ 0&0&{{m_2}}& \cdots &0&0\\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots \\ 0&0&0& \cdots &{{m_{n - 1}}}&0\\ 0&0&0& \cdots &0&{{m_n}} \end{array}} \right]\left\{ {\begin{array}{*{20}{c}} {{{\mathit{\boldsymbol{\ddot x}}}_{\rm{s}}}}\\ {{{\mathit{\boldsymbol{\ddot x}}}_{{\rm{b}}1}}}\\ {{{\mathit{\boldsymbol{\ddot x}}}_{{\rm{b}}2}}}\\ \vdots \\ {{{\mathit{\boldsymbol{\ddot x}}}_{{\rm{b}} \cdot n - 1}}}\\ {{{\mathit{\boldsymbol{\ddot x}}}_{{\rm{b}} \cdot n}}} \end{array}} \right\} + }\\ {\left[ {\begin{array}{*{20}{c}} {{k_1}}&{ - {k_1}}&0& \cdots &0&0\\ { - {k_1}}&{{k_1} + {k_2}}&{ - {k_2}}& \cdots &0&0\\ 0&{ - {k_2}}&{{k_2} + {k_3}}& \cdots &0&0\\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots \\ 0&0&0& \cdots &{{k_{n - 1}} + {k_n}}&{ - {k_n}}\\ 0&0&0& \cdots &{ - {k_n}}&{{k_n}} \end{array}} \right] \cdot }\\ {\left\{ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{x}}_{\rm{s}}}}\\ {{\mathit{\boldsymbol{x}}_{{\rm{b}}1}}}\\ {{\mathit{\boldsymbol{x}}_{{\rm{b}}2}}}\\ \vdots \\ {{\mathit{\boldsymbol{x}}_{{\rm{b}} \cdot n - 1}}}\\ {{\mathit{\boldsymbol{x}}_{{\rm{b}} \cdot n}}} \end{array}} \right\} = \left\{ {\begin{array}{*{20}{c}} 0\\ {M\ddot z}\\ 0\\ \vdots \\ 0\\ 0 \end{array}} \right\}} \end{array} (2) 将式(2)包含的所有方程累加, 可得
M{{\mathit{\boldsymbol{\ddot x}}}_{\rm{s}}} + \sum\limits_{i = 1}^n {{m_i}{{\mathit{\boldsymbol{\ddot x}}}_{{\rm{b}}i}}} = M\ddot z (3) 式中, i为各质量块编号。式(3)可写为
{{\mathit{\boldsymbol{\ddot x}}}_{\rm{s}}} + \frac{{\ddot w}}{\lambda } = \ddot z (4) 式中:\ddot w = \sum\limits_{i = 1}^n {{m_i}} {\mathit{\boldsymbol{\ddot x}}_i}/\sum\limits_{i = 1}^n {{m_i}} ;\lambda = M/\sum\limits_{i = 1}^n {{m_i}} , 为大质量和船体结构质量总和的比值, 定义为“大质量比”。可见, 当\ddot w = \sum\limits_{i = 1}^n {{m_i}} {\mathit{\boldsymbol{\ddot x}}_i}/\sum\limits_{i = 1}^n {{m_i}} 不变时, 随着质量比λ不断的增加, 有{\mathit{\boldsymbol{\ddot x}}_{\rm{s}}} \to \ddot z。因为船体结构的总质量保持不变, 故可通过增加设备质量M来实现激励加速度值和基底表面加速度值的近似。
1.2 引入阻尼项的LMM修正
在上节讨论的基础上加入对阻尼的考虑。将式(2)写成矩阵形式, 并加入阻尼项, 有
\mathit{\boldsymbol{M\ddot x}} + \mathit{\boldsymbol{C\dot x}} + \mathit{\boldsymbol{Kx}} = \mathit{\boldsymbol{F}} (5) 式中:C为阻尼矩阵; K为刚度矩阵; M为质量矩阵, \mathit{\boldsymbol{M}} = \left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{M}}_{\rm{b}}}}&0\\ 0&{{\mathit{\boldsymbol{M}}_{\rm{s}}}} \end{array}} \right], 其中Ms为响应输入对象的分块质量矩阵, 即船体的分块质量矩阵, Mb为设备大质量点的分块质量矩阵; F为激振力矩阵。按分块矩阵形式展开, 有
\begin{array}{*{20}{c}} {\left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{M}}_{\rm{b}}}}&0\\ 0&{{\mathit{\boldsymbol{M}}_{\rm{s}}}} \end{array}} \right]\left\{ {\begin{array}{*{20}{l}} {{{\mathit{\boldsymbol{\ddot x}}}_{\rm{b}}}}\\ {{{\mathit{\boldsymbol{\ddot x}}}_{\rm{s}}}} \end{array}} \right\} + \left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{C}}_{11}}}&{{\mathit{\boldsymbol{C}}_{12}}}\\ {{\mathit{\boldsymbol{C}}_{21}}}&{{\mathit{\boldsymbol{C}}_{22}}} \end{array}} \right]\left\{ {\begin{array}{*{20}{l}} {{{\mathit{\boldsymbol{\dot x}}}_{\rm{b}}}}\\ {{{\mathit{\boldsymbol{\dot x}}}_{\rm{s}}}} \end{array}} \right\} + }\\ {\left[ {\begin{array}{*{20}{l}} {{\mathit{\boldsymbol{K}}_{11}}}&{{\mathit{\boldsymbol{K}}_{12}}}\\ {{\mathit{\boldsymbol{K}}_{21}}}&{{\mathit{\boldsymbol{K}}_{22}}} \end{array}} \right]\left\{ {\begin{array}{*{20}{l}} {{{\mathit{\boldsymbol{\ddot x}}}_{\rm{b}}}}\\ {{{\mathit{\boldsymbol{\ddot x}}}_{\rm{s}}}} \end{array}} \right\} = \left\{ {\begin{array}{*{20}{c}} 0\\ {{\mathit{\boldsymbol{M}}_{\rm{b}}}\ddot z} \end{array}} \right\}} \end{array} (6) 式中:C11, C12(C21), C22分别为大质量阻尼分块矩阵、船体与大质量的耦合阻尼分块矩阵及船体阻尼分块矩阵; K11, K12, K21, K22为相应的刚度分块矩阵。将式(6)中包含的矩阵方程展开并累加, 有
\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{M}}_{\rm{b}}}{{\mathit{\boldsymbol{\ddot x}}}_{\rm{b}}} + {\mathit{\boldsymbol{M}}_{\rm{s}}}{{\mathit{\boldsymbol{\ddot x}}}_{\rm{s}}} + {\mathit{\boldsymbol{C}}_{11}}{{\mathit{\boldsymbol{\ddot x}}}_{\rm{b}}} + {\mathit{\boldsymbol{C}}_{12}}{{\mathit{\boldsymbol{\dot x}}}_{\rm{s}}} + {\mathit{\boldsymbol{C}}_{21}}{{\mathit{\boldsymbol{\dot x}}}_{\rm{b}}} + {\mathit{\boldsymbol{C}}_{22}}{{\mathit{\boldsymbol{\dot x}}}_{\rm{s}}} + }\\ {{\mathit{\boldsymbol{K}}_{11}}{\mathit{\boldsymbol{x}}_{\rm{b}}} + {\mathit{\boldsymbol{K}}_{12}}{\mathit{\boldsymbol{x}}_{\rm{s}}} + {\mathit{\boldsymbol{K}}_{21}}{\mathit{\boldsymbol{x}}_{\rm{b}}} + {\mathit{\boldsymbol{K}}_{22}}{\mathit{\boldsymbol{x}}_{\rm{s}}} = {\mathit{\boldsymbol{M}}_{\rm{b}}}\ddot z} \end{array} (7) 对式(7)等式两边左乘Mb-1, 可见Mb-1的对角线元素1/mbj→0, 式(7)左侧除{\mathit{\boldsymbol{M}}_{\rm{b}}}{\mathit{\boldsymbol{\ddot x}}_{\rm{b}}}项外, 其余项左乘Mb-1后趋近为零矩阵。然而, 当考虑瑞利阻尼时,
\mathit{\boldsymbol{C}} = \left[ {\begin{array}{*{20}{l}} {{\mathit{\boldsymbol{C}}_{11}}}&{{\mathit{\boldsymbol{C}}_{12}}}\\ {{\mathit{\boldsymbol{C}}_{21}}}&{{\mathit{\boldsymbol{C}}_{22}}} \end{array}} \right] = \alpha \left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{M}}_{\rm{b}}}}&0\\ 0&{{\mathit{\boldsymbol{M}}_{\rm{s}}}} \end{array}} \right] + \beta \left[ {\begin{array}{*{20}{l}} {{\mathit{\boldsymbol{K}}_{11}}}&{{\mathit{\boldsymbol{K}}_{12}}}\\ {{\mathit{\boldsymbol{K}}_{21}}}&{{\mathit{\boldsymbol{K}}_{22}}} \end{array}} \right] (8) 式中, α和β为瑞利阻尼系数。将式(8)代入式(7), 等式两边左乘Mb-1后有
{{\mathit{\boldsymbol{\ddot x}}}_{\rm{b}}} + \alpha {{\mathit{\boldsymbol{\ddot x}}}_{\rm{b}}} \approx \ddot z (9) α并非无穷小量, 故将\alpha {\mathit{\boldsymbol{\dot x}}_{\rm{b}}}作为干扰项, 会影响LMM计算结果的精度。因此, 需增加补偿力α\ddot z加以修正。质量相关系数α可参照文献[13]进行计算:
\left\{ {\begin{array}{*{20}{l}} {\xi = \frac{\alpha }{{2{\omega _{\rm{a}}}}} + \frac{{\beta {\omega _{\rm{a}}}}}{2}}\\ {\xi = \frac{\alpha }{{2{\omega _{\rm{b}}}}} + \frac{{\beta {\omega _{\rm{b}}}}}{2}} \end{array}} \right. (10) 式中:ξ为阻尼比, 一般在0.02~0.05内取值, 本文取ξ=0.05[14]; ωa与ωb均为结构特征频率, 前者取结构基频, 后者在对动力响应有较大贡献的高阶振型频率中选取[15]。以某近海油船(图 2(a))为研究对象, 研究工况为空载返航, 此时开启的主要动力设备为1台主机和1台柴油发电机组。根据船体设计图纸, 建立了考虑流体耦合的全船有限元模型, 并在油船双层底处建立了主机(带齿轮箱)和柴油发电机组的基座, 如图 2(b)所示。以大质量单元模拟设备, 大质量设为船体总质量的107倍, 随后, 将大质量单元与机脚建立刚性域连接。
提取主机基座与柴油发电机组基座的有限元模型, 在对基座底部进行固支约束后再分别进行振型计算。经过计算, 对柴油发电机组的基座取ωa=94.18 rad/s, ωb=447.13 rad/s(第8阶固有频率); 对主机及齿轮箱基座取ωa=38.89 rad/s, ωb=101.73rad/s(第25阶固有频率)。由式(10)计算得到柴油发电机组的瑞利阻尼系数α1=7.78, 主机及齿轮箱的瑞利阻尼系数α2=2.81。
为准确计算设备机脚载荷力, 还必须获取机脚加速度及速度。鉴于所研究油船的运行工况, 将主机、柴油发电机组以及齿轮箱作为模拟的激励源设备, 其机脚加速度通过实船测取。将各激励设备的机脚加速度实测值在时域上取平均, 用以表示各设备机脚的实际加速度值。基于MATLAB平台写入机脚加速度时域信号, 调用cumtrapz函数进行积分, 得到相应的速度时域信号如图 3所示。
按\ddot z = \ddot x + \alpha \dot x计算得到修正的设备机脚加速度, 经傅里叶变换, 得到20~200 Hz频段内的加速度频谱如图 4所示。由图可见, 修正后, 各设备的机脚振动加速度值在激励峰值频率点处改变较小, 但在其余频率点处其加速度数值均有不同程度的提高。
2. 船体声振计算分析
2.1 船体振动数值分析
对各个设备的大质量单元施加修正后的力载荷频谱, 在20~200 Hz频段内对全船进行谐响应计算, 并提取主要测点的振动数据, 用于与实船实测的振动数据进行对比。
测点选择在可能出现较大振幅的位置。参考《船上振动控制指南》中的测点布置方案, 并考虑测量现场情况, 选择如图 5所示的测点位置。舱内的振动测量则以餐厅、机舱和驾驶室为对象。
为保证仿真计算结果的准确性, 在有限元模型对应测点附近取4个节点, 提取测点及其相邻节点的振动位移数据, 并按下式计算测点的平均振动位移级:
{{\bar L}_{\rm{a}}} = 10\lg \left[ {\frac{1}{N}\sum\limits_{i = 1}^N 1 {0^{{L_{{\rm{a}}i}}/10}}} \right] (11) 式中, Lai为各节点的振动位移级。将计算得到的各测点的平均振动位移级{\bar L_{\rm{a}}}转化为振动加速度级。限于篇幅, 本文仅展示了部分测点实测值与仿真值的对比, 结果如图 6所示。
由图可见, 仿真值与实测值的频谱趋势一致性较好, 且各测点的主要振动峰值频率基本吻合, 可以较为准确地反映船体的振动状态。为了更加直观地判断M-LMM的准确性, 将本文基于M-LMM计算得到的船体各部位的振动响应与文献[12]的计算结果进行了比较。后者使用未经修正的LMM对同一艘油船进行了谐响应分析。对比参数为船体主要位置处测点在20~200 Hz(频率间隔2 Hz)内合成后得到的带宽内总体振动加速度级。对比结果如表 1所示。
表 1 各测点的总体振动加速度级Table 1. The overall vibration acceleration levels of measure points测点 总体振动加速度/dB 误差/% 实测值 LMM M-LMM LMM M-LMM 机舱侧壁 105.45 100.52 103.22 -4.67 -2.11 上甲板 83.70 77.76 82.44 -7.10 -1.51 上甲板舷侧 96.31 102.38 103.19 6.30 7.14 餐厅侧壁 97.57 95.71 96.41 -1.91 -1.19 驾驶室地板 88.26 86.69 89.59 -1.78 1.51 由表 1可见, 除上甲板舷侧的测点以外, 其余测点上M-LMM的仿真精度较之LMM均有较大改善。在所有测点中, 机舱侧壁由于靠近激振源, 其总体振动加速度级最大。上甲板舷侧因加筋较弱, 也表现出了较大的振动量, 这可能是由于船体建模误差、测点距离激振源较远、仿真计算偏差等原因所致, M-LMM在精度上的改进在该测点没有得到体现。综合考虑各个测点的仿真精度, LMM的平均误差为4.35%, M-LMM的平均误差则为2.70%, 精度提高明显。
2.2 船体声辐射计算分析
将水线以下船体外壳的振动响应结果作为声学边界条件, 应用间接边界元法进行声辐射计算:将谐响应计算的rst文件导入声学分析软件Virtual.Lab, 设置外场流体属性, 取海水密度ρw=1 025 kg/m3, 声速cw=1 500 m/s, 将船体湿表面的位移响应映射到边界元模型上, 最后进行辐射声功率的计算。基于2种LMM计算得到的声功率频谱曲线如图 7所示。
在170~200 Hz频段内, M-LMM得到的声功率辐射峰值略大于LMM所得。对声功率进行带宽内合成, 基于M-LMM的总体辐射声功率为140.547 dB, 基于LMM的则为139.021 dB。因载荷力修正所引起的船体响应数值的提高, 在声辐射的计算上也得到了体现。
船舶作为结构复杂的噪声源, 其外场噪声辐射有着很明显的方向性, 其中垂直声场指向性可用于描述位于垂直于船体纵轴线的横剖面内噪声辐射的方向性。仍以上述油船作为对象, 对其水下声场垂直指向性进行研究, 分析模型如图 8所示。
图 9设置了3个垂直于船体纵向轴线的横剖面, 分别用于显示油船艉部、舯部及艏部的声场指向性。以船体中心点所在纵轴线为中心, 分别在上述3个横剖面内设置半径为50 m的平面场点, 计算20~200 Hz频段内船体的总声压级, 得到相应场点上的声压指向性曲线如图 9所示。
由图可见, 油船的总体声压以艉部最大, 艏部最小, 舯部居间。本文所讨论的激励设备均布置在艉部机舱, 因此, 可以认为艉部的水下辐射声压主要受动力设备激励的影响, 而舯部和艏部因离振动激励源较远, 故各方向的噪声辐射水平均有所下降。除个别方向以外, 采用M-LMM计算所得的场点声压级均大于LMM的计算所得。随着截取的横剖面逐渐远离船艉激励源, 这种数值差异方面的不对称将逐渐消失, 2种LMM计算所得的垂直指向性曲线也趋于吻合。
3. 结论
本文基于理论推导分析了使用LMM时会产生的误差, 通过引入质量相关系数, 对LMM进行了修正。通过和实船实测的振动数据进行对比, 对M-LMM计算振动响应的准确性进行了验证, 同时基于谐响应计算结果, 进一步计算了船舶的水下辐射声功率与垂直指向性。综上分析, 得到如下主要结论:
1) 较之LMM, M-LMM的仿真结果能更准确地反映船体结构的振动响应水平。在20~200 Hz分频段内, M-LMM的仿真值与实测值的频谱曲线较为吻合; 对于总体振动加速度级这一评判数值, M-LMM的平均误差为2.70%, 较之LMM的平均误差4.35%有较大的提高。
2) M-LMM与LMM计算所得的辐射功率曲线相近, 前者的总体辐射声功率数值略高于后者(1.526 dB); 在以船体纵轴线为中心的各个方向上, M-LMM计算所得的场点声压级基本高于LMM的计算所得; 当远离艉部激励源时, 二者在声压级数值方面上差距逐渐减小, 指向性曲线趋于一致。
本文在进行声辐射计算时采用了有限元—边界元(FE-BEM)法, 由于M-LMM计算得到的船体振动响应结果更接近于实测值, 能够为声辐射计算提供更为精确的边界条件, 故理论上和LMM相比能更准确地模拟油船声辐射水平。但限于设备条件, 未能对声辐射计算的准确性进行验证, 这是本文今后有待拓展的方面。
-
表 1 各测点的总体振动加速度级
Table 1 The overall vibration acceleration levels of measure points
测点 总体振动加速度/dB 误差/% 实测值 LMM M-LMM LMM M-LMM 机舱侧壁 105.45 100.52 103.22 -4.67 -2.11 上甲板 83.70 77.76 82.44 -7.10 -1.51 上甲板舷侧 96.31 102.38 103.19 6.30 7.14 餐厅侧壁 97.57 95.71 96.41 -1.91 -1.19 驾驶室地板 88.26 86.69 89.59 -1.78 1.51 -
[1] 王樑.外激励下星-箭界面载荷的分析[D].长沙: 湖南大学, 2016. http://cdmd.cnki.com.cn/Article/CDMD-10532-1016252004.htm Wang L. Analysis of interface force between satellite and launch vehicle under external excitation[D]. Changsha: Hunan University, 2016(in Chinese). http://cdmd.cnki.com.cn/Article/CDMD-10532-1016252004.htm
[2] Léger P, Idé I M, Paultre P. Multiple-support seismic analysis of large structures[J]. Computers and Structures, 1990, 36(6):1153-1158. doi: 10.1016/0045-7949(90)90224-P
[3] 王有智, 夏卫华.用大质量法计算东风EQ140汽车前面罩的强迫运动响应[J].力学与实践, 1998, 20(5):32-34. http://www.cnki.com.cn/Article/CJFDTotal-LXYS805.009.htm Wang Y Z, Xia W H. The calculation of enforced motion response of the front face of Dongfeng EQ140 truck by the large mass method[J]. Mechanics and Engineering, 1998, 20(5):32-34(in Chinese). http://www.cnki.com.cn/Article/CJFDTotal-LXYS805.009.htm
[4] 李增光, 马忠俊.多支点处不同冲击作用下的舰船推进轴系冲击计算及分析[J].武汉理工大学学报(交通科学与工程版), 2013, 37(3):566-569. doi: 10.3963/j.issn.2095-3844.2013.03.028 Li Z G, Ma Z J. Anti-shock analysis of the propulsion shafting of warship under different multiple support-points load[J]. Journal of Wuhan University of Technology(Transportation Science & Engineering), 2013, 37(3):566-569(in Chinese). doi: 10.3963/j.issn.2095-3844.2013.03.028
[5] 周国良, 李小军, 刘必灯, 等.大质量法在多点激励分析中的应用、误差分析与改进[J].工程力学, 2011, 28(1):48-54. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gclx201101008 Zhou G L, Li X J, Liu B D, et al. Error analysis and improvements of large mass method used in multi-support seismic excitation analysis[J]. Engineering Mechanics, 2011, 28(1):48-54(in Chinese). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gclx201101008
[6] Paraskevopoulos E A, Panagiotopoulos C G, Manolis G D. Imposition of time-dependent boundary conditions in FEM formulations for elastodynamics:critical assessment of penalty-type methods[J]. Computational Mechanics, 2010, 45(2/3):157-166. http://cn.bing.com/academic/profile?id=a186da64b71c0878e9e44b34693c8a68&encoded=0&v=paper_preview&mkt=zh-cn
[7] Kim Y W, Jhung M J. A study on large mass method for dynamic problem of multiple degree-of-freedom system excited by ground acceleration time history?[J]. Journal of Mechanical Science and Technology, 2014, 28(1):25-41. doi: 10.1007/s12206-013-0941-x
[8] 周国良, 鲍叶欣, 李小军, 等.结构动力分析中多点激励问题的研究综述[J].世界地震工程, 2009, 25(4):25-32. http://d.old.wanfangdata.com.cn/Periodical/sjdzgc200904005 Zhou G L, Bao Y X, Li X J, et al. Review on dynamic analyses of structures under multi-support excitation[J]. World Earthquake Engineering, 2009, 25(4):25-32(in Chinese). http://d.old.wanfangdata.com.cn/Periodical/sjdzgc200904005
[9] 许玲玲, 叶继红, 薛素铎.大质量法在基于离散元法多点激励分析中的应用和误差修正[J].土木工程学报, 2016, 49(增刊1):32-36. Xu L L, Ye J H, Xue S D. Application and error improvement of large mass method used in DEM-based multi-support excitation analysis[J]. China Civil Engineering Journal, 2016, 49(Supp 1):32-36(in Chinese).
[10] 黄炳臣, 沈伟, 卢智.瑞利阻尼对大质量直接积分时程计算方法的影响[J].核动力工程, 2015, 36(2):21-23. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hdlgc201502006 Huang B C, Shen W, Lu Z. Effect of Rayleigh damping on result of large mass direct integration method[J]. Nuclear Power Engineering, 2015, 36(2):21-23(in Chinese). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hdlgc201502006
[11] 王迪, 朱翔, 李天匀, 等.基于虚拟质量法的水面船舶辐射噪声分析[C]//中国造船工程学会船舶力学学术委员会水下噪声学组成立三十周年第十五届船舶水下噪声学术讨论会论文集.郑州: 中国造船工程学会, 2015. Wang D, Zhu X, Li T Y, et al. The analysis on radiation of surface ship based on virtual mass method[C]//The 30th Anniversary for Founding Underwater Accoustic Research Group of Marine Mechanics Acedemic Committee of the Chinese Society of Naval Architecture Engineering & 15th Conference of Ship Unerwater Noise. Zhengzhou: Chinese Society of Shipbuilding Engineering, 2015(in Chinese).
[12] Zhang B, Xiang Y, He P, et al. Study on prediction methods and characteristics of ship underwater radiated noise within full frequency[J]. Ocean Engineering, 2019, 174:61-70. doi: 10.1016/j.oceaneng.2019.01.028
[13] Yoshida N, Kobayashi S, Suetomi I, et al. Equivalent linear method considering frequency dependent characteristics of stiffness and damping[J]. Soil Dynamics and Earthquake Engineering, 2002, 22(3):205-222. doi: 10.1016/S0267-7261(02)00011-8
[14] 胡成宝, 王云岗, 凌道盛.瑞利阻尼物理本质及参数对动力响应的影响[J].浙江大学学报(工学版), 2017, 51(7):1284-1290. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zjdxxb-gx201707003 Hu C B, Wang Y G, Ling D S. Physical essence and influence of model parameters on dynamic response of Rayleigh damping[J]. Journal of Zhejiang University (Engineering Science), 2017, 51(7):1284-1290(in Chinese). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zjdxxb-gx201707003
[15] Clough R, Joseph P. Dynamics of structures[M]. 3rd ed. New York:McGraw-Hill, 2003:234-236.
-
期刊类型引用(4)
1. 刘西安,杨德庆,刘见华. 加速度载荷信息部分缺失下船舶机械噪声快速预报. 振动与冲击. 2022(12): 247-253 . 百度学术
2. 易跃峰,杨德庆,刘西安,陈林. 船舶振动声学优化设计中设备运动激励的等效反演力方法. 计算力学学报. 2022(04): 397-403 . 百度学术
3. 王霂,王洪波. 输入测点误差对结构振声预报影响的实验研究. 科技创新与应用. 2021(02): 8-13 . 百度学术
4. 唐瑞鹏,孙国玺,张锋. 大数据背景下船舶移动网络信息安全风险预测模型. 舰船科学技术. 2021(20): 163-165 . 百度学术
其他类型引用(6)