CN112149226A - 一种基于局部无网格基本解法的车内噪声预测方法 - Google Patents

Abstract

本发明公开了一种基于局部无网格基本解法的车内噪声预测方法，包括以下步骤：首先，建立车腔的几何模型；其次，对所建立的几何模型进行布点，根据实际工况配置边界数据；然后，建立局部无网格基本解法数值模型，计算车内声压和声压级的分布；最后，输出车内声压和声压级的预测值，并对预测结果进行性能评估。本发明具有数学理论简单，无需网格划分和数值积分，处理速度快，预测精度高，易于大尺度和高频声波预测等优点。本发明可用于分析不同车型不同工况下车内噪声的分布规律，能够缩短车辆测试调校周期，降低车辆NVH研发成本。

Description

一种基于局部无网格基本解法的车内噪声预测方法

技术领域：

本发明属于汽车车内噪声预测领域，涉及一种汽车噪音的预测方法，特别是一种基于局部无网格基本解法的车内噪声预测方法。

背景技术：

随着科学技术的快速发展和人民生活水平的不断提高，汽车已经成为每个家庭的必备代步工具，中国也成为了全球最大的汽车消费国。伴随着国外汽车品牌战略重心向中低档车型市场转移及《中国制造2025》的提出，我国要实现从汽车大国向汽车强国转变，就必须致力于汽车质量和性能的提高。据世界知名汽车质量评估机构J.D.Power统计，汽车整体性能指标中有近三分之一与汽车的振动、噪声及声振粗糙度(NVH)直接相关。汽车车内噪声是影响乘坐舒适性的重要因素之一。高噪声既能损害乘员的身心健康，又能导致驾驶员迅速疲劳，对汽车行驶安全构成了极大的威胁。消费者在购买汽车时往往会首先感受乘坐或驾驶汽车时的振动和噪声状态。因此，汽车车内噪声预测与控制方法的研究对汽车的设计和制造具有重要的理论意义和实用价值。

近年来，随着计算机软硬件技术的日臻成熟与不断完善，CAD/CAM/CAE技术已在汽车研发生产中起到重要的作用，特别是汽车的振动噪声分析及预测。作为商用CAE软件的核心模块之一，以有限元和边界元为主体的数值计算方法对于仿真的精度和效率而言重要性可见一斑。然而，有限元网格划分成本过高，复杂几何模型边界离散误差过大。边界元方法将问题降一维求解，无需划分区域网格，但是必须进行数值积分，特别是近奇异和超奇异积分的计算好坏直接关系到最终结果的可靠性。因此，对于复杂几何模型和高频噪声分析等问题，他们的计算结果与实际结果之间仍存在着较大差距，因而还不能很好地指导汽车的振动、噪声开发设计。

作为一类新型数值计算方法，近年来发展起来的局部无网格基本解算法具有无需网格剖分和数值积分，数学理论简单，易于编程，计算精度高，适用于大尺度复杂区域问题等优点(见文献1.W.Qu,C.M.Fan,Y.Gu.Localized method of fundamental solutionsfor interior Helmholtz problems with high wave number[J].Engineering Analysiswith Boundary Elements,2019,107:25-32.文献2.W.Qu,C.M.Fan,Y.Gu,et al.Analysisof three-dimensional interior acoustic fields by using the localized methodof fundamental solutions[J].Applied Mathematical Modelling,2019,76:122-132.)。克服了传统有限元和边界元方法的自身不足，在计算力学和数值仿真中具有重要的应用前景，也为准确高效地进行汽车噪声预测提供了可能，因此，本发明寻求设计提供一种基于局部无网格基本解法的车内噪声预测方法，准确、快速、有效地预测汽车车内噪声。

发明内容：

本发明的目的在于针对现有技术中存在的上述问题，设计提供一种基于局部无网格基本解法的车内噪声预测方法，能准确、快速、有效地预测汽车车内噪声。

为了实现上述目的，本发明涉及的基于局部无网格基本解法的车内噪声预测方法通过如下技术方案实现：

一种基于局部无网格基本解法的车内噪声预测方法，具体操作步骤按照如下方式进行：

S1、根据需预测车辆的几何参数，建立车腔的几何模型，并在车腔的边界及内部配置节点；

S2、利用车腔内空气的密度、声音在空气中的传播速度等已知参数，结合相应的边界条件，建立车内声学分析的局部无网格基本解法数值模型；

S3、利用步骤S2所建立的数值模型，计算步骤S1中所配置节点上的声压和声压级，对车内噪声进行预测；

S4、输出待预测汽车车腔内的噪声分布，并对预测结果进行性能评估；

进一步的，步骤S1中车腔的边界及内部配置节点包括：

在车腔边界上配置n_b个节x⁽ⁱ⁾,i＝1,2,...,n_b；

在车腔内部配置n_i个节x⁽ⁱ⁾,i＝1,2,...,n_i。

进一步的，步骤S2中所采用的声学控制方程为：

其中

为拉普拉斯算子，Ω表示需预测车内区域，x为空间坐标，p(x)为x处的声压幅值，ω为角频率，c为声速。

边界条件为下列三种中的一种或多种：

p＝p_D,在Γ_D上

在Γ_N上

在Γ_R上

其中

ρ为空气的密度，n为边界点上的单位外法向量，p_D为本质边界Γ_D上的声压，v_n为自然边界Γ_N上的法向速度，A_n为吸收边界Γ_R上的阻尼系数；

步骤S2中建立局部无网格基本解法数值模型步骤包括：

(1)对任意节点x⁽ⁱ⁾，搜索其m个邻近节点

确定其所对应的局部子区域

并在局部子区域外一定距离d处配置相应数量的源点

(2)构建如下形式的局部矩阵H，并求其逆矩阵H^-1；

其中，G(r_kj)为声学控制方程的基本解，

为局部节点与源点之间的距离；

(3)利用给定边界条件建立相应的代数方程：对于内部节点，满足下列方程

本质边界节点，满足下列方程

p(x⁽ⁱ⁾)＝p_D

对于自然边界节点，满足下列方程

对于吸收边界节点，满足下列方程

上式中，

h⁽ⁱ⁾＝[G(r_i1),G(r_i2),...,G(r_im)]，

(4)联立(3)中的各方程，获得如下的最终稀疏矩阵方程；

Ap＝b

其中，A为(3)中联立方程所形成的系数矩阵，b为(3)中联立方程所形成的右端项，

为所有配置节点上的声压幅值；

进一步的，本发明中步骤S3中将步骤S2中所建立的线性方程组进行求解，获得所有配置节点上的声压幅值，然后利用如下公式预测各节点处的声压级：

SPL＝20log₁₀[p/p_ref]

其中p_ref为参考声压，取值2×10^-5Pa。

进一步的，本发明中步骤S4中车腔内的噪声衡量指标主要是声压级的分布，将步骤S3中获得的声压级分布绘制云图，观察最大声压级的分布情况，以及驾驶员和乘客耳朵附近的声压级大小，根据环境噪声标准对汽车噪声进行性能评估，声压级越小，则性能越优；声压级越大，则性能越差。

本发明与现有技术相比，取得的有益效果如下：本发明基于声学控制方程的基本解建立局部无网格基本解方法数值模型，将所建立的模型应用于汽车车内噪声的预测。所建立的无网格数值模型于此处最先提出，与传统的CAE软件核心算法有限元和边界元相比，具有无需网格划分，计算速度快，噪声预测精度高等优势。本发明由于利用了声学问题的半解析基本解，并且用任意分布的节点代替有限元网格，因此减轻了汽车CAE从业人员在前处理阶段的网格剖分负担、减小数值模型的离散误差并提高了汽车车内噪声预测的精度和效率。其主体构思巧妙，计算方式简便，能够简单、快捷、稳定、精确预测各类车腔噪声分布，为车内噪声预测提供了新的、简单高效的技术路线。应用环境友好，市场前景广阔。

附图说明：

图1为本发明涉及的基于局部无网格基本解法的车内噪声预测方法的流程图。

图2为本发明涉及的局部无网格基本解方法的示意图。

图3为本发明涉及的某轿车的二维声腔示意图。

图4为本发明涉及的局部基本解法的节点分布。

图5为本发明涉及的无网格局部基本解法与有限元法所得路径ab上的声压预测结果对比。

图6为本发明涉及的无网格局部基本解法与有限元法所得路径ab上的声压级预测结果对比。

图7为本发明涉及的无网格局部基本解法与有限元法所得驾驶员耳朵处的频率响应曲线对比。

具体实施方式：

下面结合附图和具体实施例，进一步阐明本发明，应理解此处实施例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。

实施例1：

如图1所示，本实施例涉及一种基于局部无网格基本解法的车内噪声预测方法，具体操作步骤按照如下方式进行：

S1、根据需预测车辆的几何参数，建立车腔的几何模型，并在车腔的边界及内部配置节点；

S2、利用车腔内空气的密度(ρ，kg/m³)、声音在空气中的传播速度(c，343m/s)、发动机的振动速度(v_n，m/s)、车顶吸声材料的声阻抗系数(A_n，m³/(Pa·s))，结合相应的边界条件，建立车内声学分析的局部无网格基本解法数值模型，其步骤包括：

(1)如图2所示，对任意节点x⁽ⁱ⁾，搜索其m个邻近节点

确定其所对应的局部子区域

并在局部子区域外一定距离d处配置相应数量的源点

(2)构建如下形式的局部矩阵H，并求其逆矩阵H^-1；

其中，G(r_kj)为声学控制方程的基本解，

为局部节点与源点之间的距离；

(3)利用给定边界条件建立相应的代数方程：对于内部节点，满足下列方程

本质边界节点，满足下列方程

p(x⁽ⁱ⁾)＝p_D

对于自然边界节点，满足下列方程

对于吸收边界节点，满足下列方程

上式中，

h⁽ⁱ⁾＝[G(r_i1),G(r_i2),...,G(r_im)]，

(4)联立(3)中的各方程，获得如下的最终稀疏矩阵方程；

Ap＝b

其中，A为(3)中联立方程所形成的系数矩阵，b为(3)中联立方程所形成的右端项，

为所有配置节点上的声压幅值。

S3、求解步骤S2中所建立的线性方程组，获得所有配置节点上的声压幅值，然后利用如下公式预测各节点处的声压级：

SPL＝20log₁₀[p/p_ref]

其中p_ref为参考声压，取值2×10^-5Pa。

S4、输出车腔内各节点处声压和声压级的预测值，绘制云图，观察最大声压级的分布情况，以及驾驶员和乘客耳朵附近的声压级大小。根据环境噪声标准对汽车噪声进行性能评估，声压级越小，则性能越优；声压级越大，则性能越差。

实施例2：

在本实施例中，基于局部无网格基本解法的车内噪声预测方法被用于某一轿车的车腔噪声预测。图3为某款小轿车二维车内声腔模型，其水平方向的长度L_x＝2.664m，竖直方向的高度L_y＝1.121m。由于发动机振动引起的噪声占车内噪声的绝大部分，因而这里假设汽车的前围板上承受一法向振动速度v_n＝0.01m/s。乘员舱的顶部为吸声材料，其吸声材料的声阻抗系数为A_n＝0.00144m³/(Pa·s)。图3也给出了驾驶员耳朵所在位置，以及需要重点预测的路径ab。

根据实施例1中一种基于局部无网格基本解法的车内噪声预测方法的步骤，在车腔内部及边界上分别布置2551和233个离散的节点，如图4所示。然后选取局部节点数目为m＝18，虚假边界距离为d＝5R_i。利用所预测车腔的已知边界条件，建立局部无网格基本解法数值模型，并对该车腔噪声分布进行预测。图5和图6分别给出了500Hz频率下本发明所提方法与有限元法所得路径ab上的声压与声压级预测结果的对比。图7给出了驾驶员耳朵处的频率响应曲线。这里参考解是有限元方法采用非常细密的单元计算得到的，其网格模型的节点数为21819，单元数为43636。有限元的预测结果是利用1392个节点和2784个单元得到的。由于参考解的节点数远远大于数值模型中的节点数，因而将其作为对比的依据是合理可信的。能够得出，有限元与参考解的误差较大，而局部无网格基本解法的预测结果与参考解十分吻合。这里局部无网格基本解法所采用的节点数目与有限元所采用的单元数目相同，但是前者的预测精度高。另外需要指出，本实施例涉及的方法无需繁琐的网格划分，仅需随机分布的离散节点，这大大降低了网各类方法的前处理成本，提高了车腔噪声预测的效率。

综上，实施例1车内噪声预测方法以局部无网格基本解法为基础，直接利用声学问题的半解析基本解作为核函数，基于局部算法理论建立无网格方法用于车内噪声预测。该发明的特点是数学理论简单，易于实施，仅需离散节点，无需网格划分和数值积分，有效避免了前处理阶段的网格剖分困难、减小了数值模型的离散误差并提高了汽车车内噪声预测的精度和效率。另，实施例1的方法作为一种局部算法，相较于传统的全局算法更易于复杂几何形状及其大规模问题的噪声预测分析。实施例1可用于分析不同车型不同工况下车内噪声的分布规律，能够缩短车辆测试调校周期，降低车辆NVH研发成本。实验对比表明，应用实施例1的方法进行车内噪声预测，在采用相同单元或节点数目条件下，预测精度远远高于传统的有限元方法。

Claims (5)

1.一种基于局部无网格基本解法的车内噪声预测方法，其特征在于，包括以下步骤：

S1、根据需预测车辆的几何参数，建立车腔的几何模型，并在车腔的边界及内部配置节点；

S2、利用车腔内空气的密度、声音在空气中的传播速度等已知参数，结合相应的边界条件，建立车内声学分析的局部无网格基本解法数值模型；

S3、利用步骤S2所建立的数值模型，计算步骤S1中所配置节点上的声压和声压级，对车内噪声进行预测；

S4、输出待预测汽车车腔内的噪声分布，并对预测结果进行性能评估。

2.根据权利要求1所述的一种基于局部无网格基本解法的车内噪声预测方法，其特征在于，步骤S1中车腔的边界及内部配置节点包括：

在车腔边界上配置n_b个节x⁽ⁱ⁾,i＝1,2,...,n_b；

在车腔内部配置n_i个节x⁽ⁱ⁾,i＝1,2,...,n_i。

3.根据权利要求1所述的一种基于局部无网格基本解法的车内噪声预测方法，其特征在于，步骤S2中所采用的声学控制方程为：

其中

为拉普拉斯算子，Ω表示需预测车内区域，x为空间坐标，p(x)为x处的声压幅值，ω为角频率，c为声速；

边界条件为下列三种中的一种或多种：

p＝p_D,在Γ_D上

在Γ_N上

在Γ_R上

其中

ρ为空气的密度，n为边界点上的单位外法向量，p_D为本质边界Γ_D上的声压，v_n为自然边界Γ_N上的法向速度，A_n为吸收边界Γ_R上的阻尼系数；

步骤S2中建立局部无网格基本解法数值模型步骤包括：

(1)对任意节点x⁽ⁱ⁾，搜索其m个邻近节点

确定其所对应的局部子区域

并在局部子区域外一定距离d处配置相应数量的源点

(2)构建如下形式的局部矩阵H，并求其逆矩阵H^-1；

其中，G(r_kj)为声学控制方程的基本解，

为局部节点与源点之间的距离；

(3)利用给定边界条件建立相应的代数方程：对于内部节点，满足下列方程

本质边界节点，满足下列方程

p(x⁽ⁱ⁾)＝p_D

对于自然边界节点，满足下列方程

对于吸收边界节点，满足下列方程

上式中，

h⁽ⁱ⁾＝[G(r_i1),G(r_i2),...,G(r_im)]，

联立步骤S3下(3)中的各方程，获得如下的最终稀疏矩阵方程；

Ap＝b

其中，A为(3)中联立方程所形成的系数矩阵，b为(3)中联立方程所形成的右端项，

为所有配置节点上的声压幅值。

4.根据权利要求1所述的一种基于局部无网格基本解法的车内噪声预测方法，其特征在于，步骤S3中将步骤S2中所建立的线性方程组进行求解，获得所有配置节点上的声压幅值，然后利用如下公式预测各节点处的声压级：

SPL＝20log₁₀[p/p_ref]

其中p_ref为参考声压，取值2×10^-5Pa。

5.根据权利要求1-4中任一权利要求所述的一种基于局部无网格基本解法的车内噪声预测方法，其特征在于，步骤4中车腔内的噪声衡量指标主要是声压级的分布，将步骤S3中获得的声压级分布绘制云图，观察最大声压级的分布情况，以及驾驶员和乘客耳朵附近的声压级大小，根据环境噪声标准对汽车噪声进行性能评估，声压级越小，则性能越优；声压级越大，则性能越差。

Images