[发明专利]一种水下声辐射预报方法、系统、计算机设备、存储介质

权利要求书

1.一种水下声辐射预报方法，其特征在于，所述水下声辐射预报方法包括：基于径向点插值法和改进的Dirichlet-to-Neumann边界条件进行纯无网格弱式的水下结构声辐射预报。

2.如权利要求1所述水下声辐射预报方法，其特征在于，所述水下声辐射预报方法包括以下步骤：

步骤一，采用人工边界截断结构外部无限问题域获得有限计算域；

步骤二，在RPIM-MDtN中，采用RPIM构造声形函数，在人工边界处施加MDtN边界条件，进行水下声辐射预报。

3.如权利要求2所述水下声辐射预报方法，其特征在于，所述采用RPIM构造声形函数包括：

N^T(x)＝[N₁(x),N₂(x),…,N_n(x)]；

其中，N^T(x)表示采用RPIM构造的n个场节点的形函数；p_s＝[p₁,p₂,…,p_n]^T表示声压值组成的向量；p^h(x)表示声压p(x)的试函数；p^T＝[p₁(x),p₂(x),…,p_m(x)]。

4.如权利要求2所述水下声辐射预报方法，其特征在于，所述在RPIM-MDtN中，采用RPIM构造声形函数，在人工边界处施加MDtN边界条件，进行水下声辐射预报还包括：

在加权残差公式中，采用RPIM构造的形函数N(x)作为检验函数W(x)，则RPIM构造声形函数表示为：

建立如下矩阵形式的离散系统方程

[K-k²M+K^b]{p}＝F；

式中：{p}表示声压值的矢量，并且C_n(x)＝[cosnθsinnθ]；

基于得到的人工边界处的声压值，利用下式计算得到V_O中任意场点处的声压值：

其中，H_n⁽¹⁾表示第一类n阶的汉克尔函数，p(R,θ)表示Γ_R上的声压，求和符号后的“'”表示对应n＝0的项乘0.5。

5.如权利要求2所述水下声辐射预报方法，其特征在于，所述MDtN边界条件包括：

其中，

其中，上“'”表示对kR求导。

6.一种接收用户输入程序存储介质，所存储的计算机程序使电子设备执行权利要求1-5任意一项所述水下声辐射预报方法包括下列步骤：

步骤一，采用人工边界截断结构外部无限问题域获得有限计算域；

步骤二，在RPIM-MDtN中，采用RPIM构造声形函数，在人工边界处施加MDtN边界条件，进行水下声辐射预报。

7.一种存储在计算机可读介质上的计算机程序产品，包括计算机可读程序，供于电子装置上执行时，提供用户输入接口以实施如权利要求1-5任意一项所述水下声辐射预报方法。

8.一种实施权利要求1-5任意一项所述水下声辐射预报方法的水下声辐射预报系统，其特征在于，所述水下声辐射预报系统包括：

有限计算域获得模块，用于采用人工边界截断结构外部无限问题域获得有限计算域；

水下声辐射预报模块，用于在RPIM-MDtN中，采用RPIM构造声形函数，在人工边界处施加MDtN边界条件，进行水下声辐射预报。