CN105092194A - 一种基于逆子结构技术的传递路径分析方法 - Google Patents
一种基于逆子结构技术的传递路径分析方法 Download PDFInfo
- Publication number
- CN105092194A CN105092194A CN201510553126.6A CN201510553126A CN105092194A CN 105092194 A CN105092194 A CN 105092194A CN 201510553126 A CN201510553126 A CN 201510553126A CN 105092194 A CN105092194 A CN 105092194A
- Authority
- CN
- China
- Prior art keywords
- point
- signal
- path analysis
- analysis method
- frequency response
- 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
Links
Abstract
本发明公开了一种基于逆子结构技术的传递路径分析方法,包括以下步骤:步骤1、测取激励信号和响应信号;步骤2、测取频响函数;步骤3、计算动刚度;步骤4、计算载荷;步骤5、计算传递路径贡献量。具有克服了专用设备、测试结果与实车状态存在差异的不足,使计算所得传递路径贡献量更接近整车实际状态等优点。
Description
技术领域
本发明涉及一种振动噪声传递路径分析技术,具体为一种基于逆子结构技术的传递路径分析方法。
背景技术
车辆,尤其是混合动力汽车,振动源众多,有路面激励、轮胎激励、发动机激励、电机激励、排气系统激励等,这些振动源对车内等处的振动噪声都有影响,判断这些振动源对目标点振动噪声贡献量的大小是减振降噪的前提。目前常见用于分析振动源贡献量的传递路径分析方法主要有传统传递路径分析方法、工况传递路径分析方法,其它方法大多基于这两种方法改进。传统传递路径分析方法具有较高的精度,但分析过程耗时较多,一般需要30个工作日,工况传递路径分析方法则相反。
传递路径分析方法基于这样一个理论:目标位置的总响应由来自不同路径的贡献量线性叠加而得,比如,影响某个位置振动的传递路径贡献总量表达如式(1)
式中,ω为频率;Y(ω)为传递路径贡献总量;Hp(ω)为频响;fp(ω)为载荷;n为激励传递路径数量;p为传递路径序数。由式(1)可见,传统传递路径分析方法计算贡献总量需要得到频响函数和载荷。目前计算载荷的方法有逆矩阵法和悬置刚度法,对于悬置刚度法,获得悬置刚度的方法是试验实测。试验实测需要专用的测试设备,且测试结果受设备精度、操作人员业务水平、周边环境等因素影响,结果误差影响贡献量的准确性。试验时悬置被从实车中分离,与实车状态刚度特性有一定差距。
发明内容
本发明的目的在于克服现有技术的缺点与不足,提供一种基于逆子结构技术的传递路径分析方法,该方法是一种基于逆子结构技术计算悬置动刚度的传递路径分析方法。
本发明的目的可以通过以下技术方案实现:一种基于逆子结构技术的传递路径分析方法,包括以下步骤:
步骤1、测取激励信号和响应信号;
步骤2、测取频响函数;
步骤3、计算动刚度;
步骤4、计算载荷;
步骤5、计算传递路径贡献量。
所述步骤2包括以下步骤:
步骤21、将拟进行传递路径分析的系统从用弹性件连接的部位划分为2个子结构;
步骤22、确定激励点和响应点,一般取其中一个子结构上弹性连接部位为激励点;
步骤23、给2个子结构每一个弹性连接点确定互相垂直的X、Y和Z三个方向,且这些连接点的三个方向相同;
步骤24、在每个连接点布置加速度传感器,逐一在这些连接点用力锤从X、Y和Z三个方向施力,记录所有传感器的信号;
步骤25、从步骤22确定的激励点处施力,测取从此激励点到响应点的频响函数。
所述步骤3包括以下步骤:
步骤31、抽取步骤24所述力信号和加速度响应信号;
步骤32、由力信号和加速度信号计算频响函数;
步骤33、由频响函数计算动刚度。
本发明的目的也可以通过以下技术方案实现:一种基于逆子结构技术的传递路径分析方法,在进行传递路径分析测试时,首先将车辆运行在设置的工况,测取选定的激励点和目标点振动信号。接着根据需要,测取整车系统水平的频响函数。第3步按式(3)方法进行动刚度计算。第4步,得到了动刚度,依据式(4)即可计算载荷。最后,依照式(1)即可进行传递路径贡献量计算。
以下述方法计算悬置动刚度。一个系统分解为子结构A和子结构B,两个子结构之间通过弹性元件,比如发动机悬置,进行多点耦合,它们存在以下动力学关系
经演算,得到以下动刚度[Kc]计算式
式中,H代表频响函数;下标S代表此频响函数处于整车系统状态测量;下标c(a)、c(b)代表子结构A和B的耦合点。
相对于现有技术,本发明具有如下的优点与有益效果:
1、本发明的传递路径分析方法是先基于在线测量所得频响函数使用逆子结构技术计算悬置动刚度,再基于悬置主动端和被动端的振动量计算载荷,然后由频响函数和载荷计算传递路径贡献量。通常,动刚度获得方法为由试验实测,实测方法需要专用设备、测试结果与实车状态存在差异等,该发明的动刚度计算方法克服了这些不足,且计算所得传递路径贡献量更接近整车实际状态。
2、本发明利用逆子结构技术在线测量整车状态下悬置动刚度,进而计算传递路径贡献量,克服了传统传递路径分析方法需要离线测试悬置动刚度和整个分析过程耗时较多的缺陷。
附图说明
图1逆子结构传递路径分析方法流程。
图2逆子结构动力学关系示意图。
图3整车振动模型;图中,1代表发动机,2代表变速器,3代表发动机前悬置,4代表发动机后悬置,5代表变速器悬置,6代表麦克风。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例
如图1所示,为本发明的逆子结构传递路径分析方法流程;进行传递路径分析测试时,首先将车辆运行在设置的工况,测取选定的激励点和目标点振动信号。接着根据需要,测取整车系统水平的频响函数。第3步按式(3)方法进行动刚度计算。第4步,得到了动刚度,依据式(4)即可计算载荷。最后,依照式(1)即可进行传递路径贡献量计算。
以下述方法计算悬置动刚度。一个系统分解为子结构A和子结构B,两个子结构之间通过弹性元件,比如发动机悬置,进行多点耦合,如图2所示,存在以下动力学关系
经演算,得到以下动刚度[Kc]计算式
式中:H代表频响函数;下标S代表此频响函数处于整车系统状态测量;下标c(a)、c(b)代表子结构A和B的耦合点。
如图3所示,建立整车振动系统模型,整车系统划分为2个子结构,发动机为子结构A,车身等为子结构B,两个子结构之间通过悬置连接,整车系统包括:发动机1、变速器2、发动机前悬置3、发动机后悬置4、变速器悬置5和麦克风6;激励点为3个发动机悬置被动端(车身侧),响应点为驾驶员耳旁噪声。每个悬置分析整车坐标系上的X、Y和Z三个振动方向,有3×3=9条路径将振动传递到驾驶员耳旁,产生噪声。按照以下步骤进行传递路径分析。
(1)在每个悬置的主动、被动侧各布置1个振动传感器,共布置6个传感器,选定一个行驶工况,比如匀速40km/h,汽车以这个工况行驶,测取各个激励点振动信号,以及响应点驾驶员耳旁噪声信号;
(2)保留第(1)工况测试时布置的振动传感器,以力锤激励法逐一在每个传感器附近对3个方向施加激励,记录每个振动传感器的信号,得到式(3)的频响函数矩阵,各个矩阵均为9×9规模。在整车状态下,以力锤激励法测取每个悬置被动端3个方向到驾驶员耳旁的频响函数,得一个9×1规模矩阵;
(3)以式(3)计算3个悬置的动刚度;
(4)基于第(1)步所得振动加速度,以式(4)计算载荷,得一个1×9规模矩阵:
式中,ω代表频率;
(5)基于第(2)步所得悬置到驾驶员耳旁的频响函数和第(4)步所得载荷使用式(1)计算各条传递路径贡献量。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (3)
1.一种基于逆子结构技术的传递路径分析方法,其特征在于,包括以下步骤:
步骤1、测取激励信号和响应信号;
步骤2、测取频响函数;
步骤3、计算动刚度;
步骤4、计算载荷;
步骤5、计算传递路径贡献量。
2.根据权利要求1所述的基于逆子结构技术的传递路径分析方法,其特征在于,所述步骤2包括以下步骤:
步骤21、将拟进行传递路径分析的系统从用弹性件连接的部位划分为2个子结构;
步骤22、确定激励点和响应点,一般取其中一个子结构上弹性连接部位为激励点;
步骤23、给2个子结构每一个弹性连接点确定互相垂直的X、Y和Z三个方向,且这些连接点的三个方向相同;
步骤24、在每个连接点布置加速度传感器,逐一在这些连接点用力锤从X、Y和Z三个方向施力,记录所有传感器的信号;
步骤25、从步骤22确定的激励点处施力,测取从此激励点到响应点的频响函数。
3.根据权利要求1所述的基于逆子结构技术的传递路径分析方法,其特征在于,所述步骤3包括以下步骤:
步骤31、抽取步骤24所述力信号和加速度响应信号;
步骤32、由力信号和加速度信号计算频响函数;
步骤33、由频响函数计算动刚度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510553126.6A CN105092194A (zh) | 2015-08-31 | 2015-08-31 | 一种基于逆子结构技术的传递路径分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510553126.6A CN105092194A (zh) | 2015-08-31 | 2015-08-31 | 一种基于逆子结构技术的传递路径分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105092194A true CN105092194A (zh) | 2015-11-25 |
Family
ID=54573125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510553126.6A Pending CN105092194A (zh) | 2015-08-31 | 2015-08-31 | 一种基于逆子结构技术的传递路径分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105092194A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106644512A (zh) * | 2016-12-22 | 2017-05-10 | 北京汽车股份有限公司 | 基于动力总成载荷的噪声分析方法和系统 |
CN106996875A (zh) * | 2017-06-02 | 2017-08-01 | 广西科技大学 | 一种车辆驾驶室传递函数获取方法 |
CN107389288A (zh) * | 2017-08-18 | 2017-11-24 | 常州市计量测试技术研究所 | 一种工业机器人的振动路径分析系统及方法 |
CN108846147A (zh) * | 2018-04-16 | 2018-11-20 | 清华大学苏州汽车研究院(相城) | 一种车辆工作时各振动激励系统贡献量计算分析方法 |
CN109211503A (zh) * | 2018-09-13 | 2019-01-15 | 中车唐山机车车辆有限公司 | 测试装置 |
CN110175392A (zh) * | 2019-05-22 | 2019-08-27 | 中国民航大学 | 基于otpa法和实物建模的航空发动机振动传递路径分析方法 |
CN110210179A (zh) * | 2019-06-28 | 2019-09-06 | 合肥工业大学 | 一种弹性元件动刚度及其两端载荷的获取方法 |
CN110220697A (zh) * | 2019-04-25 | 2019-09-10 | 重庆青山工业有限责任公司 | 一种汽车变速箱传递路径贡献量测试系统 |
CN110749406A (zh) * | 2018-07-24 | 2020-02-04 | 上汽通用五菱汽车股份有限公司 | 一种车身振动传递路径分析方法 |
CN115096582A (zh) * | 2022-07-08 | 2022-09-23 | 上海交通大学 | 一种基于原位测量的传递路径分析与齿轮故障溯源方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101271022A (zh) * | 2008-05-15 | 2008-09-24 | 上海交通大学 | 车辆系统结构振动和噪声的传递路径检测系统 |
DE102008023094A1 (de) * | 2007-05-21 | 2008-11-27 | Avl List Gmbh | Verfahren zur Berechnung von richtungskorrigierten Übertragungsfunktionen und/oder richtungskorrigierten Impedanzgrößen in einer Transferpfadanalyse einer schwingenden Struktur |
CN102519580A (zh) * | 2011-12-20 | 2012-06-27 | 重庆长安汽车股份有限公司 | 一种基于传递路径分离变速器敲击结构噪声和空气噪声的方法 |
CN103292894A (zh) * | 2013-05-17 | 2013-09-11 | 同济大学 | 一种结构振动传递函数的获取方法 |
CN103439119A (zh) * | 2013-09-04 | 2013-12-11 | 暨南大学 | 一种车辆产品机械装配的动态故障诊断方法 |
CN104215454A (zh) * | 2014-09-19 | 2014-12-17 | 武汉理工大学 | 一种汽车后桥主减速器噪音性能测试方法 |
-
2015
- 2015-08-31 CN CN201510553126.6A patent/CN105092194A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102008023094A1 (de) * | 2007-05-21 | 2008-11-27 | Avl List Gmbh | Verfahren zur Berechnung von richtungskorrigierten Übertragungsfunktionen und/oder richtungskorrigierten Impedanzgrößen in einer Transferpfadanalyse einer schwingenden Struktur |
CN101271022A (zh) * | 2008-05-15 | 2008-09-24 | 上海交通大学 | 车辆系统结构振动和噪声的传递路径检测系统 |
CN102519580A (zh) * | 2011-12-20 | 2012-06-27 | 重庆长安汽车股份有限公司 | 一种基于传递路径分离变速器敲击结构噪声和空气噪声的方法 |
CN103292894A (zh) * | 2013-05-17 | 2013-09-11 | 同济大学 | 一种结构振动传递函数的获取方法 |
CN103439119A (zh) * | 2013-09-04 | 2013-12-11 | 暨南大学 | 一种车辆产品机械装配的动态故障诊断方法 |
CN104215454A (zh) * | 2014-09-19 | 2014-12-17 | 武汉理工大学 | 一种汽车后桥主减速器噪音性能测试方法 |
Non-Patent Citations (2)
Title |
---|
莫愁等: "扩展工况传递路径分析方法改进", 《振动与冲击》 * |
莫愁等: "逆子结构传递路径分析方法", 《吉林大学学报(工学版)》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106644512A (zh) * | 2016-12-22 | 2017-05-10 | 北京汽车股份有限公司 | 基于动力总成载荷的噪声分析方法和系统 |
CN106644512B (zh) * | 2016-12-22 | 2019-04-30 | 北京汽车股份有限公司 | 基于动力总成载荷的噪声分析方法和系统 |
CN106996875A (zh) * | 2017-06-02 | 2017-08-01 | 广西科技大学 | 一种车辆驾驶室传递函数获取方法 |
CN107389288A (zh) * | 2017-08-18 | 2017-11-24 | 常州市计量测试技术研究所 | 一种工业机器人的振动路径分析系统及方法 |
CN108846147A (zh) * | 2018-04-16 | 2018-11-20 | 清华大学苏州汽车研究院(相城) | 一种车辆工作时各振动激励系统贡献量计算分析方法 |
CN108846147B (zh) * | 2018-04-16 | 2022-07-22 | 清华大学苏州汽车研究院(相城) | 一种车辆工作时各振动激励系统贡献量计算分析方法 |
CN110749406A (zh) * | 2018-07-24 | 2020-02-04 | 上汽通用五菱汽车股份有限公司 | 一种车身振动传递路径分析方法 |
CN109211503A (zh) * | 2018-09-13 | 2019-01-15 | 中车唐山机车车辆有限公司 | 测试装置 |
CN110220697A (zh) * | 2019-04-25 | 2019-09-10 | 重庆青山工业有限责任公司 | 一种汽车变速箱传递路径贡献量测试系统 |
CN110220697B (zh) * | 2019-04-25 | 2023-10-24 | 重庆青山工业有限责任公司 | 一种汽车变速箱传递路径贡献量测试系统 |
CN110175392A (zh) * | 2019-05-22 | 2019-08-27 | 中国民航大学 | 基于otpa法和实物建模的航空发动机振动传递路径分析方法 |
CN110175392B (zh) * | 2019-05-22 | 2022-12-13 | 中国民航大学 | 基于otpa法和实物建模的航空发动机振动传递路径分析方法 |
CN110210179A (zh) * | 2019-06-28 | 2019-09-06 | 合肥工业大学 | 一种弹性元件动刚度及其两端载荷的获取方法 |
CN110210179B (zh) * | 2019-06-28 | 2023-03-24 | 合肥工业大学 | 一种弹性元件动刚度及其两端载荷的获取方法 |
CN115096582A (zh) * | 2022-07-08 | 2022-09-23 | 上海交通大学 | 一种基于原位测量的传递路径分析与齿轮故障溯源方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105092194A (zh) | 一种基于逆子结构技术的传递路径分析方法 | |
CN101393079B (zh) | 汽车车身结构疲劳寿命预测系统 | |
US10753827B2 (en) | Method and system for measuring vertical wheel impact force in real-time based on tire pressure monitoring | |
CN105222886A (zh) | 一种改进的提高精度的工况传递路径分析的方法 | |
CN104006979A (zh) | 转向架悬挂系统参数测试试验装置及测试方法 | |
CN110298125A (zh) | 一种基于虚拟试验场技术的疲劳分析方法 | |
CN103761363A (zh) | 乘用车副车架强度和疲劳分析方法 | |
CN106033027A (zh) | 重型汽车动力总成悬置系统的刚体模态测试装置及其方法 | |
CN110119577A (zh) | 一种提取纯电动卡车电池包支架结构道路载荷谱的方法 | |
CN106644512A (zh) | 基于动力总成载荷的噪声分析方法和系统 | |
CN103900826A (zh) | 实时监测汽车底盘结构疲劳损伤的方法 | |
CN112595528B (zh) | 车辆行驶工况下动力总成悬置主动侧激励力的识别方法 | |
CN104296924A (zh) | 动力总成惯性参数测试方法 | |
CN106769093A (zh) | 基于车轮六分力载荷的路噪预测分析方法及系统 | |
CN108170643B (zh) | 一种免拆分的子结构频响函数识别方法 | |
CN110210179A (zh) | 一种弹性元件动刚度及其两端载荷的获取方法 | |
CN115270296A (zh) | 一种商用车驾驶室疲劳耐久性分析方法及系统 | |
JPH1151820A (ja) | 車輌のスピンドルに連結されたシミュレータで使用するタイヤをモデリングする方法および装置 | |
CN105651473A (zh) | 用于自动确定物体的动态刚度的方法 | |
CN106596008B (zh) | 一种模拟实际驾驶工况的车载燃料电池振动试验系统 | |
CN102680214B (zh) | 波纹管隔振性能试验装置与试验方法 | |
CN108228994A (zh) | 在越野道路随机激励下车辆及工作装置的应力计算方法 | |
CN109426691A (zh) | 车辆运动部件仿真方法及仿真系统 | |
Forsén | Parameter identification in heavy vehicle simulation | |
CN113656994B (zh) | 一种汽车悬置系统的悬置力获取方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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: 20151125 |