基于时频分析的高温振动环境螺栓连接件松动判别_曹芝腑.pdf

 振动与冲击 第 38 卷第 17 期JOURNAL OF VIBRATION AND SHOCKVol．38 No．17 2019 基金项目 国家自然科学基金 11572086; 11602112 收稿日期 2018 －01 －23修改稿收到日期 2018 －05 －15 第一作者 曹芝腑 男， 博士生， 1990 年生 通信作者 费庆国 男， 博士， 教授， 1977 年生 基于时频分析的高温振动环境螺栓连接件松动判别 曹芝腑1， 2，谭志勇3，姜东1， 4，何顶顶1， 2，费庆国1 1． 东南大学 空天机械动力学研究所， 南京211189; 2． 东南大学 工程力学系， 南京210096; 3． 北京临近空间飞行器系统工程研究所， 北京100076; 4． 南京林业大学 机械电子工程学院， 南京 210037 摘 要 针对某型复合材料螺栓连接件， 研究其在高温环境下的振动特性， 为表征其松动特性提供理论和试验依 据。利用振动控制系统， 温度控制系统和非接触式激光测量系统实现了 1 000 ℃高温环境下的复合材料螺栓连接件振动 试验。并利用短时傅里叶方法提取高温环境下的振动响应特征量， 结合复合材料连接件的振动基频， 二倍频幅值和二倍 频与基频幅值的比值对连接结构的松动状态进行综合评判。试验结果表明， 在连接件宏观松动现象非明显情况下， 所提 评判方法能够综合评估复合材料螺栓连接件在高温振动环境下的松动特性。 关键词 复合材料; 高温振动; 连接松动; 时频分析 中图分类号 TH131． 9; O327文献标志码 ADOI10． 13465/j． cnki． jvs． 2019． 17． 028 Loosening discrimination of a bolted connector under high- temperature vibration environment based on time- frequency analysis CAO Zhifu1， 2，TAN Zhiyong3，JIANG Dong1， 4，HE Dingding1， 2，FEI Qingguo1 1． Institute of Aerospace Machinery and Dynamics，Southeast University，Nanjing 211189，China; 2． Department of Engineering Mechanics，Southeast University，Nanjing 210096，China; 3． Beijing Institute of Nearspace Vehicle’ s Systems Engineering，Beijing 100076，China; 4． College of Mechanical and Electronic Engineering，Nanjing Forestry University，Nanjing 210037，China AbstractHere，vibration features of a certain type composite bolted connector were studied under high- temperature vibration environment to provide theoretical and test basis for judging its loosening features． Using a vibration control system，a temperature control system and a non- contact laser measurement system，vibration tests of the connector under 1 000 ℃ high- temperature environmentwere conducted． The short- time Fourier transation was applied to extract its vibration response characteristic quantities under high- temperature environment．The fundamental natural frequency，the amplitude of the second natural frequency component and the amplitude ratio between those of the second one and the fundamental one of the bolted connecter were used to judge its loosening state． The test results showed that under the condition of its macroscopic loosening phenomenon being not obvious，the proposed can be used to comprehensively uate loosening features of the composite bolted connector under high- temperature vibration environment． Key words composites; high temperature vibration; connection loosening; time- frequency analysis 航空航天、 能源等诸多领域由于其高温高压及强 振动的严苛工作环境 ［1 ］， 需要连接结构在此环境下仍 能保持足够的刚度、 强度、 稳定性及耐久性， 传统材料 很难满足这些要求。因此复合材料由于其比重小、 高 温环境下比刚度高、 比强度大及服役寿命久等诸多优 点， 得到越来越多国内外工程研究人员的广泛重视和 研究， 并在相关领域得到广泛使用和发展 ［2- 5 ］。随着复 合材料在工程中的广泛应用， 复合材料结构设计越来 越复杂， 加之由于成型工艺水平、 设计和使用维护等方 面的限制， 工程中往往需要将复合材料结构件用连接 技术连为整体， 复合材料结构的连接问题不可避免。 对于复合材料连接结构设计， 连接部位是一个相对薄 弱的环节， 结构的连接破坏是导致结构失效的主要原 因之一， 尤其是某些重要连接部位的失效破坏， 有可能 会危及整体结构的安全性和可靠性。 复合材料结构的连接形式主要有 3 种 ［6 ］ 胶接连 ChaoXing 接、 机械连接和混合连接。机械连接有螺栓连接、 销钉 连接和铆钉连接等方式， 其中螺栓连接具有可靠性高、 承载能力强、 便于重复拆装及使用维护简便等优点 ［7 ］， 是航空航天领域中复合材料结构的主要机械连接方 式。基于复合材料螺栓连接的优点和存在的缺点， 有 必要对其连接特性进行深入分析， 改善不利因素的影 响， 提高连接效率 ［8 ］。 松动失效是最常见的连接结构失效模式， 近年来 有不少国内外学者致力于连接结构松动诱致失效、 故 障的研究。研究主要分为两方面， 一方面建立连接结 构松动动力学模型， 分析连接结构松动对系统动力学 特性的影响 ［9- 10 ］; 另一方面是基于结构动响应分析的方 法 ［11- 12 ］， 利用动力学实验中结构非线性动响应， 结合非 线性动力学理论提取描述结构松动状态的特征参数， 建立结构松动状态与非线性特征参量之间的关系， 进 行连接结构松动状态监测和辨识。 L 等 ［13 ］通过试验建立螺栓残余应力与温度和初 始预载荷之间的关系， 对复合材料- 复合材料， 复合材 料- 金属和金属- 金属等不同形式连接结构的松动行为 进行分析。该方法能够利用短期试验数据构造半解析 的时间- 温度相关蠕变率方程， 并利用该方程对螺栓连 接结构的预紧力松弛现象进行长期预测。赵登峰等 ［14 ］ 建立了螺纹连接松动高度简化的动力学模型， 通过数 值仿真和实验相结合的方法， 研究了振动环境中螺栓 连接松动过程的动响应特性， 提出将谐波失真度作为 螺栓连接松动失效的预警指标。Li 等 ［15 ］针对具有螺 栓松动特征的卫星结构， 在有限传感器的条件下， 利用 二阶频谱对其松动故障进行识别。Nguyen 等 ［16 ］利用 结构健康监测的方法， 利用不同螺栓松动工况下的振 动响应数据， 对风力发电机塔结构进行螺栓松动损伤 识别。Amerini 等 ［17 ］使用健康监测系统研究了螺栓结 构在松动和收紧状态时连接结构接触面的线性及非线 性声学参数， 使用一阶声学矩作为螺栓松紧状态的指 数， 这一方法能够有效评估螺栓连接的健康状态。另 外， 研究了单频和多频激励下， 高次谐波产生及边带调 制的非线性方法， 其中在单频激励下， 结构信号的二次 谐波幅值与载体信号的基频振幅的比值可以作为一个 可靠的指标， 有效评估螺栓的装配状态。 本文进行了高温振动环境下复合材料连接件螺栓 松动试验研究， 并使用短时傅里叶变换的时频分析技 术对高温振动响应进行分析， 得到了复合材料连接结 构在各试验条件下的时频结果。采用二倍频与基频振 动对应的幅值比作为松动特征值， 对复合材料螺栓连 接结构在高温环境下的振动松动状态进行了综合 评估。 1松动特征量 图 1 所示为螺栓连接件结构， 在高温环境下， 结构 除产生 x 方向伸缩变形外， 还会产生 z 向变形。同时由 于螺栓的约束和接触面的粗糙 ［18 ］， 使得结构连接区域 的相互接触刚度具有非线性特性。一种可行的处理方 案是将连接位置的非线性接触关系简化为非线性集中 力 ［19 ］， 利用非线性弹簧进行简化建模， 如图 2 所示。 考虑如图 2 所示的均质等截面简化梁模型， 其在 位移为 z t 的基础振动下的动力学微分方程为 2w x， t t 2 EI ρS 4w x， t x 4 1 ρS fnlδ x， l 2z t t 2 1 式中 w 为悬臂梁挠度， 由线性项 wlin和非线性项 wnlin 组合而成， 即 w x， t wlin x， t wnlin x， t ，E 为弹 性模量， I 为梁横截面惯性矩， ρ 为密度， S 为横截面积， δ x， l 为狄拉克函数， fnl为非线性集中恢复力， z t 为 基础振动位移。针对均质等截面梁， ρS 为常数， 则式 1 可简化为 2w x， t t 2 EI ρS 4w x， t x 4 FNL 2 式中 FNL表示 t 时刻的等效非线性载荷， FNL 2z t t 2 x∈ 0， l ∪ l， 2l 2z t t 2 － k1w l， t k2w l， t 2 x {l 3 其中 k1和 k2为由于温度变化引起的螺栓连接界 面处的一次和二次刚度系数， 是温度 T 的函数 k1 f T ， k2 g T 。 图 1螺栓连接结构 Fig． 1Bolted connection structure 图 2等效非线性弹簧- 梁模型 Fig． 2Equivalent nonlinear spring- beam model 文献［ 20］ 指出在该等效建模前提下， 结构的非线 性动响应包含二次谐波项。在频域下， 定义无量纲量 r Z 为二次谐波幅值与基频幅值的比， 作为螺栓连接 结构的松动特征指标， 对松动特性进行评判 r Z Z2 2ωn Z1 ω n 4 602振 动 与 冲 击2019 年第 38 卷 ChaoXing 式中， Z 为振动幅值 加速度， 速度， 位移 。 2时频域分析松动判别 时频域分析方法是非平稳信号处理的一个重要分 支， 它是利用时间和频率的联合函数来表示非平稳信 号， 并对其进行分析和处理。时频分析的主要任务是 描述信号的频谱含量随时间的变化规律。时频分析的 最终目的是建立一种分布， 以便在时间和频率上同时 表示信号的能量或强度， 得到这种规律后， 就可以对连 接结构振动信号进行分析和处理， 从而提取信号中所 包含的特征信息。 傅里叶变换难以反映信号频率成分的时间特性。 但是， 如果沿着时间轴把信号在时域上加以分段， 每一 段作傅里叶变换， 计算其频率， 则可以从各段频谱特性 随时间变化上看出信号的时变特性。这就是目前在时 变信号分析中广泛使用的分析方法 短时傅里叶变换， 因其处理过程简单、 高效， 得到了广泛的研究和使用。 短时傅里叶变换 ［21 ］的概念最早由 Gabor 于 1946 年提出的， 其计算公式为 STFT t， f∫ ∞ －∞ x τ h * τ － t e －j2πfτ dτ 5 式中 h * τ － t 是滑移时窗 h τ － t 的共轭， 滑移时窗 用来对信号进行分段截取。 3试验研究 3． 1试验件 图 3 和图 4 分别是复合材料螺栓连接组件的几何 尺寸和实物图， 由同种复合材料板通过复合材料螺栓 进行连接， 构成复合材料螺栓连接试验件。 图 3复合材料螺栓连接组件几何结构 mm Fig． 3Geometric structure of composite bolted connection mm 图 4复合材料螺栓搭接板 Fig． 4Bolted composite plate 3． 2试验系统 高温振动试验系统示意图和装置图如图 5 所示， 常温振动试验系统示意图如图6 a 所示。常温试验系 统主要包括振动控制系统和接触式振动测量系统， 高 温试验系统主要包括振动控制系统、 温度控制系统和 非接触式振动测量系统。其中 振动控制系统主要装置为电动振动台 型号 D- 600- 5; 频率范围 5 ～3 000 Hz; 额定加速度 100g ; 温度控制系统由 1 000 ℃ 高温环境箱 控制精度 1 ℃ 及双铂铑型 K 型 热电偶 测温精度 1． 5 ℃ 组成; 振动测量系统由德国 Micro- Epsilon 公司的 ODC- 1200- 10 光幕式激光位移传感器及基于 Labview 的非接 触式激光振动实时测量软件组成。针对常温振动试 验， 采用接触式测量方法， 如图 6 b 所示。其中传感 器为三向加速度传感器， 型号为 LC0122C- 50; 频率范 围 1 ～4 kHz; 量程 50g; 灵敏度 6． 05 pC/g。 a 高温振动试验系统 b 试验系统布置 1 － 温度控制系统; 2 － 振动测量系统; 3 － 振动控制系统 图 5高温振动试验系统和试验系统布置 Fig． 5Experimental system and arrangement of high- temperature vibration 夹具如图 6 b 左端所示， 该夹具利用 4 根不锈钢 丝杆配合夹块紧固， 夹块材料为 45 号钢。为保证夹持 力均匀一致， 对夹紧螺母施加相同大小的拧紧力矩。 3． 3试验流程 利用上述试验件及试验系统分别进行常温和高温 环境下的振动防松试验， 试验流程图如图 3 所示。常 温振动防松试验流程为 常温环境下， 首先对试件进行 扫频试验和分析， 得到分析结果 S1; 然后进行结构振动 试验 5 min 低量级， 3 min 高量级， 15 min 低量级 ， 并 在每次振动试验结束后对构件进行扫频试验和分析， 得到分析结果 S3， S4， S5。 702第 17 期曹芝腑等 基于时频分析的高温振动环境螺栓连接件松动判别 ChaoXing a 振动试验系统 b 传感器布置 图 6螺栓搭接板常温振动试验系统和传感器布置 Fig． 6Experimental system and sensor arrangement of bolted plate 与常温振动防松试验相比， 高温振动防松试验在 整个流程中会利用温度控制系统对试件进行加热， 调 控试件在振动过程中的环境温度。针对高温振动防松 试验， 增加高温和冷却至常温后的扫频试验步骤， 得到 分析结果 S2和 S6。结构振动试验 5 min 低量级， 3 min 高量级， 15 min 低量级 及后续的扫频试验均在高温环 境下进行， 即图 7 中虚线框部分。低量级和高量级的 载荷谱分别如图8和图9所示。 扫频试验中的扫频速 图 7常/高温振动试验流程图 Fig． 7Flow chart of high temperature vibrating experiment 图 8低量级载荷谱 Grms4． 7g Fig． 8Random spectrum of the low level 图 9高量级载荷谱 Grms 14． 48g Fig． 9Random spectrum of the high level 度为 0． 1 OCT/min， 扫频中振动台加速度值恒定为 0． 5g。 4试验结果与分析 按图 7 所示试验流程完成高温振动试验， 利用非 接触式振动测量系统获取响应数据。但由于工频成分 的存在， 影响试验结果的分析。本文利用一维数字滤 波方法， 滤除工频和倍频干扰对试验的影响 a1y N b1x N b2x N －1 bnb1x N － nb － a2y N －1 － － ana1y N － na 6 式中 x 为输入响应， y 为输出响应， ai和 bj分别为滤波 系数， na和 nb分别为反馈滤波器和前馈滤波器阶数。 为了实现对工频干扰信号的滤除， 采用二阶数字滤波 器来消除工频响应。滤波系数向量 a 和 b 分别为［ 1， －2cos ω0， 1］ 和［ 1，－2rcos ω0，r 2］ 。其中 ω 02πf0/ fs， f0为滤波频率， fs为采样频率， r 为常数， 取 0． 9。 图 10滤波前 S2时频分析结果 Fig． 10Time- frequency plot before filter of case S2 high- temperature 利用短时傅里叶变换时频分析技术对试验结果进 行分析， 得到高温振动环境下各试验流程中的时频图， 802振 动 与 冲 击2019 年第 38 卷 ChaoXing 如图 11 所示。图 11 中 Z 值最大的点 图 11 a ～ 图 11 f 中的 A 点 即为结构振动基频的振动振幅。然后 选取二倍频所对应的 B 点， 利用式 4 计算得到松动特 征指标计算结果 r Z ZB 2ωn /ZA ωn 。常高温环 境下， 扫频振动后计算得到的松动特征指标结果如图 12 所示。 a S1 b S2 c S3 d S4 e S5 f S6 图 11高温振动环境下各试验流程时频图 Fig． 11Time- frequency plot with all kinds of cases in high temperature vibrating 图 12常温和高温下各扫频试验后松动指标对比图 Fig． 12Comparison of loosening index with all kinds of normal and high temperature cases 1 图 11 a 对应于试验流程中常温环境下扫频 振动后的分析结果 S1。从图中可以看出， 二倍频成分 不明显。结合图 12 可知， 其初始松动判别指标为 6． 05 10 －3， 结构在常温状态下未出现连接松动。 2 随着温度的升高， 结构振动基频增加， 且螺栓 连接处接触状态发生改变， 导致图 11 b 呈现多频成 分， 松动特征值量级增加， 综合判定连接结构在 1 000 ℃高温环境中开始出现松动迹象。 3 图11 c ～ 图11 e 分别为1 000 ℃高温环境 下， 5 分钟低量级， 3 分钟高量级和 15 分钟低量级振动 后的扫频试验结果。从图中可以看出， 松动特征值在 3 分钟高量级振动后达到最大， 是常温状态下松动特征 值的 50． 9 倍。且图中二倍频成分均能清晰显现， 说明 连接件在高温振动环境下松动特征明显。 4 图 13 所示为试验前、 后螺栓连接局部对比 图， 在恢复到室温后， 连接件未见明显变形。图 11 f 为试验件恢复到常温后扫频试验分析结果 S6， 与试验 前比较， 二倍频成分明显。图 12 所示， 在高温环境下， S6对应的松动特征指标为 3． 37 10 －2， 为高温初始状 态的 15． 18 倍。而常温振动环境下， S5对应的松动指 标为 1． 3 10 －2， 为常温初始状态的 2． 15 倍。比较可 知， 高温振动环境下， 连接结构更容易产生松动。 a 试验前 b 试验后 图 13高温振动试验前后螺栓局部对比 Fig． 13Comparison of local bolt after high- temperature vibration 5结论 本文针对某型复合材料螺栓连接结构进行不同量 级下的高温振动防松试验， 利用时频分析方法提取振 动特征响应， 通过定义的高温松动判别指标 r Z ， 实 现了高温振动环境下复合材料螺栓连接结构的松动特 性判别及分析。 1 时频分析能够清晰表明复合材料螺栓连接件 在试验过程中的频率成分及变化， 为特征响应提取提 供准确的分析结果。 2 结合基频变化， 二倍频振幅和二倍频与基频 的振动幅值比值作为松动特性评判的特征量， 能够有 效评判复合材料螺栓连接结构在常/高温环境下的振 动松动特性。 3 高温对连接件结构材料属性的影响， 振动对 902第 17 期曹芝腑等 基于时频分析的高温振动环境螺栓连接件松动判别 ChaoXing 连接处接触状态的改变， 会导致连接结构在高温振动 后更容易出现连接组件的松动。 参 考 文 献 ［1］ 杜善义． 先进复合材料与航空航天［ J］ ． 复合材料学报， 2007， 24 1 1- 12． DU Shanyi．Advanced composite materials and aerospace engineering［ J］ ． Acta Materiae Compositae Sinica， 2007， 24 1 1- 12． ［2］ 姜东，吴邵庆，费庆国，等． 蜂窝夹层复合材料不确定性 参数识别方法［ J］ ． 振动与冲击， 2015， 34 2 14- 19． JIANG Dong，WU Shaoqing，FEI Qingguo，et al． Parameter identification approach of honeycomb sandwich composite with uncertainties［ J］ ． Journal of Vibration and Shock，2015，34 2 14- 19． ［3］ 雷宝灵，易茂中，徐惠娟． C/C 复合材料飞机刹车盘的三 维温度场［ J］ ． 复合材料学报， 2009， 26 1 113- 117． LEI Baoling，YI Maozhong，XU Huijuan． 3- D temperature field for C/C composite braking discs ［J］ ． Acta Materiae Compositae Sinica， 2009， 26 1 113- 117． ［4］ 刘冬欢，郑小平，王飞，等． 内置高温热管 C/C 复合材料 热防护结构热力耦合机制［J］ ． 复合材料学报，2010，27 3 43- 49． LIU Donghuan， ZHENG Xiaoping， WANG Fei， et al． Mechanism of thermomechanical coupling of high temperature heat pipe cooled C/C composite material thermal protection structure ［J］ ． Acta Materiae Compositae Sinica，2010，27 3 43- 49． ［5］ PRI J N，CALLOCH S，CLUZEL C，et al． Analysis of a multiaxial test on a C/C composite by using digital image correlationandadamagemodel ［J］ ．Experimental Mechanics， 2002， 42 3 318- 328． ［6］ CHOWDHURY N，CHIU W K，WANG J． Review on the fatigue of composite hybrid joints used in aircraft structures ［ J］ ． Advanced Materials Research Trans Tech Publications， 2014， 891 1591- 1596． ［7］ THOPPUL S D，FINEGAN J，GIBSON R F． Mechanics of mechanically fastened joints in polymer- matrix composite structures- a review［ J］ ． Composites Science and Technology， 2009， 69 3 301- 329． ［8］ 王志强． 复合材料层压板螺栓连接性能分析［D］ ． 哈尔 滨 哈尔滨工程大学， 2010． ［9］ 陈学前，杜强，冯加权． 螺栓连接非线性振动特性研究 ［ J］ ． 振动与冲击， 2009， 28 7 196- 198． CHEN Xueqian， DU Qiang， FENG Jiaquan．Nonlinear vibrational characteristicofbolt- joints ［J］ ．Journalof Vibration and Shock， 2009， 28 7 196- 198． ［ 10］ 郝淑英，陈予恕，张琪昌． 连接结构松动对系统非线性动 力学特性的影响［J］ ． 天津大学学报，2001，34 4 452- 454． HAO Shuying，CHEN Yushu，ZHANG Qichang． Effects of connect loosing and sliding on dynamic characteristics ［J］ ． Journal of Tianjin University， 2001， 34 4 452- 454． ［ 11］ 董广明，陈进，雷宣扬，等． 基于小波包- 神经网络方法的 支撑座连接螺栓松动损伤诊断的实验研究［J］ ． 机械科 学与技术， 2006， 25 1 102- 106． DONG Guangming， CHEN Jin， LEI Xuanyang， et al． Experimentalstudyondiagnosingtheattachmentbolt looseness in a clamping support based on wavelet packet transation and neural network ［J］ ． Mechanical Science and Technology， 2006， 25 1 102- 106． ［ 12］ PAI P F，PALAZOTTO A N． HHT- based nonlinear signal processingforparametricandnon- parametric identification of dynamical systems［ J］ ． International Journal of Mechanical Sciences， 2008， 50 12 1619- 1635． ［ 13］ LV J X，XIAO Y，XUE Y D． Time- temperature- dependent response and analysis of preload relaxation in bolted composite joints［J］ ．Journal of Reinforced Plastic and Composites， 2018， 37 7 460- 474． ［ 14］ 赵登峰，曾国英． 振动环境中螺纹联结松动过程的研究 ［ J］ ． 振动与冲击， 2010， 29 10 175- 178． ZHAO Dengfeng，ZENG Guoying． Unintentional unscrewing of screw joints under vibration condition ［J］ ．Journal of Vibration and Shock， 2010， 29 10 175- 178． ［ 15］ LI Q K，JING X J． A second- order output spectrum approach for fault detection of bolt loosening in a satellite- like structure with a sensor chain［ J］ ． Nonlinear Dynamics， 2017， 89 587- 606． ［ 16］ NGUYEN T C，HUYNH T C，YI J H，et al． Hybrid bolt- loosening detection in wind turbine tower structures by vibration and impedance responses［ J］ ． Wind and Structures， 2017， 24 4 385- 403． ［ 17］ AMERINI F，MEO M． Structural health monitoring of bolted joints using linear and nonlinear acoustic/ultrasound s ［ J］ ． Structural Health Monitoring， 2011， 10 6 659- 672． ［ 18］ KIM J Y，BALTAZAR A，HU J W，et al． Hysteretic linear and nonlinear acoustic response from pressed interfaces［J］ ． International Journal of Solids and Structures，2006，43 21 6436- 6452． ［ 19］ DRINKWATER B W，DWYER- JOYCE R S，CAWLEY P． A study of the interaction between ultrasound and a partially contacting solid- solid interface［ J］ ． Proceedings of the Royal Society of London A Mathematical， Physical and Engineering Sciences， TheRoyalSociety， 1996， 452 1955 2613- 2628． ［ 20］ AMERINI F，BARBIERI E，MEO M，et al．Detecting loosening/tightening of clamped structures using nonlinear vibration techniques［J］ ．Smart Materials and Structures， 2010， 19 8 085013- 1- 9． ［ 21］ GABOR D． Theory of communication． Part I the analysis of ination［J］ ．Journal of the Institution of Electrical Engineers- Part III Radio and Communication Engineering， 1946， 93 26 429- 441． 012振 动 与 冲 击2019 年第 38 卷 ChaoXing