CN110263408B - 一种利用bni曲线评估ntf风险的方法 - Google Patents
一种利用bni曲线评估ntf风险的方法 Download PDFInfo
- Publication number
- CN110263408B CN110263408B CN201910511121.5A CN201910511121A CN110263408B CN 110263408 B CN110263408 B CN 110263408B CN 201910511121 A CN201910511121 A CN 201910511121A CN 110263408 B CN110263408 B CN 110263408B
- Authority
- CN
- China
- Prior art keywords
- curve
- ntf
- bni
- risk
- vehicle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M17/00—Testing of vehicles
- G01M17/007—Wheeled or endless-tracked vehicles
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供了一种利用BNI曲线评估NTF风险的方法,包括如下步骤:1)选取若干参考车型,利用TPA方法获得参考车型底盘安装点点的载荷数据;2)根据得到的载荷数据,获得研发车型各安装点的NTF曲线的加权系数;3)根据研发车型在研发过程中得到的各安装点的NTF曲线,对NTF曲线进行加权求和得到BNI曲线;4)当研发车型的NTF曲线部分频率超过目标线,则进一步判断得到的BNI曲线在此频率范围是否存在峰值,如果不存在峰值,则风险很低,如果存在峰值,则必须进行优化使NTF达到目标线要求。本发明利用BNI曲线对研发车型的未达到目标线的NTF曲线进行二次评价,可以减少未达标NTF曲线二次风险评估对工程师经验的依赖性,评估结果更加客观。
Description
技术领域
本发明属于汽车NVH领域,尤其是涉及一种利用BNI曲线评估NTF风险的方法。
背景技术
汽车NVH(Noise、Vibration、Harshness)性能开发,即是降低车辆运行中的噪声、振动,提高乘坐舒适性。在竞争激烈的汽车市场上,消费者对乘坐舒适性越来越关注,提高车辆NVH水平已成为各大主机厂重要竞争和技术发展方向。在此背景下,车辆的NVH性能正在成为汽车研发过程中最为重要的性能之一。
NTF(Noise Transfer Function)即激励点到响应点的声振传递函数,是车身声学响应对激励点敏感度的反映。在实际的研发和应用过程中,通过底盘与车身各安装点的NTF曲线,评价车身声学响应灵敏度。整备车身的NTF曲线是否合理,是车辆NVH性能开发过程重要评价点。
评价NTF曲线时,一般先通过标杆车或者数据库中同类车型NTF数据设定研发车型各安装点NTF目标曲线,再将各安装点NTF曲线与目标线对比,超出目标线部分认为未达到设计标准,需要进行优化。
在实际开发过程中,少有车型所有NTF曲线可以都达到设定的目标值,这和设定目标值合理性以及所研发车型的结构特点有关,对于未达标NTF曲线,还需要进行风险评估,评估其设计风险。对于风险评估,与工程师对NTF的理解、个人经验等因素有关,存在一定的主观因素。
发明内容
有鉴于此,本发明旨在提出一种利用BNI曲线评估NTF风险的方法,以解决未达标NTF曲线风险评估受工程师主观因素影响问题。
为达到上述目的,本发明的技术方案是这样实现的:
一种利用BNI曲线评估NTF风险的方法,其特征在于包括如下步骤:
1)选取若干参考车型,利用TPA方法获得参考车型底盘安装点的载荷数据;
2)根据得到的载荷数据,获得研发车型各安装点的NTFi,j曲线的加权系数ai,j;
3)根据研发车型在研发过程中得到的各安装点的NTFi,j曲线,结合如下公式获得研发车型的BNI曲线;
式中:NTFi,j为第i个激励点对应的j方向产生的NTF曲线,安装点即为激励点;
4)当研发车型的NTF曲线部分频率超过目标线,则进一步判断得到的BNI曲线在此频率范围是否存在峰值,如果不存在峰值,则风险很低,如果存在峰值,则必须进行优化使NTF达到目标线要求。
进一步的,在步骤2中,将步骤1得到的载荷数据fi,j,n按照从大到小进行排序,其中,fi,j,n表示第n辆车在第i个激励点j方向的激励载荷;取上四分位数得到Fi,j,找出最大的上四分位数Fmax,将Fi,j除以Fmax做归一化处理,通过如下公式得到加权系数ai,j:
进一步的,在步骤1中利用TPA方法获得参考车型底盘安装点的载荷数据的方法包括如下步骤:
11)利用力锤激励法获取结构声学传递函数;
12)用逆矩阵方法求取载荷。
相对于现有技术,本发明具有以下优势:
本发明利用BNI曲线对研发车型的未达到目标线的NTF曲线进行二次评价,可以减少未达标NTF曲线二次风险评估对工程师经验的依赖性,评估结果更加客观。BNI曲线可以综合所有安装点对车内噪声进行评价,而不是单单从某一个安装点角度进行评价,BNI曲线相比车辆NTF曲线趋势上更加接近路噪实测曲线。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为本发明实施例得到的BNI曲线;
图2为本发明研发车型的BNI曲线、NTF曲线和路噪实测曲线对比图。
具体实施方式
需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。
下面将参考附图并结合实施例来详细说明本发明。
本实施例一种利用BNI曲线评估NTF风险的方法,以SUV车型为例,包括如下步骤:
步骤1,选取若干参考车型,利用TPA方法获得参考车型底盘安装点的载荷数据;
步骤2,根据得到的载荷数据,获得研发车型各安装点的NTF曲线的加权系数;
步骤3,根据研发车型在研发过程中得到的各安装点的NTFi,j曲线,结合如下公式获得研发车型的BNI曲线;
式中:NTFi,j为第i个激励点对应的j方向产生的NTF曲线,安装点即为激励点;
步骤4,当研发车型的NTF曲线部分频率超过目标线,则进一步判断得到的BNI曲线在此频率范围是否存在峰值,如果不存在峰值,则风险很低,如果存在峰值,则必须进行优化使NTF达到目标线要求。
在步骤1中:
首先选取参考车型,选取参考车型的目的在于获取其底盘安装点的载荷数据,通过载荷数据获得BNI曲线的加权系数。选取满足以下条件的SUV车型:
(1)参考车型与研发车型整车重量相差5%以内
(2)参考车型前后悬架及副车架结构形式与研发车型基本一致;
(3)参考车型轴距与研发车型相差5%以内;
(4)至少6款参考车型。
之后,获得参考车型底盘安装点的载荷数据,本实施例利用TPA技术获取参考车型的底盘安装点的载荷数据。
(1)TPA的基本理论:TPA(Transfer Path Analysis)传递路径分析方法,分析车内总体噪声与每条路径贡献关系。车内所接受的噪声是结构声和空气声之和。将所有的结构声和空气声传递路径的贡献和车内总的噪声绘制在一张图上,就可以知道每条路径的具体贡献量。本申请以分析结构声传递路径贡献量为例进行方法说明:
任意一个振动源(fi)激励车身上的某一点i,都会在车内产生结构噪声。对固定的激励点和固定的声音接受点,其车内噪声与激励源振动之间的传递函数是确定的。第i个振动源在j方向的力fi,j产生车内声音响应与激励和传递函数之间的关系为:
车身上的N个振动源作用在车身N个不同的位置或方向,都会产生结构声。多个振动源通过N条传递路径产生的结构声的总和,表达为:
式中PSB(ω)表示车内结构噪声响应,Hi,j SB(ω)表示结构声振传递函数,fi,j(ω)表示激励载荷。
(2)测量结构声振传递函数:
利用力锤激励法获取结构声振传递函数。用带有力传感器的力锤敲击激励点,将传声器放置于响应点用来测量声音响应,将测量的声音和振动信号进行处理后,即可以得到结构声振传递函数。
(3)获取载荷
用逆矩阵方法求取载荷,逆矩阵法主要应用于线性系统,有响应Xi到输入Fj的频响函数Hij=Xi/Fj,故激励力可用矩阵求逆法求得:
上式(3)中:
{FN}为耦合激励力向量;{XM}为响应点上的工作响应向量;HMN -1为由输入FN到响应XM的频响函数逆矩阵。
采用逆矩阵法测量激励力时,必须要实测路况工作条件下测量指示点X/Y/Z三个方向的加速度。
需要注意的是,根据式(3)计算工作激励力过程中为了得到唯一解,需要响应数量(m)至少需要与估计的输入数量(n)相等。并且在测量过程中,可以在被动侧安装多个加速度传感器,所以在接受测可以测量更多数量的响应(m>n),即允许方程组超定解,我们使用最小二乘法获得工作力的估计值,一般使用超定因子2是结果获得更高的置信度,即在被动侧加速度传感器至少安装两个以上。
还有,因为参考车型试验周期长、费用高,车辆开发过程中获得所有参考车的载荷数据往往时间不允许可以通过NVH咨询类公司获得某些参考车型载荷数据。
在步骤2中:
加权系数ai,j,是算法BNI的核心。ai,j的获取基于所选择参考车型底盘安装点载荷数据,安装点即是激励点。根据载荷数据确定第i个激励点j方向下NTFi,j的系数ai,j,具体算法如下:
表1:
如表1所示,n表示参考车型数量,fi,j,n表示第n辆车在i激励点j方向的激励载荷。将载荷数据按照从大到小进行排序,取上四分位数得到Fi,j。找出最大的上四分位数Fmax,将Fi,j除以Fmax做归一化处理,使获得的加权系数ai,j在0和1之间;
由上式(5)得到研发车型NTFi,j曲线的加权系数ai,j。
在步骤3中:
获得设研发车型的BNI曲线:
在研发过程中获取的研发车型各安装点的NTF曲线,根据式(4)对NTFi,j曲线进行加权求和,得到BNI曲线,如图1所示。
式(4)中,ai,j为加权系数,NTFi,j为第i个激励点对应的j方向产生的NTF曲线。
本申请提出用BNI曲线二次评价未达标NTF曲线的方法。对未达标的NTF曲线二次评价时,如果满足BNI曲线评价标准,则认为设计风险低,反之,则认为设计风险高,需要优化NTF曲线至满足BNI曲线评价标准。BNI曲线是通过多辆参考车型安装点载荷数据统计规律及加权算法得到的客观的评价指标,所以通过BNI曲线对未达标NTF曲线进行二次评价来代替工程师主观评价可以规避工程师个人主观因素对评价结果的影响。
图2为本发明研发车型的BNI曲线、NTF曲线和路噪实测曲线对比图,相比NTF曲线,BNI曲线与路噪实测曲线趋势吻合的更好。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种利用BNI曲线评估NTF风险的方法,其特征在于包括如下步骤:
1)选取若干参考车型,利用TPA方法获得参考车型连接点的载荷数据;
2)根据得到的载荷数据,获得研发车型各安装点的NTFi,j曲线的加权系数ai,j;其中,将步骤1得到的载荷数据fi,j,n按照从大到小进行排序,其中,fi,j,n表示第n辆车在第i个激励点j方向的激励载荷;取上四分位数得到Fi,j,找出最大的上四分位数Fmax,将Fi,j除以Fmax做归一化处理,通过如下公式得到加权系数ai,j:
3)根据研发车型在研发过程中得到的各安装点的NTFi,j曲线,结合如下公式获得研发车型的BNI曲线;
式中:NTFi,j为第i个激励点对应的j方向产生的NTF曲线,安装点即为激励点;
4)当研发车型的NTF曲线部分频率超过目标线,则进一步判断得到的BNI曲线在此频率范围是否存在峰值,如果不存在峰值,则风险很低,如果存在峰值,则必须进行优化使NTF达到目标线要求。
2.根据权利要求1所述的一种利用BNI曲线评估NTF风险的方法,其特征在于:在步骤1中利用TPA方法获得参考车型连接点的载荷数据的方法包括如下步骤:
11)利用力锤激励法获取结构声学传递函数;
12)用逆矩阵方法求取载荷。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910511121.5A CN110263408B (zh) | 2019-06-13 | 2019-06-13 | 一种利用bni曲线评估ntf风险的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910511121.5A CN110263408B (zh) | 2019-06-13 | 2019-06-13 | 一种利用bni曲线评估ntf风险的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110263408A CN110263408A (zh) | 2019-09-20 |
CN110263408B true CN110263408B (zh) | 2023-04-25 |
Family
ID=67918064
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910511121.5A Active CN110263408B (zh) | 2019-06-13 | 2019-06-13 | 一种利用bni曲线评估ntf风险的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110263408B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113656934B (zh) * | 2021-07-02 | 2024-06-21 | 中汽研(天津)汽车工程研究院有限公司 | 一种基于传递路径分析理论的车内空气噪声目标分解方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101384463A (zh) * | 2006-03-31 | 2009-03-11 | 丰田自动车株式会社 | 车辆动态控制系统和对车辆动态进行控制的方法 |
CN104786957A (zh) * | 2015-04-29 | 2015-07-22 | 上汽通用五菱汽车股份有限公司 | 汽车阻尼板结构优化设计方法 |
CN105004416A (zh) * | 2015-07-02 | 2015-10-28 | 哈尔滨工程大学 | 基于逆边界元法机械噪声远场声压预测方法 |
CN105799487A (zh) * | 2016-03-15 | 2016-07-27 | 北京理工大学 | 一种机电复合传动履带车辆多目标优化方法 |
CN105825014A (zh) * | 2016-03-17 | 2016-08-03 | 中铁大桥勘测设计院集团有限公司 | 一种基于车桥耦合分析的车桥安全健康评估系统及方法 |
CN107085768A (zh) * | 2017-04-25 | 2017-08-22 | 交通运输部公路科学研究所 | 一种用于评价汽车使用可靠性的系统及方法 |
CN107132052A (zh) * | 2017-05-05 | 2017-09-05 | 广州汽车集团股份有限公司 | 电动车辆车内振动、噪声预测方法和装置 |
CN108509754A (zh) * | 2018-05-15 | 2018-09-07 | 北京航空航天大学 | 一种高密度封装集成电路键合丝触碰风险评估方法 |
CN109086536A (zh) * | 2018-08-13 | 2018-12-25 | 天津龙创世纪汽车设计有限公司 | 一种基于形貌优化的ntf改进方法 |
CN109753722A (zh) * | 2018-12-29 | 2019-05-14 | 北京长城华冠汽车技术开发有限公司 | 一种提高噪声传递函数的优化设计方法和系统 |
-
2019
- 2019-06-13 CN CN201910511121.5A patent/CN110263408B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101384463A (zh) * | 2006-03-31 | 2009-03-11 | 丰田自动车株式会社 | 车辆动态控制系统和对车辆动态进行控制的方法 |
CN104786957A (zh) * | 2015-04-29 | 2015-07-22 | 上汽通用五菱汽车股份有限公司 | 汽车阻尼板结构优化设计方法 |
CN105004416A (zh) * | 2015-07-02 | 2015-10-28 | 哈尔滨工程大学 | 基于逆边界元法机械噪声远场声压预测方法 |
CN105799487A (zh) * | 2016-03-15 | 2016-07-27 | 北京理工大学 | 一种机电复合传动履带车辆多目标优化方法 |
CN105825014A (zh) * | 2016-03-17 | 2016-08-03 | 中铁大桥勘测设计院集团有限公司 | 一种基于车桥耦合分析的车桥安全健康评估系统及方法 |
CN107085768A (zh) * | 2017-04-25 | 2017-08-22 | 交通运输部公路科学研究所 | 一种用于评价汽车使用可靠性的系统及方法 |
CN107132052A (zh) * | 2017-05-05 | 2017-09-05 | 广州汽车集团股份有限公司 | 电动车辆车内振动、噪声预测方法和装置 |
CN108509754A (zh) * | 2018-05-15 | 2018-09-07 | 北京航空航天大学 | 一种高密度封装集成电路键合丝触碰风险评估方法 |
CN109086536A (zh) * | 2018-08-13 | 2018-12-25 | 天津龙创世纪汽车设计有限公司 | 一种基于形貌优化的ntf改进方法 |
CN109753722A (zh) * | 2018-12-29 | 2019-05-14 | 北京长城华冠汽车技术开发有限公司 | 一种提高噪声传递函数的优化设计方法和系统 |
Non-Patent Citations (1)
Title |
---|
整备车身的CAE 解析与试验;徐仰汇等;《第六届中国CAE工程分析技术年会论文集》;20100728;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110263408A (zh) | 2019-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10013967B2 (en) | Method and system for selecting sensor locations on a vehicle for active road noise control | |
CN112697448B (zh) | 车辆怠速工况下动力总成悬置主动侧激励力的识别方法 | |
CN112595528B (zh) | 车辆行驶工况下动力总成悬置主动侧激励力的识别方法 | |
CN112100816A (zh) | 一种基于电机声学模型的电动车车内噪声预测方法 | |
CN110263408B (zh) | 一种利用bni曲线评估ntf风险的方法 | |
CN110688956A (zh) | 一种用于汽车道路噪声主动控制的参考信号选择方法 | |
Park et al. | Coherence-based sensor set expansion for optimal sensor placement in active road noise control | |
Kim et al. | Optimization of the lower arm of a vehicle suspension system for road noise reduction by sensitivity analysis | |
CN115186377A (zh) | 一种商用车悬置系统的稳健性优化方法 | |
CN115270296A (zh) | 一种商用车驾驶室疲劳耐久性分析方法及系统 | |
CN113886974B (zh) | 一种车内结构声路噪预测方法 | |
EP2574951A2 (en) | A method for analyzing sound transmission paths in a system | |
CN114329788A (zh) | 一种基于田口法和熵权灰色关联分析的车门优化设计方法 | |
Tang et al. | Operational transfer path analysis with regularized total least-squares method | |
CN113515808A (zh) | 一种动力总成的噪声预测方法 | |
Peng et al. | Research on the virtual reality of vibration characteristics in vehicle cabin based on neural networks. | |
KR100929589B1 (ko) | Mts 기법을 이용한 음질 평가 방법 | |
Dom et al. | Transfer path analysis: accurate load prediction beyond the traditional mount stiffness and matrix inversion methods | |
CN106777781B (zh) | 隔振件的隔振性能分析方法、刚度确定方法及装置 | |
CN115168999A (zh) | 一种路噪tpa优化方法及装置 | |
Janssens et al. | A novel transfer path analysis method delivering a fast and accurate noise contribution assessment | |
Lennström | In-situ characterization of vibrations from a door mounted loudspeaker | |
van der Seijs et al. | Road noise: embedding suspension test benches in sound & vibration design using virtual points and the transfer path analysis framework | |
CN112528419B (zh) | 基于相关度的目标结构件与外部激励载荷数据库构建方法 | |
Kim et al. | A study on optimal selection of reference accelerometers for active road noise control |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |