CN114777661A - Tunnel section convergence deformation quantitative calculation method based on high-density measuring point strain - Google Patents
Tunnel section convergence deformation quantitative calculation method based on high-density measuring point strain Download PDFInfo
- Publication number
- CN114777661A CN114777661A CN202210369310.5A CN202210369310A CN114777661A CN 114777661 A CN114777661 A CN 114777661A CN 202210369310 A CN202210369310 A CN 202210369310A CN 114777661 A CN114777661 A CN 114777661A
- Authority
- CN
- China
- Prior art keywords
- section
- convergence deformation
- tunnel
- strain
- deformation
- 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/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
- G01B11/18—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge using photoelastic elements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
Abstract
A tunnel section convergence deformation quantitative calculation method based on high-density measuring point strain belongs to the field of tunnel section convergence deformation monitoring. The invention constructs a conversion relation model of the circumferential strain and the convergence deformation of high-density measuring points of the tunnel section, provides a calculation method of the additional horizontal convergence deformation amount of the tunnel section by considering the extrusion effect and the shearing effect, establishes a conversion model of the strain and the horizontal convergence deformation of the tunnel section by utilizing the artificial neural network theory, and provides a calculation method of the horizontal convergence deformation of multiple sections of the tunnel. The method is based on high-density measuring point circumferential strain monitoring data provided by a distributed optical fiber sensing technology, a tunnel section horizontal convergence deformation calculation method considering a longitudinal bending additional effect is established by relying on a conjugate curved beam theory and an artificial neural network, a tunnel section strain-horizontal convergence deformation conversion model is established, the problem of performing tunnel multi-section horizontal convergence deformation calculation by using distributed optical fiber strain monitoring data is solved, and the method has good robustness.
Description
Technical Field
The invention belongs to the field of monitoring of convergence deformation of a tunnel section, and particularly relates to a quantitative calculation method for convergence deformation of a tunnel section based on high-density measuring point strain.
Background
The development and the evolution of the horizontal convergence deformation of the tunnel section can be accurately monitored, and the operation safety of the tunnel structure can be effectively guaranteed. When the full-distributed optical fiber sensing technology is adopted to monitor the horizontal convergence deformation of the section of the tunnel, an ideal optical cable arrangement scheme is to arrange the sensing optical cables along the circumferential direction of the section of the tunnel; however, considering the driving safety of the subway train, the sensing optical cable cannot be arranged in the vault area of the shield tunnel, so that the distributed sensing optical cable arranged on the section cannot adopt a closed arrangement mode along the section. In addition, the existing method for calculating the horizontal convergence deformation of the tunnel section by using the strain of the optical fiber measuring point generally cannot consider the additional horizontal convergence deformation caused by the longitudinal bending of the shield tunnel, so that a certain error exists in the deformation calculation result.
Disclosure of Invention
The invention aims to solve the problems and provides a tunnel section convergence deformation quantitative calculation method based on high-density measuring point strain, which utilizes high-density measuring point annular strain monitoring data provided by a tunnel structure safety distributed optical fiber monitoring system, can effectively consider the influence of additional horizontal convergence deformation, realizes multi-section horizontal convergence deformation calculation, and improves the calculation precision and the calculation efficiency of the tunnel section horizontal convergence deformation.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
the tunnel section convergence deformation quantitative calculation method based on the high-density measuring point strain comprises the following steps:
the method comprises the following steps: constructing a tunnel section circular curved beam model based on a conjugate curved beam theory by using tunnel section high-density measuring point circular strain monitoring data acquired by a full-distributed optical fiber sensor, establishing a tunnel section high-density measuring point circular strain and section convergence deformation conversion model, and calculating to obtain a tunnel section convergence deformation amount;
step two: constructing a basic structure of a tunnel section additional horizontal convergence deformation calculation force method, and calculating the additional horizontal convergence deformation of the section caused under the longitudinal bending effect of the tunnel; on the basis, the additional horizontal convergence deformation of the tunnel section under the extrusion effect caused by longitudinal bending is calculated by considering the constraint action of the actual tunnel section connection rigidity and the foundation reaction force;
step three: constructing a shear effect equation caused by longitudinal bending deformation, calculating the additional horizontal convergence deformation of the tunnel section under the shear effect, and establishing the calculation results of the horizontal convergence deformation of the first step and the second step in parallel to obtain the horizontal convergence deformation of the tunnel section considering the additional deformation influence;
step four: and (4) establishing a tunnel section strain-horizontal convergence deformation conversion model based on the artificial neural network by utilizing the section horizontal convergence deformation obtained by calculation in the step three, and calculating to obtain the tunnel multi-section horizontal convergence deformation.
Compared with the prior art, the invention has the beneficial effects that: the invention constructs a conversion model of the circumferential strain and the section convergence deformation of high-density measuring points of the tunnel section, provides a calculation method of the additional horizontal convergence deformation amount of the tunnel section by considering the extrusion effect and the shearing effect, establishes a conversion model of the strain-horizontal convergence deformation of the tunnel section by utilizing the artificial neural network theory, and provides a calculation method of the horizontal convergence deformation of multiple sections of the tunnel. The method is based on high-density measuring point circumferential strain monitoring data provided by a distributed optical fiber sensing technology, a tunnel section horizontal convergence deformation calculation method considering a longitudinal bending additional effect is established by relying on a conjugate curved beam theory and an artificial neural network, a tunnel section strain-horizontal convergence deformation conversion model is established, the problem of performing tunnel multi-section horizontal convergence deformation calculation by using distributed optical fiber strain monitoring data is solved, the robustness is good, and compared with the existing method, the tunnel section horizontal convergence deformation calculation precision and calculation efficiency based on the strain monitoring data are improved.
Drawings
FIG. 1 is a flow chart of a tunnel section horizontal convergence deformation calculation method based on high-density measuring point strain in the invention;
FIG. 2 is a graph of an example optical fiber strain profile for a tunnel monitoring section;
FIG. 3 is a diagram illustrating the measurement result of the horizontal convergence deformation engineering of the left line of the tunnel cross section according to the embodiment;
FIG. 4 is a diagram illustrating the measurement result of the tunnel cross section right line horizontal convergence deformation engineering of the embodiment;
FIG. 5 is a comparison graph of the calculation results of the horizontal convergence deformation of the left line of the tunnel cross section in the embodiment;
FIG. 6 is a graph comparing the calculation results of the horizontal convergence deformation of the right line of the tunnel cross section in the embodiment;
FIG. 7 is a graph showing the calculation result of the horizontal convergence deformation of the left line of the first ring section of the tunnel according to the embodiment;
FIG. 8 is a graph showing the calculation result of the horizontal convergence deformation of the right line of the first ring section of the tunnel according to the embodiment;
FIG. 9 is an iterative diagram of an embodiment artificial neural network training;
FIG. 10 is a comparison graph of the calculation result and the actual measurement of the horizontal convergence deformation of the multi-section left line of the tunnel according to the embodiment;
FIG. 11 is a graph comparing the calculation results and actual measurement results of the horizontal convergence deformation of the right line of the tunnel multi-section according to the embodiment.
Detailed Description
The technical solution of the present invention is further described below with reference to the embodiments and the drawings, but the present invention is not limited thereto, and any modification or equivalent replacement of the technical solution of the present invention without departing from the spirit of the technical solution of the present invention shall be covered within the protection scope of the present invention.
The first specific implementation way is as follows: the embodiment describes a method for quantitatively calculating convergence deformation of a tunnel section based on strain of high-density measuring points, which comprises the following steps:
the method comprises the following steps: constructing a tunnel section circular curved beam model based on a conjugate curved beam theory by using tunnel section high-density measuring point circular strain monitoring data acquired by a full-distributed optical fiber sensor, establishing a tunnel section high-density measuring point circular strain and section convergence deformation conversion model, and calculating to obtain a tunnel section convergence deformation amount;
step two: constructing a basic structure of a tunnel section additional horizontal convergence deformation calculation force method, and calculating the additional horizontal convergence deformation of the section caused under the longitudinal bending effect of the tunnel; on the basis, considering the constraint action of the actual tunnel section connection rigidity and the foundation reaction force, calculating to obtain the additional horizontal convergence deformation of the tunnel section under the extrusion effect caused by longitudinal bending;
step three: constructing a shear effect equation caused by longitudinal bending deformation, calculating the additional horizontal convergence deformation of the tunnel section under the shear effect, and establishing the calculation results of the horizontal convergence deformation of the first step and the second step in parallel to obtain the horizontal convergence deformation of the tunnel section considering the additional deformation influence;
step four: and (4) establishing a tunnel section strain-horizontal convergence deformation conversion model based on the artificial neural network by using the section horizontal convergence deformation calculated in the third step, and calculating to obtain the tunnel multi-section horizontal convergence deformation.
The invention provides a tunnel section horizontal convergence deformation calculation method based on high-density measuring point strain, which aims at the problem of how to calculate the horizontal convergence deformation of a tunnel section by using distributed optical fiber strain test data distributed annularly. Firstly, constructing a circular curved beam model of a tunnel section by relying on a conjugate curved beam theory, and providing a conversion algorithm of high-density measuring point circular strain and section convergence deformation of the tunnel section; secondly, analyzing the influence rule of the additional internal force caused by the longitudinal bending of the tunnel structure on the horizontal convergence deformation of the section, and providing a calculation method for the additional horizontal convergence deformation of the shield tunnel section under the consideration of the longitudinal bending effect; finally, a multi-section strain-horizontal convergence deformation conversion model based on an artificial neural network is provided by combining a tunnel section 'n' -shaped distributed optical cable arrangement mode; the method has good robustness and anti-noise capability, and can accurately calculate the horizontal convergence deformation of the tunnel section based on the strain monitoring data of the high-density measuring points.
The second embodiment is as follows: in a tunnel section convergence deformation quantitative calculation method based on high-density measurement point strain according to a specific embodiment, in the first step, the calculation method of the tunnel section convergence deformation amount is as follows:
(1) acquiring strain monitoring data of high-density measuring points by using a fully-distributed sensing optical cable distributed along the circumferential direction through length of a tunnel section, and calculating the convergence deformation of the circumferential any position coordinates (p, lambda) (p is 1,2, …, n, and lambda is more than or equal to 0 and less than 1) of a tunnel section duct piece caused by bending moment according to the theory of conjugate curved beams by the following formula, namely
In the formula (I), the compound is shown in the specification,
representing the convergent deformation at the azimuthal coordinate (p, λ) caused by the bending moment; beta represents a corresponding central angle between two circumferential strain measuring points according to the optical fiber; r represents the radius of a segment of the section of the tunnel; h represents the height of the segment of the tunnel section;
respectively indicate the breaking of the fiber distribution positions of i point and p point of the coordinateCircumferential strain of the face; n represents the total number of the optical fiber strain measuring points; p represents the length of the tunnel section between two strain measuring points according to the optical fiber; lambda represents the azimuth coordinate on the length of the tunnel section between two strain measurement points according to the optical fiber; i is a calculation symbol in the summation operation and represents the ith measuring point in the total n measuring points;
(2) by using the theory of conjugate curved beams and the solution of the elastic mechanics ring to the uniformly distributed pressure, the convergence deformation of the section segment ring at any azimuth coordinate (p, lambda) caused by the axial force is calculated by the following formula, namely
In the formula (I), the compound is shown in the specification,
represents a convergent deformation at the azimuth coordinate (p, λ) caused by the axial force;
representing the strain monitoring data of the circumferential measuring point of the optical fiber;
(3) by using the conjugate beam theory, the convergence deformation at any point position of the tunnel section is calculated by the following formula, namely
In the formula, omegap,λThe convergence distortion at the position coordinates (p, λ) of the tunnel cross section is represented.
The third concrete implementation mode: in a specific embodiment, the tunnel section convergence deformation quantitative calculation method based on the high-density measurement point strain includes, in step two, a tunnel section additional horizontal convergence deformation calculation method under the squeezing effect:
(1) the additional horizontal convergence deformation of the section caused by the extrusion effect is influenced by the thickness of the duct piece, the radius of the duct piece, the inertia moment of the duct piece and the longitudinal bending curvature of the tunnel, the connection rigidity of the section is considered according to the action mechanism of the extrusion effect caused by the longitudinal bending, and the additional horizontal convergence deformation of the section under the extrusion effect is calculated by the following formula, namely
In the formula (I), the compound is shown in the specification,
representing the additional horizontal convergence deformation of the tunnel section under the extrusion effect; c represents the structural thickness of the tunnel segment; k represents the longitudinal bending curvature of the tunnel structure; i isoRepresenting the inertia moment of the cross section of the tunnel;
(2) when the horizontal convergence deformation occurs on the section of the tunnel, the resistance of the foundation can restrict the development of the deformation, and the additional horizontal convergence deformation after the reaction force of the foundation is considered is calculated by adopting the following formula, namely
In the formula (I), the compound is shown in the specification,
representing the additional horizontal convergence deformation of the tunnel section under the extrusion effect; k represents a ground coefficient.
(1) The horizontal convergence deformation directly generated under the action of the horizontal extrusion is calculated, and the foundation reaction force in the step (2) is generated by the horizontal extrusion in the step (1) and can also influence the horizontal convergence deformation under the action of the foundation reaction force.
The fourth concrete implementation mode: in a specific embodiment, a method for quantitatively calculating convergence deformation of a tunnel cross section based on strain at a high-density measuring point includes the following steps:
(1) when the shield tunnel structure is longitudinally bent and deformed, additional shearing force can be generated at two sides of the tunnel section due to uneven dislocation, the additional horizontal convergence deformation of the tunnel section caused by the shearing effect is calculated by the following formula,
in the formula (I), the compound is shown in the specification,
representing the additional horizontal convergence deformation of the tunnel section caused by the shearing effect; e represents the elastic modulus of the tunnel segment structure; qsRepresenting the shear load acting on the cross section of the tunnel due to longitudinal bending deformation;
(2) the additional horizontal convergence deformation of the shield tunnel section under the longitudinal bending effect is calculated by the following formula, namely
In the formula, omegaQRepresenting the additional horizontal convergence deformation of the shield tunnel section under the longitudinal bending effect;
(3) on the basis of obtaining the section convergence deformation by using the circumferential strain of the high-density measuring points, the horizontal convergence deformation of the section after the additional deformation caused by the extrusion effect and the shearing effect is considered to be calculated by the following formula, namely
ω＝ωp,1/2+ωQ (8)
In the formula, ωp,1/2The convergence deformation of the horizontal position of the section of the tunnel is represented and determined by a formula (3); and omega represents horizontal convergence deformation of the shield tunnel section.
The fifth concrete implementation mode is as follows: in a detailed embodiment, a tunnel section convergence deformation quantitative calculation method based on high-density measurement point strain includes, in step four, a tunnel multi-section horizontal convergence deformation calculation method that includes:
(1) the horizontal convergence deformation of the first ring section of the tunnel distributed by the annular full-length of the fully-distributed optical fiber is obtained by the formula (8) in the step three, and the horizontal convergence deformation of the first ring section at multiple moments is calculated by the following formula, namely
V＝g(ε) (9)
In the formula, g (-) represents the transformation process of shield tunnel section strain-horizontal convergence deformation; v represents a calculated value of horizontal convergence deformation of the lower section at m moments; epsilon represents strain monitoring data of circumferential measuring points of the lower section at m moments;
(2) taking the side wall strain of a first ring section of a tunnel with fully distributed optical fibers as an input layer and the horizontal convergence deformation as an output layer, training a strain-horizontal convergence deformation correlation model by using an artificial neural network, and calculating according to the following formula, namely
Youtput＝Z(Xinput) (10)
In the formula, Z (-) represents a neural network hidden layer relation model; xinputRepresenting a training set input layer of shield tunnel section side wall strain monitoring data; y isoutputRepresenting a shield tunnel section horizontal convergence deformation training set output layer;
(3) constructing a strain-horizontal convergence deformation correlation model by using an artificial neural network, and calculating the data of the input layer in the training set by the following formula, namely
Xinput＝[ε1,ε2,…,εt,…,εm],t＝1,2,…,m (11)
In the formula, m represents that the corresponding relation of m moments is selected as a training set; epsilontThe strain vector of the side wall of the shield tunnel section at the t-th moment is represented by the following formula,
in the formula (I), the compound is shown in the specification,
representing strain monitoring data of a rho-th measuring point of the lower section side wall at the t-th moment; rho represents the total number of strain measurement points of the side wall of the distributed optical fiber;
(4) constructing a strain-horizontal convergence deformation correlation model by using an artificial neural network, and calculating the data of the output layer in the training set by the following formula, namely
Youtput＝[ω1,ω2,…,ωt,…,ωm] (13)
In the formula, ωtRepresenting horizontal convergence deformation of the shield tunnel section at the t moment;
(5) taking side wall strain monitoring data of the fully distributed optical fiber circumferential 'n-shaped' distribution section of the shield tunnel as input, substituting the data into a strain-horizontal convergence correlation model trained by a formula (10), and calculating the horizontal convergence deformation of multiple sections according to the following formula, namely
Vt＝g(εt) (14)
In the formula, VtRepresenting the horizontal convergence deformation of the multi-section at the time t; epsilontAnd (4) representing the side wall strain monitoring data vector distributed in a shape like a Chinese character ji of the fully distributed optical fiber at t moment.
Example 1:
in this embodiment, a certain actual tunnel engineering is selected as an example, and the proposed tunnel section horizontal convergence deformation calculation method based on the high-density measuring point strain is verified. The resolution of distributed optical fiber measuring points applied in practical engineering is 0.2m, namely, each 20cm section has one strain measuring point, the distributed optical fiber distribution and monitoring range of SK6+ 400.216-SK 6+412.816 sections in a tunnel structure is selected, 157 effective structural strain measuring points are shared in the section, and distributed optical fiber monitoring data at a certain time are shown in figure 2.
In order to verify the effectiveness of the algorithm in the calculation of the horizontal convergence deformation of the tunnel section, stable monitoring data of the horizontal convergence deformation of the section of the total station instrument for seven days are selected, as shown in fig. 3 and 4. The comparison of the horizontal convergence deformation and the actual measurement deformation results obtained by the algorithm of the present invention with the existing algorithm is shown in fig. 5 and 6.
From the comparison between fig. 5 and fig. 6, it can be known that the horizontal convergence deformation of the tunnel section is basically in the same change trend along with time, the maximum relative errors of the horizontal convergence deformation of the left line and the right line in 7 days obtained by using the existing algorithm are 71.2% and 59.7%, while the relative errors of the horizontal convergence deformation of the left line and the right line in the algorithm provided by the present invention are about 13.6% and 15%, and it can be known that the calculation accuracy of the horizontal convergence deformation of the algorithm provided by the present invention is obviously improved compared with the existing algorithm, and the effectiveness of the algorithm provided by the present invention is verified.
The horizontal convergence deformation of the first full-length distribution section of the distributed optical fiber monitoring section calculated by the algorithm provided by the invention is shown in fig. 7 and 8. Taking the first ring section side wall strain as an input layer, taking the horizontal convergence deformation of the graph 7 and the graph 8 as an output layer, and training by using an artificial neural network, wherein the training iterative process is as shown in the graph 9. And (3) bringing the side wall strain of the fully distributed optical fiber n-shaped distribution section in the tunnel monitoring section into a strain-horizontal convergence correlation model to obtain the calculation results of the horizontal convergence deformation of the left line and the right line of the multi-section tunnel, as shown in fig. 10 and fig. 11.
From the comparison between fig. 10 and fig. 11, it can be seen that the horizontal convergence deformation obtained by engineering measurement has the same tendency as the multi-section horizontal convergence deformation obtained based on the strain-convergence correlation model, and the method provided by the present invention can reflect the convergence deformation degree in the monitoring interval, and verify the effectiveness of the calculation method for the multi-section horizontal convergence deformation.
Claims (5)
1. The tunnel section convergence deformation quantitative calculation method based on the high-density measuring point strain is characterized by comprising the following steps of: the method comprises the following steps:
the method comprises the following steps: the method comprises the steps of utilizing the circumferential strain monitoring data of high-density measuring points of the tunnel section acquired by a fully-distributed optical fiber sensor to construct a circular curved beam model of the tunnel section based on a conjugate curved beam theory, establishing a conversion model of the circumferential strain and the section convergence deformation of the high-density measuring points of the tunnel section, and calculating to obtain the convergence deformation of the tunnel section;
step two: constructing a basic structure of a tunnel section additional horizontal convergence deformation calculation force method, and calculating the additional horizontal convergence deformation of the section caused under the longitudinal bending effect of the tunnel; on the basis, the additional horizontal convergence deformation of the tunnel section under the extrusion effect caused by longitudinal bending is calculated by considering the constraint action of the actual tunnel section connection rigidity and the foundation reaction force;
step three: constructing a shear effect equation caused by longitudinal bending deformation, calculating the additional horizontal convergence deformation of the tunnel section under the shear effect, and establishing the calculation results of the horizontal convergence deformation of the first step and the second step in parallel to obtain the horizontal convergence deformation of the tunnel section considering the additional deformation influence;
step four: and (4) establishing a tunnel section strain-horizontal convergence deformation conversion model based on the artificial neural network by using the section horizontal convergence deformation calculated in the third step, and calculating to obtain the tunnel multi-section horizontal convergence deformation.
2. The tunnel section convergence deformation quantitative calculation method based on the high-density measuring point strain as claimed in claim 1, wherein: in the first step, the method for calculating the convergence deformation of the tunnel section comprises the following steps:
(1) acquiring strain monitoring data of high-density measuring points by using a fully-distributed sensing optical cable distributed along the circumferential direction through length of a tunnel section, and calculating the convergence deformation of the circumferential any position coordinates (p, lambda) (p is 1,2, …, n, and lambda is more than or equal to 0 and less than 1) of a tunnel section duct piece caused by bending moment according to the theory of conjugate curved beams by the following formula, namely
In the formula (I), the compound is shown in the specification,
representing the convergent deformation at the azimuthal coordinate (p, λ) caused by the bending moment; beta represents a corresponding central angle between two circumferential strain measuring points according to the optical fiber; r represents the radius of the segment of the tunnel section; h represents the height of the segment of the tunnel section;
respectively representing the circumferential strain of the section at the fiber distribution positions of the i point and the p point of the coordinate; n represents the total number of the optical fiber strain measuring points; p represents the length of the tunnel section between two strain measuring points according to the optical fiber; lambda represents the azimuth coordinate on the length of the tunnel section between two strain measurement points according to the optical fiber; i is a calculation symbol in the summation operation and represents the ith measuring point in the total n measuring points;
(2) by utilizing the conjugate curved beam theory and the solution of the elastic mechanical ring under the uniformly distributed pressure, the convergence deformation of the section duct piece at the annular position of any azimuth coordinate (p, lambda) caused by the axial force is calculated by the following formula, namely
In the formula (I), the compound is shown in the specification,
represents a convergent deformation at the azimuth coordinate (p, λ) caused by the axial force;
representing the strain monitoring data of the circumferential measuring point of the optical fiber;
(3) by using the conjugate beam theory, the convergence deformation at any point position of the tunnel section is calculated by the following formula, namely
In the formula, omegap,λThe convergence distortion at the position coordinates (p, λ) of the tunnel cross section is represented.
3. The tunnel section convergence deformation quantitative calculation method based on the high-density measuring point strain as claimed in claim 1, wherein: in the second step, the method for calculating the additional horizontal convergence deformation of the tunnel section under the extrusion effect comprises the following steps:
(1) according to the action mechanism of the extrusion effect caused by longitudinal bending, considering the connection rigidity of the cross section, the additional horizontal convergence deformation of the cross section under the extrusion effect is calculated by the following formula, namely
In the formula (I), the compound is shown in the specification,
representing the additional horizontal convergence deformation of the tunnel section under the extrusion effect; c represents the structural thickness of the tunnel segment; k represents the longitudinal bending curvature of the tunnel structure; I.C. AoRepresenting the inertia moment of the cross section of the tunnel;
(2) when the horizontal convergence deformation occurs on the section of the tunnel, the resistance of the foundation can restrict the development of the deformation, and the additional horizontal convergence deformation after the reaction force of the foundation is considered is calculated by adopting the following formula, namely
In the formula (I), the compound is shown in the specification,
representing the additional horizontal convergence deformation of the tunnel section under the extrusion effect; k represents a ground coefficient.
4. The tunnel section convergence deformation quantitative calculation method based on the high-density measuring point strain as claimed in claim 1, wherein: in the third step, the method for calculating the horizontal convergence deformation of the tunnel section considering the additional deformation influence comprises the following steps:
(1) when the shield tunnel structure is longitudinally bent and deformed, additional shearing force can be generated at two sides of the tunnel section due to uneven dislocation, the additional horizontal convergence deformation of the tunnel section caused by the shearing effect is calculated by the following formula,
in the formula (I), the compound is shown in the specification,
representing the additional horizontal convergence deformation of the tunnel section caused by the shearing effect; e represents the elastic modulus of the tunnel segment structure; qsIndicating due to longitudinal bending deformationActing on the shear load of the cross section of the tunnel;
(2) the additional horizontal convergence deformation of the shield tunnel section under the longitudinal bending effect is calculated by the following formula, namely
In the formula, ωQRepresenting the additional horizontal convergence deformation of the shield tunnel section under the longitudinal bending effect;
(3) on the basis of obtaining the section convergence deformation by using the high-density measuring point annular strain, the horizontal convergence deformation of the section after considering the additional deformation caused by the extrusion effect and the shearing effect is calculated by the following formula, namely
ω＝ωp,1/2+ωQ (8)
In the formula, omegap,1/2The convergence deformation at the horizontal position of the section of the tunnel is represented and determined by a formula (3); and omega represents horizontal convergence deformation of the shield tunnel section.
5. The tunnel section convergence deformation quantitative calculation method based on the high-density measuring point strain as claimed in claim 1, characterized in that: in the fourth step, the method for calculating the horizontal convergence deformation of the tunnel multi-section comprises the following steps:
(1) the horizontal convergence deformation of the first ring section of the tunnel distributed by the annular full length of the fully distributed optical fiber is obtained by the formula (8) in the step three, and the horizontal convergence deformation of the first ring section at a plurality of moments is calculated by the following formula, namely
V＝g(ε) (9)
In the formula, g (-) represents the transformation process of shield tunnel section strain-horizontal convergence deformation; v represents a calculated value of horizontal convergence deformation of the lower section at m moments; epsilon represents strain monitoring data of circumferential measuring points of the lower section at m moments;
(2) taking the side wall strain of a first ring section of a tunnel with fully distributed optical fibers as an input layer and the horizontal convergence deformation as an output layer, training a strain-horizontal convergence deformation correlation model by using an artificial neural network, and calculating according to the following formula, namely
Youtput＝Z(Xinput) (10)
In the formula, Z (-) represents a neural network hidden layer relation model; xinputRepresenting a training set input layer of shield tunnel section side wall strain monitoring data; y isoutputRepresenting the output layer of the shield tunnel section horizontal convergence deformation training set;
(3) constructing a strain-horizontal convergence deformation correlation model by using an artificial neural network, and calculating the data of the input layer in the training set by the following formula, namely
Xinput＝[ε1,ε2,…,εt,…,εm],t＝1,2,…,m (11)
In the formula, m represents that the corresponding relation of m moments is selected as a training set; epsilontThe strain vector of the side wall of the shield tunnel section at the t-th moment is represented, and the definition is shown as the following formula,
in the formula (I), the compound is shown in the specification,
representing strain monitoring data of a rho-th measuring point of the side wall of the lower section at the t moment; rho represents the total number of strain measuring points of the side wall of the distributed optical fiber;
(4) constructing a strain-horizontal convergence deformation correlation model by using an artificial neural network, and calculating the data of an output layer in a training set by the following formula, namely
Youtput＝[ω1,ω2,…,ωt,…,ωm] (13)
In the formula, omegatRepresenting horizontal convergence deformation of the shield tunnel section at the t moment;
(5) taking side wall strain monitoring data of the fully distributed optical fiber circumferential 'n-shaped' distribution section of the shield tunnel as input, substituting the data into a strain-horizontal convergence correlation model trained by a formula (10), and calculating the horizontal convergence deformation of multiple sections according to the following formula, namely
Vt＝g(εt) (14)
In the formula, VtShowing the horizontal convergence deformation of multiple sections at the time t; epsilontAnd (4) representing the side wall strain monitoring data vector distributed in a shape like a Chinese character ji of the fully distributed optical fiber at t moment.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210369310.5A CN114777661B (en) | 2022-04-08 | 2022-04-08 | Tunnel section convergence deformation quantitative calculation method based on high-density measuring point strain |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210369310.5A CN114777661B (en) | 2022-04-08 | 2022-04-08 | Tunnel section convergence deformation quantitative calculation method based on high-density measuring point strain |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114777661A true CN114777661A (en) | 2022-07-22 |
| CN114777661B CN114777661B (en) | 2022-10-21 |
Family
ID=82430108
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210369310.5A Active CN114777661B (en) | 2022-04-08 | 2022-04-08 | Tunnel section convergence deformation quantitative calculation method based on high-density measuring point strain |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114777661B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116776416A (en) * | 2023-05-11 | 2023-09-19 | 西南交通大学 | Method and system for calculating internal force of cross-activity fault tunnel by considering fracture zone width |
| CN117113728A (en) * | 2023-10-20 | 2023-11-24 | 湖南大学 | Shield tunnel convergence monitoring method and system based on bolt strain |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002156229A (en) * | 2000-11-17 | 2002-05-31 | Kajima Corp | Mobile displacement measuring method and device for structure |
| US20050128196A1 (en) * | 2003-10-08 | 2005-06-16 | Popescu Voicu S. | System and method for three dimensional modeling |
| CN102345460A (en) * | 2011-06-21 | 2012-02-08 | 北京交通大学 | Settlement control method for shallow large-section loess tunnel |
| WO2013135033A1 (en) * | 2012-03-16 | 2013-09-19 | 同济大学 | Tunnel deformation online monitoring system based on image analysis and application thereof |
| CN103591982A (en) * | 2013-08-21 | 2014-02-19 | 国家电网公司 | Method for monitoring electric power tunnel structure problem |
| CN106989685A (en) * | 2017-03-21 | 2017-07-28 | 山西省交通科学研究院 | A kind of duct pieces of shield tunnel overall deformation monitoring device and its monitoring method based on distribution type fiber-optic |
| CN108050952A (en) * | 2018-01-16 | 2018-05-18 | 陕西高速星展科技有限公司 | A kind of tunnel cross-section DEFORMATION MONITORING SYSTEM and monitoring method |
| CN108180856A (en) * | 2018-01-30 | 2018-06-19 | 中国地质大学（武汉） | A kind of tunnel deformation monitoring method, equipment and storage device based on laser data |
| CN207879337U (en) * | 2018-02-24 | 2018-09-18 | 黄河勘测规划设计有限公司 | The shield formula TBM of system is monitored with country rock convergent deformation |
| CN109724530A (en) * | 2019-01-28 | 2019-05-07 | 上海同岩土木工程科技股份有限公司 | A kind of tunnel convergence deformation on-line monitoring method based on image analysis |
| US20200064119A1 (en) * | 2018-08-22 | 2020-02-27 | Government Of The United States Of America, As Represented By The Secretary Of Commerce | Non-contact coordinate measuring machine using a noncontact metrology probe |
| US20200182615A1 (en) * | 2018-12-05 | 2020-06-11 | Chang'an University | Quick fixing device for measuring tunnel peripheral convergence and application method thereof |
| CN111609805A (en) * | 2020-04-23 | 2020-09-01 | 哈尔滨工业大学 | Tunnel structure state diagnosis method based on full-distribution strain measurement point section curvature |
| CN113094778A (en) * | 2021-03-11 | 2021-07-09 | 中国铁路兰州局集团有限公司兰州工程建设指挥部 | High-ground-stress interbed soft rock tunnel failure mechanism and construction control research method |
-
2022
- 2022-04-08 CN CN202210369310.5A patent/CN114777661B/en active Active
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002156229A (en) * | 2000-11-17 | 2002-05-31 | Kajima Corp | Mobile displacement measuring method and device for structure |
| US20050128196A1 (en) * | 2003-10-08 | 2005-06-16 | Popescu Voicu S. | System and method for three dimensional modeling |
| CN102345460A (en) * | 2011-06-21 | 2012-02-08 | 北京交通大学 | Settlement control method for shallow large-section loess tunnel |
| WO2013135033A1 (en) * | 2012-03-16 | 2013-09-19 | 同济大学 | Tunnel deformation online monitoring system based on image analysis and application thereof |
| CN103591982A (en) * | 2013-08-21 | 2014-02-19 | 国家电网公司 | Method for monitoring electric power tunnel structure problem |
| CN106989685A (en) * | 2017-03-21 | 2017-07-28 | 山西省交通科学研究院 | A kind of duct pieces of shield tunnel overall deformation monitoring device and its monitoring method based on distribution type fiber-optic |
| CN108050952A (en) * | 2018-01-16 | 2018-05-18 | 陕西高速星展科技有限公司 | A kind of tunnel cross-section DEFORMATION MONITORING SYSTEM and monitoring method |
| CN108180856A (en) * | 2018-01-30 | 2018-06-19 | 中国地质大学（武汉） | A kind of tunnel deformation monitoring method, equipment and storage device based on laser data |
| CN207879337U (en) * | 2018-02-24 | 2018-09-18 | 黄河勘测规划设计有限公司 | The shield formula TBM of system is monitored with country rock convergent deformation |
| US20200064119A1 (en) * | 2018-08-22 | 2020-02-27 | Government Of The United States Of America, As Represented By The Secretary Of Commerce | Non-contact coordinate measuring machine using a noncontact metrology probe |
| US20200182615A1 (en) * | 2018-12-05 | 2020-06-11 | Chang'an University | Quick fixing device for measuring tunnel peripheral convergence and application method thereof |
| CN109724530A (en) * | 2019-01-28 | 2019-05-07 | 上海同岩土木工程科技股份有限公司 | A kind of tunnel convergence deformation on-line monitoring method based on image analysis |
| CN111609805A (en) * | 2020-04-23 | 2020-09-01 | 哈尔滨工业大学 | Tunnel structure state diagnosis method based on full-distribution strain measurement point section curvature |
| CN113094778A (en) * | 2021-03-11 | 2021-07-09 | 中国铁路兰州局集团有限公司兰州工程建设指挥部 | High-ground-stress interbed soft rock tunnel failure mechanism and construction control research method |
Non-Patent Citations (3)
| Title |
|---|
| BAHMAN MEYGHANI等: "Progress in Thermomechanical Analysis of Friction Stir Welding", 《CHINESE JOURNAL OF MECHANICAL ENGINEERING》 * |
| 夏才初等: "基于GZZ强度准则考虑应变软化特性的深埋隧道弹塑性解", 《岩石力学与工程学报》 * |
| 沈圣等: "基于分布式光纤应变传感技术的盾构隧道横截面收敛变形监测方法", 《土木工程学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116776416A (en) * | 2023-05-11 | 2023-09-19 | 西南交通大学 | Method and system for calculating internal force of cross-activity fault tunnel by considering fracture zone width |
| CN116776416B (en) * | 2023-05-11 | 2023-12-29 | 西南交通大学 | Method and system for calculating internal force of cross-activity fault tunnel by considering fracture zone width |
| CN117113728A (en) * | 2023-10-20 | 2023-11-24 | 湖南大学 | Shield tunnel convergence monitoring method and system based on bolt strain |
| CN117113728B (en) * | 2023-10-20 | 2024-02-02 | 湖南大学 | Shield tunnel convergence monitoring method and system based on bolt strain |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114777661B (en) | 2022-10-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114777661B (en) | Tunnel section convergence deformation quantitative calculation method based on high-density measuring point strain | |
| CN110127078B (en) | Optical fiber monitoring method for strain-deflection-bending moment state of helicopter blade structure | |
| CN111609805B (en) | Tunnel structure state diagnosis method based on full-distribution strain measurement point section curvature | |
| CN105956216B (en) | Correction method for finite element model greatly across steel bridge based on uniform temperature response monitor value | |
| CN109632963A (en) | It is a kind of based on when invariant features signal building structural damage four-dimensional imaging method | |
| CN114548375B (en) | Cable-stayed bridge girder dynamic deflection monitoring method based on two-way long-short-term memory neural network | |
| Luo et al. | Data missing mechanism and missing data real-time processing methods in the construction monitoring of steel structures | |
| CN110162895A (en) | A kind of two stage high energy efficiency ship form optimization design method | |
| CN114858126B (en) | Tunnel longitudinal settlement deformation quantitative calculation method based on high-density measuring point strain | |
| CN110887448B (en) | Beam structure form reconstruction method based on optical fiber strain measurement | |
| CN116362083A (en) | Rail beam finite element model parameter joint correction method based on response surface | |
| CN111475940B (en) | Flexible baseline dynamic prediction method based on fiber bragg grating sensor and wing mode | |
| CN112989265A (en) | Pile body horizontal displacement calculation method based on optical fiber sensing | |
| CN114858079B (en) | Tunnel deformation monitoring method based on distributed multipoint strain and displacement conversion network | |
| CN114741923B (en) | Large-span bridge sling force identification method based on high-density measuring point strain | |
| KR101657425B1 (en) | Stress distribution estimation system for structure using of terrestrial laser scanning, thereof method and health monitoring system using it | |
| JP2001221616A (en) | Method for detecting strained position of structure using optical fiber strain distribution measuring instrument | |
| CN116428977B (en) | Main cable strand positioning and adjusting method, system and storage medium based on laser array | |
| CN102735471A (en) | Strain monitoring-based damaged cable recognition method in support generalized displacement and temperature change | |
| KR102332188B1 (en) | A method of determining the measurement location and location of a structure using a genetic algorithm and AI technology for precise strain-displacement prediction | |
| CN117571252A (en) | Vibration mode correction method for wind vibration response of gas-elastic model of power transmission tower | |
| CN109163651B (en) | Cantilever member disturbance degree measuring device and method based on strain | |
| CN116296164A (en) | Device and method for testing arch ring deflection of arch bridge based on trabecula strain and inclination angle | |
| Campos et al. | Uncertainty Propagation in the Estimation of Bending Stress in ACSR Conductors | |
| CN116252975A (en) | Large-scale space structure thermal stability test precision quantitative evaluation method |
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 |