摘要
为研究温度对混凝土超声测试尾波信号的影响规律;将信号向量间的归一化夹角作为波动指标,反映温度效应引起的信号变化;通过向量残差矩阵SVD获得表征温度效应大小的特征向量,建立向量空间映射和温度差的数学关系。在实验室采集混凝土梁尾波信号进行验证,结果表明,随温度升高尾波信号的波形发生后移,文中方法可分段线性量化温度效应;基于量化结果,得到常温下超声尾波信号最敏感的温度区间;任意4.5 ℃范围内,可去除74%~90%的温度效应。
超声波在非均匀介质中传播时会发生散射和衍射,使得一部分波的传播路径更长,称这一部分为尾波。Larose
超声信号易受环境温度影响这一问题在尾波信号上体现得尤为明显,无论在实验室还是实际工程中,区分温度和应力引起的声学特征变化较困难。Snieder
为实现混凝土材料超声测试温度效应的量化,研究温度效应量化与去除的内在联系,将超声信号视为高维空间中的向量,考察信号向量在高维空间的方向变化规律。基于向量残差SVD提取表征温度效应大小的特征向量,在高维空间建立几何映射关系,实现温度效应的量化和去除。在实验室开展混凝土梁足尺试件超声测试实验,采集自重状态的超声尾波信号,验证文中方法的可行性。
将尾波信号视为高维空间中的向量,对于能量归一化的单位信号向量(当未特别说明时,下文的向量包括矩阵中的列向量均为单位向量),其方向变化由温度效应引起。利用SVD提取主成分的能力,锁定温度效应对测试信号产生影响的向量,实现温度效应的量化和去除。
为衡量样本信号与参考信号间的差异性,设样本信号与参考信号均为单位列向量,计算2信号向量间的归一化夹角,即波动指标
, | (1) |
式中:波动指标的取值为[0,1],当波动指标为0时,表示2向量方向一致;当波动指标为1时,表示2向量相互正交,即2信号不存在相关性。
1) 信号向量的残差
设信号样本为矩阵,在样本信号中选定1个作为基准信号,所有信号向量减去在基准信号向量上的投影,即
, | (2) |
式中:为信号向量的残差矩阵;T表示转置运算。的任一列向量在平均信号上的投影都为0,即与垂直,包含了温度效应的总体影响(中的列向量为残差向量不可单位化)。
2) 残差矩阵的SVD
残差矩阵为m×n阶矩阵,由奇异值分解可得
, | (3) |
式中:为m×m阶酉矩阵;为m×n阶对角矩阵;为n×n阶酉矩阵,是的共轭转置。其中,
, | (4) |
式中:为对角阵,对角线上的元素均非负且按递减顺序排列;r为非零元的个数,;的对角元素称为矩阵Da的奇异值;的列向量称为的左奇异向量;的列向量称为的右奇异向量。
3) 温度效应量化与去除
信号受影响的因素包括环境温度、环境湿度、随机噪声、系统误差等,视环境温度效应为主要影响,建立与温度差的关系反映测试信号的温度效应。
, | (5) |
式中:为效应指标;为环境温度差。样本信号偏离平均方向的程度可由残差矩阵的奇异向量反映。
, | (6) |
式中:是样本信号;可表征温度效应的大小,取值为[-1,1];为温度效应的特征向量。当=1(或-1)时,信号与温度效应方向一致(或相反),即与真实信号正交,效应达到正向(或反向)最大;当=0时,信号与真实信号方向一致,没有温度效应。
由温度效应量化的假设,利用
。 | (7) |
在实验室开展混凝土工字梁自重状态下的超声测试实验,同步记录环境温度,实验系统如

图1 超声测试系统示意图(单位:cm)
Fig. 1 Test system
制作了2片同尺寸的钢筋混凝土工字梁(记作1#、2#梁),尺寸如

图2 实验梁配筋图(单位:cm)
Fig. 2 Reinforcement diagram of test beam
超声实验梁 | 钢筋配置 | 直径为10 mm的HPB300级钢筋,箍筋间距15 cm。 |
---|---|---|
C30砼 配合比 | 水泥∶细骨料∶粗骨料∶水∶粉煤灰∶矿物 =1.00∶2.50∶3.85∶0.55∶0.22∶0.16。 | |
骨料性能 | 细骨料细度模数2.7,MB值1.1; 粗骨料为5~10 mm的碎石,含泥量为0.3%。 | |
实验梁支座 | 素混凝土试块,尺寸为:55 cm15 cm55 cm。 | |
超声测试仪 | 跨距 | 声程为实验梁长2 000 mm。 |
增益 | 增益不改变超声信息,取经验值80 dB。 | |
脉宽 | 为获得最大超声能量,取最大值5 μs。 | |
采样间隔 | 多次调试后取经验值8 μs。 | |
触发延时 | 为采集到尾波信号取仪器最大值9 999 μs。 | |
采样长度 | 采样长度越长信息越丰富,取最大值1 024。 | |
压电换能器 | 直径40 mm,中心频率50 kHz。 | |
环境温度计 | 采样间隔为5 min/次;量程:-30~65 ℃;精度:±0.3 ℃;分辨率:0.1 ℃。 |
超声测试使用RSM-SY5(T)非金属超声检测仪,电源为充电式电池,已验证电池更换和导线插拔对实验无明显影响。换能器采用JHP01型压电陶瓷换能器,使用环氧树脂耦合剂将换能器永久固定在梁端横截面上,一端发射一端接收。使用YEM-70L环境记录仪自动采集温湿度。
数据处理流程如

图3 数据处理流程
Fig. 3 Data processing flow
超声信号标准化处理后得到一级样本信号,按温度从小到大排序,依次进行下面2种平均,以消除权重对后续处理分析的影响和减小测试误差。
1)以最低温度对应的信号为基准,依次划分0.1 ℃为间隔的温度区格,认为各区格内的信号处于同一温度,将区格内的温度和对应的信号进行平均,记为均匀样本信号。
2)基于均匀样本信号,以每个温度点为中心把±范围内的温度和对应信号平均,即连续移动平均,得到二级样本信号。
为直观体现环境温度变化对超声测试尾波信号的影响规律,选择最低温度对应的信号为参考信号,3种样本信号分别与该参考信号做波动指标。图

图4 尾波信号波动指标与温度的关系
Fig. 4 Code wave signal fluctuation indicators as a function of temperature
根据前期的实验分析,测试仪自身温度对信号波动无明显影响,其他非随机影响因素可能是湿度效应引起
为分析温度对信号波形特征的影响规律,平均[,](t取整数,℃)范围内的二级样本信号作为t ℃的代表信号,样本信号能量归一化后的幅值记为归一化幅值。

图5 尾波波形与温度的关系
Fig. 5 The form of Code wave versus temperature
通过分析波形图可以发现,温度变化的主要效应是拉伸或压缩信号,次要效应是变形形状,尾波波形随环境温度上升发生规律性后移。尾波信号可近似看作某种受扰动的衰减正弦信号,而正弦信号的自相关函数为余弦函数。波动指标是该衰减正弦信号的自相关函数的映射,所以
研究表明,在2~3 ℃的环境温度变化时,能显著影响测试信
利用
1) 温度效应量化的最优区间
线性回归程度可用拟合优度表征。

图6 温度效应在不同温度区间的量化效果
Fig. 6 Quantitative results of temperature effect in different temperature ranges
2种情形随着拟合区间增大,拟合优度总体逐渐减小;分别取±2.0 ℃、±2.5 ℃及±3.0 ℃时,随着拟合区间增大,拟合优度变化曲线光滑下降,可取这3种量化温度效应。为了获得拟合优度较高的线性量化结果,综合比较,取=±2.0 ℃,=4.5 ℃为宜。
2) 温度效应量化
基于上述分析,将一级样本同一温度(±0.05 ℃)平均后,在±2.0 ℃范围内连续移动平均得到二级样本,并在4.5 ℃范围内量化温度效应。在选定的温度区间内,二级样本去除其在平均信号上的投影,得到信号残差矩阵。提取信号残差矩阵的奇异向量(效应的特征向量),二级样本在上的投影值即效应指标。

图7 温度效应量化结果
Fig. 7 Quantifying temperature effect
3) 温度效应去除
二级样本去除残差矩阵奇异向量上的投影,对比分析区间内去效应前后二级样本的波动指标。4个区间内,去效应前后波动指标均值明显减小,去效应比率分别为83%、77%、81%和85%,表明奇异向量上包含了大部分影响量,即温度效应。这样就得到了温度线性影响的方向向量。温度效应去除结果,如

图8 温度效应去除结果
Fig. 8 Temperature effect removal results
4) 温度效应量化和去除的关系
为分析奇异向量量化和去除温度效应的一般效果,将

图9 拟合优度与去效应比率的关系
Fig. 9 Relationship between goodness of fit and effect reduction ratio
由
上述方法得到残差矩阵时所使用的(二级)样本数量较大,实际工程应用时往往没有条件在某4.5 ℃范围内把全部温度对应的信号采集到,这就需要寻找一种可替代性的残差矩阵,期望以少量样本信号获得等价的奇异向量。
温度区间 | 7~11.5 ℃ | 11.5~16 ℃ | 16~20.5 ℃ | 20.5~25 ℃ |
---|---|---|---|---|
0.990 5 | 0.996 9 | 0.998 5 | 0.997 3 |
由
在实验室采集自重状态的梁片混凝土梁尾波信号,将信号向量间的归一化夹角作为波动指标,与环境温度变化建立联系。通过向量残差矩阵SVD提取温度产生影响的特征向量,线性量化和去除温度效应。得到结论如下。
1) 环境温度对两片梁的尾波信号有相同的影响规律,尾波波形随温度升高发生规律性后移;信号波动指标随环境温度变化近似呈类正弦变化,且同一温度下的波动指标存在幅度相近的带宽,表明除温度外还存在其他非随机影响量。
2) 可在分段区间线性量化温度效应,文中取7~11.5 ℃、11.5~16 ℃、16~20.5 ℃和20.5~25 ℃,4个区间量化时,效应指标与温度差线性拟合优度分别为0.90、0.89、0.92和0.90;11.5~16 ℃区间内效应指标灵敏度最大,尾波信号对温度效应最敏感。
3) 温度效应线性量化和去除效果与所选区间大小有关。选定相同基准温度不同温度跨度,温度效应线性量化拟合优度越高,去效应比率越大;当确定温度跨度,不同温度区间的去效应比率在一定范围内波动,任意4.5 ℃跨度内,去效应比率为74%~90%。
提高上述方法去温度效应的效果,优化该方法和量化温度以外环境因素的影响,为将来实际工程超声应力测试标定温度及其他影响提供支撑,是下一步研究的主要内容。
参考文献
Larose E, Hall S. Monitoring stress related velocity variation in concrete with a 2×10-5 relative resolution using diffuse ultrasound[J]. The Journal of the Acoustical Society of America, 2009, 125(4): 1853-1856. [百度学术]
Payan C, Garnier V, Moysan J. Determination of nonlinear elastic constants and stress monitoring in concrete by coda waves analysis[J]. European journal of environmental and civil engineering, 2008, 15(4): 519-531. [百度学术]
郭增伟, 田川, 计亲, 等. 混凝土工字梁三点弯曲荷载超声测试[J]. 重庆大学学报, 2019, 42(11): 88-97. [百度学术]
Guo, Z W, Tian C, Ji Q, et al. Three-point bending load ultrasonic test of concrete I-beam[J]. Journal of Chongqing University, 2019, 42(11): 88-97.(in Chinese) [百度学术]
Jiang H W, Zhan H Y, Zhang J Q, et al. Detecting Stress Changes and Damage in Full-Size Concrete T-Beam and Slab with Ultrasonic Coda Waves[J]. Journal of Structural Engineering, 2021, 147(9): 04021140. [百度学术]
Stähler S C, Sens-Schönfelder C, Niederleithinger E. Monitoring stress changes in a concrete bridge with coda wave interferometry[J]. The Journal of the Acoustical Society of America, 2011, 129(4): 1945-1952. [百度学术]
宋丽莉, 葛洪魁, 郭志伟, 等. 利用多次散射波监测介质性质变化的试验研究[J]. 岩石力学与工程学报, 2012, 31(4): 713-722. [百度学术]
Song L L, Ge H K, Guo Z W, et al. Experimental study of variation of media properties monitoring using multiple scattering waves[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(4): 713-722.(in Chinese) [百度学术]
Hu H W, Li D Y, Wang L, et al. An improved ultrasonic coda wave method for concrete behavior monitoring under various loading conditions[J]. Ultrasonics, 2021, 116: 106498. [百度学术]
Gondim R M L, Haach V G. Monitoring of ultrasonic velocity in concrete specimens during compressive loading-unloading cycles[J]. Construction and Building Materials, 2021, 302: 124218. [百度学术]
Finger C, Saydak L, Vu G, et al. Sensitivity of Ultrasonic Coda Wave Interferometry to Material Damage—Observations from a Virtual Concrete Lab[J]. Materials, 2021, 14(14): 4033. [百度学术]
Snieder R, Grêt A, Douma H, et al. Coda wave interferometry for estimating nonlinear behavior in seismic velocity[J]. Science, 2002, 295(5563): 2253-2255. [百度学术]
Wang X, Chakraborty J, Niederleithinger E. Noise reduction for improvement of ultrasonic monitoring using coda wave interferometry on a real bridge[J]. Journal of Nondestructive Evaluation, 2021, 40(1): 1-14. [百度学术]
Niederleithinger E, Wang X, Herbrand M, et al. Processing ultrasonic data by coda wave interferometry to monitor load tests of concrete beams[J]. Sensors, 2018, 18(6): 1971. [百度学术]
Wang X,Niederleithinger E,Hindersmann I.The installation of embedded ultrasonic transducers inside a bridge to monitor temperature and load influence using coda wave interferometry technique[J].Structural Health Monitoring,2022,21(3):913-927. [百度学术]
Larose E, De Rosny J, Margerin L, et al. Observation of multiple scattering of kHz vibrations in a concrete structure and application to monitoring weak changes[J]. Physical Review E, 2006, 73(1): 016609. [百度学术]
Lu Y H, Michaels J E. A methodology for structural health monitoring with diffuse ultrasonic waves in the presence of temperature variations[J]. Ultrasonics, 2005, 43(9): 717-731. [百度学术]
Zhang Y X, Abraham O, Tournat V, et al. Validation of a thermal bias control technique for Coda Wave Interferometry (CWI)[J]. Ultrasonics, 2012, 53(3): 658-664. [百度学术]
郑丹, 任涛. 频率及含水量对混凝土损伤超声检测的影响[J]. 水利学报, 2014, 45(S1): 90-94. [百度学术]
Zheng D, Ren T. Influence of frequency and water content on ultrasonic detection of concrete damage[J]. Journal of Hydraulic Engineering, 2014, 45(S1): 90-94.(in Chinese) [百度学术]
Croxford A J, Moll J, Wilcox P D, et al. Efficient temperature compensation strategies for guided wave structural health monitoring[J]. Ultrasonics, 2010, 50(4/5): 517-528. [百度学术]
Niederleithinger E, Wunderlich C. Influence of small temperature variations on the ultrasonic velocity in concrete[C]//AIP Conference Proceedings. Denver, Colorado, USA. AIP, 2013: 390-397. [百度学术]