机械系统的粘性或结构阻尼系数可以从空间上描述系统的阻尼分布，这对于阻尼结构减振设计具有重要的意义。本文研究基于频响函数矩阵法来辨识悬臂薄板结构系统阻尼系数。详细叙述了频响函数矩阵法的原理和测试流程。在此基础上组建了辨识结构系统阻尼系数的测试系统，对薄板结构进行了测试分析。从识别的阻尼矩阵可以看出，无论是粘性阻尼还是结构阻尼，薄板根部的数值均大于自由端的数值。比较不同频带的阻尼系数，可以看出在低频段粘性阻尼系数较大，而在高频段结构阻尼系数较大。 The distribution of system damp could be described by the internal friction and structure damped coefficient. It was important for design of reducing vibration of the damping structure. The damping coefficient of cantilever thin plate was identified by the frequency response function matrix method in this paper. The principle and test process were introduced. The test system was built to identify the damping coefficient. The thin plate was analyzed by the system. The data of the plate root were all bigger than that of freedom end no matter the damping matrix was inherent friction or the structural damping coefficient. The inner friction coefficient of low frequency band was bigger, while the structural damping coefficient of high frequency band was bigger.
朱全军1，樊习英2，曹枚根3，韩清鹏4，任建兴4，李天成4
1全球能源互联网研究院，北京
2国网山西省电力公司经济技术研究院，山西 太原
3中国电力科学研究院，北京
4上海电力学院能源与机械工程学院，上海
收稿日期：2017年5月25日；录用日期：2017年6月7日；发布日期：2017年6月14日
机械系统的粘性或结构阻尼系数可以从空间上描述系统的阻尼分布，这对于阻尼结构减振设计具有重要的意义。本文研究基于频响函数矩阵法来辨识悬臂薄板结构系统阻尼系数。详细叙述了频响函数矩阵法的原理和测试流程。在此基础上组建了辨识结构系统阻尼系数的测试系统，对薄板结构进行了测试分析。从识别的阻尼矩阵可以看出，无论是粘性阻尼还是结构阻尼，薄板根部的数值均大于自由端的数值。比较不同频带的阻尼系数，可以看出在低频段粘性阻尼系数较大，而在高频段结构阻尼系数较大。
关键词 :频响函数，矩阵法，薄板结构，阻尼系数
Copyright © 2017 by authors and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
机械系统的阻尼减振降噪技术离不开对结构的阻尼测试。阻尼测试是验证阻尼模型正确性、检验阻尼减振效果以及实施阻尼结构优化设计的重要技术基础。识别阻尼的方法包括频域法 [
传统上，通常用阻尼比、模态损耗因子、品质因子等来描述系统的阻尼特性，这些参数一般对应着某阶模态，用半功率带宽法可以求出具体的阻尼参数。但是，对于阻尼减振设计，单用阻尼比还不能满足研究需要。主要因为阻尼比不能从空间上描述系统的阻尼分布，也不能明确反映作用在系统上的阻尼机理。
要从空间上描述系统的阻尼分布，可以考虑用阻尼系数来表征。本文研究基于频响函数矩阵法来辨识悬臂薄板结构系统阻尼系数的方法，详细叙述了频响函数矩阵法的原理和测试流程，并用数值模型进行了验证。在此基础上组建了阻尼系数辨识测试系统，对实际结构进行测试，并对测试结果进行了分析，最后，对该识别算法的有效性进行了评价。
所谓频响函数矩阵法是通过测量获得结构系统的频响函数，再对频响函数应用相应的辨识算法来识别出系统的阻尼矩阵。阻尼矩阵可以描述系统的阻尼类型、阻尼大小、阻尼分布等特征，因而用阻尼矩阵来表征系统的阻尼特性，可以用于指导阻尼减振设计。另外，测试的阻尼矩阵还可以用于阻尼结构系统的动力学建模与分析。相对于质量及刚度，很难用一个合适的阻尼模型来预估出系统的阻尼特性，因而将测试获得的阻尼矩阵引入到动力学有限元分析中，可以提高结构系统动力学分析的精度。频响函数矩阵法又可分为基于动柔度矩阵 [
该识别算法是针对既具有粘性阻尼又具有结构阻尼的动力学系统进行阻尼系数辨识，其动力学方程可描述为
这里
假设
动刚度矩阵定义为
这里
实际应用中，总是先测试出频响函数，进一步通过求逆获得动刚度。
方程(3)的动刚度又可以分为实部及虚部
方程(5)和(6)又可以进一步写成
标记测量点，n个测量点，则需要
假如有必要，刚度和质量矩阵也可以被求出
假如仅用粘性阻尼来表示系统的阻尼特性，一个等价粘性阻尼矩阵Ceq求解式为
假如仅用结构阻尼来表示系统的阻尼特性，一个等价粘性阻尼矩阵
(1) 理论模型
图1为一个三自由度集中质量系统，相应的参数见表1。通过计算得出系统的频响函数，进一步利用频响函数反推出系统中的阻尼矩阵，从而检验算法的有效性。
(2) 识别算法校验
利用基于动刚度矩阵的识别算法进行校验。识别出的阻尼矩阵如表2所示。可见用频响函数识别出的阻尼矩阵与预设的相同，证明识别算法是有效的。
将TC4钛基薄板作为金属基层，其几何参数为长100 mm，宽80 mm，厚1.5 mm。材料参数为杨氏模量
搭建实验系统对涂层前后钛基悬臂薄板模态参数进行测试。分别对涂层前后钛基薄板梁进行测试，移动力锤敲击钛合金薄板，敲击顺序与电涡流传感器布置顺序一致，最后获得4 × 4的频响函数矩阵。针对测量点和敲击点进行测试，每点测试3次，取平均值作为测试结果，完成测试，如图2所示。
由于能量的衰减，力锤激励对于所研究的实验结构仅能达到3000 Hz，且仅在大约1000 Hz内能实现较高信噪比的测试。图1为试验敲击获得的频响，从中可以看出，大约1000 Hz以内存在三阶固有频率。
选择0~1000 Hz为测试范围，根据频响函数测试结果可知，该频段范围共包括了3阶固有频率。根据每阶固有频率的位置，将0~1000 Hz分成0~200 Hz、201~600 Hz、601~1000 Hz三个部分。受试验条件和试验设备的制约，试验测试所得钛合金薄板的频响函数是非对称矩阵，而理想状态下，在对应的激励点和响应点处获得的频响函数是相同的。例如：
图1. 三自由度集中质量系统
图2. 钛基薄板频响函数图
| 质量矩阵M/kg | 粘性阻尼矩阵C/Ns/m | 刚度矩阵K/N/m | 结构阻尼矩阵D/N/m |
|---|---|---|---|
|
|
|
|
|
表1. 三自由度系统参数表
频段下的阻尼系数，结果见表3、表4。
从识别结果可以看出：以悬臂薄板结构系统中，阻尼主要来自于结构阻尼，结构阻尼系数大约比粘性阻尼系数大两个数量级。结构阻尼的产生主要源于薄板装卡产生的结合部界面摩擦阻尼和材料内耗。而从实验中可以看出材料内耗阻尼份额较小，因而无法用上述识别方法比较涂层前后阻尼系数的大小。
| 识别方法 | 粘性阻尼矩阵/Ns/m | 结构阻尼矩阵/N/m |
|---|---|---|
| 动刚度矩阵法 |
|
|
表2. 两种算法求解三自由度系统的阻尼矩阵
| 频率范围/Hz | 粘性阻尼系数/Ns/m | 结构阻尼系数/N/m |
|---|---|---|
| 0~200 |
|
|
| 201~600 |
|
|
| 601~1000 |
|
|
表3. 利用对称动刚度矩阵识别未涂层钛基薄板的阻尼系数
| 频率范围/Hz | 粘性阻尼系数/Ns/m | 结构阻尼系数/N/m |
|---|---|---|
| 0~200 |
|
|
| 201~600 |
|
|
| 601~1000 |
|
|
表4. 利用对称动刚度矩阵识别涂层后钛合金薄板的阻尼系数
从识别的阻尼矩阵可以看出，无论是粘性阻尼还是结构阻尼，薄板根部的数值均大于自由端的数值。比较不同频带的阻尼系数，可以看出在低频段粘性阻尼系数较大，而在高频段结构阻尼系数较大。
机械系统的粘性或结构阻尼系数可以从空间上描述系统的阻尼分布，这对于阻尼结构减振设计具有重要的意义。本文研究基于频响函数矩阵法来辨识悬臂薄板结构系统阻尼系数。
详细叙述了频响函数矩阵法的原理和测试流程。用一个三自由度的集中质量模型，做了原理性验证。结果表明：用频响函数识别出的阻尼矩阵与预设的相同，证明识别算法是有效的。对涂层前后钛基悬臂薄板模态参数进行测试。力锤激励对于所研究的实验结构仅能达到3000 Hz，且仅在大约1000 Hz内能实现较高信噪比的测试。通过对称动刚度矩阵可以识别涂层后钛合金薄板的阻尼系数。
由国家电网科技项目(SGRI-WD-71-15-010)支持。
朱全军,樊习英,曹枚根,韩清鹏,任建兴,李天成. 基于频响函数矩阵法的薄板结构阻尼系数识别Recognition of Structural Damping Coefficient of Thin Plate Based on the Frequency Response Matrix Method[J]. 机械工程与技术, 2017, 06(02): 133-139. http://dx.doi.org/10.12677/MET.2017.62019