一种基于稀疏网格配点理论的车内随机振动噪声预测的方法
【专利摘要】本发明公开了一种基于稀疏网格配点理论的车内随机振动噪声预测的方法。首先根据工程领域的实际需求建立车内振动噪声预测的有限元模型，确定目标空间位置和目标频率范围；其次，在以随机模型实现相关不确定性的定量化后，基于稀疏网格配点理论对随机参数抽样，利用车内振动噪声预测的有限元模型计算每个随机参数样本点处的响应值；最后，根据离散格式响应值计算获得车内振动噪声响应的混沌多项式代理模型的系数矩阵，并进一步依此系数矩阵计算获得车内振动噪声的均值频响分布和方差频响分布。本发明同时考虑了外部载荷与结构材料参数、空气介质特性参数对车内振动噪声的随机效应，为车内噪声优化与控制等降噪措施的制定提供依据。
【专利说明】
-种基于稀疏网格配点理论的车内随机振动噪声预测的方法
技术领域
[0001] 本发明设及装备噪声测量的技术领域，具体设及一种基于稀疏网格配点理论的车 内随机振动噪声预测的方法，适用于汽车车内噪声分析、优化与控制领域。
【背景技术】
[0002] 随着世界汽车工业发展及市场需求的扩大，持续保持汽车产品的强势竞争力并不 断扩大市场占有份额是所有汽车生产商孜孜追求的目标。同时，随着科学技术的日新月异 及人们消费观念的改变，汽车的概念正由"代步工具"向"生活空间"发展和延伸，消费者对 汽车综合性能的要求越来越高。在考虑节油性、动力性、安全性、操纵稳定性、平顺性等性能 指标的同时，汽车噪声与振动特性(NVH)因消费者在乘坐或试驾过程中的最直观且下意识 的感受而受到越来越多的关注，特别是车内振动噪声性能正逐步成为衡量乘用车辆及其档 次的重要指标之一，对于高档及豪华车市场尤其如此。车内振动噪声是影响驾乘人员舒适 性、听觉损害程度、语言清晰度W及辨识车外各种声音信号能力的重要因素，且直接影响汽 车的耐久性和可靠性。经验表明，汽车产品在生产制造阶段的振动噪声特性的改进需付出 极其高昂的代价，因此在设计阶段就应将车内振动噪声的优化与控制作为重要内容和关键 环节之一。然而，车内振动噪声优化与控制的首要前提是准确地预测车内噪声特性。
[0003] 从辩证法的角度看，确定性是相对的，不确定性是绝对的。在车内振动噪声预测过 程中存在多源不确定性，主要表现为:其一，驾乘室的空气介质声学环境随溫度等因素的变 化而不断变化，引起空气介质密度与声速的波动；其二，作用于汽车结构的外部激励波动； 其=，汽车结构几何尺寸、材料属性等参数的波动;其四，不同边界条件(如焊接、馴接等）的 简化假设而造成的振动噪声分析模型的不确定性。运些不确定性在车内振动噪声预测过程 中广泛存在，且不确定性禪合效应会导致噪声特性出现较大波动。因此，考虑不确定性对车 内振动噪声的影响是十分必要且有意义的。在不确定参数试验数据的样本容量足够大而可 实现W高精度拟合参数的概率密度函数条件下，本发明W正态随机模型定量化不确定参数 并基于稀疏网格配点理论建立车内振动噪声的混浊多项式代理模型，进一步计算车内振动 噪声的均值与方差，提供了一种基于稀疏网格配点理论的车内随机振动噪声预测的方法。
【发明内容】
[0004] 本发明要解决的技术问题是:考虑材料参数、声场介质特性参数、外部载荷参数的 不确定性对车内振动噪声的影响，弥补当前该领域仅考虑随机载荷效应的空白，提供一种 准确的预测车内振动噪声统计特性的方法。
[0005] 本发明采用的技术方案是:一种基于稀疏网格配点理论的车内随机振动噪声预测 的方法，其实现步骤是：
[0006] 第一步:确定车内振动噪声预测问题的目标空间位置及目标频率范围，从而确定 车内振动噪声预测的响应列向量L为目标空间位置处在目标频率范围内不同频率下振动噪 声A声级所组成的向量；
[0007] 第二步:利用第一步中响应列向量L关于材料属性、外部载荷、声场介质特性参数 的灵敏度分析数据，确定车内振动噪声预测过程中不确定参数向量X所包含的参数，并根据 不确定参数试验数据Wn维随机参数行向量X定量化，其中均值行向量为ii，方差行向量为 O2，对应标准差行向量为曰；
[0008] 第=步:基于稀疏网格配点理论生成n维标准正态随机参数向量的化个样本点，并 逐行W标准随机样本点矩阵Mg存储，矩阵Mg的每行对应标准正态随机参数向量的一个具体 实现，利用第二步中均值行向量y和标准差行向量0,将标准随机样本点矩阵Mc映射至随机 参数列向量X空间内W生成随机样本点，并逐行W随机样本点矩阵Mx存储；
[0009] 第四步:将第=步获得的标准随机样本点矩阵Mc逐行代入n维二阶混浊多项式函 数组成的基底行向量r lx(n+l)(n+2)/2中，计算基底行向量r lx(n+l)(n+2)/2在Mc每行处的取值，并 逐行存储于基底矩阵H中；
[0010] 第五步：取随机参数向量X的具体实现为第=步中随机样本点矩阵Mx的第k行，并 将其代入车内振动噪声预测的有限元模型中，计算得到第一步中响应列向量L在随机参数 行向量X的第k个样本点处的值Lk,遍历随机样本点矩阵Mx的所有行，获得响应矩阵；
[0011] 第六步:根据第四步获得的基底矩阵H和第五步获得的响应矩阵化，计算响应列向 量L关于混浊多项式基底行向量r lx(n+l)(n+2)/2的展开系数矩阵C，并进一步基于系数矩阵C 计算车内振动噪声响应列向量L的均值列向量化和方差列向量O自，最终获得车内任意特定 空间位置处在给定频率范围内振动噪声A声级的频响分布。
[0012] 其中，所述方法W随机定量化模型考虑了除随机载荷效应外的汽车结构材料属性 参数、介质特性参数的不确定性对车内振动噪声的影响。
[0013] 其中，所述方法基于稀疏网格配点理论对随机参数抽样W建立车内振动噪声响应 的混浊多项式代理模型，并进一步通过混浊多项式系数获得车内振动噪声的均值频响分布 与方差频响分布。
[0014] 本发明与现有技术相比的优点在于：
[0015] (1)本发明不仅考虑了随机载荷效应，且同时考虑了工程领域内不可避免的材料 属性参数与声场介质特性参数的不确定性对车内振动噪声特性的影响。
[0016] (2)本发明基于稀疏网格配点理论对随机参数抽样W建立混浊多项式代理模型， 与随机领域其他分析方法相比，W牺牲小幅计算效率的代价获得计算精度的显著提升。
【附图说明】
[0017] 图1为基于稀疏网格配点法的车内随机振动噪声预测的原理图；
[0018] 图2为基于稀疏网格配点法的车内随机振动噪声预测的流程图；
[0019] 图3为车内振动噪声预测的空气介质有限元模型图；
[0020] 图4为车内振动噪声预测的结构有限元模型图；
[0021] 图5为车内振动噪声声压级的均值与方差的频响分布图。
【具体实施方式】
[0022] 下面结合附图W及【具体实施方式】进一步说明本发明。
[0023] 本发明一种基于稀疏网格配点理论的车内随机振动噪声预测的方法，首先建立如 图3和图4所示的车内振动噪声预测有限元模型并确定所关屯、空间位置及频率范围；其次， 在W随机模型对不确定性定量化后，基于稀疏网格配点理论对随机参数进行抽样并计算随 机参数样本点处响应值;最后基于响应的混浊多项式展开系数计算获得车内振动噪声的均 值频响分布与方差频响分布。但汽车车内声学品质有多种不同的度量与评价指标，本发明 选择但不限于车内振动噪声的A声级，其它度量与评价指标可W平行应用至本发明中。如图 2所示，其具体实施步骤是：
[0024] 第一步：如图2所示，确定车内振动噪声预测问题的目标空间位置(如人耳所处空 间位置)及目标频率范围（如75赫兹至100赫兹之间），从而确定车内振动噪声预测的响应列 向量L为目标空间位置处目标频率范围内不同频率下振动噪声A声级所组成的向量；
[0025] 第二步:利用第一步中响应列向量L关于材料属性、外部载荷、声场介质特性参数 的灵敏度分析数据，确定车内振动噪声预测过程中不确定参数向量X所包含的参数，并根据 不确定参数试验数据Wn维随机参数行向量X定量化，其中均值行向量为ii，方差行向量为 O2，对应标准差行向量为曰；
[0026] 第=步:根据随机参数的分布类型选择配点集合，本发明选择但不限于正态随机 参数为对象，WGauss-Hermite点（如表1)为配点集合，利用稀疏网格配点配点理论生成n维 标准正态随机参数向量的化个样本点，并逐行W标准随机样本点矩阵Mc存储，矩阵Mc的每行 对应标准正态随机参数向量的一个具体实现。利用第二步中均值行向量y和标准差行向量 O，将标准随机样本点矩阵Mc映射至随机参数列向量X空间内W生成随机样本点，并逐行W 随机样本点矩阵Mx存储，有：
[0027] Mx 化，：）= 〇〇Mc 化，：）+y，k=l，2，...，^ (1)
[0028] 其中（k，：）表示对应矩阵的第k行，k为矩阵行索引值，符号V'表示两个行向量的 对应元素相乘且结果仍为相同维数的行向量；
[0029] 表 1
[0030]
[0031] 第四步：在n维随机参数空间内，二阶混浊多项式函数组成的基底行向量 厂 lX(n+l)(n+2)/2为：
[0032] (2)
[0033]
[0034] 巧
[003引而二阶基底是由矩阵r海行中非零元素逐行顺序组合而得，矩阵r康示为：
[0036]
(4)
[0037] 其中. . .，Cn表示相互独立的标准正态随机变量。将第=步中标准随机样本 点矩阵Mc逐行代入基底行向量r ix(n+l)(n刮/2中W计算基底行向量r ix(n+l)(n刮/2在Mc每行处 的取值.并该斤在储干甚底巧陈H中.有，
[00；3引
（苟
[0039] 第立步：取随机参数问量X的具体实现为第S步中随机祥本点矩阵Mx的第k行，并 将其代入如图3和图4所示的车内振动噪声预测的有限元模型中，计算得到第一步中响应列 向量L在随机参数行向量X的第k个样本点处的值Lk,遍历随机样本点矩阵Mx的所有行，获得 响应矩阵化，有：
[0040]
(6)
[0041] 第六步:根据第四步中基底矩阵H和第五步中响应矩阵化，计算响应列向量L关于 混浊多项式基底行向量r lx(n+l)(n+2)/2的展开系数矩阵C，有：
[0042] C=化)-i(Ml)t (7)
[0043] 其中（?尸表示矩阵的广义逆，（?）T表示矩阵的转置矩阵。进一步，车内振动噪声 响应列向量L的均值列向量化和方差列向量G的作化十算为：
[0044] 化(l)=c(l，l) (8)
[0045]
(9)
[0046] 其中1是响应列向量L中分量的索引值，N表示响应列向量L的维数，即响应个数，分 块行向量百表示混浊多项式基底行向量rix(n+l)(n+2)/2的方差，有：
[0047] (10)
[004引 与：
[0049] 抑
[0050] 且行向量6:由矩阵化每行中非零元素逐行顺序组合而得，矩阵化表示为：
[0化1]
(12)
[0052]根据式(8)和(9)最终获得车内任意特定空间位置处在给定频率范围内振动噪声A 声级的频响分布。
[0053]本发明W图3和图4所示的车内振动噪声预测的有限元模型为对象，根据第一步至 第六步所述的一种基于系数网格配点理论的车内随机振动噪声预测方法的【具体实施方式】， 确定随机参数及其定量化结果如表2所示，利用基于Gaus S-化rmi te点的稀疏网格配点理论 生成的随机参数向量样本点如表3所示，获得驾驶者左耳位置处在75赫兹至100赫兹频率范 围内A声级的均值与方差的频响分布如图5所示。
[0化4] 表2
[0化6] 表3
[0化9]
[0060]本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
【主权项】
1. 一种基于稀疏网格配点理论的车内随机振动噪声预测的方法，其特征在于包括以下 步骤： 第一步：确定车内振动噪声预测问题的目标空间位置及目标频率范围，从而确定车内 振动噪声预测的响应列向量L为目标空间位置处在目标频率范围内不同频率下振动噪声A 声级所组成的向量； 第二步:利用第一步中响应列向量L关于材料属性、外部载荷、声场介质特性参数的灵 敏度分析数据，确定车内振动噪声预测过程中不确定参数向量X所包含的参数，并根据不确 定参数试验数据以η维随机参数行向量X定量化，其中均值行向量为μ，方差行向量为σ 2,对 应标准差行向量为σ; 第三步:基于稀疏网格配点配点理论生成η维标准正态随机参数向量的Ns个样本点，并 逐行以标准随机样本点矩阵Me存储，矩阵Me的每行对应标准正态随机参数向量的一个具体 实现，利用第二步中均值行向量μ和标准差行向量〇,将标准随机样本点矩阵Me映射至随机 参数列向量X空间内以生成随机样本点，并逐行以随机样本点矩阵M x存储； 第四步:将第三步获得的标准随机样本点矩阵Me逐行代入η维二阶混沌多项式函数组成 的基底行向量Γ 1Χ(η+1)(η+2)/2中，计算基底行向量r 1Χ(η+1)(η+2)/2在Me每行处的取值，并逐行存 储于基底矩阵Η中； 第五步:取随机参数向量X的具体实现为第三步中随机样本点矩阵Μχ的第k行，并将其代 入车内振动噪声预测的有限元模型中，计算得到第一步中响应列向量L在随机参数行向量X 的第k个样本点处的值Lk，遍历随机样本点矩阵M x的所有行，获得响应矩阵； 第六步:根据第四步获得的基底矩阵Η和第五步获得的响应矩阵Ml，计算响应列向量L关 于混沌多项式基底行向量Γ lx(n+lKn+2)/2的展开系数矩阵C，并进一步基于系数矩阵C计算车 内振动噪声响应列向量L的均值列向量R和方差列向量<，最终获得车内任意特定空间位 置处在给定频率范围内振动噪声Α声级的频响分布。2. 根据权利要求1所述的基于稀疏网格配点理论的车内随机振动噪声预测的方法，其 特征在于，所述方法以随机定量化模型考虑了除随机载荷效应外的汽车结构材料属性参 数、介质特性参数的不确定性对车内振动噪声的影响。3. 根据权利要求1所述的基于稀疏网格配点理论的车内随机振动噪声预测的方法，其 特征在于，所述方法基于稀疏网格配点理论对随机参数抽样以建立车内振动噪声响应的混 沌多项式代理模型，并进一步通过混沌多项式系数获得车内振动噪声的均值频响分布与方 差频响分布。
【文档编号】G06F17/50GK105956283SQ201610295005
【公开日】2016年9月21日
【申请日】2016年5月6日
【发明人】邱志平, 许孟辉, 王晓军, 王冲, 李云龙, 郑宇宁, 陈贤佳, 张泽晟
【申请人】北京航空航天大学