船舶力学

部分浸没充液无限长圆柱壳声辐射特性研究

PDF下载

王文涛 , 王斌 , 范军 , 周富霖

船舶力学 | 流体与结构声学 2024,28(1): 144-153

收起

船舶力学 | 流体与结构声学 2024, 28(1): 144-153

部分浸没充液无限长圆柱壳声辐射特性研究

全屏

王文涛, 王斌, 范军, 周富霖

作者信息

上海交通大学　海洋工程国家重点实验室　高新船舶与深海开发装备协同创新中心，上海　200240

王文涛（1994-），男，硕士研究生

王斌（1981-），男，博士，教授，博士生导师，通讯作者，E-mail：bin_wang@sjtu.edu.cn。

通讯作者:

通讯作者，E-mail：bin_wang@sjtu.edu.cn

Acoustic radiation characteristics of partially immersed and filled cylindrical shell

Wen-tao WANG, Bin WANG, Jun FAN, Fu-lin ZHOU

Affiliations

Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration, State Key Laboratory of Ocean Engineering, Shanghai Jiao Tong University, Shanghai 200240, China

出版时间: 2024-01-20 doi: 10.3969/j.issn.1007-7294.2024.01.013

文章导航

摘要

收起

本文旨在研究圆柱壳的浸没和充液对其声辐射特征的影响，为界面附近充液壳体目标声学性能的评估和测量提供理论依据和研究方法。通过有限元数值仿真方法计算得到半充液无限长圆柱壳辐射声压的频率-浸没深度谱和半浸没无限长圆柱壳辐射声压的频率-充液高度谱，结果均呈现出两种明显的共振现象，分别由弯曲波和低阶流体附加波激发。根据相速度频散曲线推导出共振频率的预报公式，能够准确地预报辐射声场中出现的共振现象，可为系泊状态内部液舱结构振动声辐射线谱特征控制与评估提供理论支撑。

关键词

圆柱壳 / 声辐射 / 充液 / 弹性波 / 流体附加波

Abstract

收起

This paper aims to study the influence of immersion and liquid filling on the acoustic radiation characteristics of cylindrical shells, and to provide theoretical basis and research methods for the evaluation and measurement of acoustic performance of liquid-filled shells near the interface. The frequency immersion depth spectrum of the radiated sound pressure of a semi-filled infinite cylindrical shell and the frequency liquid-filled height spectrum of a semi-immersed infinite cylindrical shell were calculated by using the finite element numerical simulation method. The results show that there are two obvious resonance phenomena, which are respectively excited by the bending waves and the low-order fluid additional waves. According to the dispersion curve of phase velocity, the prediction formula of resonance frequency was deduced, which can accurately predict the resonance phenomenon in the radiated sound field and provide theoretical support for the control and evaluation of the characteristics of the vibration acoustic radiation line spectrum of the internal tank structure in the mooring state.

Key words

cylindrical shell / sound radiation / liquid-filled shell / elastic wave / fluid additional wave

引用本文

王文涛, 王斌, 范军, 周富霖. 部分浸没充液无限长圆柱壳声辐射特性研究. 船舶力学, 2024 , 28 (1) : 144 -153 . DOI: 10.3969/j.issn.1007-7294.2024.01.013

Wen-tao WANG, Bin WANG, Jun FAN, Fu-lin ZHOU. Acoustic radiation characteristics of partially immersed and filled cylindrical shell[J]. Journal of Ship Mechanics, 2024 , 28 (1) : 144 -153 . DOI: 10.3969/j.issn.1007-7294.2024.01.013

正文

收起

0　引言

收起

壳体是工程应用中的常见结构，例如系泊的潜艇、水面舰船、油气运输管道等均可简化为部分浸没的充液圆柱壳。近些年来，单层圆柱壳、双层圆柱壳^[1]和内部加肋圆柱壳^[2]等的声振特性被广泛研究，但是大多聚焦在全无限流体域中，忽略了自由液面边界的影响，因此本文开展部分浸没充液圆柱壳的研究，重点关注浸没深度和充液高度对其辐射声场的影响。

在现有研究中，Salaün^[3]建立了半浸没状态下无限长圆柱壳声辐射的解析模型，研究表明自由液面对半浸没圆柱壳辐射声场具有重要影响；王斌等^[4]研究了半浸没无限长圆柱壳的振动和声辐射特性，通过对比激励力位于液面上方及下方时壳体的均方振速和辐射声功率，得出壳体均方振速受激励力位置影响较小的结论；叶文兵^[5]基于Flügge薄壳理论建立了半潜无限长圆柱壳振动和声辐射的解析模型，将壳体周围流体域近似为扇形，得到了部分浸没圆柱壳振动的近似理论模型，这些工作为后续研究部分浸没圆柱壳的振动和声特性奠定了理论基础。当壳体部分浸没时，壳体与流体之间的耦合面会随着浸没深度变化，严格解析解较难获得，大多采用数值仿真计算方法，因此Ergin^[6]和Adams ^[7]等利用有限元方法讨论了部分浸没在流体中的圆柱壳自由振动特性。关于部分浸没真空圆柱壳的声辐射研究已较为成熟，而针对内部充液圆柱壳的研究，Li^[8]通过波数域方法研究了半浸没全充液圆柱壳在流动的流体中的声辐射特性，并指出辐射声场与激励点位置、流体流速和壳体厚度有关；林子钦等^[9]利用附带水质量法讨论了部分浸没充液圆柱壳的自由振动特性。但是上述工作对于壳体中弹性波的研究大多聚焦在定量分析而未对其机理进行深入讨论。

本文首先利用有限元软件COMSOL建立半充液无限长圆柱壳和半浸没无限长圆柱壳的仿真模型。然后，为了研究壳体浸没深度变化对辐射特性的影响，计算了半充液无限长圆柱壳辐射声压的频率-浸没深度谱，为了研究壳体充液高度对辐射特性的影响，计算了半浸没无限长圆柱壳辐射声压的频率-充液高度谱。考虑到在后续的研究中方便与现有实际模型对比，计算参数取壳体中面半径a =0.1365 m，厚度h_s =0.004 m，材料为钢，密度ρ_s =7850 kg/m³，杨氏模量E_s =1.94201×10¹¹ Pa，泊松比σ_s =0.25；水介质密度ρ_f =1000 kg/m³，声速c_f =1500 m/s，结果中呈现出两种不同的辐射共振亮线。最后，通过对壳体中弹性环绕波的激发、传播路径以及辐射的研究，分析了共振现象的产生机理，并给出共振频率的预报公式。

1　仿真建模

收起

假设圆柱壳为弹性薄壳，沿轴向无限长，顶部施加轴向线激励F。以壳体截面圆心为原点建立柱坐标系（r，φ，z），其中φ的取值范围为0到2π。

对于部分浸没半充液无限长圆柱壳，其浸没深度动态变化，引入参量h表示壳体底部与外部流体液面间的距离，仿真模型示意图如图1所示。

为了方便对计算结果的观测，引入归一化浸没深度H，令H =（h-a）/a。当壳体下表面刚刚接触外部流体自由液面时H =-1；当壳体恰好完全浸没于壳体外部流体时H =1。因此当H的变化范围为-1~1时表示壳体连续下潜的过程。远场辐射声压观察点取在壳体正横位置的正下方100 m处（r =-100，φ =3π 2，z =0）。力源激励频率范围为20~5000 Hz。

在COMSOL软件中添加二维轴对称几何，并设置上文所述参数和相应的区域以及边界条件，建模结果如图2所示。

然后扫描归一化浸没深度和激励频率，可以得到不同浸没深度下圆柱壳远场辐射声压的频谱结果，最终可得到壳体远场辐射声压的频率-浸没深度谱。

对于半浸没部分充液无限长圆柱壳，建立相同的壳体结构坐标系，如图3所示。

壳体内部充液高度为d，为了方便对计算结果的观测，引入归一化充液高度D，令D =（d-a）/a。当壳体内部真空时D =-1，当壳体内部完全充液时D =1。因此当D的变化范围为-1~1时表示壳体连续充液的过程。

在COMSOL软件中采用类似的建模方法，得到的结果如图4所示。

通过扫描归一化充液高度和激励频率，可以得到不同充液高度下圆柱壳远场辐射声压的频谱结果，最终得到壳体远场辐射声压的频率-充液高度谱。

上述两种工况的模型计算均属于振动问题，在COMSOL计算过程中只扫描参数的变化，原理相同。

2　半充液圆柱壳浸没深度对声辐射特性的影响研究

收起

本章主要研究半充液圆柱壳辐射声场随浸没深度的变化规律，通过对辐射声压的频率-浸没深度谱的分析，基于弹性环绕波理论对结果中呈现出的共振现象做出机理解释。

基于第1章仿真模型和计算方法可以得到半充液圆柱壳辐射声压的频率-浸没深度谱，结果如图5所示。其中：颜色代表声压幅值，单位为Pa；纵轴为圆柱壳归一化浸没深度H；横轴为力源激励频率f。

图5中呈现出两种明显特征：一是存在一系列规律分布且大致平行的倾斜亮线，这是由壳体内弯曲波a₀辐射共振产生的^[10]；二是存在两段独立分布的共振亮线，如图中红色虚线框标注所示，说明壳体内还存在其他形式的波。两种特征均与激励频率以及浸没深度有关。

首先分析a₀波的辐射共振规律，a₀波在壳体中的传播速度分为两种情况：一是壳外流体液面低于充液液面，如图6所示；二是壳外流体液面高于充液液面，如图7所示。

对于情况一，圆柱壳浸没在水中，壳体外部与水面交于A₁、B₁两点，壳体内部与水面交于A₂、B₂两点。以顺时针传播为例，a₀波的路径为A₁—>A₂—>B₂—>B₁—>A₁，会途经壳体两侧均为真空的部分（A₂—>B₂）、壳体外部真空内部充液部分（A₁—>A₂、B₂—>B₁）和壳体两侧均为水的部分（B₁—>A₁）。设a₀波在其中的相速度分别为

、

和

。根据不同部分传播的相速度，可以得到共振频率的预报公式：

(1)

随着壳体浸没深度的增加，θ₁逐渐减小，共振频率f逐渐减小，所以共振峰随着浸没深度的增加向低频方向倾斜，在频率-浸没深度谱上呈现出从左上到右下的共振亮线。

对于情况二，壳体外部与水面交于A₄、B₄两点，壳体内部与水面交于A₃、B₃两点。同样以顺时针传播为例，a₀波的路径为A₄—>B₄—>B₃—>A₃—>A₄，会途经壳体两侧均为真空的部分、壳体外部浸水内部真空部分（B₄—>B₃、A₃—>A₄）和壳体两侧均为水的部分。设a₀波在B₄—>B₃、A₃—>A₄段传播的相速度为

，有

，由此可以得到共振频率的预报公式：

(2)

随着壳体浸没深度的增加，θ₂逐渐增大，共振频率f逐渐减小，共振峰分布向低频倾斜。

由于浸没状态下的充液圆柱壳频散曲线较难获得，且a₀波在圆柱壳中的传播速度与在平板中的传播速度相似，可通过相同材料、相同厚度的板中波的相速度频散曲线预报辐射共振。通过谱方法计算两侧真空的钢板、一侧真空一侧水的钢板和两侧均为水的钢板等三种工况的频散曲线，结果如图8所示。

图中横轴表示激励频率，纵轴表示不同工况钢板中a₀波波速与水中声速的比值。黑色实线表示两侧均为真空的钢板中a₀波相速度频散曲线（

值），红色虚线表示一侧真空一侧为水的钢板中a₀波的相速度频散曲线（

、

值），蓝色点划线表示两侧均为水的钢板中a₀波的相速度频散曲线（

值）。基于预报公式和频散曲线，将计算出的预报结果添加至辐射声压的频率-浸没深度谱中，结果如图9所示。

图中白色虚线即为预报结果，可以看出预报结果与数值仿真计算结果中的亮线吻合得很好，说明导出的预报公式能够准确地解释共振现象。此外，声压的频率-浸没深度谱中除了a₀波贡献的倾斜条纹外，还存在独立分布的共振亮条纹，这是由于壳体内被激发出的一系列流体附加波导致的。流体附加波只在壳体内表面流体中往返传播^[11]，路径如图10中白色曲线所示。

可见流体附加波往返传播会激发较强的辐射共振，为了导出共振频率的预报公式，需要计算半充液圆柱壳中各种波的相速度频散曲线。已知壳体浸没对流体附加波的传播影响很小^[12]，因此本章仅针对外部真空状态下的半充液圆柱壳模型展开计算，为了避免手动插值造成结果不准确，采用COMSOL仿真计算求解，结果如图11所示。

图中横轴表示激励频率，纵轴表示真空中半充液圆柱壳内流体附加波F_i的波速与水中声速的比值。在研究频段内，起主要贡献作用的是低阶流体附加波F₀和F₁。由于壳体内部恒为半充液，因此共振频率的预报公式为

(3)

式中，

表示第i阶流体附加波相速度。分别代入F₀和F₁波的波速值，并将预报结果添加至辐射声压的频率-浸没深度谱中，结果如图12所示。

图12中红色虚线表示F₀波激发的共振峰，红色“*”线表示F₁波激发的共振峰。趋势吻合说明推导的预报公式可以预报辐射声压的频率-浸没深度谱中的共振现象。但是在数值计算结果中，当H在-1~0范围内变化时，流体附加波的辐射很弱，起主要贡献作用的是F₀波，相速度几乎不变，表现为不随浸没深度变化的共振条纹；当H在0~1范围内变化时，起主要贡献作用的是F₁波，表现为随浸没深度变化的弯曲亮条纹，其中的原因有待进一步研究。

3　半浸没圆柱壳内部充液对声辐射特性的影响研究

收起

接下来研究半浸没圆柱壳辐射声场随壳体内部充液高度的变化规律，通过对辐射声压的频率-充液高度谱的分析对呈现出的共振现象做出机理解释。

采用相同计算参数，得到半浸没圆柱壳辐射声压的频率-充液高度谱，结果如图13所示。其中颜色代表声压幅值，单位为Pa；纵轴为圆柱壳归一化充液高度D；横轴为力源激励频率f。

从图13中可以观察到一系列规律分布的由a₀波激发的共振峰和一条独立的由F_i波激发的共振峰。首先分析a₀波辐射共振机理，a₀波的传播也分为两种情况：一是充液高度低于浸没深度，如图14所示；二是充液高度高于浸没深度，如图15所示。

对于情况一，圆柱壳半浸没在水中，壳体外部与水面交于C₂、D₂两点，壳体内部与水面交于C₁、D₁两点。以顺时针传播为例，a₀波的路径为C₁—>C₂—>D₂—>D₁—>C₁，会途经壳体两侧均为真空的部分、壳体外部浸水内部真空部分和壳体两侧均为水的部分。根据不同部分传播的相速度，可以得到共振频率的预报公式：

(4)

随着壳体浸没深度的增加，γ₁逐渐减小，共振频率f逐渐减小，所以共振峰随着浸没深度的增加向低频方向倾斜，在频率-浸没深度谱上呈现出从左上到右下的共振亮线。

对于情况二，壳体外部与水面交于C₃、D₃两点，壳体内部与水面交于C₄、D₄两点。同样以顺时针传播为例，a₀波的路径为C₄—>D₄—>D₃—>C₃—>C₄，会途经壳体两侧均为真空的部分、壳体外部真空内部充液部分和壳体两侧均为水的部分，由此可以得到共振频率的预报公式：

(5)

随着壳体浸没深度的增加，γ₂逐渐增大，共振频率f逐渐减小，共振峰分布向低频倾斜。

由于预报公式中需要的a₀波在不同部分的相速度

、

和

均在图5中给出，因此可直接对共振峰进行预报，将预报结果添加至辐射声压的频率-充液高度谱中，结果如图16所示。

图中白色曲线即为预报结果，可以看出预报结果与数值仿真计算结果中的亮线吻合得很好，说明了预报公式的准确性。接下来分析流体附加波的共振机理，其传播路径与上一章分析的类似，限于篇幅，在此不重复描述。计算外部真空全充液圆柱壳中各种波的相速度频散曲线，结果如图17所示。

图17中横轴表示激励频率，纵轴表示真空中的全充液圆柱壳内流体附加波F_i波速与水中声速的比值。可见在研究频段内，起主要贡献作用的只有F₀波，可得共振频率的预报公式为

(6)

式中，β表示图11中OC₁与OD₁之间的夹角。将预报结果添加至辐射声压的频率-充液高度谱中，结果如图18所示。

图中红色曲线即为F₀波共振的预报结果。由频散曲线易知，F₀波相速度在研究频段内几乎不变，表现在辐射声压的频率-充液高度谱中为近似线性分布，吻合得很好。说明了本章推导的预报公式能准确地预报辐射声压的频率-充液高度谱中的共振现象，同时解释了流体附加波对壳体辐射声场的影响规律。

4　结论

收起

本文分析了部分浸没圆柱壳和部分充液圆柱壳辐射声场的构成，利用有限元软件建立了部分浸没半充液无限长圆柱壳和半浸没部分充液无限长圆柱壳辐射的仿真模型，基于COMSOL软件计算得到辐射声压的频率-浸没深度/充液高度谱。通过壳体中各种波的相速度频散曲线解释了辐射声场中呈现的辐射共振现象。得到如下结论：

（1）部分浸没半充液无限长圆柱壳和半浸没部分充液无限长圆柱壳辐射声场均呈现出一系列规律分布的共振峰，这是由于壳体中弯曲波a₀在不同部分的传播速度不同导致的；

（2）部分浸没半充液圆柱壳的辐射声场中存在的独立分布的共振峰，是由壳体内低阶流体附加环绕波F₀和F₁激发的，其分布规律与壳体浸没深度有关；

（3）半浸没部分充液圆柱壳的辐射声场中存在的独立共振峰则主要由F₀波贡献，在壳体辐射的频率-充液高度谱中近似呈现线性分布规律；

（4）基于壳体中各种波的相速度频散曲线推导出的频率预报公式能够准确地预报声场中出现的共振现象；

（5）经过验证，当壳体的厚径比改变时，共振频点间隔随之改变，但是分布规律以及产生的机理相同，因此本文结论具有普适性。

基金

收起

国家自然科学基金资助项目(52001211)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

杨晓刚，林立，等. 径向激励作用下有限长双层圆柱壳声振传递特性研究[J]. 船舶力学，2015，19(5): 600-608.

Yang

, Lin

, el al. Sound and vibration transmission characteristic of finite length double cylindrical shells under radial excitation[J]. Journal of Ship Mechanics, 2015, 19(5): 600-608. (in Chinese)

[2]

谢官模，李军向，等. 环肋、舱壁和纵骨加强的无限长圆柱壳在水下的声辐射特性[J]. 船舶力学，2004，8(1): 101-108.

Xie

, Li

, el al. Sound radiation from fluid-loaded infinite cylindrical shells with rings, bulkheads and longitudinal stiffeners[J]. Journal of Ship Mechanics, 2004, 8(1): 101-108. (in Chinese)

[3]

Salaün

. Effect of a free surface on the far-field sound pressure radiated by a point-excited cylindrical shell[J]. Journal of the Acoustical Society of America, 1991, 90: 2173-2181.

[4]

王斌，汤渭霖. 半潜状态圆柱壳振动声辐射特性研究[C]//第十二届船舶水下噪声学术讨论会论文集，2009.

Wang

, Tang

. Study on vibration and sound radiation characteristics of semi-submerged cylindrical shell[C]//Proceedings of the 12th Symposium on Underwater Noise from Ships, 2009. (in Chinese)

[5]

叶文兵. 潜浮状态下圆柱壳结构的声振特性研究[D]. 武汉：华中科技大学，2012.

. Study on acoustic and vibration characteristics of cylindrical shell structure in submerged floating state[D]. Wuhan: Huazhong University of Science and Technology, 2012. (in Chinese)

[6]

Ergin

, Temarel

. Free vibration of a partially liquid-filled and immersed, horizontal cylindrical shell[J]. Journal of Sound and Vibration, 2002, 254: 951-965.

[7]

Adams

S D

, Craster

R V

, Williams

D P

, et al. Two-dimensional frequency domain scattering and point excitation of partially filled elastic cylinders[J]. Journal of Sound and Vibration, 2010, 329(11): 2147-2157.

[8]

H L

, Wu

, Huang

, et al. Parametric study on sound radiation from an infinite fluid-filled/semi-immersed cylindrical shell[J]. Applied Acoustics, 2003, 64(5): 495-509.

[9]

林子钦，许学强，等. 基于虚拟质量法的部分充液且部分浸没圆柱壳自振特性分析[J]. 中国水运（下半月），2018，18(12): 156-160.

Lin

, Xu

, el al. Analysis of natural vibration characteristics of partially submerged cylindrical shell filled with liquid based on virtual mass method[J]. China Water Transport (Second Half of the Month), 2018, 18(12): 156-160. (in Chinese)

[10]

赵开琦. 界面附近目标低频振动与声辐射特性研究[D]. 上海：上海交通大学，2020.

Zhao

. Research on the low frequency vibroacoustic behavior of near-surface target[D]. Shanghai: Shanghai Jiao Tong University, 2020. (in Chinese)

[11]

Tong

, Fan

, Wang

, et al. Vibroacoustic behavior of an infinitely long cylindrical shell with periodic internal lengthwise ribs[J]. Journal of the Acoustical Society of America, 2018, 144(3): 1528-1538.

[12]

范威，郑国垠，范军. 充水圆柱壳声散射的环绕波分析[J]. 声学学报，2010，35(7): 419-426.

Fan

, Zheng

, Fan

. Surround wave analysis of sound scattering from water filled cylindrical shell[J]. Acta Acustica, 2010, 35(7): 419-426. (in Chinese)

2024年第28卷第1期

PDF下载

引用本文

BibTeX

文章信息

doi: 10.3969/j.issn.1007-7294.2024.01.013

接收时间：2023-07-10
首发时间：2026-03-21
出版时间：2024-01-20

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2023-07-10

基金

国家自然科学基金资助项目(52001211)

作者信息

上海交通大学　海洋工程国家重点实验室　高新船舶与深海开发装备协同创新中心，上海　200240

通讯作者:

通讯作者，E-mail：bin_wang@sjtu.edu.cn

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/cblx/CN/10.3969/j.issn.1007-7294.2024.01.013

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT