CN113190912B - 舰艇振动壳体的辐射噪声预测方法 - Google Patents

舰艇振动壳体的辐射噪声预测方法 Download PDF

Info

Publication number
CN113190912B
CN113190912B CN202110273677.2A CN202110273677A CN113190912B CN 113190912 B CN113190912 B CN 113190912B CN 202110273677 A CN202110273677 A CN 202110273677A CN 113190912 B CN113190912 B CN 113190912B
Authority
CN
China
Prior art keywords
shell
sound
tangential
sound field
vibration
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110273677.2A
Other languages
English (en)
Other versions
CN113190912A (zh
Inventor
张兴武
何燕飞
张涛
刘金鑫
陈雪峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN202110273677.2A priority Critical patent/CN113190912B/zh
Publication of CN113190912A publication Critical patent/CN113190912A/zh
Application granted granted Critical
Publication of CN113190912B publication Critical patent/CN113190912B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
    • G01H1/12—Measuring characteristics of vibrations in solids by using direct conduction to the detector of longitudinal or not specified vibrations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Soundproofing, Sound Blocking, And Sound Damping (AREA)

Abstract

公开舰艇振动壳体的辐射噪声预测方法，方法中，采集舰艇振动壳体的振动信号，对两个切向加速度信号
Figure DDA0002972971010000011
和一个法向加速度信号
Figure DDA0002972971010000012
进行时间积分得到两个切向速度信号
Figure DDA0002972971010000013
和一个法向速度信号
Figure DDA0002972971010000014
基于声波速度和声压的控制方程空间离散化和时间积分，构造舰艇振动壳体辐射噪声预测模型，其中，声场变量包括声压大小p(x，t)、声波切向速度uτ′(x，t)和uτ″(x，t)、声波法向速度un(x，t)，x表示为声辐射场的不同位置，针对声场变量采用DRP空间离散格式和LDDRK时间积分格式构造PML的无反射预测模型；在模型中添加切向速度变量
Figure DDA0002972971010000015
Figure DDA0002972971010000016
和法向速度变量
Figure DDA0002972971010000017
以计算壳体振动情况下的声场随时间的变化情况；计算每一时间步长下，舰艇振动壳体运动对声场变量的影响以预测每一时刻舰艇振动壳体的声辐射大小。

Description

舰艇振动壳体的辐射噪声预测方法
技术领域
本发明涉及声学仿真技术领域，尤其涉及一种舰艇振动壳体的辐射噪声预测方法。
背景技术
舰艇隐身性能是舰艇生存能力和战斗力的关键特征，也是各国提升海上军事威慑力的关键技术。然而，在浩瀚的海洋水介质中，声场是唯一能够远距离传递信息的物理场。因此，舰艇水下辐射噪声是其隐身性能的核心指标。为了降低舰艇水下辐射噪声水平，在设计阶段常用缩比试制的方法测试该指标的大小。但是缩比试制方法不仅耗费大量的人力物力，而且受到尺度效应和实验环境的影响。由于数值仿真具备适应性强、不受环境影响、耗费小、周期短等优势，逐渐成为设计阶段的主要分析手段。但是传统噪声预测方法精度不高，并且容易受到边界效应的影响，不能精确有效地分析振动壳体的辐射噪声大小。
在背景技术部分中公开的上述信息仅仅用于增强对本发明背景的理解，因此可能包含不构成本领域普通技术人员公知的现有技术的信息。
发明内容
本发明的目的是提供一种舰艇振动壳体的辐射噪声预测方法。采用已知的振动壳体加速度信号，利用高精度的空间离散格式和时间积分格式，以及消除边界反射效应的边界吸收层，从而构造出振动壳体辐射噪声预测模型。在此基础上添加随时间变化的壳体运动速度边界条件，以此实现辐射噪声的精确预测。
为了实现上述目的，本发明提供如下技术方案：
本发明的一种舰艇振动壳体的辐射噪声预测方法包括：
S100：采集舰艇振动壳体的振动信号，得到壳体的切向加速度信号
Figure BDA0002972968990000011
和法向加速度信号
Figure BDA0002972968990000012
其中X(t)表示为壳体边界的位置；
S200：对所述切向加速度信号
Figure BDA0002972968990000013
和法向加速度信号
Figure BDA0002972968990000014
进行时间积分得到切向速度
Figure BDA0002972968990000015
和法向速度
Figure BDA0002972968990000016
S300：基于声波速度和声压的控制方程空间离散化和时间积分，构造舰艇振动壳体辐射噪声预测模型，其中，声场变量包括：声压大小p(x，t)、声波切向速度uτ′(x，t)和uτ″(x，t)、声波法向速度un(x，t)，其中，p(x，t)和uτ′(x，t)以及uτ″(x，t)中的x表示声辐射场的不同位置，针对所述声场变量采用DRP空间离散格式和LDDRK时间积分格式构造PML的无反射预测模型；
S400：在所述模型中添加所述切向速度变量
Figure BDA0002972968990000021
和法向速度变量
Figure BDA0002972968990000022
以计算壳体振动情况下的声场随时间的变化情况；
S500：计算每一时间步长下，舰艇振动壳体运动对声场变量的影响以预测每一时刻舰艇振动壳体的声辐射大小。
所述的方法中，步骤S100中，舰艇振动壳体上均匀布置三向加速度传感器，以采集舰艇振动壳体上多点的切向加速度信号
Figure BDA0002972968990000023
和法向加速度信号
Figure BDA0002972968990000024
所述的方法中，步骤S200中，对振动壳体的切向加速度信号
Figure BDA0002972968990000025
和法向加速度
Figure BDA0002972968990000026
进行时间积分获得振动壳体上多点的切向速度
Figure BDA0002972968990000027
和法向速度
Figure BDA0002972968990000028
Figure BDA0002972968990000029
Figure BDA00029729689900000210
Figure BDA00029729689900000211
其中，
Figure BDA00029729689900000212
是t时刻下的两个切向速度，
Figure BDA00029729689900000213
是t时刻下的法向速度；Δt是振动加速度传感器的采样时间间隔大小，其值为采样频率的倒数Δt＝1/f；f为振动加速度传感器的采样频率。
所述的方法中，所述步骤S300包括：
声波速度和声压的控制方程包括声波速度控制方程和声压的控制方程，其中，
S301：声场变量计算的控制方程：
在x、y、z三维方向上的声波速度控制方程为：
Figure BDA0002972968990000031
Figure BDA0002972968990000032
Figure BDA0002972968990000033
声压为标量，声压的控制方程为：
Figure BDA0002972968990000034
其中，
声波速度为矢量，在x、y、z笛卡尔坐标系三维方向上分别记为ux、uy、uz；声压为标量，在辐射声场中采用p表示，
该控制方程表达为矩阵形式：
Figure BDA0002972968990000035
其中，U为声场变量矩阵，为U＝[ux uy uz p]T；A、B、C为对应的系数矩阵，
S302：采用四阶DRP离散格式对控制方程中的变量进行空间离散，
S303：对控制方程进行LDDRK时间积分：
该控制方程写为矩阵形式，表达式为：
Figure BDA0002972968990000036
其中，U为声场变量矩阵，为U＝[ux uy uz p]T，F为矩阵U的控制方程表达式的简化符号，每一步的积分计算表达式如下，
Ki＝U(t)+βiΔtF(Ki-1)，i＝0，1，…s，其中，s为积分步数，βi为时间积分系数，下一时刻t+Δt的声场速度、声压变量的大小U(t+Δt)，则通过下式计算获得：
U(t+Δt)＝Ks，
当s＝5时，积分所需的系数为：
β1＝0.181575486327091，β2＝0.238260222208392，
β3＝0.330500707328000，β4＝0.500000000000000，β5＝1；
S304：构建无反射的PML边界吸收层：
在声场上下左右前后边界区域上和壳体边界区域内，修正控制方程，修正后的控制方程表达式为：
Figure BDA0002972968990000041
Figure BDA0002972968990000042
其中，q为吸收层中的辅助变量，σx、σy、σz为笛卡尔坐标系三维方向上的吸收系数；
在上下左右前后边界区域中，设置吸收系数大小的计算表达式为：
Figure BDA0002972968990000043
Figure BDA0002972968990000044
Figure BDA0002972968990000045
其中，上标β表示常值指数，进一步地，常值指数为2，σmax为最大吸收系数，x、y、z为声场上下左右前后边界区域的坐标值，xl、yl、zl为吸收层在x、y、z方向上的初始位置，Dl为吸收层的大小，
在壳体边界内部，设置吸收系数大小的计算表达式为：
Figure BDA0002972968990000046
其中，r(x，y，z)表示壳体内部任意空间点与壳体中心轴线间的距离，r0为壳体曲率半径大小，并且rl为辐射声场和边界吸收层间过渡层的大小，rl＝r0-4h，h为空间离散是的网格大小。
所述的方法中，所述步骤S400中，将随时间变化的壳体运动速度边界条件引入所述模型中：
Figure BDA0002972968990000051
Figure BDA0002972968990000052
Figure BDA0002972968990000053
其中，左端为壳体附近流体的切向和法向速度，右端为通过振动壳体上加速度传感时间积分获得的切向和法向壳体速度，其中τ′和τ″表示壳体曲面上的两个切向方向，τ′x、τ′v、τ′z分别为壳体切向向量τ′在x、y、z方向上的分量，τ″x、τ″y、τ″z分别为壳体切向向量τ″在x、y、z方向上的分量，nx、ny、nz分别为壳体法向向量n在x、y、z方向上的分量；
在每一个时间步长内，修正壳体附近流体的速度，使其与壳体运动速度保持一致，通过计算所述模型获得下一时刻的声场变量大小。
所述的方法中，所述步骤S500中，通过不断迭代添加壳体运动速度边界条件，并且计算振动壳体声辐射模型实现振动壳体辐射噪声时域预测。
所述的方法中，声场内部点的空间离散中，ux沿着x方向空间离散为：
Figure BDA0002972968990000054
其中，j为整数，h为声场网格间距大小，aj为对应的系数：
a0＝0，
a1＝-a-1＝0.77088238051822552，
a2＝-a-2＝-0.166705904414580469，
a3＝-a-3＝0.02084314277031176；
声场边界点的空间离散中，以ux沿着x方向空间离散为：
Figure BDA0002972968990000061
其中，j为整数，h为声场网格间距大小，
Figure BDA0002972968990000062
为对应系数，n为目标节点的左边节点数，m为目标节点的右边节点数，对应左边界节点
Figure BDA0002972968990000063
和右边界节点
Figure BDA0002972968990000064
的系数：
Figure BDA0002972968990000065
对应次左边界节点
Figure BDA0002972968990000066
和次右边界节点
Figure BDA0002972968990000067
的系数：
Figure BDA0002972968990000068
对应次次左边界节点
Figure BDA0002972968990000069
和次次右边界节点
Figure BDA00029729689900000610
的系数：
Figure BDA00029729689900000611
在上述技术方案中，本发明提供的一种舰艇振动壳体的辐射噪声预测方法，具有以下有益效果：舰艇振动壳体的辐射噪声预测方法能够实现辐射噪声的精确预测。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明中记载的一些实施例，对于本领域普通技术人员来讲，还可以根据这些附图获得其他的附图。
图1是本发明一个实施例的舰艇振动壳体辐射噪声预测方法的步骤示意图；
图2是本发明一个实施例的脉动圆柱壳体辐射噪声计算示意图；
图3是本发明一个实施例的加速度传感器布置图；
图4(a)至图4(c)是本发明一个实施例的不同时刻下声压登高线图；
图5(a)至图5(c)是本发明一个实施例计算得到的不同时刻下径向声压分布情况图。
具体实施方式
为使本发明实施方式的目的、技术方案和优点更加清楚，下面将结合本发明实施方式中的附图1至图5(c)，对本发明实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式是本发明一部分实施方式，而不是全部的实施方式。基于本发明中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都属于本发明保护的范围。
因此，以下对在附图中提供的本发明的实施方式的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施方式。基于本发明中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都属于本发明保护的范围。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”、“顺时针”、“逆时针”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的设备或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中，“多个”的含义是两个或两个以上，除非另有明确具体的限定。
在本发明中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或成一体；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通或两个元件的相互作用关系。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
在本发明中，除非另有明确的规定和限定，第一特征在第二特征之“上”或之“下”可以包括第一和第二特征直接接触，也可以包括第一和第二特征不是直接接触而是通过它们之间的另外的特征接触。而且，第一特征在第二特征“之上”、“上方”和“上面”包括第一特征在第二特征正上方和斜上方，或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”包括第一特征在第二特征正下方和斜下方，或仅仅表示第一特征水平高度小于第二特征。
为了使本领域的技术人员更好地理解本发明的技术方案，下面将结合附图对本发明作进一步的详细介绍。
为了更好地理解，图1为舰艇振动壳体辐射噪声预测方法的步骤示意图；如图1所示，舰艇振动壳体辐射噪声预测方法包括以下步骤：
S100：采集舰艇振动壳体的振动信号，得到壳体的切向加速度信号
Figure BDA0002972968990000081
和法向加速度信号
Figure BDA0002972968990000082
其中X(t)表示为壳体边界的位置；
S200：对所述切向加速度信号
Figure BDA0002972968990000083
和法向加速度信号
Figure BDA0002972968990000084
进行时间积分得到切向速度
Figure BDA0002972968990000085
和法向速度
Figure BDA0002972968990000086
S300：给出了声波速度和声压的控制方程，并对控制方程进行空间离散化和时间积分，从而构造处振动壳体噪声预测的计算模型；在该模型中计算的声场变量有声压大小p(x，t)、声波切向速度uτ′(x，t)和uτ″(x，t)、声波法向速度un(x，t)，x表示为声辐射场的不同位置，针对这些变量采用DRP空间离散格式和LDDRK时间积分格式，在此基础上构造了PML的无反射吸收边界层，从而实现振动壳体声辐射模型的构建；声波速度和声压的控制方程包括声波速度控制方程和声压的控制方程，其中，声波变量控制方程也称之为线性欧拉方程。舰艇振动壳体辐射噪声预测模型是采用DRP空间离散格式，LDDRK时间积分格式，以及PML吸收层计算上述声波变量控制方程的集合，具体表现形式为S300中所有的计算公式。其中DRP空间离散格式，LDDRK时间积分格式，以及PML吸收层分别为色散关系守恒(Dispersion-Relation-Preserving，DRP)的空间离散格式，低色散低耗散Runge-Kutta法(LowDissipation and low-Dispersion Runge-Kutta algorithm，LDDRK法)的时间积分方案，以及完美匹配层技术(Perfect Matched Layer，PML技术)对边界层内的干涉声波进行吸收，故而称之PML吸收层。
S400：随着时间变化，在振动壳体声辐射模型中添加壳体切向速度变量
Figure BDA0002972968990000091
Figure BDA0002972968990000092
和法向速度变量
Figure BDA0002972968990000093
以此计算壳体振动情况下，声场随时间的变化情况；
S500：计算每一时间步长下，振动壳体运动对声场变量的影响，从而预测每一时刻振动壳体的声辐射大小。
上述实施例构成了本发明的完整技术方案，与现有技术不同，上述实施例采用了已知的振动壳体加速度信号，利用了高精度的空间离散格式和时间积分格式，以及消除边界反射效应的边界吸收层，从而构造了振动壳体辐射噪声预测模型。在此基础上添加随时间变化的壳体运动速度边界条件，以此实现辐射噪声的精确预测。
图2为脉动圆柱壳体辐射噪声计算示意图，其中脉动运动的圆柱壳体位于声场计算域的中心位置，振动壳体的脉动运动速度表达式为ur＝yreiωt＝yrcos(ωt)，其中ω为脉动运动角频率，大小为31.4159rad/s。测试条件为：圆柱壳体的直径大小为D＝1m，计算域的大小为12m，声传播介质密度为1Kg/m3，声波的传播速度为1m/s，振动壳体的脉动运动幅值为1×10-4m/s，计算的时间步长大小为0.005s。
在本实施例中，步骤S100中，所述切向和法向上的振动加速度信号
Figure BDA0002972968990000094
和
Figure BDA0002972968990000095
是由三向加速度传感器采集获得。图3是振动圆柱壳体加速度传感器的布置示意图。从图3中可以看到，在每一层截面和竖直方向上，传感器均布置，以便构造振动圆柱壳体的实际运动速度情况。
本实施例中，步骤S200中，对振动壳体的切向加速度
Figure BDA0002972968990000096
和法向加速度
Figure BDA0002972968990000097
进行时间积分获得振动壳体上多点的切向速度
Figure BDA0002972968990000098
和法向速度
Figure BDA0002972968990000101
Figure BDA0002972968990000102
Figure BDA0002972968990000103
Figure BDA0002972968990000104
其中，
Figure BDA0002972968990000105
是t时刻下的两个切向速度，
Figure BDA0002972968990000106
是t时刻下的法向速度；Δt是振动加速度传感器的采样时间间隔大小，其值为采样频率的倒数Δt＝1/f；f为加速度传感器的采样频率。
本实施例中，所述步骤S300包括：
S301：声场变量计算的控制方程：
在x、y、z三维方向上的声波速度控制方程为：
Figure BDA0002972968990000107
Figure BDA0002972968990000108
Figure BDA0002972968990000109
声压为标量，故而声压的控制方程为：
Figure BDA00029729689900001010
其中声波速度为矢量，在x、y、z三维方向上分别记为ux、uy、uz；声压为标量，在辐射声场中采用p表示。
该控制方程可表达为矩阵形式：
Figure BDA00029729689900001011
其中，U为声场变量矩阵，即为U＝[ux uy uz p]T；A、B、C为对应的系数矩阵。
S302：对控制方程中的变量进行空间离散，采用四阶DRP离散格式：
声场内部点的空间离散中，ux沿着x方向空间离散为：
Figure BDA0002972968990000111
其中，j为整数，h为声场网格间距大小，aj为对应的系数：
a0＝0，
a1＝-a-1＝0.77088238051822552，
a2＝-a-2＝-0.166705904414580469，
a3＝-a-3＝0.02084314277031176；
声场边界点的空间离散中，以ux沿着x方向空间离散为：
Figure BDA0002972968990000112
其中，j为整数，h为声场网格间距大小，
Figure BDA0002972968990000113
为对应系数，n为目标节点的左边节点数，m为目标节点的右边节点数，对应左边界节点
Figure BDA0002972968990000114
和右边界节点
Figure BDA0002972968990000115
的系数：
Figure BDA0002972968990000116
对应次左边界节点
Figure BDA0002972968990000117
和次右边界节点
Figure BDA0002972968990000118
的系数：
Figure BDA0002972968990000119
对应次次左边界节点
Figure BDA0002972968990000121
和次次右边界节点
Figure BDA0002972968990000122
的系数：
Figure BDA0002972968990000123
S303：对控制方程进行LDDRK时间积分方案：
该控制方程可写为矩阵形式，表达式为：
Figure BDA0002972968990000124
其中，U包含声波速度、声压等变量的矩阵。四阶LDDRK时间积分方案中，每一步的积分计算表达式如下
Ki＝U(t)+βiΔtF(Ki-1)，i＝0，1，…s，其中，s为积分步数，βi为时间积分系数，下一时刻t+Δt的声场速度、声压变量的大小U(t+Δt)，则通过下式计算获得：
U(t+Δt)＝Ks，
当s＝5时，积分所需的系数为：
β1＝0.181575486327091，β2＝0.238260222208392，
β3＝0.330500707328000，β4＝0.500000000000000，β5＝1；
S304：构建无反射的PML边界吸收层：
其本质为在声场上下左右前后边界区域上和壳体边界区域内，修正控制方程，添加吸收系数，从而实现边界的无反射效果。修正后的控制方程表达式为：
Figure BDA0002972968990000125
Figure BDA0002972968990000126
其中，q为吸收层中的辅助变量，σx、σy、σz为吸收系数；
在上下左右前后边界区域中，设置吸收系数大小的计算表达式为：
Figure BDA0002972968990000131
Figure BDA0002972968990000132
Figure BDA0002972968990000133
其中，上标β表示常值指数，一般取值为2，σmax为最大吸收系数，x、y、z为声场上下左右前后边界区域的坐标值，xl、yl、zl为吸收层在x、y、z方向上的初始位置，Dl为吸收层的大小。
在壳体边界内部，设置吸收系数大小的计算表达式为：
Figure BDA0002972968990000134
其中，r(x，y，z)表示壳体内部任意空间点与壳体中心轴线间的距离，r0为壳体曲率半径大小，并且rl为辐射声场和边界吸收层间过渡层的大小，rl＝r0-4h。
本实施例中，步骤S400中，将随时间变化的壳体运动速度边界条件引入至声辐射计算模型中：
Figure BDA0002972968990000135
Figure BDA0002972968990000136
Figure BDA0002972968990000137
其中，左端为壳体附近流体的切向和法向速度，右端为通过振动壳体上加速度传感时间积分获得的切向和法向壳体速度，其中τ′和τ″表示壳体曲面上的两个切向方向，τ′x、τ′y、τ′z分别为壳体切向向量τ′在x、y、z方向上的分量，τ″x、τ″y、τ″z分别为壳体切向向量τ″在x、y、z方向上的分量，nx、ny、nz分别为壳体法向向量n在x、y、z方向上的分量。在每一个时间步长内，修正壳体附近流体的速度，使其与壳体运动速度保持一致。
在添加壳体运动速度边界条件的基础上，通过计算振动壳体声辐射模型，可获得下一时刻的声场变量大小，即预测了下一时刻辐射噪声；通过不断迭代添加壳体运动速度边界条件，可实现振动壳体辐射噪声时域预测。
本实施例中，步骤S500中，通过不断迭代添加壳体运动速度边界条件，并且计算振动壳体声辐射模型，可实现振动壳体辐射噪声时域预测。
图4(a)至图4(c)和图5(a)至图5(c)是本发明一个实施例通过上述步骤计算得到的，分别为t＝2、3、4s时的径向声压分布情况图，以及t＝2、3、4s时的声压等高线图。在图4(a)至图4(c)中给出了该实例的解析解，从中可以发现本发明方法计算结果与解析结果完全吻合，说明本发明方法的准确性。从图5(a)至图5(c)中可以看出，在不同时刻下，在圆柱内部的声压被完美吸收了，不会影响脉冲振动壳体带来的向外扩散的辐射噪声，可以获得干净理想的计算结果。
本方法采集壳体切向和法向加速度振动信号；计算振动壳体的切向和法向速度变量大小；通过色散关系守恒(Dispersion-Relation-Preserving，DRP)的空间离散格式，4阶6步的低色散低耗散Runge-Kutta法(Low Dissipation and low-Dispersion Runge-Kuttaalgorithm，LDDRK法)进行时间积分和PML吸收层，构建基于线性欧拉方程的振动壳体辐射噪声预测模型；随时间变化，添加对应的边界条件，实现振动壳体辐射噪声的连续预测。
最后应该说明的是：所描述的实施例仅是本申请一部分实施例，而不是全部的实施例，基于本申请中的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其它实施例，都属于本申请保护的范围。
以上只通过说明的方式描述了本发明的某些示范性实施例，毋庸置疑，对于本领域的普通技术人员，在不偏离本发明的精神和范围的情况下，可以用各种不同的方式对所描述的实施例进行修正。因此，上述附图和描述在本质上是说明性的，不应理解为对本发明权利要求保护范围的限制。

Claims (6)

1.一种舰艇振动壳体的辐射噪声预测方法，其特征在于，其包括以下步骤：
S100：采集舰艇振动壳体的振动信号，得到壳体的切向加速度信号
Figure DEST_PATH_IMAGE001
、
Figure 287594DEST_PATH_IMAGE002
和法向加速度信号
Figure DEST_PATH_IMAGE003
，其中
Figure 774070DEST_PATH_IMAGE004
表示为壳体边界的位置;
S200：对所述切向加速度信号
Figure DEST_PATH_IMAGE005
、
Figure 997198DEST_PATH_IMAGE006
和法向加速度信号
Figure DEST_PATH_IMAGE007
进行时间积分得到切向速度
Figure 670756DEST_PATH_IMAGE008
、
Figure DEST_PATH_IMAGE009
和法向速度
Figure 541760DEST_PATH_IMAGE010
；
S300：基于声波速度和声压的控制方程空间离散化和时间积分，构造舰艇振动壳体辐射噪声预测模型；
S400：在所述模型中添加所述切向速度变量
Figure DEST_PATH_IMAGE011
、
Figure 617163DEST_PATH_IMAGE012
和法向速度变量
Figure DEST_PATH_IMAGE013
以计算壳体振动情况下的声场随时间的变化情况；
S500：计算每一时间步长下，舰艇振动壳体运动对声场变量的影响以预测每一时刻舰艇振动壳体的声辐射大小；
其中，
对于所述步骤S300的模型，其中：
声场变量包括：声压大小
Figure 17052DEST_PATH_IMAGE014
、声波切向速度
Figure DEST_PATH_IMAGE015
和
Figure 177906DEST_PATH_IMAGE016
、声波法向速度
Figure DEST_PATH_IMAGE017
，其中，
Figure 321442DEST_PATH_IMAGE018
和
Figure DEST_PATH_IMAGE019
以及
Figure 251352DEST_PATH_IMAGE020
中的
Figure DEST_PATH_IMAGE021
表示声辐射场的不同位置，针对所述声场变量采用DRP空间离散格式和LDDRK时间积分格式构造PML的无反射预测模型；
声波速度和声压的控制方程包括声波速度控制方程和声压的控制方程，其中，步骤S300进一步包括如下子步骤：
S301: 获得声场变量计算的控制方程：
在x、y、z三维方向上的声波速度控制方程为：
Figure 87721DEST_PATH_IMAGE022
，
声压为标量，声压的控制方程为：
Figure DEST_PATH_IMAGE023
，其中，
声波速度为矢量，在x、y、z笛卡尔坐标系三维方向上分别记为
Figure 735871DEST_PATH_IMAGE024
、
Figure DEST_PATH_IMAGE025
、
Figure 683099DEST_PATH_IMAGE026
；声压为标量，在辐射声场中采用
Figure DEST_PATH_IMAGE027
表示，
该控制方程表达为矩阵形式：
Figure 453386DEST_PATH_IMAGE028
其中，
Figure DEST_PATH_IMAGE029
为声场变量矩阵，为
Figure 382027DEST_PATH_IMAGE030
；
Figure DEST_PATH_IMAGE031
、
Figure 704424DEST_PATH_IMAGE032
、
Figure DEST_PATH_IMAGE033
为对应的系数矩阵，
S302：采用四阶DRP离散格式对控制方程中的变量进行空间离散，
S303：对控制方程进行LDDRK时间积分：
该控制方程写为矩阵形式，表达式为：
Figure 376714DEST_PATH_IMAGE034
，其中，
Figure 405850DEST_PATH_IMAGE029
为声场变量矩阵，为
Figure DEST_PATH_IMAGE035
，
Figure 770972DEST_PATH_IMAGE036
为矩阵
Figure 518348DEST_PATH_IMAGE029
的控制方程表达式的简化符号,每一步的积分计算表达式如下，
Figure DEST_PATH_IMAGE037
，其中，s为积分步数，
Figure 994329DEST_PATH_IMAGE038
为时间积分系数，下一时刻
Figure DEST_PATH_IMAGE039
的声场速度、声压变量的大小
Figure 940288DEST_PATH_IMAGE040
，则通过下式计算获得：
Figure DEST_PATH_IMAGE041
，
当s=5时，积分所需的系数为：
Figure 210733DEST_PATH_IMAGE042
，
Figure DEST_PATH_IMAGE043
，
Figure 242142DEST_PATH_IMAGE044
，
Figure DEST_PATH_IMAGE045
，
Figure 990656DEST_PATH_IMAGE046
；
S304：构建无反射的PML边界吸收层：
在声场上下左右前后边界区域上和壳体边界区域内，修正控制方程，修正后的控制方程表达式为：
Figure DEST_PATH_IMAGE047
，
Figure 62560DEST_PATH_IMAGE048
，其中，
Figure DEST_PATH_IMAGE049
为吸收层中的辅助变量，
Figure 769485DEST_PATH_IMAGE050
、
Figure DEST_PATH_IMAGE051
、
Figure 288191DEST_PATH_IMAGE052
为笛卡尔坐标系三维方向上的吸收系数；
在上下左右前后边界区域中，设置吸收系数大小的计算表达式为：
Figure DEST_PATH_IMAGE053
，
Figure 840395DEST_PATH_IMAGE054
，
Figure DEST_PATH_IMAGE055
，
其中，上标β表示常值指数，为
Figure 495367DEST_PATH_IMAGE056
最大吸收系数，x、y、z为声场上下左右前后边界区域的坐标值，
Figure DEST_PATH_IMAGE057
、
Figure 373193DEST_PATH_IMAGE058
、
Figure DEST_PATH_IMAGE059
为吸收层在x、y、z方向上的初始位置，
Figure 379196DEST_PATH_IMAGE060
为吸收层的大小，
在壳体边界内部，设置吸收系数大小的计算表达式为：
Figure DEST_PATH_IMAGE061
，
其中，
Figure 670DEST_PATH_IMAGE062
表示壳体内部任意空间点与壳体中心轴线间的距离，
Figure DEST_PATH_IMAGE063
为壳体曲率半径大小，并且
Figure 775728DEST_PATH_IMAGE064
为辐射声场和边界吸收层间过渡层的大小，
Figure DEST_PATH_IMAGE065
, h为空间离散时的网格大小。
2.根据权利要求1所述的方法，其特征在于，步骤S100中，舰艇振动壳体上均匀布置三向加速度传感器，以采集舰艇振动壳体上多点的切向加速度信号
Figure 293297DEST_PATH_IMAGE066
、
Figure DEST_PATH_IMAGE067
和法向加速度信号
Figure 786595DEST_PATH_IMAGE068
。
3.根据权利要求1所述的方法，其特征在于，步骤S200中，对振动壳体的切向加速度信号
Figure DEST_PATH_IMAGE069
、
Figure 940321DEST_PATH_IMAGE070
和法向加速度
Figure DEST_PATH_IMAGE071
，进行时间积分获得振动壳体上多点的切向速度
Figure 569886DEST_PATH_IMAGE072
、
Figure DEST_PATH_IMAGE073
和法向速度
Figure 523935DEST_PATH_IMAGE074
：
Figure DEST_PATH_IMAGE075
，
Figure 504530DEST_PATH_IMAGE076
Figure DEST_PATH_IMAGE077
，
其中，
Figure 467807DEST_PATH_IMAGE078
、
Figure DEST_PATH_IMAGE079
是
Figure 686298DEST_PATH_IMAGE080
时刻下的两个切向速度，
Figure DEST_PATH_IMAGE081
是
Figure 76828DEST_PATH_IMAGE080
时刻下的法向速度；
Figure 482402DEST_PATH_IMAGE082
是振动加速度传感器的采样时间间隔大小，其值为采样频率的倒数
Figure DEST_PATH_IMAGE083
；f为振动加速度传感器的采样频率。
4.根据权利要求1所述的方法，其特征在于，所述步骤S400中，将随时间变化的壳体运动速度边界条件引入所述模型中：
Figure 249370DEST_PATH_IMAGE084
，
Figure DEST_PATH_IMAGE085
，
Figure 853526DEST_PATH_IMAGE086
，
其中，左端为壳体附近流体的切向和法向速度，右端为通过振动壳体上加速度传感时间积分获得的切向和法向壳体速度，其中
Figure DEST_PATH_IMAGE087
和
Figure 414958DEST_PATH_IMAGE088
表示壳体曲面上的两个切向方向，
Figure DEST_PATH_IMAGE089
、
Figure 838986DEST_PATH_IMAGE090
、
Figure DEST_PATH_IMAGE091
分别为壳体切向向量
Figure 149925DEST_PATH_IMAGE092
在x、y、z方向上的分量，
Figure DEST_PATH_IMAGE093
、
Figure 139746DEST_PATH_IMAGE094
、
Figure DEST_PATH_IMAGE095
分别为壳体切向向量
Figure 872079DEST_PATH_IMAGE096
在x、y、z方向上的分量，
Figure DEST_PATH_IMAGE097
、
Figure 48982DEST_PATH_IMAGE098
、
Figure DEST_PATH_IMAGE099
分别为壳体法向向量
Figure 892173DEST_PATH_IMAGE100
在x、y、z方向上的分量；
在每一个时间步长内，修正壳体附近流体的速度，使其与壳体运动速度保持一致，通过计算所述模型获得下一时刻的声场变量大小。
5.根据权利要求1所述的方法，其特征在于，所述步骤S500中，通过不断迭代添加壳体运动速度边界条件，并且计算振动壳体声辐射模型实现振动壳体辐射噪声时域预测。
6.根据权利要求1所述的方法，其特征在于，声场内部点的空间离散中，
Figure 939764DEST_PATH_IMAGE024
沿着x方向空间离散为：
Figure DEST_PATH_IMAGE101
其中，j为整数，h为声场网格间距大小，
Figure 842998DEST_PATH_IMAGE102
为对应的系数：
Figure DEST_PATH_IMAGE103
，
Figure 507197DEST_PATH_IMAGE104
，
Figure DEST_PATH_IMAGE105
，
Figure 419658DEST_PATH_IMAGE106
；
声场边界点的空间离散中，以
Figure DEST_PATH_IMAGE107
沿着x方向空间离散为：
Figure 852914DEST_PATH_IMAGE108
，
其中，j为整数，h为声场网格间距大小，
Figure DEST_PATH_IMAGE109
为对应系数，n为目标节点的左边节点数，m为目标节点的右边节点数，对应左边界节点
Figure 210206DEST_PATH_IMAGE110
和右边界节点
Figure DEST_PATH_IMAGE111
的系数：
Figure 96123DEST_PATH_IMAGE112
，
Figure DEST_PATH_IMAGE113
，
Figure 812275DEST_PATH_IMAGE114
，
Figure DEST_PATH_IMAGE115
，
Figure 365616DEST_PATH_IMAGE116
，
Figure DEST_PATH_IMAGE117
，
Figure 610653DEST_PATH_IMAGE118
；
对应次左边界节点
Figure DEST_PATH_IMAGE119
和次右边界节点
Figure 515024DEST_PATH_IMAGE120
的系数：
Figure DEST_PATH_IMAGE121
，
Figure 238129DEST_PATH_IMAGE122
，
Figure DEST_PATH_IMAGE123
，
Figure 380397DEST_PATH_IMAGE124
，
Figure DEST_PATH_IMAGE125
，
Figure 61914DEST_PATH_IMAGE126
，
Figure DEST_PATH_IMAGE127
；
对应次次左边界节点
Figure 922423DEST_PATH_IMAGE128
和次次右边界节点
Figure DEST_PATH_IMAGE129
的系数：
Figure 455079DEST_PATH_IMAGE130
，
Figure DEST_PATH_IMAGE131
，
Figure 186274DEST_PATH_IMAGE132
，
Figure DEST_PATH_IMAGE133
，
Figure 38693DEST_PATH_IMAGE134
，
Figure DEST_PATH_IMAGE135
，
Figure 386497DEST_PATH_IMAGE136
。
CN202110273677.2A 2021-03-12 2021-03-12 舰艇振动壳体的辐射噪声预测方法 Active CN113190912B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110273677.2A CN113190912B (zh) 2021-03-12 2021-03-12 舰艇振动壳体的辐射噪声预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110273677.2A CN113190912B (zh) 2021-03-12 2021-03-12 舰艇振动壳体的辐射噪声预测方法

Publications (2)

Publication Number Publication Date
CN113190912A CN113190912A (zh) 2021-07-30
CN113190912B true CN113190912B (zh) 2022-10-28

Family

ID=76973348

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110273677.2A Active CN113190912B (zh) 2021-03-12 2021-03-12 舰艇振动壳体的辐射噪声预测方法

Country Status (1)

Country Link
CN (1) CN113190912B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109269497A (zh) * 2018-07-31 2019-01-25 哈尔滨工程大学 基于auv切法向速度模型的多尺度无迹卡尔曼滤波估计方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3686871B2 (ja) * 2002-02-27 2005-08-24 三ツ星ベルト株式会社 歯付ベルトの騒音予測方法及び装置、並びにプログラム
US10296683B2 (en) * 2013-11-06 2019-05-21 Thornton Tomasetti, Inc. Computer implemented apparatus and method for finite element modeling using hybrid absorbing element

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109269497A (zh) * 2018-07-31 2019-01-25 哈尔滨工程大学 基于auv切法向速度模型的多尺度无迹卡尔曼滤波估计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
The fast prediction of sound radiation from typical submarine cabin based on acoustic transfer vector;Qifan Zhang等;《IEEE Xplore》;20160808;全文 *
水下圆柱壳结构表面振速插值方法及其辐射声场预测;吴海军等;《声学学报》;20150915(第05期);全文 *

Also Published As

Publication number Publication date
CN113190912A (zh) 2021-07-30

Similar Documents

Publication Publication Date Title
CN108763683B (zh) 一种三角函数框架下新weno格式构造方法
CN111307266B (zh) 声速获取方法及基于该方法的全球海洋声速场构建方法
JP4269623B2 (ja) 血流可視化診断装置
CN109446634B (zh) 基于泰勒展开边界元方法的船舶运动预报方法
CN104678384B (zh) 一种波束域的声压差互相关谱分析水下目标速度估计方法
He et al. A non-contact acoustic pressure-based method for load identification in acoustic-structural interaction system with non-probabilistic uncertainty
CN115630446B (zh) 用于潜航器结构低频辐射噪声实时仿真的快速推演方法
CN112257311B (zh) Pekeris波导下结构声振计算的FE/WSM方法
CN112577592B (zh) 基于空间傅里叶变换的有限空间平面近场声全息测量方法
CN109489796A (zh) 一种基于单元辐射法的水下复杂结构辐射噪声源定位识别与声辐射预报方法
CN113190912B (zh) 舰艇振动壳体的辐射噪声预测方法
CN115906715B (zh) 一种波浪作用下可压缩淤泥质海床土体运动速度计算方法及系统
CN105487047A (zh) 一种高参数减压阀噪声源位置及噪声指向性的确定方法
Goody et al. Separated flow surface pressure fluctuations and pressure-velocity correlations on prolate spheroid
Zhang et al. A gradient vector descent strategy for localizing acoustic emission sources in discontinuous structures with a hole
Dong et al. Simplified calculation method for hydrodynamic noise of multi-tube heat exchanger based on acoustic cavity model
CN115336432B (zh) 一种用于大型复杂结构水下辐射噪声的快速预报方法
CN117606611A (zh) 一种混合激励下的舱室结构水下辐射噪声预报方法
CN113343499B (zh) 一种舰艇壳体的流噪声预测方法
Yang et al. The Influence of Vertical Cable on Flow Field and Acoustic Analysis of A Submersible Buoy System Based on CFD
Stephen et al. Natural frequency error estimation using a patch recovery technique
Che et al. Real-time array shape estimation method of horizontal suspended linear array based on non-acoustic auxiliary sensors
CN118036387A (zh) 基于有限元及边界元法的不同潜深状态下声辐射预报方法
CN117782299A (zh) 一种基于声全息的直升机声足迹高效预测方法
CN117010293A (zh) 基于水下航行器的流固声耦合数值分析方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant