CN108760890A - 一种使用信号能谱系数监测脆性材料损伤演化过程的方法 - Google Patents

一种使用信号能谱系数监测脆性材料损伤演化过程的方法 Download PDF

Info

Publication number
CN108760890A
CN108760890A CN201810291038.7A CN201810291038A CN108760890A CN 108760890 A CN108760890 A CN 108760890A CN 201810291038 A CN201810291038 A CN 201810291038A CN 108760890 A CN108760890 A CN 108760890A
Authority
CN
China
Prior art keywords
frequency band
band signals
window
fragile material
acoustic emission
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810291038.7A
Other languages
English (en)
Other versions
CN108760890B (zh
Inventor
宁建国
任会兰
王宗炼
马天宝
李健
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201810291038.7A priority Critical patent/CN108760890B/zh
Publication of CN108760890A publication Critical patent/CN108760890A/zh
Application granted granted Critical
Publication of CN108760890B publication Critical patent/CN108760890B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/14—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object using acoustic emission techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/4409—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/46—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by spectral analysis, e.g. Fourier analysis or wavelet analysis
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/023—Solids
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/028—Material parameters
    • G01N2291/0289—Internal structure, e.g. defects, grain size, texture

Abstract

本发明涉及一种使用信号能谱系数监测脆性材料损伤演化过程的方法，属于材料损伤检测技术领域。其步骤为：①在待监测脆性材料上布置声发射传感器，采集声发射信号；②根据声发射信号的波形特征选取小波基；③对声发射信号进行2级小波包分解，分解成AA2，DA2，AD2和DD24个频段信号；④计算AA2频段信号的能谱系数P1；⑤以20个声发射信号为滑动窗口，按采集次序依次选取20个声发射信号为一个窗口，计算每个窗口中P1的平均值，并绘制P1随窗口数的分布图；⑥根据P1的变化状态监测所述待监测脆性材料损伤演化过程。本发明提出的方法可以定量获取材料内部细观结构变化的物理信息，实时监测脆性材料损伤演化过程。

Description

一种使用信号能谱系数监测脆性材料损伤演化过程的方法
技术领域
本发明涉及一种使用信号能谱系数监测脆性材料损伤演化过程的方法，属于材料损伤检测技术领域。
背景技术
脆性材料中往往存在着大量弥散的初始微裂纹，微裂纹损伤是脆性材料破坏的显著特征。微裂纹的萌生、扩展和汇合将导致材料逐渐劣化直至最终的破坏失效。材料的变形与破坏是一种能量耗散的不可逆过程，能量变化贯穿于材料变形破坏的整个过程。脆性材料中裂纹在萌生和扩展演化时，局部源快速释放能量产生瞬态弹性波，即声发射现象，声发射也称为应力波发射。声发射技术作为一种物理检测方法可以实时获取材料内部细观结构变化的物理信息。因此，声发射技术在脆性材料损伤检测领域具有广阔的应用前景。
不同的损伤模式对应的声发射信号具有不同的频率特征。声发射信号的小波包特征能谱系数能够有效地表征声发射信号的频率信息，在定量声发射分析方面扮演着重要作用。已有学者证实脆性材料损伤破坏过程中裂纹源的尺度与声发射信号的频率成反相关的关系。因此，基于声发射信号的能谱系数可以有效地判断裂纹的长度及危险程度，并监测脆性材料损伤演化过程。
发明内容
本发明的目的是提供一种使用信号能谱系数监测脆性材料损伤演化过程的方法，该方法可以定量获取材料内部细观结构变化的物理信息，实时监测脆性材料损伤演化过程。
本发明的目的是通过以下技术方案实现的。
一种使用信号能谱系数监测脆性材料损伤演化过程的方法，具体操作步骤为：
步骤一、在待监测脆性材料上布置声发射传感器，用于实时采集待监测脆性材料损伤破坏产生的声发射信号。
所述脆性材料为陶瓷类材料、岩石和混凝土。
步骤二、根据步骤一采集到的声发射信号的波形特征，选取小波基。
步骤三、使用步骤二选取的小波基，对步骤一采集到的声发射信号进行2级小波包分解，声发射信号由低到高被分解成4个频段信号，分别用符号AA2，DA2，AD2和DD2表示；其中，AA2频段信号的频率范围为DA2频段信号的频率范围为AD2频段信号的频率范围为DD2频段信号的频率范围为fmax为所述声发射信号的最大频率。
步骤四、计算AA2频段信号的能谱系数，用符号P1表示；具体步骤为：
步骤4.1：用公式(1)表示AA2频段信号的小波包分解系数；用公式将(2)表示DA2频段信号的小波包分解系数；用公式(3)表示AD2频段信号的小波包分解系数；用公式将(4)表示DD2频段信号的小波包分解系数；
S1n＝{x1,n,n＝1,2,...,i/4} (1)
其中，S1n表示AA2频段信号的小波包分解系数集合；x1,n表示AA2频段信号的小波包分解的各个系数；i表示声发射信号的采样点数；
S2n＝{x2,n,n＝1,2,...,i/4} (2)
其中，S2n表示DA2频段信号的小波包分解系数集合；x2,n表示DA2频段信号的小波包分解的各个系数；
S3n＝{x3,n,n＝1,2,...,i/4} (3)
其中，S3n表示AD2频段信号的小波包分解系数集合；x3,n表示AD2频段信号的小波包分解的各个系数；
S4n＝{x4,n,n＝1,2,...,i/4} (4)
其中，S4n表示DD2频段信号的小波包分解系数集合；x4,n表示DD2频段信号的小波包分解的各个系数；
步骤4.2：用公式(5)表示AA2频段信号的能量，用符号E1表示；用公式将(6)表示DA2频段信号的能量，用符号E2表示；用公式(7)表示AD2频段信号的能量，用符号E3表示；用公式将(8)表示DD2频段信号的能量，用符号E4表示；
步骤4.3：通过公式(9)得到AA2频段信号的能谱系数P1；
P1＝100*(E1/(E1+E2+E3+E4)) (9)
步骤五、在步骤四的操作基础上，以20个发射信号为滑动窗口，按采集次序依次选取20个声发射信号为一个窗口，计算每个窗口中P1的平均值，并绘制P1随窗口数的分布图，横坐标为窗口的序号；纵坐标为对应窗口中P1的平均值。
步骤六、根据P1随窗口数的分布图中，P1的变化状态监测所述待监测脆性材料损伤演化过程。监测所述待监测脆性材料损伤演化过程的具体方法为：当窗口M之后连续10个窗口中P1的平均值大于该窗口及之前连续10个窗口中P1的平均值，则判断待监测脆性材料在M窗口由裂纹萌生阶段进入裂纹生长阶段；当窗口N，N>M，之后连续10个窗口中P1的平均值大于该窗口及之前连续10个窗口中P1的平均值，则判断待监测脆性材料在N窗口由裂纹生长阶段进入裂纹汇合阶段，裂纹汇合阶段持续到待监测脆性材料最终的破坏失效。
有益效果
本发明提出的一种使用信号能谱系数监测脆性材料损伤演化过程的方法与已有技术相比较，其优点为：可以定量获取材料内部细观结构变化的物理信息，实时监测脆性材料损伤演化过程。
附图说明
图1为本发明具体实施方式中db3小波基函数图；
图2为本发明具体实施方式中小波包分解流程图；
图3为本发明具体实施方式中P1随窗口数的分布图。
具体实施方式
下面结合附图和具体实施例对本发明技术方案作进一步说明。
实施例1
本实施例中，待监测脆性材料为混凝土材料，对混凝土进行三点弯破坏实验。
使用本发明提出的使用信号能谱系数监测脆性材料损伤演化过程的方法监测混凝土损伤演化过程，具体操作步骤为：
步骤一、在待监测脆性材料上布置声发射传感器，用于实时采集待监测脆性材料上的声发射信号。
步骤二、根据步骤一采集到的声发射信号的波形特征，选取db3小波基，db3小波基函数图如图1所示。
步骤三、使用步骤二选取的小波基，对步骤一采集到的声发射信号进行2级小波包分解，声发射信号由低到高被分解成AA2，DA2，AD2和DD24个频段信号，小波包分解流程如图2所示。其中，AA2频段信号的频率范围为(0,125)Hz，DA2频段信号的频率范围为(125,250)Hz，AD2频段信号的频率范围为(250,375)Hz，DD2频段信号的频率范围为(375,500)Hz。
步骤四、计算AA2频段信号的能量特征值P1；具体步骤为：
步骤4.1：用公式(1)表示AA2频段信号的小波包分解系数；用公式将(2)表示DA2频段信号的小波包分解系数；用公式(3)表示AD2频段信号的小波包分解系数；用公式将(4)表示DD2频段信号的小波包分解系数公式(1)、公式(2)、公式(3)和公式(4)中的i＝1024。
步骤4.2：用公式(5)表示AA2频段信号的能量E1；用公式将(6)表示DA2频段信号的能量E2；用公式(7)表示AD2频段信号的能量E3；用公式将(8)表示DD2频段信号的能量E4。
步骤4.3：通过公式(9)得到AA2频段信号的能谱系数P1。
步骤五、在步骤四的操作基础上，以20个声发射信号为滑动窗口，按采集次序依次选取20个声发射信号为一个窗口，计算每个窗口中P1的平均值，并绘制P1随窗口数的分布图，横坐标为窗口的序号；纵坐标为对应窗口中P1的平均值；如图3所示。
步骤六、根据P1随窗口数的分布图中，P1的变化状态监测所述待监测脆性材料损伤演化过程。监测所述待监测脆性材料损伤演化过程的具体方法为：当第17个窗口之后连续10个窗口中P1的平均值大于第17个窗口及之前连续10个窗口中P1的平均值，则判断待监测脆性材料在第17个窗口由裂纹萌生阶段进入裂纹生长阶段；当第53个窗口之后连续10个窗口中P1的平均值大于第53个窗口及之前连续10个窗口中P1的平均值，则判断待监测脆性材料在第53个窗口由裂纹生长阶段进入裂纹汇合阶段，裂纹汇合阶段持续到待监测脆性材料最终的破坏失效。
本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，熟悉本领域的技术人员可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。

Claims (3)

1.一种使用信号能谱系数监测脆性材料损伤演化过程的方法，其特征在于：具体操作步骤为：
步骤一、在待监测脆性材料上布置声发射传感器，用于实时采集待监测脆性材料损伤破坏产生的声发射信号；
步骤二、根据步骤一采集到的声发射信号的波形特征，选取小波基；
步骤三、使用步骤二选取的小波基，对步骤一采集到的声发射信号进行2级小波包分解，声发射信号由低到高被分解成4个频段信号，分别用符号AA2，DA2，AD2和DD2表示；其中，AA2频段信号的频率范围为DA2频段信号的频率范围为AD2频段信号的频率范围为DD2频段信号的频率范围为fmax为所述声发射信号的最大频率；
步骤四、计算AA2频段信号的能谱系数，用符号P1表示；具体步骤为：
步骤4.1：用公式(1)表示AA2频段信号的小波包分解系数；用公式将(2)表示DA2频段信号的小波包分解系数；用公式(3)表示AD2频段信号的小波包分解系数；用公式(4)表示DD2频段信号的小波包分解系数；
S1n＝{x1,n,n＝1,2,...,i/4} (1)
其中，S1n表示AA2频段信号的小波包分解系数集合；x1,n表示AA2频段信号的小波包分解的各个系数；i表示声发射信号的采样点数；
S2n＝{x2,n,n＝1,2,...,i/4} (2)
其中，S2n表示DA2频段信号的小波包分解系数集合；x2,n表示DA2频段信号的小波包分解的各个系数；
S3n＝{x3,n,n＝1,2,...,i/4} (3)
其中，S3n表示AD2频段信号的小波包分解系数集合；x3,n表示AD2频段信号的小波包分解的各个系数；
S4n＝{x4,n,n＝1,2,...,i/4} (4)
其中，S4n表示DD2频段信号的小波包分解系数集合；x4,n表示DD2频段信号的小波包分解的各个系数；
步骤4.2：用公式(5)表示AA2频段信号的能量，用符号E1表示；用公式将(6)表示DA2频段信号的能量，用符号E2表示；用公式(7)表示AD2频段信号的能量，用符号E3表示；用公式(8)表示DD2频段信号的能量，用符号E4表示；
步骤4.3：通过公式(9)得到AA2频段信号的能谱系数P1；
P1＝100*(E1/(E1+E2+E3+E4)) (9)
步骤五、在步骤四的操作基础上，以A个发射信号为滑动窗口，即按采集次序依次选取A个声发射信号为一个窗口，计算每个窗口中P1的平均值，并绘制P1随窗口数的分布图，横坐标为窗口的序号；纵坐标为对应窗口中P1的平均值；
步骤六、根据P1随窗口数的分布图中，P1的变化状态监测所述待监测脆性材料损伤演化过程；监测所述待监测脆性材料损伤演化过程的具体方法为：当窗口M之后连续B个窗口中P1的平均值大于该窗口及之前连续B个窗口中P1的平均值，则判断待监测脆性材料在M窗口由裂纹萌生阶段进入裂纹生长阶段；当窗口N，N>M，之后连续B个窗口中P1的平均值大于该窗口及之前连续B个窗口中P1的平均值，则判断待监测脆性材料在N窗口由裂纹生长阶段进入裂纹汇合阶段，裂纹汇合阶段持续到待监测脆性材料最终的破坏失效。
2.如权利要求1所述的一种使用信号能谱系数监测脆性材料损伤演化过程的方法，其特征在于：所述脆性材料为陶瓷类材料、岩石和混凝土。
3.如权利要求1所述的一种使用信号能谱系数监测脆性材料损伤演化过程的方法，其特征在于：所述A为20；所述B为10。
CN201810291038.7A 2018-04-03 2018-04-03 一种使用信号能谱系数监测脆性材料损伤演化过程的方法 Active CN108760890B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810291038.7A CN108760890B (zh) 2018-04-03 2018-04-03 一种使用信号能谱系数监测脆性材料损伤演化过程的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810291038.7A CN108760890B (zh) 2018-04-03 2018-04-03 一种使用信号能谱系数监测脆性材料损伤演化过程的方法

Publications (2)

Publication Number Publication Date
CN108760890A true CN108760890A (zh) 2018-11-06
CN108760890B CN108760890B (zh) 2021-03-09

Family

ID=63981116

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810291038.7A Active CN108760890B (zh) 2018-04-03 2018-04-03 一种使用信号能谱系数监测脆性材料损伤演化过程的方法

Country Status (1)

Country Link
CN (1) CN108760890B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112525672A (zh) * 2021-01-21 2021-03-19 东北石油大学 基于裂纹扩展能量演化的页岩储层脆性测试装置及方法
CN113686968A (zh) * 2021-08-10 2021-11-23 北京科技大学 实时检测薄膜裂纹间距的方法和装置
CN117191952A (zh) * 2023-11-06 2023-12-08 中冶建筑研究总院有限公司 基于声发射信号小波包分解频带能量谱的疲劳损伤识别与寿命预测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102062747B (zh) * 2010-11-15 2012-11-21 湖南大学 基于压电陶瓷和小波包分析的钢管混凝土构件钢管壁剥离监测方法
CN102854252A (zh) * 2012-09-10 2013-01-02 北京理工大学 一种检测金属材料疲劳状态的方法和系统
CN104965026A (zh) * 2015-05-25 2015-10-07 北京理工大学 一种使用声发射信号能量特征值预测脆性材料临界失稳的方法
CN106370730A (zh) * 2016-08-25 2017-02-01 中国科学院武汉岩土力学研究所 一种基于声发射技术精确测量脆性材料损伤阈值的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102062747B (zh) * 2010-11-15 2012-11-21 湖南大学 基于压电陶瓷和小波包分析的钢管混凝土构件钢管壁剥离监测方法
CN102854252A (zh) * 2012-09-10 2013-01-02 北京理工大学 一种检测金属材料疲劳状态的方法和系统
CN104965026A (zh) * 2015-05-25 2015-10-07 北京理工大学 一种使用声发射信号能量特征值预测脆性材料临界失稳的方法
CN106370730A (zh) * 2016-08-25 2017-02-01 中国科学院武汉岩土力学研究所 一种基于声发射技术精确测量脆性材料损伤阈值的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张立新: "《基于小波分析及神经网络的混凝土内部缺陷智能识别》", 30 June 2016 *
梁天成 等: "利用声发射和波速变化判定岩石损伤状态", 《中国地震》 *
赖于树 等: "受载混凝土破坏全过程声发射信号频带能量特征", 《振动与冲击》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112525672A (zh) * 2021-01-21 2021-03-19 东北石油大学 基于裂纹扩展能量演化的页岩储层脆性测试装置及方法
CN113686968A (zh) * 2021-08-10 2021-11-23 北京科技大学 实时检测薄膜裂纹间距的方法和装置
CN113686968B (zh) * 2021-08-10 2022-08-02 北京科技大学 实时检测薄膜裂纹间距的方法和装置
CN117191952A (zh) * 2023-11-06 2023-12-08 中冶建筑研究总院有限公司 基于声发射信号小波包分解频带能量谱的疲劳损伤识别与寿命预测方法
CN117191952B (zh) * 2023-11-06 2024-01-23 中冶建筑研究总院有限公司 基于声发射信号小波包分解频带能量谱的疲劳损伤识别与寿命预测方法

Also Published As

Publication number Publication date
CN108760890B (zh) 2021-03-09

Similar Documents

Publication Publication Date Title
CN108760890A (zh) 一种使用信号能谱系数监测脆性材料损伤演化过程的方法
CN104965026B (zh) 一种使用声发射信号能量特征值预测脆性材料临界失稳的方法
US8275548B2 (en) Method and apparatus for detecting lightning activity
CN105676268B (zh) 一种基于声音信号波形变化特征的应变型岩爆预警方法
CN102854252B (zh) 一种检测金属材料疲劳状态的方法和系统
Lin et al. HHT-based AE characteristics of natural fatigue cracks in rotating shafts
CN103546853A (zh) 一种基于短时傅里叶变换的扬声器异常音检测方法
CN106841403A (zh) 一种基于神经网络的声学玻璃缺陷检测方法
CN103901111A (zh) 风力发电机组叶片的无损检测系统及方法
CN109085244A (zh) 一种基于压电阵列的非线性Lamb波结构疲劳损伤层析成像方法
CN102890750B (zh) 运输安全记录仪数据分析方法
CN108519436B (zh) 一种基于声发射信号统计分析的混凝土损伤评估方法
JP2007243338A (ja) ログ分析装置、ログ分析プログラム、および記録媒体
CN116257540A (zh) 一种地理环境数据的管理方法和管理系统
Wang et al. Weak fault diagnosis of rolling bearing under variable speed condition using IEWT-based enhanced envelope order spectrum
CN107255677A (zh) 一种基于压电智能传感器的大坝安全损伤监测系统
Wang et al. The optimal wavelet threshold de-nosing method for acoustic emission signals during the medium strain rate damage process of concrete
CN116796571A (zh) 基于声发射特征的岩石构造应力预测系统
Wu et al. A novel wireless acoustic emission sensor system for distributed wooden structural health monitoring
CN116519112A (zh) 一种基于激光测振的超声刀头振动能量效率评价方法
CN109782064B (zh) 一种风电场输出阻抗频率特性测试分析方法
CN103345921B (zh) 基于多特征的夜间睡眠声信号分析方法
CN106353406A (zh) 一种风电机组螺栓断裂监测装置
CN117783294B (zh) 一种基于声发射能量熵的木材损伤动态检测方法与系统
Montano Low frequency noise and infrasound: A new method to determine the specific sound from the total sound; a plausible statistical algorithm for use in Legal Noise Assessment

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant