CN110715929B - Distributed strain micro crack detection system and method based on stacking self-encoder - Google Patents
Distributed strain micro crack detection system and method based on stacking self-encoder Download PDFInfo
- Publication number
- CN110715929B CN110715929B CN201910973652.6A CN201910973652A CN110715929B CN 110715929 B CN110715929 B CN 110715929B CN 201910973652 A CN201910973652 A CN 201910973652A CN 110715929 B CN110715929 B CN 110715929B
- Authority
- CN
- China
- Prior art keywords
- strain
- encoder
- subsequences
- subsequence
- distributed
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/84—Systems specially adapted for particular applications
- G01N21/88—Investigating the presence of flaws or contamination
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种基于堆叠自编码器的分布式应变裂缝检测系统及方法,利用深度神经网络良好的特征表征能力,将裂缝检测视为一个二分类问题,构建一个基于堆叠自编码器的深度神经网络,实现结构体应变子序列的裂缝与非裂缝分类。本发明方法在实验室钢梁上可以准确并且无遗漏地检测张口宽度为32μm的微小裂缝,为结构体表面分布式应变裂缝检测提供了一种具有良好噪声鲁棒性的解决方案。
The invention discloses a distributed strain crack detection system and method based on a stacked self-encoder, which utilizes the good feature representation ability of a deep neural network, regards crack detection as a two-class problem, and constructs a depth based on the stacked self-encoder. A neural network for crack and non-crack classification of structural strain subsequences. The method of the invention can accurately and without omission detect tiny cracks with an opening width of 32 μm on a laboratory steel beam, and provides a solution with good noise robustness for the detection of distributed strain cracks on the surface of the structure.
Description
技术领域technical field
本发明属于模式识别领域,具体涉及一种基于堆叠自编码器的分布式应变微小裂缝检测系统及方法。The invention belongs to the field of pattern recognition, and in particular relates to a distributed strain micro-crack detection system and method based on a stacked self-encoder.
背景技术Background technique
裂缝检测一直是结构健康监测领域中的重要课题。裂缝检测方法包含了人工观测的方法和无损检测方法。人工观测的方法需要专门的维护人员使用专业的工具进行定期检查,该方法效率低、主观性强。无损检测方法主要通过超声波、X射线、探地雷达以及摄像机等获得的数据对结构体裂缝进行检测。这些传感器都是点对点传感器,无法针对结构体的整体数据进行测量,容易遗漏裂缝。Crack detection has always been an important topic in the field of structural health monitoring. Crack detection methods include manual observation methods and non-destructive testing methods. The manual observation method requires specialized maintenance personnel to use professional tools to conduct regular inspections, which is inefficient and highly subjective. Non-destructive testing methods mainly detect structural cracks through the data obtained by ultrasonic, X-ray, ground penetrating radar and cameras. These sensors are all point-to-point sensors, which cannot measure the overall data of the structure and are prone to missing cracks.
发明内容SUMMARY OF THE INVENTION
本发明的目的在于提供一种基于堆叠自编码器的分布式应变微小裂缝检测系统及方法,以克服上述现有技术存在的缺陷,本发明能够对结构体的整体应变数据进行有效的采集,结构体整体应变数据进行裂缝检测,显著改善了裂缝检测的正确率,提升裂缝检测的检测效果,为结构健康监测提供了一种高效的裂缝检测方案。The purpose of the present invention is to provide a distributed strain micro-crack detection system and method based on stacked self-encoders, so as to overcome the above-mentioned defects in the prior art, and the present invention can effectively collect the overall strain data of the structure. The overall strain data is used for crack detection, which significantly improves the accuracy of crack detection, improves the detection effect of crack detection, and provides an efficient crack detection solution for structural health monitoring.
为达到上述目的,本发明采用如下技术方案:To achieve the above object, the present invention adopts the following technical solutions:
一种基于堆叠自编码器的分布式应变微小裂缝检测系统,包括:应变序列采集模块,用于对结构体表面的分布式应变进行采集;应变序列预处理模块:用于将采集所得的应变序列进行z-score标准化并截取为应变子序列;基于堆叠自编码器的特征自学习与表征模块:用于提取所划分的应变子序列的特征;Softmax分类识别模块,用于对提取到的子序列特征进行二分类,判别每一个子序列属于裂缝子序列和非裂缝子序列的概率;A distributed strain micro-crack detection system based on stacked self-encoders, comprising: a strain sequence acquisition module for acquiring distributed strain on the surface of a structure; a strain sequence preprocessing module: for collecting the acquired strain sequence Standardize z-score and truncate into strain subsequences; feature self-learning and characterization module based on stacked autoencoder: used to extract the features of the divided strain subsequences; Softmax classification and recognition module, used to extract the subsequences The features are classified into two categories to determine the probability that each subsequence belongs to a cracked subsequence and a non-cracked subsequence;
进一步地,应变序列采集模块:将光纤传感器敷设于结构体表面,使用基于BOTDA的分布式光纤传感系统对结构体表面的分布式应变进行采集;应变序列预处理模块包括:Further, the strain sequence acquisition module: lays the optical fiber sensor on the surface of the structure, and uses the distributed optical fiber sensing system based on BOTDA to collect the distributed strain on the surface of the structure; the strain sequence preprocessing module includes:
进一步地,z-score标准化模块和滑窗截取模块,z-score标准化模块将应变序列标准化为0均值1标准差的数据。滑窗模块通过长度为21,步长为1的滑动窗口将标准化后的应变序列截取了一组长度都是21的应变子序列。基于堆叠自编码器的特征自学习与表征模块:由3个自动编码器模块构成,将3个自动编码器模块的编码部分用作特征表征。Further, the z-score normalization module and the sliding window interception module, the z-score normalization module normalizes the strain series to data with 0 mean and 1 standard deviation. The sliding window module intercepts a set of strain subsequences with a length of 21 from the normalized strain sequence through a sliding window with a length of 21 and a step size of 1. Feature Self-Learning and Representation Module Based on Stacked Auto-Encoder: It consists of 3 auto-encoder modules, and the encoded part of the 3 auto-encoder modules is used as feature representation.
进一步地,基于堆叠自编码器的特征自学习与表征模块:由3个自动编码器模块构成,用于提取所划分的应变子序列的特征,自动编码器模块对于输入数据x,特征h与输出之间的关系可以表示为fθ(·)和gθ'(·)两个函数,具体如下所示:Further, the feature self-learning and characterization module based on stacked autoencoders: composed of 3 autoencoder modules, used to extract the features of the divided strain subsequences, the autoencoder module for input data x, feature h and output The relationship between can be expressed as two functions f θ ( ) and g θ' ( ), as follows:
h=fθ(x)=sf(Wx+bh)h=f θ (x)=s f (Wx+b h )
其中,W和bh分别为输入数据x和特征h之间的连接矩阵和偏置向量,WT和bv分别为特征h和输出之间的连接矩阵和偏置向量,WT为W的转置矩阵。sf和sg是激活函数。Among them, W and b h are the connection matrix and bias vector between the input data x and feature h, respectively, W T and b v are the feature h and output, respectively Between the connection matrix and the bias vector, W T is the transpose matrix of W. s f and s g are activation functions.
基于堆叠自编码器的特征自学习与表征模块的过程具体如下所示:The process of feature self-learning and representation module based on stacked autoencoder is as follows:
hi=sf(Wixi+bh,i)h i =s f (W i x i +b h,i )
其中,xi,hi,分别为第i(i=1,2,…,n)个自动编码器模块的输入数据x,特征h与输出Wi和Wi T为其连接权重矩阵,bh,i和bv,i为其偏置向量,xi=hi-1,sf和sg为激活函数。Among them, x i , h i , are the input data x, feature h and output of the i-th (i=1,2,...,n) autoencoder module, respectively Wi and Wi T are their connection weight matrices, b h,i and b v,i are their bias vectors, xi =h i-1 , s f and s g are activation functions.
一种基于堆叠自编码器的分布式应变裂缝检测方法,包括以下步骤:A distributed strain crack detection method based on stacked autoencoders, comprising the following steps:
步骤1:应变序列采集;Step 1: Strain sequence acquisition;
步骤2:使用z-score标准化将采集到的应变序列进行标准化,使用长度为21,步长为1的滑动窗口截取应变序列,得到应变子序列,将应变子序列按照截取位置进行标记;Step 2: Use z-score normalization to standardize the collected strain sequence, use a sliding window with a length of 21 and a step size of 1 to intercept the strain sequence to obtain a strain subsequence, and mark the strain subsequence according to the interception position;
步骤3:使用基于堆叠自编码器的神经网络自动学习表征应变子序列的特征;Step 3: Automatically learn features representing strain subsequences using a stacked autoencoder-based neural network;
步骤4:采用Softmax分类器实现对提取到的应变子序列的特征进行二分类,完成裂缝检测;Step 4: Use the Softmax classifier to implement binary classification of the features of the extracted strain subsequences to complete crack detection;
进一步地,步骤1中应变序列采集具体过程为:将光纤传感器通过环氧树脂粘附于结构体表面,光纤的两端接于基于BOTDA的分布式光纤传感系统,基于BOTDA的分布式光纤传感系统通过两个光源——泵浦光和探测光测量光纤的布里渊频移,通过布里渊频移与应变的线性关系得到结构体表面的分布式应变。Further, the specific process of strain sequence acquisition in step 1 is as follows: the optical fiber sensor is adhered to the surface of the structure through epoxy resin, and the two ends of the optical fiber are connected to the BOTDA-based distributed optical fiber sensing system, and the BOTDA-based distributed optical fiber transmission system. The sensing system measures the Brillouin frequency shift of the fiber through two light sources—pump light and probe light, and obtains the distributed strain on the surface of the structure through the linear relationship between the Brillouin frequency shift and strain.
进一步地,步骤2中应变序列处理的具体过程为:Further, the specific process of the strain sequence processing in
步骤2.1:将采集到的应变序列减去其均值,除以其标准差,得到0均值1标准差的数据。Step 2.1: Subtract the mean value of the acquired strain series and divide it by its standard deviation to obtain data with 0 mean and 1 standard deviation.
步骤2.2:使用长度为21,步长为1的滑动窗口沿着采集到的应变序列进行滑动,将应变序列截取为一组长度为21的应变子序列。Step 2.2: Use a sliding window with a length of 21 and a step of 1 to slide along the acquired strain sequence, and cut the strain sequence into a set of strain subsequences with a length of 21.
步骤2.3:将所得应变子序列按照截取位置标记标签,将以裂缝处为中心截取的应变子序列标记为裂缝子序列,并将其左侧3个和右侧4个应变子序列标记为裂缝子序列,其余标记为非裂缝子序列。Step 2.3: Label the obtained strain subsequence according to the interception position, mark the strain subsequence intercepted at the crack as the crack subsequence, and mark the left 3 strain subsequences and the right 4 strain subsequences as crack subsequences sequences, and the rest are marked as non-cracked subsequences.
进一步地,步骤3中使用基于堆叠自编码器的神经网络自动学习表征应变子序列特征的具体过程如下:Further, the specific process of using the stacked autoencoder-based neural network to automatically learn the characteristics of the strain subsequence in
步骤3.1:模型初始化,确定模型的层数与神经元个数。随机初始化模型中的连接权重矩阵和偏置向量。输入层的神经元个数等于21,为应变子序列的长度。Step 3.1: Model initialization, determine the number of layers and neurons of the model. Randomly initialize the connection weight matrix and bias vector in the model. The number of neurons in the input layer is equal to 21, which is the length of the strain subsequence.
步骤3.2:预训练堆叠自编码器,堆叠自编码器由3个自动编码器构成,用得到的应变子序列预训练每个自动编码器。预训练自动编码器的损失函数为输入与输出之间的均方误差,具体如下:Step 3.2: Pretrain the stacked autoencoder, which consists of 3 autoencoders, and pretrain each autoencoder with the resulting strain subsequences. The loss function of the pretrained autoencoder is the mean squared error between the input and the output, as follows:
其中,x为输入的应变子序列,为自动编码器输出的重构数据,M为所有输入的应变子序列的数量,Xm、分别是输入模型的第m条应变子序列和对应的输出重构的第m条子序列。where x is the input strain subsequence, is the reconstructed data output by the autoencoder, M is the number of all input strain subsequences, X m , are the mth strain subsequence of the input model and the mth subsequence of the corresponding output reconstruction, respectively.
进一步地,步骤4中采用Softmax分类器实现对应变子序列的分类,具体方法为:Further, in step 4, the Softmax classifier is used to realize the classification of the strain subsequences, and the specific method is:
步骤4.1:构建Softmax分类器,对于给定的输入z,用假设函数hδ(z)针对每一个类别l估算出概率值p(y=l|z),l∈{0,1},假设函数hδ(z)输出一个t维的向量表示这t个估计的概率值,t=2,假设函数hδ(z)如下:Step 4.1: Build a Softmax classifier. For a given input z, use the hypothesis function h δ (z) to estimate the probability value p(y=l|z) for each category l, l∈{0,1}, suppose The function h δ (z) outputs a t-dimensional vector representing the t estimated probability values, t=2, assuming that the function h δ (z) is as follows:
其中,δ1,δ2是Softmax分类器的全部参数,z(i)为输入,y(i)为输出,Softmax分类器将z分为类别l的概率为:Among them, δ 1 , δ 2 are all parameters of the Softmax classifier, With z (i) as input and y (i) as output, the probability that the Softmax classifier classifies z into class l is:
其中,z(i)为输入,y(i)为输出;where z (i) is the input and y (i) is the output;
步骤4.2:预训练Sofmax,将应变子序列输入预训练后的堆叠自编码器,得到输出的特征z(i),以z(i)及其标签类别y(i)预训练Softmax,损失函数为交叉熵函数,具体如下:Step 4.2: Pre-train Sofmax, input the strain subsequence into the pre-trained stacked autoencoder, get the output feature z (i) , pre-train Softmax with z (i) and its label category y (i) , the loss function is Cross entropy function, as follows:
其中,为Softmax分类器的全部参数,为输出的类别概率,λ1为Softmax中连接权重矩阵和偏置向量正则项的权重系数,M为输入应变子序列的总个数,K为类别数,为2。in, are all parameters of the Softmax classifier, is the output category probability, λ 1 is the weight coefficient of the regular term connecting the weight matrix and the bias vector in Softmax, M is the total number of input strain subsequences, and K is the number of categories, which is 2.
步骤4.3:微调,堆叠自编码器的编码部分后接Softmax分类器,可以使其具有分类功能。利用预训练得到的应变子序列微调堆叠自编码器的编码部分与Softmax分类器整体结构的连接权重矩阵和偏置向量。微调时的损失函数为交叉损失函数,具体如下:Step 4.3: Fine-tuning, the encoding part of the stacked autoencoder is followed by a Softmax classifier, which can make it have a classification function. Using the pretrained strain subsequences to fine-tune the connection weight matrix and bias vector between the encoding part of the stacked autoencoder and the overall structure of the Softmax classifier. The loss function during fine-tuning is a cross loss function, as follows:
其中,ω为堆叠自编码器中的连接权重矩阵和偏置向量,Θ为ω和δ,λ2为堆叠自编码器中连接权重矩阵和偏置向量正则项的权重系数。where ω is the connection weight matrix and bias vector in the stacked autoencoder , Θ is ω and δ, and λ2 is the weight coefficient of the regular term of the connection weight matrix and bias vector in the stacked autoencoder.
步骤4.4:Softmax分类器接收堆叠自编码器输出的特征作为其输入,输出应变子序列的类别0或1,0表示非裂缝,1表示裂缝;对于堆叠自编码器输出的特征z(i),选择概率p(y(i)=l|z(i);δ)最大的类别l作为该特征对应的类别。将分布式应变子序列还原至截取的位置且将位置相邻的裂缝子序列合并为一个裂缝。Step 4.4: The Softmax classifier receives the features output by the stacked autoencoder as its input, and outputs the category 0 or 1 of the strain subsequence, where 0 means non-slit and 1 means crack; for the feature z (i) output by the stacked autoencoder, Select the category l with the largest probability p(y (i) =l|z (i) ; δ) as the category corresponding to the feature. The distributed strain subsequence is restored to the truncated position and the adjacent fracture subsequences are merged into one fracture.
与现有技术相比,本发明具有以下有益的技术效果:Compared with the prior art, the present invention has the following beneficial technical effects:
本发明实现以分布式光纤传感器采集数据,改变以往点对点传感的方式。通过标准化缩小了不同数据之间的差异。同时,以基于堆叠自编码器的方法克服了分布式光纤传感器高空间分辨率与低信噪比的矛盾。堆叠自编码器能够在低信噪比的数据中提取高鲁棒性、可辨识的特征用以分类。在裂缝检测中效果显著,能够检测微小裂缝,在微小裂缝的检测效果上得提升。The invention realizes data collection by distributed optical fiber sensors, and changes the previous point-to-point sensing method. Standardization narrows the differences between different data. At the same time, the method based on stacked autoencoder overcomes the contradiction between high spatial resolution and low signal-to-noise ratio of distributed optical fiber sensor. Stacked autoencoders can extract highly robust and identifiable features for classification in low signal-to-noise ratio data. It has a significant effect in crack detection, can detect tiny cracks, and has improved the detection effect of tiny cracks.
附图说明Description of drawings
图1是本发明系统的流程示意图;Fig. 1 is the schematic flow chart of the system of the present invention;
图2是本发明中的自动编码器的示意图;Fig. 2 is the schematic diagram of the autoencoder in the present invention;
图3是本发明中的堆叠自编码器的示意图;Fig. 3 is the schematic diagram of the stacked autoencoder in the present invention;
图4是本发明方法的过程示意图;Fig. 4 is the process schematic diagram of the inventive method;
图5是本发明中具体的预训练与微调示意图。FIG. 5 is a schematic diagram of specific pre-training and fine-tuning in the present invention.
具体实施方式Detailed ways
下面结合附图对本发明作进一步详细描述:Below in conjunction with accompanying drawing, the present invention is described in further detail:
参见图1至图5,一种基于堆叠自编码器的分布式应变裂缝检测系统,包括应变序列采集模块;应变序列预处理模块;基于堆叠自编码器的特征自学习与表征模块;Softmax分类识别模块(具体流程如图1所示)。Referring to Figures 1 to 5, a distributed strain crack detection system based on stacked autoencoders includes a strain sequence acquisition module; a strain sequence preprocessing module; a feature self-learning and characterization module based on stacked autoencoders; Softmax classification and recognition module (the specific process is shown in Figure 1).
应变序列采集模块,用于采集结构体的分布式应变,采集到的结构体分布式应变为一个一维序列;The strain sequence acquisition module is used to collect the distributed strain of the structure, and the collected distributed strain of the structure is a one-dimensional sequence;
应变序列预处理模块包括:z-score标准化模块和滑窗模块,z-score标准化模块将应变序列标准化为0均值1标准差的数据。滑窗模块通过长度为21,步长为1的滑动窗口将标准化后的应变序列截取了一组长度都是21的应变子序列。将所得应变子序列按照截取位置标记标签,将以裂缝处为中心截取的应变子序列标记为裂缝子序列,并将其左侧3个和右侧4个应变子序列标记为裂缝子序列,其余标记为非裂缝子序列。The strain sequence preprocessing module includes: z-score normalization module and sliding window module. The z-score normalization module normalizes the strain sequence to data with 0 mean and 1 standard deviation. The sliding window module intercepts a set of strain subsequences with a length of 21 from the normalized strain sequence through a sliding window with a length of 21 and a step size of 1. Label the obtained strain subsequence according to the interception position, mark the strain subsequence intercepted at the crack as the crack subsequence, and mark the left 3 strain subsequences and the right four strain subsequences as the crack subsequence, and the rest marked as non-cracked subsequences.
基于堆叠自编码器的特征自学习与表征模块包括:3个自动编码器模块如附图2,自动编码器模块对于输入数据x,特征h与输出之间的关系可以表示为fθ(·)和gθ'(·)两个函数,具体如下所示:The feature self-learning and characterization module based on stacked autoencoder includes: 3 autoencoder modules as shown in Figure 2, the autoencoder module for input data x, feature h and output The relationship between can be expressed as two functions f θ ( ) and g θ' ( ), as follows:
h=fθ(x)=sf(Wx+bh)h=f θ (x)=s f (Wx+b h )
其中,W和bh分别为输入数据x和特征h之间的连接矩阵和偏置向量,WT和bv分别为特征h和输出之间的连接矩阵和偏置向量,WT为W的转置矩阵。sf和sg是激活函数。Among them, W and b h are the connection matrix and bias vector between the input data x and feature h, respectively, W T and b v are the feature h and output, respectively Between the connection matrix and the bias vector, W T is the transpose matrix of W. s f and s g are activation functions.
基于堆叠自编码器的特征自学习与表征模块的过程具体如下所示:The process of feature self-learning and representation module based on stacked autoencoder is as follows:
hi=sf(Wixi+bh,i)h i =s f (W i x i +b h,i )
其中,xi,hi,分别为第i(i=1,2,…,n)个自动编码器模块的输入数据x,特征h与输出Wi和Wi T为其连接权重矩阵,bh,i和bv,i为其偏置向量,xi=hi-1,sf和sg为激活函数。Among them, x i , h i , are the input data x, feature h and output of the i-th (i=1,2,...,n) autoencoder module, respectively Wi and Wi T are their connection weight matrices, b h,i and b v,i are their bias vectors, xi =h i-1 , s f and s g are activation functions.
采用Softmax分类器实现对应变子序列的分类,具体方法为:The classification of strain subsequences is realized by Softmax classifier, and the specific method is as follows:
构建Softmax分类器,对于给定的输入z,用假设函数hδ(z)针对每一个类别l估算出概率值p(y=l|z),l∈{0,1},假设函数hδ(z)输出一个t维的向量表示这t个估计的概率值,t=2,假设函数hδ(z)如下:Construct a Softmax classifier. For a given input z, use the hypothesis function h δ (z) to estimate the probability value p(y=l|z) for each category l, l∈{0,1}, assuming the function h δ (z) Output a t-dimensional vector representing the t estimated probability values, t=2, assuming that the function h δ (z) is as follows:
其中,δ1,δ2是Softmax分类器的全部参数,z(i)为输入,y(i)为输出,Softmax分类器将z分为类别l的概率为:Among them, δ 1 , δ 2 are all parameters of the Softmax classifier, With z (i) as input and y (i) as output, the probability that the Softmax classifier classifies z into class l is:
其中,z(i)为输入,y(i)为输出;where z (i) is the input and y (i) is the output;
Softmax分类器接收堆叠自编码器输出的特征作为其输入,输出应变子序列的类别0或1,0表示非裂缝,1表示裂缝;对于堆叠自编码器输出的特征z(i),选择概率p(y(i)=l|z(i);δ)最大的类别l作为该特征对应的类别。将分布式应变子序列还原至截取的位置且将位置相邻的裂缝子序列合并为一个裂缝。The Softmax classifier receives the features output by the stacked self-encoder as its input, and outputs the category 0 or 1 of the strain subsequence, where 0 means no crack and 1 means crack; for the feature z (i) output by the stacked self-encoder, the probability p is chosen (y (i) =l|z (i) ; δ) The largest category l is taken as the category corresponding to this feature. The distributed strain subsequence is restored to the truncated position and the adjacent fracture subsequences are merged into one fracture.
一种基于堆叠自编码器的分布式应变裂缝检测方法,具体步骤如图5所示:A distributed strain crack detection method based on stacked autoencoders, the specific steps are shown in Figure 5:
1)、应变序列采集;1), strain sequence collection;
2)、使用z-score标准化将采集到的应变序列进行标准化,使用长度为21,步长为1的滑动窗口截取应变序列,得到应变子序列,将应变子序列按照截取位置进行标记;2), use z-score standardization to standardize the acquired strain sequence, use a sliding window with a length of 21 and a step size of 1 to intercept the strain sequence, obtain a strain subsequence, and mark the strain subsequence according to the interception position;
2.1:将采集到的应变序列减去其均值,除以其方差,得到0均值1标准差的数据。2.1: Subtract the mean value of the acquired strain series and divide it by its variance to obtain data with 0 mean and 1 standard deviation.
2.2:使用长度为21,步长为1的滑动窗口沿着采集到的应变序列进行滑动,将应变序列截取为一组长度为21的应变子序列。2.2: Use a sliding window with a length of 21 and a step size of 1 to slide along the acquired strain sequence, and cut the strain sequence into a set of strain subsequences with a length of 21.
2.3:将所得应变子序列按照截取位置标记标签,将以裂缝处为中心截取的应变子序列标记为裂缝子序列,并将其左侧3个和右侧4个应变子序列标记为裂缝子序列。2.3: Label the obtained strain subsequence according to the interception position, mark the strain subsequence intercepted at the crack as the crack subsequence, and mark the 3 left and 4 right strain subsequences as the crack subsequence .
3)、用基于堆叠自编码器的神经网络自动学习表征应变子序列的特征;3), use the neural network based on the stacked autoencoder to automatically learn the characteristics of the strain subsequence;
3.1:模型初始化,确定模型的层数与神经元个数。随机初始化模型中的连接权重矩阵和偏置向量。输入层的神经元个数等于21,为应变子序列的长度。3.1: Model initialization, determine the number of layers and neurons of the model. Randomly initialize the connection weight matrix and bias vector in the model. The number of neurons in the input layer is equal to 21, which is the length of the strain subsequence.
3.2:预训练堆叠自编码器,堆叠自编码器由3个自动编码器构成,用得到的应变子序列预训练每个自动编码器。预训练自动编码器的损失函数为输入与输出之间的均方误差,具体如下:3.2: Pre-training stacked autoencoders. The stacked autoencoders consist of 3 autoencoders, and each autoencoder is pretrained with the resulting strain subsequences. The loss function of the pretrained autoencoder is the mean squared error between the input and the output, as follows:
其中,x为输入的应变子序列,为自动编码器输出的重构数据,M为所有输入的应变子序列的数量,Xm、分别是输入模型的第m条应变子序列和对应的输出重构的第m条子序列。where x is the input strain subsequence, is the reconstructed data output by the autoencoder, M is the number of all input strain subsequences, X m , are the mth strain subsequence of the input model and the mth subsequence of the corresponding output reconstruction, respectively.
4)、采用Softmax分类器实现对提取到的应变子序列的特征进行二分类,完成裂缝检测;4), adopt the Softmax classifier to realize the two-classification of the features of the extracted strain subsequence, and complete the crack detection;
4.1:构建Softmax分类器,对于给定的输入z,用假设函数hδ(z)针对每一个类别l估算出概率值p(y=l|z),l∈{0,1},假设函数hδ(z)输出一个t维的向量表示这t个估计的概率值,t=2,假设函数hδ(z)如下:4.1: Build a Softmax classifier, for a given input z, use the hypothesis function h δ (z) to estimate the probability value p(y=l|z) for each category l, l∈{0,1}, the hypothesis function h δ (z) outputs a t-dimensional vector representing the t estimated probability values, t=2, assuming that the function h δ (z) is as follows:
其中,δ1,δ2是Softmax分类器的全部参数,z(i)为输入,y(i)为输出,Softmax分类器将z分为类别l的概率为:Among them, δ 1 , δ 2 are all parameters of the Softmax classifier, With z (i) as input and y (i) as output, the probability that the Softmax classifier classifies z into class l is:
其中,z(i)为输入,y(i)为输出;where z (i) is the input and y (i) is the output;
4.2:预训练Sofmax,将应变子序列输入预训练后的堆叠自编码器,得到输出的特征z(i),以z(i)及其标签类别y(i)预训练Softmax,损失函数为交叉熵函数,具体如下:4.2: Pre-training Sofmax, input the strain subsequence into the pre-trained stacked autoencoder, get the output feature z (i) , pre-train Softmax with z (i) and its label category y (i) , and the loss function is crossover The entropy function, as follows:
其中,为Softmax分类器的全部参数,为输出的类别概率,λ1为Softmax中连接权重矩阵和偏置向量正则项的权重系数,M为输入应变子序列的总个数,K为类别数,为2。in, are all parameters of the Softmax classifier, is the output category probability, λ 1 is the weight coefficient of the regular term connecting the weight matrix and the bias vector in Softmax, M is the total number of input strain subsequences, and K is the number of categories, which is 2.
4.3:微调,堆叠自编码器的编码部分后接Softmax分类器,可以使其具有分类功能。利用预训练得到的应变子序列微调堆叠自编码器的编码部分与Softmax分类器整体结构的连接权重矩阵和偏置向量。微调时的损失函数为交叉损失函数,具体如下:4.3: Fine-tuning, the encoding part of the stacked autoencoder is followed by a Softmax classifier, which can make it have a classification function. Using the pretrained strain subsequences to fine-tune the connection weight matrix and bias vector between the encoding part of the stacked autoencoder and the overall structure of the Softmax classifier. The loss function during fine-tuning is a cross loss function, as follows:
其中,ω为堆叠自编码器中的连接权重矩阵和偏置向量,Θ为ω和δ,λ2为堆叠自编码器中连接权重矩阵和偏置向量正则项的权重系数。where ω is the connection weight matrix and bias vector in the stacked autoencoder , Θ is ω and δ, and λ2 is the weight coefficient of the regular term of the connection weight matrix and bias vector in the stacked autoencoder.
4.4:Softmax分类器接收堆叠自编码器输出的特征作为其输入,输出应变子序列的类别0或1,0表示非裂缝,1表示裂缝;对于堆叠自编码器输出的特征z(i),选择概率p(y(i)=l|z(i);δ)最大的类别l作为该特征对应的类别。将分布式应变子序列还原至截取的位置且将位置相邻的裂缝子序列合并为一个裂缝。4.4: The Softmax classifier receives the features output by the stacked self-encoder as its input, and outputs the category 0 or 1 of the strain subsequence, 0 means non-crack and 1 means crack; for the feature z (i) output by the stacked self-encoder, choose The category l with the largest probability p(y (i) =l|z (i) ; δ) is used as the category corresponding to the feature. The distributed strain subsequence is restored to the truncated position and the adjacent fracture subsequences are merged into one fracture.
实施效果Implementation Effect
首先将光纤传感器预张紧,然后通过环氧树脂粘附于钢结构体表面。光纤传感器的两端接于基于BOTDA的分布式光纤传感系统,得到结构体表面沿光纤传感器径向分布的分布式应变数据。采用本发明基于堆叠自编码器的微小裂缝检测方法,基于采集的分布式应变数据可以准确无遗漏地检出张口宽度为32μm的微小裂缝,是一种可以用于钢结构表面微小裂缝分布式检测的有效方法。The fiber optic sensor is first pre-tensioned and then adhered to the surface of the steel structure by epoxy. The two ends of the optical fiber sensor are connected to the distributed optical fiber sensing system based on BOTDA, and the distributed strain data distributed along the radial direction of the optical fiber sensor on the surface of the structure body are obtained. By adopting the micro-crack detection method based on the stacked self-encoder of the present invention, the micro-cracks with an opening width of 32 μm can be accurately detected based on the collected distributed strain data, which is a distributed detection method for micro-cracks on the surface of steel structures. effective method.
Claims (5)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910973652.6A CN110715929B (en) | 2019-10-14 | 2019-10-14 | Distributed strain micro crack detection system and method based on stacking self-encoder |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910973652.6A CN110715929B (en) | 2019-10-14 | 2019-10-14 | Distributed strain micro crack detection system and method based on stacking self-encoder |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110715929A CN110715929A (en) | 2020-01-21 |
CN110715929B true CN110715929B (en) | 2022-03-25 |
Family
ID=69211557
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910973652.6A Expired - Fee Related CN110715929B (en) | 2019-10-14 | 2019-10-14 | Distributed strain micro crack detection system and method based on stacking self-encoder |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110715929B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111750795A (en) * | 2020-06-18 | 2020-10-09 | 哈尔滨工程大学 | A distributed creep measurement system and measurement method |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102322808B (en) * | 2011-08-09 | 2013-03-27 | 中国计量学院 | Very long range pulse coding distribution type Fiber Raman and Brillouin photon sensor |
CN202195825U (en) * | 2011-08-09 | 2012-04-18 | 中国计量学院 | Extra-long distance pulse-coding distributed optical fiber Raman and Brillouin photon sensor |
US10846447B2 (en) * | 2016-04-29 | 2020-11-24 | Exxonmobil Upstream Research Company | Method and system for stacking fracture prediction |
-
2019
- 2019-10-14 CN CN201910973652.6A patent/CN110715929B/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN110715929A (en) | 2020-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109765053B (en) | Rolling bearing fault diagnosis method using convolutional neural network and kurtosis index | |
Iyer et al. | Ultrasonic signal processing methods for detection of defects in concrete pipes | |
CN107144643B (en) | A kind of damnification recognition method of Lamb wave monitoring signals statistical parameter | |
CN102928435A (en) | Aircraft skin damage identification method and device based on image and ultrasound information fusion | |
CN114782753A (en) | Lung cancer histopathological whole section classification method based on weakly supervised learning and converter | |
CN119023677A (en) | A highly accurate intelligent detection system for wiring harness defects | |
CN108537790A (en) | Heterologous image change detection method based on coupling translation network | |
CN110715929B (en) | Distributed strain micro crack detection system and method based on stacking self-encoder | |
CN110738168A (en) | A distributed strain micro-crack detection system and method based on stacked convolutional autoencoders | |
Chen et al. | A mixed samples-driven methodology based on denoising diffusion probabilistic model for identifying damage in carbon fiber composite structures | |
CN114358189A (en) | Fault diagnosis method of hydraulic waterproof valve based on multimodal deep residual shrinkage network | |
CN119086721B (en) | Plastic product strength detection method and device based on non-contact measurement | |
TW202144763A (en) | Tomography method, system and apparatus based on time-domain spectroscopy | |
Banerjee et al. | Damage detection and localization by learning deep features of elastic waves in piezoelectric ceramic using point contact method | |
CN118447000A (en) | Welding defect detection method and system | |
Masaki et al. | Surface undulation detection system using wearable artificial skin layer with strain gauge | |
CN112446612B (en) | Assessment method of damage assessment system of soft rigid arm mooring system connection structure | |
CN116818574A (en) | Material mechanical properties measurement method based on terahertz spectrum characteristic parameter characterization | |
Liu et al. | Automatic terahertz recognition of hidden defects in layered polymer composites based on a deep residual network with transfer learning | |
Lin et al. | Automatic elimination of invalid impact-echo signals for detecting delamination in concrete bridge decks based on deep learning | |
CN119150120A (en) | Impact damage identification method based on multi-mode data fusion and ensemble learning | |
CN114104328A (en) | Aircraft state monitoring method based on deep migration learning | |
CN106679886A (en) | Nonlinear fault detecting and identifying method of self-confirming air data system | |
CN118777562B (en) | A steel quality analysis and feedback system | |
CN117274691B (en) | Hyperspectral image classification method based on space pooling transducer |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220325 |
|
CF01 | Termination of patent right due to non-payment of annual fee |