CN113418594A - 基于动态力识别的振动噪声贡献量分析方法及车辆 - Google Patents

基于动态力识别的振动噪声贡献量分析方法及车辆 Download PDF

Info

Publication number
CN113418594A
CN113418594A CN202110674846.3A CN202110674846A CN113418594A CN 113418594 A CN113418594 A CN 113418594A CN 202110674846 A CN202110674846 A CN 202110674846A CN 113418594 A CN113418594 A CN 113418594A
Authority
CN
China
Prior art keywords
suspension
vehicle
target point
point
sensor
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.)
Pending
Application number
CN202110674846.3A
Other languages
English (en)
Inventor
张健
贾文宇
马东正
范大力
李沙
宫世超
贾小利
尹冬冬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing Changan Automobile Co Ltd
Original Assignee
Chongqing Changan Automobile Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing Changan Automobile Co Ltd filed Critical Chongqing Changan Automobile Co Ltd
Priority to CN202110674846.3A priority Critical patent/CN113418594A/zh
Publication of CN113418594A publication Critical patent/CN113418594A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M17/00Testing of vehicles
    • G01M17/007Wheeled or endless-tracked vehicles
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种基于动态力识别的振动噪声贡献量分析方法,包括:S1.确定待测车辆的目标点,在目标点处布置传感器;S2.基于驾乘场景对待测车辆进行试验,并收集传感器采集的工况数据GR;S3.将待测车辆的悬架拆卸,并将由金属材质制成的悬架替换件按照设定扭力值安装到原悬架所在位置处;S4.采用锤击法对拆卸掉悬架后的待测车辆进行测试,获得悬架安装点至目标点处的频响函数HR;S5.根据工况数据GR和频响函数HR建立动态载荷目标函数,运用迭代法求解悬架接附点动态载荷F;S6.根据悬架接附点动态载荷F与频响函数HR计算各路径点贡献量,以此进行主路径识别与目标点响应预测。本发明能规避矩阵求逆带来的病态性问题,能够提高工作效率,且通用性强,可操作性高。

Description

基于动态力识别的振动噪声贡献量分析方法及车辆
技术领域
本发明属于动态载荷力识别技术领域,具体涉及基于动态力识别的振动噪声贡献量分析方法。
背景技术
汽车NVH(噪声、振动与声振粗糙度)性能不仅是乘坐舒适性的重要指标,同时也是提升汽车品质感与强化品牌力的关键因素。因此,各大汽车厂商在车型研发中投入了大量人力与物力对其进行研究。
如专利文献CN106644512A公开了一种基于动力总成载荷的噪声分析方法和系统,方法包括以下步骤:获取传声器位置信息和三向加速度传感器位置信息;基于不同工况对待测车辆进行道路试验以获取传声器和三向加速度传感器的道路试验测试数据;将待测车辆的动力总成拆卸后,将以金属零件改制的动力总成悬置安装在拆卸后的待测车辆,并保留设置的传声器和三向加速度传感器;采用锤击法对拆卸后的待测车辆进行测试以获取传声器和三向加速度传感器的频率响应函数;根据道路试验测试数据和频率响应函数进行动力总成载荷计算,以获取动力总成载荷矩阵;根据动力总成载荷矩阵预测待测车辆在实际道路工况下的车内噪声,从而能综合考虑各个方向及多个动力总成悬架相互之间的影响。
又如专利文献CN106679990A公开了一种汽车六自由度轮心力测试及振动噪声贡献率的计算方法,该方法包括:建立坐标系、布置测点及安装传感器、布置激励点及测量激励点坐标、测量频率响应函数、求解从轮毂中心到响应点的频率响应函数矩阵、求解从轮毂中心到目标点的频率响应函数向量、路试、求解激励功率谱密度矩阵、验证激励功率谱密度矩阵识别有效性、计算各自由度激励对目标点响应的能量贡献率。
以上专利文献中对动态载荷识别皆采用频响函数矩阵求逆法,解决了矩阵求逆病态性成了首要难题。矩阵求逆病态性问题的解决途径通常分为两大类,分别为参考点的三向振动传感器数量和正则化方法。矩阵求逆法计算动态载荷需在每个路径下布置三向振动传感器,并且为规避矩阵病态性问题,要求传感器的数量至少是被识别激励力数量的2倍及以上,造成传感器数量需求大,数据采集系统通道需求多。同时对三向振动传感器布置位置有一定要求,同一路径点处各参考点三向振动传感器位置要适宜,参考点三向振动传感器距离路径点太近容易获取重复振动信息,造成矩阵病态性,距离太远获取信息不纯,引入误差。对实施者而言无疑带来一定难度,要同时兼顾位置与空间布置空间。在后处理计算中,要运用正则化方法修正矩阵病态性,目前商业软件如LMS.TSET LAB测试系统大多采用奇异值截断法修正矩阵病态性,对于奇异值的设定多凭经验,且需反复调试。
因此,有必要开发一种新的基于动态力识别的振动噪声贡献量分析方法及车辆。
发明内容
本发明的目的是提供基于动态力识别的振动噪声贡献量分析方法,能规避矩阵求逆带来的病态性问题,能节省设备资源,提高工作效率,且通用性强,可操作性高。
本发明所述的基于动态力识别的振动噪声贡献量分析方法,包括以下步骤:
S1.确定待测车辆的目标点,在目标点处布置传感器,若目标点需要解决声噪问题,则在该目标点处安装麦克风传感器,若目标点处需要解决振动问题,则在该目标点处布置三向振动传感器;
S2.基于驾乘场景对待测车辆进行试验,并收集传感器采集的工况数据GR
S3.将待测车辆的悬架拆卸,并将由金属材质制成的悬架替换件按照设定扭力值安装到原悬架所在位置处;
S4.采用锤击法对拆卸掉悬架后的待测车辆进行测试,获得悬架安装点至目标点处的频响函数HR
S5.根据步骤S2所得的工况数据GR和步骤S4所得的频响函数HR建立动态载荷目标函数,运用迭代法求解悬架接附点动态载荷F;
S6.根据悬架接附点动态载荷F与步骤S4所得的频响函数HR计算各路径点贡献量,以此进行主路径识别与目标点响应预测。
可选地,所述传感器方向依据整车坐标系O-XYZ的轴线方向布置,车头至车尾方向为+X,主驾至副驾方向为+Y,竖直向上为+Z。
可选地,若目标点需要解决声噪问题,则在该目标点处安装麦克风传感器,若目标点处需要解决振动问题,则在该目标点处布置三向振动传感器。
本发明具有以下优点:利用路径点至目标点的频响函数和目标点的工况数据,通过迭代法求解路径点动态载荷。无需矩阵求逆,因此规避了矩阵求逆带来的病态性问题。无需在路径点附近布置参考点三向振动传感器,对数据采集前端通道需求少,大大节省了设备资源,同时省去了布置大量三向振动传感器工作量,提升了工作效率;此方法在静态频响函数测试时拆除了路径点主被动侧,特别适用于同一车身状态下匹配不同状态悬架的噪声振动分析,每次仅采集目标点传感器工况数据即可重复本方法,通用性强,可操作性高。
附图说明
图1为本实施例的流程图;
图2为本实施例中迭代算法的流程图。
具体实施方式
下面结合附图对本发明作可选的说明。
如图1所示,本实施例中,一种基于动态力识别的振动噪声贡献量分析方法,包括以下步骤:
S1.布置目标点传感器:确定待测车辆的目标点,在目标点处布置传感器,若目标点需要解决声噪问题,则在该目标点处安装麦克风传感器,若目标点处需要解决振动问题,则在该目标点处布置三向振动传感器。比如:若要解决乘员舱的声噪问题,可以选择在乘员舱座椅人耳附近位置布置麦克风传感器。若要解决座椅的振动问题,可以选择在座椅导轨处布置三向振动传感器。
本实施例中,传感器方向依据整车坐标系O-XYZ(车头至车尾方向为+X,主驾至副驾方向为+Y,竖直向上为+Z)的轴线方向布置。
S2.采集工况数据:基于驾乘场景对待测车辆进行试验,并收集传感器采集的工况数据GR
S3.拆卸悬架,安装悬架替换件:为了消除结构耦合问题,在进行静态频响函数测试前,将待测车辆的悬架拆卸,并将由金属材质制成的悬架替换件(比如:金属套筒)按照设定扭力值安装到原悬架所在位置处。
S4.采用锤击法对拆卸掉悬架后的待测车辆进行测试,获得悬架安装点至目标点处的频响函数HR
S5.根据步骤S2所得的工况数据GR和步骤S4所得的频响函数HR建立动态载荷目标函数,运用迭代法求解悬架接附点动态载荷F。
S6.根据悬架接附点动态载荷F与步骤S4所得的频响函数HR计算各路径点贡献量,以此进行主路径识别与目标点响应预测。
本实施例中,所述步骤S5,以解决声噪为例进行说明:
根据声音传递特性,车内声音由如下公式获得:
GR=HR·F;
对于求解动态载荷F(即求解出的动态载荷的最优值),即寻找Fd满足目标min||GR-HR·Fd||,相当于Fd满足min(GR-HR·F)2,其中,Fd为根据传感器测量的数据计算出的实际动态载荷;令U(F)=(GR-HR·F)2,即得到:
Figure BDA0003120311570000041
其中,U(F)是Rn上具有一阶连续偏导数的函数,Rn是每次迭代实际产生的误差值,选用梯度下降法对此进行迭代优化,选定适当的初值F(0),不断迭代,更新F的值,进行目标函数的极小化,直到收敛。
如图2所示,以
Figure BDA0003120311570000042
为搜索方向,k为迭代次数,λk为第k次迭代的步长,由一维搜索确定,即λk使得:
Figure BDA0003120311570000043
迭代算法流程如下:
输入:目标函数U(F),搜索方向
Figure BDA0003120311570000044
设置精度ε;
输出:U(F)的极小点F*
(5.1)取初始值F(0)∈Rn,置k=0;
(5.2)计算U(F(k));
(5.3)计算梯度
Figure BDA0003120311570000045
Figure BDA0003120311570000046
时,停止迭代,令F*=F(k);否则,以
Figure BDA0003120311570000047
为搜索方向,使
Figure BDA0003120311570000048
(5.4)置
Figure BDA0003120311570000049
计算U(F(k+1)),当||U(F(k+1))-U(F(k))||<ε或||F(k+1)-F(k)||<ε,停止迭代,令F*=F(k+1),即F=F*;否则,置k=k+1,转入(5.3)。
根据悬架接附点动态载荷F与静态测试频响函数HR计算各路径点至车内的贡献量,以此进行主路径识别,为问题整改与优化提供参考。此外,还可对悬架接附点动态载荷F与静态测试频响函数HR进行调制,预测目标点响应。

Claims (3)

1.一种基于动态力识别的振动噪声贡献量分析方法,其特征在于:包括以下步骤:
S1.确定待测车辆的目标点,在目标点处布置传感器;
S2.基于驾乘场景对待测车辆进行试验,并收集传感器采集的工况数据GR
S3.将待测车辆的悬架拆卸,并将由金属材质制成的悬架替换件按照设定扭力值安装到原悬架所在位置处;
S4.采用锤击法对拆卸掉悬架后的待测车辆进行测试,获得悬架安装点至目标点处的频响函数HR
S5.根据步骤S2所得的工况数据GR和步骤S4所得的频响函数HR建立动态载荷目标函数,运用迭代法求解悬架接附点动态载荷F;
S6.根据悬架接附点动态载荷F与步骤S4所得的频响函数HR计算各路径点贡献量,以此进行主路径识别与目标点响应预测。
2.根据权利要求1所述的基于动态力识别的振动噪声贡献量分析方法,其特征在于:所述传感器方向依据整车坐标系O-XYZ的轴线方向布置,车头至车尾方向为+X,主驾至副驾方向为+Y,竖直向上为+Z。
3.根据权利要求1或2所述的基于动态力识别的振动噪声贡献量分析方法,其特征在于:若目标点需要解决声噪问题,则在该目标点处安装麦克风传感器,若目标点处需要解决振动问题,则在该目标点处布置三向振动传感器。
CN202110674846.3A 2021-06-17 2021-06-17 基于动态力识别的振动噪声贡献量分析方法及车辆 Pending CN113418594A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110674846.3A CN113418594A (zh) 2021-06-17 2021-06-17 基于动态力识别的振动噪声贡献量分析方法及车辆

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110674846.3A CN113418594A (zh) 2021-06-17 2021-06-17 基于动态力识别的振动噪声贡献量分析方法及车辆

Publications (1)

Publication Number Publication Date
CN113418594A true CN113418594A (zh) 2021-09-21

Family

ID=77789012

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110674846.3A Pending CN113418594A (zh) 2021-06-17 2021-06-17 基于动态力识别的振动噪声贡献量分析方法及车辆

Country Status (1)

Country Link
CN (1) CN113418594A (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104537251A (zh) * 2015-01-06 2015-04-22 西安交通大学 一种风机叶片的冲击载荷识别方法
CN106644512A (zh) * 2016-12-22 2017-05-10 北京汽车股份有限公司 基于动力总成载荷的噪声分析方法和系统
KR102211079B1 (ko) * 2019-09-23 2021-02-02 넥센타이어 주식회사 타이어에 의한 차량의 실내 소음 예측 장치 및 방법

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104537251A (zh) * 2015-01-06 2015-04-22 西安交通大学 一种风机叶片的冲击载荷识别方法
CN106644512A (zh) * 2016-12-22 2017-05-10 北京汽车股份有限公司 基于动力总成载荷的噪声分析方法和系统
KR102211079B1 (ko) * 2019-09-23 2021-02-02 넥센타이어 주식회사 타이어에 의한 차량의 실내 소음 예측 장치 및 방법

Similar Documents

Publication Publication Date Title
CN107393522B (zh) 选择车辆上主动道路噪声控制的传感器位置的方法和系统
Abdullah et al. Dynamics properties of a Go-kart chassis structure and its prediction improvement using model updating approach
CN112597595B (zh) 汽车车内结构噪声诊断及优化方法
CN111444623B (zh) 一种阻尼非线性商用车悬架动力学的协同优化方法及系统
CN106066921A (zh) 一种整车路面振动噪声的cae仿真预测方法
CN111680409B (zh) 一种汽车结构耐久程序的试验场关联方法
CN112595528A (zh) 车辆行驶工况下动力总成悬置主动侧激励力的识别方法
CN112697448A (zh) 车辆怠速工况下动力总成悬置主动侧激励力的识别方法
CN110715783A (zh) 一种长头载重汽车驾驶室振动特性的检测方法
Kim et al. Optimization of the lower arm of a vehicle suspension system for road noise reduction by sensitivity analysis
CN115270296A (zh) 一种商用车驾驶室疲劳耐久性分析方法及系统
CN113886974B (zh) 一种车内结构声路噪预测方法
CN113987678B (zh) 一种基于虚拟轮心力的路噪混合预测及优化方法
CN113515808B (zh) 一种动力总成的噪声预测方法
CN113418594A (zh) 基于动态力识别的振动噪声贡献量分析方法及车辆
Song et al. A methodolgy for evaluating the structure-borne road noise prior to a prototype vehicle using direct force measured on a suspension rig
Olatunbosun et al. Dynamic performance analysis of a light van body-in-white structure
CN114936424A (zh) 一种车辆橡胶件多自由度分析方法和相关装置
Minervini et al. FBS Decoupling at Suspension Level for Road Noise Applications
CN110263408B (zh) 一种利用bni曲线评估ntf风险的方法
CN111241724B (zh) 一种宽体矿用车车架疲劳寿命预测方法
CN113591228B (zh) 一种噪声传递路径分析及优化方法
CN114896708B (zh) 一种车架质量属性量化方法及系统
CN116628855A (zh) 基于对标的车身悬置设计方法、装置、系统、设备及介质
Buchmann Krivachy et al. Virtual operational strength analysis for commercial vehicles

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210921