大型客滚船首部砰击引发舱室噪声的数值模拟.pdf

第45卷第2 0 期2023年10 月舰船科学技术SHIP SCIENCEANDTECHNOLOGYVol.45,No.20Oct.,2023大型客滚船首部砰击引发舱室噪声的数值模拟何丽丝，王德禹(上海交通大学船舶海洋与建筑工程学院，上海2 0 0 2 40)摘要：噪声规范中以校核机械振动产生的噪声为主，但高幅值的砰击力周期性地与船体相互作用，已成为船舶首部舱室的重要噪声源，故开展由击引发的舱室噪声问题数值模拟研究。以某大型客滚船为例，采用CFD软件STAR-CCM+对船体首部的砰击现象进行数值模拟，编制时域-频域转化程序，并借助声学分析软件VAOne，将砰击力作为外载荷计人声学分析模型，获得船体典型舱室内的噪声水平。结果表明，该计算方法能够快速有效地得到由砰击引发的船舶舱内噪声水平。针对某大型客滚船首部区域，每增加一级海况，船体首部由于砰击引发的声压级提高约3dB左右。该计算流程实现了流体力学与声学计算一体化，给出一种关于砰击现象引发舱室内噪声的数值预报方法，为实船降噪提供了可靠的保障，对相关船型的设计工作具有参考价值。关键词：砰击；声学分析；计算流体力学；时域-频域转化；噪声评价中图分类号：U662.2文章编号：16 7 2-7 6 49(2 0 2 3)2 0-0 0 2 5-0 7Numerical simulation of the bow slamming induced cabin noise of(School of Naval Architecture,Ocean&Civil Engineering,Shanghai Jiaotong University,Shanghai 200240,China)Abstract:The existing regulation of cabin noise is the noise of this focus due to the mechanical vibration.Since theperiodical interaction between the high amplitude slamming forces at the ship bow may be an important noise source of thecabins,numerical investigations on the bow slamming induced cabin noise were conducted.Taking a large RO-RO ship as acomputational example,numerical studies were performed by combing the computational fluid dynamics(CFD)and acous-tic simulation.The slamming phenomena were simulated using the CFD code STAR-CCM+.The computed time-domainslamming forces were then converted into the frequency-domain using the self-developed signal converter before being sub-stituted into the acoustic analysis as the external excitations.The acoustic simulation yielded finally the noise level of theoverall ship hull and some typical cabins.Numerical results show that the sound pressure level caused by the slamming in-creases about 3 dB with each increment of the sea state level.A numerical prediction method for cabin noise caused by slam-ming is presented,which provides a useful tool for noise performance analysis of large passenger ships.The presented nu-merical results are of reference values to the design of related ship types.Key words:slamming;acoustic analysis;computational fluid dynamics;time domain-frequency domain trans-form;noise rating0引言船舶舱室内的噪声问题是船舶工程中的重要研究领域之一。传统上，船舶噪声问题研究多着眼于由主机、辅机、泵、螺旋桨、空调通风设备等各种机械设备激发船体结构振动所引发的噪声!。现阶段的噪声收稿日期：2 0 2 2-0 9-2 0基金项目：教育部、财政部重大科研项目（教技函2 0 13 35号）作者简介：何丽丝（198 8-），女，硕士，工程师，研究方向为船舶结构的可靠性和安全性。文献标识码：Aa large RO-RO passenger shipHE Li-si,WANG De-yudoi:10.3404/j.issn.1672-7649.2023.20.005相关规范也仅限于校核船舶机械振动引起的噪声。而近年来，随着国际海事组织（IMO）及各国船级社不断出台更为严格的船舶噪声控制标准，对于船舶噪声问题的研究开始朝着更加系统化、复杂化的方向发展2-3。事实上，海上航行船舶在遭遇较大风浪时容易26产生大幅度的垂荡和纵摇运动，这些运动使得船首部分与海水在很短时间内发生大幅度的相对运动，形成所谓的“砰击”现象。砰击力作用在船体表面，可能成为船体首部舱室的重要噪声源。尽管砰击现象引发舱室噪声的物理现象比较直观，但对于船舶首部砰击引发舱室噪声的相关研究却较少见诸文献。究其原因，很大程度上是由于砰击问题本身是一个瞬时发生的强非线性问题，要准确获得船体首部不同位置处的砰击压力是困难的。早期对于物体入水砰击引起的噪声问题研究多采用实验方法4-9。近年来，借助于计算流体力学（CFD）技术的快速发展，对船体砰击问题的数值模拟已经取得了很大的进展。借助高精度的CFD计算，开展砰击引起噪声问题研究逐渐成为可能。邓方舟7 开展了低速物体人水声的特性研究，Barjasteh等8 则对楔形体入水砰击造成的水下噪声进行了数值和实验研究。以1艘大型客滚船为对象，开展大型客滚船首部砰击引发舱室噪声的数值模拟。由于客滚船运载旅客，对舱室内噪声控制标准就比较严苛，但在相关规范中引起噪声的原因未涉及首部砰击。另一方面，客滚船车辆甲板为前后贯通的有大型舱室，因此对于噪声的抑制难度更高。为了评价船舶在航行过程中由于砰击引发的舱室噪声，首先采用计算流体力学软件STAR-CCM+对船舶的砰击现象进行数值模拟，并通过在船体首部不同位置建立监测点，获得典型位置上时域状态下的砰击压力和船体的入水速度。编制时域-频域转化软件，借助声学分析软件VAOne将压力和入水速度作为边界条件（频域）输人VAOne软件中，反映不同强度的砰击对舱室内噪声的影响。建立数学模型采用的数值方法，对于计算船体由于入水砰击产生的舱室噪声问题具有有效性，得到的数值结果也能为相关船型设计提供参考数据。1计算方法及船型选择1.1砰击引发舱室噪声的计算方法采用CFD软件STAR-CCM+对船舶在不同海况下航行时的砰击现象进行模拟，取随船运动坐标系。坐标系的远点取在船舶静止时水线面的形心，x轴指向船首，轴指向左舫，z轴铅直向上。将船体沿x、J、z轴的线运动分别称为纵荡、横荡和垂荡，绕x、y、z轴的角运动则称为横摇、纵摇和首摇。在STAR-CCM+软件平台中建立了虚拟波浪水池。为简化问题，只考虑规则入射波的情况，在STAR-CCM+平台中已经内置了常用的人射波形式，因此只需设置数值波舰船科学技术浪水池的区域和边界条件，即可进行人射波的数值模拟。参考之前的波浪-船舶相互作用计算案例，使用如图1所示的计算区域及边界条件，具体为：流域左侧为速度入口边界，右侧为压力出口边界，两侧和底部均为无滑移壁面边界，顶部为速度入口边界。其中，上游和下游均设置3倍波长（3L），两侧设置2.5倍波长(2.5L)。入射波图1数值波浪水池立体示意图Fig,1 Schematic diagram of the numerical wave tank以VOF（V o l u m e o f Fl u i d）流体体积法捕捉自由液面的变化。时间离散采用隐式非定常时间推进格式。所有模拟均采用RealizableK-Epsilon2层湍流模型，应用两层全y+壁面处理模型。使用SIMPLE（Se m i-I m p l i-cit Method for Pressure-Linked Equations）算法求解压力场和速度场。船舶在波浪上的运动情况采用STAR-CCM+中的DFBI（D y n a mi c Fl u i d-Bo d y In t e r a c t i o n）模型，对联立后的方程组进行求解。CFD计算得到的是砰击力的时域信号，而目前常用的噪声分析数值工具多采用频域方法，两者之间存在数据转换的障碍。砰击噪声主要考虑的是船舶在高海况下，由于波浪诱导的大幅度运动而导致的首部砰击问题，尽管该物理过程存在强烈的非线性效应，但船体收到的砰击力本质上仍属于周期性激励。基于此，采用离散傅里叶变换算法，自主开发船舶水砰击力信号的时域-频域转化软件。以船舶的砰击激励力时域信号为输人，通过内置数值计算和转换工具，实现时域信号到频域信号的快速计算。舱室噪声计算采用声学分析软件VAOne中的统计能量法（SEA）进行，Burroughs等9 对该方法进行详细的阐述。在VAOne软件平台中应用SEA方法进行船舶舱室内的噪声分析时，首先需要建立目标船的声学分析统计能量模型（简称“声学模型”），而后需要将噪声的激励源以外激励的方式施加在声学模型上，即可通过求解统计能量方程获得各声腔中的声压级分布情况。第45卷顶部（速度入口）波面底部（无滑移壁面）第45卷1.2目标船型的选取选择总长为2 0 0 m的某大型客滚船作为目标船型开展相关计算。目标船型自下而上共有9层甲板，第3、第5甲板为贯穿全船的大型车道舱室，第2 甲板高度9.3m，第5甲板高度15m。第7、8、9 甲板为载客甲板，其距离基线高度分别为2 3.5m、2 6.5m和2 9.2 5m。其中，第1甲板为双层底，第4和第6 甲板为船中部甲板室内的平台。对于所考虑的舱室噪声问题而言，这几层甲板的影响较低，因此在后续声学模型的建立中将其忽略。目标船的总布置情况如图2 所示。何丽丝，等：大型客滚船首部砰击引发舱室噪声的数值模拟274326422540244123(a)160m(FR202)和17 3m(FR219)处的砰击压力监测点27.441910图2 目标船型的布置情况13Fig.2General arrangement of the subject ship14(b)183.6m(FR236)和191.4m(FR249)处的砰击压力监测点2模型的建立图4砰击压力监测点的布置情况2.1CFD计算模型的建立对于涉及的CFD计算部分而言，只需船体的几何外形即可。为此，根据船体的肋骨型线图进行了船体表面的几何建模，并对船体的上层建筑进行简化后，得到的CFD计算模型如图3所示。图3简化后用于CFD计算的三维模型Fig.3Simplified ship model for the CFD computation为了获得船首在砰击现象发生时各部分所受的压力情况，在船体首部总共设置45个压力监测点，具体分布如图4所示。这些监测点位于船体表面曲率变化的显著位置，能够描述不同情况下，船体外部水压力对船体结构产生的激励情况，2.2声学模型的建立在VAOne软件平台中，根据船体的肋骨型线图和基本结构图，建立目标船的SEA声学模型，主要使用SEA平板（SEAPlate)对船体中部及上层建筑进行建模，船舶首尾曲率较大处则采用SEA单曲率曲面板壳（SEA s in g le c u r v e d s h e ll）进行表达。通过依次连接各横剖面上的节点，最终形成船体表面的板壳模型如图5所示。在完成全船SEA建模之后，即可对船体进行SEA声腔（Cavity）的划分。由于本研究主要考虑船舶首部发生砰击所诱发的舱室噪声情况，根据船体的结Fig.4Arrangement of pressure monitoring point图5VAOne中建立的全船SEA模型分解视图Fig.5 Shrink view of the SEA model built in VA One构形式及布置情况，在船首区域内声腔的划分较细，船舶中部、尾部由于距离碎击发生处相对较远，噪声的水平不高，可适当增加声腔的尺寸以简化建模和计算的过程。最终全船划分成为7 0 个声腔，具体的分布情况如图6 所示。图6 VAOne中建立的全船声腔模型分解视图Fig.6Shrink view of the SEA cavities built in VA One根据CFD中压力监测点的布置情况，首先将船体首部左右两舫的SEA板壳各划分为15个单元如图7 所示。假定每个单元上各点的压力变化情况相同，则将各单元对应监测点监测到的砰击压力乘以单元面积后，便可得到不同单元上的砰击压力情况。考虑到船首15单元上的压力变化剧烈，将每个单元上对应的压力监测信息进行平均，作为该单元上的砰击力载荷。28（图中方形标志为压力监测点的位置，锥形箭头表示SEA模型力载荷的位置）图7 船体首部砰击压力的计算单元分布Fig.7Computation panels for slamming pressure at the bow确定单元分布后，将点力载荷（PointForce）施加于每个单元的形心位置，即可完成力载荷的施加。SEA模型通过力的频谱定义载荷的强度，因此在力的载荷布置完成之后，每一个载荷都关联一个频谱。后续计算中，只需根据不同海况下的砰击压力更新此频谱，便可模拟同一位置不同载荷引起的噪声问题，对载荷的施加位置不再进行调整，3计算流程和目标区域3.1时域-频域转化软件采用自编程的时域-频域信号转换软件，实现信号由时域到频域的转换，如图8 所示。入读时域信号后，软件自动读取信号的长度和采用间隔。考虑到CFD计算在初始时刻存在一定的数值误差，软件中允许用户自定义进行傅里叶编号的初始信号序列。以图8 信号为例，该信号的采样间隔为0.0 2 s（即CFD计算的时间步长），指定从第6 0 0 个采样点开始，选取后续的10 2 4时域信号4.445F2.96314820-1.482-2.963-4445频域信号2.490F199214949964980Fig.8 The functional interface of the software舰船科学技术个信号进行离散傅里叶变换，可以快速得出具体的频谱及其对应的数值结果。3.2噪声监测位置的确定11聚焦于船体首部发生砰击而引发的舱室内噪声问12题，尽管船首各声腔距离砰击点的位置最近，但这些声腔实际上多为封闭舱室和压载舱室，对噪声问题的13研究意义不大。根据目标船的实际布置情况，选择如148510期域计算结果：2.043.292.535.963.128.244.020.595.014.206.38.768.04.6610.02.2212.50.8816.00.2420.00.0525.00.0131.50.0040.00.0050.00.0063.00.0025图8 软件功能界面第45卷图9 所示的5个舱室，作为噪声监测的重点位置。将其分别命名为第1 第5目标声腔。第1、第2 目标声腔分别位于第3甲板和第5甲板前部，是客滚船车辆停放的区域，其舱室空间大，距离砰击点相对较近，可能发生较大的砰击噪声响应。第3 第5目标声腔位于上层建筑甲板室的前部，是旅客活动区域，对噪声水平有较高的要求。图9舱室噪声监测位置选择情况Fig.9 Selected cabin for noise monitoring4计算结果分析4.1砰击问题的数值模拟结果考虑目标船在58 级海况下迎浪航行的情况，采用工程中常用的“设计波”简化方法，即以波能相当的规则波代替不规则波，反映船舶在不同海况中的运动响应情况。计算中不同海况下实际使用的设计波波幅及周期如表1所示10 1，计算航速根据船舶的设计航速，取2 8.5kn。计算得到不同海况下船舶在波浪上航行时的气-液2.000两相流体的流动情况，其中每种海况间隔的时长为半个遭遇周期。结果表明，5级海况下船体首部相对于流体的运动并不明显，但海况增加以后，船首人水的情况越发明显。在8 级海况下，船首已经出现大幅度的出水和人水情况，其伴随的砰击和上浪现象非常显著。图10 给出了8 级海况下，几个典型压力监测点监测到的压力随时间变换情况。可以看出，各压力监测71080.00.00100.00.00125.00.00160.0-0.00点得出的压力结果并非完全规则的正弦余弦曲线，而是呈现出类似脉冲信号的变化规律，这些类似于压力脉冲的信号反映出砰击现象发生时的压力变化。在船第45卷表1不同海况下设计波的波幅和周期Tab.1Amplitude and period of designed waves indifferent sea states海况(S.S.)5波高（即2 倍波幅）/m3.25波周期/s9.7240000200000160000120000出80000400000-400000(a)压力监测点6，9，12，13，14结果240000监测点2 3200000监测点2 4监测点2 5160000监测点2 6120000监测点2 7力80000出400000-400000200000160000120000ed/80000400000-400000图10 八级海况(S.S.8)下典型位置压力监测点监测到的Fig.10Time history of the monitored pressure under Sea State 8体进行大幅度垂荡-纵摇耦合运动，部分监测点会不断离开并回到水中，因此该点的压力不断在空气压力和水压力之间跳跃。在监测点接触水面的短暂瞬间，压力会以砰击现象的发生瞬间达到一个较高水平，而后随着船体入水，压力水平逐渐降低。当监测点离开水体接触空气时，压力则接近于0。这样的过程不断反复，便形成了对舱室内噪声的激发效果。4.2砰击引发舱室内噪声的数值结果为计算不同海况下砰击现象引发的舱室噪声情况，假定图7 中定义的每个单元上各点的压力变化情况相同。在每种海况下，首先将各单元上包含的压力监测点所测得的压力时历曲线进行平均，再与单元的面积相乘后，以获得得到该单元时域的激励力信号。何丽丝，等：大型客滚船首部砰击引发舱室噪声的数值模拟6757.512.415监测点6监测点9监测点12监测点13-监测点142040时间/sH20(b)压力监测点2 3 2 7 结果2040时间/s(c)压力监测点40 44结果压力变化时历结果29通过时域-频域转化软件，即可得到该单元上的激励力频谱。图11展示了8 级海况(S.S.8)第2 个单元上的砰击力频谱。可以看出，尽管CFD计算中只考虑了规则8入水波的情况，但由于砰击问题的非线性特性，得到11.5的能量频谱还是覆盖了一定的频率区间，这也符合16.4SEA方法对载荷频谱特性的要求。由于声学计算模型在船体左右两各划分了15个单元，因此每个海况下都需要通过时域-频域转化软件求得15个激励力频谱。将这些频谱作为外激励加载到VAOne软件中，即可得到舱室内的声学水平。10:00010006010010N/半楼10.10.010.0010.00011e-51e-614060时间/s监测点40监测点41监测点42监测点43监测点446010频率/Hz图11八级海况(S.S.8)第2 个单元上的砰击力频谱Fig.11Computed slamming force spectrum for the panelNo.2 under sea state 8图12 给出了不同海况下船体总体声压的计算结果。可以看出，由于砰击力载荷作用于首部，在球鼻首及其附近产生的声压级最高，58 海况下对应的最高声压分别为10 6.2 dB、110 d B、118.4d B和12 1.6 dB。而随着舱室的位置逐渐远离砰击发生的位置，各舱室内的声压水平显著下降。船首部分的声压级虽高，但由于其远离人员活动区域，并不会造成明显的危害。相比较而言，人员活动区域的声压水平更受关注。为此，将图9中5个目标声腔中的噪声水平进行汇总，并采用国际上常用的NR曲线方法，对噪声的水平进行评价。NR曲线噪声评价方法采用“相切法”（tangencymethod）确定噪声的NR值，具体如下：先测量各个倍频带声压级，再把倍频带噪声谱合在NR曲线上，与频谱相切的最高NR曲线为该噪声的NR噪声等级。ANSI/ASAS12.2-2008标准11推荐的不同场合NR曲线如表2 所示。将各目标声腔在不同海况等级下的声压级计算曲线分别绘人NR曲线图，所得到的结果如图13所示。表3则列出了各舱室在不同海况下的最终NR等级。目标舱室1和2 分别位于第3和第5甲板的前端，这是客滚船装载车辆的区域，在57 级海况之下，这2个舱室的噪声水平在NR-20上下，属于静谧的环10030图12 各海况下船体总体声压分布计算结果Fig.12Computed overall sound level of the shipat different Sea State表2 各类建筑物室内允许噪声级NR曲线The allowable NR level for various room noisesTab.2NR曲线建筑物类别广播录音室、播音室音乐厅、剧院、电视演播室电影院、演讲厅、会议厅办公室、阅览室、设计室餐厅、宴会厅、体育馆、商场候机厅、候车厅、候船厅洁净车间、带机械设备的办公室舰船科学技术1001062102.601.7695.3884.5281477.28973.6670.0466.4262.8591855.5651.94vAOne(a)5级海况11010696.1310219417902286.2678.35K2.3174.466.4970.4562.5458.5854.6350.67vVAOne(b)6级海况n118.4114100-41124102598:4894.590.5386.552.578674.6370.656.6862.758.72vAOne(c)7级海况12126118:11146110.6104.1100.697.0493.590.086.579.4782.9972.4575.9668:94vAOner(d)8级海况1520202525303035354540455060第45卷S.S.5S.S.6NR-9080S.S.73.S.8NR-80NR-7060NR-60NR-50租4020031.510080604020031.5100808P/田604020031.5100808P/60租4020031.5100808P/王604020031.5图13目标舱室的NR曲线评定结果Fig.13 Evaluated NR-Level of the target cabinNR-40NR-30NR-20NR-1063125250频率/Hz(a)第1目标声腔S.S.53:6NR-90S.S.7S.S.8NR-80NR-60NR-70NR-50NR-40NR-30NR-20NR-10631252505001000频率/Hz(b)第2 目标声腔S.S.5S.S.6NR-90S.S.7S.S.8NR-80NR-70NR-60NR-50NR-40NR-30NR-20NR-10631252505001000频率/Hz()第3目标声腔S.S.5S.S.6NR-90S.S.7S.S.8NR-80NR-70NR-60NR-50NR-40NR-30NR-20NR-10631252505001000频率/Hz(d)第4目标声腔S.S.53:6NR-90S.S.7S.S.8NR-80NR-70NR-60NR-50NR-40NR-30NR-20NR-10631252505001000频率/Hz(e)第5目标声腔5001000第45卷表3各目标舱室的NR曲线评级结果Tab.3 Summary of the evaluated NR level of selected cabins5级海况8级海况目标舱室1NR-19目标舱室2NR-17目标舱室3NR-10目标舱室4NR-10目标舱室5NR-10境。目标舱室35属于人员活动的区域，在57 级海况之下，这些舱室内由于砰击现象造成的噪声水平低于NR-10，说明此时的砰击现象不会激起上层建筑内的噪声情况。8级海况下，各舱室内的噪声水平均发生显著提升。这是由于此海况下船体首部已出现了大范围的人水现象，此时船体首部所受的激励力与57 级海况下相比已出现明显变化。在8 级海况下，目标舱室1舱室3的NR等级在50 左右，目标舱室4内的噪声等级NR-45，目标舱室5内则为NR-20。对标表2 中给出的噪声参考标准可以看出，此时目标舱室1舱室4中的噪声已值得关注，特别是在目标舱室4内，对于人员活动的区域，较高的噪声等级可能影响旅客的乘坐舒适度，并可能导致一定的安全隐患。船舶在8 级海况或更高级海况航行下，可通过主动降低航速的方法，减少在高海况下大幅度垂荡和纵摇运动的产生。垂荡和纵摇耦合运动的减弱会带来砰击力的显著降低，相应地，由于砰击产生的噪声水平也会明显下降。5结语针对船舶在波浪中航行时，由于受到波浪的作用所发生的入水碎击现象，以及砰击现象诱发的舱室噪声进行了数值研究。以一艘大型客滚船为对象，对其在58级海况下迎浪航向时的首部砰击现象进行了模拟，并对5个典型舱室内由于砰击引发的噪声进行了计算。砰击现象的模拟采用计算流体力学（CFD）软件STAR-CCM+进行，声学分析则借助VAOne软件中的SEA模块进行。通过将CFD计算得到的结果经过FFT转换为频谱后输入声学分析软件中，获得了砰击引发的舱室噪声数值预报结果。通过计算发现，目标船型航行时海况等级每提升一级，砰击引起的声学能量激励增加约一个数量级，相应地，全船总体的噪声水平提升3dB左右。为评价船体由于砰击而产生的噪声水平，在船体前部的5个典型位置处选择了5个目标声腔，其中第1和第2 目标声腔则分别位于第1层和第2 层车辆甲板靠近船首何丽丝，等：大型客滚船首部砰击引发舱室噪声的数值模拟6级海况7级海况NR-19NR-21NR-16NR-18NR-10NR-10NR-10NR-10NR-10NR-1031的区域，第3、4、5、目标声腔则位于甲板室前部靠近船首的位置。引人国际上常用的NR曲线方法对各目标舱室的噪声等级进行评价后发现，57 级海况NR-48下，船体5个目标声腔的噪声水平都比较低，但在NR-508级海况下，各目标舱室内的噪声相当可观。造成这NR-49种现象的原因是因为8 级海况下，船首已出现了显著NR-45的入水情况，各压力监测点的出现强非线性的变化，NR-20导致砰击力激励频谱的幅值显著提升，能量分布的频率区间也更大。要降低砰击引发的舱室噪声水平，改善船舶的航向状态，尽量选择较低的海况以减少砰击载荷是一种最为简单可行的方案。实现了流体力学CFD计算与声学统计能量SEA计算的结合，船舶流体力学与声学计算一体化流程，给出了一种简便可行的预报船舶在波浪上航行时由于砰击引发舱内噪声的计算方法。为大型客船的设计和噪声性能分析提供了有效实用的工具，同时对相关噪声规范中未涉及的噪声计算起到指导作用。参考文献：1陈小剑.舰船噪声控制技术M.上海：上海交通大学出版社,2 0 132 International Maritime Organization,Code on noise levels onboard shipsS.2014.【3】李倩，刘志坚.吉春正.MLC2006公约新造船验证分析船舶,2 0 14,2 0 14(5):114-1147LI Q,LIU Z J,JI C Z.Analysis of verifi cation of newly builtships following MLC2006J.Ship,2014,2014(5):114-1147.4 RICHARDSON E G.The sounds of impact of a solid on a liq-uid surfaceJ.Proceedings of the Physical Society Section B,1955,61(4):352.5鲍筱玲，汤渭霖.击水声特性研究.哈尔滨船舶工程学院学报,198 1,(1):2 7-39.BAO X L,TANG W L.Investigations of the impact noiseJ.Journal of Harbin shipbuilding Engineering Institute,1981,(1):27-39.6汤渭霖，俞孟萨，王斌,水动力噪声理论M.北京：科学出版社,2 0 19.【7】邓方舟.低速物体入水声特性研究D.哈尔滨：哈尔滨工程大学.2 0 2 0.8 BARJASTEH M,ZERAATGAR H.A contribution to under-water acoustic problem in water entry of wedge-shapedbodiesJ.Applied Ocean Research,2021.111:p.102649.9 BURROUGHS C,FISCHER R,KERN F,An introduction tostatistical energy analysisJ,Journaal of Acoustical Society ofAmerican,1997.10 FALTINSEN,O M.Sea loads on ships and offshorestructuresMJ.Cambridge:Cambridge University Press.1990.11 American National Standard.Criteria for evaluating roomnoiseS.ANSI/ASA S12.2,2008.