CN114777661B - 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
- CN114777661B CN114777661B CN202210369310.5A CN202210369310A CN114777661B CN 114777661 B CN114777661 B CN 114777661B CN 202210369310 A CN202210369310 A CN 202210369310A CN 114777661 B CN114777661 B CN 114777661B
- Authority
- CN
- China
- Prior art keywords
- section
- convergence deformation
- tunnel
- deformation
- strain
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Operations Research (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Length Measuring Devices By Optical Means (AREA)
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 the high-density measuring points of the tunnel section, provides a calculation method of the additional horizontal convergence deformation amount of the tunnel section 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 the 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 tunnel section convergence deformation quantitative calculation method based on high-density measuring point strain.
Background
The development and the evolution of the horizontal convergence deformation of the section of the tunnel 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 cannot generally consider the additional horizontal convergence deformation caused by the longitudinal bending of the shield tunnel, so that a certain error exists in a 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: 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.
Compared with the prior art, the invention has the beneficial effects that: the invention constructs a conversion model of the high-density measuring point circumferential strain and the section convergence deformation of the tunnel section, provides a calculation method of the tunnel section additional horizontal convergence deformation amount considering the extrusion effect and the shearing effect, establishes a tunnel section strain-horizontal convergence deformation conversion model by utilizing the artificial neural network theory, and provides a calculation method of the tunnel multi-section horizontal convergence deformation. 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 results of the engineering of horizontal convergence and deformation of the right line of the tunnel cross section according to 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 comparison graph of the calculation results of the horizontal convergence deformation of the right line of the tunnel cross section according to 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 embodiment artificial neural network training iteration graph;
FIG. 10 is a comparison graph of the calculated 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 calculated result of the horizontal convergence deformation of the right line of the tunnel multi-section according to the embodiment with the actual measurement.
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 modifications or equivalent substitutions may be made to the technical solution of the present invention without departing from the spirit of the technical solution of the present invention, and the technical solution of the present invention is covered by 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: 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 section additional horizontal convergence deformation 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.
The invention provides a tunnel section horizontal convergence deformation calculation method based on high-density measuring point strain, aiming at the problem of how to calculate the horizontal convergence deformation of the 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 specific embodiment, a method for quantitatively calculating convergence deformation of a tunnel cross section based on strain of a high-density measuring point includes:
(1) Acquiring strain monitoring data of high-density measuring points by using a fully-distributed sensing optical cable distributed along the circumferential direction of a tunnel section, and calculating the convergence deformation of the circumferential random position coordinates (p, lambda) (p =1,2, \ 8230;, n; 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 arrangement 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 distance between two strain measurement points according to the optical fiberA tunnel segment length; 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 convergence deformation at the azimuth coordinates (p, λ) caused by the axial force;
representing the strain monitoring data of the circumferential measuring points of the optical fiber;
(3) By using the conjugate beam theory, the convergence deformation of the tunnel section at any point position is calculated by the following formula, namely
In the formula, ω p,λ Representing the convergent deformation at the coordinates (p, λ) of the location of the tunnel section.
The third concrete implementation mode: in a specific embodiment, the method for quantitatively calculating the convergence deformation of the tunnel section based on the strain of the high-density measuring points includes the following steps:
(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 is o Representing 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; q s Representing 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, ω Q Representing 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, omega p,1/2 The 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 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 distributed with fully distributed optical fibers as an input layer and the horizontal convergence deformation thereof as an output layer, training a strain-horizontal convergence deformation correlation model by utilizing an artificial neural network, and calculating by the following formula, namely
Y output ＝Z(X input ) (10)
In the formula, Z (-) represents a neural network hidden layer relation model; x input Representing a training set input layer of shield tunnel section side wall strain monitoring data; y is output Representing 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
X input ＝[ε 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; epsilon t The 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 side wall of the lower section at the t 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
Y output ＝[ω 1 ,ω 2 ,…,ω t ,…,ω m ] (13)
In the formula, ω t Indicates the t-th timeThe section of the lower shield tunnel is horizontally converged and deformed;
(5) Taking side wall strain monitoring data of a shield tunnel full-distributed optical fiber annular 'n-shaped' distribution section 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 by the following formula, namely
V t ＝g(ε t ) (14)
In the formula, V t Representing the horizontal convergence deformation of the multi-section at the time t; epsilon t And (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 has one strain measuring point, the distributed optical fiber arrangement and monitoring range of an SK6+ 400.216-SK 6+412.816 section 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 fig. 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 calculated by the existing algorithm and the algorithm provided by the present invention is shown in fig. 5 and fig. 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 invention can reflect the convergence deformation degree in the monitoring interval, thereby verifying 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 section additional horizontal convergence deformation 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 utilizing the section horizontal convergence deformation obtained by calculation in the step three, 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 of a tunnel section, and calculating the convergence deformation of the circumferential random position coordinates (p, lambda) (p =1,2, \ 8230;, n; 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 of the tunnel section at any point position is calculated by the following formula, namely
In the formula, omega p,λ Representing the convergent deformation at the coordinates (p, λ) of the location of the tunnel section.
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. A o Representing 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, characterized in that: 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; q s The shear load acting on the cross section of the tunnel caused by longitudinal bending deformation is shown;
(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, omega Q Representing 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, ω p,1/2 The 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, wherein: 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 distributed with fully distributed optical fibers as an input layer and the horizontal convergence deformation thereof as an output layer, training a strain-horizontal convergence deformation correlation model by utilizing an artificial neural network, and calculating by the following formula, namely
Y output ＝Z(X input ) (10)
In the formula, Z (-) represents a neural network hidden layer relation model; x input Representing a shield tunnel section side wall strain monitoring data training set input layer; y is output Representing 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
X input ＝[ε 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; epsilon t The 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 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
Y output ＝[ω 1 ,ω 2 ,…,ω t ,…,ω m ] (13)
In the formula, ω t Representing horizontal convergence deformation of the shield tunnel section at the t moment;
(5) Taking side wall strain monitoring data of a shield tunnel full-distributed optical fiber annular 'n-shaped' distribution section 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 by the following formula, namely
V t ＝g(ε t ) (14)
In the formula, V t Representing the horizontal convergence deformation of the multi-section at the time t; epsilon t And (3) 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 time.
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 CN114777661A (en) | 2022-07-22 |
| CN114777661B true 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) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| CN117113728B (en) * | 2023-10-20 | 2024-02-02 | 湖南大学 | Shield tunnel convergence monitoring method and system based on bolt strain |
Citations (11)
| 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 |
| 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 |
| 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 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7747067B2 (en) * | 2003-10-08 | 2010-06-29 | Purdue Research Foundation | System and method for three dimensional modeling |
| US11073373B2 (en) * | 2018-08-22 | 2021-07-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 |
| CN109405815B (en) * | 2018-12-05 | 2024-03-19 | 长安大学 | Quick fixing device and method for measuring peripheral convergence in tunnel |
-
2022
- 2022-04-08 CN CN202210369310.5A patent/CN114777661B/en active Active
Patent Citations (11)
| 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 |
| 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 |
| 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 |
|---|
| Progress in Thermomechanical Analysis of Friction Stir Welding;Bahman Meyghani等;《Chinese Journal of Mechanical Engineering》;20200215(第01期);全文 * |
| 基于GZZ强度准则考虑应变软化特性的深埋隧道弹塑性解;夏才初等;《岩石力学与工程学报》;20180825(第11期);全文 * |
| 基于分布式光纤应变传感技术的盾构隧道横截面收敛变形监测方法;沈圣等;《土木工程学报》;20130915(第09期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114777661A (en) | 2022-07-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114777661B (en) | Tunnel section convergence deformation quantitative calculation method based on high-density measuring point strain | |
| CN105956216B (en) | Correction method for finite element model greatly across steel bridge based on uniform temperature response monitor value | |
| CN111609805B (en) | Tunnel structure state diagnosis method based on full-distribution strain measurement point section curvature | |
| CN105956218B (en) | Steel bridge correction method for finite element model based on non-uniform temperature response monitor value | |
| Huang et al. | A novel two-stage structural damage identification method based on superposition of modal flexibility curvature and whale optimization algorithm | |
| CN102156796B (en) | Device and method for determining wind conditions using multiple wind resource grids | |
| CN110162895A (en) | A kind of two stage high energy efficiency ship form optimization design method | |
| CN113609749B (en) | Current calculation method applicable to multiple scenes and based on magnetic field signals | |
| CN114858126B (en) | Tunnel longitudinal settlement deformation quantitative calculation method based on high-density measuring point strain | |
| Simiu et al. | Advances in the design of high-rise structures by the wind tunnel procedure: Conceptual framework | |
| CN116362083A (en) | Rail beam finite element model parameter joint correction method based on response surface | |
| CN113188462B (en) | Fiber grating sensing rod for three-dimensional deformation measurement and measurement method thereof | |
| KR102435166B1 (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 | |
| CN113255035B (en) | Method for fusing monitoring data of small and medium-sized bridges | |
| CN114858079B (en) | Tunnel deformation monitoring method based on distributed multipoint strain and displacement conversion network | |
| CN111475940B (en) | Flexible baseline dynamic prediction method based on fiber bragg grating sensor and wing mode | |
| CN114722455A (en) | Three-dimensional engineering control network construction method combining total station and laser tracker | |
| CN112989265A (en) | Pile body horizontal displacement calculation method based on optical fiber sensing | |
| Peng et al. | Dynamic strain signal monitoring and calibration with neural network based on hierarchical orthogonal artificial bee colony | |
| CN111561907A (en) | Tower drum uneven settlement monitoring method based on plane dip angle measurement | |
| CN116070533B (en) | Neutron energy spectrum determination method | |
| CN114182770B (en) | Digital twinning-based marine booster station jacket foundation corrosion prediction method and system | |
| CN116186859B (en) | Wind vibration power calculation method and system for large-span cable roof surface component area | |
| CN116467878B (en) | Method, system, equipment and medium for detecting safety and health of assembled building | |
| CN117571252A (en) | Vibration mode correction method for wind vibration response of gas-elastic model of power transmission tower |
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 |