A method, system, electronic device and storage medium for estimating the bending moment of an enclosure structure

Abstract

A method, system, electronic device and storage medium for estimating the bending moment of an enclosure structure, wherein the method includes: establishing a coordinate system YOZ; determining a plurality of measuring points on the enclosure structure and obtaining the deformation measurement values at each measuring point; assuming that the fitting function of the horizontal displacement of the enclosure structure changing with depth is #imgabs0# calculating the horizontal displacement values at each measuring point according to the fitting function; based on the deformation measurement values and horizontal displacement values at each measuring point, obtaining the cumulative residual function #imgabs1# according to the constraint rules, obtaining a preset number of constraint functions; based on the cumulative residual function and the preset number of constraint functions, constructing a Lagrangian optimization function #imgabs2# according to the Lagrangian multiplication, taking the partial derivatives of F with respect to a i and λ k and setting the partial derivatives to zero, calculating the values of each coefficient of u(z); based on the values of each coefficient of u(z), obtaining the function of the bending moment of the enclosure structure changing with depth according to the relationship between the displacement and bending moment of the enclosure structure #imgabs3# The present invention overcomes the defects of the existing methods for estimating the bending moment of the enclosure structure and is reasonable and accurate.

Description

Technical Field

[0001] The present invention belongs to the technical field of foundation pits, and particularly relates to a method, a system, an electronic device and a storage medium for estimating the bending moment of a retaining structure. Background Art

[0002] Comprehensively considering the deep deformation and internal force of the retaining structure can better evaluate the safety state of the foundation pit. The deep deformation of the retaining structure is usually measured by an inclinometer. This monitoring technology is well-developed, the measurement accuracy meets the engineering requirements, and the deformation of the retaining structure can be densely obtained along the depth direction at low cost. The bending moment of the retaining structure is indirectly calculated by measuring the steel bar stress of the steel bar stress gauge buried in it. However, the survival rate of ordinary steel bar stress gauges during the embedding and use processes is relatively low. If high-performance steel bar stress gauges are used, the engineering monitoring cost will increase. In addition, if dense bending moment data needs to be obtained along the depth direction, the arrangement of the steel bar stress gauges along the depth direction is too dense, which will cause problems such as difficult element embedding and a significant increase in cost. Therefore, when considering the engineering monitoring plan, engineering technicians generally arrange fewer or no steel bar stress measurement points considering factors such as economic cost and practicability. There are also no mandatory requirements for the internal force monitoring of the retaining structure of conventional foundation pits in the current monitoring specifications, while deformation monitoring is generally mandatory. However, in order to better judge the safety state of the foundation pit, the bending moment is often estimated using the measured value of the deformation of the retaining structure.

[0003] Currently, the main common methods for estimating the bending moment of a retaining structure using the measured value of its deformation are as follows: First, the fixed circle method: The relationship between the horizontal displacement of the retaining structure and its deformation curvature is theoretically deduced. However, due to the errors generated during instrument measurement, the curvature calculated by this method oscillates and has a large error. Second, the method of fitting a five-point curve with a cubic function: This method requires selecting points for fitting in segments, and the fitting result is directly related to the selected points in segments, and its result has a certain randomness. In addition, dividing the deformation measurement points of the retaining structure into multiple segments is relatively complex and cumbersome to use in engineering. Third, the polynomial fitting method and the smoothing spline curve method: This method is simple to use, and the obtained bending moment of the retaining wall is relatively accurate in the middle of the retaining wall, but the end constraint conditions of the retaining structure are not considered, resulting in a large error in the bending moment estimation results of the two at the end part of the retaining structure, and the closer to the end part, the larger the error in the estimated bending moment.

[0004] In view of the above defects, there is an urgent need to provide a simple and easy-to-use, and more reasonable and accurate method for estimating the bending moment of a retaining structure. Summary of the Invention

[0005] The purpose of the present invention is to provide a method, a system, an electronic device and a storage medium for estimating the bending moment of a retaining structure, which overcome the defects of the existing methods for estimating the bending moment of a retaining structure and are reasonable and accurate.

[0006] The present invention is realized through the following technical solutions:

[0007] A method for estimating the bending moment of an enclosure structure includes the following steps:

[0008] Taking the centroid of the top section of the enclosure structure as the coordinate origin O, establish a coordinate system YOZ;

[0009] Determine multiple measuring points at different depths on the enclosure structure, and obtain the deformation measurement values of the enclosure structure at each measuring point;

[0010] Assume that the fitting function of the horizontal displacement of the enclosure structure along the depth in the coordinate system YOZ is Calculate the horizontal displacement values at each measuring point according to the fitting function, where z is the distance between the measuring point and the top section of the enclosure structure, u(z) is the horizontal displacement value of the enclosure structure at the depth z, a i is the coefficient of the fitting function, i = 0, 1,..., m;

[0011] Based on the deformation measurement values and horizontal displacement values at each measuring point, obtain the cumulative residual function where n is the number of measuring points, u(z i ) is the horizontal displacement value at the i-th measuring point, y i is the deformation measurement value at the i-th measuring point;

[0012] According to the constraint rules, obtain a preset number of constraint functions;

[0013] Based on the cumulative residual function and a preset number of constraint functions, according to the Lagrange multiplier method, introduce a preset number of Lagrange factors to construct a Lagrange optimization function where l is the preset number, λ k is the k-th Lagrange multiplier factor, is the k-th constraint function, where k = 1, 2,..., l;

[0014] Take the partial derivatives of F with respect to a i and λ k and set the partial derivatives to zero to calculate the values of each coefficient of u(z);

[0015] Based on the calculated values of each coefficient of u(z), obtain the function of the bending moment of the enclosure structure changing along the depth according to the relationship between the displacement and bending moment of the enclosure structure where E is the elastic modulus of the enclosure structure, I is the cross-sectional moment of inertia of the enclosure structure, is the value of the coefficient a i in u(z).

[0016] Furthermore, the steps of obtaining a preset number of constraint functions according to the constraint rules include:

[0017] Judge whether there is a support structure at the top of the retaining structure;

[0018] If there is no support structure at the top of the retaining structure, obtain the constraint function and

[0019] If there is a support structure at the top of the retaining structure, confirm whether the support structure is a hinged support or a fixed support;

[0020] If the support structure is a hinged support, obtain the constraint function

[0021] Based on the deformation measurement values at each measuring point, judge whether the bottom end of the retaining structure deforms;

[0022] If the bottom end of the retaining structure does not deform, obtain the constraint function where h is the height of the retaining structure;

[0023] If the bottom end of the retaining structure deforms, obtain the constraint function

[0024] Furthermore, the steps of obtaining a preset number of constraint functions according to the constraint rules include:

[0025] Select a preset number of target measuring points from multiple measuring points, and bury moment monitors at each target measuring point;

[0026] Based on the moment monitoring values measured by the moment monitors, obtain a preset number of constraint functions where z k is the distance between the i-th target measuring point and the top end of the maintenance structure, and b k is the moment monitoring value measured by the moment monitor at the i-th target measuring point.

[0027] Furthermore, the preset number is greater than or equal to 2.

[0028] The present invention also provides a retaining structure moment estimation system, including:

[0029] A building module, used to establish a coordinate system YOZ with the centroid of the top cross-section of the retaining structure as the coordinate origin O;

[0030] A determination module, used to determine multiple measuring points at different depths on the retaining structure and obtain the deformation measurement values of the retaining structure at each measuring point;

[0031] An assumption module, used to assume that the fitting function of the horizontal displacement of the retaining structure along the depth in the coordinate system YOZ is The horizontal displacement values at each measurement point are calculated according to the fitting function. Here, z is the distance between the measurement point and the top section of the retaining structure, u(z) is the horizontal displacement value of the retaining structure at the depth z, and a i is the coefficient of the fitting function, where i = 0, 1, …, m;

[0032] A first obtaining module, configured to obtain a cumulative residual function based on the deformation measurement values and horizontal displacement values at each measurement point where n is the number of measurement points, u(z i ) is the horizontal displacement value at the i-th measurement point, and y i is the deformation measurement value at the i-th measurement point;

[0033] An obtaining module, configured to obtain a preset number of constraint functions according to the constraint rules;

[0034] A second obtaining module, configured to introduce a preset number of Lagrange factors based on the cumulative residual function and the preset number of constraint functions, and construct a Lagrange optimization function according to the Lagrange multiplication where l is the preset number, and λ k is the k-th Lagrange multiplier factor, is the k-th constraint function, where k = 1, 2, …, l;

[0035] A calculation module, configured to take the partial derivatives of F with respect to a i and λ k and set the partial derivatives to zero, and calculate the values of each coefficient of u(z);

[0036] A third obtaining module, configured to obtain a function of the bending moment of the retaining structure varying with depth based on the calculated values of each coefficient of u(z where E is the elastic modulus of the retaining structure, I is the sectional moment of inertia of the retaining structure, is the value of the coefficient a i in u(z).

[0037] The present invention also discloses an electronic device, including a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, the steps of any one of the above methods are implemented.

[0038] The present invention also discloses a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of any one of the above methods are implemented.

[0039] Compared with the prior art, the beneficial effects of the present invention are as follows: simple and easy to use, more reasonable and accurate. Based on the measured deformation data of the retaining structure, the horizontal displacement fitting function and the constraint function of the retaining structure are organically combined, and the Lagrange multiplier factor is introduced to finally obtain the function of the bending moment of the retaining structure changing with depth, realizing the estimation of the bending moment at any depth of the retaining structure, so as to realize the automatic real-time estimation of the bending moment of the retaining structure. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] Figure 1 It is a flowchart of the steps of the method for estimating the bending moment of the retaining structure of the present invention;

[0041] Figure 2 It is a schematic diagram of the deformation monitoring of the retaining structure on the coordinate system YOZ;

[0042] Figure 3 It is a schematic diagram of the modules of the system for estimating the bending moment of the retaining structure of the present invention;

[0043] Figure 4 It is a schematic block diagram of the structure of an embodiment of the electronic device of the present invention;

[0044] Figure 5 It is a schematic diagram of the structure of an embodiment of the computer-readable storage medium of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0045] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. Usually, the components of the embodiments of the present invention described and illustrated herein can be arranged and designed in various different configurations.

[0046] Therefore, the following detailed description of the embodiments of the present invention provided in the drawings is not intended to limit the scope of the claimed present invention, but merely represents selected embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts fall within the scope of protection of the present invention.

[0047] It should be noted that: similar reference numerals and letters denote similar items in the following drawings. Therefore, once an item is defined in one drawing, it does not need to be further defined and explained in subsequent drawings. At the same time, in the description of the present invention, the terms "first", "second", etc. are only used for distinguishing descriptions and cannot be understood as indicating or implying relative importance.

[0048] It should be noted that in this article, relational terms such as first and second are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the term "comprising", "including" or any other variant thereof is intended to cover non-exclusive inclusion, so that a process, method, article or device comprising a series of elements not only includes those elements, but also includes other elements not expressly listed, or also includes elements inherent to such process, method, article or device. Without further limitation, an element defined by the statement "comprising an..." does not exclude the presence of additional identical elements in the process, method, article or device comprising the said element.

[0049] In the description of the present invention, it should be noted that the orientation or positional relationship indicated by the terms "upper", "lower", "inner", "outer", etc. is based on the orientation or positional relationship shown in the drawings, or the orientation or positional relationship in which the invention product is usually placed during use. It is only for the convenience of describing the present invention and simplifying the description, and does not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore should not be construed as a limitation of the present invention.

[0050] Please refer to Figure 1 and Figure 2 , Figure 1 is the flowchart of the steps of the method for estimating the bending moment of the enclosure structure of the present invention. Figure 2 is the schematic diagram of the deformation monitoring of the enclosure structure on the coordinate system YOZ. A method for estimating the bending moment of an enclosure structure includes the following steps:

[0051] S1. Taking the centroid of the top section of the enclosure structure as the coordinate origin O, establish a coordinate system YOZ;

[0052] S2. Determine a plurality of measuring points at different depths on the enclosure structure, and obtain the deformation measurement values of the enclosure structure at each measuring point;

[0053] S3. Assume that the fitting function of the horizontal displacement of the enclosure structure along the depth in the coordinate system YOZ is Calculate the horizontal displacement values at each measuring point according to the fitting function, where z is the distance between the measuring point and the top section of the enclosure structure, u(z) is the horizontal displacement value of the enclosure structure at the depth z, a i is the coefficient of the fitting function, i = 0, 1,..., m;

[0054] S4. Based on the deformation measurement values and horizontal displacement values at each measuring point, obtain the cumulative residual function where n is the number of measuring points, u(z i ) is the horizontal displacement value at the i-th measuring point, yi is the deformation measurement value at the i-th measurement point;

[0055] S5. According to the constraint rules, obtain a preset number of constraint functions;

[0056] S6. Based on the cumulative residual function and a preset number of constraint functions, according to the Lagrange multiplier method, introduce a preset number of Lagrange multipliers to construct a Lagrangian optimization function where l is the preset number, and λ k is the k-th Lagrange multiplier factor, is the k-th constraint function, where k = 1, 2,..., l;

[0057] S7. Take the partial derivatives of F with respect to a i and λ k , and set the partial derivatives to zero to calculate the values of each coefficient of u(z);

[0058] S8. Based on the values of each coefficient of u(z) calculated, obtain the function of the change of the bending moment of the retaining structure along the depth according to the relationship between the displacement and the bending moment of the retaining structure where E is the elastic modulus of the retaining structure, I is the cross-sectional moment of inertia of the retaining structure, is the value of the coefficient a i in u(z).

[0059] In the above step S1, taking the centroid of the top section of the retaining structure as the coordinate origin, the vertically downward direction as the positive direction of the z-axis, and the direction pointing to the inside of the foundation pit surrounded by the retaining structure as the positive direction of the y-axis, establish a plane rectangular coordinate system YOZ, as Figure 2 shown.

[0060] In the above step S2, select n measurement points on the retaining structure, and the measurement points are arranged at intervals along the height direction of the retaining structure so that the depths of the measurement points are different. The measurement point located at the topmost is the top of the retaining structure, and they are numbered in ascending order of depth from bottom to top. Then, obtain the deformation measurement values of the retaining structure at each measurement point through an inclinometer, so as to obtain a set of coordinates (z1, y1), (z2, y2), (z3, y3),..., (z n , y n ) in the plane rectangular coordinate system YOZ, where z1, z2, z3,..., z n are the distance values of the first measurement point, the second measurement point, the third measurement point,..., the n-th measurement point from the top of the retaining structure respectively, and y1, y2, y3,..., y n are the deformation measurement values at the first measurement point, the second measurement point, the third measurement point,..., the n-th measurement point respectively.

[0061] In the above step S3, according to mechanical theory, earth pressure is the fourth derivative of the displacement curve. Based on Terzaghi's research on the foundation pit of the Berlin subway and the recommended distribution form of earth pressure in the "Japanese Building Structure Foundation Design Code", it is considered that the distribution pattern of earth pressure is similar to a parabola. Therefore, it is assumed that the fitting function of the horizontal displacement of the retaining structure along the depth in the coordinate system YOZ is a polynomial of order m: Then the coefficient of the fitting function u(z) is a i , i = 0, 1, ……, m. For example, the fitting function is assumed to be a sixth-order polynomial: u(z) = a0 + a1z 1 + a2z 2 + a3z 3 + a4z 4 + a5z 5 + a6z 6 .

[0062] In the above step S4, for each measurement point, the difference function between the deformation measurement value and the horizontal displacement value at this measurement point is calculated, and n difference functions are obtained: L1 = u(z1) - y1, L2 = u(z2) - y2, ……, L n = u(z n ) - y n ; The sum of the squares of the n obtained difference functions is accumulated to obtain the cumulative residual function

[0063] In the above step S5, based on the measured data of the deformation of the retaining structure, the type of the retaining structure and the actual situation of the project, the constraint conditions of the retaining structure are determined. Then, based on the constraint conditions, a preset number of constraint functions are obtained according to the constraint rules k = 1, ……, l, is the distance of the constraint position from the top of the retaining structure. In the calculation process, the more constraint conditions are considered, the more constraint functions are obtained, and the more accurate the estimation is. Therefore, to ensure the accuracy of the calculation, at least 2 constraint conditions should be considered, so the preset number is greater than or equal to 2.

[0064] Furthermore, in step S5, the steps of obtaining a preset number of constraint functions according to the constraint rules include:

[0065] S51. Judge whether there is a support structure at the top of the retaining structure;

[0066] S52. If there is no support structure at the top of the retaining structure, then the constraint functions and

[0067] S53. If there is a support structure at the top of the retaining structure, confirm whether the support structure is a hinged support or a fixed support;

[0068] S54. If the support structure is a hinged support, the constraint function is obtained.

[0069] S55. Based on the deformation measurement values at each measuring point, determine whether the bottom end of the retaining structure has deformed.

[0070] S56. If the bottom end of the retaining structure has not deformed, the constraint function is obtained. where h is the height of the retaining structure.

[0071] S57. If the bottom end of the retaining structure has deformed, the constraint function is obtained.

[0072] In the above steps S51 to S57, when confirming the constraint conditions, the boundary constraints of the retaining structure can be considered, and the constraint conditions of the retaining structure are determined by using the deformation measurement values of the retaining structure measured by the inclinometer and the type of the retaining structure. Specifically, the form of the support structure at the top of the retaining structure directly determines the internal force boundary conditions at the top of the retaining structure. When there is no support structure at the top end of the retaining structure, it means that there is no constraint at the top of the retaining structure, and both the shear force and the bending moment at the top of the retaining structure are 0. And the depth at the top of the retaining structure is 0. Then, according to the relationship between the deformation and the internal force of the cross-section of the retaining structure in Table 1, and Therefore, two constraint functions can be obtained. and That is, two constraint functions are obtained. and When there is a support structure at the top of the retaining structure and the support structure is a hinged support, the bending moment at the hinged point is 0. Therefore, the bending moment at the top of the retaining structure is 0, and the shear force is not 0. According to Table 1, Thus, a constraint function can be obtained. That is, the constraint function is obtained. where the hinged support is a steel support, a cable anchor, etc. When there is a support structure at the top of the retaining structure and the support structure is a fixed support, both the shear force and the bending moment at the top of the retaining structure are not equal to 0. Therefore, the shear force and the bending moment at the top of the retaining structure cannot be determined and cannot be used as constraint conditions. Therefore, the constraint function cannot be obtained.

[0073] Table 1 Relationship between the deformation and the internal force of the cross-section of the retaining structure

[0074]

[0075] Whether the bottom end of the retaining structure is deformed determines the internal force boundary conditions at the bottom end of the retaining structure. According to the deformation measurement values of the retaining structure at each measuring point obtained by using an inclinometer in step S2, combined with the existing finite element numerical analysis method, the horizontal displacement of the retaining structure at different depths is calculated. Whether the bottom end of the retaining structure is deformed can be judged based on the horizontal displacement of the retaining structure at different depths. If there is no horizontal displacement at the bottom end of the retaining structure, it means that the bottom end of the retaining structure has not deformed, indicating that the constraint at the bottom end of the retaining structure is strong and the bending moment is large. Then, it can be considered that both the displacement and rotation angle at the bottom end of the retaining structure are 0, and the distance from the bottom end of the retaining structure to the top of the retaining structure is the height h of the retaining structure. Therefore, according to Table 1, Thus, a constraint function can be obtained That is, the constraint function is obtained If there is a horizontal displacement at the bottom end of the retaining structure, it means that the bottom end of the retaining structure has deformed, indicating that the constraint at the bottom end of the retaining structure is weak and the bending moment is small. At the same time, due to the end effect of the retaining structure, it can be considered that the bending moment at the bottom end of the retaining structure is 0. Therefore, according to Table 1, Thus, a constraint function can be obtained That is, the constraint function is obtained

[0076] Furthermore, to facilitate obtaining a preset number of constraint functions, in step S5, according to the constraint rules, the steps for obtaining a preset number of constraint functions include:

[0077] S58. Select a preset number of target measuring points from multiple measuring points and install bending moment monitors at each target measuring point;

[0078] S59. Based on the bending moment monitoring values measured by the bending moment monitors, obtain a preset number of constraint functions where z k is the distance between the i-th target measuring point and the top end of the retaining structure, and b k is the bending moment monitoring value measured by the bending moment monitor at the i-th target measuring point.

[0079] In the above steps S58 to S59, bending moment monitors are installed at each selected target measuring point to obtain the bending moment values of the retaining structure at each target measuring point through the bending moment monitors. The bending moment values measured by the bending moment monitors are actual values, so they can be used as constraint conditions. According to Table 1, Thus, Furthermore, steps S58 to S59 can be executed independently, or when the number of constraint functions obtained in steps S51 to S57 is less than the preset number, steps S58 to S59 are executed.

[0080] In the above step S6, the Lagrange multiplier method is a method for finding the extreme value of a multivariate function whose variables are restricted by one or more conditions. After obtaining a preset number of constraint functions, a preset number of Lagrange factors λ1, λ2, ……, λ l are introduced, and then combined with the cumulative residual function obtained in step S4 to construct a Lagrangian optimization function

[0081] In the above step S7, take the partial derivatives of F with respect to a i and λ k , and set the partial derivatives to zero to obtain n + 1 + l equations:

[0082]

[0083]

[0084]

[0085]

[0086]

[0087] Solve the n + 1 + l equations simultaneously to obtain the coefficients a i of u(z)

[0088] In the above step S8, substitute the value of the coefficient a i of u(z) obtained into u(z) to obtain According to the relationship between the displacement and bending moment of the retaining structure, obtain the bending moment function of the bending moment of the retaining structure varying with depth Substitute into to obtain where E is the elastic modulus of the retaining structure, and the elastic modulus E of the retaining structure material can be obtained through indoor mechanical tests. I is the moment of inertia of the cross-section of the retaining structure, which is related to the cross-section properties of the retaining pile. The retaining structure takes 1 linear meter to calculate the moment of inertia I of the cross-section. The specific calculation method can refer to the principles of theoretical mechanics and structural mechanics. This is the prior art and will not be elaborated here. Then calculate the flexural rigidity EI of the cross-section of the retaining structure according to the elastic modulus and moment of inertia of the cross-section. At this time, input any depth on the retaining structure into the bending moment function M(z) to obtain the estimated value of the bending moment at that depth, thus realizing the automatic real-time estimation of the bending moment of the retaining structure.

[0089] The following describes the application of the bending moment estimation method of the retaining structure of the present invention in an actual project:

[0090] A certain foundation pit project uses a diaphragm wall with a thickness of t and a depth of h as the retaining structure, and the top support of the retaining structure uses a steel pipe support.

[0091] Taking the centroid of the top section of the retaining structure as the origin of the coordinate system, a plane rectangular coordinate system YOZ is established. In this embodiment, it is assumed that the fitting function of the horizontal displacement of the retaining structure along the depth in the coordinate system YOZ is a sixth-order polynomial:

[0092] Then, multiple measuring points are determined from bottom to top in the retaining structure, and the distances between each measuring point and the top of the retaining structure are z1, z2, z3,..., z n , and the deformation measured values y1, y2, y3,..., y n of each measuring point are measured from bottom to top using an inclinometer. Substitute the depths z1, z2, z3,..., z n of each measuring point into the fitting function u(z) to obtain the horizontal displacement values at each measuring point.

[0093] Then calculate the difference function between the deformation measured value and the horizontal displacement value at each measuring point, and accumulate the squared values of the difference function to obtain the cumulative residual function

[0094] The support structure used for the retaining structure is a steel pipe support, which belongs to a hinged support. And it can be seen from the inclinometer data that there is deformation at the bottom of the retaining structure. Therefore, the boundary constraint conditions of this retaining structure are that the bending moments at both ends are 0. According to the boundary constraint conditions of the retaining structure, two constraint functions can be obtained from the constraint rules and

[0095] There are two boundary constraint conditions in this embodiment. According to the Lagrange multiplier method, two Lagrange multiplier factors λ1 and λ2 are introduced to construct the Lagrangian function Take the partial derivatives of F with respect to the coefficients a0, a1,..., a6 and the Lagrange multiplier factors λ1 and λ2, and set the partial derivatives to 0 to obtain nine equations, and construct the following system of equations:

[0096]

[0097] Solving the above system of equations can obtain the values of the coefficients a0, a1,..., a6 of the fitting function u(z) Substitute into the fitting function u(z) to obtain Substitute the fitting function u(z) into to obtain the bending moment function:

[0098]

[0099] By inputting the moment function M(z) at any depth of the retaining structure, the estimated moment value of the retaining structure at that depth can be obtained.

[0100] Please refer to Figure 3 , Figure 3 , which is a schematic diagram of the modules of the moment estimation system for the retaining structure of the present invention. The present invention also provides a moment estimation system for a retaining structure, including:

[0101] A establishing module 1, configured to establish a coordinate system YOZ with the centroid of the top cross-section of the retaining structure as the coordinate origin O;

[0102] A determining module 2, configured to determine a plurality of measuring points at different depths on the retaining structure and obtain the deformation measurement values of the retaining structure at each measuring point;

[0103] A hypothesizing module 3, configured to hypothesize that the fitting function of the horizontal displacement of the retaining structure along the depth in the coordinate system YOZ is The horizontal displacement values at each measuring point are calculated according to the fitting function, where z is the distance between the measuring point and the top cross-section of the retaining structure, u(z) is the horizontal displacement value of the retaining structure at the depth z, a i is the coefficient of the fitting function, and i = 0, 1,..., m;

[0104] A first obtaining module 4, configured to obtain a cumulative residual function based on the deformation measurement values and horizontal displacement values at each measuring point where n is the number of measuring points, u(z i ) is the horizontal displacement value at the i-th measuring point, and y i is the deformation measurement value at the i-th measuring point;

[0105] An obtaining module 5, configured to obtain a preset number of constraint functions according to the constraint rules;

[0106] A second obtaining module 6, configured to introduce a preset number of Lagrange factors based on the cumulative residual function and the preset number of constraint functions, and construct a Lagrange optimization function according to the Lagrange multiplication where l is the preset number, and λ k is the k-th Lagrange multiplier factor, is the k-th constraint function, where k = 1, 2,..., l;

[0107] A calculating module 7, configured to take the partial derivatives of F with respect to a i and λ k , and set the partial derivatives to zero, and calculate the values of each coefficient of u(z);

[0108] The third obtaining module 8 is configured to obtain a function of the moment of the retaining structure varying with depth based on the values of each coefficient of the calculated u(z). Where E is the elastic modulus of the retaining structure, and I is the cross-sectional moment of inertia of the retaining structure. is the coefficient a in u(z) i value.

[0109] The establishing module 1 takes the centroid of the top section of the retaining structure as the coordinate origin, the vertically downward direction as the positive direction of the z-axis, and the direction pointing to the inside of the foundation pit surrounded by the retaining structure as the positive direction of the y-axis to establish a plane rectangular coordinate system YOZ, as Figure 2 shown.

[0110] The determining module 2 selects n measuring points on the retaining structure. The measuring points are arranged at intervals along the height direction of the retaining structure so that the depths of the respective measuring points are different. The measuring point located at the uppermost position is the top of the retaining structure, and the measuring points are numbered in order of depth from bottom to top. Then, the deformation measurement values of the retaining structure at each measuring point are obtained by an inclinometer, so that a set of coordinates (z1, y1), (z2, y2), (z3, y3),..., (z n , y n ) in the plane rectangular coordinate system YOZ can be obtained. z1, z2, z3,..., z n are the distance values of the first measuring point, the second measuring point, the third measuring point,..., the nth measuring point from the top end of the retaining structure respectively, and y1, y2, y3,..., y n are the deformation measurement values at the first measuring point, the second measuring point, the third measuring point,..., the nth measuring point respectively.

[0111] The assuming module 3, according to the mechanical theory, the earth pressure is the fourth derivative of the displacement curve. According to the research on the foundation pit of the Berlin subway by Terzaghi and the recommended distribution form of the earth pressure in the "Japanese Building Structure Foundation Design Code", it is considered that the distribution form of the earth pressure is similar to a parabola. Therefore, it is assumed that the fitting function of the horizontal displacement of the retaining structure varying with depth in the coordinate system YOZ is a polynomial of order m: Then the coefficient of the fitting function u(z) is a i , i = 0, 1,..., m. For example, the fitting function is assumed to be a polynomial of order 6: u(z) = a0 + a1z 1 + a2z 2 + a3z 3 + a4z 4 + a5z 5 + a6z 6 .

[0112] The first obtaining module 4 calculates the difference function between the deformation measurement value and the horizontal displacement value at each measurement point for each measurement point, and obtains n difference functions L1 = u(z1) - y1, L2 = u(z2) - y2, ……, L n = u(z n ) - y n ; The sum of the squared values of the n obtained difference functions is accumulated to obtain the cumulative residual function

[0113] The obtaining module 5 determines the constraint conditions of the retaining structure based on the measured deformation data of the retaining structure, the type of the retaining structure, and the actual situation of the project, and thus obtains a preset number of constraint functions according to the constraint rules based on the constraint conditions is the distance from the constraint position to the top of the retaining structure. In the calculation process, the more constraint conditions are considered, the more constraint functions are obtained, and the more accurate the estimation is. Therefore, to ensure accurate calculation, at least 2 constraint conditions should be considered, so the preset number is greater than or equal to 2

[0114] Furthermore, the obtaining module 5 includes:

[0115] The first judgment unit is used to judge whether there is a support structure at the top of the retaining structure;

[0116] The first obtaining unit is used to obtain the constraint function if the first judgment unit judges that there is no support structure at the top of the retaining structure and

[0117] The confirmation unit is used to confirm whether the support structure is a hinged support or a fixed support if the first judgment unit judges that there is a support structure at the top of the retaining structure;

[0118] The second obtaining unit is used to obtain the constraint function if the confirmation unit judges that the support structure is a hinged support

[0119] The second judgment unit is used to judge whether the bottom end of the retaining structure deforms based on the deformation measurement values at each measurement point;

[0120] The third obtaining unit is used to obtain the constraint function if the second judgment unit judges that the bottom end of the retaining structure does not deform where h is the height of the retaining structure;

[0121] The fourth obtaining unit is used to obtain the constraint function if the second judgment unit judges that the bottom end of the retaining structure deforms

[0122] When confirming the constraint conditions, the boundary constraints of the retaining structure can be considered, and the constraint conditions of the retaining structure can be determined by using the measured deformation values of the retaining structure measured by the inclinometer and the type of the retaining structure. Specifically, the form of the support structure at the top of the retaining structure directly determines the internal force boundary conditions at the top of the retaining structure. When the first judgment unit judges that there is no support structure at the top of the retaining structure, it means that there is no constraint at the top of the retaining structure, and both the shear force and the bending moment at the top of the retaining structure are 0, while the depth at the top of the retaining structure is 0. Then, according to the relationship between the deformation and the internal force of the cross-section of the retaining structure in Table 1, and Therefore, the first obtaining unit can obtain two constraint functions and That is, two constraint functions are obtained and When the first judgment unit judges that there is a support structure at the top of the retaining structure and the confirmation unit judges that the support structure is a hinged support, the bending moment at the hinged place is 0. Therefore, the bending moment at the top of the retaining structure is 0, and the shear force is not 0. According to Table 1, Thus, the second obtaining unit can obtain a constraint function That is, the constraint function is obtained Among them, the hinged support is a steel support, a cable anchor, etc. When there is a support structure at the top of the retaining structure and the support structure is a fixed support, both the shear force and the bending moment at the top of the retaining structure are not equal to 0. Therefore, the shear force and the bending moment at the top of the retaining structure cannot be determined and cannot be used as constraint conditions, so no constraint function can be obtained.

[0123] Whether the bottom end of the retaining structure deforms determines the internal force boundary conditions at the bottom end of the retaining structure. The second judgment unit calculates the horizontal displacement of the retaining structure at different depths by using the measured deformation values of the retaining structure at each measuring point obtained by the inclinometer in step S2 and combining the existing finite element numerical analysis means. According to the horizontal displacement of the retaining structure at different depths, it can be judged whether the bottom end of the retaining structure deforms. If there is no horizontal displacement at the bottom end of the retaining structure, it means that the bottom end of the retaining structure has not deformed, indicating that the constraint at the bottom end of the retaining structure is strong and the bending moment is large. Then, it can be considered that both the displacement and the rotation angle at the bottom end of the retaining structure are 0, and the distance from the bottom end of the retaining structure to the top of the retaining structure is the height h of the retaining structure. Therefore, according to Table 1, Thus, the third obtaining unit can obtain a constraint function That is, the constraint function is obtained If the bottom end of the retaining structure has a horizontal displacement, it means that the bottom end of the retaining structure deforms, indicating that the constraint at the bottom end of the retaining structure is weak and the bending moment is small. At the same time, due to the end effect of the retaining structure, it can be considered that the bending moment at the bottom end of the retaining structure is 0. Therefore, according to Table 1, Thus, the fourth obtaining unit can obtain a constraint function The constraint function can be obtained

[0124] Further, to facilitate obtaining a preset number of constraint functions, the acquisition module 5 includes:

[0125] A burying unit, configured to select a preset number of target measuring points from multiple measuring points and bury a bending moment monitor at each target measuring point;

[0126] A fifth obtaining unit, configured to obtain a preset number of constraint functions based on the bending moment monitoring values measured by the bending moment monitor where z k is the distance between the i-th target measuring point and the top of the retaining structure, and b k is the bending moment monitoring value measured by the bending moment monitor at the i-th target measuring point.

[0127] The burying unit buries a bending moment monitor at each selected target measuring point to obtain the bending moment value of the retaining structure at each target measuring point through the bending moment monitor. The bending moment value measured by the bending moment monitor is the actual value, so it can be used as a constraint condition. According to Table 1, it can be obtained so that the fifth obtaining unit can obtain Further, the burying unit and the fifth obtaining unit can be executed separately, or when the number of constraint functions obtained after executing the first judgment unit to the fourth obtaining unit is less than the preset number, the burying unit and the fifth obtaining unit are executed.

[0128] Lagrange multiplication is a method for finding the extreme value of a multi-variable function where the variables are restricted by one or more conditions. After obtaining a preset number of constraint functions, the second obtaining module 6 introduces a preset number of Lagrange factors λ1, λ2, ……, λ l , and then combines with the cumulative residual function obtained by the first obtaining module 4 to construct a Lagrangian optimization function

[0129] The calculation module 7 calculates the partial derivatives of F with respect to a i and λ k , and sets the partial derivatives to zero to obtain n + 1 + l equations:

[0130]

[0131]

[0132]

[0133]

[0134]

[0135] Solving the system of \(n + 1 + l\) equations simultaneously to obtain the coefficients \(a\) of \(u(z)\) i value

[0136] The third obtaining module 8 will obtain the coefficients \(a\) of \(u(z)\) i value Substitute into \(u(z)\) to obtain According to the relationship between the displacement and bending moment of the retaining structure, obtain the bending moment function of the bending moment of the retaining structure varying with depth Substitute into to obtain Among them, \(E\) is the elastic modulus of the retaining structure, and the elastic modulus \(E\) of the retaining structure material can be obtained through indoor mechanical tests. \(I\) is the cross-sectional moment of inertia of the retaining structure, which is related to the cross-sectional properties of the retaining pile. The retaining structure takes 1 linear meter to calculate the cross-sectional moment of inertia \(I\). The specific calculation method can refer to the principles of theoretical mechanics and structural mechanics. This is the prior art and will not be elaborated here. Then, calculate the cross-sectional flexural stiffness \(EI\) of the retaining structure according to the elastic modulus and cross-sectional moment of inertia. At this time, input the bending moment function \(M(z)\) at any depth on the retaining structure to obtain the bending moment estimated value at that depth, so as to realize the automatic real-time estimation of the bending moment of the retaining structure.

[0137] Please refer to Figure 4 , Figure 4 The structural schematic block diagram of an embodiment of the electronic device of the present invention. An embodiment of the present invention also proposes an electronic device 1001, including a memory 1003 and a processor 1002. The memory 1003 stores a computer program 1004. When the processor 1002 executes the computer program 1004, it realizes the steps of any one of the above methods for enhancing the terminal identifier, including: S1. Taking the centroid of the top cross-section of the retaining structure as the coordinate origin \(O\), and establishing a coordinate system \(YOZ\); S2. Determining multiple measuring points at different depths on the retaining structure, and obtaining the deformation measurement values of the retaining structure at each measuring point; S3. Assuming in the coordinate system \(YOZ\) that the fitting function of the horizontal displacement of the retaining structure varying with depth is Calculating the horizontal displacement values at each measuring point according to the fitting function, where \(z\) is the distance between the measuring point and the top cross-section of the retaining structure, \(u(z)\) is the horizontal displacement value of the retaining structure at depth \(z\), and \(a\) i is the coefficient of the fitting function, \(i = 0, 1,\cdots,m\); S4. Based on the deformation measurement values and horizontal displacement values at each measuring point, obtain the cumulative residual function where \(n\) is the number of measuring points, \(u(z i ) is the horizontal displacement value at the \(i\)-th measuring point, and \(y iis the deformation measurement value at the i-th measurement point; S5. According to the constraint rules, obtain a preset number of constraint functions; S6. Based on the cumulative residual function and the preset number of constraint functions, according to the Lagrange multiplication, introduce a preset number of Lagrange factors, and construct a Lagrange optimization function where l is the preset number, and λ k is the k-th Lagrange multiplier factor, is the k-th constraint function, where k = 1, 2,..., l; S7. Take the partial derivatives of F with respect to a i and λ k and set the partial derivatives to zero, and calculate the values of each coefficient of u(z); S8. Based on the values of each coefficient of u(z) calculated, obtain the function of the bending moment of the retaining structure varying with depth according to the relationship between the displacement and bending moment of the retaining structure where E is the elastic modulus of the retaining structure, I is the sectional moment of inertia of the retaining structure, is the value of the coefficient a i in u(z).

[0138] Please refer to Figure 5 . Figure 5 is the structural schematic block diagram of an embodiment of the computer-readable storage medium of the present invention. An embodiment of the present invention also provides a computer-readable storage medium 2001, on which a computer program 1004 is stored. When the computer program 1004 is executed by a processor 1002, the steps of the method for enhancing the terminal identifier described in any one of the above are implemented, including: S1. Taking the centroid of the top cross-section of the retaining structure as the coordinate origin O, establish a coordinate system YOZ; S2. Determine a plurality of measurement points at different depths on the retaining structure, and obtain the deformation measurement values of the retaining structure at each measurement point; S3. Assume that the fitting function of the horizontal displacement of the retaining structure varying with depth in the coordinate system YOZ is Calculate the horizontal displacement values at each measurement point according to the fitting function, where z is the distance between the measurement point and the top cross-section of the retaining structure, u(z) is the horizontal displacement value of the retaining structure at the depth z, and a i is the coefficient of the fitting function, i = 0, 1,..., m; S4. Based on the deformation measurement values and horizontal displacement values at each measurement point, obtain a cumulative residual function where n is the number of measurement points, u(z i ) is the horizontal displacement value at the i-th measurement point, and y i is the deformation measurement value at the i-th measurement point; S5. According to the constraint rules, obtain a preset number of constraint functions; S6. Based on the cumulative residual function and the preset number of constraint functions, according to the Lagrange multiplication, introduce a preset number of Lagrange factors, and construct a Lagrange optimization function where l is the preset number, and λ k is the k-th Lagrange multiplier factor, is the k-th constraint function, where k = 1, 2, ……, l; S7. Take the partial derivatives of F with respect to a i and λ k and set the partial derivatives to zero to calculate the values of each coefficient of u(z); S8. Based on the calculated values of each coefficient of u(z), obtain the function of the variation of the bending moment of the retaining structure along the depth according to the relationship between the displacement and the bending moment of the retaining structure where E is the elastic modulus of the retaining structure, I is the sectional moment of inertia of the retaining structure, is the value of the coefficient a i in u(z).

[0139] Compared with the prior art, the beneficial effects of the present invention are: simple to use and more reasonable and accurate. Based on the measured deformation data of the retaining structure, the horizontal displacement fitting function of the retaining structure and the constraint function are organically combined, and the Lagrange multiplier factor is introduced to finally obtain the function of the variation of the bending moment of the retaining structure along the depth, realizing the estimation of the bending moment at any depth of the retaining structure, so as to realize the automatic real-time estimation of the bending moment of the retaining structure.

[0140] The above is only a preferred embodiment of the present invention, and does not impose any form of limitation on the present invention. Therefore, any simple modification, equivalent change and modification made to the above embodiment according to the technical essence of the present invention without departing from the technical solution of the present invention still fall within the scope of the technical solution of the present invention.

Claims

1. A method for estimating the bending moment of an enclosure structure, characterized in that, It includes the following steps: Taking the centroid of the top cross-section of the retaining structure as the coordinate origin O, establish a coordinate system YOZ; Determine multiple measuring points at different depths on the retaining structure, and obtain the deformation measurement values of the retaining structure at each measuring point; In the coordinate system YOZ, assume that the fitting function of the horizontal displacement of the enclosure structure varying with depth is The horizontal displacement values at each measuring point are calculated according to the fitting function, where z is the distance between the measuring point and the top section of the enclosure structure, u(z) is the horizontal displacement value of the enclosure structure at depth z, a i is the coefficient of the fitting function, and i = 0, 1, ……, m; Based on the deformation measurement values and horizontal displacement values at each of the measurement points, a cumulative residual function is obtained where n is the number of the measurement points, u(z i ) is the horizontal displacement value at the i-th measurement point, and y i is the deformation measurement value at the i-th measurement point; According to the constraint rules, construct a preset number of constraint functions; Based on the cumulative residual function and a preset number of constraint functions, according to the Lagrange multiplier method, introduce the preset number of Lagrange factors to construct a Lagrangian optimization function where l is the preset number, and λ k is the k-th Lagrange multiplier factor, is the k-th constraint function, where k = 1,..., l; Find a with respect to F i and λ k Take the partial derivatives and set them to zero to calculate the values of each coefficient of u(z); Based on the values of each coefficient of the calculated u(z), a function of the variation of the bending moment of the retaining structure along the depth is obtained according to the relationship between the displacement and the bending moment of the retaining structure where E is the elastic modulus of the retaining structure and I is the sectional moment of inertia of the retaining structure, is the coefficient a in u(z) i value.

2. The method for estimating the bending moment of the enclosure structure according to claim 1, wherein The step of obtaining a preset number of constraint functions according to the constraint rules includes: Judge whether there is a support structure at the top of the retaining structure; If there is no support structure at the top of the enclosure structure, the constraint function is obtained and If there is a support structure at the top of the retaining structure, confirm whether the support structure is a hinged support or a fixed support; If the support structure is a hinged support, the constraint function is obtained Based on the deformation measurement values at each measuring point, judge whether the bottom end of the retaining structure deforms; If the bottom end of the retaining structure does not deform, the constraint function is obtained where h is the height of the retaining structure; If the bottom end of the enclosure structure is deformed, the constraint function is obtained 3. The method for estimating the bending moment of the enclosure structure according to claim 1, wherein The step of obtaining a preset number of constraint functions according to the constraint rules includes: Select a preset number of target measuring points from the multiple measuring points, and bury moment monitors on each of the target measuring points; Based on the bending moment monitoring values measured by the bending moment monitor, the preset number of constraint functions are obtained where z k is the distance between the i-th target measuring point and the top of the retaining structure, and b k is the bending moment monitoring value measured by the bending moment monitor at the i-th target measuring point.

4. The method for estimating the bending moment of the enclosure structure according to claim 1, wherein The preset number is greater than or equal to 2.

5. A bending moment estimation system for an enclosure structure, characterized in that It includes: A establishing module, which is used to take the centroid of the top cross-section of the retaining structure as the coordinate origin O and establish a coordinate system YOZ; A determining module, which is used to determine multiple measuring points at different depths on the retaining structure and obtain the deformation measurement values of the retaining structure at each measuring point; A hypothesis module, which is used to assume that the fitting function of the horizontal displacement of the enclosure structure varying with depth under the coordinate system YOZ is The horizontal displacement values at each measuring point are calculated according to the fitting function, where z is the distance between the measuring point and the top section of the enclosure structure, u(z) is the horizontal displacement value of the enclosure structure at the depth z, and a i is the coefficient of the fitting function, and i = 0, 1,..., m; The first obtaining module is configured to obtain a cumulative residual function based on the deformation measurement values and the horizontal displacement values at each of the measurement points where n is the number of the measurement points, u(z i ) is the horizontal displacement value at the i-th measurement point, and y i is the deformation measurement value at the i-th measurement point; An obtaining module, which is used to obtain a preset number of constraint functions according to the constraint rules; A second obtaining module, configured to introduce the preset number of Lagrange factors according to the Lagrange multiplication based on the cumulative residual function and the preset number of constraint functions, and construct a Lagrange optimization function where l is the preset number, and λ k is the k-th Lagrange multiplier factor, is the k-th constraint function, where k = 1, 2,..., l; A calculation module for calculating the partial derivatives of F with respect to a i and λ k and setting the partial derivatives to zero to calculate the values of each coefficient of u(z); A third obtaining module, configured to obtain a function of the bending moment of the enclosure structure varying along the depth based on the value of each coefficient of the calculated u(z). where E is the elastic modulus of the enclosure structure, and I is the cross-sectional moment of inertia of the enclosure structure is the coefficient a in u(z) i value.

6. An electronic device, comprising a memory and a processor, the memory storing a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1-4.

7. A computer-readable storage medium storing a computer program thereon, characterized in that, When the computer program is executed by the processor, it implements the steps of the method according to any one of claims 1-4.

Images