基于总体最小二乘法的工况传递路径分析

桂俊涛; 樊晓波; 满石林; 吴凇

doi:10.11993/j.issn.2096-3920.2025-0024

基于总体最小二乘法的工况传递路径分析

doi: 10.11993/j.issn.2096-3920.2025-0024

中国船舶集团有限公司第七〇五研究所, 云南昆明, 650106

详细信息

作者简介:
桂俊涛, 男, 硕士, 工程师, 主要研究方向为振动噪声控制

中图分类号: TJ630; U661.44
计量
- 文章访问数: 3
- HTML全文浏览量: 2
- PDF下载量: 1
- 被引次数: 0
出版历程
- 收稿日期: 2025-02-13
- 修回日期: 2025-09-02
- 录用日期: 2025-09-08
- 网络出版日期: 2026-03-27

Operational transfer path analysis with total least-squares method

The 705 Research Institute, China Shipbuilding Industry Corporation Limited, Kunming 650106, China

摘要

摘要: 工况传递路径分析(OTPA)通过不同工况下的响应数据, 实现对振动噪声的分解和预测, 因此广泛应用于各类工程领域中。但振动噪声的响应数据不可避免包含误差, 严重影响OTPA的准确性。为减小误差的影响及提高传递率函数矩阵的准确性, 采用总体最小二乘法来估计传递率函数矩阵, 再进行OTPA。相较传统方法, 总体最小二乘法同时考虑了目标点与指示点数据的误差。在数值模型与测试模型中开OTPA, 分别应用最小二乘模型和总体最小二乘模型获取各路经贡献量。仿真结果显示总体最小二乘法识别的贡献量与经典传递路径分析贡献量吻合度更高, 表明相比正则化最小二乘法, 总体最小二乘法在OTPA中适用性更好, 有效提高了OTPA的准确度。
- 工况传递路径分析 /
- 观测误差 /
- 总体最小二乘法 /
- 振动噪声
Abstract: The operational transfer path analysis(OTPA) utilizes response data under different working conditions to decompose and predict vibration and noise, and is therefore widely used in various engineering fields. However, the response data of vibration noise inevitably contains errors, which seriously affect the accuracy of the transfer path analysis. To reduce the impact of errors and improve the accuracy of the transfer rate function matrix, the total least squares method is used to estimate the transfer rate function matrix, and then the analysis of the working condition transfer path is carried out. Compared to traditional methods, the overall least squares method takes into account the errors in both the target point and indicator point data. In numerical model simulations and experimental model of open condition transmission path analysis, least squares method and overall least squares method were applied to obtain the contribution of each path. The simulation results show that the contribution identified by the overall least squares method is more consistent with the contribution of classical transmission path analysis, indicating that compared with the regularized least squares method, the total least-squares method has better applicability in the analysis of transmission paths under working conditions, effectively improving the accuracy of transmission path analysis under working conditions.
- operational transfer path analysis /
- observation error /
- total least-squares method /
- vibration and noise

HTML全文

图 1 主被动线性系统

Figure 1. Active and passive linear systems

下载: 全尺寸图片幻灯片

图 2 九自由度模型

Figure 2. The nine-degree-of-freedom model system

下载: 全尺寸图片幻灯片

图 3 拆解九自由度模型

Figure 3. The split nine-degree-of-freedom model system

下载: 全尺寸图片幻灯片

图 4 合成响应对比

Figure 4. Comparison of synthetic responses

下载: 全尺寸图片幻灯片

图 5 路径35贡献量对比

Figure 5. Comparison of Contribution of Path 35

下载: 全尺寸图片幻灯片

图 6 路径46贡献量对比

Figure 6. Comparison of contribution of path 46

下载: 全尺寸图片幻灯片

图 7 路径78贡献量对比

Figure 7. Comparison of contribution of path 78

下载: 全尺寸图片幻灯片

图 8 实验装置

Figure 8. The experimental setup

下载: 全尺寸图片幻灯片

图 9 频响函数H_ya和H_yb

Figure 9. The FRF H_yaand H_yb

下载: 全尺寸图片幻灯片

图 10 TPA合成响应与实测响应对比

Figure 10. The Comparison of TPA synthetic response and test response

下载: 全尺寸图片幻灯片

图 11 O-RLS合成响应与实测响应对比

Figure 11. Comparison of o-rls synthetic response and test response

下载: 全尺寸图片幻灯片

图 12 O-TLS合成响应与实测响应对比

Figure 12. Comparison of o-tls synthetic response and measured response

下载: 全尺寸图片幻灯片

图 13 路径1贡献量

Figure 13. The individual contribution of path 1

下载: 全尺寸图片幻灯片

图 14 路径2贡献量

Figure 14. The individual contribution of path 2

下载: 全尺寸图片幻灯片

图 15 路径贡献量云图

Figure 15. Path contribution cloud chart

下载: 全尺寸图片幻灯片

图 16 FRAC值对比

Figure 16. The comparison of FRAC value

下载: 全尺寸图片幻灯片

表 1 九自由度模型参数列表

Table 1. Parameter list of the 9-DOF model

参数类型	参数标识及数值	单位
质量	M₁=2; M₂=1; M₃=5.5; M₄=3; M₅=6; M₆=8; M₇=2.5; M₈=3.5; M₉=4;	kg
刚度	K₀₁=0.1; K₁₂=0.4; K₁₇=0.5; K₂₃=0.35; K₂₄=0.15; K₃₅=0.13; K₄₆=0.6; K₅₉=0.27; K₆₉=0.39; K₇₈=0.48; K₈₉=0.29;	kN/mm
阻尼	C₀₁=1.1; C₁₂=3.6; C₁₇=6.2; C₂₃=4.5; C₂₄=9.0; C₃₅=8.3; C₄₆=11.3; C₅₉=14; C₆₉=16.3; C₇₈=5.2; C₈₉=7.0;	kg/s

下载: 导出CSV

表 2 实验仪器参数表

Table 2. Experimental Instrument parameters

仪器	数量	标称灵敏度	采样频率/Hz
B&K加速度计(4524 B)	3	100 mv/g	0.25 ~3 000
CL-YD力传感器	2	3.25 N/Pc	1 000
JZK-2 激振器	1	-	15 000
JZK-20 激振器	1	-	2 000
B&K数据采集处理系统	1	-	25 600

下载: 导出CSV

表 3 实验工况

Table 3. Experimental conditions

工况	点A激励	点B激励
1	随机	随机
2	随机	正弦扫频
3	正弦扫频	随机
4	正弦扫频	正弦扫频

下载: 导出CSV

参考文献(25)

[1]	王路, 尹韶平, 曹小娟, 等. 基于Workbench响应谱分析的鱼雷振动传递优化[J]. 鱼雷技术, 2016, 24(1): 13-17. Wang L, Yin S P, Cao X J, et al. Vibration transfer optimization of torpedo based on workbench response spectrum[J]. Journal of Unmanned Undersea Systems, 2016, 24(1): 13-17.
[2]	王红瑞, 曹小娟, 尹韶平, 等. 鱼雷振动试验系统结构传递特性影响因素分析[J]. 水下无人系统学报, 2019, 27(5): 574-579. doi: 10.11993/j.issn.2096-3920.2019.05.014 Wang H R, Cao X J, Yin S P, et al. Analysis on the factors influencing structural transfer characteristic of torpedo vibration test system[J]. Journal of Unmanned Undersea Systems, 2019, 27(5): 574-579. doi: 10.11993/j.issn.2096-3920.2019.05.014
[3]	Yang S, Chang H, Wang Y, et al. A phononic crystal suspension for vibration isolation of acoustic loads in underwater gliders[J]. Applied Acoustics, 2024, 216: 109731. doi: 10.1016/j.apacoust.2023.109731
[4]	Sun Q, Yang Y, Wu P, et al. Study on the vibration-damping mechanism of a new phononic crystal suspension equipped on underwater gliders[J]. Journal of Marine Science and Engineering, 2024, 12(11): 2088. doi: 10.3390/jmse12112088
[5]	孙旭阳, 周景军, 王谦, 等. 基于声学超材料的鱼雷动力舱段减振方法[J]. 水下无人系统学报, 2024, 32(6): 1072-1081. Sun X Y, Zhou J J, Wang Q, et al. Vibration reduction method for power cabin of torpdoes based on acoustic metamaterials[J]. Journal of Unmanned Undersea Systems, 2024, 32(6): 1072-1081.
[6]	马皋, 陈亚明, 沈德明, 等. 基于工况传递路径的旋转机械振动源分析[J]. 振动与冲击, 2022, 41(17): 276-281. doi: 10.13465/j.cnki.jvs.2022.17.034 Ma G, Chen Y M, Shen D M, et al. Vibrtion source analysis of rotating machinery based on working condition transmission path[J]. Journal of Vibrtion and Shock, 2022, 41(17): 276-281. doi: 10.13465/j.cnki.jvs.2022.17.034
[7]	Acri A, Offner G, Nijman E, et al. Substructuring of multibody systems for numerical transfer path analysis in internal combustion engines[J]. Mechanical Systems and Signal Processing, 2016, 79: 254-270. doi: 10.1016/j.ymssp.2016.02.034
[8]	Oktav A, Yılmaz Ç, Anlaş G. Transfer path analysis: Current practice, trade-offs and consideration of damping[J]. Mechanical Systems and Signal Processing, 2017, 85: 760-772. doi: 10.1016/j.ymssp.2016.09.013
[9]	Thite A N, Thompson D J. The quantification of structure-borne transmission paths by inverse methods. part 2: Use of regularization techniques[J]. Journal of Sound and Vibration, 2003, 264(2): 433-451. doi: 10.1016/S0022-460X(02)01203-8
[10]	Ye S, Hou L, Zhang P, et al. Transfer path analysis and its application in low-frequency vibration reduction of steering wheel of a passenger vehicle[J]. Applied Acoustics, 2020, 157: 107021. doi: 10.1016/j.apacoust.2019.107021
[11]	Maarten V, Dennis D K, Daniel J. R. General framework for transfer path analysis: History, theory and classification of techniques[J]. Mechanical Systems and Signal Processing, 2016, 68-69: 217-244.
[12]	Lage Y E, Neves M M, Maia N M M, et al. Force transmissibility versus displacement transmissibility[J]. Journal of Sound and Vibration, 2014, 333(22): 5708-5722. doi: 10.1016/j.jsv.2014.05.038
[13]	Lage Y E, Maia N M M, Neves M M, et al. Force identification using the concept of displacement transmissibility[J]. Journal of Sound and Vibration, 2013, 332(7): 1674-1686. doi: 10.1016/j.jsv.2012.10.034
[14]	De Sitter G, Devriendt C, Guillaume P, et al. Operational transfer path analysis[J]. Mechanical Systems and Signal Processing, 2010, 24(2): 416-431. doi: 10.1016/j.ymssp.2009.07.011
[15]	Gajdatsy P, Janssens K, Desmet W, et al. Application of the transmissibility concept in transfer path analysis[J]. Mechanical Systems and Signal Processing, 2010, 24(7): 1963-1976. doi: 10.1016/j.ymssp.2010.05.008
[16]	Klerk D, Ossipov A. Operational transfer path analysis: Theory, guidelines and tire noise application[J]. Mechanical Systems and Signal Processing, 2010, 24(7): 1950-1962. doi: 10.1016/j.ymssp.2010.05.009
[17]	Choi H G, Thite A N, Thompson D J. Comparison of methods for parameter selection in tikhonov regularization with application to inverse force determination[J]. Journal of Sound and Vibration, 2007, 304(3-5): 894-917. doi: 10.1016/j.jsv.2007.03.040
[18]	Cheng W, Lu Y, Zhang Z. Tikhonov regularization-based operational transfer path analysis[J]. Mechanical Systems and Signal Processing, 2016, 75: 494-514. doi: 10.1016/j.ymssp.2015.12.025
[19]	Liu Y, Shepard W S. Dynamic force identification based on enhanced least squares and total least-squares schemes in the frequency domain[J]. Journal of Sound and Vibration, 2005, 282(1-2): 37-60. doi: 10.1016/j.jsv.2004.02.041
[20]	Jia Y, Yang Z, Guo N, et al. Random dynamic load identification based on error analysis and weighted total least squares method[J]. Journal of Sound and Vibration, 2015, 358: 111-123. doi: 10.1016/j.jsv.2015.07.035
[21]	鲁铁定, 周世健. 总体最小二乘的迭代解法[J]. 武汉大学学报(信息科学版), 2010, 35(11): 1351-1354. Lu T D, Zhou S J. An iterative algorithm for total least squares estimation[J]. Geomatics and Information Science of Wuhan University, 2010, 35(11): 1351-1354.
[22]	Wang Z, Zhu P. A system response prediction approach based on global transmissibilities and its relation with transfer path analysis methods[J]. Applied Acoustics, 2017, 123: 29-46. doi: 10.1016/j.apacoust.2017.02.016
[23]	Vaitkus D, Tcherniak D, Brunskog J. Application of vibro-acoustic operational transfer path analysis[J]. Applied Acoustics, 2019, 154: 201-212. doi: 10.1016/j.apacoust.2019.04.033
[24]	唐中华, 昝鸣, 张志飞, 等. 广义Tikhonov正则化工况传递路径分析[J]. 振动与冲击, 2022, 41(24): 270-277. Tang Z H, Zan M, Zhang Z F, et al. Operational transfer path analysis with generalized Tikhonov regularization method[J]. Journal of Vibrtion and shock, 2022, 41(24): 270-277.
[25]	Shin K. An alternative approach to measure similarity between two deterministic transient signals[J]. Journal of Sound and Vibration, 2016, 371: 434-445. doi: 10.1016/j.jsv.2016.02.037