CN108170643B - 一种免拆分的子结构频响函数识别方法 - Google Patents

一种免拆分的子结构频响函数识别方法 Download PDF

Info

Publication number
CN108170643B
CN108170643B CN201711426851.2A CN201711426851A CN108170643B CN 108170643 B CN108170643 B CN 108170643B CN 201711426851 A CN201711426851 A CN 201711426851A CN 108170643 B CN108170643 B CN 108170643B
Authority
CN
China
Prior art keywords
frequency response
response function
mechanical system
passive
point
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.)
Expired - Fee Related
Application number
CN201711426851.2A
Other languages
English (en)
Other versions
CN108170643A (zh
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.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong University
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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201711426851.2A priority Critical patent/CN108170643B/zh
Publication of CN108170643A publication Critical patent/CN108170643A/zh
Application granted granted Critical
Publication of CN108170643B publication Critical patent/CN108170643B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • 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
    • G06F17/17Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method

Landscapes

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

Abstract

一种基于原位测量频响函数且免拆分的子结构频响函数识别方法,首先将机械系统中的部件分类为主动件、被动件和弹性元件并建立耦合频响函数矩阵,然后基于被动件频响函数预测公式得到解耦机械系统的频响函数。本发明基于“源‑路径‑接受体”模型,采用耦合机械频响函数预测解耦机械频响函数,克服现有技术无法计算机械子结构任意自由度之间频响函数的问题,提高频响函数的计算精度,并为机械系统振动传递路径分析奠定了基础。

Description

一种免拆分的子结构频响函数识别方法
技术领域
本发明涉及的是一种机械振动分析检测领域的技术,具体是一种基于原位测量频响函数且免拆分的子结构频响函数识别方法。
背景技术
传递路径分析方法(TPA)广泛用于分析和处理复杂机械系统的振动与噪声问题,通过TPA能够识别和量化激励源,分析能量从激励源传递至目标点的路径,准确评估和排序不同传递路径对目标点的贡献量,通过控制和改进这些路径可以使噪声和振动控制在预定的目标值内。经典TPA因考虑信息全面、分析精度高等优点已成为汽车NVH领域中的标准TPA。经典TPA主要包含两个步骤:工况力的识别和被动部件频响函数的测量。其中被动部件频响函数的测量需要拆掉主动部件,这导致经典TPA消耗大量的人力和时间。后续发展的TPA虽然可以提高分析效率,但是以牺牲精度为代价。如何既能实施经典TPA分析,又能避免拆分系统一直是TPA的一个研究难点。
发明内容
本发明针对现有技术存在的上述不足,提出一种基于原位测量频响免拆分的子结构频响函数识别方法,基于“源-路径-接受体”模型,采用耦合机械频响函数预测解耦机械频响函数,克服现有技术无法计算机械子结构任意自由度之间频响函数的问题,提高频响函数的计算精度,并为机械系统振动传递路径分析奠定了基础。
本发明是通过以下技术方案实现的:
本发明首先将机械系统中的部件分类为主动件、被动件和弹性元件并建立耦合机械系统的频响函数矩阵,然后基于被动件频响函数预测公式得到解耦机械系统的频响函数。
本发明具体包含以下步骤:
步骤一:机械系统划分,具体为:根据机械系统的结构特点和受力情况,将直接受到激励作用的部件定为主动件,非直接受到激励作用的定为被动件,主动件和被动件之间通过弹性元件连接。
步骤二:测量耦合机械系统的频响函数矩阵,具体为:根据关注的响应位置和振动噪声频率范围,确定被动件上的目标点及分析频率范围,以弹性元件两端点和目标点为测量点,测量弹性元件的主动件一侧端点、弹性元件的被动件一侧端点及目标位置点之间的频响函数,组成耦合机械系统的频响函数矩阵。
步骤三:解耦机械系统频响函数预测,具体为:以测量得到的耦合机械系统频响函数为输入,根据子结构频响函数预测公式得到解耦机械系统的频响函数。
所述的耦合机械系统是指包含主动件和被动件的机械系统。
所述的解耦机械系统是指去除主动件后,只包含被动件的机械系统。
所述的解耦机械系统的频响函数是指被动件的频响函数。
所述的子结构频响函数预测公式为:Hd,tp=Hc,tp-Hc,ta(Hc,aa-Hc,pa)-1(Hc,ap-Hc,pp),其中:Hc,aa为弹性元件的主动件一侧端点频响函数矩阵,Hc,pp为弹性元件的被动件一侧端点频响函数矩阵,Hc,ap为弹性元件的被动件一侧端点到主动件一侧端点的频响函数矩阵,Hc,pa为弹性元件的主动件一侧端点到被动件一侧端点的频响函数矩阵,Hc,tp为弹性元件的被动件一侧端点到目标点的频响函数矩阵,Hc,ta为弹性元件的主动件一侧端点到目标点的频响函数矩阵,Hd,tp为待求的弹性元件的被动件一侧端点到目标点的解耦频响函数矩阵。
技术效果
与现有技术相比,本发明的技术效果包括:
1)本发明基于子结构频响函数预测公式,避开了对解耦机械系统频响函数的直接实验测量,实现了解耦机械系统任意自由度之间频响函数的预测,提高了机械系统振动传递路径分析效率。
2)本发明所提的子结构频响函数预测公式,包含矩阵求逆操作少,受噪声等因素影响较小,提高了解耦频响函数的预测精度。
3)本发明的工程应用性强。利用本发明计算出的解耦机械系统频响函数,结合系统工况数据,可以很方便计算出悬置安装点的工况耦合力,得到振动传递路径贡献度,大大提高了传递路径分析方法在机械振动问题中的可操作性和分析效率,虽然本发明是针对解耦频响函数提出的,但本发明的理论成果具有很强的应用前景,稍加扩展便可应用于悬置动刚度的识别与优化。
附图说明
图1为本发明方法流程图;
图2为离散机械系统图;
图3为被动件解耦频响函数预测值和真实值的比较图;
图4为简易车身实验测试图;
图5为悬置示意图;
图6为简易车身解耦频响函数预测值和真实值的比较图。
具体实施方式
实施例1
如图1所示,本实施例包括如下步骤:
步骤一、以图2所示的离散机械系统,其包含主动件和被动件,主动件由4个质量块M5、M6、M8、M9组成,被动件由5个质量块M1~M4、M7组成,主动件与被动件之间有3条传递路径K35、K46、K78,质量块M3、M4和M7为被动件一侧耦合点,而质量块M5、M6、M8为主动件一侧耦合点,质量块M1的位移为目标响应,一共有7个测量点,分析频率范围为1-250Hz。
步骤二、测量耦合机械系统频响函数:测量耦合点主动件一侧频响函数矩阵Hc,aa、耦合点被动件一侧频响函数矩阵Hc,pp、耦合点被动件一侧到主动件一侧频响函数矩阵Hc,ap、耦合点主动件一侧到被动件一侧频响函数矩阵Hc,pa,耦合点被动件一侧到目标点频响函数矩阵Hc,tp,耦合点主动件一侧到目标点频响函数矩阵Hc,ta
以上耦合机械系统的频响函数矩阵展开为:
Figure GDA0002815498650000031
Figure GDA0002815498650000032
Hc,tp=[Hc,13 Hc,14 Hc,17],Hc,ta=[Hc,15 Hc,16 Hc,18],其中:Hc,ij为施力点i到响应点j的耦合频响函数。
步骤三、解耦机械系统频响函数预测:以步骤二中测量的耦合机械系统的频响函数矩阵作为输入,用被动件频响函数预测公式计算解耦频响函数。计算结果如图3所示,由图可知,预测得到的3个解耦频响函数与它们的真实值分别相同,验证了本方法的理论正确性。
实施例2
如图4所示,本实施例展示了简易车身物理模型,该模型包含车身(被动件,如图4(a)所示)和“发动机”支架(主动件,如图4(c)所示),支架通过3个橡胶悬置与车身连接,实验装置如图4(b)所示。
本实施例包括如下步骤:
步骤一、确定测量点和分析频率段:将主动件一侧悬置点记为a1、a2和a3,被动件一侧悬置点记为p1、p2和p3,悬置示意图及全局坐标系如图5所示。目标点为车身右侧一点(单向)振动响应,记为t。只考虑悬置点的平动自由度,因此共有19个自由度需要关注,分析频率范围为20-400Hz。
步骤二、测量耦合机械系统频响函数:在悬置点用3D加速度传感器拾取加速度信号,在目标点用1D加速度传感器拾取加速度。在同一激励点连续锤击5次,通过LMS系统采集、分析力和加速度信号,取5次测量的平均值,得到耦合系统频响函数矩阵。
为了方便比较,将“发动机”支架去掉,进行车身解耦频响函数测试,如图4(d)所示。
步骤三、解耦机械系统频响函数预测:以步骤二中测得的耦合系统频响函数作为输入,用被动件频响函数预测公式预测被动件解耦频响函数。简易车身解耦频响函数的预测值和测量值比较结果如图6所示,图6(a)为前悬置下点X方向到目标点的解耦频响函数,图6(b)为前悬置下点Y方向到目标点的解耦频响函数,图6(c)为前悬置下点Z方向到目标点的解耦频响函数。由图可知,简易车身解耦频响函数的预测值与测量值在整体上比较吻合。
上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式对其进行局部调整,本发明的保护范围以权利要求书为准且不由上述具体实施所限,在其范围内的各个实现方案均受本发明之约束。

Claims (1)

1.一种基于原位测量频响函数且免拆分的子结构频响函数识别方法,其特征在于,首先将机械系统中的部件分类为主动件、被动件和弹性元件并建立耦合机械系统的频响函数矩阵,然后基于被动件频响函数预测公式得到解耦机械系统的频响函数;
所述的分类是指:根据机械系统的结构特点和受力情况,将直接受到激励作用的部件定为主动件,非直接受到激励作用的定为被动件,主动件和被动件之间通过弹性元件连接;
所述的解耦机械系统的频响函数,以耦合机械系统的频响函数矩阵作为输入,根据被动件频响函数预测公式预测解耦机械系统的频响函数,即去除主动件后的机械系统的频响函数,其中被动件频响函数预测公式为:
Figure 755119DEST_PATH_IMAGE002
,其中:Hc,aa为弹性元件的主动件一侧端点频响函数矩阵,Hc,pp为弹性元件的被动件一侧端点频响函数矩阵,Hc,ap为弹性元件的被动件一侧端点到主动件一侧端点的频响函数矩阵,Hc,pa为弹性元件的主动件一侧端点到被动件一侧端点的频响函数矩阵,Hc,tp为弹性元件的被动件一侧端点到目标点的频响函数矩阵,Hc,ta为弹性元件的主动件一侧端点到目标点的频响函数矩阵,Hd,tp为待求的弹性元件的被动件一侧端点到目标点的解耦频响函数矩阵;
所述的耦合机械系统的频响函数矩阵,根据关注的响应位置和振动噪声频率范围,确定被动件上的目标点及分析频率范围,以弹性元件两端点和目标点为测量点,测量悬置的主动件一侧安装点、悬置的被动件一侧安装点及目标位置点之间的频响函数,组成耦合机械系统的频响函数矩阵;
所述的耦合机械系统是指包含主动件和被动件的机械系统。
CN201711426851.2A 2017-12-26 2017-12-26 一种免拆分的子结构频响函数识别方法 Expired - Fee Related CN108170643B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711426851.2A CN108170643B (zh) 2017-12-26 2017-12-26 一种免拆分的子结构频响函数识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711426851.2A CN108170643B (zh) 2017-12-26 2017-12-26 一种免拆分的子结构频响函数识别方法

Publications (2)

Publication Number Publication Date
CN108170643A CN108170643A (zh) 2018-06-15
CN108170643B true CN108170643B (zh) 2021-02-09

Family

ID=62520920

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711426851.2A Expired - Fee Related CN108170643B (zh) 2017-12-26 2017-12-26 一种免拆分的子结构频响函数识别方法

Country Status (1)

Country Link
CN (1) CN108170643B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110749406B (zh) * 2018-07-24 2021-06-18 上汽通用五菱汽车股份有限公司 一种车身振动传递路径分析方法
CN110210179B (zh) * 2019-06-28 2023-03-24 合肥工业大学 一种弹性元件动刚度及其两端载荷的获取方法
CN114778047B (zh) * 2022-06-17 2022-09-02 中国飞机强度研究所 一种飞机减振构件的动刚度和一致性测试装置
CN114778048B (zh) * 2022-06-17 2022-09-02 中国飞机强度研究所 一种飞机减振构件的动刚度和一致性测试方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102689229A (zh) * 2012-05-04 2012-09-26 华中科技大学 基于响应耦合的刀尖点频响函数获取方法
CN102880803A (zh) * 2012-09-29 2013-01-16 西安交通大学 一种复杂机械结构的转动自由度频率响应函数计算方法
CN103308333A (zh) * 2013-05-17 2013-09-18 同济大学 一种隔振器动刚度测试方法
KR101406778B1 (ko) * 2013-04-02 2014-06-17 (주)대주기계 가상 주파수 응답함수 기반의 기계이상 진단방법 및 이를 이용한 진단 시스템
CN104180896A (zh) * 2014-09-11 2014-12-03 江南大学 一种间接测定极脆部件非耦合界面频率响应函数的方法
CN104880248A (zh) * 2015-05-07 2015-09-02 中国船舶重工集团公司第七一二研究所 一种电机结构噪声激励源的贡献量定量识别方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102689229A (zh) * 2012-05-04 2012-09-26 华中科技大学 基于响应耦合的刀尖点频响函数获取方法
CN102880803A (zh) * 2012-09-29 2013-01-16 西安交通大学 一种复杂机械结构的转动自由度频率响应函数计算方法
KR101406778B1 (ko) * 2013-04-02 2014-06-17 (주)대주기계 가상 주파수 응답함수 기반의 기계이상 진단방법 및 이를 이용한 진단 시스템
CN103308333A (zh) * 2013-05-17 2013-09-18 同济大学 一种隔振器动刚度测试方法
CN104180896A (zh) * 2014-09-11 2014-12-03 江南大学 一种间接测定极脆部件非耦合界面频率响应函数的方法
CN104880248A (zh) * 2015-05-07 2015-09-02 中国船舶重工集团公司第七一二研究所 一种电机结构噪声激励源的贡献量定量识别方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Response prediction for modified mechanical systems based on in-situ frequency response functions: Theoretical and numerical studies;Zengwei Wang et al.;《Journal of Sound and Vibration》;20170424;第417-441页 *
Structural modifications using frequency response functions;H. NEVZAT;《Mechanical Systems and Signal Processing》;19901231;第53-63页 *
基于RCSA的深孔内圆磨床主轴端点频响函数预测;李孝茹 等;《中国机械工程》;20151031;第26卷(第19期);第2652-2661页 *
频率响应函数估计方法综述;段虎明 等;《振动与冲击》;20080531;第48-51页 *

Also Published As

Publication number Publication date
CN108170643A (zh) 2018-06-15

Similar Documents

Publication Publication Date Title
CN108170643B (zh) 一种免拆分的子结构频响函数识别方法
Janssens et al. OPAX: A new transfer path analysis method based on parametric load models
WO2016197552A1 (zh) 基于模型识别与等效简化的高速平台运动参数自整定方法
EP2113755B1 (en) Vibrational and/or acoustic transfer path analysis
CN112595528B (zh) 车辆行驶工况下动力总成悬置主动侧激励力的识别方法
CN105222886A (zh) 一种改进的提高精度的工况传递路径分析的方法
JP2008134182A5 (zh)
CN106033027A (zh) 重型汽车动力总成悬置系统的刚体模态测试装置及其方法
CN112697448A (zh) 车辆怠速工况下动力总成悬置主动侧激励力的识别方法
EP3397935B1 (en) Vibration and noise mapping system and method
CN106644512A (zh) 基于动力总成载荷的噪声分析方法和系统
CN110210179B (zh) 一种弹性元件动刚度及其两端载荷的获取方法
JP5493373B2 (ja) 伝達経路毎の成分を算定するための方法
CN106908205A (zh) 一种发动机悬置动刚度的测试装置
CN108593092A (zh) 一种基于传递路径分析的高端纺织机械噪声识别方法
EP2574951A2 (en) A method for analyzing sound transmission paths in a system
CN102270249B (zh) 识别零部件特征频率的方法
CN114491825B (zh) 一种汽车翼子板振动强度分析方法
Kim et al. Local stiffness control for reducing vehicle interior noise by using FRF-based synthesis method
JP2009059094A (ja) 振動解析システム及び振動解析方法
Janssens et al. A novel transfer path analysis method delivering a fast and accurate noise contribution assessment
Wagner et al. High frequency source characterization of an e-motor using component-based TPA
KR101449331B1 (ko) 엔진 주변부품의 가속 진동 수명시험방법 및 그 장치
CN113297907A (zh) 一种脉冲激励下基于数据驱动的非线性阻尼识别方法
Wu et al. Comparison of classical and operational transfer path analysis in vehicle NVH improvement

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210209

Termination date: 20211226