Error correction method for pipeline morphology inversion based on strain correction coefficient matrix optimization
Through the strain correction coefficient matrix optimization method, the problem of cumulative error in pipeline morphology inversion is solved, and higher precision pipeline morphology monitoring is achieved, and the error is significantly reduced.
Patent Information
- Application Number
- CN202310167388.3
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2023-02-27
- Publication Date
- 2025-08-29
- Estimated Expiration
- 2043-02-27
AI Technical Summary
The prior art fails to effectively correct the cumulative error during the pipeline morphology inversion process, resulting in large deviations in the pipeline morphology results, affecting monitoring accuracy.
The strain correction coefficient matrix optimization method is adopted, and the strain correction coefficient matrix is defined by design variables, the objective function is optimized, and the fiber sensor measurement data and boundary conditions are combined, and the strain matrix is corrected by an optimization algorithm to reduce errors.
The pipeline morphology inversion error is significantly reduced, the deflection error accumulation effect is reduced, and the monitoring accuracy is improved. The error is reduced from 1.243mm to 0.257mm and 0.417mm.
Smart Images
Figure CN116384057B_ABST
Abstract
Description
Technical Field
[0001] The present invention belongs to the field of structural health monitoring based on optical fiber sensors, and in particular relates to a pipeline morphology inversion error correction method based on strain correction coefficient matrix optimization. Background Art
[0002] Due to the strong adaptability of optical fiber to harsh monitoring environments, its small size, and the application of time division multiplexing, wavelength division multiplexing, and space division multiplexing technologies in the optical fiber field, optical fiber sensors have good engineering applicability for pipeline morphology inversion.
[0003] In terms of morphological inversion algorithms, Frenet's moving frame theory is an important method for morphological reconstruction based on fiber optic sensors. Wu Yaxing and others conducted in-depth research on the pipeline morphology inversion method based on the minimum rotation frame, which effectively overcomes the shortcomings of the Frenet moving frame theory in terms of instantaneous inflection points. The overall frame transition tends to be smooth, and it can better describe the morphological characteristics of spatial curves. However, when this method is applied to a single-core optical fiber, it requires the use of multiple sensing paths, which increases the complexity of the sensing network. Therefore, it is urgent to study new strain sensing layout forms suitable for pipeline configurations to improve the efficiency of sensing strain information on the surface of the pipeline structure. Current research on pipeline morphology monitoring and inversion rarely involves the correction of accumulated errors in the pipeline morphology reconstruction process. Considering that the accumulated errors corresponding to the actual inversion process are unavoidable, if they are not corrected, it will lead to large deviations in the pipeline morphology results. Therefore, reducing and correcting the accumulated errors in the pipeline morphology reconstruction process has become an urgent problem that researchers need to solve. Summary of the Invention
[0004] Purpose of the invention: The technical problem to be solved by the present invention is to address the deficiencies of the prior art and provide a pipeline morphology inversion error correction method based on strain correction coefficient matrix optimization, comprising the following steps:
[0005] Step 1: The design variables are defined as the strain correction coefficient matrix, and the strain correction coefficient matrix optimization objective function is determined;
[0006] Step 2: collect constraint information for the strain correction coefficient matrix optimization method;
[0007] Step 3: Carry out morphological inversion error correction based on the strain correction coefficient matrix optimization method.
[0008] Step 1 includes:
[0009] The strain measured by the distributed optical fiber sensor is the strain at high-density discrete measuring points. The strain measured at n discrete measuring points is set to be The true strain values corresponding to n discrete measuring points are ε1,ε2,…,ε n , respectively define the measurement strain matrix and the true strain matrix
[0010]
[0011]
[0012] Where n is the number of strain sensing points of the distributed optical fiber sensor. Under ideal conditions, the relationship is as follows:
[0013]
[0014] The design variables are defined as the strain correction coefficient matrix K:
[0015] K=diag(k1,k2,…,k n )
[0016] Among them, k i Represents the i-th strain correction coefficient, which is currently an undetermined coefficient. The diag function is used to construct a diagonal matrix;
[0017] Construct an approximate strain correction matrix Used to replace the true strain matrix
[0018]
[0019] By optimizing the strain correction coefficient matrix K, With the measured strain matrix With the best similarity, on this basis, the objective function F of strain correction coefficient matrix optimization is defined K for:
[0020]
[0021] Where min refers to finding the minimum value of the function.
[0022] Step 2 includes:
[0023] If there is a point on the pipeline whose deflection does not change with the change of the overall deflection of the pipeline, the point is called an anchor point. If the structure measured by the inclinometer tube has a rigid support, the deflection at the support is always 0, and the support is the anchor point.
[0024] Based on the deflection and the pipeline constraints, the pipeline anchor points are set as points A, B, ..., X from left to right, and the boundary conditions of the inclinometer pipe with anchor points are determined as follows:
[0025]
[0026] Where w(x) is the deflection of the pipe at a distance x from the left end, x A 、xB ,…,x X are the distances between points A, B, ..., X and the left end of the pipeline respectively;
[0027] Based on the rotation angle, when the anchor point feature is unknown, the rotation angle measured by the accelerometer is used. Assuming that the accelerometer setting points are A, B, ..., X from left to right, the constraint condition expression of the inclinometer casing is obtained and regarded as the boundary condition:
[0028]
[0029] Where θ(x) is the angle of rotation between the pipe and the left end at a distance of x, θ A ,θ B ,…,θ X are the rotation angles at points A, B, ..., and X, measured by the accelerometer.
[0030] Step 3 includes:
[0031] Case 1: Obtaining boundary conditions based on the inclinometer casing with anchor points;
[0032] Case 2: Obtaining the boundary conditions of the inclinometer casing based on unknown anchor points.
[0033] In step 3, the step of obtaining the boundary conditions of the inclinometer casing with anchor points specifically includes: using the constraint conditions of the optimization algorithm to equivalently replace the boundary conditions of the inclinometer casing:
[0034]
[0035] w(x) is the pipe deflection curve equation, is the deflection equation obtained by inversion after correction;
[0036] Assume W(x,K) is the correction equation for the deflection of the inclinometer tube with anchor point:
[0037]
[0038] in, Refers to the corrected relationship between pipeline strain and x, Refers to Functional relationship between the fitted corrected strain and x;
[0039] Substitute x = x A 、x=x B ,…,x=x C , establish the following mathematical model for the strain correction coefficient matrix optimization problem:
[0040]
[0041] Solving the mathematical model to obtain the optimal strain correction coefficient matrix K;
[0042] according to Obtain the optimized strain correction matrix On this basis, based on the strain correction matrix The modified deflection curve equation obtained by inversion using the strain-curvature-deflection multiple integration algorithm is
[0043]
[0044] Where dx is the differential of x, C and D are constants determined by the boundary conditions, and R is the radius of the pipe.
[0045] In step 3, the acquisition of the inclinometer casing boundary conditions based on the unknown anchor point specifically includes:
[0046] The constraints of the optimization algorithm are used to replace the boundary conditions of the inclinometer casing:
[0047]
[0048] Among them, θ(x) is the angle equation of the pipeline, θ(x A ),θ(x B ),…,θ(x X ) is measured using an accelerometer, is the rotation angle equation obtained by inversion after correction;
[0049] Assume Theta(x,K) is the correction equation for the inclinometer tube rotation angle when the anchor point characteristics are unknown:
[0050]
[0051] Where x represents the distance from the left end of the pipe, K is the strain correction coefficient matrix, and R is the pipe radius. Refers to the corrected relationship between pipeline strain and x, Refers to Functional relationship between the fitted corrected strain and x;
[0052] Substitute x = x A 、x=x B ,…,x=x C , establish the following mathematical model of strain correction coefficient matrix optimization problem:
[0053]
[0054] Solve the above mathematical model to obtain the optimal strain correction coefficient matrix K;
[0055] according to Obtain the optimized strain correction matrix On this basis, based on the strain correction matrix The modified multi-sensor fusion inclinometer tube deflection curve equation is obtained by inversion using the strain-curvature-deflection multiple integration algorithm.
[0056]
[0057] The present invention also provides a storage medium storing a computer program or instruction. When the computer program or instruction is executed, the pipeline morphology inversion error correction method based on strain correction coefficient matrix optimization is implemented.
[0058] The present invention has the following beneficial effects: Strain measurement error is the main source of morphological inversion error. Therefore, this method proposes to use the strain correction coefficient matrix optimization method to correct it. This method uses the constraint conditions of the optimization algorithm to equivalently replace the boundary conditions of the measured structure. At the same time, it is necessary to establish a strain matrix that is consistent with the measured strain matrix. Strain correction matrix with optimal similarity On this basis, based on the strain-curvature-deflection multi-integral method, a more accurate deflection equation can be inverted Used to describe the pipeline morphology. Combined with the results of two numerical simulations of a double-end-clamped, indeterminate pipeline solid finite element model, the root mean square error of the corrected pipeline morphology inversion based on the strain correction coefficient matrix optimization was reduced from 1.243mm to 0.257mm and from 0.733mm to 0.417mm, respectively. The research results show that the strain correction coefficient matrix optimization method helps to reduce the pipeline morphology inversion error. Compared with the inclinometer tube morphology inversion method based only on fiber optic sensors and without correction, the inclinometer tube morphology error correction method based on the strain correction coefficient matrix optimization method can effectively reduce the cumulative effect of deflection errors in the morphology inversion process of the conventional strain-curvature-deflection multiple integration method, resulting in a significant reduction in the pipeline morphology inversion error, proving that this method has extremely high application value. BRIEF DESCRIPTION OF THE DRAWINGS
[0059] The present invention will be further described below in conjunction with the accompanying drawings and specific embodiments, and the above and / or other advantages of the present invention will become more apparent.
[0060] Figure 1a It is a schematic diagram of the boundary conditions of the pipeline with anchor points and their positional relationships.
[0061] Figure 1b It is a schematic diagram of the applied load and its position relationship.
[0062] Figure 2a Schematic diagram of pipeline boundary conditions and applied loads.
[0063] Figure 2b It is a schematic diagram of the positions of known corner points A, B, and C.
[0064] Figure 3a It is a schematic diagram of the inversion morphological error correction effect with anchor points.
[0065] Figure 3b It is a schematic diagram of the inversion morphological error correction effect when the anchor point position is unknown.
[0066] Figure 4 This is a flow chart of the pipeline morphology inversion error correction method based on strain correction coefficient matrix optimization. DETAILED DESCRIPTION
[0067] like Figure 4 As shown, the present invention provides a pipeline morphology inversion error correction method based on strain correction coefficient matrix optimization, comprising the following steps:
[0068] Step 1: The design variables are defined as the strain correction coefficient matrix, and the strain correction coefficient matrix optimization objective function is determined;
[0069] Step 2: collect constraint information for the strain correction coefficient matrix optimization method;
[0070] Step 3: Carry out morphological inversion error correction based on the strain correction coefficient matrix optimization method.
[0071] Step 1 includes:
[0072] The strain measured by the distributed optical fiber sensor is the strain at high-density discrete measuring points. The strain measured at n discrete measuring points is set to be The true strain values corresponding to n discrete measuring points are ε1,ε2,…,ε n , respectively define the measurement strain matrix and the true strain matrix
[0073]
[0074]
[0075] Where n is the number of strain sensing points of the distributed optical fiber sensor. Under ideal conditions, the relationship is as follows:
[0076]
[0077] The design variables are defined as the strain correction coefficient matrix K:
[0078] K=diag(k1,k2,…,k n )
[0079] Among them, k i Represents the i-th strain correction coefficient, which is currently an undetermined coefficient. The diag function is used to construct a diagonal matrix;
[0080] Construct an approximate strain correction matrix Used to replace the true strain matrix
[0081]
[0082] By optimizing the strain correction coefficient matrix K, With the measured strain matrix With the best similarity, on this basis, the objective function F of strain correction coefficient matrix optimization is defined K for:
[0083]
[0084] Where min refers to finding the minimum value of the function.
[0085] Step 2 includes:
[0086] If there is a point on the pipeline whose deflection does not change with the change of the overall deflection of the pipeline, the point is called an anchor point. If the structure measured by the inclinometer tube has a rigid support, the deflection at the support is always 0, and the support is the anchor point.
[0087] Based on the deflection and the pipeline constraints, the pipeline anchor points are set as points A, B, ..., X from left to right, and the boundary conditions of the inclinometer pipe with anchor points are determined as follows:
[0088]
[0089] Where w(x) is the deflection of the pipe at a distance x from the left end, x A 、x B ,…,x X are the distances between points A, B, ..., X and the left end of the pipeline respectively;
[0090] Based on the rotation angle, when the anchor point feature is unknown, the rotation angle measured by the accelerometer is used. Assuming that the accelerometer setting points are A, B, ..., X from left to right, the constraint condition expression of the inclinometer casing is obtained and regarded as the boundary condition:
[0091]
[0092] Where θ(x) is the angle of rotation between the pipe and the left end at a distance of x, θ A ,θ B ,…,θ X are the rotation angles at points A, B, ..., and X, measured by the accelerometer.
[0093] Step 3 includes:
[0094] Case 1: Obtaining boundary conditions based on the inclinometer casing with anchor points;
[0095] Case 2: Obtaining the boundary conditions of the inclinometer casing based on unknown anchor points.
[0096] In step 3, the step of obtaining the boundary conditions of the inclinometer casing with anchor points specifically includes: using the constraint conditions of the optimization algorithm to equivalently replace the boundary conditions of the inclinometer casing:
[0097]
[0098] w(x) is the pipe deflection curve equation, is the deflection equation obtained by inversion after correction;
[0099] Assume W(x,K) is the correction equation for the deflection of the inclinometer tube with anchor point:
[0100]
[0101] in, Refers to the corrected relationship between pipeline strain and x, Refers to Functional relationship between the fitted corrected strain and x;
[0102] Substitute x = x A 、x=x B ,…,x=x C , establish the following mathematical model for the strain correction coefficient matrix optimization problem:
[0103]
[0104] Solving the mathematical model to obtain the optimal strain correction coefficient matrix K;
[0105] according to Obtain the optimized strain correction matrix On this basis, based on the strain correction matrix The modified deflection curve equation obtained by inversion using the strain-curvature-deflection multiple integration algorithm is
[0106]
[0107] Where dx is the differential of x, C and D are constants determined by the boundary conditions, and R is the radius of the pipe.
[0108] In step 3, the acquisition of the inclinometer casing boundary conditions based on the unknown anchor point specifically includes:
[0109] The constraints of the optimization algorithm are used to replace the boundary conditions of the inclinometer casing:
[0110]
[0111] Among them, θ(x) is the angle equation of the pipeline, θ(x A ),θ(x B ),…,θ(x X ) is measured using an accelerometer, is the rotation angle equation obtained by inversion after correction;
[0112] Assume Theta(x,K) is the correction equation for the inclinometer tube rotation angle when the anchor point characteristics are unknown:
[0113]
[0114] Where x represents the distance from the left end of the pipe, K is the strain correction coefficient matrix, and R is the pipe radius. Refers to the corrected relationship between pipeline strain and x, Refers to Functional relationship between the fitted corrected strain and x;
[0115] Substitute x = x A 、x=x B ,…,x=x C , establish the following mathematical model of strain correction coefficient matrix optimization problem:
[0116]
[0117] Solve the above mathematical model to obtain the optimal strain correction coefficient matrix K;
[0118] according to Obtain the optimized strain correction matrix On this basis, based on the strain correction matrix The modified multi-sensor fusion inclinometer tube deflection curve equation is obtained by inversion using the strain-curvature-deflection multiple integration algorithm.
[0119]
[0120] The present invention also provides a storage medium storing a computer program or instruction. When the computer program or instruction is executed, the pipeline morphology inversion error correction method based on strain correction coefficient matrix optimization is implemented.
[0121] Example
[0122] For the finite element model of a left-end clamped inclinometer tube with an outer diameter of 70 mm, an inner diameter of 60 mm, a length of l = 3000 mm, a material elastic modulus of 3 GPa, and a Poisson's ratio of 0.33, the anchor points are set at 750 mm, 1500 mm, and 2250 mm from the clamped end, respectively. Concentrated forces F1 = 400 N and F2 = 200 N are applied at distances X1 = 500 mm and X2 = 2000 mm from the clamped end, respectively.
[0123] like Figure 4 As shown, this embodiment provides a pipeline morphology inversion error correction method based on strain correction coefficient matrix optimization, including the following steps:
[0124] Step 1: The design variable is the strain correction coefficient matrix, and then the strain correction coefficient matrix optimization objective function is determined:
[0125] The strain measured by the distributed optical fiber sensor is the strain at high-density discrete measuring points. Assume that the strain measured at each discrete measuring point is The true strain values corresponding to each discrete measuring point are ε1,ε2,…,ε n , then the measurement strain matrix can be defined separately and the true strain matrix
[0126]
[0127]
[0128] Where n is the number of strain sensing points of the distributed optical fiber sensor. Under ideal conditions, the relationship between the two is:
[0129]
[0130] The design variables are defined as the strain correction coefficient matrix K:
[0131] K=diag(k1,k2,…,k n )
[0132] Among them, k i represents the i-th strain correction coefficient, which is currently an undetermined coefficient. The diag function is used to construct a diagonal matrix. In this embodiment, m = 101;
[0133] Construct an approximate strain correction matrix Used to replace the true strain matrix
[0134]
[0135] By optimizing the strain correction coefficient matrix K, With the measured strain matrix With the best similarity, on this basis, the objective function F of strain correction coefficient matrix optimization is defined K for:
[0136]
[0137] Where min refers to finding the minimum value of the function.
[0138] Step 2: Collect constraint information for the strain correction coefficient matrix optimization method:
[0139] Method 1: Based on deflection and pipeline constraints, the boundary conditions of the inclinometer pipe with anchor points are determined as follows:
[0140]
[0141] Method 2: Based on the rotation angle, when the anchor point characteristics are unknown, the rotation angle measured by the accelerometer is used to obtain the inclinometer tube constraint expression and treat it as a boundary condition:
[0142]
[0143] Step 3: Carry out morphological inversion error correction based on the strain correction coefficient matrix optimization method:
[0144] Case 1: Boundary conditions based on an inclinometer casing with anchor points
[0145] The constraints of the optimization algorithm are used to replace the boundary conditions of the inclinometer casing:
[0146]
[0147] in, is the strain distribution after correction by the strain correction coefficient matrix, w(x) is the pipe deflection curve equation, x A , x B , x C are the distances between points A, B, and C and the fixed support end, R is the radius of the inclinometer tube, They represent the corrected deflections of points A, B, and C. In this embodiment, x A =750mm, x B =1500mm, x C =2250mm, such as Figure 1a 、 Figure 1b shown.
[0148] Assume W(x,K) is the correction equation for the deflection of the inclinometer tube with anchor point:
[0149]
[0150] Where x represents the position relative to the fixed end, and K is the strain correction coefficient matrix.
[0151] Based on the above content, a mathematical model for the strain correction coefficient matrix optimization problem is established:
[0152]
[0153] Solve the above model and get the optimal strain correction coefficient matrix K. Obtain the optimized strain correction matrix On this basis, based on the strain correction matrix The modified deflection curve equation obtained by inversion using the strain-curvature-deflection multiple integration algorithm is like Figure 3a shown.
[0154] Case 2: Boundary conditions of the inclinometer casing with unknown anchor points
[0155] The constraints of the optimization algorithm are used to replace the boundary conditions of the inclinometer casing:
[0156]
[0157] in, is the strain distribution after correction by the strain correction coefficient matrix, R is the radius of the inclinometer tube, such as Figure 2a 、 Figure 2b shown.
[0158] Assume Theta(x,K) is the correction equation for the inclinometer tube rotation angle when the anchor point characteristics (including number and position) are unknown:
[0159]
[0160] Where x represents the position relative to the fixed end, and K is the strain correction coefficient matrix.
[0161] According to the above content, the mathematical model of strain correction coefficient matrix optimization problem is established:
[0162]
[0163] Solving the above model, we can get the optimal strain correction coefficient matrix K. Obtain the optimized strain correction matrix
[0164] On this basis, based on the strain correction matrix The modified multi-sensor fusion inclinometer tube deflection curve equation is obtained by inversion using the strain-curvature-deflection multiple integration algorithm. like Figure 3b shown.
[0165] In a specific implementation, the present application provides a computer storage medium and a corresponding data processing unit, wherein the computer storage medium is capable of storing a computer program that, when executed by the data processing unit, executes the invention disclosure of the pipeline morphology inversion error correction method based on strain correction coefficient matrix optimization provided by the present invention, as well as some or all of the steps in each embodiment. The storage medium may be a magnetic disk, an optical disk, a read-only memory (ROM), or a random access memory (RAM).
[0166] Those skilled in the art can clearly understand that the technical solutions in the embodiments of the present invention can be implemented by means of computer programs and their corresponding general hardware platforms. Based on this understanding, the technical solutions in the embodiments of the present invention, in essence or in other words, the part that contributes to the prior art, can be embodied in the form of a computer program, i.e., a software product. The computer program software product can be stored in a storage medium and includes several instructions for enabling a device including a data processing unit (which can be a personal computer, a server, a single-chip microcomputer, a MUU, or a network device, etc.) to execute the methods described in various embodiments of the present invention or certain parts of the embodiments.
[0167] The present invention provides a pipeline morphology inversion error correction method based on strain correction matrix optimization. There are numerous methods and approaches for implementing this technical solution. The above description is merely a preferred embodiment of the present invention. It should be noted that those skilled in the art may make various improvements and modifications without departing from the principles of the present invention, and such improvements and modifications are also within the scope of protection of the present invention. Any components not specified in this embodiment may be implemented using existing technologies.
Claims
1. A pipeline morphology inversion error correction method based on strain correction coefficient matrix optimization is characterized by: The following steps are involved: Step 1: The design variables are defined as the strain correction coefficient matrix, and the strain correction coefficient matrix optimization objective function is determined; Step 2: collect constraint information for the strain correction coefficient matrix optimization method; Step 3: Carry out morphological inversion error correction based on the strain correction coefficient matrix optimization method; Step 1 includes: The strain measured by the distributed optical fiber sensor is the strain at high-density discrete measuring points. The strain measured at n discrete measuring points is set to be The true strain values corresponding to n discrete measuring points are ε1,ε2,…,ε n , respectively define the measurement strain matrix and the true strain matrix Where n is the number of strain sensing points of the distributed optical fiber sensor. Under ideal conditions, the relationship is as follows: The design variables are defined as the strain correction coefficient matrix K: K=diag(k1,k2,…,k n ) Among them, k i Represents the i-th strain correction coefficient, which is currently an undetermined coefficient. The diag function is used to construct a diagonal matrix; Construct an approximate strain correction matrix Used to replace the true strain matrix By optimizing the strain correction coefficient matrix K, With the measured strain matrix With the best similarity, on this basis, the objective function F of strain correction coefficient matrix optimization is defined K for: Where min refers to finding the minimum value of the function; Step 2 includes: If there is a point on the pipeline whose deflection does not change with the change of the overall deflection of the pipeline, the point is called an anchor point. If the structure measured by the inclinometer tube has a rigid support, the deflection at the support is always 0, and the support is the anchor point. Based on the deflection and the pipeline constraints, the pipeline anchor points are set as points A, B, ..., X from left to right, and the boundary conditions of the inclinometer pipe with anchor points are determined as follows: Where w(x) is the deflection of the pipe at a distance x from the left end, x A 、x B ,…,x X are the distances between points A, B, ..., X and the left end of the pipeline respectively; Based on the rotation angle, when the anchor point feature is unknown, the rotation angle measured by the accelerometer is used. Assuming that the accelerometer setting points are A, B, ..., X from left to right, the constraint condition expression of the inclinometer casing is obtained and regarded as the boundary condition: Where θ(x) is the angle of rotation between the pipe and the left end at a distance of x, θ A ,θ B ,…,θ X are the rotation angles at points A, B, ..., and X, respectively, measured by the accelerometer; Step 3 includes: Case 1: Obtaining boundary conditions based on the inclinometer casing with anchor points; Case 2: Obtaining the boundary conditions of the inclinometer casing based on unknown anchor points; In step 3, the acquisition of the inclinometer casing boundary conditions based on the unknown anchor point specifically includes: The constraints of the optimization algorithm are used to replace the boundary conditions of the inclinometer casing: Among them, θ(x) is the angle equation of the pipeline, θ(x A ),θ(x B ),…,θ(x X ) is measured using an accelerometer, is the rotation angle equation obtained by inversion after correction; Assume Theta(x,K) is the correction equation for the inclinometer tube rotation angle when the anchor point characteristics are unknown: Where x represents the distance from the left end of the pipe, K is the strain correction coefficient matrix, and R is the pipe radius. Refers to the corrected relationship between pipeline strain and x, Refers to Functional relationship between the fitted corrected strain and x; Substitute x = x A 、x=x B ,…,x=x C , establish the following mathematical model of strain correction coefficient matrix optimization problem: Solve the above mathematical model to obtain the optimal strain correction coefficient matrix K; according to Obtain the optimized strain correction matrix On this basis, based on the strain correction matrix The modified multi-sensor fusion inclinometer tube deflection curve equation is obtained by inversion using the strain-curvature-deflection multiple integration algorithm.
2. The method according to claim 1, characterized in that In step 3, the step of obtaining the boundary conditions of the inclinometer casing with anchor points specifically includes: using the constraint conditions of the optimization algorithm to equivalently replace the boundary conditions of the inclinometer casing: w(x) is the pipe deflection curve equation, is the deflection equation obtained by inversion after correction; Assume W(x,K) is the correction equation for the deflection of the inclinometer tube with anchor point: in, Refers to the corrected relationship between pipeline strain and x, Refers to Functional relationship between the fitted corrected strain and x; Substitute x = x A 、x=x B ,…,x=x C , establish the following mathematical model for the strain correction coefficient matrix optimization problem: Solving the mathematical model to obtain the optimal strain correction coefficient matrix K; according to Obtain the optimized strain correction matrix On this basis, based on the strain correction matrix The modified deflection curve equation obtained by inversion using the strain-curvature-deflection multiple integration algorithm is Where dx is the differential of x, C and D are constants determined by the boundary conditions, and R is the radius of the pipe.
3. A storage medium, characterized in that: A computer program or instruction is stored, and when the computer program or instruction is executed, the method according to any one of claims 1 to 2 is implemented.
Citation Information
Patent Citations
Deep soft soil metro foundation pit soil body parameter inversion analyzing method
CN104077451A
Optical fiber monitoring method for strain-deflection-bending moment state of helicopter blade structure
CN110127078A