CN114923420A - Crack diagnosis method and system based on fiber Bragg grating and storage medium - Google Patents
Crack diagnosis method and system based on fiber Bragg grating and storage medium Download PDFInfo
- Publication number
- CN114923420A CN114923420A CN202210544130.6A CN202210544130A CN114923420A CN 114923420 A CN114923420 A CN 114923420A CN 202210544130 A CN202210544130 A CN 202210544130A CN 114923420 A CN114923420 A CN 114923420A
- Authority
- CN
- China
- Prior art keywords
- crack
- fiber bragg
- length
- bragg grating
- reflection spectrum
- 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
Links
Images
Classifications
-
- 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/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
-
- 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
- G01B11/165—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge by means of a grating deformed by the object
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
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
The invention relates to the technical field of train structure crack monitoring, in particular to a crack diagnosis method, a crack diagnosis system and a storage medium based on fiber Bragg gratings, wherein the method comprises the steps of determining material parameters of a manufacturing material of a structure to be analyzed; establishing a finite element simulation model, simulating the crack propagation condition based on the finite element simulation model, and acquiring structural strain data of the crack propagation detected by the fiber Bragg grating sensor to different lengths; reconstructing the structural strain data into a reflectance spectrum; extracting a damage sensitive characteristic value and a reference signal in the reflection spectrum, wherein the damage sensitive characteristic value is used for indicating the length of the crack, and the reference signal is used for indicating that no crack is expanded; and taking the damage sensitivity characteristic value as an input, taking the crack length as an output to construct a crack length regression model, and diagnosing the crack based on the crack length regression model. The method can solve the problems that the traditional finite element method is complex in calculation for acquiring the strain output, complicated in steps and incomplete in quantitative diagnosis and monitoring of the cracks based on a single characteristic value.
Description
Technical Field
The invention relates to the technical field of train structure crack monitoring, in particular to a crack diagnosis method and system based on fiber Bragg gratings and a storage medium.
Background
A large number of porous structures contained in rail vehicles, mechanical equipment and the like are easy to crack under the action of cyclic load, and the crack is allowed to expand to cause structural function failure and serious safety accidents. The existing crack quantitative diagnosis method based on the fiber Bragg grating reflection spectrum generally utilizes a single characteristic value represented by reflection spectrum central wavelength and broadening to establish a crack quantitative diagnosis model. As the reflection spectrum of the fiber Bragg grating depends on the accuracy of a strain field to a great extent, the traditional finite element method has complex calculation and complex steps for acquiring strain output, and the quantitative diagnosis and monitoring of the cracks based on a single characteristic value are not comprehensive enough.
Disclosure of Invention
The invention provides a crack diagnosis method, a crack diagnosis system and a storage medium based on fiber Bragg gratings, and aims to solve the problems that the traditional finite element method is complex in calculation for obtaining strain output, complicated in steps and incomplete in quantitative crack diagnosis and monitoring based on a single characteristic value.
In order to achieve the purpose, the invention is realized by the following technical scheme:
in a first aspect, the present invention provides a crack diagnosis method based on fiber bragg grating, including:
determining material parameters of a manufacturing material of a structure to be analyzed;
establishing a finite element simulation model according to the material parameters, simulating the crack propagation condition based on the finite element simulation model, and acquiring structural strain data of the crack propagation detected by the fiber Bragg grating sensor under different lengths;
reconstructing the structural strain data into a reflection spectrum by adopting a transmission matrix method;
extracting a damage sensitivity characteristic value and a reference signal in the reflection spectrum, wherein the damage sensitivity characteristic value is used for representing the length of the crack, and the reference signal is used for representing no expansion crack;
and taking the damage sensitivity characteristic value as an input, taking the crack length as an output to construct a crack length regression model, and diagnosing the crack based on the crack length regression model.
In a second aspect, the present application provides a fiber bragg grating based crack diagnosis system, including M fiber bragg grating sensors disposed at a target structure of a train, where M is a positive integer, and a processing center connected to the M sensors, where the processing center is configured to:
determining material parameters of a manufacturing material of a structure to be analyzed;
establishing a finite element simulation model according to the material parameters, simulating the crack propagation condition based on the finite element simulation model, and acquiring structural strain data of the crack propagation detected by the fiber Bragg grating sensor under different lengths;
reconstructing the structural strain data into a reflection spectrum by adopting a transmission matrix method;
extracting damage sensitive characteristic values and reference signals in the reflection spectrum, wherein the damage sensitive characteristic values are used for representing the length of the crack, and the reference signals are used for representing no expansion crack;
and taking the damage sensitivity characteristic value as an input, taking the crack length as an output to construct a crack length regression model, and diagnosing the crack based on the crack length regression model.
In a third aspect, the present application provides a fiber bragg grating based crack diagnosis system, which includes a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor implements the steps of the method according to the first aspect when executing the computer program.
In a fourth aspect, the present application provides a computer readable storage medium having stored thereon a computer program which, when executed by a processor, performs the method steps as set forth in the first aspect.
Has the advantages that:
the crack diagnosis method based on the fiber Bragg grating provided by the invention simulates the crack expansion process under the condition of cyclic loading by using an expansion finite element method, reconstructs the fiber Bragg grating reflection spectrum by using a transmission matrix method, further determines the action mechanism of crack expansion change and the fiber Bragg grating reflection spectrum during crack expansion, extracts a plurality of damage sensitive characteristic values of the reflection spectrum, and constructs a crack length regression model of the plurality of damage sensitive characteristic values and the crack length by using a support vector regression method, thus, the accuracy of the constructed model can be higher, the diagnosis result of the model is close to the actual crack length, and overcomes the problem that the traditional finite element method has complex calculation in the crack propagation simulation, in addition, the method can comprehensively monitor the structure to be analyzed based on a plurality of damage sensitivity characteristic values.
According to the crack diagnosis system based on the fiber Bragg grating, provided by the invention, the fiber Bragg grating sensor arranged at the key position of the structure can accurately sense the structural strain change caused by the crack, and a crack quantitative monitoring model can be established by extracting the characteristic value capable of representing the crack length in the reflection spectrum of the fiber Bragg grating, so that the real-time monitoring on the crack damage is realized.
Drawings
Fig. 1 is a flowchart of a crack diagnosis method based on fiber bragg grating according to a preferred embodiment of the present invention;
FIG. 2 is a schematic diagram of the dimensions of a simulation test piece according to the preferred embodiment of the present invention;
FIG. 3 is a schematic reflection spectrum of the FBG3 sensor location in accordance with the preferred embodiment of the present invention;
FIG. 4 is a schematic diagram illustrating the variation of damage characteristics in simulation according to the preferred embodiment of the present invention;
FIG. 5 is a schematic diagram of the diagnostic results of the preferred embodiment of the present invention.
Detailed Description
The technical solutions of the present invention are described clearly and completely below, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without inventive step based on the embodiments of the present invention, are within the scope of protection of the present invention.
Unless defined otherwise, technical or scientific terms used herein shall have the ordinary meaning as understood by one of ordinary skill in the art to which this invention belongs. The use of "first," "second," and similar terms in the present application do not denote any order, quantity, or importance, but rather the terms are used to distinguish one element from another. Also, the use of the terms a, an, etc. do not denote a limitation of quantity, but rather denote the presence of at least one. The terms "connected" or "coupled" and the like are not restricted to physical or mechanical connections, but may include electrical connections, whether direct or indirect. "upper", "lower", "left", "right", and the like are used only to indicate relative positional relationships, and when the absolute position of the object to be described is changed, the relative positional relationships are changed accordingly.
It should be noted that, because the operating environment of the high-speed rail is complex, the situation of foreign object impact often occurs in each structure, and such an impact event easily causes structural damage and threatens driving safety, so it is very necessary to monitor the impact event in each structure of the high-speed rail. Because the reflection spectrum of the fiber Bragg grating depends on the accuracy of a strain field to a great extent, the traditional finite element method has complex calculation and complex steps for acquiring strain output, and the quantitative diagnosis and monitoring of cracks based on a single characteristic value are not comprehensive enough. Based on the method, the application provides a crack diagnosis method based on the fiber Bragg grating.
Referring to fig. 1, an embodiment of the present application provides a crack diagnosis method based on a fiber bragg grating, including:
determining material parameters of a manufacturing material of a structure to be analyzed;
establishing a finite element simulation model according to the material parameters of the damage sensitivity characteristic value, simulating the crack propagation condition based on the finite element simulation model of the damage sensitivity characteristic value, and acquiring structural strain data of the crack propagation detected by the fiber Bragg grating sensor under different lengths;
reconstructing the damage sensitive characteristic value structure strain data into a reflection spectrum by adopting a transmission matrix method;
extracting a damage sensitive characteristic value and a reference signal in a damage sensitive characteristic value reflection spectrum, wherein the damage sensitive characteristic value of the damage sensitive characteristic value is used for indicating the length of the crack, and the reference signal of the damage sensitive characteristic value is used for indicating that no crack is expanded;
and taking the damage sensitive characteristic value as an input, taking the crack length as an output to construct a crack length regression model, and diagnosing the crack based on the crack length regression model of the damage sensitive characteristic value.
In this embodiment, a Fiber Bragg Grating (FBG) sensor is used as the detection sensor. That is, an optical fiber sensor using a bragg grating as a sensing element can directly measure temperature and strain and can indirectly measure physical quantities related to temperature and strain. Thus, the sensitivity of the fiber bragg grating reflectance spectrum to the effects of strain fields is extremely high. Therefore, in the embodiment, accurate and reliable extension monitoring is carried out by analyzing the rule between the change of the reflection spectrum and the crack length, so that the failure rate of mechanical equipment can be reduced, and property loss can be reduced.
The crack diagnosis method based on the fiber Bragg grating simulates the crack expansion process under the condition of cyclic loading by using the finite element expansion method, reconstructs the fiber Bragg grating reflection spectrum by using the transmission matrix method, further determines the action mechanism of the crack expansion change and the fiber Bragg grating reflection spectrum during the crack expansion, extracts a plurality of damage sensitive characteristic values of the reflection spectrum, and constructs a crack length regression model of the plurality of damage sensitive characteristic values and the crack length by using the support vector regression method, thus, the accuracy of the constructed model can be higher, the diagnosis result of the model is close to the actual crack length, and overcomes the problem of complex calculation in crack propagation simulation by the traditional finite element method, in addition, the method can comprehensively monitor the structure to be analyzed based on a plurality of damage sensitivity characteristic values.
The reconstructing the structural strain data into a reflection spectrum by adopting a transmission matrix method comprises the following steps:
assuming that the length of the uniform grating is L, when calculating the fiber Bragg grating reflection spectrum (FBG reflection spectrum) under the action of non-uniform strain, the non-uniform grating is uniformly divided into N small segments, wherein N is a positive integer, and simultaneously the requirement of calculating the FBG reflection spectrum under the action of non-uniform strain is metλ B For the center wavelength, the expression is:
λ B =2n eff Λ;
wherein Λ is the period of change of the refractive index of the grating, n eff Is the effective refractive index.
Regarding the average period of each section in the N small sections as the equivalent period of the section, and regarding the refractive index of each section as the equivalent refractive index of the section;
substituting the parameters of each section into a coupling equation to carry out iterative calculation so as to obtain the reflection spectrum of the whole FBG, wherein the calculation process is as follows:
Λ i =Λ 0 (1+aε zz );
in the formula, Λ 0 Initial grating period, ε zz Is the axial average strain of the ith segment, a is the grating strain coefficient, and the expression is as follows:
in the formula, n eff0 Is the average effective refractive index of FBG in free state, v is the modulation depth, p 11 And p 12 Respectively setting an optical transfer matrix of each grating segment to generate a 2 x 2T matrix T for effective optical stress tensor components based on modal coupling theory i The following:
in the formula, R i And S i Is the amplitude of the forward transmission mode and the amplitude of the backward transmission mode of the i-th section, T i The expression satisfies the following relationship:
wherein, Delta z is the length of the grating segment,and κ is the direct current self-coupling coefficient and the alternating current self-coupling coefficient of the i-th segment, respectively, and the expressions are:
in the formula, δ n eff Defining gamma as the mean of the refractive index variation over the grating periodObtaining T matrix T ═ T of all grating segments through iterative computation N T N-1 T N-2 …T 2 T 1 Then, the following are obtained:
in the formula, R 0 、S 0 Is the amplitude of the forward transmission mode and the amplitude of the backward transmission mode of the first section, R L 、S L The amplitude of the forward transmission mode and the amplitude of the backward transmission mode of the tail section;
the reflectivity r of the FBG for each wavelength is calculated as follows:
where ρ is the reflection coefficient of FBG, T 21 And T 11 The intermediate calculation results for calculating the reflectivity r are shown as the 1 st number in the 2 nd row and the 1 st number in the 1 st row in the T matrix, respectively.
And obtaining the FBG reflection spectrum in the whole wavelength interval through iterative calculation, and calculating the corresponding FBG reflection spectrum under different crack lengths to obtain the reflection spectrum of the whole crack propagation process.
In the embodiment, simulation can be subdivided into two steps, firstly, crack propagation is simulated through abaqus software, and strain data at the position along the sensor is obtained; and then, converting the strain into a reflection spectrum by using a transmission matrix method, further analyzing the action mechanism of the strain on the reflection spectrum, extracting characteristics and the like. Therefore, reflection spectrum data under multiple groups of crack length samples can be obtained, and the method is low in cost, simple, efficient and feasible compared with the test.
Wherein the extracting of the damage sensitive characteristic value and the reference signal in the reflection spectrum comprises:
the change of the reflection spectrum along with crack propagation is represented by wavelength shift, the wavelength shift can be obtained by subtracting the wavelength of a reference spectrum from the central wavelength under the damage of cracks, and the central wavelength lambda of the reflection spectrum is extracted c Satisfies the following relation:
wherein λ represents a reflectance spectrum wavelength, and R (λ) represents a reflectance;
the variation of the different crack length spreads during crack propagation is characterized by the full width of the reflection line at a reflectance value L, wherein the spread b satisfies the following relationship:
b=|λL1-λ Lend |;
in the formula, λ L1 Is the wavelength corresponding to the first limit of reflectivity L, and λ Lend The wavelength corresponding to the last reflectivity limit value L;
the wave crest number is used for representing the change of the wave crest numbers of different crack lengths in the crack propagation process;
and representing the area change of the reflection spectrum of different crack lengths in the crack propagation process by using the area of the shape surrounded by the reflection spectrum line and the horizontal axis of the coordinate, wherein the area S expression of the normalized reflection spectrum meets the following relational expression:
in the formula, λ 1 Is the first wavelength, lambda, of the wavelength range of the reflectance spectrum en d is the last wavelength of the wavelength range of the reflectance spectrum,represents a reflection line;
representing the change of the overlapping area of different crack lengths in the crack propagation process by the overlapping part of the crack-free reference reflection spectral line and each reflection spectrum in the crack propagation process, wherein the overlapping area S is normalized c The expression satisfies the following relation:
in the formula of lambda s Is the first wavelength, λ, of the coincidence range of the wavelengths of the reflection lines e Is the last wavelength of the coincidence range of the wavelengths of the reflection lines,representing coincident reflection lines.
Representing the change of correlation coefficients of different crack lengths in the crack propagation process by using the crack-free reference reflection spectral line and the damage spectrum in the crack propagation process, wherein the correlation coefficient C m Satisfies the following relational expression:
In the formula, ρ 0 As a reference spectral reflectance vector, ρ m Is the damage spectral reflectance vector, N is the length of the reflectance vector,for any wavelength shift from the damage spectrum, k is the number of reflectivities in the reflectivity vector, k is 1,2 R ;
Characterizing the significance degree of reflection spectrum chirp phenomena of different crack lengths in a crack propagation process by a fractal dimension, wherein the fractal dimension FD is calculated by a box counting method, and the expression of the FD satisfies the following relational expression:
wherein r is the length of the grid side dividing the image, and r is 2 i And i is 0,1,2, 10, m is the number of grids, a line graph of log (1/r) -log (m) is drawn finally, and the fractal dimension FD is calculated by taking two points in the log graph to calculate the slope.
It should be noted that the fractal dimension serves as a measure of the irregularity of complex geometric figures. In the embodiment, the box counting method is selected to calculate the fractal dimension of the reflection spectrum, which mainly reflects the significance degree of the chirp phenomenon of the reflection spectrum, wherein the box counting method is meaningful and intuitive, and the calculation is simple and efficient.
When the box counting method is adopted for calculation, firstly, the proportion of a horizontal axis and a vertical axis of coordinates is set as 1, and a square reflection spectrum is drawn; then removing the coordinate axis and the frame of the image to obtain a pure reflection spectrum graph; dividing the whole square image by a small grid with the side length of r, and calculating the number m of grids in which the ratio of the included parts to the non-included parts of the reflection spectrum curve in the grids is larger than a certain value; finally, a line plot of log (1/r) -log (m) is plotted. The fractal dimension is calculated by taking two points in a logarithmic graph to calculate the slope, and consideringUntil the grid size is smaller, the grid division is finer, the fractal description of the reflection spectrum is more accurate, and therefore, r is 2 in the embodiment 0 And r is 2 1 The slope is calculated at two points, so that the fractal dimension of the reflection spectrum is calculated by adopting a box counting method, and the method is simple and direct and has high calculation efficiency.
The method for constructing the crack length regression model by taking the damage sensitivity characteristic value as an input and the crack length as an output comprises the following steps:
acquiring crack length sample data set { (x) i , y i 1,2, …, n being the number of samples of different crack lengths, x i ,y i Respectively an input quantity and an output quantity, where x i ={x i1 ,x i2 ,…,x ip To influence crack length y i Damage sensitive eigenvalue of, y i Is the true crack length for the ith sample, p is the effect y i Wherein the estimation function of the crack length satisfies the following relation:
f(x)=ωφ(x)+b;
in the formula, omega is weight, phi (x) is high-dimensional nonlinear function, and b is offset;
solving ω and b requires minimizing the optimization model and introducing the relaxation factor ξ as follows:
y i -ωx i -b≤ε+ξ i ;ωx i -y i +b≤ε+ξ i * ;
wherein C is a penalty factor xi i 、ξ i * Introducing alpha for relaxation variables, epsilon for loss function and N for number of training samples i ,The crack length regression model obtained by establishing a Lagrange function by the Lagrange operator to solve the dual problem meets the following relational expression:
in the formula, alpha i 、Lagrangian for the ith sample, K (x) i ,x j ) Is a radial basis kernel function, x i ,x j Represents any two samples in the training samples, the training samples (x) for a given N crack lengths i I ═ 1, 2., N), each training sample corresponds to a feature vector composed of seven features (wavelength shift, broadening, reflection spectrum area, 1-coincidence area, number of peaks, 1-correlation coefficient, fractal dimension), and any two crack length samples { x } in the training samples are processed i , x j 1,2, ·, N; j ═ 1,2,. N } (with x) j Representing another crack length sample) into the radial basis kernel function computation K (x) i ,x j )。
Wherein the method further comprises: and optimizing the penalty factor C and the variance g of the radial basis kernel function by using a cross-validation method to optimize the crack length regression model. Finally, with R 2 And (4) making an evaluation standard, and then applying a crack length regression model to diagnose the crack length.
In conclusion, the quantitative relation model between the characteristics and the crack length is established by extracting the characteristic values. Specifically, when a support vector regression method is used for establishing a quantitative relation between multiple features and crack length, a part of samples are selected from all sample data to serve as training samples, a multiple feature-crack length support vector regression model is trained, training accuracy depends on selection of a penalty factor C and a variance g of a radial basis kernel function, C and g can be selected iteratively one by using a cross-validation method to obtain optimal model parameters, and meanwhile, the fitting goodness R2 is used for evaluating the goodness and the badness of the training model. Finally, the remaining test sample can be used to diagnose the crack length. Thus, the diagnosis efficiency and the diagnosis accuracy can be improved.
The method steps of the present application are further explained below with a certain simulation example.
In this example, the strain at the crack propagation extraction sensor length location was simulated using the Abaqus finite element simulation software. The mechanical property parameters of the test piece material are shown in table 1, and the sensor parameters are shown in table 2.
TABLE 1 mechanical Property parameters of test piece materials
TABLE 2 sensor parameters
Referring to fig. 2, fig. 2 is a schematic diagram of the size and position of the test piece sensor, where in fig. 2, 1 indicates FBG 1; 2 denotes an FBG 2; FBG3 is denoted 3; 4 denotes an FBG 4; 5 denotes an FBG 5; and 6 denotes an FBG 6. Fatigue crack propagation was triggered by introducing a 3mm pre-crack through a 10mm diameter hole in the center of the slab. And applying a constant-amplitude load to the test piece, wherein the maximum loading amplitude is 80MPa, the minimum amplitude is 8MPa, the loading frequency is 5Hz, the boundary of the top of the test piece is fixed, and the bottom of the test piece is applied with an even tensile load of 80 MPa.
In the simulation, the grid formed by linear triangular elements is discretized in the whole test piece plate range, the grid is refined in the right center part of the test piece with expected crack propagation, and the initial crack size is 3 mm. The strain distribution along the FBG1-6 sensor for various crack sizes was calculated using higher order propagation finite elements, and then the strain was reconstructed as a reflection spectrum using a transmission matrix method.
Referring to fig. 3, fig. 3 shows a reflection spectrum of the FBG3 sensor position, wherein a total of 46 sets of crack sample data are extracted, and as the crack propagates, the reflection spectrum generates a position shift chirp phenomenon. Referring to fig. 4, fig. 4 shows the change of the extracted damage features, and it can be determined from fig. 4 that the rule of the change of the features along with the crack length is more obvious, and the features are different, so that the single features are combined to obtain a more comprehensive monitoring result.
And (3) taking the characteristic value of the sample reflection spectrum as the input of the neural network and the corresponding crack length as the output, establishing the relation between the multi-damage sensitive characteristic and the crack length, and realizing the monitoring of the crack length. For each FBG, 36 groups of samples are randomly selected from 46 groups of crack sample data to be used as training samples, each group of samples contains 7 characteristic values which are used as x i Input, x i The input quantity is 7 input quantities of the corresponding central wavelength, broadening, wave peak number, reflection spectrum area, coincidence area of damage spectrum and health spectrum, fractal dimension and correlation coefficient of each sensor, y i Inputting 1 output quantity corresponding to the crack length of the sample into an SVR neural network for training, and optimizing a penalty factor C and a kernel function g by using a cross-validation method. The remaining 10 groups were used as test samples, and the crack length diagnostic values and errors for all sensor fusions are shown in table 3, which is schematically shown in fig. 5:
TABLE 3 crack Length diagnostic values and errors
As can be determined from table 3 and fig. 5, the results of the diagnosis crack length and the actual crack length are very similar, and it can be seen that the diagnosis result of the present application is relatively accurate.
The application also provides a crack diagnostic system based on fiber bragg grating, including setting up in M fiber bragg grating sensors and the processing center of the target structure department of train, M is the positive integer, the processing center with M sensor connection, the processing center is used for:
determining material parameters of a manufacturing material of a structure to be analyzed;
establishing a finite element simulation model according to the material parameters, simulating the crack propagation condition based on the finite element simulation model, and acquiring structural strain data of the crack propagation detected by the fiber Bragg grating sensor under different lengths;
reconstructing the structural strain data into a reflection spectrum by adopting a transmission matrix method;
extracting damage sensitive characteristic values and reference signals in the reflection spectrum, wherein the damage sensitive characteristic values are used for representing the length of the crack, and the reference signals are used for representing no expansion crack;
and taking the damage sensitivity characteristic value as an input, taking the crack length as an output to construct a crack length regression model, and diagnosing the crack based on the crack length regression model.
According to the crack diagnosis system based on the fiber Bragg grating, the fiber Bragg grating sensors arranged at the key positions of the structure can accurately sense the structural strain change caused by the crack, and a crack quantitative monitoring model can be established by extracting the characteristic value capable of representing the crack length in the reflection spectrum of the fiber Bragg grating, so that the real-time monitoring on the crack damage is realized.
The present application further provides a fiber bragg grating based crack diagnosis system, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of the method when executing the computer program. The crack diagnosis system based on the fiber bragg grating can realize each embodiment of the crack diagnosis method based on the fiber bragg grating, and can achieve the same beneficial effects, and the details are not repeated here.
Embodiments of the present application also provide a computer-readable storage medium, on which a computer program is stored, which when executed by a processor implements the method steps as described above. The readable storage medium can implement the embodiments of the method described above, and can achieve the same beneficial effects, which are not described herein again.
The foregoing detailed description of the preferred embodiments of the invention has been presented. It should be understood that numerous modifications and variations could be devised by those skilled in the art in light of the present teachings without departing from the inventive concepts. Therefore, the technical solutions that can be obtained by a person skilled in the art through logical analysis, reasoning or limited experiments based on the prior art according to the concepts of the present invention should be within the scope of protection determined by the claims.
Claims (8)
1. A crack diagnosis method based on fiber Bragg grating is characterized by comprising the following steps:
determining material parameters of a manufacturing material of a structure to be analyzed;
establishing a finite element simulation model according to the material parameters, simulating the crack propagation condition based on the finite element simulation model, and acquiring structural strain data of the crack propagation detected by the fiber Bragg grating sensor under different lengths;
reconstructing the structural strain data into a reflection spectrum by adopting a transmission matrix method;
extracting damage sensitive characteristic values and reference signals in the reflection spectrum, wherein the damage sensitive characteristic values are used for representing the length of the crack, and the reference signals are used for representing no expansion crack;
and taking the damage sensitivity characteristic value as an input, taking the crack length as an output to construct a crack length regression model, and diagnosing the crack based on the crack length regression model.
2. The fiber bragg grating based crack diagnosis method according to claim 1, wherein the reconstructing the structural strain data into a reflection spectrum by using a transmission matrix method comprises:
setting the length of the uniform grating as L, and when calculating the reflection spectrum of the fiber Bragg grating under the action of non-uniform strain, uniformly dividing the non-uniform grating into N small sections, wherein N is a positive integer and simultaneously satisfies the requirementλ B For the center wavelength, the expression satisfies the following relation:
λ B =2n eff Λ;
wherein Λ is the period of change of the refractive index of the grating, n eff Is the effective refractive index;
regarding the average period of each segment in the N small segments as the equivalent period of the segment, and regarding the refractive index of each segment as the equivalent refractive index of the segment;
substituting the parameters of each section into a coupling equation to carry out iterative computation, wherein the process of computing the fiber Bragg grating reflection spectrum comprises the following steps:
Λ i =Λ 0 (1+aε zz );
in the formula, Λ 0 Is the initial grating period, epsilon zz And (b) taking the axial average strain of the ith section, wherein a is a grating strain coefficient, and the expression is as follows:
in the formula, n eff0 Is the average effective refractive index of FBG in free state, v is the modulation depth, p 11 And p 12 Respectively setting an optical transfer matrix of each grating segment to generate a 2 x 2T matrix T for effective optical stress tensor components based on modal coupling theory i The following were used:
in the formula, R i And S i Is the amplitude of the forward transmission mode and the amplitude of the backward transmission mode of the i-th section, T i The expression satisfies the following relationship:
where Δ z is the length of the grating segment and γ is an intermediate calculation parameter,and κ is the dc self-coupling coefficient and ac self-coupling coefficient of the i-th segment, respectively, the expressions are:
wherein π is constant and δ n eff Defining gamma as the mean of the refractive index variation over the grating periodObtaining T matrix T ═ T of all grating segments by iterative computation N T N-1 T N-2 …T 2 T 1 Then, the following are obtained:
in the formula, R 0 、S 0 Is the amplitude of the forward transmission mode and the amplitude of the backward transmission mode of the first section, R L 、S L The amplitude of the forward transmission mode and the amplitude of the backward transmission mode of the tail section;
the reflectivity r of the FBG for each wavelength is calculated as follows:
where ρ is the reflection coefficient of the fiber Bragg grating, T 21 And T 11 Is the intermediate calculation result for calculating the reflectivity r, which respectively represents the 1 st number of the 2 nd row and the 1 st number of the 1 st row in the T matrix;
and obtaining the fiber Bragg grating reflection spectrum in the whole wavelength interval through iterative calculation, and calculating the corresponding fiber Bragg grating reflection spectrum under different crack lengths to obtain the reflection spectrum of the whole crack expansion process.
3. The method for diagnosing fiber bragg grating-based cracks according to claim 1, wherein the extracting damage sensitive characteristic values and reference signals in the reflection spectrum comprises:
representing the change of the reflection spectrum along with crack propagation by wavelength shift, wherein the wavelength shift is obtained by subtracting the central wavelength of a reference signal from the central wavelength under crack damage, and extracting the central wavelength lambda of the reflection spectrum c Satisfies the following relation:
wherein λ represents a reflectance spectrum wavelength, and R (λ) represents a reflectance;
the variation of the different crack length spreads during crack propagation is characterized by the full width of the reflection line at a reflectance value L, wherein the spread b satisfies the following relationship:
b=|λ L1 -λ Lend |;
in the formula, λ L1 Is the wavelength corresponding to the first limit of reflectivity L, and λ Lend The last wavelength for which the reflectivity is a limit value L;
the wave crest number is used for representing the change of the wave crest numbers of different crack lengths in the crack propagation process;
and representing the area change of the reflection spectrums with different crack lengths in the crack propagation process by using the area of the shape surrounded by the reflection spectrum line and the horizontal axis of the coordinate, wherein the area S expression of the normalized reflection spectrums meets the following relational expression:
in the formula, λ 1 Is the first wavelength, λ, of the wavelength range of the reflectance spectrum end Is to turn overThe last wavelength of the spectral wavelength range,represents a reflection line;
representing the change of the overlapping area of different crack lengths in the crack propagation process by using the overlapping part of the crack-free reference reflection spectral line and each reflection spectrum in the crack propagation process, wherein the overlapping area S is normalized c The expression satisfies the following relation:
in the formula of lambda s Is the first wavelength, λ, of the coincident range of the wavelengths of the reflection lines e Is the last wavelength of the coincidence range of the wavelengths of the reflection lines,representing coincident reflection lines;
representing the change of correlation coefficients of different crack lengths in the crack propagation process by using the crack-free reference reflection spectral line and the damage spectrum in the crack propagation process, wherein the correlation coefficient C m Satisfies the following relation:
in the formula, ρ 0 As a reference spectral reflectance vector, ρ m As the spectral reflectance vector of the damage, N R Is the length of the reflectivity vector and,for any wavelength shift from the damage spectrum, k represents the number of reflectivities in the reflectivity vector, k being 1,2 R ;
Characterizing the significance degree of reflection spectrum chirp phenomena of different crack lengths in a crack propagation process by a fractal dimension, wherein the fractal dimension FD is calculated by a box counting method, and an expression formula meets the following relational expression:
wherein r is the side length of the grid dividing the image, wherein r is 2 i I is 0,1,2, 10, m is the number of grids.
4. The fiber bragg grating based crack diagnosis method as claimed in claim 1, wherein the constructing a crack length regression model with the damage sensitivity characteristic value as an input and the crack length as an output includes:
acquiring a crack length sample data set { (x) i ,y i 1,2, …, n is the number of samples with different crack lengths, x i ,y i Respectively an input quantity and an output quantity, where x i ={x i1 ,x i2 ,…,x ip To influence crack length y i Damage sensitive eigenvalue of, y i Is the true crack length for the ith sample, p is the effect y i Wherein the estimation function of the crack length satisfies the following relation:
f(x)=ωφ(x)+b;
in the formula, omega is weight, phi (x) is high-dimensional nonlinear function, and b is offset;
solving ω and b requires minimizing the optimization model and introducing the relaxation factor ξ as follows:
y i -ωx i -b≤ε+ξ i ;ωx i -y i +b≤ε+ξ i * ;
wherein C is a penalty factor xi i 、ξ i * Introducing alpha for relaxation variable, epsilon is loss function, N is number of training samples i ,The crack length regression model obtained by establishing a Lagrange function by the Lagrange operator to solve the dual problem meets the following relational expression:
5. The fiber bragg grating based crack diagnosis method as claimed in claim 4, further comprising: and optimizing the penalty factor C and the variance g of the radial basis kernel function by using a cross-validation method to optimize the crack length regression model.
6. A crack diagnosis system based on fiber bragg grating, comprising M fiber bragg grating sensors arranged at a target structure of a train, wherein M is a positive integer, and a processing center connected to the M sensors, wherein the processing center is configured to:
determining material parameters of a manufacturing material of a structure to be analyzed;
establishing a finite element simulation model according to the material parameters, simulating the crack propagation condition based on the finite element simulation model, and acquiring structural strain data of the crack propagation detected by the fiber Bragg grating sensor under different lengths;
reconstructing the structural strain data into a reflection spectrum by adopting a transmission matrix method;
extracting damage sensitive characteristic values and reference signals in the reflection spectrum, wherein the damage sensitive characteristic values are used for representing the length of the crack, and the reference signals are used for representing no expansion crack;
and taking the damage sensitivity characteristic value as an input, taking the crack length as an output to construct a crack length regression model, and diagnosing the crack based on the crack length regression model.
7. A fiber bragg grating based crack diagnosis system comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the steps of the method as claimed in any one of the preceding claims 1 to 6 when executing the computer program.
8. A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the method steps of any one of claims 1 to 6.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210544130.6A CN114923420B (en) | 2022-05-18 | 2022-05-18 | Crack diagnosis method and system based on fiber Bragg grating and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210544130.6A CN114923420B (en) | 2022-05-18 | 2022-05-18 | Crack diagnosis method and system based on fiber Bragg grating and storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114923420A true CN114923420A (en) | 2022-08-19 |
CN114923420B CN114923420B (en) | 2022-12-27 |
Family
ID=82808881
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210544130.6A Active CN114923420B (en) | 2022-05-18 | 2022-05-18 | Crack diagnosis method and system based on fiber Bragg grating and storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114923420B (en) |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001059797A (en) * | 1999-08-24 | 2001-03-06 | Mitsubishi Heavy Ind Ltd | Optical fiber underlaying method and distortion detecting device through use of optical fiber |
JP2005274200A (en) * | 2004-03-23 | 2005-10-06 | Nippon Steel Weld Prod & Eng Co Ltd | Optical fiber strain sensor |
EP1826545A2 (en) * | 2006-02-27 | 2007-08-29 | Fuji Jukogyo Kabushiki Kaisha | Damage detection system using optical fibre sensors |
CN101105404A (en) * | 2006-07-11 | 2008-01-16 | 上海理工大学 | Method for synchronously measuring temperature and stress using single optical grating |
CN102313779A (en) * | 2011-07-22 | 2012-01-11 | 北京航空航天大学 | Acoustic emission signal sensing system based on fiber Bragg grating |
JP2013213710A (en) * | 2012-04-02 | 2013-10-17 | Kawasaki Heavy Ind Ltd | Minute crack sensor of concrete |
CN103528599A (en) * | 2013-09-29 | 2014-01-22 | 哈尔滨理工大学 | Spectral discrimination method based on non-uniform fiber Bragg grating distributed sensing |
CN103575664A (en) * | 2013-10-08 | 2014-02-12 | 南京航空航天大学 | Structural multiscale health monitoring device and structural multiscale health monitoring method on basis of novel integrated type composite sensor array |
CN105372182A (en) * | 2015-11-03 | 2016-03-02 | 中国科学院半导体研究所 | Railway rail transverse crack monitoring method based on fiber Bragg grating (FGB) array |
CN105954293A (en) * | 2016-05-03 | 2016-09-21 | 北京航空航天大学 | Fiber Bragg grating sensor-based local crack real-time monitoring method |
CN108333043A (en) * | 2017-12-14 | 2018-07-27 | 中国特种飞行器研究所 | A kind of composite repair structural crack method of real-time based on FBG sensor |
CN108562490A (en) * | 2018-04-13 | 2018-09-21 | 武汉理工大学 | Crack Damage monitoring based on FBG sensor array spectral and strain field measurement method |
CN109916741A (en) * | 2019-01-18 | 2019-06-21 | 昆明理工大学 | A kind of composite impact crack detection algorithm based on optical fiber grating sensing |
CN110095470A (en) * | 2019-04-15 | 2019-08-06 | 北京航空航天大学 | A kind of Crack Damage quantitative approach based on fiber-optic grating sensor |
CN110567378A (en) * | 2019-09-04 | 2019-12-13 | 中南大学 | optical fiber displacement sensor and measuring method thereof |
-
2022
- 2022-05-18 CN CN202210544130.6A patent/CN114923420B/en active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001059797A (en) * | 1999-08-24 | 2001-03-06 | Mitsubishi Heavy Ind Ltd | Optical fiber underlaying method and distortion detecting device through use of optical fiber |
JP2005274200A (en) * | 2004-03-23 | 2005-10-06 | Nippon Steel Weld Prod & Eng Co Ltd | Optical fiber strain sensor |
EP1826545A2 (en) * | 2006-02-27 | 2007-08-29 | Fuji Jukogyo Kabushiki Kaisha | Damage detection system using optical fibre sensors |
CN101105404A (en) * | 2006-07-11 | 2008-01-16 | 上海理工大学 | Method for synchronously measuring temperature and stress using single optical grating |
CN102313779A (en) * | 2011-07-22 | 2012-01-11 | 北京航空航天大学 | Acoustic emission signal sensing system based on fiber Bragg grating |
JP2013213710A (en) * | 2012-04-02 | 2013-10-17 | Kawasaki Heavy Ind Ltd | Minute crack sensor of concrete |
CN103528599A (en) * | 2013-09-29 | 2014-01-22 | 哈尔滨理工大学 | Spectral discrimination method based on non-uniform fiber Bragg grating distributed sensing |
CN103575664A (en) * | 2013-10-08 | 2014-02-12 | 南京航空航天大学 | Structural multiscale health monitoring device and structural multiscale health monitoring method on basis of novel integrated type composite sensor array |
CN105372182A (en) * | 2015-11-03 | 2016-03-02 | 中国科学院半导体研究所 | Railway rail transverse crack monitoring method based on fiber Bragg grating (FGB) array |
CN105954293A (en) * | 2016-05-03 | 2016-09-21 | 北京航空航天大学 | Fiber Bragg grating sensor-based local crack real-time monitoring method |
CN108333043A (en) * | 2017-12-14 | 2018-07-27 | 中国特种飞行器研究所 | A kind of composite repair structural crack method of real-time based on FBG sensor |
CN108562490A (en) * | 2018-04-13 | 2018-09-21 | 武汉理工大学 | Crack Damage monitoring based on FBG sensor array spectral and strain field measurement method |
CN109916741A (en) * | 2019-01-18 | 2019-06-21 | 昆明理工大学 | A kind of composite impact crack detection algorithm based on optical fiber grating sensing |
CN110095470A (en) * | 2019-04-15 | 2019-08-06 | 北京航空航天大学 | A kind of Crack Damage quantitative approach based on fiber-optic grating sensor |
CN110567378A (en) * | 2019-09-04 | 2019-12-13 | 中南大学 | optical fiber displacement sensor and measuring method thereof |
Also Published As
Publication number | Publication date |
---|---|
CN114923420B (en) | 2022-12-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pan et al. | Moving force identification based on redundant concatenated dictionary and weighted l1-norm regularization | |
US20150057988A1 (en) | Method for Identifying Stochastic Information of Heterogeneous Materials | |
Catbas et al. | Predictive analysis by incorporating uncertainty through a family of models calibrated with structural health-monitoring data | |
KR102023701B1 (en) | Apparatus and Method for Measuring Distribution of Structure Strain with Combination of Discrete Sensor and Distributed Sensor | |
Sankararaman et al. | Statistical inference of equivalent initial flaw size with complicated structural geometry and multi-axial variable amplitude loading | |
Gan et al. | Estimation of remaining fatigue life under two-step loading based on kernel-extreme learning machine | |
CN111024821A (en) | Composite material storage box health monitoring system and method | |
CN117540174B (en) | Building structure multi-source heterogeneous data intelligent analysis system and method based on neural network | |
Yuan et al. | Equivalence of partial-least-squares SEM and the methods of factor-score regression | |
CN108846200B (en) | Quasi-static bridge influence line identification method based on iteration method | |
Du et al. | Three-dimensional stress wave imaging of wood internal defects using TKriging method | |
Sierra-Pérez et al. | Damage detection in composite materials structures under variable loads conditions by using fiber Bragg gratings and principal component analysis, involving new unfolding and scaling methods | |
Lin et al. | An integrated neuro-genetic approach incorporating the Taguchi method for product design | |
CN114896712A (en) | Structural damage identification method based on reference Bayesian principle and sparse regularization | |
CN116703323A (en) | Bridge cantilever construction monitoring system optimization method based on grey prediction theory | |
Lu et al. | Structural response estimation method based on particle swarm optimisation/support vector machine and response correlation characteristics | |
Zhu et al. | Structural damage detection of the bridge under moving loads with the quasi-static displacement influence line from one sensor | |
Seemab et al. | Automated detection of propagating cracks in RC beams without shear reinforcement based on DIC-controlled modeling of damage localization | |
CN114923420B (en) | Crack diagnosis method and system based on fiber Bragg grating and storage medium | |
Shen et al. | Prediction of critical strains of flexible pavement from traffic speed deflectometer measurements | |
CN113723496A (en) | Reactor iron core loosening fault diagnosis method, device and equipment | |
CN116227045B (en) | Local stress strain field construction method and system for structural test piece | |
Yang et al. | A multi-sensor mapping Bi-LSTM model of bridge monitoring data based on spatial-temporal attention mechanism | |
CN113973403B (en) | Temperature-induced strain field redistribution intelligent sensing method based on structure discrete measurement point topology | |
CN114254533B (en) | Method for examining influence and prediction of fatigue vibration on fixed angle of product group component |
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 |