2. 西安交通大学 机械结构强度与振动国家重点实验室,陕西 西安 710049
2. State Key Laboratory for Strength and Vibration of Mechanical Structures, Xi′an Jiaotong University, Xi'an 710049, China
裂纹不仅会引起结构局部刚度的变化,还会影响整个结构的机械特性,裂纹的出现与发展是整个机械结构故障的潜在危险.因此,研究裂纹结构的相关动力特性是十分必要的.在不损坏整体结构的前提下,可以通过裂纹结构的振动特性来进行诊断[1-2].Peng等[3]从频域角度提出了非线性输出频响函数 (NOFRFs) 的概念, 并将其用于裂纹结构的无损检测.此外,结构中裂纹的变化也会对阻尼产生一定的影响,即裂纹结构的阻尼特性可以对裂纹的严重程度进行标识和诊断.Panteliou等[4]从热力学的角度,研究了受交变轴向应力的单边裂纹梁的阻尼及阻尼比与系统振动的关系.Rezaee等[5]建立了一个单裂纹的简支梁模型,研究了系统阻尼比随着裂纹深度和裂纹位置的变化规律.Bovsunovsky等[6]建立了一种呼吸裂纹梁的有限元模型,研究了超谐共振下,阻尼对系统非线性响应的影响.裂纹的张开或呼吸会引起梁系统阻尼的增加,阻尼比是用于描述系统阻尼的一个特性参数.因此,可以通过系统阻尼比的变化来判断系统裂纹参数的变化,进而实现对裂纹结构的故障诊断.
本文基于文献[7]中模型的材料和几何参数,在ANSYS软件中建立了一种混合单元的直裂纹悬臂梁模型,并与文献[7]中实验测定的固有频率进行对比,验证了模型的有效性.在此模型的基础上,研究了冲击和正弦激励下,裂纹处界面摩擦对系统的影响;并且在摩擦可以忽略不计的情况下,分析了系统阻尼比和固有频率随裂纹深度及裂纹位置的变化规律.
1 直裂纹梁模型的建立本文研究了一个直裂纹悬臂梁,悬臂梁的长度L=300 mm,裂纹位置距离梁固定端的距离为d,直裂纹的深度为a,梁的截面尺寸为b×h=20 mm×20 mm.材料为结构钢12CrNi3A[8],弹性模量E=206.6 GPa,泊松比为0.3,密度为7 850 kg/m3,加载点如图 1b所示.考虑到全部使用平面单元建模会降低计算效率,因此,本文基于ANSYS软件,采用平面单元Plane183模拟裂纹区域,如图 1a所示,裂纹区的局部细化可以提高计算精度;采用梁单元Beam188模拟远离裂纹的区域,可以提高计算效率.在Plane183与Beam188的连接位置采用绑定接触来保证结构的整体性.需要注意的是裂纹区域宽度w及比例阻尼项的选取,本文通过对开裂纹梁混合模型和平面单元模型的一阶弯曲固有频率的对比,选择裂纹区宽度w=20 mm,二者频率吻合较好 (见表 1);瞬态计算设置比例阻尼项α=50, β=0[8].
呼吸裂纹的张开、闭合状态会导致裂纹梁的非线性现象.在一些实际应用中,系统会被认为是双线性的.当接近裂纹尖端的应力为正时,裂纹处于张开状态,否则,裂纹被视为闭合状态.无裂纹梁的固有频率为f1,开裂纹梁的固有频率为f2,呼吸裂纹梁的双线性固有频率f0为[9]
(1) |
当a=10 mm,d=140 mm时,将平面单元模型、混合单元模型与实验数据[7]的前3阶固有频率进行对比 (见表 1),结果表明平面单元模型、混合单元模型与实验数据吻合很好,这在一定程度上验证了混合单元模型的有效性.裂纹深度以无量纲参数s=a/h来表示,裂纹位置以p=d/L来表示.
3 裂纹悬臂梁系统阻尼特性分析本文讨论裂纹位置与裂纹深度对阻尼比的影响,通过施加2 kN的冲击力,冲击时间为0.02T0[10],T0=1/f0,得到加载点的y向位移自由衰减振动曲线,见图 2a,并通过式 (2) 和式 (3) 确定呼吸裂纹梁的阻尼比ζ:
(2) |
(3) |
式中:δ为对数衰减率;Ax为一选定振动波形峰值;Ax+n为推后n个周期后的波形峰值.
3.1 摩擦对裂纹梁振动响应的影响基于所建的裂纹梁模型,在表 2所示的两种工况下,其时域、频域和接触应力如图 2和图 3所示.其中,t为载荷作用时间.计算100个周期,每个周期划分60个载荷步.值得指出的是,图 2b中二阶固有频率两侧出现了调制均匀的边频.这是由裂纹导致的周期性脉冲冲击力,在频谱上体现为固有频率两侧出现调制均匀的边频[11].
图 2及图 3表明,针对裂纹悬臂梁,不论是正弦激励还是冲击激励,摩擦对其影响甚小,可以忽略不计.因此,后续的研究中都不考虑摩擦.图 2c和图 3c, 3d给出了无摩擦时,裂纹尖端、中间及末端的单元接触应力曲线及局部放大图 (单元位置见图 1a).
3.2 裂纹深度对阻尼比和固有频率的影响在p=0.27的情况下,研究s对阻尼比的影响.阻尼比的计算见式 (2) 和 (3)[12].其中双线性频率的计算见式 (1),呼吸频率是通过冲击加载获得的.图 4表明呼吸频率与公式计算所得的双线性频率误差很小,因此,双线性频率可以近似代替呼吸频率.对于冲击激励下的自由衰减振动,关于峰值点的选取目前还没有固定标准,本文对峰值点的选取进行了讨论 (见图 5).
首先,由于自由衰减振动初始时,振动波形不稳定,因此本文建议尽量推后几个周期开始选取峰值点,本文从第4个周期以后开始选取峰值点.图 5a表明,当选取相邻两个周期的峰值点 (第四周期和第六周期峰值点) 计算阻尼比时,阻尼比呈现较强的分散性,这是由于两个相近峰值相差较小,因此取点和计算值的误差较大,因此,适当地多间隔几个周期选取峰值点,计算得出的阻尼比的点集中度较好,如图 5b,5c所示.本文得到阻尼比的一些数据点,为了更好地预测阻尼比与裂纹深度的变化趋势,可以采用最小二乘多项式拟合得到阻尼比数据点的拟合曲线.
进行数据拟合时需要保证两点:①误差平方和越接近0,拟合曲线与数据点的误差越小;②确定系数越接近1,拟合度越好.据此,采用最小二乘法得到的三阶多项式来进行曲线拟合,拟合公式为y=a0+a1x+a2x2+a3x3,拟合曲线见图 5.图 5表明,随着裂纹深度的增加,系统阻尼呈现上升趋势,其变化情况与文献[4]的实验结果类似.
3.3 裂纹位置对阻尼比和固有频率的影响图 6给出了s=0.5时,裂纹梁的双线性频率和呼吸频率随裂纹位置p的变化曲线.从图中可以看出,双线性频率和呼吸频率相差不大.这再一次证明了用计算得到的双线性频率代替呼吸频率的合理性.阻尼比随裂纹位置的变化规律仍采用对数衰减率的方式计算,峰值点以及拟合曲线的选取见3.2所述.本节中,拟合曲线采用四阶多项式:y=a0+a1x+a2x2+a3x3+a4x4.图 7再一次表明适当地多选几个周期计算阻尼比,其点的分布更集中,效果更好.且随着裂纹位置的偏移,阻尼比呈现先增加后减小的趋势.这是由于:无量纲裂纹位置为0时,代表无裂纹健康梁的情况,随着裂纹的产生,其阻尼比增加,然后,随着裂纹位置越来越远离固定端,裂纹梁的状态越来越接近于无裂纹的情况,阻尼比也会随之减小.
1) 无论是冲击激励还是正弦激励下,摩擦对系统振动响应的影响很小,研究中可以忽略不计.通过冲击激励确定的系统固有频率和计算得到的双线性固有频率吻合很好.
2) 随着裂纹严重程度的增加,系统的阻尼比呈现出递增趋势.利用呼吸裂纹梁的自由衰减振动曲线进行阻尼比的计算时,适当地增加衰减的周期数可以更好地计算系统阻尼比.采用多项式曲线拟合数据,可以更好地描述系统阻尼比随裂纹深度和位置的变化规律.
[1] | Chati M, Rand R, Mukherjee S J. Modal analysis of a cracked beam[J]. Journal of Sound and Vibration, 1997, 207(2): 249–270. DOI:10.1006/jsvi.1997.1099 |
[2] | Andreaus U, Baragatti P. Cracked beam identification by numerically analyzing the nonlinear behavior of the harmonically forced response[J]. Journal of Sound and Vibration, 2011, 330(4): 721–742. DOI:10.1016/j.jsv.2010.08.032 |
[3] | Peng Z K, Lang Z Q, Billings S A. Crack detection using nonlinear output frequency response functions[J]. Journal of Sound and Vibration, 2007, 301(3): 777–788. |
[4] | Panteliou S D, Chondros T G, Argyrakis V C. Damping factor as an indicator of crack severity[J]. Journal of Sound and Vibration, 2001, 241(2): 235–245. DOI:10.1006/jsvi.2000.3299 |
[5] | Rezaee M, Hassannejad R. Free vibration analysis of simply supported beam with breathing crack using perturbation method[J]. Acta Mechanica Solida Sinica, 2010, 23(5): 459–470. DOI:10.1016/S0894-9166(10)60048-1 |
[6] | Bovsunovsky A P, Surace C. Considerations regarding superharmonic vibrations of a cracked beam and the variation in damping caused by the presence of the crack[J]. Journal of Sound and Vibration, 2005, 288(4/5): 865–886. |
[7] | Rizos P F, Aspragathos N, Dimarogonas A D. Identification of crack location and magnitude in a cantilever beam from the vibration modes[J]. Journal of Sound and Vibration, 1990, 138(3): 381–388. DOI:10.1016/0022-460X(90)90593-O |
[8] |
崔韦, 王建军.
共振条件下的裂纹梁振动与裂纹扩展耦合分析[J]. 推进技术, 2014, 35(10): 1404–1411.
( Cui Wei, Wang Jian-jun. Coupling analysis of vibration and crack propagation for a cracked beam at reesonant state[J]. Journal of Propulsion Technology, 2014, 35(10): 1404–1411. ) |
[9] | Chu Y C, Shen M H H. Analysis of forced bilinear oscillators and the application to cracked beam dynamics[J]. Journal of American Institute of Aeronautics and Astronautics, 1992, 30(10): 2512–2519. DOI:10.2514/3.11254 |
[10] | Andreaus U, Casini P, Vestroni F. Non-linear dynamics of a cracked cantilever beam under harmonic excitation[J]. International Journal of Non-linear Mechanics, 2007, 42(3): 566–575. DOI:10.1016/j.ijnonlinmec.2006.08.007 |
[11] | Cheng S M, Wu X J, Wallace W. Vibration response of a beam with a breathing crack[J]. Journal of Sound and Vibration, 1999, 225(1): 201–208. DOI:10.1006/jsvi.1999.2275 |
[12] | Andreaus U, Baragatti P. Fatigue crack growth, free vibrations and breathing crack detection of aluminum alloy and steel beams[J]. The Journal of Strain Analysis for Engineering Design, 2009, 44(7): 595–608. DOI:10.1243/03093247JSA527 |