作者简介:孙伟(1975-), 男, 副教授,博士生导师.研究方向:机械系统动力学及振动控制.E-mail:weisun@mail.neu.edu.cn
针对硬涂层复合结构共振衰减的时域响应,研究了基于Hilbert变换辨识硬涂层复合结构的非线性刚度及阻尼的方法。首先,确定了硬涂层复合结构共振状态下自由衰减的运动方程,进而推导了辨识该复合结构非线性刚度及阻尼参数的原理性公式;其次,用数值算例证明了辨识算法的正确性;然后,提出了将辨识算法应用于硬涂层复合结构非线性刚度及阻尼参数辨识的流程;最后,进行了实验验证,辨识出单面涂敷NiCrAlY+YSZ硬涂层的悬臂钛板非线性刚度及阻尼参数,并与未涂层试件进行了比对。结果表明,涂覆硬涂层后钛板的刚度非线性增强,阻尼增大,硬涂层具有减振效果。
Based on the resonance decay response of time domain, a method to identify the nonlinear stiffness and damping was proposed using Hilbert transform. First, the free decay motion equation of the hard-coating composite structure under resonance status was determined and the formulas to identify nonlinear stiffness and damping parameters were derived. A numerical example was used to verify the correctness of the proposed identification algorithm. Then, the procedure of identifying the above parameters was presented. Finally, a titanium plate coated with NiCrAIY + YSZ hard coating was chosen as a study case, the nonlinear stiffness and damping parameters of the composite plate were identified and compared with the results of uncoated plate. It is shown that the hard coating can enhance the nonlinearity and increasing the damping of the structure, thus, improving the effect of vibration reduction.
硬涂层是指用金属或陶瓷制成的涂层材料, 由于硬涂层可以在高温、高腐蚀环境对薄壳结构进行减振, 因而近年来受到了越来越广泛的关注[1, 2, 3]。硬涂层材料具有非线性特点, 其杨氏模量及耗因子都会随着系统激励幅度的不同而发生改变[4, 5]。正因为硬涂层材料具有上述特点, 致使涂敷硬涂层的复合结构表现出变刚度及变阻尼的非线性振动特性。
有效辨识出这些非线性刚度及阻尼参数是硬涂层复合结构动力学分析以及硬涂层减振有效性评价的基础, 但是目前还没有提出有效的辨识算法来实现这一目标。Hilbert变换广泛应用于结构动力学特性参数的识别, 例如, 赵岩等[6]建立了基于扭摆法的时变转动惯量测量模型, 提出了基于Hilbert变换辨识时变转动惯量计算方法。Feldman[7, 8]推导了基于Hilbert变换和经验模态分解(Empirical mode decomposition, EMD)辨识动力学系统振动特性参数的方法。由此可以推断, Hilbert变换技术也可应用于硬涂层复合结构非线性刚度及阻尼参数的辨识。
本文针对硬涂层复合结构共振衰减的时域响应, 研究基于Hilbert变换辨识涂层复合结构非线性刚度和阻尼参数的方法。首先, 确定了硬涂层复合结构共振状态下自由衰减的运动方程, 进而推导了辨识该复合结构非线性刚度及阻尼参数的原理性公式; 其次, 用数值算例证明了辨识算法的正确性; 然后, 提出了将辨识算法应用于硬涂层复合结构非线性刚度及阻尼参数辨识的流程; 最后, 进行了实验验证, 辨识出单面涂敷NiCrAlY+YSZ硬涂层悬臂钛板的非线性刚度及阻尼参数, 并与未涂层薄板进行了比对。
结构处于共振状态时, 当各阶模态耦合不严重时, 其动力学特性可由单一模态来确定[9]。当切断激振力时, 由于存在阻尼, 结构上某点的响应信号会逐渐衰减。对于硬涂层复合结构, 由于其刚度及阻尼参数均随着系统的振幅而发生改变, 因而对于这一衰减过程, 可以认为复合结构的刚度及阻尼随着时间发生改变。这样, 可将上述过程用如下单自由度方程来描述:
式中: 为复合结构的模态质量c(t)和k(t)分别表示硬涂层复合结构的时变模态阻尼及刚度q(t); 为对应模态坐标的振动衰减响应。以下基于Hilbert变换技术推导辨识非线性刚度及阻尼的辨识公式。
对衰减响应信号
式中:H(· )表示Hilbert变换;
由变换后的信号和原始衰减响应信号组成分析信号, 表达为:
则分析信号
对瞬时相位
对式(1)进行Hilbert变换, 可以得到:
需要说明的是在进行Hilbert变换时, 这里假设时变阻尼
将式(7)两端乘以j加到式(1), 可得到用分析信号表达的运动方程, 即:
利用式(3)(4)(5), 可将式(10)中分析信号一阶及二阶导数
将式(11)(12)代入到式(10), 并整理得:
令方程(13)的实部和虚部都等于零, 可获得非线性阻尼及非线性刚度的表达式分别如下:
需要说明的是上述刚度及阻尼辨识原理主要适用于具有低模态密度的结构系统, 即各模态之间耦合效应小, 因而可用单自由度方程来描述某一阶次的运动状态。如果结构具有密集频率特性, 模态耦合严重, 则所提出的辨识方法失效。
以下用数值算例校验辨识算法的正确性, 设时变刚度(单位:N/m)及阻尼(单位:(N· s)/m)的表达式分别为
首先, 利用Newton Raphson法, 并设初始条件为
进一步对q(t)进行Hilbert变换, 并获得包络线A(t)和即时频率ω(t),如图2所示。最后利用式(14)(15)辨识得出系统的时变刚度和时变阻尼, 辨识结果与已知值的对比见图3。
从图3可以看出, 最终辨识的时变刚度及阻尼与预设值差距不大, 从而证明了上述辨识算法的正确性。
按照辨识算法, 测试硬涂层复合结构非线性刚度及阻尼, 可参照如下流程进行。
(1)测试硬涂层复合结构件的共振衰减响应
在找到硬涂层复合结构各阶共振频率的前提下, 利用稳态激励装置(例如电磁振动台)以共振频率作为激振频率激励复合结构件, 待达到稳态振动后, 切断激励装置, 记录衰减位移响应。
(2)计算拾振点位置处的模态质量
不同的拾振点对应结构的模态质量不同, 需要利用结构件的模态振型计算出该拾振点处的模态质量。对于单面全部涂敷硬涂层的悬臂板(见图4), 可以认为涂层前后模态振型没有改变。按照双向梁函数法悬臂薄板的振型函数可表示为[10]:
式中:
式中: S为薄板的面积;ρd 、ρc分别为基体和硬涂层的密度;Hd 、Hc 分别为基体和硬涂层的厚度。
(3)对获取的振动响应进行EMD分解
Hilbert变换需要信号在任意时刻具有单频性, 而测试获得的信号在某一时刻不能完全保证处于单一频率, 因而需要对原始信号进行EMD分解。此外, 对原始测得的信号进行EMD分解也能用于判断测得的衰减响应是否为共振衰减响应, 因为实践发现结构件的共振与非共振衰减响应经EMD分解后存在明显的差异。例如, 图5为对硬涂层板共振及非共振响应进行EMD分解获得的响应信号, 可发现由共振响应EMD分解后得到的是单一频率的主体信号和一些振动很小的扰乱信号, 而由非共振响应EMD分解后得到的是多个频率的振动信号。
(4)获得模态坐标下的振动响应
可按照模态理论将真实坐标下的响应变换为模态坐标下的响应。同样以硬涂层板为例, 变换公式为:
式中:
(5)进行Hilbert变换, 获取
对模态坐标下的振动响应
(6)带入公式获得时变刚度及阻尼
将获得包络线A(t)和即时频率ω(t)带入到式(14)(15)则可获得被测硬涂层复合结构的非线性阻尼及刚度。
需要说明的是, 上述非线性刚度及阻尼辨识流程同样也适用于未涂层结构。受到夹持边界条件的影响, 未涂层结构也会产生上述变刚度及变阻尼的非线性现象, 只是非线性的程度会有所减弱。
以未涂层的和单面涂敷NiCrAlY+YSZ硬涂层的悬臂钛板为研究对象, 来检验本文方法的可行性, 同时对比涂层前后薄板动力学参数的变化。具体试验件见图6, 硬涂层及薄板的材料及相关尺寸见表1。需要说明的是, 硬涂层材料的杨氏模量不是一个固定值, 其大小随激励幅度的不同而改变, 表1中给定值是参照文献[11]给出的大致范围。
所组配的实验系统见图7, 板以悬臂的方式固定在夹具上(夹持区为20 mm), 测试装置包括:金盾EM-1000F振动台、LMS SCADAS采集前端、金盾EM-1000F电磁振动台和Polytec PLD-100激光测振仪等。
这里以在1.5 g激励下的钛板第6阶共振特性为例对本文算法进行验证。经扫频测试, 未涂层钛板第6阶的固有频率为1297.88 Hz, 涂层后第6阶的固有频率为1262.25 Hz。分别按以上固有频率定频激励涂层前后的钛板, 待达到稳定后切掉激励源, 获得的共振状态的自由振动衰减响应见图8。
按照第3节所述的流程辨识该硬涂层薄板的时变刚度及阻尼。其中对应于拾振点处(x=0.029 m, y=0.041 m), 由式(17)计算的涂层前第6阶的模态质量为0.1250 kg, 涂层后的模态质量为0.1617 kg。根据式(18), 指定点处模态坐标与真实坐标之间变换系数为-0.9266。
首先对测得的共振衰减响应信号进行EMD分解, 提取主体信号。进而用获得的变换系数将真实坐标下的响应变换为模态坐标下的响应。然后, 对模态坐标响应进行Hilbert变换, 利用式(4)和式(5), 识别出复合结构的包络线和即时频率, 分别见图9和图10。
最后由式(14)和式(15)辨识出硬涂层薄板的非线性阻尼和刚度, 由测得的阻尼可进一步得出模态阻尼比, 分别如图11和图12所示。
从测试结果可以看出, 从0.3 s以后, 辨识结果发散, 结果的可信度降低。发散的原因主要来自于测试获得的时域信号在高响应幅度时信噪比高, 而在低响应幅度时容易受到外界因素(例如电磁噪声)干扰而使信噪比降低。低信噪比信号将失去原始结构的力学信息, 导致辨识结果发散。因而后续需要进一步研究低响应幅度时, 针对时域衰减信号的降噪技术。
大约在0~0.3 s的时间范围内, 结果较为稳定, 可以反应出一定的规律, 因而本文的测试也具有一定的意义。对比涂层前后0~0.3 s时间范围内钛板的刚度及阻尼参数变化情况可以看出:①未涂层板在这一时间段内刚度和即时频率基本不变, 而涂层板的刚度和即时频率发生显著变化, 具体为随响应幅度的减小而快速增大; ②未涂层板和涂层板在这一时间段内的阻尼都发生变化, 表明阻尼与响应幅度相关, 此外涂层板的阻尼大于未涂层板的阻尼。
以上现象说明, 涂覆硬涂层后钛板的刚度非线性增强, 阻尼增大, 硬涂层具有减振效果。实际上, 从共振衰减响应(见图8)中也可以看出, 硬涂层具有减振效果, 涂层后共振响应的幅值显著降低。
(1)将硬涂层复合结构简化为单自由度分析模型, 并基于Hilbert变换可辨识出硬涂层复合结构的非线性刚度及阻尼。基于本文提出的辨识流程, 可有效测试出涂层前后钛板的时变刚度和时变阻尼。
(2)辨识算法对信号的信噪比要求较高, 对于本文所研究的实例, 在0.3 s以后的数据结果发散, 可信度降低。因而后续需要进一步研究低响应幅度时, 针对时域衰减信号的降噪技术。
(3)对比钛板涂覆NiCrAlY+YSZ硬涂层前后的时变刚度和时变阻尼的变化可以看出:硬涂层使钛板的刚度非线性增强, 阻尼增大, 并有效抑制了复合结构的共振响应。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|