Python-Fish-based double-line shield tunnel parameterization automatic simulation method and product

Abstract

The invention provides a Python-Fish-based double-line shield tunnel parameterization automatic simulation method and product, and relates to the technical field of tunnel engineering numerical simulation. According to the method, an efficient Python parameterized modeling framework is provided, parameterized control over the geometric model of the double-line shield tunnel can be achieved based on the geometric parameters of the target double-line shield tunnel, and therefore automatic grid generation under different layouts of the double-line shield tunnel is achieved. Multi-parameter driving simulation based on Fish is provided, calculation simulation parameters of the double-line shield tunnel can be determined on the basis of engineering working conditions actually needing to be simulated, and automatic simulation and multi-working-condition batch simulation of different layouts and different construction working conditions of the double-line shield tunnel are achieved on the basis of the calculation simulation parameters of the double-line shield tunnel. And parameterized geometric modeling and multi-parameter driven automatic simulation of the double-line shield tunnel under complex working conditions are realized.

Description

TECHNICAL FIELD

[0001] The embodiment of the application relates to the technical field of tunnel engineering numerical simulation, and particularly relates to a double-line shield tunnel parameterized automatic simulation method based on Python-Fish and a product. BACKGROUND

[0002] In the field of tunnel engineering, shield tunnel construction is a widely used efficient construction method. With the continuous deepening of urban underground space development, the demand for double-line shield tunnel construction is increasing. In order to ensure construction safety and optimize design schemes, numerical simulation technology is increasingly important in double-line shield tunnel construction. Fine three-dimensional numerical simulation can accurately simulate the stress and deformation distribution of surrounding soil and tunnel structure during double-line shield tunnel construction. Through numerical simulation, the influence of double-line shield tunnel spatial layout on the deformation of surrounding soil and tunnel structure can be studied to provide a scientific basis for optimizing design. The influence of different construction parameters (such as face pressure, grouting pressure, and jack thrust) on tunnel construction can be studied to determine the optimal combination of construction parameters. The influence of double-line shield tunnel construction sequence and lag distance on surrounding soil and existing structures can be studied to determine reasonable construction sequence and lag distance. The interaction mechanism and construction disturbance mechanism of double-line shield tunnel can be deeply understood to provide strong support for the optimization design and construction safety control of double-line shield tunnel.

[0003] Currently, FLAC3D has the following shortcomings when modeling and analyzing double-line shield tunnels: (1) Geometric modeling is difficult. FLAC3D has weak pre-processing capabilities. Geometric modeling using Fish command streams is limited by the stacking of built-in basic models. For complex geometric structures such as double-line shield tunnels with arbitrary layout, the modeling process is complex and prone to errors, making it difficult to automatically generate mesh models and requiring a lot of manual intervention, resulting in low modeling efficiency; (2) Parameter adjustment is not flexible. Existing methods are mostly based on specific engineering simulation under specific conditions. After adjusting engineering parameters or condition parameters, the process and simulation need to be started from scratch. It is difficult to achieve rapid adjustment and automatic simulation process of multiple parameters, and it is difficult to achieve batch simulation of multiple conditions for parameter sensitivity analysis and optimization. SUMMARY

[0004] The embodiment of the application provides a double-line shield tunnel parameterized automatic simulation method and product based on Python-Fish, which efficiently and automatically performs modeling, construction simulation, and batch simulation of multiple conditions of double-line shield tunnels.

[0005] The first aspect of the embodiment of the application provides a double-line shield tunnel parameterized automatic simulation method based on Python-Fish, which comprises the following steps: The geometric parameters of the target double-line shield tunnel are imported into a numerical simulation grid generator, an automatic modeling Python script is run, a parameterized double-line shield tunnel geometric model is established, three-dimensional grid division is performed, and a grid model file is exported; In FLAC3D, the grid model file is imported by using a Fish command; According to the target engineering working condition, double-line shield tunnel calculation simulation parameters are assigned, including constitutive model parameters, shield construction parameters, support parameters, and double-line shield tunnel lag distance. Based on the double-line shield tunnel calculation simulation parameters and the grid model file, double-line shield tunnel excavation cycle simulation is performed to generate a simulation result file; The simulation result file is imported into FLAC3D for viewing and analyzing the stress and deformation of the soil and structure.

[0006] Optionally, in FLAC3D, the grid model file is imported by using a Fish command, including: The grid model file is imported into FLAC3D, and the geometric model is grouped to divide the leading tunnel T1, the trailing tunnel T2, and the surrounding soil S into three groups.

[0007] Optionally, the geometric parameters of the target double-line shield tunnel include model size L, width W, height H, tunnel diameter D, tunnel burial depth C, double-line tunnel center distance B, double-line tunnel offset angle θ, and grid size Δ. L, W, and H are used to limit the external boundary of the double-line shield tunnel geometric model. D and C are used to limit the diameter and vault burial depth of the double-line shield tunnel geometric model, respectively. B and θ limit the position of the two tunnels in the double-line shield tunnel geometric model, and the center coordinates of the two tunnels are determined by the following formula: ; ; ; ; Δ is used to control the size of three-dimensional grid division.

[0008] Optionally, the soil constitutive model uses a small-strain plastic hardening model; the shield construction parameters include face pressure, grouting pressure, and jack thrust; the support parameters include the stiffness of the lining segment and the elastic modulus of the grouting; and the double-line shield tunnel lag distance includes the distance between the leading tunnel and the trailing tunnel face.

[0009] Optionally, the double-line shield tunnel excavation cycle simulation is performed based on the double-line shield tunnel calculation simulation parameters and the grid model file, and the double-line shield tunnel excavation cycle simulation comprises the following steps: According to the constitutive model parameters, the physical and mechanical parameters of the three-dimensional numerical model are automatically assigned, and the initial stress balance of the model is completed; according to the shield construction parameters, the support parameters and the double-line shield tunnel lag distance, the double-line shield tunnel cycle excavation simulation is automatically performed to generate a simulation result file.

[0010] Optionally, according to the shield construction parameters, the support parameters and the double-line shield tunnel lag distance, the double-line shield tunnel cycle excavation simulation is automatically performed, and the double-line shield tunnel cycle excavation simulation comprises the following steps: The face pressure is determined as the horizontal stress of the trapezoidal distribution acting on the face vertically, and increases linearly along the depth; The shield machine is simulated by using a shell element, the stiffness decreases from front to back to simulate the taper effect, and the self-weight pressure is the vertical stress acting on the length and 90° cross-section range of the shield machine; The jack thrust acts on the newly spliced ring segment as a uniformly distributed nodal concentrated force along the circumference of the segment; The grouting pressure acts within the width range of the newly spliced ring segment as the normal pressure applied to the surrounding rock and the segment, and increases linearly along the depth; The trolley pressure is determined as the vertical pressure acting on the two parallel lines of the segment; The lining segment is simulated by using a Liner element, and the joint connection between the segments and the interaction between the segments and the grouting are considered; The grouting body is simulated by using a linear elastic solid element, the hardening process of the grouting body is considered, the grouting body of the first ring after the shield tail is dragged out is fresh grout, and lower strength parameters are assigned to the fresh grout; the fresh grout is converted into hardened grouting body after two excavation steps, and higher strength parameters are assigned to the hardened grouting body.

[0011] The second aspect of the embodiment of the application provides a double-line shield tunnel parameterized automatic simulation device based on Python-Fish, and the device comprises: A model establishment module is configured to import the geometric parameters of a target double-line shield tunnel into a numerical simulation grid generator, run an automatic modeling Python script, establish a parameterized control double-line shield tunnel geometric model, perform three-dimensional grid division and export a grid model file; An import module is configured to import the grid model file by using a Fish command in FLAC3D; An assignment module is configured to assign double-line shield tunnel calculation simulation parameters according to a target engineering working condition, wherein the double-line shield tunnel calculation simulation parameters comprise constitutive model parameters, shield construction parameters, support parameters and double-line shield tunnel lag distance; An analog module is configured to perform a double-line shield tunnel cyclic excavation simulation based on the double-line shield tunnel calculation simulation parameters and the grid model file, and generate an analog result file. A display module is configured to import the analog result file into FLAC3D, and view and analyze the stress and deformation of the soil body and the structure.

[0012] In a third aspect, an electronic device is provided, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the double-line shield tunnel parameterized automatic simulation method based on Python-Fish when executed.

[0013] In a fourth aspect, a computer readable storage medium is provided, which stores a computer program, and the program is executable on a processor to implement the double-line shield tunnel parameterized automatic simulation method based on Python-Fish.

[0014] In a fifth aspect, a computer program product is provided, which includes computer programs / instructions, and the computer programs / instructions are executable on a processor to implement the steps in the double-line shield tunnel parameterized automatic simulation method based on Python-Fish.

[0015] In the embodiments of the present application, an efficient Python parameterized modeling framework is provided, which can realize parameterized control of the double-line shield tunnel geometric model based on the geometric parameters of the target double-line shield tunnel, thereby realizing automatic mesh generation of the double-line shield tunnel under different layouts.

[0016] In the embodiments of the present application, Fish-based multi-parameter driven simulation is provided, which can determine the double-line shield tunnel calculation simulation parameters based on the actual engineering working conditions to be simulated, and realize automatic simulation and multi-working condition batch simulation of the double-line shield tunnel under different layouts and different construction working conditions based on the double-line shield tunnel calculation simulation parameters.

[0017] Therefore, in the embodiments of the present application, the automatic mesh generation of the double-line shield tunnel with complex geometric structure is realized through the parameterized control of the Python script and the powerful modeling capability of the numerical simulation mesh generator, which significantly improves the modeling efficiency. The dynamic adjustment and automatic simulation of multiple parameters are realized based on the double-line shield tunnel calculation simulation parameters through the parameterized control of Fish, which can quickly generate simulation results under different engineering working conditions and improve the simulation efficiency. It is convenient to realize rapid and large numerical simulation results for the double-line shield tunnel, which is beneficial to data mining and optimization design, and improves the scientificity and economy of the construction scheme.

[0018] In summary, the embodiment of the present application realizes parameterized geometric modeling and automatic simulation driven by multiple parameters of a double-line shield tunnel under complex working conditions. BRIEF DESCRIPTION OF DRAWINGS

[0019] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings needed to be used in the description of the embodiments of the present application will be briefly introduced as follows. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can also be obtained by those skilled in the art without creative labor on the basis of these drawings.

[0020] Figure 1 is a step flow chart of the double-line shield tunnel parameterized automatic simulation method based on Python-Fish provided by the embodiment of the present application; Figure 2 is an implementation flow chart of the double-line shield tunnel parameterized automatic simulation method based on Python-Fish provided by the embodiment of the present application; Figure 3 is an exemplary schematic diagram of geometric parameter control of a double-line shield tunnel model in the double-line shield tunnel parameterized automatic simulation method based on Python-Fish provided by the embodiment of the present application; Figure 4 is an exemplary geometric parameter input diagram in the double-line shield tunnel parameterized automatic simulation method based on Python-Fish provided by the embodiment of the present application; Figure 5 is an exemplary grid model diagram automatically generated in the double-line shield tunnel parameterized automatic simulation method based on Python-Fish provided by the embodiment of the present application; Figure 6 is an exemplary input diagram of double-line shield tunnel calculation simulation parameters in the double-line shield tunnel parameterized automatic simulation method based on Python-Fish provided by the embodiment of the present application; Figure 7 is a shield-stratum interaction model schematic diagram used in the double-line shield tunnel parameterized automatic simulation method based on Python-Fish provided by the embodiment of the present application; Figure 8 is an exemplary result diagram of the double-line shield tunnel parameterized automatic simulation method based on Python-Fish provided by the embodiment of the present application; Figure 9 is a hardware structure diagram of any device with data processing capability on which the double-line shield tunnel parameterized automatic simulation device based on Python-Fish provided by the embodiment of the present application is located. DETAILED DESCRIPTION

[0021] In order to make the above objectives, characteristics and advantages of the present application more obvious and easy to understand, the present application will be further described in detail below with reference to the drawings and specific embodiments.

[0022] As shown in Figure 1 The method comprises the following steps: S101, the geometric parameters of the target double-line shield tunnel are imported into a numerical simulation grid generator, an automatic modeling Python script is run, a parameterized control double-line shield tunnel geometric model is established, three-dimensional grid division is performed, and a grid model file is exported.

[0023] In the embodiment of the present application, an existing numerical simulation grid generator (for example, Gmsh library) can be used to run an automatic modeling Python script to generate a parameterized control double-line shield tunnel geometric model based on the geometric parameters of the target double-line shield tunnel.

[0024] FLAC3D is a professional geotechnical engineering numerical analysis software that integrates a Python interpreter. With this feature, users can use Python language to write scripts to control various functions of FLAC3D.

[0025] Gmsh is an open source CAD and mesh generation tool that provides a Python API for creating complex geometric models and performing mesh division. In the embodiment of the present application, Gmsh can be used in the Python environment of FLAC3D.

[0026] The construction of the double-line shield tunnel geometric model depends on the geometric parameters of the target double-line shield tunnel, mainly including: Geometric parameters of a single tunnel, such as the outer diameter, inner diameter of the tunnel, and the thickness of the lining, etc.

[0027] Spatial position parameters of the tunnel, such as the horizontal and vertical spacing between the two tunnels, and the burial depth of the tunnel, etc.

[0028] Tunnel orientation parameters, including the length, slope, and turning radius of the tunnel, etc.

[0029] Other parameters, such as segment length, shape of the excavation face, etc.

[0030] In the Python script, these parameters can be defined as variables for easy modification, thereby achieving parameterized control of the model.

[0031] In the embodiment of the application, parameterized control means that the double-line shield tunnel geometric model can be flexibly adjusted by modifying parameters, without the need to recompile the entire script, and the combined parameters can be defined as variables or function parameters to easily realize parameterized control.

[0032] In the embodiment of the application, Gmsh provides a plurality of meshing algorithms, and a suitable algorithm can be selected according to the characteristics of the model. Before meshing, a mesh size parameter needs to be set.

[0033] Specifically, in the embodiment of the application, the geometric parameters of the target double-line shield tunnel include: model size L, W, H; tunnel diameter D, tunnel burial depth C, double-line tunnel center distance B, double-line tunnel offset angle θ and mesh size Δ. L, W and H are used to limit the external boundary of the double-line shield tunnel geometric model. D and C are respectively used to limit the diameter and vault burial depth of the double-line shield tunnel geometric model. B and θ limit the positions of the two tunnels in the double-line shield tunnel geometric model, and the center coordinates of the two tunnels are determined by the following formula: ; ; ; ; Δ is used to control the size of three-dimensional meshing.

[0034] In the embodiment of the application, the double-line shield tunnel geometric model can be flexibly adjusted by modifying the above parameters, without the need to recompile the entire script.

[0035] S102, in FLAC3D, the mesh model file is imported by using a Fish command.

[0036] Fish is a programming language built in FLAC3D, which can realize automatic operation, parameterized modeling and complex numerical calculation. Fish language has many functions, such as variable definition, function writing, loop control structure, and can call various commands of FLAC3D. To execute Fish code, it can be directly input in the command line of FLAC3D, or the code can be written in a file, and then the code is run by using the program call command.

[0037] In step S101, the mesh model file (.inp format) has been generated by using Gmsh, and then the file can be imported in FLAC3D by using a Fish command.

[0038] In geotechnical engineering analysis, in order to give different material properties to different regions, set boundary conditions or carry out result analysis, it is necessary to group and name entities (such as regions, contact surfaces, nodes, etc.) in the grid model. Through reasonable grouping and naming of the grid model, the management of the model can be more clear, and the analysis and calculation can be more efficient.

[0039] Specifically, step S102 includes: importing the grid model file into FLAC3D, grouping the geometric model, and dividing the first tunnel T1, the second tunnel T2 and the surrounding soil S into three groups.

[0040] In the embodiment of the application, the gmsh.model.geo.addPoint(), gmsh.model.geo.addLine(), gmsh.model.geo.addCircleArc(), gmsh.model.geo.addPlaneSurface() commands can be used to generate point, straight line, arc line and surface models respectively, the gmsh.model.geo.extrude() command can be used to stretch the surface model into a body model, and the gmsh.model.mesh.generate() command can be used to automatically divide the grid.

[0041] Specifically, as shown in Figure 2 As shown in the flowchart of the implementation of the double-line shield tunnel parameterization automatic simulation method based on Python-Fish provided by the embodiment of the application, in the embodiment of the application, based on Python control, based on the necessary libraries Gmsh, itasca, math and model set parameters, the model construction script model.py is executed, the required geometric model is automatically generated, the grid division is automatically performed, the required grid model is generated, and the grid model "model.inp" file is obtained.

[0042] In the embodiment of the application, the "zone import" command can be used to import the grid file into FLAC3D, and the "zone group range" command can be used to group the model.

[0043] S103, according to the target engineering working condition, the double-line shield tunnel calculation simulation parameters are valued, and the double-line shield tunnel calculation simulation parameters include: constitutive model parameters, shield construction parameters, support parameters and double-line shield tunnel lag distance.

[0044] In the embodiment of the present application, the soil constitutive model adopts a small-strain plastic hardening model, and specific parameters can be selected according to an actual survey report; shield construction parameters include face pressure, grouting pressure and jack thrust; support parameters include the rigidity of a lining segment and the elastic model of grouting; and the double-line shield tunnel lag distance includes the distance between the faces of the preceding tunnel and the following tunnel.

[0045] Specifically, as shown in the figure, Figure 2 In the embodiment of the present application, based on Fish control, the grid model "model.inp" can be imported into FLAC3D, and simulation calculation parameters, constitutive model parameters, shield construction parameters (face pressure, grouting pressure and jack thrust), support parameters (lining rigidity and grouting rigidity) and double-line shield tunnel lag distance are defined.

[0046] In the embodiment of the present application, the above parameters can be freely changed through parameterization processing to adapt to different geological conditions and construction schemes.

[0047] In the embodiment of the present application, based on the double-line shield tunnel calculation simulation parameters, the dynamic adjustment and automatic simulation of multiple parameters are realized by using the parameterization control of Fish, the simulation results under different engineering conditions can be quickly generated, and the simulation efficiency is improved.

[0048] S104, based on the double-line shield tunnel calculation simulation parameters and the grid model file, a double-line shield tunnel excavation cycle simulation is performed to generate a simulation result file.

[0049] The numerical simulation of the shield tunnel excavation process is a dynamic process, mainly including the following aspects: Step-by-step excavation: the tunnel excavation process is divided into multiple steps according to a certain length (for example, the length of each ring segment) and is performed in sequence.

[0050] Stress release: excavation will cause changes in the stress of the surrounding soil, and this stress release process needs to be simulated.

[0051] Support construction: as the tunnel is excavated, lining structures need to be constructed in time to provide support.

[0052] Soil-structure interaction: the interaction between soil and lining is simulated.

[0053] In the embodiment of the present application, a shield-stratum interaction model is used to carry out double-line shield tunnel excavation cycle simulation.

[0054] In the numerical simulation of shield tunnel construction, the shield-stratum interaction model is used to simulate the mechanical behavior and mutual influence between the shield machine, the lining structure and the surrounding rock-soil mass.

[0055] During shield tunnel construction, the shield-ground interaction mainly includes the following aspects: Shield thrust force: The shield machine provides thrust force through jacks to overcome the resistance of the surrounding soil and advance forward.

[0056] Shield shell friction: There is friction between the shield shell and the surrounding soil, which affects the stress state and deformation of the soil.

[0057] Excavation face stability: The excavation face of the shield machine needs to maintain a certain pressure to balance the pressure of the soil in front and prevent soil collapse or heave.

[0058] Lining interaction with soil: As the shield machine advances, the lining structure (such as segments) is installed behind the shield tail and interacts with the surrounding soil.

[0059] Shield tail gap: After the shield machine passes, a shield tail gap is formed between the lining and the soil, which needs to be filled with grouting in a timely manner.

[0060] In numerical simulation, the shield-ground interaction model needs to consider the following key elements: Geometric model: Shield machine: Usually simplified as a rigid cylinder, considering its outer diameter, length and shape.

[0061] Lining structure: Can be simulated as an elastic body or a combined structure (such as the joint between segments).

[0062] Surrounding soil: Determine the calculation range and boundary conditions according to the actual engineering situation.

[0063] Material constitutive model: Soil: Can choose linear elastic model, elastoplastic model (such as Mohr-Coulomb model, Drucker-Prager model) or more complex constitutive model (such as plastic hardening model, critical state model).

[0064] Lining: Generally uses linear elastic model, but also can consider the nonlinear characteristics of concrete.

[0065] Contact surface: Simulate the relative sliding and separation between the shield machine, lining and soil, usually use Coulomb friction model.

[0066] Load and boundary conditions: Excavation face pressure: Simulate the support pressure applied by the shield machine excavation face.

[0067] Shield shell friction: Friction distributed along the shield shell surface.

[0068] Grouting pressure: Pressure generated by grouting in the shield tail gap.

[0069] Boundary conditions: Displacement constraints or stress boundary conditions are usually imposed on the model boundaries.

[0070] In the FLAC3D, a script is written by using a Fish language to realize the excavation cycle simulation of the shield tunnel.

[0071] Specifically, the step S104 comprises: S1041, automatically assigning the physical and mechanical parameters of the three-dimensional numerical model according to the constitutive model parameters, and completing the initial ground stress balance of the model.

[0072] In the embodiment of the application, the mechanical behavior of the rock-soil body is described by the constitutive model, and a specific constitutive model can be selected according to actual needs.

[0073] Through the automatic assignment of parameters and the initial ground stress balance, the numerical model can accurately reflect the actual mechanical state of the project, and provide a reliable foundation for subsequent shield construction simulation.

[0074] S1042, automatically performing the double-line shield tunnel cycle excavation simulation according to the shield construction parameters, the support parameters and the double-line shield tunnel lag distance, and generating a simulation result file.

[0075] Specifically, the step S1042 comprises: The face pressure is determined as the horizontal stress of the trapezoidal distribution acting on the face, and increases linearly along the depth.

[0076] The shield machine is simulated by using a shell element, the stiffness increases from front to back to simulate the taper effect, and the self-weight pressure is the vertical stress acting on the length and 90° cross-section range of the shield machine.

[0077] The jack thrust acts on the newly spliced ring segment, and is a uniformly distributed nodal concentrated force along the circumference of the segment.

[0078] The grouting pressure acts within the width range of the newly spliced ring segment, and is the normal pressure applied to the surrounding rock and the segment, and increases linearly along the depth.

[0079] The trolley pressure is determined as the vertical pressure acting on the two parallel lines of the segment.

[0080] The lining segment is simulated by using a Liner element, and the joint connection between the segments and the interaction between the segments and the grouting are considered.

[0081] The grouting body is simulated by using a linear elastic solid element, the grouting hardening process is considered, the first ring after the shield tail is dragged out is fresh grout, and is given a lower strength parameter, and the fresh grout is converted into hardened grouting body after two excavation steps, and is given a higher strength parameter.

[0082] AsFigure 2 As shown, in the embodiment of the present application, simulation calculation script simulation.f3dat can be executed based on Fish control to automatically complete step-by-step simulation, and save result file result.sav.

[0083] S105, import the simulation result file to FLAC3D, and view and analyze the stress and deformation of the soil body and structure.

[0084] In the embodiment of the present application, the simulation result file is imported to FLAC3D, and the stress and deformation of the soil body and structure can be viewed and analyzed.

[0085] For ease of understanding, an exemplary case is provided to illustrate the Python-Fish-based double-line shield tunnel parameterized automatic simulation method. Specifically, taking the double-line shield tunnel project of Huangniupu-Huangjiang Central Station of Dongguan Rail Transit Line 1 as an example, the specific steps are as follows: 1) Import Gmsh library by using the Python module of FLAC3D to build an environment for subsequent automatic modeling. Input the geometric parameters of the double-line shield tunnel model, and the required geometric control parameters and their meanings are as shown in the table. Figure 3 As shown. Figure 4 The geometric parameters input in this embodiment. Run the automatic modeling Python script, and the script execution steps are as follows: a. According to L, W, and H, the outer boundary of the model is limited, and according to B and θ, the coordinates of the two tunnels are determined:

[0086]

[0087]

[0088]

[0089] According to D and C, the diameter and vault depth of the tunnel are determined.

[0090] Use the gmsh.model.geo.addPoint() and gmsh.model.geo.addLine() commands to generate point and straight line models to create the outer boundary of the model; use the gmsh.model.geo.addPoint() and gmsh.model.geo.addCircleArc() commands to generate point and arc line models to create the tunnel boundary; use the gmsh.model.geo.addPlaneSurface() command to generate surface models; and use the gmsh.model.geo.extrude() command to stretch the surface models into body models.

[0091] The grid size is determined according to delta, and the grid is automatically divided by using the gmsh.model.mesh.generate() command to generate a double-line shield tunnel grid model.

[0092] The double-line shield tunnel grid model is exported by using the gmsh.write('''{}\model.inp'''.format(cur_dir)) command to obtain a grid model file "model.inp".

[0093] 2) The grid file is imported into FLAC3D by using the "zone import'model.inp' format abaqus" command, and the model is grouped by using the "zone group range" command to divide the first tunnel T1, the second tunnel T2 and the surrounding soil S into three groups. The grid model diagram completed by importing and grouping in this embodiment is shown in Figure 5

[0094] 3) The double-line shield tunnel calculation simulation parameters are input, and the calculation simulation parameters input in this embodiment are shown in Figure 6 After the calculation simulation parameters are input, the Fish script is run to automatically complete the initial ground stress balance and double-line shield tunnel excavation cycle simulation. The shield tunnel excavation cycle simulation is carried out by using a shield-ground interaction model, and the schematic diagram of the shield-ground interaction model used in the shield tunnel excavation cycle simulation is shown in Figure 7 ​As shown, specifically: 1) the face pressure is the horizontal stress acting on the face in trapezoidal distribution, which increases linearly with depth, and in this embodiment, it is 413.22 kPa; 2) the shield machine is simulated by shell elements, and the stiffness decreases from front to back to simulate the taper effect, and the self-weight pressure is the vertical stress acting on the length and 90° cross-sectional range of the shield machine; 3) the jack thrust is the nodal concentrated force acting on the newly spliced ring segment, which is uniformly distributed along the circumference of the segment, and in this embodiment, it is 2000 kN; 4) the grouting pressure acts on the width range of the newly spliced ring segment, which is the normal pressure applied to the surrounding rock and the segment, and increases linearly with depth, and in this embodiment, it is 688.80 kPa; 5) the trolley pressure is the vertical pressure acting on the two parallel lines of the segment; 6) the lining segment is simulated by Liner elements, and the joint connection between the segments and the interaction between the segments and the grouting are considered, and in this embodiment, the lining stiffness is 31 GPa; 7) the grouting body is simulated by linear elastic solid elements, considering the grouting hardening process, the first ring after the shield tail is dragged out, the grouting body is fresh grout, and is given a lower strength parameter, and the fresh grout is converted into hardened grouting body after two excavation steps, and is given a higher strength parameter, and in this embodiment, it is 18 MPa. The simulation results are saved in the simulation result file "result.sav".

[0095] 4) FLAC3D imports the simulation result file "result.sav", and views and analyzes the simulation results of the double-line shield tunnel construction. As shown in Figure 8 the upper part is the ground settlement, and the lower part is the maximum bending moment of the lining.

[0096] Based on the same inventive concept, the embodiments of the present application also provide a Python-Fish-based double-line shield tunnel parameterized automatic simulation device, which comprises: A model establishment module is configured to import the geometric parameters of a target double-line shield tunnel into a numerical simulation grid generator, run an automatic modeling Python script, establish a parameterized control double-line shield tunnel geometric model, perform three-dimensional grid division, and export a grid model file. An import module is configured to import the grid model file in FLAC3D by using a Fish command. A value assignment module is configured to assign double-line shield tunnel calculation simulation parameters according to a target engineering working condition, wherein the double-line shield tunnel calculation simulation parameters comprise: constitutive model parameters, shield construction parameters, support parameters, and double-line shield tunnel lag distance. A simulation module is configured to perform double-line shield tunnel cyclic excavation simulation based on the double-line shield tunnel calculation simulation parameters and the grid model file, and generate a simulation result file. A display module is configured to import the simulation result file into FLAC3D, and view and analyze the stress and deformation of the soil body and the structure.

[0097] Optionally, the import module is configured to import the grid model file into FLAC3D, and group the geometric model to divide the preceding tunnel T1, the following tunnel T2 and the surrounding soil S into three groups.

[0098] Optionally, the geometric parameters of the target double-line shield tunnel include model size L, W and H, tunnel diameter D, tunnel burial depth C, double-line tunnel center distance B, double-line tunnel offset angle θ and grid size Δ. L, W and H are used to limit the outer boundary of the double-line shield tunnel geometric model. D and C are used to limit the diameter and vault burial depth of the double-line shield tunnel geometric model, respectively. B and θ limit the position of the two tunnels in the double-line shield tunnel geometric model, and the center coordinates of the two tunnels are determined by the following formula: ; ; ; ; Δ is used to control the size of the three-dimensional grid division.

[0099] Optionally, the soil constitutive model adopts a small-strain plastic hardening model; the shield construction parameters include face pressure, grouting pressure and jack thrust; the support parameters include the stiffness of the lining segment and the elastic model of the grouting; and the double-line shield tunnel lag distance includes the distance between the faces of the preceding tunnel and the following tunnel.

[0100] Optionally, the simulation module is configured to: automatically assign the physical and mechanical parameters of the three-dimensional numerical model according to the constitutive model parameters, complete the initial stress balance of the model, and automatically execute the double-line shield tunnel cyclic excavation simulation according to the shield construction parameters, the support parameters and the double-line shield tunnel lag distance, to generate a simulation result file.

[0101] Optionally, the simulation module is configured to: determine the face pressure as a trapezoidal distribution of horizontal stress acting vertically on the face, which increases linearly along the depth; simulate the shield machine by using a shell element, the stiffness of which increases from front to back to simulate the taper effect, and the self-weight pressure of which is the vertical stress acting on the length and 90° cross-section range of the shield machine; act the jack thrust on the latest spliced ring segment as a uniformly distributed nodal concentrated force along the segment circumference; The grouting pressure is applied to the latest spliced annular segment in the range of the segment width, and the normal pressure applied to the surrounding rock and the segment linearly increases with the depth; The trolley pressure is determined as the vertical pressure applied to the two parallel lines of the segment; The liner unit is used to simulate the lining segment, and the joint connection between the segments and the interaction between the segments and the grouting are considered; The linear elastic solid element is used to simulate the grouting body, the grouting hardening process is considered, the first ring after the shield tail is dragged out, the grouting body is fresh grout, and lower strength parameters are given, and the fresh grout is converted into hardened grouting body after two excavation steps, and higher strength parameters are given.

[0102] The embodiment of the device provided by the application can be applied to any device with data processing capability, such as a computer or the like. The device embodiment can be realized by software, or by hardware or a combination of software and hardware. Taking software realization as an example, as a logical device, it is formed by reading the corresponding computer program instructions in the non-volatile memory into the memory and running by the processor of the device with data processing capability. From the hardware level, as shown in Figure 9 As shown in the figure, in addition to the processor, the memory, the network interface, and the non-volatile memory shown in the figure, the device in the embodiment is usually provided with other hardware according to the actual functions of the device with data processing capability, and details are not repeated. Figure 9

[0103] Based on the same inventive concept, the embodiment of the application further provides an electronic device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the steps in the Python-Fish-based double-line shield tunnel parameterized automatic simulation method according to any of the above embodiments when executed.

[0104] Based on the same inventive concept, the embodiment of the application further provides a computer readable storage medium, which stores a computer program, and the program implements the steps in the Python-Fish-based double-line shield tunnel parameterized automatic simulation method according to any of the above embodiments when executed by a processor.

[0105] ​Based on the same inventive concept, the embodiments of the present application provide a computer program product comprising computer programs / instructions which, when executed by a processor, implement the steps of the Python-Fish-based double-line shield tunnel parameterized automatic simulation method according to any of the above embodiments.

[0106] Each of the embodiments in the specification is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other.

[0107] Those skilled in the art should understand that the embodiments of the present application can be provided as a method, device, or computer program product. Therefore, the embodiments of the present application can be in the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the embodiments of the present application can be in the form of a computer program product implemented on one or more computer usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) containing computer usable program code.

[0108] The embodiments of the present application are described with reference to flowcharts and / or block diagrams according to the method, terminal device (apparatus), and computer program product of the embodiments of the present application. It should be understood that each flow and / or block in the flowcharts and / or block diagrams, and the combination of the flows and / or blocks in the flowcharts and / or block diagrams can be implemented by computer program instructions. These computer program instructions can be provided to a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable terminal device to produce a machine, so that the instructions executed by the processor of the computer or other programmable terminal device produce a device that implements the functions specified in the flowcharts and / or block diagrams. Figure 1 The functions specified in one flow or multiple flows and / or blocks Figure 1 The functions specified in one flow or multiple flows and / or blocks

[0109] These computer program instructions can also be stored in a computer readable storage medium that can guide the computer or other programmable terminal device to work in a specific manner, so that the instructions stored in the computer readable storage medium produce a product including instruction devices that implement the functions specified in the flowcharts and / or block diagrams. Figure 1 The functions specified in one flow or multiple flows and / or blocks Figure 1 The functions specified in one flow or multiple flows and / or blocks

[0110] These computer program instructions can also be loaded into a computer or other programmable terminal device, so that a series of operation steps are performed on the computer or other programmable terminal device to produce a computer-implemented process, so that the instructions executed on the computer or other programmable terminal device provide a process for implementing the functions specified in the flowcharts and / or block diagrams. Figure 1 The functions specified in one flow or multiple flows and / or blocksFigure 1 the steps of the function specified in the one or more blocks.

[0111] While preferred embodiments of the application have been described, those skilled in the art will appreciate that other modifications and variations to the preferred embodiments are possible in light of the above teachings. It is, therefore, to be understood that within the scope of the appended claims, the application can be practiced otherwise than as specifically described.

[0112] Finally, it should be noted that the terms "first", "second", and the like, herein do not denote any order, quantity, combination, or importance, but rather are used to distinguish one element from another, and do not imply singular or plural. Moreover, the term "include", "have", or "contain" or any other variant thereof, are intended to cover a non-exclusive inclusion, such that processes, methods, articles, or apparatuses that comprise a set of elements not expressly listed are also within the scope of the present application. Further, any statement of a means plus function format is intended to cover structures, apparatuses, and / or methods for performing the recited function and / or achieving the recited result.

[0113] The above describes in detail a Python-Fish-based double-line shield tunnel parameterization automatic simulation method provided by the present application, and the principles and implementation manners of the present application are described by using specific examples. The above description of the embodiments is only used to help understand the method of the present application and its core idea; meanwhile, for those skilled in the art, the specific implementation manners and application ranges will be changed according to the idea of the present application. In conclusion, the content of the specification should not be understood as a limitation of the present application.

Claims

1. A Python-Fish-based double-line shield tunnel parameterized automatic simulation method, characterized in that, The method comprises: The geometric parameters of the target double-line shield tunnel are imported into a numerical simulation grid generator, an automatic modeling Python script is run, a parameterized double-line shield tunnel geometric model is established, three-dimensional grid division is performed, and a grid model file is exported; In FLAC3D, the grid model file is imported by using a Fish command; According to the target engineering working condition, double-line shield tunnel calculation simulation parameters are assigned, and the double-line shield tunnel calculation simulation parameters comprise: constitutive model parameters, shield construction parameters, support parameters, and double-line shield tunnel lag distance; Based on the double-line shield tunnel calculation simulation parameters and the grid model file, double-line shield tunnel excavation cycle simulation is performed to generate a simulation result file; The simulation result file is imported into FLAC3D to view and analyze the stress and deformation of the soil body and the structure.

2. The Python-Fish based double-line shield tunnel parameterized automatic simulation method according to claim 1, characterized in that, In FLAC3D, the grid model file is imported by using a Fish command, which comprises: The grid model file is imported into FLAC3D, the geometric model is grouped, and the first tunnel T1, the second tunnel T2, and the surrounding soil S are divided into three groups.

3. The Python-Fish based double-line shield tunnel parameterized automatic simulation method according to claim 1, characterized in that, The geometric parameters of the target double-line shield tunnel comprise: model size L, width W, height H; tunnel diameter D, tunnel burial depth C, double-line tunnel center distance B, double-line tunnel offset angle θ, and grid size Δ; L, W, and H are used to limit the external boundary of the double-line shield tunnel geometric model; D and C are used to limit the diameter and vault burial depth of the double-line shield tunnel geometric model, respectively; B and θ limit the positions of the two tunnels in the double-line shield tunnel geometric model, and the center coordinates of the two tunnels are determined by the following formula: ； ； ； ； Δ is used to control the size of the three-dimensional grid division.

4. The Python-Fish based double shield tunnel parameterized automatic simulation method according to any one of claims 1-3, characterized in that, The constitutive model of the soil body adopts a small-strain plastic hardening model; the shield construction parameters comprise face pressure, grouting pressure, and jack thrust; the support parameters comprise the stiffness of the lining segment and the elastic modulus of the grouting; and the double-line shield tunnel lag distance comprises the distance between the faces of the first tunnel and the second tunnel.

5. The Python-Fish based double-line shield tunnel parameterized automatic simulation method according to claim 4, characterized in that, Based on the double-line shield tunnel calculation simulation parameters and the grid model file, double-line shield tunnel cycle excavation simulation is performed, which comprises: According to the constitutive model parameters, automatic assignment of the physical and mechanical parameters of the three-dimensional numerical model is completed, and the initial stress balance of the model is completed; according to the shield construction parameters, the support parameters, and the double-line shield tunnel lag distance, automatic double-line shield tunnel excavation cycle simulation is performed to generate a simulation result file.

6. The Python-Fish based double-line shield tunnel parameterized automatic simulation method according to claim 5, characterized in that, According to the shield construction parameters, the support parameters, and the double-line shield tunnel lag distance, automatic double-line shield tunnel excavation cycle simulation is performed, which comprises: The face pressure is determined as a trapezoidal distribution of horizontal stress acting vertically on the face, which increases linearly along the depth; The shield machine is simulated by using a shell element, the stiffness decreases from front to back to simulate the taper effect, and the self-weight pressure is the vertical stress acting on the length and 90° cross-section range of the shield machine; The jack thrust is applied to the latest spliced ring segment as a uniformly distributed nodal concentrated force along the segment circumference; The grouting pressure is applied within the width range of the latest spliced ring segment as a normal pressure applied to the surrounding rock and the segment, which increases linearly along the depth; The trolley pressure is determined as a vertical pressure acting on two parallel lines of the segment; The liner element is used to simulate the lining segment, and the joint connection between the segments and the interaction between the segments and the grouting are considered; The linear elastic solid element is used to simulate the grouting body, the hardening process of the grouting body is considered, the first ring after the shield tail is dragged out, the grouting body is fresh grout, and low strength parameters are given, and the fresh grout is converted into hardened grouting body after two excavation steps, and high strength parameters are given.

7. A device for automatic simulation of double-line shield tunnel parameterization based on Python-Fish, characterized in that, The device comprises: A model establishing module is configured to import geometric parameters of a target double-line shield tunnel into a numerical simulation grid generator, run an automatic modeling Python script, establish a parameterized double-line shield tunnel geometric model, perform three-dimensional grid division, and export a grid model file; An importing module is configured to import the grid model file into FLAC3D by using a Fish command; A value assigning module is configured to assign double-line shield tunnel calculation simulation parameters according to a target engineering working condition, wherein the double-line shield tunnel calculation simulation parameters include constitutive model parameters, shield construction parameters, support parameters, and double-line shield tunnel lag distance; A simulation module is configured to perform double-line shield tunnel cyclic excavation simulation based on the double-line shield tunnel calculation simulation parameters and the grid model file, and generate a simulation result file; A display module is configured to import the simulation result file into FLAC3D, and view and analyze the stress and deformation of the soil body and the structure.

8. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor executes the computer program to realize the Python-Fish-based double-line shield tunnel parameterized automatic simulation method in any one of claims 1-6.

9. A computer-readable storage medium having stored thereon a computer program, characterized in that, The computer program is executed by the processor to realize the Python-Fish-based double-line shield tunnel parameterized automatic simulation method in any one of claims 1-6.

10. A computer program product comprising computer programs / instructions, characterized in that, The computer program / instructions are executed by the processor to realize the steps in the Python-Fish-based double-line shield tunnel parameterized automatic simulation method in any one of claims 1-6.

Images