CN108195937B - 一种基于导波的损伤概率成像方法 - Google Patents
一种基于导波的损伤概率成像方法 Download PDFInfo
- Publication number
- CN108195937B CN108195937B CN201711230877.XA CN201711230877A CN108195937B CN 108195937 B CN108195937 B CN 108195937B CN 201711230877 A CN201711230877 A CN 201711230877A CN 108195937 B CN108195937 B CN 108195937B
- Authority
- CN
- China
- Prior art keywords
- damage
- monitoring
- imaging
- monitoring path
- path
- 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.)
- Active
Links
Images
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/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于导波的损伤概率成像方法,涉及结构健康监测技术领域,其包括下列步骤:(1)在待测结构上布置压电片组成N条压电监测路径;(2)获取每条监测路径的基准信号;(3)获取每条监测路径的当前信号;(4)计算每条监测路径的损伤指数;(5)对监测区域建立坐标系,划分成像网格,并计算每个成像网格的损伤概率成像值,其中损伤概率成像值越大,表明了损伤存在的概率越大。本发明提供的损伤概率成像方法,具有:损伤成像的图像清晰,损伤定位精度高,可以对多损伤进行成像识别,方法简单易行,可应用于复杂结构的结构健康监测。
Description
技术领域
本发明属于结构健康监测技术领域,涉及一种基于导波的损伤概率成像方法。
背景技术
近年来,基于导波的结构健康监测技术受到了广泛的关注,成为了一个研究热点。这是由于导波能够在结构中传播较长的距离,并且对裂纹、腐蚀和分层等损伤敏感。在基于导波的结构健康监测中,损伤成像算法由于能用图像直观地识别和定位损伤,从而受到了重视,并得到了广泛研究。其中,基于椭圆分布概率的损伤概率成像方法由于能利用稀疏传感器网络产生较好质量的损伤图像,且适用于复杂结构,而得到了很多研究。在此方法中,对于每个激励-传感路径,损伤存在分布概率是一个线性衰减的椭圆分布概率,即成像点离激励- 传感路径越远,损伤存在的概率越小。并且,对于每个激励-传感路径,此方法的损伤存在概率,无论损伤处于何位置,总是在激励-传感的直达路径上最大。这样的损伤分布概率是一种不准确的分布概率,与事实不符合,会造成损伤定位精度的降低。另一方面,控制椭圆分布概率影响区域大小的比例因子,需要靠经验进行选择确定。以上,限制了此方法在实际监测中的应用。
发明内容
本发明的目的是:为了克服基于椭圆分布概率的损伤概率成像方法的不足,本发明提供了一种基于椭圆环分布概率的损伤概率成像方法,可以改善损伤存在分布概率的准确性,提高损伤定位的精度。
本发明的技术方案是:一种基于导波的损伤概率成像方法,包括下列步骤: (1)在待测结构上布置压电片组成N条压电监测路径;(2)获取每条监测路径的基准信号;(3)获取每条监测路径的当前信号;(4)计算每条监测路径的损伤指数;(5)对监测区域建立坐标系,划分成像网格,并计算每个成像网格的损伤概率成像值,其中损伤概率成像值越大,表明了损伤存在的概率越大。
本发明具体包括下列步骤:
(1)在待测结构4上布置压电片阵列3组成N条压电监测路径;
(2)获取每条监测路径的基准信号Sbi (t)(i=1,2,...,N):由波形发生器1通过如下公式产生波形,并通过功率放大器2激励压电片阵列3组成的监测路径中的激励器;由每条监测路径中的传感器接收信号,并通过电荷放大器5、数据采集处理系统6存储每条监测路径的基准信号Sbi (t);
式中:
A—信号的幅度调制;
fc—信号中心频率;
n—信号波峰个数;
H—Heaviside阶梯函数;
t—时间。
(3)与步骤(2)相同的方法,获取每条监测路径的当前信号Sci (t)(i=1,2,...,N);
(4)计算每条监测路径的损伤指数DIi(i=1,2,...,N);
式中:
ti—积分起始时间;
tf—积分结束时间;
(5)在监测区域建立笛卡尔平面坐标系,并将监测区域划分为均匀大小的成像网格,得到压电片坐标和各网格中心的坐标。对于中心坐标为(x,y)的网格,其损伤概率成像值为:
其中DIi是第i条监测路径的损伤指数,N是监测路径的总数,Di(x,y)是中心坐标为(x,y)的网格分别到第i条监测路径的激励器和传感器的距离之和,可以表示为:
Di(x,y)=Da,i(x,y)+Ds,i(x,y)
其中Da,i(x,y)和Ds,i(x,y)分别是中心坐标为(x,y)的网格到第i条监测路径的激励器和传感器的距离。
Wi[Di(x,y)]是第i条监测路径的损伤分布概率的权重函数,可以表示为:
其中D1i是第i条监测路径的估计的最近损伤位置,可以表达为:
D1i=c1ti
其中c1为估计的导波波速下限,ti为第i条监测路径的损伤散射信号的到达时间。D2i是第i条监测路径的估计的最远损伤位置,可以表达为:
D2i=c2ti
其中c2为估计的导波波速上限。
上述得到的各成像网格的损伤概率成像值,其中成像值越大,表明该网格处存在损伤的可能性越大。
本发明相比现有技术的优点在于:损伤成像的图像清晰,损伤定位精度高,可以对多损伤进行成像识别,方法简单易行。
附图说明
图1是基于导波的损伤概率成像方法实施例示意图;
图中,1-波形发生器,2-功率放大器,3-压电片阵列,4-待测结构,5-电荷放大器,6-数据采集处理系统。
图2是本发明一具体实施例中板尺寸、压电片布置及坐标系示意图。
图3为使用本发明的损伤概率成像图。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步说明。
参照图1,现有基于导波的结构损伤监测基本由波形发生器1,功率放大器 2,压电片阵列3,待测结构4,电荷放大器5,数据采集处理系统6组成。波形发生器1通过导线与功率放大器2连接;功率放大器2通过导线与压电片阵列3 组成的监测路径中的激励器连接;监测路径中的传感器通过导线与电荷放大器5 连接;电荷放大器5通过导线与数据采集处理系统连接。在监测过程中,由系统在待测结构4中产生导波,并采集在待测结构4监测区域传播的导波,对采集到的导波信号进行处理,从而得到损伤概率成像图。具体步骤为:
(1)在待测结构4上布置压电片阵列3组成N条压电监测路径;
(2)获取每条监测路径的基准信号Sb (t):由波形发生器1产生波形,并通过功率放大器2放大后,激励压电片阵列3组成的监测路径中的激励器;由每条监测路径中的传感器接收信号,并通过电荷放大器5、数据采集处理系统6存储每条监测路径的基准信号Sb(t);
(3)与步骤(2)相同的方法,获取每条监测路径的当前信号Sc (t);
(4)计算每条监测路径的损伤指数DI;
(5)对监测区域建立坐标系,划分成像网格,并计算每个成像网格的损伤概率成像值。
按照上述步骤,对一碳纤维复合材料板进行损伤监测。复合材料板的材料为:T700/BA9916,铺层为:[45/0/-45/90/0/45/0/-45/0/45/90/-45]S,尺寸为: 450mm×500mm×3mm。在板上布置了8个压电片。复合材料板尺寸,压电片布置及坐标系建立,如图2所示。由于导波传播的互易性,此压电片网络共组成28条压电监测路径。信号激励频率为90KHz,信号采样率为10MHz。监测开始前,先采集了一组基准信号。然后通过冲击引入了两处损伤,再采集了一组当前信号。利用基准信号和当前信号,计算了每条监测路径的损伤指数。然后通过损伤概率成像公式,得到了损伤成像图,如图3所示,其中损伤成像图采用了阈值化处理,对小于阈值1.05的成像值设置为了零。
Claims (1)
1.一种基于导波的损伤概率成像方法,其特征为所述方法包括下列步骤:
1.1)在待测结构上布置压电片组成N条压电监测路径;
1.2)获取每条监测路径的基准信号;
1.3)获取每条监测路径的当前信号;
1.4)计算每条监测路径的损伤指数;
1.5)对监测区域建立坐标系,划分成像网格,并计算每个成像网格的损伤概率成像值,其中损伤概率成像值越大,表明了损伤存在的概率越大;所述方法包括下列步骤:
2.1)在待测结构(4)上布置压电片阵列(3)组成N条压电监测路径;
2.2)获取每条监测路径的基准信号Sbi (t)(i=1,2,...,N):由波形发生器(1)通过如下公式产生波形,并通过功率放大器(2)激励压电片阵列(3)组成的监测路径中的激励器;由每条监测路径中的传感器接收信号,并通过电荷放大器(5)、数据采集处理系统(6)存储每条监测路径的基准信号Sbi (t);
式中:
A—信号的幅度调制;
fc—信号中心频率;
n—信号波峰个数;
H—Heaviside阶梯函数;
t—时间;
2.3)与步骤2.2)相同的方法,获取每条监测路径的当前信号Sci (t)(i=1,2,...,N);
2.4)计算每条监测路径的损伤指数DIi(i=1,2,...,N);
式中:
ti—积分起始时间;
tf—积分结束时间;
2.5)在监测区域建立笛卡尔平面坐标系,并将监测区域划分为均匀大小的成像网格,得到压电片坐标和各网格中心的坐标;对于中心坐标为(x,y)的网格,其损伤概率成像值为:
其中DIi是第i条监测路径的损伤指数,N是监测路径的总数,Di(x,y)是中心坐标为(x,y)的网格分别到第i条监测路径的激励器和传感器的距离之和,表示为:
Di(x,y)=Da,i(x,y)+Ds,i(x,y)
其中Da,i(x,y)和Ds,i(x,y)分别是中心坐标为(x,y)的网格到第i条监测路径的激励器和传感器的距离;
Wi[Di(x,y)]是第i条监测路径的损伤分布概率的权重函数,表示为:
其中D1i是第i条监测路径的估计的最近损伤位置,表达为:
D1i=c1ti
其中c1为估计的导波波速下限,ti为第i条监测路径的损伤散射信号的到达时间;D2i是第i条监测路径的估计的最远损伤位置,表达为:
D2i=c2ti
其中c2为估计的导波波速上限;
上述得到的各成像网格的损伤概率成像值,其中成像值越大,表明该网格处存在损伤的可能性越大。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711230877.XA CN108195937B (zh) | 2017-11-29 | 2017-11-29 | 一种基于导波的损伤概率成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711230877.XA CN108195937B (zh) | 2017-11-29 | 2017-11-29 | 一种基于导波的损伤概率成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108195937A CN108195937A (zh) | 2018-06-22 |
CN108195937B true CN108195937B (zh) | 2020-09-18 |
Family
ID=62573396
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711230877.XA Active CN108195937B (zh) | 2017-11-29 | 2017-11-29 | 一种基于导波的损伤概率成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108195937B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109900804B (zh) * | 2019-03-25 | 2021-11-09 | 中国特种设备检测研究院 | 一种基于超声导波的金属材料裂纹定量监测方法 |
CN110849724A (zh) * | 2019-11-23 | 2020-02-28 | 福州大学 | 一种用于装配式混凝土剪力墙损伤识别的概率成像方法 |
CN111208206B (zh) * | 2020-02-20 | 2022-11-25 | 北京博感科技有限公司 | 一种温度变化环境下的复合材料超声导波结构健康监测方法 |
CN112179990A (zh) * | 2020-09-15 | 2021-01-05 | 昆明理工大学 | 一种基于ToF损伤因子的碳纤维复合材料疲劳损伤概率成像方法 |
CN112660417B (zh) * | 2020-12-25 | 2022-03-25 | 湖南航天机电设备与特种材料研究所 | 一种用于飞行器结构部件的结构损伤诊断方法及诊断系统 |
CN113702505A (zh) * | 2021-08-19 | 2021-11-26 | 河北工程大学 | 一种垃圾填埋场hdpe膜损伤定位系统及方法 |
CN114184680A (zh) * | 2021-12-21 | 2022-03-15 | 中国飞机强度研究所 | 一种基于权重分布修正的导波损伤概率成像方法 |
CN114460175B (zh) * | 2022-02-28 | 2024-03-15 | 西北工业大学 | 一种薄壁结构损伤检测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102043016A (zh) * | 2010-11-05 | 2011-05-04 | 上海交通大学 | 基于兰姆波的自主式损伤识别成像方法 |
CN102331455A (zh) * | 2011-08-18 | 2012-01-25 | 中国飞机强度研究所 | 一种基于主动Lamb波损伤指数的工程结构损伤监测方法 |
CN102928511A (zh) * | 2012-10-24 | 2013-02-13 | 西安交通大学 | 一种基于rapid层析的机械结构损伤无损识别方法 |
CN105488795A (zh) * | 2015-11-26 | 2016-04-13 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种复合材料损伤识别方法 |
-
2017
- 2017-11-29 CN CN201711230877.XA patent/CN108195937B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102043016A (zh) * | 2010-11-05 | 2011-05-04 | 上海交通大学 | 基于兰姆波的自主式损伤识别成像方法 |
CN102331455A (zh) * | 2011-08-18 | 2012-01-25 | 中国飞机强度研究所 | 一种基于主动Lamb波损伤指数的工程结构损伤监测方法 |
CN102928511A (zh) * | 2012-10-24 | 2013-02-13 | 西安交通大学 | 一种基于rapid层析的机械结构损伤无损识别方法 |
CN105488795A (zh) * | 2015-11-26 | 2016-04-13 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种复合材料损伤识别方法 |
Non-Patent Citations (4)
Title |
---|
Lamb波谱元法建模及多损伤定量化监测方法研究;刘梦龙;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20150215;全文 * |
便携式主动Lamb波结构健康监测集成系统设计;王星等;《压电与声光》;20150228;第37卷(第1期);138-141 * |
基于加权PCR证据融合的多传感器Lamb波损伤成像;刘学君等;《南京航空航天大学学报》;20170831;第49卷(第4期);479-486 * |
基于多维阵列和空间滤波器的损伤无波速成像定位方法;刘彬等;《复合材料学报》;20140630;第31卷(第3期);835-844 * |
Also Published As
Publication number | Publication date |
---|---|
CN108195937A (zh) | 2018-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108195937B (zh) | 一种基于导波的损伤概率成像方法 | |
CN106525968B (zh) | 基于子区域的损伤概率成像定位方法 | |
CN102323337B (zh) | 一种采用合成波阵面激励的工程结构损伤主动监测方法 | |
Su et al. | Predicting delamination of composite laminates using an imaging approach | |
CN102043016B (zh) | 基于兰姆波的自主式损伤识别成像方法 | |
CN111122700B (zh) | 一种提高激光超声saft缺陷定位速度的方法 | |
CN104965025A (zh) | 一种基于Lamb波信号相关系数的多区域损伤检测方法 | |
CN102353718A (zh) | 用于复合材料板结构损伤监测的Lamb波损伤概率成像方法 | |
CN111208206B (zh) | 一种温度变化环境下的复合材料超声导波结构健康监测方法 | |
CN106770664A (zh) | 一种基于全聚焦成像算法改进边缘缺陷检测的方法 | |
CN105784853A (zh) | 一种工程结构方向性损伤主动监测与评估方法 | |
CN109632963B (zh) | 一种基于时不变特征信号构建的结构损伤四维成像方法 | |
CN112179990A (zh) | 一种基于ToF损伤因子的碳纤维复合材料疲劳损伤概率成像方法 | |
CN109696480B (zh) | 一种基于改进时间反转算法的玻璃纤维复合材料声发射源定位成像方法 | |
CN110687208A (zh) | 一种基于双曲线定位的无基准Lamb波损伤监测方法 | |
CN111521691B (zh) | 一种基于时间反转加权分布的复合材料Lamb波损伤成像方法 | |
CN102998369B (zh) | 一种二维损伤定量化监测方法 | |
CN105372327B (zh) | 基于十字阵和波数扫描滤波器的在线损伤成像方法 | |
CN105488795A (zh) | 一种复合材料损伤识别方法 | |
CN103792287A (zh) | 一种基于Lamb波的大面积结构损伤检测方法 | |
CN104181235A (zh) | 一种基于虚拟时间反转的声发射被动成像方法和装置 | |
CN110702792A (zh) | 一种基于深度学习的合金组织超声检测分类方法 | |
CN104730152A (zh) | 一种基于分形维数的复合材料结构件裂纹损伤监测方法 | |
CN103412053A (zh) | 一种基于双声发射传感阵列和波束形成的无需波速的声发射源定位方法 | |
He et al. | A fast acoustic emission beamforming localization method based on Hilbert curve |
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 |