一种基于步进差分声时差的空耦超声高空间分辨率应力测量方法。根据倾角组装测量装置；设置激励空耦超声换能器与接收空耦换能器的使Lamb波的传播距离为L0，测量此时的声时差TL0并代入利用标定的应力系数K1,L0获得L0段的平均应力σL0；保持激励空耦超声换能器位置不变，采用步进方式使接收空耦超声换能器移动ΔL，使Lamb波的传播距离增加至L；使L＝L0+ΔL；Lamb波的传播距离为L的应力系数为K1,L＝L/L0·K1,L0；获得L段的平均应力σL；通过σL0、σL以及ΔL段的应力值σΔL的关系得到。用以解决空耦超声Lamb波应力测量方法存在测量精度与空间分辨率的矛盾问题及耦合材料对待测件带来的二次污染问题。
1.一种基于步进差分声时差的空耦超声高空间分辨率应力测量方法，其特征在于，所述测量方法具体为，在碳纤维复合材料板(3)内激发出相对纯净的单一Lamb波模态从而确定激励空耦超声换能器(1)发出的激励与接收空耦换能器(2)的倾角，再根据倾角组装测量装置；
设置激励空耦超声换能器(1)与接收空耦换能器(2)的Lamb波传播距离为L0，测量此时的声时差TL0，利用标定的应力系数K1,L0获得L0段的平均应力σL0，将声时差TL0带入平均应力σL0；
保持激励空耦超声换能器(1)位置不变，采用步进方式使接收空耦超声换能器移动ΔL，使Lamb波的传播距离增加至L；使L＝L0+ΔL；
Lamb波的传播距离为L的应力系数为K1,L＝L/L0·K1,L0；
通过σL0、σL以及ΔL段的应力值σΔL的关系得到，σΔL＝(TL‑TL0)/(K1,L0ΔL)＝ΔT/(K1,L0ΔL) (19)。
2.根据权利要求1所述空耦超声高空间分辨率应力测量方法，其特征在于，所述组装测量装置具体为，所述空耦换能器(1)与碳纤维复合材料板(3)的夹角为θT，所述接收空耦换能器(2)与碳纤维复合材料板(3)的夹角为θR；
在包含应力集中区域的碳纤维复合材料板(3)的一侧，按照确定好的倾角放置激励空耦超声换能器(1)与接收空耦换能器(2)；
所述碳纤维复合材料板(3)的另一侧设置接收空耦换能器(2)，所述接收空耦换能器(2)与空耦换能器(1)在碳纤维复合材料板(3)的同一侧。
3.根据权利要求2所述空耦超声高空间分辨率应力测量方法，其特征在于，确定激励信号周期为N，并选择合适空耦超声换能器的中心频率f，采用任意信号函数发生器产生N周期、Hanning窗调制、频率为f的正弦脉冲信号作为激励信号，经过50欧姆负载进行阻抗匹配，并经过低通滤波器后，施加在激励空耦超声换能器(1)上。
4.根据权利要求1所述空耦超声高空间分辨率应力测量方法，其特征在于，当Lamb波的传播距离为L，其应力系数为K1,L＝L/L0·K1,L0时，测量此时的声时差TL，利用应力系数K1,L获得L段的平均应力σL，将声时差TL带入平均应力σL。
5.根据权利要求2所述空耦超声高空间分辨率应力测量方法，其特征在于，所述Lamb波超声应力检测具体为，根据线弹性理论，小振幅超声波在处于自然状态的介质中传播时的波动方程为，设谐振平面波为，
ui＝Aiexp[i(kmxm‑ωt)] (2)式中，Ai＝Aαi，A为谐振波合振幅，αi为质点位移的方向余弦，km为单位波矢量；
式中，v为声速，δil为Kronecher函数，Γil＝Cijklljlk为Christoffel声张量，lj、lk为波
前法向量的方向余弦，满足ui＝ulδil、kj＝klj、kk＝klk、v＝ω/k。
6.根据权利要求5所述空耦超声高空间分辨率应力测量方法，其特征在于，当弹性波在处于变形形态的受力介质中传播时，变形状态坐标系下的波动方程更改为令Aijkl＝Cijkl+tjlδik为受力介质的等效刚度系数，得到各向异性介质在受力变形状态时的Christoffel方程，式(3)和式(5)存在非零解的条件为
det|ρv0δil‑Cijklljlk|＝0 (6)2
det|ρvδil‑Aijklljlk|＝0 (7)式中，v0为各向异性介质不受力时的声速，v为各向异性介质受力时的声速，故2
由于速度在受力前后的变化往往是微弱的，故v≈v0，联立式(8)和式(9)得到由上式可知，在受力的各向异性介质中，超声波的速度变化与应力存在近似线性关系。
7.根据权利要求6所述空耦超声高空间分辨率应力测量方法，其特征在于，已知两个主应力对Lamb波相速度的影响及两个正应力的Lamb声弹性效应满足线性叠加关系，Δv1＝K11σ11+K12σ22 (11)Δv2＝K21σ11+K22σ22 (12)式中，Δv1和Δv2分别是在x1和x2方向上有无应力前后的Lamb波相速度变化，σ11和σ22分别为x和y方向上的两个主应力，K11、K12、K21、K22是声弹性系数对于x1轴上的单轴应力场，式(11)和式(12)可以简化为Δv1＝K11σ11 (13)
通常在实际测量过程中，声速的变化很微弱，直接测量声速会引入更多误差；为了测量更加直接，引入物理量声时差T，将声速的测量转化为对声时差的测量；
将式(13)代入式(14)即可得到声时差T与应力的关系，式中，K1为应力系数；
由式(15)可知，通过实验标定获得声弹性系数K1，即改变σ11的大小，获得x1方向上Lamb波在相同传播距离下的声时差T1，即可实现碳纤维复合材料板x1轴上的单轴应力场的应力测量。
8.根据权利要求1所述空耦超声高空间分辨率应力测量方法，其特征在于，公式(19)具体为，σL0、σL以及ΔL段的应力值σΔL满足以下关系式，σL0＝TL0/K1,L0 (16)
联立式(16)～(18)得ΔL段的应力值σΔL为，σΔL＝(TL‑TL0)/(K1,L0ΔL)＝ΔΤ/(K1,L0ΔL)(19)式中，ΔΤ＝ΤL‑TL0，称之为步进差分声时差，为步进ΔL前后两声时差TL0和TL之差。
基于步进差分声时差的空耦超声高空间分辨率应力测量方法
技术领域
[0001] 本发明属于超声波检测领域，具体涉及基于步进差分声时差的空耦超声高空间分辨率应力测量方法。
背景技术
[0002] 2010年左右，波音787问世，随后是空客A350，两者都有超过50％的飞机结构重量由复合材料制成，相当于体积的80％。在飞机中广泛使用复合材料的主要驱动因素是在保持所需的承载能力的同时显著降低了结构重量，同时兼顾了经济和环境目标。碳纤维增强树脂基复合材料具有性能可设计、多功能兼容、可整体成型等特点，以及高比强度和比刚度、疲劳性能好、耐腐蚀等特点，广泛应用于新型飞机的机身、机翼和机舱以及汽车、船舶等领域。复合材料中的应力主要表现为残余应力和结构承力，其中有些可以弥合部分分层、裂纹等损伤而增加结构的强度和稳定性。然而，有些则会使分层、裂纹等损伤加剧，进而造成结构强度和稳定性大幅下降，从而影响其应用领域的安全性。因此，发展适用于各向异性碳纤维复合材料板的应力测量方法，对其应力状态及时进行快速、精确的测量，对于实现其应用领域的健康监测和评价具有重要的指导意义。
[0003] 根据是否使用耦合剂(如水或凝胶)，超声检测技术分为接触式超声检测和非接触式超声检测。基于声弹性效应的接触式超声检测技术需要使用耦合剂实现换能器和待检测样本之间的有效耦合。然而，耦合剂的厚度不均匀会直接影响声时差的测量从而引入应力测量误差。此外，耦合剂的存在会影响测量的灵活性，因此，非接触式超声检测更适合用于实现碳纤维复合材料板的在役快速应力测量。非接触式超声检测技术主要包括电磁超声技术、激光超声技术以及空气耦合超声技术。电磁超声技术使用磁场在试样中激发超声波，要求试样由导电材料制成。故电磁超声仅限于金属，不适用于碳纤维复合材料。激光超声检测技术使用短激光脉冲瞬间加热试样表面的一个点，导致快速热膨胀，从而激发超声波。激光超声检测技术受限于能量，存在烧蚀试样的风险。空气耦合超声检测技术以空气作为耦合剂，具有完全非接触和无污染的特点。尤其对于复合材料，其声阻抗更接近于空气的声阻抗，更容易在复合材料板中激励产生模态单一的Lamb波。因此，基于Lamb波的空气耦合超声检测技术对于实现碳纤维复合材料板完全非接触的在役应力测量是可行的。
[0004] 传统基于声时差的应力测量方法，是获得的固定传播距离下的平均应力，应力测量的空间分辨率直接取决于传播距离的大小。然而，为了实现碳纤维复合材料板的原位应力测量，空耦超声换能器对需采用同侧俯仰式进行Lamb波的激励和接收。为了消除直达波对Lamb波信号的干扰，提高应力测量精度，两空耦超声换能器之间需要有足够长的距离，从而降低了应力测量的空间分辨率。因此，空耦超声Lamb波应力测量方法存在测量精度与空间分辨率的矛盾问题。
发明内容
[0005] 本发明提供基于步进差分声时差的空耦超声高空间分辨率应力测量方法，用以解决空耦超声Lamb波应力测量方法存在测量精度与空间分辨率的矛盾问题及耦合材料对待测件带来的二次污染问题。
[0006] 本发明通过以下技术方案实现：
[0007] 一种基于步进差分声时差的空耦超声高空间分辨率应力测量方法，所述测量方法具体为，在碳纤维复合材料板3内激发出相对纯净的单一Lamb波模态从而确定激励空耦超声换能器1发出的激励与接收空耦换能器2的倾角，再根据倾角组装测量装置；
[0008] 设置激励空耦超声换能器1与接收空耦换能器2的Lamb波的传播距离为L0，测量此时的声时差TL0，利用标定的应力系数K1,L0获得L0段的平均应力σL0，将声时差TL0带入平均应力σL0；
[0009] 保持激励空耦超声换能器1位置不变，采用步进方式使接收空耦超声换能器移动ΔL，使Lamb波的传播距离增加至L；使L＝L0+ΔL；
[0010] Lamb波的传播距离为L的应力系数为K1,L＝L/L0·K1,L0；获得L段的平均应力σL；
[0011] 通过σL0、σL以及ΔL段的应力值σΔL的关系得到，
[0012] σΔL＝(TL‑TL0)/(K1,L0ΔL)＝ΔT/(K1,L0ΔL) (19)。
[0013] 进一步的，所述组装测量装置具体为，所述空耦换能器(1)与碳纤维复合材料板3的夹角为θT，所述接收空耦换能器2与碳纤维复合材料板3的夹角为θR；
[0014] 在包含应力集中区域的碳纤维复合材料板3的一侧，按照确定好的倾角放置激励空耦超声换能器1与接收空耦换能器2；
[0015] 所述碳纤维复合材料板3的另一侧设置接收空耦换能器2，所述接收空耦换能器2与空耦换能器1在碳纤维复合材料板3的同一侧。
[0016] 进一步的，确定激励信号周期为N，并选择合适空耦超声换能器的中心频率f，采用任意信号函数发生器产生N周期、Hanning窗调制、频率为f的正弦脉冲信号作为激励信号，经过50欧姆负载进行阻抗匹配，并经过低通滤波器后，施加在激励空耦超声换能器1上。
[0017] 进一步的，当Lamb波的传播距离为L，其应力系数为K1,L＝L/L0·K1,L0时，测量此时的声时差TL，利用应力系数K1,L获得L段的平均应力σL，将声时差TL带入平均应力σL。
[0018] 进一步的，所述Lamb波超声应力检测具体为，根据线弹性理论，小振幅超声波在处于自然状态的介质中传播时的波动方程为，
[0019]
[0020] 设谐振平面波为，
[0021] ui＝Aiexp[i(kmxm‑ωt)] (2)
[0022] 式中，Ai＝Aαi，A为谐振波合振幅，αi为质点位移的方向余弦，km为单位波矢量。
[0023] 通过公式(1)和(2)联立可得
[0024] (ρv2δil‑Γil)ul＝0 (3)
[0025] 式中，v为声速，δil为Kronecher函数，Γil＝Cijklljlk为Christoffel声张量，lj、lk2 2 2
为波前法向量的方向余弦，满足ui＝ulδil、kj＝klj、kk＝klk、v＝ω/k。
[0026] 进一步的，当弹性波在处于变形形态的受力介质中传播时，变形状态坐标系下的波动方程更改为
[0027]
[0028] 令Aijkl＝Cijkl+tjlδik为受力介质的等效刚度系数，得到各向异性介质在受力变形状态时的Christoffel方程，
[0029] (ρv2δil‑Aijklljlk)ul＝0 (5)
[0030] 式(3)和式(5)存在非零解的条件为
[0031] det|ρv02δil‑Cijklljlk|＝0 (6)
[0032] det|ρv2δil‑Aijklljlk|＝0 (7)
[0033] 式中，v0为各向异性介质不受力时的声速，v为各向异性介质受力时的声速，故[0034] ρv02δil＝Cijklljlk (8)
[0035] ρv2δil＝Aijklljlk (9)
[0036] 由于速度在受力前后的变化往往是微弱的，故v≈v0，联立式(8)和式(9)得到[0037]
[0038] 由上式可知，在受力的各向异性介质中，超声波的速度变化与应力存在近似线性关系。
[0039] 进一步的，已知两个主应力对Lamb波相速度的影响及两个正应力的Lamb声弹性效应满足线性叠加关系，
[0040] Δv1＝K11σ11+K12σ22 (11)
[0041] Δv2＝K21σ11+K22σ22 (12)
[0042] 式中，Δv1和Δv2分别是在x1和x2方向上有无应力前后的Lamb波相速度变化，σ11和σ22分别为x和y方向上的两个主应力，K11、K12、K21、K22是声弹性系数
[0043] 对于x1轴上的单轴应力场，式(11)和式(12)可以简化为
[0044] Δv1＝K11σ11 (13)
[0045] 通常在实际测量过程中，声速的变化很微弱，直接测量声速会引入更多误差。为了测量更加直接，引入物理量声时差T，将声速的测量转化为对声时差的测量；
[0046]
[0047] 将式(13)代入式(14)即可得到声时差T与应力的关系，
[0048]
[0049] 式中，K1为应力系数
[0050] 由式(15)可知，通过实验标定(即改变σ11的大小，获得x1方向上Lamb波在相同传播距离下的声时差T1)获得声弹性系数K1，即可实现碳纤维复合材料板x1轴上的单轴应力场的应力测量。
[0051] 进一步的，公式(19)具体为，σL0、σL以及ΔL段的应力值σΔL满足以下关系式，[0052] σL0＝TL0/K1,L0 (16)
[0053] σL＝TL/K1,L (17)
[0054] σΔLΔL+σL0L0＝σLL (18)
[0055] 联立式(31)～(33)得ΔL段的应力值σΔL为，
[0056] σΔL＝(TL‑TL0)/(K1,L0ΔL)＝ΔΤ/(K1,L0ΔL) (19)
[0057] 式中，ΔΤ＝ΤL‑TL0，称之为步进差分声时差，为步进ΔL前后两声时差TL0和TL之差。
[0058] 本发明的有益效果是：
[0059] 本发明通过步进前后声时差之差作为应力表征指标，在消除直达波干扰保证应力测量精度的同时，实现高空间分辨率应力测量。
[0060] 本发明采用空气耦合Lamb波超声检测，在检测过程中以空气作为传输媒质代替了传统超声无损检测中的耦合剂。同时，Lamb波在传播过程中能量衰减小且传播距离远，能够大范围、高效率地实现应力的测量。
[0061] 本发明在检测过程中具有完全无接触、无侵入和无损害的优势，也可大大延长空耦超声换能器的使用寿命，增强了检测的灵活性、提高了检测效率，更好地实现碳纤维复合材料板的原位应力测量。
附图说明
[0062] 图1是本发明的结构示意图。
具体实施方式
[0063] 下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0064] 一种基于步进差分声时差的空耦超声高空间分辨率应力测量方法，所述测量方法具体为，在碳纤维复合材料板3内激发出相对纯净的单一Lamb波模态从而确定激励空耦超声换能器1发出的激励与接收空耦换能器2的倾角，再根据倾角组装测量装置；
[0065] 以消除直达波影响为依据，设置激励空耦超声换能器1与接收空耦换能器2的Lamb波的传播距离为L0，测量此时的声时差TL0，利用标定的应力系数K1,L0获得L0段的平均应力σL0，将声时差TL0带入平均应力σL0；
[0066] 保持激励空耦超声换能器1位置不变，采用步进方式使接收空耦超声换能器移动ΔL，使Lamb波的传播距离增加至L；使L＝L0+ΔL；
[0067] 由于应力系数与传播距离成线性关系，故Lamb波的传播距离为L的应力系数为K1,L＝L/L0·K1,L0；
[0068] 测量此时的声时差TL并代入利用应力系数K1,L获得L段的平均应力σL；
[0069] 通过σL0、σL以及ΔL段的应力值σΔL的关系得到，
[0070] σΔL＝(TL‑TL0)/(K1,L0ΔL)＝ΔT/(K1,L0ΔL) (19)。
[0071] 所述空耦超声高空间分辨率应力测量方法，所述组装测量装置具体为，为在碳纤维复合材料板内激发出相对纯净的单一Lamb波模态，需要利用对应Lamb波模态的相速度以及空气声速并结合Snell定律确定激励与接收空耦超声换能器的倾角；
[0072] 所述空耦换能器1与碳纤维复合材料板3的夹角为θT，所述接收空耦换能器2与碳纤维复合材料板3的夹角为θR；
[0073] 在包含应力集中区域的碳纤维复合材料板3的一侧，按照确定好的倾角放置激励空耦超声换能器1与接收空耦换能器2；
[0074] 所述碳纤维复合材料板3的另一侧设置接收空耦换能器2，所述接收空耦换能器2与空耦换能器1在碳纤维复合材料板3的同一侧。
[0075] 所述空耦超声高空间分辨率应力测量方法，为了保证声束具有足够能量，确定激励信号周期为N，并选择合适空耦超声换能器的中心频率f，采用任意信号函数发生器产生N周期、Hanning窗调制、频率为f的正弦脉冲信号作为激励信号，经过50欧姆负载进行阻抗匹配，并经过低通滤波器后，施加在激励空耦超声换能器1上。
[0076] 所述空耦超声高空间分辨率应力测量方法，当Lamb波的传播距离为L，其应力系数为K1,L＝L/L0·K1,L0时，测量此时的声时差TL，利用应力系数K1,L获得L段的平均应力σL，将声时差TL带入平均应力σL。
[0077] 所述空耦超声高空间分辨率应力测量方法，所述Lamb波超声应力检测具体为，根据线弹性理论，小振幅超声波在处于自然状态的介质中传播时的波动方程为，[0078]
[0079] 设谐振平面波为，
[0080] ui＝Aiexp[i(kmxm‑ωt)] (2)
[0081] 式中，Ai＝Aαi，A为谐振波合振幅，αi为质点位移的方向余弦，km为单位波矢量。
[0082] 通过公式(1)和(2)联立可得
[0083] (ρv2δil‑Γil)ul＝0 (3)
[0084] 式中，v为声速，δil为Kronecher函数，Γil＝Cijklljlk为Christoffel声张量，lj、lk2 2 2
为波前法向量的方向余弦，满足ui＝ulδil、kj＝klj、kk＝klk、v＝ω/k。
[0085] 所述空耦超声高空间分辨率应力测量方法，当弹性波在处于变形形态的受力介质中传播时，变形状态坐标系下的波动方程更改为
[0086]
[0087] 令Aijkl＝Cijkl+tjlδik为受力介质的等效刚度系数，得到各向异性介质在受力变形状态时的Christoffel方程，
[0088] (ρv2δil‑Aijklljlk)ul＝0 (5)
[0089] 式(3)和式(5)存在非零解的条件为
[0090] det|ρv02δil‑Cijklljlk|＝0 (6)
[0091] det|ρv2δil‑Aijklljlk|＝0 (7)
[0092] 式中，v0为各向异性介质不受力时的声速，v为各向异性介质受力时的声速，故[0093] ρv02δil＝Cijklljlk (8)
[0094] ρv2δil＝Aijklljlk (9)
[0095] 由于速度在受力前后的变化往往是微弱的，故v≈v0，联立式(8)和式(9)得到[0096]
[0097] 由上式可知，在受力的各向异性介质中，超声波的速度变化与应力存在近似线性关系。
[0098] 所述空耦超声高空间分辨率应力测量方法，已知两个主应力对Lamb波相速度的影响及两个正应力的Lamb声弹性效应满足线性叠加关系，
[0099] Δv1＝K11σ11+K12σ22 (11)
[0100] Δv2＝K21σ11+K22σ22 (12)
[0101] 式中，Δv1和Δv2分别是在x1和x2方向上有无应力前后的Lamb波相速度变化，σ11和σ22分别为x和y方向上的两个主应力，K11、K12、K21、K22是声弹性系数
[0102] 对于x1轴上的单轴应力场，式(11)和式(12)可以简化为
[0103] Δv1＝K11σ11 (13)
[0104] 通常在实际测量过程中，声速的变化很微弱，直接测量声速会引入更多误差。为了测量更加直接，引入物理量声时差T，将声速的测量转化为对声时差的测量；
[0105]
[0106] 将式(13)代入式(14)即可得到声时差T与应力的关系，
[0107]
[0108] 式中，K1为应力系数
[0109] 由式(15)可知，通过实验标定(即改变σ11的大小，获得x1方向上Lamb波在相同传播距离下的声时差T1)获得声弹性系数K1，即可实现碳纤维复合材料板x1轴上的单轴应力场的应力测量。
[0110] 所述空耦超声高空间分辨率应力测量方法，公式19)具体为，σL0、σL以及ΔL段的应力值σΔL满足以下关系式，
[0111] σL0＝TL0/K1,L0 (16)
[0112] σL＝TL/K1,L (17)
[0113] σΔLΔL+σL0L0＝σLL (18)
[0114] 联立式(31)～(33)得ΔL段的应力值σΔL为，
[0115] σΔL＝(TL‑TL0)/(K1,L0ΔL)＝ΔΤ/(K1,L0ΔL) (19)
[0116] 式中，ΔΤ＝ΤL‑TL0，称之为步进差分声时差，为步进ΔL前后两声时差TL0和TL之差。
