一种基于代理模型的车内振动噪声预测的顶点分析方法
【专利摘要】本发明一种基于代理模型的车内振动噪声预测的顶点分析方法，首先根据工程领域的实际需求建立车内振动噪声预测的有限元模型，通过给定目标空间位置和目标频率范围而确定目标响应；其次，以区间模型定量化不确定参数，利用高斯积分点对区间参数抽样，通过车内振动噪声预测的有限元模型计算区间参数样本点处的真实响应；最后，基于高维模型表征理论以离散格式真实响应值建立目标响应的代理模型，并基于顶点分析理论计算其在超正方体所有顶点处取值，最终获得车内振动噪声区间界限的频响分布。本发明克服了车内振动噪声预测对参数大样本容量试验数据的依赖与小区间参数的适用性限制，为车内噪声优化与控制等降噪措施的制定提供依据。
【专利说明】
一种基于代理模型的车内振动噪声预测的顶点分析方法
技术领域
[0001] 本发明涉及装备噪声测量的技术领域，具体涉及一种基于代理模型的车内振动噪 声预测的顶点分析方法，适用于汽车车内噪声分析、优化与控制领域。
【背景技术】
[0002] 随着科学技术的进步与经济水平的提高，世界汽车市场的竞争日趋激烈，在每个 市场层次上都充斥着各大汽车生产制造商不断推出的在功能上大同小异的类似车型。随着 道路状况的改善及车辆振动平顺性研究成果的应用，车内噪声在评价汽车乘坐舒适性中起 着越来越重要的作用。在这种背景下，作为汽车品质最直观体验的车内声学环境已然成为 影响消费者购买意愿的重要决策因素。同时，作为汽车零部件机械疲劳重要成因之一的车 内噪声也直接影响汽车的耐久性和可靠性。以改善车内声学品质、降低车内噪声水平为目 标，车内噪声优化与控制作为汽车设计制造领域的重要课题之一而受到研究人员的日益重 视。振动噪声作为车内噪声的主要来源，其准确预测是进行车内噪声优化与控制的重要前 提。
[0003] 然而，多源不确定性广泛存在于车内振动噪声预测过程中，如车身壁板厚度等几 何参数与结构材料属性等材料参数因生产制造工艺限制而波动、空气介质特性因环境温度 变化而变化、外部载荷因路面平整状况变化或发动机转速变化而不同、噪声分析模型的边 界条件等简化假设而引入误差，这些因素均诱导了车内振动噪声的不确定性。当前不确定 性声学分析领域主要通过随机有限元方法考虑了随机载荷对声学响应特性的影响规律，相 应研究成果可应用于车内随机振动噪声预测。然而，随机有限元方法是基于不确定参数的 大样本容量试验数据拟合而得的概率密度函数或分布函数。由于试验条件或经济代价等客 观条件的限制，工程领域内获得大容量试验数据对某些参数难以实现，在不确定参数试验 数据有限而难以高精度拟合其概率密度函数或分布函数条件下，区间模型定量化方法因可 以克服对大样本容量试验数据的限制而成为车内振动噪声预测的可行解决途径。在这种背 景下，该领域当前主要提出了区间摄动分析方法，但潜在保守估计或区间平移等限制。因 此，本发明基于代理模型发明了一种车内振动噪声预测的顶点分析方法。
【发明内容】
[0004] 本发明要解决的技术问题是:克服车内随机振动噪声预测方法对参数的大样本容 量试验数据的依赖，克服已有方法仅适用于小区间参数的适用性限制，提供一种具有满意 计算精度与计算效率的车内振动噪声波动范围的预测方法。
[0005] 本发明采用的技术方案是:一种基于代理模型的车内振动噪声预测的顶点分析方 法，其实现步骤是：
[0006] 第一步:建立车内振动噪声预测的有限元模型，确定目标空间位置及目标频率范 围，进一步确定车内振动噪声目标响应列向量P为目标空间位置处目标频率范围内不同频 率下振动噪声声压级组成的向量，维数为Np;
[0007] 第二步:根据第一步中响应列向量P关于材料属性、空气介质属性、外部载荷I参 数的灵敏度分析数据，给定不确定参数向量z所包含的具体参数，并基于参数的试验数据以 区间模型将其定量化为η维区间参数列向量z 1，其下界列向量与上界列向量分别表示为/和 ζυ，中点列向量和半径列向量分别表示为F和
[0008] 第三步:根据第一步中响应列向量P关于区间参数列向量Z1每个分量的非线性程 度，确定响应列向量P在区间参数空间内每个维度上的最佳平方逼近函数的阶数N，并依此 确定高斯积分点的个数m，进一步确定高斯积分点列向量为XC=[X1，X2, . . .，Xm]T，利用第二 步获得的中点列向量Zc"和半径列向量Zlr对区间参数列向量Z抽样，将样本点存储于分块形式 的输入样本点矩阵Minput中；
[0009] 第四步:将第三步中输入样本点矩阵Minput逐行输入车内振动噪声有限元模型中， 计算第一步中响应列向量P在每个样本点处的取值，并以第三步中输入样本点矩阵M input的 分块格式将响应列向量P的值存储于输出样本点矩阵Mciutput中；
[0010] 第五步:根据第四步中输出样本点矩阵Mciutput，以Chebyshev多项式为基底函数，建 立响应列向量P的任意第1个分量关于区间参数向量Z 1的任意第i个分量的最佳平方逼近函 数户；)，进一步将最佳平方逼近函数P(//) h ；)顺序组成最佳平方逼近函数矩阵PMat(x);
[0011] 第六步:将区间参数列向量Z1的具体实现取为中点值列向量z%利用车内振动噪声 预测的有限元模型计算获得响应列向量名义值P%根据高维模型表征理论与第五步中最佳 平方逼近函数矩阵P Mat(x)，建立响应列向量P的代理模型
[0012] 第七步：计算由η个标准区间[_1，1]张成的η维超正方体所有顶点的空间位置坐 标，并逐列存储于顶点矩阵Mv中；
[0013] 第八步:将第七步中顶点矩阵Mv逐列代入第六步中响应列向量代理模型Is(X)中， 计算获得响应矩阵，利用顶点分析方法计算获得响应列向量P的区间界限Ρ1，即车内振 动噪声区间界限的频响分布。
[0014] 其中，所述方法以区间模型实现车内振动噪声不确定参数的定量化，并利用高斯 积分点生成区间参数样本点。
[0015] 其中，所述方法利用基于Chebyshev多项式的最佳平方逼近理论与高维模型表征 理论建立车内振动噪声预测的目标响应向量的代理模型。
[0016] 其中，所述方法利用顶点分析方法计算车内振动噪声响应的区间界限。
[0017] 本发明与现有技术相比的优点在于：
[0018] (1)本发明以区间模型定量化不确定参数而克服了车内振动噪声预测对参数大样 本容量试验数据的依赖，实现在有限试验数据条件下车内振动噪声波动范围的预测；
[0019] (2)本发明克服了区间摄动分析方法仅适用于小区间参数车内振动噪声预测的适 用性限制，并进一步避免了该方法在大区间参数输入条件下车内振动噪声预测的可观计算 代价。
【附图说明】
[0020] 图1为基于代理模型的车内振动噪声顶点分析方法的原理图；
[0021] 图2为基于代理模型的车内振动噪声顶点分析方法的流程图；
[0022] 图3为车内振动噪声预测的空气介质有限元模型图；
[0023] 图4为车内振动噪声预测的汽车结构有限元模型图；
[0024] 图5为车内振动噪声声压级区间界限的频响分布图。
【具体实施方式】
[0025]下面结合附图以及【具体实施方式】进一步说明本发明。
[0026]本发明一种基于代理模型的车内振动噪声预测的顶点分析方法，首先根据工程领 域的实际需求建立车内振动噪声预测的有限元模型，通过给定目标空间位置和目标频率范 围而确定目标响应;其次，以区间模型定量化不确定参数，利用高斯积分点对区间参数抽 样，通过车内振动噪声预测的有限元模型计算区间参数样本点处的真实响应;最后，基于高 维模型表征理论以离散格式真实响应值建立目标响应的代理模型，并基于顶点分析理论计 算其在超正方体所有顶点处取值，最终获得车内振动噪声区间界限的频响分布。如图2所 示，其具体实施步骤是：
[0027]第一步:建立如图3和图4所示的车内振动噪声预测的有限元模型，确定目标空间 位置(如人耳所处空间位置)及目标频率范围（如75赫兹至100赫兹之间），进一步确定车内 振动噪声目标响应列向量P为目标空间位置处目标频率范围内不同频率下振动噪声声压级 组成的向量，维数为Np;
[0028] 第二步:根据第一步中响应列向量P关于材料属性、空气介质属性、外部载荷等参 数的灵敏度分析数据，给定不确定参数向量z所包含的具体参数(如材料弹性模量与外部载 荷），并基于参数的试验数据以区间模型将其定量化为η维区间参数列向量z 1，其下界列向 量与上界列向量分别表示为ζ4Ρζυ，中点列向量Zcl和半径列向量Z lr分别计算为：
[0029]
(1)
[0030] 第三步:根据第一步中响应列向量P关于区间参数列向量Z1每个分量的非线性程 度，确定响应列向量P在区间参数空间内每个维度上的最佳平方逼近函数的阶数N，该最佳 平方逼近函数以Chebyshev正交多项式为基底，依此确定高斯积分点的个数m，满足关系：
[0031] m 彡 2N+1 (2)
[0032] 进一步确定高斯积分点列向量为XC= [X1，X2, ...，Xm]T，利用第二步获得的中点列 向量Zc=和半径列向量Zlr对区间参数列向量z抽样，将样本点存储于分块形式的输入样本点矩 阵Minput中，有：
[0033]
(3)
[0034] 其中任意第i个分块输入样本点矩阵的确定过程为:利用高斯积分点列向量XG 对区间参数列向量z的第i个分量抽样，获得m个样本点为..,Zw,将此m个样本点逐 一与其它n-1个区间参数的中点值2丨（./V /j:组合形成行向量，组合过程以各参数在区间参数 列向量Z1的顺序为标准，并逐行排列而确定取，即有：
[0035]
(4)
[0036] 其中（：，j)表示对应矩阵的第j列元素，（j)表示对应向量的第j个元素，3^为 Kronecker符号，满足：
[0037]
(5)
[0038]其中i，j表不参数在区间参数列向量Z1中的索引值，二者相等时Sij为1，否则为〇。 [0039]第四步:将第三步中输入样本点矩阵Minput逐行输入图3和图4所示的车内振动噪声 有限元模型中，计算第一步中响应列向量P在每个样本点处的取值，并以第三步中输入样本 点矩阵M input的分块格式将响应列向量P的值存储于输出样本点矩阵Mciutput中，有：
[0040]
(6)
[0041] 其中任意第i个分块输出样本点矩阵兪;对应于第三步中第i个分块输入样本点矩 阵1?，且M;表示为：
[0042]
(7)
[0043] 其中Pf1表示对应于区间参数列向量Z1第i个分量的第j个输入样本点处的响应列 向量P的取值，即将图3和图4所示的车内振动噪声预测有限元模型的输入参数取为输入样 本点矩阵^^的第j行而计算获得的响应列向量P。
[0044] 第五步:根据第四步中输出样本点矩阵Mcmtput，以Chebyshev多项式为基底函数，建 立响应列向量P的任意第1个分量关于区间参数向量Z 1的任意第i个分量的最佳平方逼近函 数戶"'l.v,),有：
[0045]
[0046]
[0047]
[0048] 其中xG(k)表示第三步中高斯积分点列向量XG的第k个分量，Ρ^/；)表示第四步中 响应列向量的第1个分量。将最佳平方逼近函数顺序组成如下格式的最佳平方 逼近函数矩阵PMat(X)，有：
[0049]
[0050]其中X=[X1，X2,…，Xn]是由位于标准区间[-1,1]内的变量X1，X2,…，Xn组成的向 量。
[0051]第六步:将区间参数列向量Z1的具体实现取为中点值列向量z%利用车内振动噪声 预测的有限元模型计算获得响应列向量名义值Pe。根据高维模型表征理论与第五步中最佳 平方逼近函数矩阵PMat(x)，建立响应列向量P的代理模型尹⑷，有：
[0052]
(Ii)
[0053] 其中1。。1表示Np维元素均为1的列向量，有：
[0054] Ι〇〇ι = [1,1,...,1]τ (12)
[0055] 第七步：计算由η个标准区间[_1，1]张成的η维超正方体所有顶点的空间位置坐 标，并逐列存储于顶点矩阵Mv中。以如下格式遍历i和j，即：
[0059] 第八步:将第七步中顶点矩阵Mv逐列代入第六步中响应列向量代理模型中， 计算获得响应矩阵Mjs，有：
[0056]
[0057]
[0058] (14)：
[0060]
(15)
[0061] 利用顶点分析方法和响应矩阵计算获得响应列向量P的区间界限Pi，其下界 列向量炒和上界列向量Pu分别为：
[0062] (I6)
[0063] 分别表不对2n个向量的对应位置分 量取最小值和最大值而获得的列向量。最终获得的响应列向量区间界限P1即为车内振动噪 声区间界限的频响分布。
[0064] 本发明以图3和图4所示的车内振动噪声预测的有限元模型为对象，根据第一步至 第八步所述的基于代理模型的车内振动噪声预测的顶点分析方法的【具体实施方式】，确定区 间参数及其定量化结果如表1所示，计算获得驾驶者左耳位置处在75赫兹至100赫兹频率范 围内车内振动噪声的声压级区间界限的频响分布如图5所示。
[0065] 表 1
LOO67」本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
【主权项】
1. 一种基于代理模型的车内振动噪声预测的顶点分析方法，其特征在于包括以下步 骤： 第一步:建立车内振动噪声预测的有限元模型，确定目标空间位置及目标频率范围，进 一步确定车内振动噪声目标响应列向量P为目标空间位置处目标频率范围内不同频率下振 动噪声声压级组成的向量，维数为NP; 第二步:根据第一步中响应列向量P关于材料属性、空气介质属性、外部载荷等参数的 灵敏度分析数据，给定不确定参数向量z所包含的具体参数，并基于参数的试验数据以区间 模型将其定量化为n维区间参数列向量z 1，其下界列向量与上界列向量分别表示为zl和z11， 中点列向量和半径列向量分别表示为一和^; 第三步:根据第一步中响应列向量P关于区间参数列向量z1每个分量的非线性程度，确 定响应列向量P在区间参数空间内每个维度上的最佳平方逼近函数的阶数N，并依此确定高 斯积分点的个数m，进一步确定高斯积分点列向量为XC=[X1，X2, . . .，Xm]T，利用第二步获得 的中点列向量z""和半径列向量Z1对区间参数列向量Z抽样，将样本点存储于分块形式的输入 样本点矩阵Minput中； 第四步:将第三步中输入样本点矩阵Mlnput逐行输入车内振动噪声有限元模型中，计算 第一步中响应列向量P在每个样本点处的取值，并以第三步中输入样本点矩阵Minput的分块 格式将响应列向量P的值存储于输出样本点矩阵M?t_中； 第五步:根据第四步中输出样本点矩阵Mcmtput，以Chebyshev多项式为基底函数，建立响 应列向量P的任意第1个分量关于区间参数向量z1的任意第i个分量的最佳平方逼近函数 户"，1 (x,)，进一步将最佳平方逼近函数戶叫(x,)顺序组成最佳平方逼近函数矩阵PMat(x); 第六步:将区间参数列向量z1的具体实现取为中点值列向量匕利用车内振动噪声预测 的有限元模型计算获得响应列向量名义值P%根据高维模型表征理论与第五步中最佳平方 逼近函数矩阵PMat(x)，建立响应列向量P的代理模型P(x): 第七步:计算由n个标准区间[_1，1]张成的n维超正方体所有顶点的空间位置坐标，并 逐列存储于顶点矩阵Mv中； 第八步:将第七步中顶点矩阵Mv逐列代入第六步中响应列向量代理模型中，计算 获得响应矩阵Mf，利用顶点分析方法计算获得响应列向量P的区间界限P1，即车内振动噪 声区间界限的频响分布。2. 根据权利要求1所述的基于代理模型的车内振动噪声预测的顶点分析方法，其特征 在于，所述方法以区间模型实现车内振动噪声不确定参数的定量化，并利用高斯积分点生 成区间参数样本点。3. 根据权利要求1所述的基于代理模型的车内振动噪声预测的顶点分析方法，其特征 在于，所述方法利用基于Chebyshev多项式的最佳平方逼近理论与高维模型表征理论建立 车内振动噪声预测的目标响应向量的代理模型。4. 根据权利要求1所述的基于代理模型的车内振动噪声预测的顶点分析方法，其特征 在于，所述方法利用顶点分析方法计算车内振动噪声响应的区间界限。
【文档编号】G06F19/00GK106055861SQ201610291147
【公开日】2016年10月26日
【申请日】2016年5月5日
【发明人】邱志平, 许孟辉, 王晓军, 王冲, 李云龙, 郑宇宁, 王磊, 陈贤佳
【申请人】北京航空航天大学