Study on Ultrasonic Guided Wave Propagation Characteristics and Damage Imaging for Composite Structures Under Variable Temperature Field
-
摘要: 基于Lamb波的损伤监测方法已在复合材料结构健康监测中得到了广泛的应用。然而，复合材料结构的服役环境复杂多变，超声导波信号极易被外界因素干扰。本文研究了玻璃纤维层合板在变温场下的超声导波特性及损伤成像修正方法。首先，利用红外热成像仪研究持续加热下超声导波信号幅值和相位的变化规律，并构建因变温引起的幅相误差矩阵。其次，将误差矩阵代入多重信号分类（2D-MUSIC）算法中修正导向矢量，并建立代价函数迭代出准确损伤位置。玻璃纤维复合材料层合板的模拟损伤实验表明该方法有效提高了2D-MUSIC算法在变温场下的损伤定位分辨率和精度。
-
关键词:
- 变温场 /
- 超声导波 /
- 2D-MUSIC算法 /
- 损伤成像
Abstract: The Lamb-wave-based damage location method has been widely used for health monitoring of composite structures. However, it is easily interfered by external factors because its service environment is complex and changeable. To determine the influence of the temperature field on the propagation of Lamb waves in composite structures, in this study we first investigated the propagation characteristics of ultrasonic guided waves on glass fiber laminates under a temperature field using an infrared thermal imager. Subsequently, a corrected model was established using the extracted amplitude attenuation and phase delay errors. Consequently, a modified multiple signal classification(2D-MUSIC) algorithm-based damage imaging method is proposed for composite structures. The experimental results on glass fiber composite laminates show that the proposed method can effectively improve the resolution and accuracy of the original algorithm under a variable temperature field. -
0. 引言
回热器是回热式低温制冷机中的关键部件，冷热流体工质在其中交变流动，直接与回热器填料进行热交换。根据计算回热器1%的效率损失会引起接近10%的理论制冷量损失，当回热器效率为90%时，理论制冷量损失会达到100%[1]，由此可见回热器效率的重要性。因此对回热器的研究是制冷机研究中的一大重点。根据填料填充方式的不同，回热器可分为径向填充（层叠金属丝网式，堆叠金属珠式，金属丝网和金属珠混合式），轴向填充（平行金属丝式，平行金属片式），随机填充3种类型[2]，以上3种填充结构都为多孔介质结构。当前对回热器的研究方法主要有实验研究和数值模拟，不管是实验研究方法还是数值模拟方法都未能从微观结构上对回热器进行研究，因此本文采用分形理论对回热器进行表征。
1. 分形表征理论基础
多孔介质是指由许多骨架形成大量微小缝隙的物质，丝网填充的回热器是一种典型的多孔介质。回热器中的丝网一般为100目~500目，将丝网裁剪后填入回热器外管，形成回热器的基本结构。
图 1为不锈钢网片堆叠的示意图，图 2为不锈钢网片和不锈钢毡的显微照片，除图 1和图 2所示的径向填充回热器填充结构，还有轴向填充、随机填充等填充结构的回热器，这些回热器都是多孔结构。一般实际多孔介质都具有统计自相似的分形特征。已有文献报道[3-5]，多孔介质的孔隙大小分布满足分形标度律。多孔介质中孔隙直径大于或等于λ的孔隙数目和大小满足如下分形标度关系[6]：
![]()
图 2 不锈钢网片（左）和不锈钢毡（右）Figure 2. Stainless steel mesh(left) & stainless steel felt(right)$$N\left( {L \geqslant \lambda } \right) = {\left( {{\lambda _{{\rm{max}}}}/\lambda } \right)^{{D_{\rm{f}}}}}$$ (1) 式中：N为孔隙数目；L为孔隙直径；Df为分形维数；λmax为最大孔隙直径。当λ取最小值λmin时，即可得孔隙总数目为：
$${N_{\rm{t}}} = {({\lambda _{\max }}/{\lambda _{\min }})^{{D_{\rm{f}}}}}$$ (2) 需要注意的是自然界和人造材料的分形结构一般是统计自相似的，因此上式的适用范围会在一定的尺度区间内。λmin~λmax内分形维数相同，超出这个区间分形维数可能不同也可能不存在。式(2)即为多孔介质中孔隙数目的分布。
对于精确自相似分形，其分形维数可以直接通过计算得出[7-8]。对于统计自相似分形，因其并不是由数学迭代生成，无法知道其具体结构。通常用其他方法确定分形维数。
实验方法依据分形标度关系得到分形维数，Yu[9]使用分形分析推导得出统计自相似分形多孔介质的分形维数解析表达式：
$$\mathit{\Phi } = {\left( {{\lambda _{{\rm{max}}}}/{\lambda _{{\rm{min}}}}} \right)^{2 - {D_{\rm{f}}}}}$$ (3) 式中：Φ为多孔介质的孔隙率，Φ=孔隙体积/总体积。
式(3)建立了多孔介质孔隙率，最大最小孔隙和分形维数的关系。该式中的2指的是结构拓扑维数为2，将该式的维度推广，可得：
$${D_{\rm{f}}} = {d_{\rm{E}}} - \frac{{{\rm{ln}}\mathit{\Phi }}}{{{\rm{ln}}\left( {{\lambda _{{\rm{max}}}}/{\lambda _{{\rm{min}}}}} \right)}}$$ (4) 式中：三维空间dE为3，二维空间dE为2，同时也表明若二维空间与三维空间孔隙率相同，最大孔径与最小孔径比也相等，则分形维数相差1。在实际中可通过测定孔隙率和最大最小孔径来得到结构的分形维数。如采用气体吸附法和压汞法[10]。
孔隙率是多孔介质结构固有的参数，对于多孔介质的孔隙，通常可以用实验方法得出。显微分析法是一种粗略的测量分析方法，它用二维平面测量结果近似代替三维空间指标。在这种孔隙率测试方法中，先用显微镜观察出多孔样品截面的总面积和其中包含的孔隙面积，然后直接通过两者比值计算出多孔体的孔隙率。除此之外还有质量-体积计算法，浸泡介质法等[10]。
在分形分析中，可以推导得出孔隙率与微观结构之间关系的表达式[11-12]：
$${L_0} = {\left[ {\frac{{1 - \mathit{\Phi } }}{\mathit{\Phi } }\frac{{{\rm{ \mathit{ π} }}{D_{\rm{f}}}{\lambda _{{\rm{max}}}}^2}}{{4\left( {2 - {D_{\rm{f}}}} \right)}}} \right]^{1/2}}$$ (5) 方程建立了分形多孔介质分形维数、孔隙率、介质微结构之间的定量关系式。其中L0指研究对象的边长。
基于现有的实验材料与实验条件，拟采用压汞法对回热器孔隙分布进行实验。压汞法的原理为非浸润液体注入多孔介质时，表面张力阻止液体进入孔隙，需额外施加压力才能将液体注入多孔介质孔隙中，施加的压力等于毛细压力[13]。由于不同孔隙大小产生的毛细压力不同，因此可以建立外加压力和孔隙大小的关系。汞进入孔隙所需压力取决于汞的表面张力、接触角和孔径，测试不同压力下进入多孔试样中汞的体积量，就可以计算出相应压力下大于某个半径值的孔隙的体积，从而计算出多孔介质的孔隙尺寸分布和比表面积[10]。
2. 分形压汞实验
本次压汞实验使用Quantachrome公司POREMASTER系列压汞仪，实验压汞仪如图 3。该压汞仪含低压和高压站口2个，低压范围1.5~350 kPa，高压范围140 kPa~420 Mpa，可测量直径为0.0035~400 μm范围内变化的孔容。该压汞仪测试原理如图 4，将待测样品放入样品管中，随着压力提高汞液体会逐渐进入到样品孔隙中，图中横轴为汞液体压力，纵轴为进入孔隙的体积，由此可以绘制出原始的体积-压力曲线，由该曲线可以进一步分析得出其他孔隙参数。
本次实验采用的样品为不锈钢丝网填充的回热器，5个样品填充参数如表 1所示。
表 1 实验样品参数Table 1. Sample parametersNumber 1# 4# 7# 10# 13# mesh number 500 420 220 420 500 mass/g 3.53 2.64 2.15 2.27 3.53 Filling method Manual Manual Manual Manual Mechanical 以上5个样品主要区别在于填充的网片目数和质量不同，同时1#、4#、7#、10#样品为手工填充，而13#样品为机器填充。
对以上5个样品进行压汞测试，初始压力-体积曲线实验结果如图 5所示，其中横坐标为压力，纵坐标为体积。
3. 回热器分形特征分析
Washburn方程描述了将液态汞压入特定直径的孔所需要的压力，方程如下：
$$D = \frac{{ - 4\gamma \left( {cos\theta } \right)}}{P}$$ (6) 式中：P为压力；D为孔隙直径；θ为汞和孔隙壁面的接触角，在本次测试中为140°；γ为汞的表面张力，取0.48 N/m。
根据式(5)以及实验测得的体积-压力关系，即可得出样品的孔隙分布与孔隙率等参数。
图 6即为5个样品的孔隙分布情况，其中横坐标为孔隙直径，纵坐标为体积对直径对数的微分。由该图可以看出对于使用220目~500目网片填充的回热器其孔径分布在20~100 μm之间，在多孔介质分类中属于大孔[14]。样品1#、4#、7#目数分别为500、420和220目，目数逐渐减小同时由图中可以看出孔径分布逐渐增大，与实际经验相符。其中1#样品500目孔径为20~60 μm，4#样品420目孔径为20~50 μm，7#样品220目孔径为40~100 μm。同时可以看出由于目数的减少样品孔隙率逐渐升高。10#样品目数也为420目，但质量相较于4#样品较小，两个样品孔径分布大致相同。13#样品为500目相较于其他样品不同在于使用机器填充，孔径分布与手工填充大致相同，但出现了小部分50~100 μm孔径，考虑该部分孔隙是由于填充过程中不均匀造成的。同时由图 6孔隙分布情况可以看出，除13#样品外，其余样品均为出现孔隙分布分段情况，最大、最小孔之间为连续分布。
5个样品孔隙率如表 2。
表 2 样品孔隙率及分形维数Table 2. Porosity & DfNumber 1# 4# 7# 10# 13# Porosity 63.94% 69.41% 75.15% 72.44% 68.72% Df 2.592 2.639 2.763 2.681 2.351 样品1#、4#、7#目数逐渐减小，同时由表中可以看出孔隙率逐渐增大。1#与13#目数都为500目，但13#比1#孔隙率高5%，机填相对于手工填充孔隙率更大。4#与10#样品网片目数都为420目，但孔隙率有一定差别，可以看出由手工填充的网片因填充工艺限制会造成孔隙率变化。
根据式(4)可对5个回热器样品进行分形维数计算，计算结果见表 2。
由表 2可知，1#、4#、7#、10#、13#四个样品分形维数在2~3之间，且随着网片目数的减小分形维数增大，主要是由于试样孔隙率的增加。13#样品分形维数计算得2.351，结合13#样品的孔隙分布情况，推测是由于样品制备过程中压力不均匀，导致出现两段孔隙分布，使得后续计算结果出现误差。由式(4)可知随着孔隙率增加分形维数会增大，同时最大、最小孔径之差增大会导致分形维数增大。
4. 结论
1）通过对回热器进行压汞实验，得到回热器中不同孔径的孔隙体积分数分布情况，得出回热器中孔隙尺寸为大孔，以及不同目数网片填充得到的孔隙率为60%~75%之间。
2）根据多孔介质分形维数的计算模型，得出回热器分形维数数值在2和3之间，说明回热器具有分形特征，相比其他具有分形结构特征的多孔介质，不锈钢网片填充的回热器内部结构较为均匀，且最大、最小孔隙差值较小，并未出现孔隙分布分段的现象。
-
![]()
图 3 变温场下玻璃纤维环氧树脂板的实验测试系统
Figure 3. Experimental setup of inconsistent temperature field
![]()
图 4 传感器PZT S0超声导波响应信号在局部加热300 s前后的幅值相位变化
Figure 4. The direct wave of PZT S0 gain-phase contrast after heating 300 s
![]()
图 5 局部持续加热状态下传感器PZT S0超声导波信号幅值衰减与相位延迟曲线
Figure 5. Gain-phase errors contrast of PZT S0 under variable temperature field
表 1 不同变温度环境下损伤定位结果及误差统计
Table 1 Damage localization results and errors under variable temperature
No. Temperature/℃ ${\hat r_1}$/cm ${\hat \theta _1}$/° E1r E1θ 1 21.87 25.7 93 0.7 3 2 29.54 24.3 93 0.7 1 3 39.80 24.4 93 0.6 3 4 49.86 25.9 93 0.9 3 5 60.65 24.6 91 0.4 1 6 70.77 23.8 93 1.2 3 7 78.67 25.7 93 0.7 3 
下载: 导出CSV
-
[1] WANG X Y, HE J J, GAO W H, et al. Three-dimensional damage quantification of low velocity impact damage in thin composite plates using phased-array ultrasound[J]. Ultrasonics, 2020, 110: 106264.
[2] 王志凌, 袁慎芳, 邱雷, 等. 基于压电超声相控阵方法的结构多损伤监测[J]. 振动测试与诊断, 2014(5): 796-801. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDCS201405002.htm WANG Z L, YUAN S F, QIU L, et al. Multi-damage monitoring method based on piezoelectric ultrasonic phased array[J]. Journal of Vibration, Measurement & Diagnosis, 2014(5): 796-801. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDCS201405002.htm
[3] YANG H, LEE Y J, Lee S K. Impact source localization in plate utilizing multiple signal classification[C]//Proceedings of the Institution of Mechanical Engineers Part C-Journal of Mechanical Engineering Science, 2013, 227: 703-713.
[4] 苏永振, 袁慎芳, 王瑜. 基于多重信号分类算法的复合材料冲击定位[J]. 复合材料学报, 2010(3): 105-110. https://www.cnki.com.cn/Article/CJFDTOTAL-FUHE201003016.htm SU Y Z, YUAN S F, WANG Y. Impact localization in composite using multiple signal classification method[J]. Acta Materiae Compositae Sinica, 2010(3): 105-110. https://www.cnki.com.cn/Article/CJFDTOTAL-FUHE201003016.htm
[5] ZHONG Y T, XIANG J W, CHEN X. et al. Multiple signal classification-based impact localization in composite structures using optimized ensemble empirical mode decomposition[J]. Appl. Sci. , 2018, 8(9): 1447. DOI: 10.3390/app8091447
[6] ZUO H, YANG Z, Xu C, et al. Damage identification for plate-like structures using ultrasonic guided wave based on improved MUSIC method[J]. Composite Structure, 2018, 203: 164-171. DOI: 10.1016/j.compstruct.2018.06.100
[7] Radecki R, Staszewski W J, Uhl T. Impact of changing temperature on lamb wave propagation for damage detection[J]. Key Engineering Materials, 2014, 588: 140-148. https://www.scientific.net/AMM.105-107.621
[8] Sikdar S, Fiborek P, Kudela P, et al. Effects of debonding on Lamb wave propagation in a bonded composite structure under variable temperature conditions[J]. Smart Materials & Structures, 2019, 28: 015021.
[9] Francesco Lanza di Scalea, Salvatore Salamone. Temperature effects in ultrasonic lamb wave structural health monitoring systems[J]. Journal of the Acoustical Society of America, 2008, 161: https://doi.org/10.1121/1.2932071.
[10] BAO Q, YUAN S F, WANG Y W, et al. Anisotropy compensated MUSIC algorithm based composite structure damage imaging method[J]. Composite Structure, 2019, 214: 293-303.
[11] ZHANG Z H, ZHONG Y T, XIANG J W, et al. Phase correction improved multiple signal classification for impact source localization under varying temperature condition[J]. Measurement, 2020, 152: 107374.
-
期刊类型引用(0)
其他类型引用(1)
计量
- 文章访问数: 118
- HTML全文浏览量: 30
- PDF下载量: 21
- 被引次数: 1
