疲劳裂纹闭合的非线性超声在线测量.DOC

1、1疲劳裂纹闭合的非线性超声在线测量 贾俊 1，陶利民 2，胡海峰 1，胡永攀 3，吕伟 1（1. 国防科学技术大学 装备综合保障技术重点实验室，湖南 长沙 410073；2.钦州学院 机械与船舶海洋工程学院，广西 钦州 535011;3.后勤工程学院 国家救灾应急装备工程技术研究中心，重庆 沙坪坝 401331）摘要：建立了外部载荷作用下用非线性振动-声场调制指数表征裂纹尺寸的数学模型，分析了外载荷与塑性致闭应力共同作用下裂纹面的开闭响应，提出了裂纹张开和闭合应力的在线测量方法。在一块含中心穿透裂纹的铝合金薄板试件上进行了裂纹闭合参数的测量实验，实验结果表明使用振动-声场调制技术在线测量裂纹闭

2、合参数是可行且有效的，且该方法还能用于分析小裂纹的特殊开闭行为。关键词：裂纹闭合；张开应力；振动-声场调制；结构健康监测；非线性超声检测中图分类号：TP277 文献标志码：A 文章编号： On-line fatigue crack closure determination using nonlinear ultrasound testingJIA Jun1, TAO Limin2, HU Haifeng1， HU Yongpan3，LV Wei 1(1. Science and Technology on Integrated Logistics Support Laboratory，Nat

3、ional University of Defense Technology, Changsha 410073, China；2.College of Mechanical and Marine Engineering, Qinzhou University, Qinzhou 535011，China ；3. National Engineering Research Center for Disaster and Emergency Rescue Equipment, Logistical Engineering University, Chongqing, 401331, China) A

4、bstract: A mathematical model which reveals the relation between the vibro-acoustic modulation index and the crack contact area is developed. Basing on the analysis of crack faces response to the resultant of applied load and plasticity induced crack closure stress, a procedure of on-line approximat

5、ing the opening and closure stress is proposed. To study the performance of this measuring method experimentally, on-line vibro-acoustic modulation tests were carried out on a central cracked aluminum alloy fatigue sample. The experimental estimations match well with other researches. Analysis of th

6、e experimental results suggests that the vibro-acoustic modulation technology is not only capable of measuring the fatigue crack opening and closure stresses but also capable of analyzing the crack closure behavior, even for physically short cracks. Keywords: crack closure; opening stress; vibro-aco

7、ustic modulation; structural health monitoring; nonlinear acoustic testing收稿日期：基金项目：国家自然科学基金资助项目（51305449）作者简介：贾俊（1985），男，湖南常德人，博士研究生，E-mail：junjia_；陶利民（通信作者），男，研究员，博士，博士生导师，E-mail：受损结构在完全卸载之前（即在某一大于零的拉伸载荷下） ，疲劳裂纹面出现部分接触的现象，称为裂纹闭合。Elber 1于 1971 年首次在平面应力试件拉-拉疲劳裂纹扩展试验中观察到裂纹闭合现象。裂纹闭合，改变了人们对疲劳裂纹扩展真正动因的认

8、识，被认为是影响疲劳裂纹扩展特性的关键因素。疲劳裂纹的张开和闭合应力是裂纹闭合理论的重要参数。对其准确地在线测量，是在线预测构件剩余寿命的前提。张开应力和闭合应力近似相等且测量方法类似，因此本文之后如不特殊说明，均只讨论张开应力。测定张开应力的方法很多，如柔度法 2、电位法 3及超声波方法 4, 5等等。其中最可靠且应用最广的是利用裂纹张开位移测量闭合应力的柔度法。但该方法是一种离线检测方法，需要在裂纹附近位置2安装引伸计，因此并不适合在线应用。超声方法实施方便，灵敏度较高，是一种极具潜力的无损检测方法。相对柔度法，由于其能对目标构件进行全域检测，因此对损伤的先验知识需求较少；相对电位法，由于

9、其无需局限于导电材料，因此适应性较强。超声检测方法的这些优良特性，使其适用于张开应力的在线测量。Rokhlin 等 4使用一次反射和板底反射波幅度为表征，实现了 2024-T3 铝合金板内裂纹张开应力的在线测量，并观察到了轻微的滞回效应。Singh 等 5使用 Lamb 波回波脉冲能量为表征，实现了 7075-T6 铝合金板内穿透裂纹张开应力的在线测量，并将测量结果与Newman 的数值计算结果进行了对比验证，研究结果表明使用超声方法用于裂纹张开应力的在线测量及构件剩余寿命的在线评估是可行的。基于反射回波检测等传统线性超声方法的缺点是灵敏度相对较低，且易受孔洞等几何构形干扰。振动-声场调制（V

10、ibro-Acoustic Modulation， VAM）是一种新颖的非线性超声方法，其使用高低频两种声场同时激励受损构件，并以散射声场的频域畸变来表征损伤程度。相对线性和其他非线性检测方法，其抗干扰能力强，对微小裂纹的灵敏度高，适用于裂纹的早期探测 6。大量应用研究证实了 VAM 技术在损伤探测和扩展跟踪方面的先进性和可靠性。Klepka 等 7使用 VAM 技术实现了对复合材料冲击损伤的检测，并分析了激励声场及试样振动模态对调制强度的影响。Lim 等 8基于 VAM方法使用空气耦合换能器（Air-Coupled Transduces, ACTs）实现了铝板疲劳裂纹的非接触式探测与成像。G

11、ian 等 9分析了裂纹扩展时VAM 调制指数与试样剩余寿命的对应关系，并实验验证了用 VAM 方法估测结构剩余寿命的可行性。国内的胡海峰等 10和刘斌等 11分别使用希尔伯特变换和同步解调等时频分析方法提出了改进的 VAM 损伤表征参量。然而，使用 VAM 技术测量疲劳裂纹张开应力还未见公开报道。1 理论与方法1.1 外载荷作用下的 VAM 调制指数如图 1 所示，考虑一长度为 ，厚度为L受损薄板，在远端疲劳载荷 的作()HL用下板内出现一穿透型裂纹，裂纹面垂直于轴，位置为 。设有两列沿 轴正方向x(,)cxzx传播的应力波，其中低频振动应力波为弯曲波，在杆内激发弯曲振动模态，频率为 。另一

12、激1f励波为高频超声纵波，频率为 。根据 VAM2f理论，裂纹等接触类损伤将会导致散射声场发生非线性畸变，即出现频率为 的和差频21f分量。z xH0cxczL1f2f,Ub ,U 裂纹图 1 薄板内裂纹的散射声场Fig.1 An illustration of the ultrasonic field scattered by a fatigue crack in a thin plateDonskoy 等 12基于裂纹的非线性弹簧模型，给出了构件自由状态（即 ）下散射声场0调制边带幅值为：* 20101020sin()cos()coABUzkDSkLLxxMERGEFORMAT (1)式中

13、， 为裂纹的等效非线性弹簧系数； 为0D3裂纹面的平均距离（或裂纹体的等效厚度） ；为裂纹面的实际接触面积；S 为薄板的横截cs面积；A 和 B 分别为两入射波的幅度； 、1k和 分别为两入射波和调制边带分量的波数。2kDonskoy 等 12还进一步推证了 VAM 调制指数（Modulation Index, MI）与裂纹接触面积 的cs正比关系：* 12MIcAsMERGEFORMAT (2)式中： 、 和 分别为两入射波和调制边1A带分量的频谱幅值。MI 作为损伤指数评价结构的损伤严重程度已得到广泛应用 9, 13。然而，当考虑构件外部载荷及疲劳裂纹闭合应力的影响时，裂纹面微凸体的接触状

14、态和非线性特性将发生改变，MI 也将随之变化。此时 MI 的推导应作适当修正。 Drr图 2 受载裂纹体的等效厚度Fig.2 The effective thickness of a loaded crack首先，如图 2 所示，设在外部拉应力 与裂纹闭合应力 的合成应力场 的作用下裂纹r接触面的平均距离由 变为 。此时接触面非0D线性弹簧模型可表示为：* 201020()()KMERGEFORMAT (3)式中： 为自由状态下的外部应力； 为自由0 0D状态下裂纹接触面的平均距离； 和 分别为1K2接触面一阶（线性）刚度与二阶（非线性）刚度，且有 14：* 21100200(),()mmKD

15、CKDCMERGEFORMAT (4)式中：C 为与粗糙度等表面参数有关的常数。Biwa 等 14指出，对微凸体高度完全随机分布的裂纹面可取 ，而对具有相同高度的微1凸体表面 m 可取接近于零的小值。考虑到疲劳裂纹扩展过程中因裂纹面张开闭合而导致的“撕裂”和“压平”效应，现取 。结合0.5m式* MERGEFORMAT (3)与式* MERGEFORMAT (4)，并考虑到 ，可解0得考虑外部应力时裂纹面的平均距离为：* 02DCMERGEFORMAT (5)其次，考虑外部应力后，裂纹面的“应力突变”边界条件应修正为：* 32ccxxFSMERGEFORMAT (6)式中： 表示和差频散射场在

16、裂纹处的应力幅度； 为因裂纹开闭运动产生的频率为和差频3F的总压力；S 为薄板的横截面积。经过以上两处修正后，可求得和差频分量的反射和透射场非线性系数（损伤指数）分别为：* , 20,2MI()bffc rfrsaDCAB MERGEFORMAT (7)式中：* 222112(1)cos()insicoi4b cf ccvkLxaECkxkbSMERGEFORMAT (8)上两式中，上标 和 分别表示后向和前向散bf4射场；除 和 外，其它所有符号意义与文献csk12相同。1.2 疲劳裂纹的开闭行为及张开应力的估计裂纹的闭合效应，主要表现为裂纹面存在分布的闭合压力，且在外载荷与闭合压力共同作用

17、下裂纹出现部分张开和部分闭合。目前认为在平面应力状态时塑性致闭（Plasticity Induced Crack Closure, PICC）效应是导致闭合效应的主要原因 15。PICC 理论认为产生闭合压力的原因是裂纹尖端附近材料的塑性变形,随着裂纹的扩展这种塑性变形将会形成裂纹扩展尾迹，导致闭合压力的产生。 xy张开 闭合闭合裂纹闭合压力图 3 裂纹在外载荷和 PICC 应力作用下部分张开Fig.3 A partially open crack subjecting to the resultant of tensil load and PICC closure stress 图 3 所示

18、为一疲劳裂纹在外部拉应力和 PICC 闭合应力的作用下部分张开。当拉应力从零逐渐增大时，部分闭合压应力被抵消，相应部分裂纹面张开，剩余部分仍闭合。如图4 所示，当拉应力增大到等于裂纹张开应力时，裂纹面完全张开。此后随着拉应力继续增大，裂纹开始扩展。若不考虑滞回效应，卸载时裂纹面的响应与加载时类似。M I0Ics部分张开完全张开加载卸载图 4 用调制指数变化曲线估计张开应力示例Fig.4 An example of opening stress approximating根据第 1.1 节的讨论，调制指数 MI 与裂纹的实际接触面积成正比，因此裂纹的张开部分对 MI 并无贡献，即可用 MI 跟踪

19、裂纹闭合面积的变化。当裂纹在拉应力作用下完全张开时，闭合面积为零，此时 MI 将减小至某一接近于零的值，且随后将不再随外部应力变化。因此，如图 4 所示，可用 MI 曲线的拐点对应的外部载荷 作为张开应力 的估计值，即 。opop闭合应力 的估计可按同样方法在裂纹卸载过cl程得到。需要说明的是，当裂纹完全张开后，MI 并不为零的原因是材料晶格非线性及试件夹持边界非线性的影响 16。3 实验验证3.1 疲劳试件与载荷序列疲劳试件为一块尺寸为 的34012m6061-T6 铝合金薄板。板中部开有一长度为10mm 的裂纹扩展槽。为避免试件两端的加持部分产生压痕及减小夹持边界非线性效应 16，试件两端

20、粘有四块保护垫片，见图 5。5C H 2C H 3C H 4C H 1C H 2C H 1C h a n e l 1 C h a n e l 2tyA m p l i t u d e : 2 VF r e q u e n c y :3 K H zA m p l i t u d e : 1 0 VF r e q u e n c y :5 0 K H zAAA - A疲劳试件压电放大器滤波放大器 示波器信号发生器计算机图例 ：加强垫片 低频压电陶瓷片 压电传感器S 1V AU ES 2高频压电陶瓷片扩展槽图 5 疲劳试样与 VAM 实验装置Fig.5 Fatigue sample and VAM

21、testing setup时间0. . .载荷序列 1V R T准备阶段裂纹长度测量. . . . .序列 5V S R疲劳V A M疲劳V A Mmaxmin裂纹长度测量图 6 疲劳载荷谱Fig.6 Fatigue load spectrum疲劳试验在 MTS 810 材料性能试验机上进行。如图 6 所示，疲劳载荷分为一个准备阶段和五个加载序列，均由一个拉- 拉周期疲劳载荷段及一个准静态载荷段组成。准备阶段目的是系统初始化、减少板内残余应力以及确定低频激励信号频率，主要进行振动去应力（Vibratory Stress Relief, VSR）和振动响应测试（Vibration Respons

22、e Tests, VRT） 。五个加载序列主要进行裂纹的疲劳扩展及 VAM 实验，且在每个序列的周期疲劳段结束后测量裂纹的长度。除 VRT 外，其余准静态载荷段均由一个加载过程和一个卸载过程组成，分别用于张开和闭合应力的测量。为防止薄板试样发生压缩失稳，整个实验过程只施加拉伸载荷，疲劳载荷的最大值 为 70MPa，应力比 R 为 0.1；max静态加载时最小值为 0MPa，最大值为70MPa，步进值 10MPa。3.2 VAM 实验为保证在 VAM 试验时裂纹面在不同的静态载荷下，其振动模态和幅值尽量相同，在准备阶段进行了 VRT 实验（即不同载荷下裂纹附近处的模态测试实验） ，用于选择 VA

23、M 的低频激励频率，测试结果如表 1 所示。模态选取的原则是选择模态振幅在载荷变程中的相对稳定者。表 1 VAM 实验低频激励信号频率Tab.1 Pumping signal frequence for VAM testing静态载荷 (MPa)模态频率(Hz)静态载荷 (MPa)模态频率(Hz)0 3110 40 312410 3113 50 312720 3117 60 313030 3122 70 3135VAM 实验装置如图所示。4 个压电陶瓷片6粘贴于试件的右半部分，VA 用于低频激励；S2 靠近裂纹萌生处，作为 VRT 测量传感器；S1 用于接收散射声场；UE 为高频激励元件。VA

24、M 实验的高频激励信号频率为 50kHz，低频激励信号频率按表 1 设置，且两激励信号均为幅值为 20 的正弦波p-V3.3 实验结果与分析3.3.1 张开应力测量图 7 所示为裂纹长度为 2.46mm-5.56mm 时，VAM 调制指数 MI 随外部拉应力 变化的归一化测量结果。因当裂纹较小（1.06mm）时，其开闭特性差异较大，故其测量结果将另述于3.3.2 小节。图 7 表明，MI 随拉应力变化的趋势与图 4 是总体一致的，且根据 1.2 节中对 MI曲线拐点的讨论，张开应力和闭合应力的估计结果在 30MPa 到 40MPa 之间。若缩小外部拉应力的采样间隔，则估计精度还可进一步提高。T

25、ong 和 Wu17使用一种改进的有限元分析方法给出了裂纹张开应力的估计式：* 23max01max(),.40opyARA且MERGEFORMAT (9)其中：* 1max20 max120133(.9580.416.8)cos().)yyAAA MERGEFORMAT (10)式中， 为约束因子，平面应力状态时取 ，=1平面应变状态时取 ； 为材料的拉伸强度=3y与屈服强度的平均值。对本实验，使用式（9）与式（10）计算张开应力的结果为 39.2MPa，与测量结果吻合。图 7 还表明，裂纹长度为 2.46mm-5.56mm时，张开应力测量值均在同一范围，表明张开应力在裂纹扩展过程中较为稳定

26、。该结论同样与其它研究结果吻合 1870 20 40 6000.20.40.60.81一一一一一一一 MPa0 20 40 6000.20.40.60.81一一一一一一一 MPa0 20 40 6000.20.40.60.81一一一一一一一 MPa0 20 40 6000.20.40.60.81一一一一一一一 MPa一一一一一一一一一一一一一一一一(c) c= 4.48mm (d) c= 5.56mm(a) c= 2.46mm (b) c= 3.52mm图 7 裂纹较长时调制指数的归一化测量值Fig.7 Normalized MI against applied load of long cr

27、acks 3.3.2 小裂纹开闭特性图 8 所示为当裂纹长度为 1.06mm 时，MI随拉应力的变化情况。对比图 7 和图 8 可以看出，小裂纹表现出与大裂纹明显不同的开闭特性。MI 开始随拉应力增大而增大，当拉应力超过 40MPa 后，MI 随拉应力的增大而减小，且直到最大拉应力时也未出现拐点。这表明当外部拉应力低于 40MPa 时，裂纹一直处于完全闭合状态，此后裂纹面在拉应力作用下开始部分张开，且一直到最大拉应力时裂纹面始终未完全张开，即当裂纹长度为 1.06mm 时，张开应力的估计值不小于 70MPa。0 20 40 6000.20.40.60.81一一一一一一一 MPa一一一一c=1.

28、06mm图 8 裂纹较小时调制指数的归一化测量值Fig.8 Normalized MI against applied load of a short crack另外，当裂纹面完全闭合（即 ）40MPa时，MI 随外部拉应力 的增大而增大，这与式（7）吻合，进一步验证了受载裂纹 VAM 模型的正确性。小裂纹在拉应力作用下仍然完全闭合的主要原因，是当裂纹开始在扩展槽根部萌生时，材料塑性变形较为剧烈，应力集中现象明显，PICC 闭合应力较大 19，因此当 较小时不足以克服 PICC 闭合应力而使裂纹面张开。4 结论根据 VAM 调制指数对裂纹接触面积敏感的基本特性，提出了使用 VAM 技术在线测量

29、裂纹闭合参数的方法。该方法通过 VAM 调制指数变化曲线拐点对应的外部载荷来估计裂纹的张开和闭合应力。6061-T6 铝合金板上裂纹张开应力的测量数据表明，当裂纹相对较大时，张开应力的估计值介于 30 Mpa 到 40Mpa 之间，且随着裂纹的扩展而较为稳定，实验结果与理论计算值吻合。当裂纹较小时，VAM 调制指8数能正确反映裂纹的闭合状态。因此，VAM技术不仅能较为准确地在线测量疲劳裂纹的张开和闭合应力，还能用于分析小裂纹的特殊闭合行为。参考文献（References ）1 Elber W. The significance of fatigue crack closureJ. Damage

30、 Tolerance in Aircraft Structures, 1971, 486: 230-242. 2 Song J H, Chung Y I. A review of crack closure measurement by compliance technique and the normalized-extended ASTM method as a currently most refined, practical and simple oneJ. Procedia Engineering, 2010, 2(1): 777-786. 3 Andersson M, Persso

31、n C, Melin S. Experimental and numerical investigation of crack closure measurements with electrical potential drop techniqueJ. International Journal of Fatigue, 2006, 28(9): 1059-1068. 4 Rokhlin S I, Kim J Y. In situ ultrasonic measurement of crack closureJ. International Journal of Fatigue, 2003,

32、25(1): 51-58. 5 Singh D S, Srivastav A, Gupta S, et alUltrasonic measurement of crack opening load for life-extending control of mechanical structuresC/ American Control Conference, 2009: 210-215.6 Zagrai A, Donskoy D. Micro/meso scale fatigue damage accumulation monitoring using nonlinear acoustic

33、vibro-modulation measurementsC/ Nondestructive Evaulation for Health Monitoring and Diagnostics, 2006: 617506-617512.7 Klepka A, Pieczonka L, Staszewski W J, et al. Impact damage detection in laminated composites by non-linear vibro-acoustic wave modulationsJ. Composites Part B Engineering, 2014, 65

34、(10): 99-108. 8 Lim H J, Song B, Park B, et al. Noncontact fatigue crack visualization using nonlinear ultrasonic modulationJ. NDT&E International, 2015, 73: 8-14. 9 Fierro G P M, Meo M. Residual fatigue life estimation using a nonlinear ultrasound modulation methodJ. Smart Materials and Structures,

35、 2015, 24(2): 1473-1484. 10 胡海峰，胡茑庆，秦国军非线性振动声调制信号耦合特征分析J机械工程学报， 2010，46(23)：68-76.HU Haifeng, HU Niaoqing, QIN Guojun. Coupling analysis of nonlinear vibro-acoustic modulation signalsJ. Journal of Mechanical Engineering, 201046(23)68-76. (in Chinese) 11 刘斌，刚铁，万楚豪，等金属杆内疲劳裂纹与振动/超声相互作用及其定量表征J声学学报，2016(

36、4)：507-514.LIU Bin, GANG Tie, WAN Chuhao, et al. Interaction of fatigue crack with vibration and ultrasound in metallic rod and its use for quantitative characterizationJ. ACTA ACUSTICA, 2016(4)507-514.12 Donskoy D, Sutin A, Ekimov A. Nonlinear acoustic interaction on contact interfaces and its use

37、for nondestructive testingJ. NDT&E International, 2001, 34(4): 231-238.13 Li N, Sun J J, Jiao J P, et al. Quantitative evaluation of micro-cracks using nonlinear ultrasonic modulation methodJ. NDT&E International, 2016, 79: 63-72. 14 Biwa S, Nakajima S, Ohno N. On the acoustic nonlinearity of solid-

38、solid contact with pressure-dependent interface stiffnessJ. Journal of Applied Mechanics, 2004, 71(4): 508-515. 15 Newman J C，Elber W. Mechanics of fatigue crack closureM. Philadelphia: ASTM, 1988.16 Aymerich F, Staszewski W J. Impact damage detection in composite laminates using nonlinear acoustics

39、J. Composites Part A Applied Science & Manufacturing, 2010, 41(41): 1084-1092. 17 Tong D H, Wu X R. Analysis of crack opening stresses for center- and edge-crack tension specimensJ. Chinese Journal of Aeronautics, 2014, 27(2): 291-298. 18 Antunes F V, Chegini A G, Correia L, et al. Numerical study of contact forces for crack closure analysisJ. International Journal of Solids & Structures, 2014, 51(6): 13301339. 19 Williams J F, Stouffer D C. An estimate of the residual stress distribution in the vicinity of a 9propagating fatigue crackJ. Engineering Fracture Mechanics, 1979, 11(3): 547-557.