一种基于高维模型表征与顶点分析策略的空腔噪声预测方法
【专利摘要】本发明公开了一种基于高维模型表征与顶点分析策略的空腔噪声预测方法。首先，将空腔噪声预测过程中的不确定性以相互独立的区间模型定量化，利用空腔噪声响应关于参数非线性程度确定其最佳平方逼近函数的阶数及高斯积分点。其次，利用高斯积分点对区间参数进行抽样，计算每个区间参数样本点处空腔噪声响应的离散值，进一步利用高维模型表征理论建立基于逐维最佳平方逼近函数的代理模型。最后，计算由所有区间参数张成的超立方体顶点的空间位置坐标及相应顶点处代理模型的取值，进一步确定区间参数影响下空腔噪声响应的区间波动范围。本发明考虑了不确定性对空腔噪声预测的影响，与实际工程领域具有良好吻合性。
【专利说明】
-种基于高维模型表征与顶点分析策略的空腔噪声预测方法
技术领域
[0001] 本发明设及装备噪声测量的技术领域，具体设及一种基于高维模型表征与顶点分 析策略的空腔噪声预测的方法，适用于汽车、飞机、轮船、潜艇等具有空腔结构的内部噪声 预测。
【背景技术】
[0002] 由弹性壳体围合而成的空腔结构广泛存在于诸如汽车、飞机、轮船、潜艇、航天器 等海陆空天各类军用与民用运输装备中，弹性壳体在外界激励作用下振动而向空腔内福射 噪声，即结构振动声福射。科学技术的日新月异让人们追求高品质生活环境、高舒适性交通 工具与高安全性武器装备的梦想照进了现实，并因此强烈凸显了结构振动噪声问题的重要 性。过高水平的噪声不仅严重影响相关装备内部乘员的身屯、健康与舒适性体验，而且极大 程度地妨碍甚至危害装备综合性能或内部关键设备的正常运行状态。因此，各类运输装备 都面临越来越多的如何减少甚至避免弹性壳体振动噪声危害的问题，而空腔噪声的准确快 速预测是相关装备声学品质优化改善的首要前提。
[0003] 然而，不确定性广泛存在于空腔噪声预测过程中，如:海水环境变化、路面平整状 况差异、阵风载荷及相关装备发动机不同转速等引起的外部激励载荷的不确定性;环境变 化(如溫度）引起的声场流体介质特性(如声速与质量密度)的波动;生产工艺限制等引起的 系统几何参数加工误差与结构材料属性差异;研究人员主观认知水平限制、制定相关规范 与执行标准等诱导的认知不确定性。研究表明：小幅不确定性在结构与声场禪合特性影响 下导致声学特性大幅偏差，甚至出现反相现象。因此，研究不确定性对空腔噪声的影响规律 对空腔噪声预测具有显著工程价值。相关不确定参数的定量化是进行空腔噪声预测的必要 前提，基于概率与统计理论的定量化方法需要不确定参数的大量统计数据，而由于试验条 件等限制导致统计数据所拟合的参数概率密度函数等精度有限。在工程领域内，当不确定 参数试验数据的样本容量十分有限时，可行途径是利用区间模型定量化不确定参数，即区 间参数。在运种背景下，当前空腔噪声分析主要W区间摄动分析方法完成。针对区间摄动分 析方法的限制，本发明基于高维模型表征与顶点分析策略发明了一种空腔噪声预测的方 法。
【发明内容】
[0004] 本发明要解决的技术问题是：克服当前方法仅可应用于小区间参数的适用性限 审IJ，克服当前方法在大区间参数输入条件下的计算效率限制，提供一种快速准确预测空腔 噪声波动范围的方法。
[0005] 本发明采用的技术方案是一种基于高维模型表征与顶点分析策略的空腔噪声预 测的方法，其实现步骤是：
[0006] 第一步:根据工程需求确定空腔噪声分析的响应列向量P所包含的响应分量(如空 腔内特定位置处不同频率点处的声压级等），根据研究人员或者工程经验数据确定不确定 参数列向量Z所包含的不确定参数(如结构材料弹性模量、结构材料质量密度、结构几何参 数、空腔内流体介质参数等），基于不确定参数列向量Z的试验数据，利用区间模型将其定量 化为区间参数列向量zi，其下界列向量为/，上界列向量为Z"，中点列向量为Z可P半径列向 量为zT。
[0007]第二步:根据第一步中响应列向量P关于区间参数列向量zi的非线性程度，确定响 应列向量P关于每个区间参数最佳平方逼近函数的阶数N、高斯积分点个数m及高斯积分点 列向量
。利用第一步确定的区间参数列向量zl的中点列向量2；\半径列向 量Zt及高斯积分点列向量S，对区间参数列向量Z进行抽样，将样本点存储于输入样本点矩 阵化中。
[000引第=步:将第二步中所得到的输入样本点矩阵Mi逐行代入空腔噪声分析的有限元 模型中，计算响应列向量P在每个样本点处的值，并存储于输出样本点矩阵Mo。
[0009] 第四步:根据第=步所获得输出样本点矩阵Mo建立响应列向量P的第1个分量关于 第i个区间参数的最佳平方逼近函数aU'i>(x)。
[0010] 第五步:将第一步中区间参数向量zi的中点列向量代入空腔噪声分析的有限元 模型中，计算响应列向量值名义值为pe。基于高维模型表征理论，根据第四步所获得最佳平 方逼近函数aU'i>(x)，建立响应列向量P的第1个分量在n维区间参数空间内的代理模型 (X)。
[0011] 第六步:根据顶点分析策略确定n维超正方体所有顶点的空间位置坐标，并存储于 顶点矩阵Mv中。
[0012] 第屯步:基于第五步所获得的第1个响应分量的代理模型AW(X)与第六步所获得 的顶点矩阵Mv,将顶点矩阵Mv的第k列代入代理模型A^(X)中，计算获得n维超正方体第k个 顶点处响应列向量P第1个分量的离散值，进一步根据顶点分析策略确定响应列向量P第 1个分量的区间界限巧。
[0013] 第八步:利用第四步至第六步遍历响应分量索引值1，可W计算得到响应列向量P 的下界列向量PL和上界列向量PU，最终可W确定空腔噪声响应的波动范围为pi = [pL，pU]。
[0014] 本发明与现有技术相比的优点在于：
[0015] (1)本发明考虑了空腔噪声分析过程中不确定性对空腔噪声预测的影响，与工程 实际领域具有良好吻合性。
[0016] (2)本发明克服了区间摄动分析方法仅适用于小区间参数空腔噪声分析问题的限 审IJ，具有更广泛的适用性。
[0017] (3)本发明对大区间参数空腔噪声分析问题具有满意的计算效率，有效避免了区 间摄动分析方法在该类分析问题中可观的计算代价。
【附图说明】
[0018] 图1为基于高维模型表征与顶点分析策略的空腔噪声预测方法的原理图；
[0019] 图2为基于高维模型表征与顶点分析策略的空腔噪声预测方法的流程图；
[0020] 图3为圆柱壳体的几何构型图；
[0021] 图4为频率为180赫兹的圆柱形空腔内中线Z向各位置处声压级区间界限；
[0022] 图5为频率为196赫兹的圆柱形空腔内中线Z向各位置处声压级区间界限；
[0023] 图6为频率为220赫兹的圆柱形空腔内中线Z向各位置处声压级区间界限。
【具体实施方式】
[0024] 下面结合附图W及【具体实施方式】进一步说明本发明。
[0025] 本发明一种准确的空腔噪声预测的方法，首先基于高维模型表征理论建立空腔噪 声响应关于区间参数的代理模型，进一步利用顶点分析策略确定区间参数影响下空腔噪声 响应的波动范围。如图2所示，其具体实施步骤是：
[0026] 第一步：如图2所示，根据工程需求确定空腔噪声分析的响应列向量P所包含的响 应分量(如空腔内特定位置处不同频率点处的声压级等），根据研究人员经验确定空腔噪声 分析的不确定参数列向量Z所包含的不确定参数（如结构材料弹性模量、结构材料质量密 度、结构几何参数、空腔内流体介质参数等），基于不确定参数列向量Z的试验数据，利用区 间模型将其定量化为区间参数列向量zi，其下界列向量为/，上界列向量为Z"，中点列向量 为Z。和半巧列向畳为zT,滿巧羊系，
[0027]
[002引 巧
[0029] 第二步:如图2所示，根据第一步确定的响应列向量P和区间参数列向量zi，基于研 究人员经验与参数灵敏度分析数据等确定响应列向量P关于区间参数列向量zi的非线性程 度，从而确定响应列向量P关于每个区间参数最佳平方逼近函数的阶数N、高斯积分点个数m 及高斯积分点列向量i:。.利用第一步确定的中点列向量(如式（1 ))、半径列向量(如式 (2))及高斯积分点列向量& =巧，起，...，对区间参数列向量Z进行抽样，将样本点存储于 输入样本点矩阵化中，输入样本点矩阵化W分块矩阵形式表示为：
[0030] Mi=[S(i);S(2); . . . ;S(n)] (3)
[0031] 其中关于第i个区间参数的输入样本点矩阵SW的维数为mXn，其具体计算为：
[0032]
{]}
[0033] 其中Sy表示Kronecker符号，满足：
[0034；
(5)
[0035] 其中i，j表示参数在区间参数列向量zi中的索引值，二者相等时Kronecker符号取 值为1，否则取值为0。
[0036] 第S步：如图2所示，将第二步中式(3)所表示的输入样本点矩阵Mi逐行代入空腔 噪声分析的有限元模型中，计算响应列向量P在每个区间参数列向量的样本点处的值，并存 储于输出样本点矩阵Mo,输出样本点矩阵MoW分块矩阵形式表示为：
[0037]
銜 [003引其中关于第i个区间参数的输出样本点矩阵Mf*表示为：
[0039]
(I)
[0040] 输出样本点矩阵的第j列if表示对应第i个区间参数的第j个输入样本点处的 响应值，即堪输入样本点矩阵第j行处计算的响应列向量P的值。
[0041] 第四步:如图2所示，根据第S步所获得输出样本点矩阵Mo建立响应列向量P的第1 个分量关于第i个区间参数的最佳平方逼近函数aU'i>(x)，有：
[0042]
(8)
[0043] 其中Lk(X)表示第k阶勒让德多项式，最佳平方逼近函数的系数可W计算为：
[0044]
蹲
[0045] 其中考表示第二步中高斯积分点S的第j个分量，P!'i表示第S步中所获得对应于 第i个区间参数的第j个输入样本点处的响应列向量，IZm(X)表示第m阶勒让德多项式的导 函数。
[0046] 第五步:如图2所示，将第一步中区间参数向量zi的中点列向量代入空腔噪声分 析的有限元模型中，计算响应列向量值名义值为pe。基于高维模型表征理论，根据第四步所 获得最佳平方逼近函数aU'i>(x)，建立响应列向量P的第1个分量在n维区间参数空间内的代 理模型A^(X),有：
[0047]
(峭
[004引其中n维列向量X= [X1，X2, . . .，Xn]T的第i个分量对应第i个区间参数，列向量X每 个元素的取值范围为[-1,1]，pe(l)表示响应列向量名义值F的第1个分量。
[0049] 第六步：如图2所示，根据顶点分析策略确定n维超正方体所有顶点的空间位置坐 标，并存储于顶点矩阵Mv中。利用如下格式遍历i和j，即：
[0050] i = l ,2, . . .，n且j = l，2, . . .，21-1 (11)
[0化1 ]顶点矩阵Mv的第i行第k列元素 Vik计算为：
[0化2]
(12)
[0053] 顶点矩阵Mv的每个列向量对应n维超正方体的一个顶点，对应于n维超正方体的2" 个顶点共计有2"个列向量。
[0054] 第屯步:如图2所示，基于第五步所获得的第1个响应分量的代理模型AW(X)与第 六步所获得的顶点矩阵Mv,将顶点矩阵Mv的第k列代入代理模型AW(X)中，计算获得n维超正 方体第k个顶点处响应列向量P第1个分量的离散值，有：
[0化5]
(13)
[0化6]根据顶点分析策略，响应列向量P第1个分量的区间界限好可W计算为：
[0化7]
I 14；
[005引第八步:如图2所示，利用第四步至第六步遍历响应分量索引值1，可W计算得到响 应列向量P的下界列向量片和上界列向量PU，最终可W确定空腔噪声响应的波动范围为Pi = [pL，pU]。
[0059] W图3所示的圆柱壳体为对象，由其围合而成的圆柱形空腔内部充满一定溫度范 围内的均匀空气介质，空腔两个端面为刚性壁面，壳体环向壁面为弹性圆柱壳体结构与空 气介质声场的禪合界面，边界条件为圆柱体壳两端固支，载荷条件为圆柱体壳周向外侧施 加简谐激励场。根据工程需求确定圆柱壳体空腔噪声分析的响应列向量P为给定频率点 (180赫兹、196赫兹和220赫兹)处空腔中线沿Z向各位置处的声压级，所考虑的区间参数列 在表1中，利用本发明【具体实施方式】的第一步至第八步可W预测圆柱形空腔噪声波动范围 如图4、图5和图6所示。
[0060] 表 1
[0062]本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
【主权项】
1. 一种基于高维模型表征与顶点分析策略的空腔噪声预测方法，其特征在于包括以下 步骤： 第一步:根据工程需求确定空腔噪声分析的响应列向量P所包含的响应分量，根据研究 人员或者工程经验数据确定不确定参数列向量Z所包含的不确定参数，基于不确定参数列 向量Z的试验数据，利用区间模型将其定量化为区间参数列向量Z 1，其下界列向量为Z1，上界 列向量为Λ中点列向量为2°和半径列向量为 第二步:根据第一步中响应列向量P关于区间参数列向量Z1的非线性程度，确定响应列 向量P关于每个区间参数最佳平方逼近函数的阶数N、高斯积分点个数m及高斯积分点列向 量i = 利用第一步确定的区间参数列向量Z1的中点列向量^、半径列向量f及 高斯积分点列向量_i :，.对区间参数列向量z进行抽样，将样本点存储于输入样本点矩阵Mi中； 第三步:将第二步中所得到的输入样本点矩阵此逐行代入空腔噪声分析的有限元模型 中，计算响应列向量P在每个样本点处的值，并存储于输出样本点矩阵Mo; 第四步:根据第三步所获得输出样本点矩阵Mo建立响应列向量P的第1个分量关于第i个 区间参数的最佳平方逼近函数Ααι)(χ); 第五步:将第一步中区间参数向量Z1的中点列向量f代入空腔噪声分析的有限元模型 中，计算响应列向量值名义值为PS基于高维模型表征理论，根据第四步所获得最佳平方逼 近函数A(1，U(X)，建立响应列向量P的第1个分量在η维区间参数空间内的代理模型A (1) (X); 第六步:根据顶点分析策略确定η维超正方体所有顶点的空间位置坐标，并存储于顶点 矩阵Mv中； 第七步:基于第五步所获得的第1个响应分量的代理模型Α(1)(χ)与第六步所获得的顶 点矩阵Μν，将顶点矩阵Mv的第k列代入代理模型Α(1)(χ)中，计算获得η维超正方体第k个顶点 处响应列向量P第1个分量的离散值i?^，进一步根据顶点分析策略确定响应列向量P第1个 分量的区间界限Pi1; 第八步:利用第四步至第六步遍历响应分量索引值1，可以计算得到响应列向量P的下 界列向量P1和上界列向量Pu，最终可以确定空腔噪声响应的波动范围为P1 = P,Pu]。2. 根据权利要求1所述的基于高维模型表征与顶点分析策略的空腔噪声预测方法，其 特征在于，所述方法以区间模型实现了空腔噪声预测过程中不确定性的定量化。3. 根据权利要求1所述的基于高维模型表征与顶点分析策略的空腔噪声预测方法，其 特征在于，所述方法基于最佳平方逼近理论与高维模型表征理论建立了空腔噪声响应关于 区间参数的代理模型。4. 根据权利要求1所述的基于高维模型表征与顶点分析策略的空腔噪声预测方法，其 特征在于，所述方法构造了任意η维超正方体所有顶点的空间坐标计算方式。5. 根据权利要求1所述的基于高维模型表征与顶点分析策略的空腔噪声预测方法，其 特征在于，所述方法利用基于高维模型表征的代理模型与顶点矩阵计算空腔噪声响应的区 间波动范围。
【文档编号】G06F19/00GK106021847SQ201610280321
【公开日】2016年10月12日
【申请日】2016年4月28日
【发明人】邱志平, 许孟辉, 王晓军, 王冲, 李云龙, 王磊, 仇翯辰, 陈贤佳, 郑宇宁
【申请人】北京航空航天大学