A method for embedding and fixing a stiff column of an underground facility integrated with an aerial rail

Abstract

The application discloses a kind of air rail ground integration underground facilities stiff column embedded construction methods, it is related to underground deep foundation pit engineering technical field, including the following steps: obtaining the surface point cloud data of stiff column, and remapping generates stiff column feature point cloud;Stiff column is suspended and implanted pile hole and grouting;Stiff column stress distribution field is constructed, and initial stress point cloud is generated in combination with stiff column feature point cloud;Based on grouting pressure field, variation iteration corrects initial stress point cloud, and obtains stiff column stress point cloud, the variation iteration includes eliminating the fluid-solid coupling physical residual term based on energy gradient construction;Stiff column feature point cloud and stiff column stress point cloud are compared, and the pose deviation of stiff column is obtained, and according to the deviation difference, multi-channel grouting rate is regulated and controlled, until stiff column embedding is embedded.The application is used to solve the rectification optimization problem of stiff column embedding process.

Description

Technical Field

[0001] This invention relates to the field of underground deep foundation pit engineering technology, and more specifically, to a method for reinforcing the stiffening columns of an integrated underground facility with air-rail-ground structure. Background Technology

[0002] With the deep integration of urban rail transit and underground space development, integrated air-rail-ground projects have become the mainstream model of modern urban construction. In such projects, in order to effectively block the upward transmission of vibration and secondary noise generated by subway train operation, the superstructure usually passes through the subway station level or depot through deep stiffened steel columns, and high-precision vibration isolation supports are set at the top as structural transition nodes. This places extremely high installation requirements on the planar position and axial torsional accuracy of the interface at the top of the stiffened column.

[0003] Construction monitoring technology for such deep-stiffened columns mainly relies on optical equipment such as plumb bobs and total stations to measure the exposed parts at the borehole opening, or on distributed fiber optic sensing technology attached to the surface of the steel column. For fiber optic monitoring, the existing common processing method is to obtain strain data based on Rayleigh scattering or Brillouin scattering principles, and then invert the column morphology through tangent recursion or curvature integral algorithms using purely geometric paths.

[0004] In terms of construction control technology, existing correction operations mainly rely on mechanical adjustment frames or hydraulic correction platforms installed on the ground at the pile hole opening. These devices clamp the top of the steel column, apply horizontal thrust or adjust the position of the suspension point to achieve passive mechanical adjustment of the verticality and center position of the steel column, maintaining the posture of the steel column until the concrete hardens and forms.

[0005] The existing technology has at least the following problems: Existing technologies mostly rely on curvature integral algorithms based on pure geometric paths. Due to the lack of physical and mechanical boundary conditions in the calculation process, they cannot truly reflect the downhole attitude of the stiffening column. Furthermore, in existing construction methods, the stiffening column is easily subject to the random effects of fluid pressure during mud injection, resulting in passive deviations and making it difficult to achieve accurate embedding at full depth.

[0006] To address the above problems, this invention proposes a solution. Summary of the Invention

[0007] To overcome the aforementioned deficiencies of the prior art, embodiments of the present invention provide a method for the construction of stiffened column embedding in an integrated air-rail-ground underground facility. By constructing a feature point cloud of the stiffened column, fusing strain data, and using fluid-structure interaction physical residual terms as constraints for variational iterative correction, a stress point cloud of the stiffened column is obtained. The pose deviation of the generated point cloud is used to control the multi-channel grouting rate, thereby solving the problem of correction and optimization in the stiffened column embedding process.

[0008] To achieve the above objectives, the present invention provides the following technical solution: A method for embedding stiffening columns in an integrated air-rail-ground underground facility includes the following steps: acquiring surface point cloud data of the stiffening column and remapping it to generate a characteristic point cloud of the stiffening column; suspending the stiffening column and inserting it into the pile hole and grouting it; constructing a stress distribution field of the stiffening column and generating an initial stress point cloud by combining it with the characteristic point cloud of the stiffening column; using variational iteration to correct the initial stress point cloud based on the grouting pressure field to obtain the stress point cloud of the stiffening column, wherein the variational iteration includes eliminating fluid-structure interaction physical residual terms constructed based on energy gradients; comparing the characteristic point cloud of the stiffening column and the stress point cloud of the stiffening column to obtain the posture deviation of the stiffening column, and adjusting the multi-channel grouting rate according to the deviation difference until the stiffening column is embedded.

[0009] In a preferred embodiment, the step of acquiring the surface point cloud data of the stiffening column and remapping it to generate the feature point cloud of the stiffening column includes: acquiring the surface point cloud data of the stiffening column and remapping the surface point cloud data using a preset local coordinate system of the column; dividing the remapped surface point cloud data into equidistant sections along the axial direction to obtain a section point cloud; and constructing a feature vector using the section center and origin of the section point cloud to generate the feature point cloud of the stiffening column.

[0010] In a preferred embodiment, the step of constructing the stress distribution field of the stiffening column and generating an initial stress point cloud by combining the characteristic point cloud of the stiffening column includes: acquiring strain data of the stiffening column and fitting it to generate a linear strain distribution field for each section of the stiffening column; based on the elastic modulus of the stiffening column material, converting the strain values ​​in the linear strain distribution field into stress values ​​to construct the stress distribution field of the stiffening column; and substituting the characteristic point cloud of the stiffening column into the stress distribution field of the stiffening column for numerical calculation to generate the initial stress point cloud.

[0011] In a preferred embodiment, acquiring stiffener strain data and fitting it to generate a linear strain distribution field for each section of the stiffener includes: acquiring stiffener strain data through distributed fiber optic sensors and associating the strain data with the spatial coordinates of each sensor's section; performing least-squares plane fitting on the spatial coordinates and stiffener strain data based on the plane section assumption to construct a strain distribution plane equation for the current section; analyzing the strain distribution plane equation to extract the axial reference strain and bending curvature of the section and generate a linear strain distribution field for each section of the stiffener.

[0012] In a preferred embodiment, the variational iterative correction of the initial stress point cloud specifically involves: mapping the initial stress point cloud to a state vector; constructing the total energy functional of the fluid-structure interaction system, and combining it with the state vector to obtain the fluid-structure interaction physical residual term; performing variational iteration on the state vector along the opposite direction of the gradient with the goal of eliminating the residual term, to obtain a converged state vector; and performing inverse mapping based on the converged state vector to obtain the corrected stiffener stress point cloud.

[0013] In a preferred embodiment, the process of constructing the total energy functional of the fluid-structure interaction system, combined with the state vector, to obtain the fluid-structure interaction physical residual term includes: performing a volume integral operation on the volume domain of the stiffening column based on the initial stress point cloud to obtain the elastic strain energy; acquiring the fluid pressure distribution of the grouting environment and performing an integral operation on the side surface of the stiffening column to obtain the fluid pressure potential energy; constructing a total energy functional that includes the elastic strain energy and the fluid pressure potential energy, and calculating the energy gradient of the total energy functional with respect to the state vector; constructing a global gradient vector based on the energy gradient, and using the global gradient vector as the fluid-structure interaction physical residual term.

[0014] In a preferred embodiment, comparing the feature point cloud and stress point cloud of the stiffening column to obtain the pose deviation of the stiffening column specifically involves: extracting the center coordinates and axial angle of the feature point cloud of the stiffening column to determine the target pose of the stiffening column; extracting the center coordinates and axial angle of the stress point cloud of the stiffening column to determine the real-time pose of the stiffening column; calculating the spatial deviation of the real-time pose relative to the target pose; and decomposing the spatial deviation along the local coordinate system of the column to obtain the pose deviation of the stiffening column.

[0015] In a preferred embodiment, the step of differentially controlling the multi-channel grouting rate based on the deviation includes: calculating the target correction force required for correction based on the stiffener's positional deviation; constructing a kinematic Jacobian matrix based on the coupling relationship between the multi-channel pumping rate and the fluid pressure on the stiffener's side surface; and differentially controlling the pumping rate required for each channel based on the kinematic Jacobian matrix according to the target correction force.

[0016] In a preferred embodiment, the step of differentially adjusting the pumping rate of each channel according to the target correction force specifically involves: calculating the generalized inverse matrix of the kinematic Jacobian matrix; using the generalized inverse matrix to perform a linear mapping solution on the target correction force to obtain the pumping rate adjustment vector corresponding to each channel; and superimposing the pumping rate adjustment vector onto the current pumping rate of each channel to perform asymmetric control of the fluid pressure distribution in the grouting environment.

[0017] The technical effects and advantages of the present invention regarding the construction method for reinforcing stiffened columns in integrated underground facilities with air-rail-ground systems are as follows: This invention constructs a stiffening column feature point cloud by remapping the surface point cloud data obtained from scanning, making the stiffening column feature point cloud more accurately reflect the geometric shape of the stiffening column before grouting. The stiffening column feature point cloud is fused with strain data to generate an initial stress point cloud. The initial stress point cloud is then iteratively corrected using fluid-structure interaction physical residual terms as constraints to obtain the stiffening column stress point cloud. By analyzing the mechanical field of the grouting fluid and the stiffening column solid, the physical changes under grouting conditions are considered. The stiffening column stress point cloud accurately characterizes the physical shape of the stiffening column during the actual grouting process. The stiffening column pose deviation is generated using the stiffening column feature point cloud and the stiffening column stress point cloud. The pose deviation is then used to control the multi-channel grouting rate, enhancing the targeting and accuracy of the control and effectively solving the problem of correction and optimization during the stiffening column embedding process. Attached Figure Description

[0018] Figure 1 This is a schematic diagram of a construction method for embedding stiffening columns in an integrated underground facility with air-rail-ground connection, provided as an embodiment of the present invention. Detailed Implementation

[0019] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of the present invention.

[0020] Example 1, Figure 1 This invention discloses a method for reinforcing stiffened columns in an integrated underground facility with air-rail-ground connection, comprising the following steps: S1, acquire the surface point cloud data of the stiffening column, and remap to generate the feature point cloud of the stiffening column; S2, the stiffening column is suspended and inserted into the pile hole and grout is injected; S3, construct the stress distribution field of the stiffened column, and generate the initial stress point cloud by combining the characteristic point cloud of the stiffened column; S4. Based on the grouting pressure field, the initial stress point cloud is modified by variational iteration to obtain the stress point cloud of the stiffened column. The variational iteration includes eliminating the fluid-structure interaction physical residual term constructed based on the energy gradient. S5. Compare the feature point cloud and stress point cloud of the stiffening column to obtain the posture deviation of the stiffening column, and adjust the multi-channel grouting rate according to the deviation difference until the stiffening column is embedded.

[0021] In this embodiment, S1 includes: S11, acquire the surface point cloud data of the stiffened column, and remap the surface point cloud data using the preset column local coordinate system; S12, obtain the cross-sectional point cloud by dividing and remapping the surface point cloud data along the axial equidistant section; S13: Construct feature vectors using the cross-sectional center and origin of the cross-sectional point cloud to generate the stiffened column feature point cloud.

[0022] It should be noted that a vibration isolation flange is installed at the top of the stiffening column.

[0023] In this embodiment, the method for constructing the preset local coordinate system of the cylinder in step S11 includes: Feature point cloud of vibration isolation flange is extracted from surface point cloud data of stiffened column; The origin is the geometric center of the feature point cloud of the vibration isolation flange; The least squares method is used to fit the plane normal vector of the feature point cloud of the vibration isolation flange as the Z-axis; Extract the surface point cloud data of the stiffened column, and take the characteristic direction of the top plane of the stiffened column as the X-axis; The Y-axis is determined using the right-hand rule, and a local coordinate system for the cylinder is established.

[0024] In S13, a feature vector is constructed using the cross-sectional center and the origin of the cross-sectional point cloud, including the lateral displacement deviation of the cross-sectional center relative to the origin and the axial torsion angle formed with the Z-axis.

[0025] It should be noted that in this embodiment, equal-interval slicing is performed along the axial direction in a normalized local coordinate system to generate a cross-sectional point cloud for accurately extracting lateral displacement at each depth.

[0026] In this embodiment, S3 includes: S31, Obtain strain data of stiffened column and fit to generate linear strain distribution field of each section of stiffened column; S32, based on the elastic modulus of the stiffened column material, converts the strain values ​​in the linear strain distribution field into stress values ​​to construct the stiffened column stress distribution field; S33, substitute the characteristic point cloud of the stiffening column into the stress distribution field of the stiffening column for numerical solution to generate the initial stress point cloud.

[0027] It should be noted that in step S32, for the stiffened column of the composite component, if the cross-section contains multiple materials, the material elastic modulus of the material region where the feature point is located is called through the feature point cloud of the stiffened column in step S1 to perform numerical calculation, and the strain value in the linear strain distribution field is converted into stress value.

[0028] It should be noted that in S33, the generated initial stress point cloud has been updated to the actual stress data after being affected by grouting in terms of physical properties, but the initial stiffening column feature point cloud is still temporarily used in terms of spatial properties. It is precisely because this geometric shape lags behind the misalignment of the stress state that there is a non-zero energy gradient inside the fluid-structure interaction system composed of stiffening columns and grout. This gradient is the direct source of the fluid-structure interaction physical residual term in subsequent steps.

[0029] In this embodiment, S31 includes: S311, acquire strain data of stiffened column through distributed optical fiber sensors, and associate the strain data with the spatial coordinates of the cross section where each sensor is located; S312, Based on the plane section assumption, least squares plane fitting is performed on the spatial coordinates and stiffened column strain data to construct the strain distribution plane equation of the current section; S313, Analyze the strain distribution plane equation, extract the axial reference strain and bending curvature of the cross section, and generate the linear strain distribution field of each cross section of the stiffened column.

[0030] In this embodiment, S312 specifically refers to: Based on the plane section assumption, a general linear equation for the strain distribution of the current section is established; Substitute the spatial coordinates and stiffened column strain data as sample data into the general linear equation to obtain the optimal characteristic coefficients; Substitute the optimal characteristic coefficients back into the general linear equation to construct the strain distribution plane equation for the current section.

[0031] It should be noted that in S312, the plane section assumption is based on the Bernoulli-Euler beam theory, which assumes that the cross section of the stiffened column remains plane after deformation, thereby defining that the strain value of any node in the section and the spatial coordinates of the node satisfy a three-dimensional linear equation relationship, and establishing a general linear equation for the strain distribution of the current section.

[0032] In this embodiment, S313 specifically refers to: Based on the cross-sectional center of the cross-sectional point cloud, the axial reference strain is obtained by substituting it into the strain distribution plane equation. The linear coefficients of the strain distribution plane equation are extracted, and the modulus of the linear coefficients is calculated by vector synthesis to obtain the bending curvature. Substituting the axial reference strain and bending curvature back into the strain distribution plane equation generates the linear strain distribution field of each section of the stiffened column.

[0033] In this embodiment, in step S32, based on the elastic modulus of the stiffening column material, the strain values ​​in the linear strain distribution field are converted into stress values, specifically as follows: Based on the preset elastic modulus of the stiffening column material, Hooke's law is used to linearly scale the strain distribution values ​​in the linear strain distribution field, converting the strain values ​​in the linear strain distribution field into stress values.

[0034] It should be noted that this embodiment transforms sparse distributed optical fiber data into a continuous stress distribution field of stiffened columns by assuming a plane section and using the elastic modulus, and further maps it onto a high-density feature point cloud, enabling the system to acquire the stress state of the entire cross section of the stiffened column, which greatly improves the spatial resolution and completeness of the monitoring data.

[0035] In this embodiment, step S4, variational iterative correction of the initial stress point cloud, specifically involves: S41, map the initial stress point cloud into a state vector; S42, construct the total energy functional of the fluid-structure interaction system, and combine it with the state vector to obtain the fluid-structure interaction physical residual term; S43, with the goal of eliminating the residual term, perform variational iteration on the state vector in the opposite direction of the gradient to obtain a convergent state vector; S44, based on the convergent state vector, performs inverse mapping to obtain the corrected stress point cloud of the stiffening column.

[0036] Preferably, before calculating the energy gradient in step S42, the coordinate components of the corresponding feature points at the bottom of the stiffened column in the state vector need to be locked according to the construction state of the stiffened column, and the displacement increment of the stiffened column is forced to zero during the variational iteration process to eliminate the rigid body displacement mode and ensure that the total energy functional has a unique minimum solution.

[0037] In this embodiment, the variational and iterative process in S43 is specifically as follows: The opposite direction of the energy gradient obtained in step S42 is used as the search direction for energy descent, and the variational increment of the state vector is calculated in combination with the preset step size factor. The current state vector is updated by accumulating the variational increment, and the updated fluid-structure interaction physical residual term is monitored to see if it meets the convergence threshold. The converged state vector is then output.

[0038] It should be noted that in S42, the fluid-structure interaction system is composed of a stiffening column and uncured grouting material inside the pile hole. In this system, the stiffening column and the grouting material interact mechanically through the side surface of the stiffening column, i.e., the coupling interface.

[0039] In this embodiment, S42 includes: S421, based on the initial stress point cloud, performs volume integral operation on the volume domain of the stiffened column to obtain the elastic strain energy; S422, obtain the fluid pressure distribution in the grouting environment, perform integral calculation on the side surface of the stiffened column, and obtain the fluid pressure potential energy; S423, construct a total energy functional that includes elastic strain energy and fluid pressure potential energy, and calculate the energy gradient of the total energy functional with respect to the state vector; S424 constructs a global gradient vector based on the energy gradient and uses the global gradient vector as a physical residual term of fluid-structure interaction.

[0040] In step S421, a volume integral operation is performed on the stiffened column volume domain to obtain the elastic strain energy, specifically: Based on the spatial distribution density of the feature point cloud of the stiffened column, the total volume of the stiffened column is discretized into volume micro-elements corresponding to the feature point cloud. Based on the theory of linear elasticity, the characteristic points and corresponding stress values ​​in the initial stress point cloud are read, and the elastic strain energy density at the characteristic points is calculated by combining the elastic modulus of the stiffened column material. The elastic strain energy of the stiffened column is obtained by weighted summation of the strain energy density and volumetric infinitesimal elements at all characteristic points using a numerical integration method.

[0041] It should be noted that the volume element is defined as the ratio of the total volume of the stiffened column to the total number of feature points, and the formula is: , In the formula, It is a volume infinitesimal element. The total volume of the stiffened column. The total number of feature points in the feature point cloud of the stiff columnar feature point cloud.

[0042] It should be noted that the formula for calculating the elastic strain energy of the stiffened column is as follows: , In the formula, The elastic strain energy of a stiffened column. This is the axial stress value. The elastic modulus of the stiffened column material.

[0043] Optionally, the calculation process for obtaining the elastic strain energy of the stiffened column can be reduced to: , In the formula, Let be the axial stress value at the i-th feature point. The overall value represents the elastic strain energy density at the i-th feature point.

[0044] In step S422, the fluid pressure potential energy is obtained by integrating the surface of the stiffening column, specifically as follows: Based on the feature point cloud of the stiffened column, the side surface of the stiffened column is discretized into a set of area micro-elements with unit normal vectors; Obtain the current grouting parameters and construct the fluid pressure vector field acting on the area micro-element; The work done by the fluid pressure vector field on each area element is calculated, and the fluid pressure potential energy is obtained by global integration and accumulation.

[0045] It should be noted that the grouting parameters include the density of the concrete, the pumping rate, and the additional pressure at the pump inlet.

[0046] It should be noted that in step S423, by using the state vector as the independent variable of the fluid-structure interaction system, the elastic strain energy of step S421 and the fluid pressure potential energy of step S422 are used as scalar functions of the state vector to construct the total energy functional.

[0047] In this embodiment, S423 specifically refers to: Performing a first-order partial differential operation on the total energy functional with respect to the state vector yields the partial derivative of the total energy functional with respect to each element in the state vector, i.e., the energy gradient. All energy gradients obtained from partial differential operations are arranged and combined according to the index order of the state vector to form the global gradient vector.

[0048] It should be noted that the specific formula for calculating the fluid-structure interaction physical residual term is as follows: , In the formula, This is the physical residual term of fluid-structure interaction. For state vectors, It is the total energy functional.

[0049] It should be noted that in this embodiment, a total energy functional containing the elastic strain energy and fluid pressure potential energy of the stiffening column is constructed. Combined with the state vector, a fluid-structure interaction physical residual term is obtained. This term is used to reflect the misalignment of the characteristic point cloud of the stiffening column with the initial stress point cloud in terms of geometric shape and stress state. Then, the fluid-structure interaction physical residual term is eliminated through variational iteration to obtain the physical shape of the stiffening column in the actual grouting process.

[0050] In this embodiment, in step S5, the stiffener feature point cloud and the stiffener stress point cloud are compared to obtain the stiffener pose deviation, specifically as follows: S51, extract the center coordinates and axial angle of the feature point cloud of the stiffening column to determine the target pose of the stiffening column; S52, extract the center coordinates and axial angle of the stress point cloud of the stiffening column to determine the real-time pose of the stiffening column; S53, calculate the spatial deviation of the real-time pose relative to the target pose; S54, decompose the spatial deviation along the local coordinate system of the column to obtain the stiffening column position deviation.

[0051] Optionally, in S51, the method for calculating the axial angle in this embodiment adopts principal component analysis. By constructing the covariance matrix of the feature point cloud, the matrix is ​​decomposed into eigenvalues, and the eigenvector corresponding to the largest eigenvalue is selected as the target axial vector of the stiffening column.

[0052] Furthermore, by calculating the difference vector between the real-time center and the target center, a spatial translation vector is obtained, and the angle between the real-time axial vector and the target axial vector is calculated to characterize the overall tilt of the stiffening column.

[0053] In this embodiment, step S5, adjusting the multi-channel grouting rate according to the deviation difference, includes: S55, calculates the target correction force required for correction based on the stiffening column's positional deviation; S56, based on the coupling relationship between the multi-channel pumping rate and the fluid pressure on the stiffening column, constructs the kinematic Jacobian matrix; S57, based on the kinematic Jacobian matrix, differentially adjusts the pumping rate required for each channel according to the target correction force.

[0054] It should be noted that in step S55, the target correction force required for correction is calculated by combining the stiffener posture deviation output in step S54 with the lateral bending stiffness coefficient of the stiffener, and according to Hooke's Law, the reverse force required to counteract the current deviation is calculated.

[0055] In this embodiment, the kinematic Jacobian matrix is ​​intended to map a linear relationship from grouting rate to force space.

[0056] It should be noted that, in this embodiment, the intertwined pressure effects between each pumping port are mathematically calculated by establishing a kinematic Jacobian matrix. This makes the system no longer blindly inject grout at full speed, but generates a resultant force with controllable direction and adjustable magnitude through asymmetric pressure ratio, correcting the posture deviation of the stiffening column and effectively solving the problem of correction and optimization in the stiffening column embedding process.

[0057] In this embodiment, step S57, which involves differentially adjusting the pumping rate of each channel based on the target correction force, includes: S571, Calculate the generalized inverse matrix of the kinematic Jacobian matrix; S572, the target correction force is solved by linear mapping using the generalized inverse matrix to obtain the pumping rate adjustment vector corresponding to each channel; S573, the pumping rate adjustment vector is superimposed on the current pumping rate of each channel to perform asymmetric control of the fluid pressure distribution in the grouting environment.

[0058] It should be noted that in S572, the pumping rate adjustment increment vector corresponding to each channel should be limited to a small step size range. Since the fluid-structure interaction system of grouting fluid and stiffening column has a high degree of nonlinearity, and a single linear Jacobian matrix can only represent the local linear relationship, it is necessary to limit the step size and use a high-frequency measurement, calculation and fine-tuning process to make the position deviation of the stiffening column gradually converge to zero, so as to avoid system oscillation caused by excessive adjustment amplitude in a single operation.

[0059] The above formulas are all dimensionless calculations. The formulas are derived from software simulations based on a large amount of collected data to obtain the most recent real-world results. The preset parameters in the formulas are set by those skilled in the art according to the actual situation.

[0060] The above embodiments can be implemented, in whole or in part, by software, hardware, firmware, or any other combination thereof. When implemented using software, the above embodiments can be implemented, in whole or in part, in the form of a computer program product.

[0061] Those skilled in the art will recognize that the modules and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0062] In addition, the functional modules in the various embodiments of this application can be integrated into one processing module, or each module can exist physically separately, or two or more modules can be integrated into one module.

[0063] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

[0064] In conclusion, the above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for reinforcing stiffened columns in an integrated underground facility with air-rail-ground connection, characterized in that, Includes the following steps: Acquire the surface point cloud data of the stiffened column and remap it to generate the feature point cloud of the stiffened column; The stiffened column is suspended and inserted into the pile hole and grout is injected; Construct the stress distribution field of the stiffened column, and generate the initial stress point cloud by combining the characteristic point cloud of the stiffened column; Based on the grouting pressure field, the initial stress point cloud is corrected by variational iteration to obtain the stress point cloud of the stiffened column. The variational iteration includes eliminating the fluid-structure interaction physical residual term constructed based on the energy gradient. By comparing the feature point cloud and stress point cloud of the stiffening column, the positional deviation of the stiffening column is obtained, and the multi-channel grouting rate is adjusted according to the deviation difference until the stiffening column is embedded.

2. The method for reinforcing rigid columns in an integrated underground facility with air-rail-ground connection as described in claim 1, characterized in that, The step of acquiring the surface point cloud data of the stiffening column and remapping it to generate the feature point cloud of the stiffening column includes: Acquire the surface point cloud data of the stiffened column, and remap the surface point cloud data using a preset local coordinate system of the column; The surface point cloud data after remapping is divided into equidistant sections along the axial direction to obtain the section point cloud; Feature vectors are constructed using the cross-sectional center and origin of the cross-sectional point cloud to generate the stiffened column feature point cloud.

3. The method for reinforcing rigid columns in an integrated underground facility with air-rail-ground connection according to claim 2, characterized in that, The construction of the stiffened column stress distribution field, combined with the stiffened column characteristic point cloud, to generate an initial stress point cloud includes: Obtain strain data of stiffened columns and fit to generate linear strain distribution fields for each section of the stiffened column; Based on the elastic modulus of the stiffened column material, the strain values ​​in the linear strain distribution field are converted into stress values ​​to construct the stiffened column stress distribution field; The characteristic point cloud of the stiffened column is substituted into the stress distribution field of the stiffened column for numerical solution to generate the initial stress point cloud.

4. The method for reinforcing stiffened columns in an integrated underground facility with air-rail-ground connection as described in claim 3, characterized in that, The process of acquiring strain data of the stiffened column and fitting it to generate a linear strain distribution field for each section of the stiffened column includes: Strain data of the stiffened column is acquired by distributed fiber optic sensors, and the strain data is correlated with the spatial coordinates of the cross section where each sensor is located; Based on the plane section assumption, the spatial coordinates and stiffened column strain data are fitted with least squares plane to construct the strain distribution plane equation of the current section. The strain distribution plane equation is analyzed, the axial reference strain and bending curvature of the cross section are extracted, and the linear strain distribution field of each cross section of the stiffened column is generated.

5. The method for reinforcing stiffened columns in an integrated underground facility with air-rail-ground connection according to claim 4, characterized in that, The variational iterative correction of the initial stress point cloud is specifically as follows: The initial stress point cloud is mapped to a state vector; Construct the total energy functional of the fluid-structure interaction system, and combine it with the state vector to obtain the fluid-structure interaction physical residual term; With the goal of eliminating the residual term, variational iteration is performed on the state vector in the opposite direction of the gradient to obtain a convergent state vector; Based on the convergent state vector, an inverse mapping is performed to obtain the corrected stress point cloud of the stiffened column.

6. The method for reinforcing stiffened columns in an integrated underground facility with air-rail-ground connection according to claim 5, characterized in that, The total energy functional of the constructed fluid-structure interaction system, combined with the state vector, yields the fluid-structure interaction physical residual term, including: Based on the initial stress point cloud, a volume integral operation is performed on the volume domain of the stiffened column to obtain the elastic strain energy; The fluid pressure distribution in the grouting environment is obtained, and the fluid pressure potential energy is obtained by integral calculation on the side surface of the stiffened column. Construct a total energy functional that includes elastic strain energy and fluid pressure potential energy, and calculate the energy gradient of the total energy functional with respect to the state vector; Based on the energy gradient, a global gradient vector is constructed and used as a physical residual term in fluid-structure interaction.

7. The method for reinforcing rigid columns in an integrated underground facility with air-rail-ground connection according to claim 6, characterized in that, The comparison of the feature point cloud and stress point cloud of the stiffening column yields the pose deviation of the stiffening column, specifically: Extract the center coordinates and axial angle of the feature point cloud of the stiffening column to determine the target pose of the stiffening column; Extract the center coordinates and axial angle of the stress point cloud of the stiffening column to determine the real-time pose of the stiffening column; Calculate the spatial deviation of the real-time pose relative to the target pose; The spatial deviation is decomposed along the local coordinate system of the cylinder to obtain the stiffening cylinder pose deviation.

8. The method for reinforcing stiffened columns in an integrated underground facility with air-rail-ground connection according to claim 7, characterized in that, The method of adjusting the multi-channel grouting rate based on the deviation difference includes: The target correction force required for correction is calculated based on the stiffening column's posture deviation. Based on the coupling relationship between the multi-channel pumping rate and the fluid pressure on the side surface of the stiffened column, a kinematic Jacobian matrix is ​​constructed. Based on the kinematic Jacobian matrix, the pumping rate required for each channel is differentially adjusted according to the target correction force.

9. The method for reinforcing rigid columns in an integrated underground facility with air-rail-ground connection according to claim 8, characterized in that, The specific method for adjusting the pumping rate of each channel based on the target correction force differential is as follows: Calculate the generalized inverse of the kinematic Jacobian matrix; The target correction force is solved by linear mapping using the generalized inverse matrix to obtain the pumping rate adjustment vector corresponding to each channel. The pumping rate adjustment vector is superimposed on the current pumping rate of each channel to perform asymmetric control of the fluid pressure distribution in the grouting environment.

Images