CN108267221B - 一种消除串扰的振动传递路径测试方法 - Google Patents
一种消除串扰的振动传递路径测试方法 Download PDFInfo
- Publication number
- CN108267221B CN108267221B CN201710002316.8A CN201710002316A CN108267221B CN 108267221 B CN108267221 B CN 108267221B CN 201710002316 A CN201710002316 A CN 201710002316A CN 108267221 B CN108267221 B CN 108267221B
- Authority
- CN
- China
- Prior art keywords
- vibration source
- vibration
- transmission function
- vector
- broad sense
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种消除串扰的振动传递路径测试方法,包括:采用测量传递函数或非正交向量的形式获得解除串扰传递矩阵,并计算解串扰广义参考信号或受挡力;基于广义参考信号或受挡力对工况传递路径进行分析,获得传递系数;根据广义参考信号或受挡力,以及传递系数,计算各个振源在一响应位置处的贡献量。
Description
技术领域
本发明涉及一种振动源的测量识别方法,特别涉及一种消除串扰的振动传递路径测试方法。
背景技术
复杂结构的振动包括三部分:振动源、振动传递路径及振动接受结构,若采用一定技术手段发现主要振源和传递路径,针对主要振源和传递路径采取减振降噪办法将能有效降低振动与噪声,在非关键部位采用的减振降噪手段起不到主要作用。减振降噪的前提也是传递路径分析和振源大小的辨识问题。
传递路径分析需要解除串扰激励源,即拆卸发动机等振源获得传递函数,很多情况不易于实现。在参考文献1“A.T.Moorhouse,A.S.Elliott,T.A.Evans.In situmeasurement of the blocked force of structure-borne sound sources.Journal ofSound and Vibration,2009(325):679–685”中提出了一种基于受挡力传递路径分析法,其中的受挡力是指将机器拆离固定于刚性基础上施加给基础的力。如:一振源结构如图1所示,受到力f1作用下,自由状态2点响应为a2;若将2点安装在固定地面上,使该点响应a2=0,这时作用在地面2点的力为f2 B,称该力为受挡力。根据传递函数法,可以列出1、2点间的方程为:
假设a2=0,推导出因此受挡力完全与力源和力源结构有关,受挡力代表了振源结构的固有特性。
对于实际结构,可通过不拆卸机器获得受挡力,如:一个受到三个源的复合系统,在不拆离机器的情况下,依次敲击机器隔振器根部位置,例如:在图2中的步骤(1)中,敲击A隔振器根部,测量A,B,C根部响应然后计算传递函数式中上标表示实验步骤,下标表示位置。同理步骤(2),敲击B隔振器根部获得传递函数步骤(3)时,敲击C隔振器根部获得传递函数
因此,获取上述传递函数后可以通过下式计算获得受挡力:
三个机器对任何一个响应点0的传递函数为:则各振源对该点贡献为:
该方法的缺点是需要在各振源处激励获得复合结构的传递函数,但实际源处空间有限,无法布置激振器,应用存在困难;另外也无法反映运行工况下的结构特性。
另外一种传递路径分析法是OTPA(operational transfer path analysis)传递路径测试法,参见参考文献2“D.deKlerk,A.Ossipov.Operational transfer pathanalysis:Theory,guidelines and tire noise application.Mechanical Systems andSignal Processing,2010(24):1950–1962”,其原理是在激励源处布置传感器的测量响应信号作为参考信号(如下方程的a1,a2,...am),将关心响应点处的响应作为观察点(如下方程中的ay),通过多工况运行获得以下方程:
则可解出传递率系数:
一般效工况数s应大于或等于m,上式右边的求逆运算唯一,可获得传递率系数矩阵,计算某个工况下源的贡献大小为:
该方法不需要拆卸机器,可在线工况测量,测试结果反映实际工况情况。但表征振源的参考响应信号存在串扰,如图3所示:用来代替振源的振动响应不仅仅来源于该源,还来源于其它源,串扰造成源分析的错误。
综上所述,现有的传递路径分析方法各自存在一定的缺陷,需要一种易于应用、且能够有效消除串扰的路径测试方法。
发明内容
本发明的目的在于克服现有的传递路径分析方法所存在的缺陷,提供一种易于应用、且能够有效消除串扰的路径测试方法。
为了实现上述目的,本发明提供了一种消除串扰的振动传递路径测试方法,包括:
步骤1)、采用测量传递函数或非正交向量的形式获得解除串扰传递矩阵,并计算解串扰广义参考信号或受挡力;
步骤2)、基于步骤1)得到的广义参考信号或受挡力对工况传递路径进行分析,获得传递系数;
步骤3)、根据步骤1)得到的广义参考信号或受挡力,以及步骤2)得到的传递系数,计算各个振源在一响应位置处的贡献量。
上述技术方案中,在步骤1)中,在测量传递函数或非正交向量时分为三种情况,包括:能够分步运转的情况;无法分步运转的情况;部分能够分步运转、部分无法分步运转的情况;其中,
对于能够分步运转的情况,在测量时,在各振源处布置参考点响应传感器,同时分步运转各独立激励源;假设共n个振源,m个参考点响应传感器,各点的测量响应构成坐标向量其中上标表示第i个振源独立运行的工况,下标表示参考点响应传感器编号,即表示第i个振源独立作用时在j位置的响应频谱;
参考点响应传感器个数应多于振源个数,即m≥n,且尽量靠近振源处;n个振源单步运转获得n个独立参考信号向量;利用这些向量构建解耦矩阵,获得广义参考向量,由于广义参考向量与受挡力在本申请方法中的作用没有本质的区别,所以后面对广义参考向量和受挡力都用符号表示;广义参考向量的计算过程为:
其中:上标rpm(j)表示第j个转速工况,“+”表示求广义逆;
对于无法分步运转的情况,在测量时,通过力锤激励的方式获得传递函数,进而得到解耦受挡力;具体包括:
采用力锤在激励源处一次激励起结构对所有参考点的传递函数,其传递函数形式是通常参考点的个数m大于振源个数,对于n个振源各独立激励获得所有传递函数:
采用力锤激励的方式获得传递函数后,通过下面的公式(2)获得解耦受挡力:
对于部分能够分步运转、部分无法分步运转的情况,假设前k个能够分步运转,后k+1到n个振源无法分步运转,则构建以下非正交分解向量矩阵,以及获取广义参考向量的计算如公式(3):
上述技术方案中,对于无法分步运转的情况,若无法在振源处激励,则在参考点激励,并利用互易原理获得传递函数;具体包括:
若在一参考点j处激励,能够一次性获得j点对所有参考点的传递函数:
通过对m个参考点激励,通常m≥n,获得所有传递函数其中任何一行是一非正交传递向量;
采用互易法获得传递函数后,同样通过下面的公式(4)获得解耦受挡力:
上述技术方案中,所述步骤2)包括:构建基于受挡力的工况传递路径分析法,如下:
其中表示第i个工况的响应信号,通过求广义逆获得传递系数:
其中:s为总工况数,通常s≥n。
上述技术方案中,所述步骤3)采用如下公式计算各个振源的贡献量:
本发明的优点在于:
本发明的方法精度很高,且易于实现。
附图说明
图1是振源结构的示意图;
图2是现有技术中梁上含三振源受挡力测量示意图;
图3是串扰对工况传递路径分析影响示意图;
图4是互易示意图;
图5是采用本发明方法的实验装置示意图;
图6是采用现有技术中的OTPA方法的分析结果和真实结果的对比图;
图7是采用本发明的方法的分析结果和真实结果的对比图。
具体实施方式
现结合附图对本发明作进一步的描述。
本发明的消除串扰的振动传递路径测试方法利用互易原理及增加虚拟振源点分步运转法解决了传递函数测试的困难,利用解除串扰矩阵求振源的广义参考向量或受挡力,利用广义参考向量或受挡力参与工况传递路径分析。
本发明的方法具体包括以下步骤:
步骤1.采用测量传递函数或非正交向量的形式获得解除串扰传递矩阵,并计算解串扰广义参考信号或受挡力。
具体的说,可进一步分为三种情况,一是可以分步运转的情况,二是无法分步运转的情况,三是部分可以分步运转、部分无法分步运转的情况。下面分别予以说明。
步骤1.1对可以分步运转的情况,一般可在各振源处布置参考点响应传感器,同时分步运转各独立激励源。
假设共n个振源,m个参考点响应传感器,各点的测量响应构成坐标向量其中上标表示第i个振源独立运行的工况,下标表示参考点响应传感器编号,即表示第i个源独立作用时在j位置的响应频谱,该值为某个频率点下的复数频谱值,分析时通常一个一个频率点进行。
参考点响应传感器个数应多于振源个数,即m≥n,且尽量靠近振源处。n个振源单步运转可以获得n个独立参考信号向量,这些向量之间一般不正交,但是反映了各振源独立运行时的状态。利用这些向量可以构建解耦矩阵,获得广义参考向量,由于广义参考向量与受挡力在本申请方法中的作用没有本质的区别,所以后面对广义参考向量和受挡力都用符号表示,广义参考向量的计算过程为:
其中:上标rpm(j)表示第j个转速工况,“+”表示求广义逆。
其中参考点的选择对广义参考信号的影响比较大,可以通过优化选择测点降低矩阵的条件数,从而减小求解过程对测试误差的敏感程度。
步骤1.2对于无法分步运转的情况,可以通过力锤激励的方式获得传递函数,进而得到解耦受挡力。
具体的说,采用力锤可在激励源处一次激励起结构对所有参考点的传递函数,其传递函数形式是通常参考点的个数m可大于振源个数,对于n个振源各独立激励可以获得所有传递函数:
作为一种优选实现方式,当有些结构空间比较狭窄,无法布放激振器时,可采用互易法获得非正交传递向量形式的传递函数,进而再得到解耦受挡力。
具体的说,对无法在振源处激励的情形,可在参考点激励,利用互易原理获得以上传递函数,力学互易原理如图4所示,对一结构板,为了求2点对1点的传递函数,通常需要在1位置施加力同时测量2位置响应计算得到传递函数但是由于1位置没有足够空间布置激振器施加激励,则可以在2位置施加激励测量1位置的响应同时利用互易原理获得传递函数:
所以在参考点j激励,可一次性获得j点对所有参考点的传递函数:通过对m个参考点激励,通常m≥n,可获得所有传递函数其中任何一行是一非正交传递向量。
采用力锤激励的方式获得传递函数或采用互易法获得传递函数后,通过下面的公式(2)获得解耦受挡力:
步骤1.3对部分可以分步运转的情况,部分无法分步运转的情况,应该尽量采用分步运转获得的非正交向量。分步运转向量获取方式参考之前的步骤1.1,传递函数向量获取方式参考之前的步骤1.2。
假设前k个源可分步运转,后k+1到n个源无法分步运转,则可以构建以下非正交分解向量矩阵,以及获取广义参考向量的计算如公式(3)。
步骤2.基于步骤1所得到的广义参考向量或受挡力,对工况传递路径进行分析,获得传递系数。
构建基于受挡力的工况传递路径分析法,如下:
其中表示第i个工况的响应信号,通过求广义逆获得传递系数:
其中:s为总工况数,通常s≥n。
3.计算各个振源在一响应位置处的贡献量。
在分析响应位置j处的工况时,可以通过以下过程计算各振源贡献量。
3.1计算参考信号
根据响应位置j处的实际情况,采用公式(1)或(2)或(3)计算广义参考向量或受挡力。
3.2由参考信号以及传递系数,计算各个振源的贡献量
其中加法中的每个项代表各振源的贡献量。
参与工况传递路径分析的量是广义参考向量或受挡力,即实测响应通过特定运算获得的参考量,采用广义参考信号和受挡力构建多工况方程获得传递率系数,同时对某被分析工况求得广义参考信号,将其值与传递率系数做乘法获得各源贡献,消除了串扰影响,精度更高。
为了证明本申请方法的实际效果,申请人在一个试验中设计了传递路径实验台架,并进行了相关试验,实验装置如图5所示。
在该试验中,采用两个力源激励结构,分别采用现有技术中的工况传递路径(OTPA)方法和本申请的方法来计算源的贡献量。两个激振器可以独立控制输入信号电压大小,是两个可独立运行的源。在平板上共同布置4个加速度测量点,其中1、2号点在激励源附近,3、4号点在离开激励源的位置。采用采集设备采集四个点的振动加速度响应,先采用OTPA分析方法,然后采用基于本申请的串扰消除的方法。分析时,采用1、2、4号点为参考点,将点3号点为响应点,分别计算响应源对总响应的贡献。各工况安装状况不变,通过调节输入信号电压大小调节工况状态,各工况输入信号电压表如下:
采用工况(1)、(2)获得广义参考信号的解耦矩阵,采用工况(3~8)进行基于广义参考信号的工况传递路径分析,采用(9)进行验证实验。图6是采用现有技术中的OTPA方法分析结果和真实结果的对比图,图7是采用本申请方法分析结果和真实结果对比图。可以看出,采用本申请方法所得到的结果几乎和真实结果完全重合,而采用OTPA部分峰值频率点误差在3dB以上,因此本申请方法的精度远高于OTPA方法。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (4)
1.一种消除串扰的振动传递路径测试方法,包括:
步骤1)、采用测量传递函数或非正交向量的形式获得解除串扰传递矩阵,并计算解串扰广义参考信号或受挡力;
在步骤1)中,在测量传递函数或非正交向量时分为三种情况,包括:能够分步运转的情况;无法分步运转的情况;部分能够分步运转、部分无法分步运转的情况;其中,
对于能够分步运转的情况,在测量时,在各振源处布置参考点响应传感器,同时分步运转各独立激励源;假设共n个振源,m个参考点响应传感器,各点的测量响应构成坐标向量其中上标表示第i个振源独立运行的工况,下标表示参考点响应传感器编号,即表示第i个振源独立作用时在j位置的响应频谱;
参考点响应传感器个数应多于振源个数,即m≥n,且尽量靠近振源处;n个振源单步运转获得n个独立参考信号向量;利用这些向量构建解耦矩阵,获得广义参考向量,广义参考向量的计算过程为:
其中:上标rpm(j)表示第j个转速工况,“+”表示求广义逆;
对于无法分步运转的情况,在测量时,通过力锤激励的方式获得传递函数,进而得到解耦受挡力;具体包括:
采用力锤在激励源处一次激励起结构对所有参考点的传递函数,其传递函数形式是参考点的个数m大于振源个数,对于n个振源各独立激励获得所有传递函数:
采用力锤激励的方式获得传递函数后,通过下面的公式(2)获得解耦受挡力:
对于部分能够分步运转、部分无法分步运转的情况,假设前k个能够分步运转,后k+1到n个振源无法分步运转,则构建以下非正交分解向量矩阵,以及获取广义参考向量的计算如公式(3):
步骤2)、基于步骤1)得到的广义参考信号或受挡力对工况传递路径进行分析,获得传递系数;
步骤3)、根据步骤1)得到的广义参考信号或受挡力,以及步骤2)得到的传递系数,计算各个振源在一响应位置处的贡献量。
2.根据权利要求1所述的消除串扰的振动传递路径测试方法,其特征在于,对于无法分步运转的情况,若无法在振源处激励,则在参考点激励,并利用互易原理获得传递函数;具体包括:
若在一参考点j处激励,能够一次性获得j点对所有参考点的传递函数:
通过对m个参考点激励,m≥n,获得所有传递函数其中任何一行是一非正交传递向量;
采用互易法获得传递函数后,同样通过下面的公式(4)获得解耦受挡力:
3.根据权利要求2或1所述的消除串扰的振动传递路径测试方法,其特征在于,所述步骤2)包括:构建基于受挡力的工况传递路径分析法,如下:
其中表示第i个工况的响应信号,通过求广义逆获得传递系数:
其中:s为总工况数,s≥n。
4.根据权利要求3所述的消除串扰的振动传递路径测试方法,其特征在于,所述步骤3)采用如下公式计算各个振源的贡献量:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710002316.8A CN108267221B (zh) | 2017-01-03 | 2017-01-03 | 一种消除串扰的振动传递路径测试方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710002316.8A CN108267221B (zh) | 2017-01-03 | 2017-01-03 | 一种消除串扰的振动传递路径测试方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108267221A CN108267221A (zh) | 2018-07-10 |
CN108267221B true CN108267221B (zh) | 2019-09-24 |
Family
ID=62770646
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710002316.8A Active CN108267221B (zh) | 2017-01-03 | 2017-01-03 | 一种消除串扰的振动传递路径测试方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108267221B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110749406B (zh) * | 2018-07-24 | 2021-06-18 | 上汽通用五菱汽车股份有限公司 | 一种车身振动传递路径分析方法 |
CN111353329B (zh) * | 2018-12-20 | 2023-05-19 | 中国科学院声学研究所 | 一种传感器组合源表征的振源串扰消除方法 |
CN110175392B (zh) * | 2019-05-22 | 2022-12-13 | 中国民航大学 | 基于otpa法和实物建模的航空发动机振动传递路径分析方法 |
CN112432702B (zh) * | 2020-11-09 | 2021-10-15 | 大连理工大学 | 一种基于离心泵振动传递路径重合的振动源识别方法 |
CN115096582B (zh) * | 2022-07-08 | 2023-06-02 | 上海交通大学 | 一种基于原位测量的传递路径分析与齿轮故障溯源方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103630232A (zh) * | 2013-10-29 | 2014-03-12 | 南车青岛四方机车车辆股份有限公司 | 一种高速动车组噪声源识别测试方法 |
US8885843B2 (en) * | 2010-04-15 | 2014-11-11 | Hitachi, Ltd. | Sound generator |
CN104880248A (zh) * | 2015-05-07 | 2015-09-02 | 中国船舶重工集团公司第七一二研究所 | 一种电机结构噪声激励源的贡献量定量识别方法 |
-
2017
- 2017-01-03 CN CN201710002316.8A patent/CN108267221B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8885843B2 (en) * | 2010-04-15 | 2014-11-11 | Hitachi, Ltd. | Sound generator |
CN103630232A (zh) * | 2013-10-29 | 2014-03-12 | 南车青岛四方机车车辆股份有限公司 | 一种高速动车组噪声源识别测试方法 |
CN104880248A (zh) * | 2015-05-07 | 2015-09-02 | 中国船舶重工集团公司第七一二研究所 | 一种电机结构噪声激励源的贡献量定量识别方法 |
Non-Patent Citations (1)
Title |
---|
工况传递路径分析法原理及其应用;伍先俊 等;《噪声与振动控制》;20140228;第34卷(第1期);第28-31页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108267221A (zh) | 2018-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108267221B (zh) | 一种消除串扰的振动传递路径测试方法 | |
CN107256204B (zh) | 基于传递函数的多点振动响应频域预测的实验装置及方法 | |
US7756649B2 (en) | Method for fault detection and diagnosis of a rotary machine | |
CN109902408B (zh) | 一种基于数值运算和改进的正则化算法的载荷识别方法 | |
JP5493373B2 (ja) | 伝達経路毎の成分を算定するための方法 | |
CN106528955A (zh) | 基于阻尼奇异值分解的运行工况传递路径分析方法 | |
Korkmaz et al. | Dynamic characterization of multi-axis dynamometers | |
Gajdatsy et al. | Critical assessment of Operational Path Analysis: mathematical problems of transmissibility estimation | |
Gajdatsy et al. | Critical assessment of Operational Path Analysis: effect of coupling between path inputs | |
Peukert et al. | Dynamic interaction between precision machine tools and their foundations | |
US9947361B2 (en) | Active vibration control device and design method therefor | |
CN107257914A (zh) | 用于检查惯性传感器的设备和方法 | |
CN110096779A (zh) | 一种伺服机构动特性分析方法 | |
Wang et al. | Fuzzy norm method for evaluating random vibration of airborne platform from limited PSD data | |
Cheng et al. | Transfer path analysis and contribution evaluation using SVD-and PCA-based operational transfer path analysis | |
Barton et al. | In-situ transfer path analysis of multiple vibration sources in a complex source-receiver assembly | |
Kim et al. | Indirect input identification by modal filter technique | |
Banwell et al. | Assessment of experimental techniques to characterise the vibration source strength of a motor radially mounted with resilient elements | |
Urgueira | On the rigid body properties estimation from modal testing | |
CN112883607A (zh) | 一种用于结构损伤测试的损伤模拟方法与系统 | |
Bianciardi | Critical assessment of OPA: effect of coherent path inputs and SVD truncation | |
CN111353329B (zh) | 一种传感器组合源表征的振源串扰消除方法 | |
Cheng et al. | AR model-based crosstalk cancellation method for operational transfer path analysis | |
Chen et al. | Application of Blocked Force Methodology in NVH development of Electrical Machines | |
Chakraborty et al. | Estimation of in-plane elastic parameters and stiffener geometry of stiffened plates |
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 |