Formation loss simulation method, device, equipment and storage medium

Abstract

The present invention provides a formation loss simulation method, device, equipment and storage medium, which relate to the technical field of internal force analysis of underground structure engineering. The formation loss simulation method includes: performing an initial stress state analysis on the finite element model of the underground structure through the start command of the simulation system; performing an excavation analysis on the finite element model after the initial stress state analysis through the removal command of the simulation system to obtain the reaction force at the top of the excavated soil; gradually reducing the reaction force to zero, and outputting the deformation amount of the top of the rigid body inside the excavated soil; when the deformation amount of the top of the rigid body inside the excavated soil exceeds a preset deformation amount, obtaining the formation loss; the preset deformation amount is the product of a predetermined formation loss rate and a station height. The formation loss simulation method, device, equipment and storage medium provided by the present invention can accurately simulate the process of passive settlement of soil into arches, generate accurate formation losses, and thus improve the accuracy of underground dynamic structure analysis.

Description

Technical Field

[0001] The present invention relates to the technical field of internal force analysis of underground structural engineering, and in particular to a formation loss simulation method, device, equipment and storage medium. Background Art

[0002] There are currently three ways to generate stratum loss: one is the ground stress release rate control method: by adjusting the stress release rate of the tunnel excavation boundary, the convergence of the tunnel over-excavation boundary and the stress release process are simulated. One is the radial displacement control method: similar to the ground stress release rate control method, it actively simulates the convergence process of the over-excavation stratum by setting a specific radial displacement at the tunnel excavation boundary. Another is the equal-generation layer method: setting an equal-thickness, homogeneous weak layer between the tunnel and the surrounding rock to simulate the excavation gap, thereby realizing the simulation of tunnel over-excavation, in which the thickness and material parameters of the equal-generation layer need to be finely calibrated according to the actual settlement conditions.

[0003] However, the current way of generating stratum loss is to actively control the boundary of excavated soil to shrink it, which may exceed the shrinkage parameters required for the soil to form an arch. It does not passively simulate the process of soil settlement and arching, but is actively set. This will affect the initial stress state of the structure and lead to inaccurate simulation of the response of the underground structure. Summary of the invention

[0004] The present invention provides a formation loss simulation method, device, equipment and storage medium, which are used to solve the technical problem of low accuracy of underground dynamic structure analysis in the prior art.

[0005] In a first aspect, the present invention provides a formation loss simulation method, comprising the following steps.

[0006] An initial stress state analysis of the finite element model of the underground structure is performed by means of a start command of the simulation system;

[0007] The excavation analysis is performed on the finite element model after the initial stress state analysis by removing the command of the simulation system to obtain the reaction force on the top of the excavated soil;

[0008] gradually reducing the reaction force to zero, and outputting the deformation of the top of the rigid body inside the excavated soil;

[0009] When the deformation amount of the top of the rigid body inside the excavated soil exceeds the preset deformation amount, the formation loss is obtained; the preset deformation amount is the product of the predetermined formation loss rate and the station height.

[0010] In some embodiments, before obtaining the formation loss, the method further includes:

[0011] Activate the rigid body unit inside the excavated soil; wherein the rigid body unit has the same weight as the excavated soil; the four side surfaces and the bottom surface of the rigid body unit are connected through hard contact; and the top surface normal of the rigid body unit is connected using the GAP unit.

[0012] In some embodiments, after obtaining the formation loss, the method further comprises:

[0013] removing the rigid body unit and releasing the contacts on the six faces of the rigid body unit;

[0014] The initial stress state analysis is ended by an end command, and the rigid body units inside the excavated soil and the output deformation are reset.

[0015] In some embodiments, the performing of an initial stress state analysis on the finite element model of the underground structure by using a start command of the simulation system includes:

[0016] Performing elastic gravity analysis on the finite element model through a start command of a simulation system, and performing linear gravity analysis on the finite element model;

[0017] The soil in the elastic gravity analysis process is linear elastic; the soil in the linear gravity analysis process is nonlinear.

[0018] In some embodiments, before performing an initial stress state analysis on the finite element model of the underground structure through a start command of the simulation system, the method further includes:

[0019] Horizontal constraints are imposed on the four side boundaries of the finite element model, and fixed constraints are imposed on the bottom boundary of the finite element model; the top boundary of the finite element model is a free surface.

[0020] In some embodiments, after obtaining the formation loss, the method further comprises:

[0021] Removing constraints on four side boundaries and a bottom boundary of the finite element model;

[0022] Activate the relevant elements and nodes of the absorbing layer, and apply fixed constraints with six degrees of freedom to the four side boundaries and the bottom boundary of the absorbing layer.

[0023] In some embodiments, after obtaining the formation loss, the method further comprises:

[0024] Apply gravity to stations and perform station gravity analysis;

[0025] The seismic wave time history is converted into nodal forces and regional seismic wave simulation and structural dynamic analysis are performed based on the nodal forces.

[0026] In a second aspect, the present invention provides a formation loss simulation device, comprising the following modules.

[0027] A first analysis module is used to perform an initial stress state analysis on a finite element model of the underground structure through a start command of a simulation system;

[0028] The second analysis module is used to perform excavation analysis on the finite element model after the initial stress state analysis by removing commands of the simulation system to obtain the reaction force at the top of the excavated hole soil;

[0029] An output module, used to gradually reduce the reaction force to zero and output the deformation of the top of the rigid body inside the excavated soil;

[0030] The acquisition module is used to obtain the stratum loss when the deformation of the top of the rigid body inside the excavated soil exceeds a preset deformation; the preset deformation is the product of the predetermined stratum loss rate and the station height.

[0031] In a third aspect, the present invention provides an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the formation loss simulation method as described in the first aspect above is implemented.

[0032] In a fourth aspect, the present invention further provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the formation loss simulation method as described in the first aspect above.

[0033] In a fifth aspect, the present invention further provides a computer program product, comprising a computer program, which, when executed by a processor, implements the formation loss simulation method as described in the first aspect above.

[0034] The stratum loss simulation method, device, equipment and storage medium provided by the present invention perform an initial stress state analysis on the finite element model of the underground structure through the start command of the simulation system, perform an excavation analysis on the finite element model after the initial stress state analysis through the removal command of the simulation system, obtain the reaction force at the top of the excavated soil, gradually reduce the reaction force to zero, output the deformation of the top of the rigid body inside the excavated soil, and form a stratum loss after the deformation of the top of the rigid body inside the excavated soil exceeds the preset deformation, wherein the preset deformation is the product of the predetermined stratum loss rate and the station height. The stratum loss simulation method, device, equipment and storage medium provided by the present invention can accurately simulate the process of passive settlement of soil to form an arch, generate accurate stratum loss, rather than forcibly shrinking the excavation boundary or simply releasing the ground stress around the hole, avoid the set stratum loss rate exceeding the stratum loss rate required for this type of soil to form a complete soil arch, thereby improving the accuracy of underground dynamic structure analysis. BRIEF DESCRIPTION OF THE DRAWINGS

[0035] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0036] Figure 1 It is a schematic flow chart of the formation loss simulation method provided by the present invention.

[0037] Figure 2 It is a schematic diagram of the functional relationship between the reaction force provided by the present invention and the deformation amount of the top of the rigid body inside the excavated soil.

[0038] Figure 3 It is a schematic diagram of the structure of the free field boundary unit provided by the present invention.

[0039] Figure 4 It is a schematic diagram of the process of underground dynamic structure analysis provided by the present invention.

[0040] Figure 5 It is a schematic diagram of changes in soil stress field for formation loss analysis provided by the present invention.

[0041] Figure 6 It is a structural schematic diagram of the formation loss simulation device provided by the present invention.

[0042] Figure 7 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0043] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0044] Figure 1 Schematic diagram of the process of the formation loss simulation method provided by the present invention, such as Figure 1 As shown, the present invention provides a formation loss simulation method. The method comprises the following steps:

[0045] Step 101: Perform an initial stress state analysis on the finite element model of the underground structure through a start command of the simulation system.

[0046] Specifically, the simulation system may use OpenSees, which is an open system for earthquake engineering simulation. The start command may be an "InitialStateAnalysis on" command provided in OpenSees.

[0047] An initial stress state analysis is performed on the finite element model of the underground structure through a start command of the simulation system, for example, a gravity load is applied to the finite element model of the underground structure to analyze the stress of the finite element model.

[0048] Step 102: Perform excavation analysis on the finite element model after the initial stress state analysis by using the removal command of the simulation system to obtain the reaction force at the top of the excavated soil.

[0049] Specifically, the removal command of the simulation system may be a "remove" command provided by OpenSees. According to the removal command, the excavation process is performed on the finite element model, and the reaction force (ReactionForce) at the top of the excavated soil is obtained through excavation analysis.

[0050] Step 103: gradually reduce the reaction force to zero, and output the deformation of the top of the rigid body inside the excavated soil.

[0051] Specifically, after completing the initial stress state analysis and excavation analysis, the stratum loss analysis is performed to gradually reduce the reaction force to zero, so that the units of the internal excavated soil and the soil outside the hole invade each other and produce deformation, and the deformation of the top of the rigid body inside the excavated soil is output.

[0052] Step 104: when the deformation of the top of the rigid body inside the excavated soil exceeds a preset deformation, obtain the stratum loss; the preset deformation is the product of a predetermined stratum loss rate and a station height.

[0053] Specifically, after the reaction force is reduced to zero, if the deformation amount of the top of the rigid body inside the excavated soil exceeds the preset deformation amount, stratum loss is formed.

[0054] Among them, the preset deformation amount (i.e., GAP value) is equal to the product of the predetermined ground loss rate and the station height, which is publicly expressed as GAP=Ground loss ratio×H, where Ground loss ratio is the specified ground loss rate and H is the station height.

[0055] The stratum loss simulation method provided in the embodiment of the present application performs an initial stress state analysis (i.e., soil gravity analysis) on the finite element model of the underground structure through the start command of the simulation system, performs excavation analysis on the finite element model after the initial stress state analysis through the removal command of the simulation system, obtains the reaction force at the top of the excavated soil, and outputs the deformation of the top of the rigid body inside the excavated soil by gradually reducing the reaction force to zero. After the deformation of the top of the rigid body inside the excavated soil exceeds the preset deformation, stratum loss is formed, wherein the preset deformation is the product of the predetermined stratum loss rate and the station height. It can accurately simulate the process of passive settlement of soil to form an arch, generate accurate stratum loss, rather than forcibly shrinking the excavation boundary or simply releasing the ground stress around the hole, avoid the set stratum loss rate exceeding the stratum loss rate required for this type of soil to form a complete soil arch, and thus improve the accuracy of underground dynamic structure analysis.

[0056] In some embodiments, before performing an initial stress state analysis on the finite element model of the underground structure through a start command of the simulation system, the method further includes:

[0057] Horizontal constraints are imposed on the four side boundaries of the finite element model, and fixed constraints are imposed on the bottom boundary of the finite element model; the top boundary of the finite element model is a free surface.

[0058] Specifically, before performing the initial stress state analysis on the finite element model, the boundaries of each surface of the finite element model are first constrained:

[0059] Apply fixed constraints to the bottom boundary of the finite element model, that is, constrain the translational and rotational degrees of freedom in three directions; apply horizontal constraints to the four lateral boundaries, that is, horizontally constrain the degrees of freedom in the X direction on the left and right sides, and horizontally constrain the degrees of freedom in the Y direction on the front and back sides; set the top boundary of the finite element model as a free surface.

[0060] In some embodiments, the performing of an initial stress state analysis on the finite element model of the underground structure by using a start command of the simulation system includes:

[0061] Performing elastic gravity analysis on the finite element model through a start command of a simulation system, and performing linear gravity analysis on the finite element model;

[0062] The soil in the elastic gravity analysis process is linear elastic; the soil in the linear gravity analysis process is nonlinear.

[0063] Specifically, the initial stress state analysis of the finite element model of the underground structure is performed through the start command of the simulation system. First, a gravity load is applied to the model for elastic gravity analysis, and the stage-0 stage of the "PIMY" material is activated. At this stage, the soil is linear elastic (the bulk modulus is the small strain bulk modulus, and the shear modulus is the small strain shear modulus). Secondly, a nonlinear gravity analysis is performed on the model, and the stage-1 stage of the "PIMY" material is activated. At this stage, the soil is nonlinear.

[0064] The "PIMY" material refers to the Pressure Independ Multi Yield material, which is an elastoplastic material that exhibits plasticity only in the deviatoric stress-strain response. The volumetric stress-strain response is linear elastic. This material is used to simulate the monotonic or cyclic response of materials whose shear behavior is insensitive to changes in constraints.

[0065] The command used in the update / activation stage is "updateMaterialStage". When the "PIMY" material is in stage-0, the material is a linear elastic material; when it is updated to stage-1, the nonlinearity of the material is activated.

[0066] After the nonlinearity of the "PIMY" material is activated, there will be a modulus reduction curve. The small strain modulus corresponds to the maximum modulus of the initial input of the material, that is, the maximum elastic modulus and the maximum shear modulus.

[0067] The formation loss simulation method provided in the embodiment of the present application constrains the boundaries of each surface of the finite element model and performs an initial stress state analysis on the finite element model through two steps of elastic gravity analysis and linear gravity analysis, thereby achieving accurate initial stress state analysis, which is conducive to accurately simulating formation loss and thus improving the accuracy of underground dynamic structure analysis.

[0068] In some embodiments, performing excavation analysis on the finite element model by removing commands from a simulation system to obtain a reaction force at the top of the excavated soil includes:

[0069] Removing relevant elements and relevant nodes of the excavated soil from the finite element model through a removal command of the simulation system to obtain an excavation analysis result;

[0070] The reaction force at the top of the excavation soil is obtained from the excavation analysis result.

[0071] Specifically, after the initial stress state analysis is performed, the finite element model is subjected to excavation analysis through the simulation system, such as removing the excavation soil related elements and nodes from the model according to the "remove" command provided by OpenSees, and extracting the ReactionForce (i.e., reaction force) at the top of the excavation soil from the analysis results.

[0072] In some embodiments, before obtaining the formation loss, the method further includes:

[0073] Activate the rigid body unit inside the excavated soil; wherein the rigid body unit has the same weight as the excavated soil; the four side surfaces and the bottom surface of the rigid body unit are connected through hard contact; and the top surface normal of the rigid body unit is connected using the GAP unit.

[0074] Specifically, after the initial stress state analysis and excavation analysis, the ground loss analysis is performed: first, the rigid body elements inside the excavated soil are activated.

[0075] The rigid body unit has the same weight as the excavated soil.

[0076] The four sides and bottom of the rigid body unit are connected by hard contact, which can be achieved by the ZeroLengthContactASDimplex element (the connection unit is compressed but not pulled, and there is friction under compression), and the friction coefficient is 0.

[0077] The top surface normal of the rigid body element is connected using the GAP element (also known as the gap element). Figure 2 Schematic diagram of the functional relationship between the reaction force provided by the present invention and the deformation of the top of the rigid body inside the excavated soil, such as Figure 2 As shown in the figure, when the top surface normal of the rigid unit is connected with the GAP unit, it is realized by the 0 length unit, and no unit is set in the tangent direction, and E is the stiffness. When analyzing the formation loss, the reaction force (ReactionForce) is gradually reduced to 0, so that the units of the internal excavation soil and the soil outside the hole invade each other, and contact each other after exceeding the GAP value, thereby forming formation loss, where the GAP value is calculated by the specified formation loss rate.

[0078] The formation loss simulation method provided in the embodiment of the present application improves the accuracy of formation loss simulation by using GAP units to connect the top surface normals of rigid units and combining preset GAP values ​​to form formation loss, thereby improving the accuracy of underground dynamic structure analysis.

[0079] In some embodiments, after obtaining the formation loss, the method further comprises:

[0080] removing the rigid body unit and releasing the contacts on the six faces of the rigid body unit;

[0081] The initial stress state analysis is ended by an end command, and the rigid body units inside the excavated soil and the output deformation are reset.

[0082] Specifically, after obtaining the formation loss, a ground stress balance operation is performed: the rigid body unit is removed, and the contact on the six faces of the rigid body unit is released, and the initial stress state analysis is ended using an end command (such as the "InitialStateAnalysis off" command), and the rigid body unit inside the excavated soil and the output deformation are reset.

[0083] In some embodiments, after obtaining the formation loss, the method further comprises:

[0084] Removing constraints on four side boundaries and a bottom boundary of the finite element model;

[0085] Activate the relevant elements and nodes of the absorbing layer, and apply fixed constraints with six degrees of freedom to the four side boundaries and the bottom boundary of the absorbing layer.

[0086] Specifically, the ground stress balance operation also includes: removing all constraints on the four side boundaries and the bottom boundary of the input layer finite element model, activating the relevant elements and nodes of the absorption layer, and applying fixed constraints of six degrees of freedom on the four side boundaries and the bottom boundary of the absorption layer.

[0087] Among them, the input layer is used to convert seismic waves into equivalent nodal force time history data and input it into the model to restore the seismic wave field in the internal soil domain. The absorption layer (absorption boundary) is used to absorb seismic waves propagating outward from the internal soil domain to prevent the influence of seismic wave reflection on the wave field in the internal soil domain. Figure 3 is a schematic diagram of the structure of the free field boundary unit provided by the present invention, such as Figure 3 As shown, nodes 1 to 4 form an absorbing unit, where nodes 1 and 4 are connected to nodes 2 and 3 through springs and dampers ( Figure 3 Nodes 1 and 4 are not coupled to nodes 2 and 3 because the seismic waves are only transmitted from the input layer to the inner domain, preventing the seismic waves from propagating to the outer domain.

[0088] The formation loss simulation method provided in the embodiment of the present application performs a ground stress balance process after obtaining the formation loss, thereby providing an accurate initial stress state for subsequent underground structure dynamic analysis.

[0089] In some embodiments, after obtaining the formation loss, the method further comprises:

[0090] Apply gravity to stations and perform station gravity analysis;

[0091] The seismic wave time history is converted into nodal forces and regional seismic wave simulation and structural dynamic analysis are performed based on the nodal forces.

[0092] Specifically, after performing ground stress balance, gravity is applied to the station and station gravity analysis is performed. Then, the seismic wave time history is converted into nodal forces and based on the nodal forces, regional seismic wave simulation and structural dynamic analysis are performed.

[0093] When performing gravity analysis of the station, first activate the station and the related elements and nodes of the side contact. The contact uses a 2-node zero-length element "ZeroLengthContactASDimplex" to simulate the potential friction and separation phenomenon of the stratum-tunnel interface. This element adopts a penalty formula and follows Coulomb friction behavior. The friction coefficient μ is the control parameter, which can be set to μ=0.4. Then apply gravity to the station and perform analysis.

[0094] Before conducting earthquake dynamic analysis, a pre-dynamic analysis of 1000s is first performed to eliminate the initial numerical oscillation. Then, the earthquake wave time history is converted into node forces and input into the Input layer to simulate regional earthquake waves and perform dynamic analysis of the structure.

[0095] The formation loss simulation method provided in the embodiment of the present application improves the accuracy of underground dynamic structure analysis by using 2-node zero-length elements to perform station gravity analysis and eliminating initial numerical oscillations before performing dynamic analysis.

[0096] Figure 4 It is a schematic diagram of the process of underground dynamic structure analysis provided by the present invention, such as Figure 4 As shown, the specific steps of underground dynamic structure analysis include: (a) soil gravity analysis, (b) excavation analysis, (c) stratum loss analysis, (d) ground stress balance, (e) station gravity analysis, and (f) earthquake dynamic analysis.

[0097] First, (a) soil gravity analysis is performed to obtain the initial stress field of the soil. On this basis, (b) excavation analysis of the soil is performed to remove the excavated soil. Then, (c) stratum loss analysis is performed to deflect the stress around the excavated soil. Figure 5 Schematic diagram of the change of soil stress field in the stratum loss analysis provided by the present invention, such as Figure 5 As shown, Figure 5 (1) is the schematic diagram of the initial soil stress field (in MPa). Figure 5(2) is a schematic diagram of the soil stress field that produces the soil arch unloading effect during the stratum loss analysis; (d) ground stress balance is performed, that is, the soil stress field at this moment is maintained and the deformation of all units is restored to 0; the station and its contact units are activated, and the (e) gravity analysis of the station is performed on the basis of the soil arch effect produced by the soil; at this point, the initial stress of the soil and station considering the stratum loss is ready, and the prepared seismic wave is input into the input layer to perform the (f) seismic dynamic analysis of the station-soil.

[0098] Figure 6 Schematic diagram of the structure of the formation loss simulation device provided by the present invention. Figure 6 As shown, the present invention provides a formation loss simulation device, including a first analysis module 601, a second analysis module 602, an output module 603 and an acquisition module 604.

[0099] The first analysis module 601 is used to perform an initial stress state analysis on the finite element model of the underground structure through a start command of the simulation system.

[0100] The second analysis module 602 is used to perform excavation analysis on the finite element model after the initial stress state analysis by removing the command of the simulation system, and obtain the reaction force on the top of the excavated soil;

[0101] The output module 603 is used to gradually reduce the reaction force to zero and output the deformation of the top of the rigid body inside the excavated soil;

[0102] The acquisition module 604 is used to acquire the stratum loss when the deformation of the top of the rigid body inside the excavated soil exceeds a preset deformation; the preset deformation is the product of the predetermined stratum loss rate and the station height.

[0103] In some embodiments, the second analysis module includes:

[0104] A first acquisition unit is used to remove the relevant elements and relevant nodes of the excavated soil from the finite element model through a removal command of the simulation system to obtain an excavation analysis result;

[0105] The second acquisition unit is used to acquire the reaction force of the top of the excavation hole soil from the excavation analysis result.

[0106] In some embodiments, it also includes:

[0107] The first activation module is used to activate the rigid body unit inside the excavated soil; wherein the rigid body unit has the same weight as the excavated soil; the four side surfaces and the bottom surface of the rigid body unit are connected through hard contact; and the top surface normal of the rigid body unit is connected using the GAP unit.

[0108] In some embodiments, it also includes:

[0109] A first removal module, used to remove the rigid body unit and release the contact on the six faces of the rigid body unit;

[0110] The reset module is used to end the initial stress state analysis through an end command, and reset the rigid body units inside the excavated soil and the output deformation.

[0111] In some embodiments, the first analysis module includes:

[0112] an analysis unit, configured to perform elastic gravity analysis on the finite element model and linear gravity analysis on the finite element model through a start command of a simulation system;

[0113] The soil in the elastic gravity analysis process is linear elastic; the soil in the linear gravity analysis process is nonlinear.

[0114] In some embodiments, it also includes:

[0115] The constraint applying module is used to apply horizontal constraints to the four side boundaries of the finite element model and to apply fixed constraints to the bottom boundary of the finite element model; the top boundary of the finite element model is a free surface.

[0116] In some embodiments, it also includes:

[0117] A second removal module, used to remove the constraints of the four side boundaries and the bottom boundary of the finite element model;

[0118] The second activation module is used to activate the relevant elements and nodes of the absorption layer and impose fixed constraints of six degrees of freedom on the four side boundaries and the bottom boundary of the absorption layer.

[0119] In some embodiments, it also includes:

[0120] The third analysis module is used to apply gravity to the station and perform station gravity analysis;

[0121] The fourth analysis module is used to convert the seismic wave time history into nodal forces and simulate regional seismic waves and perform dynamic analysis of the structure based on the nodal forces.

[0122] Specifically, the above-mentioned formation loss simulation device provided by the present invention can implement all the method steps implemented by the above-mentioned formation loss simulation method embodiment, and can achieve the same technical effect. The parts and beneficial effects of this embodiment that are the same as the method embodiment will not be described in detail here.

[0123] It should be noted that the division of units / modules in the above-mentioned embodiments of the present invention is schematic and is only a logical function division. There may be other division methods in actual implementation. In addition, the functional units in the various embodiments of the present application may be integrated into a processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated units may be implemented in the form of hardware or in the form of software functional units.

[0124] Figure 7 is a schematic diagram of the structure of the electronic device provided by the present invention, such as Figure 7 As shown, the electronic device may include: a processor 701, a communications interface 702, a memory 703 and a communication bus 704, wherein the processor 701, the communications interface 702 and the memory 703 communicate with each other via the communication bus 704. The processor 701 may call the logic instructions in the memory 703 to execute the formation loss simulation method, which includes:

[0125] An initial stress state analysis of the finite element model of the underground structure is performed by means of a start command of the simulation system;

[0126] The excavation analysis is performed on the finite element model after the initial stress state analysis by removing the command of the simulation system to obtain the reaction force on the top of the excavated soil;

[0127] gradually reducing the reaction force to zero, and outputting the deformation of the top of the rigid body inside the excavated soil;

[0128] When the deformation amount of the top of the rigid body inside the excavated soil exceeds the preset deformation amount, the formation loss is obtained; the preset deformation amount is the product of the predetermined formation loss rate and the station height.

[0129] In some embodiments, before obtaining the formation loss, the method further includes:

[0130] Activate the rigid body unit inside the excavated soil; wherein the rigid body unit has the same weight as the excavated soil; the four side surfaces and the bottom surface of the rigid body unit are connected through hard contact; and the top surface normal of the rigid body unit is connected using the GAP unit.

[0131] In some embodiments, after obtaining the formation loss, the method further comprises:

[0132] removing the rigid body unit and releasing the contacts on the six faces of the rigid body unit;

[0133] The initial stress state analysis is ended by an end command, and the rigid body units inside the excavated soil and the output deformation are reset.

[0134] In some embodiments, the performing of an initial stress state analysis on the finite element model of the underground structure by using a start command of the simulation system includes:

[0135] Performing elastic gravity analysis on the finite element model through a start command of a simulation system, and performing linear gravity analysis on the finite element model;

[0136] The soil in the elastic gravity analysis process is linear elastic; the soil in the linear gravity analysis process is nonlinear.

[0137] In some embodiments, before performing an initial stress state analysis on the finite element model of the underground structure through a start command of the simulation system, the method further includes:

[0138] Horizontal constraints are imposed on the four side boundaries of the finite element model, and fixed constraints are imposed on the bottom boundary of the finite element model; the top boundary of the finite element model is a free surface.

[0139] In some embodiments, after obtaining the formation loss, the method further comprises:

[0140] Removing constraints on four side boundaries and a bottom boundary of the finite element model;

[0141] Activate the relevant elements and nodes of the absorbing layer, and apply fixed constraints with six degrees of freedom to the four side boundaries and the bottom boundary of the absorbing layer.

[0142] In some embodiments, after obtaining the formation loss, the method further comprises:

[0143] Apply gravity to stations and perform station gravity analysis;

[0144] The seismic wave time history is converted into nodal forces and regional seismic wave simulation and structural dynamic analysis are performed based on the nodal forces.

[0145] Specifically, the processor 701 may be a central processing unit (CPU), an application specific integrated circuit (ASIC), a field programmable gate array (FPGA) or a complex programmable logic device (CPLD), and the processor may also adopt a multi-core architecture.

[0146] The logic instructions in the memory 703 can be implemented in the form of software functional units and can be stored in a processor-readable storage medium when sold or used as an independent product. Based on this understanding, the technical solution of the present application is essentially or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product. The computer software product is stored in a storage medium, including several instructions to enable a computer device (which can be a personal computer, server, or network device, etc.) or a processor to perform all or part of the steps of the method described in each embodiment of the present application. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM), random access memory (RAM), disk or optical disk and other media that can store program code.

[0147] In some embodiments, a computer program product is further provided. The computer program product includes a computer program. The computer program can be stored in a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the formation loss simulation method provided by each of the above method embodiments. The method includes:

[0148] An initial stress state analysis of the finite element model of the underground structure is performed by means of a start command of the simulation system;

[0149] The excavation analysis is performed on the finite element model after the initial stress state analysis by removing the command of the simulation system to obtain the reaction force on the top of the excavated soil;

[0150] gradually reducing the reaction force to zero, and outputting the deformation of the top of the rigid body inside the excavated soil;

[0151] When the deformation amount of the top of the rigid body inside the excavated soil exceeds the preset deformation amount, the formation loss is obtained; the preset deformation amount is the product of the predetermined formation loss rate and the station height.

[0152] Specifically, the above-mentioned computer program product provided in the embodiment of the present application can implement all the method steps implemented by the above-mentioned method embodiments, and can achieve the same technical effect. The parts and beneficial effects of this embodiment that are the same as the method embodiment will not be described in detail here.

[0153] In some embodiments, a computer-readable storage medium is further provided, wherein the computer-readable storage medium stores a computer program, wherein the computer program is used to enable a computer to execute the formation loss simulation method provided by each of the above method embodiments, the method comprising:

[0154] An initial stress state analysis of the finite element model of the underground structure is performed by means of a start command of the simulation system;

[0155] The excavation analysis is performed on the finite element model after the initial stress state analysis by removing the command of the simulation system to obtain the reaction force on the top of the excavated soil;

[0156] gradually reducing the reaction force to zero, and outputting the deformation of the top of the rigid body inside the excavated soil;

[0157] When the deformation amount of the top of the rigid body inside the excavated soil exceeds the preset deformation amount, the formation loss is obtained; the preset deformation amount is the product of the predetermined formation loss rate and the station height.

[0158] Specifically, the above-mentioned computer-readable storage medium provided by the present invention can implement all the method steps implemented by the above-mentioned method embodiments, and can achieve the same technical effects. The parts and beneficial effects that are the same as the method embodiments in this embodiment will not be described in detail here.

[0159] It should be noted that the computer-readable storage medium can be any available medium or data storage device that can be accessed by the processor, including but not limited to magnetic storage (such as floppy disks, hard disks, magnetic tapes, magneto-optical disks (MO), etc.), optical storage (such as CD, DVD, BD, HVD, etc.), and semiconductor storage (such as ROM, EPROM, EEPROM, non-volatile memory (NAND FLASH), solid-state drive (SSD)), etc.

[0160] It should also be noted that the terms "first", "second", etc. in the present invention are used to distinguish similar objects, and are not used to describe a specific order or sequence. It should be understood that the terms used in this way can be interchangeable under appropriate circumstances, so that the embodiments of the present invention can be implemented in an order other than those illustrated or described herein, and the objects distinguished by "first" and "second" are usually of the same type, and the number of objects is not limited. For example, the first object can be one or more.

[0161] In the present invention, "determine B based on A" means that the factor A should be considered when determining B. It is not limited to "B can be determined based on A alone", but should also include: "determine B based on A and C", "determine B based on A, C and E", "determine C based on A, and further determine B based on C", etc. In addition, it can also include taking A as a condition for determining B, for example, "when A meets the first condition, use the first method to determine B"; for another example, "when A meets the second condition, determine B", etc.; for another example, "when A meets the third condition, determine B based on the first parameter", etc. Of course, it can also be a condition that takes A as a factor for determining B, for example, "when A meets the first condition, use the first method to determine C, and further determine B based on C", etc.

[0162] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present application may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage and optical storage, etc.) containing computer-usable program codes.

[0163] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer executable instructions. These computer executable instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0164] These processor executable instructions may also be stored in a processor readable memory that can direct a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the processor readable memory produce an article of manufacture including an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0165] These processor-executable instructions may also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for executing on the computer or other programmable device to implement the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0166] Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include these modifications and variations.

Claims

1. A formation loss simulation method, characterized in that: include: An initial stress state analysis of the finite element model of the underground structure is performed by means of a start command of the simulation system; The excavation analysis is performed on the finite element model after the initial stress state analysis by removing the command of the simulation system to obtain the reaction force on the top of the excavated soil; gradually reducing the reaction force to zero, and outputting the deformation of the top of the rigid body inside the excavated soil; When the deformation amount of the top of the rigid body inside the excavated soil exceeds the preset deformation amount, the formation loss is obtained; the preset deformation amount is the product of the predetermined formation loss rate and the station height.

2. The formation loss simulation method according to claim 1, characterized in that: Before obtaining the formation loss, the method further includes: Activate the rigid body unit inside the excavated soil; wherein the rigid body unit and the excavated soil have the same weight; the four side surfaces and the bottom surface of the rigid body unit are connected through hard contact; the top surface normal of the rigid body unit is connected using the GAP unit.

3. The formation loss simulation method according to claim 2, characterized in that: After obtaining the formation loss, the method further includes: Removing the rigid body unit and releasing the contact on the six faces of the rigid body unit; The initial stress state analysis is ended by an end command, and the rigid body units inside the excavated soil and the output deformation are reset.

4. The formation loss simulation method according to claim 1, characterized in that: The initial stress state analysis of the finite element model of the underground structure is performed by using the start command of the simulation system, including: Performing elastic gravity analysis on the finite element model through a start command of a simulation system, and performing linear gravity analysis on the finite element model; The soil in the elastic gravity analysis process is linear elastic; the soil in the linear gravity analysis process is nonlinear.

5. The formation loss simulation method according to claim 1, characterized in that: Before the initial stress state analysis of the finite element model of the underground structure is performed by the start command of the simulation system, the method further includes: Horizontal constraints are imposed on the four side boundaries of the finite element model, and fixed constraints are imposed on the bottom boundary of the finite element model; the top boundary of the finite element model is a free surface.

6. The formation loss simulation method according to claim 5, characterized in that: After obtaining the formation loss, the method further includes: Removing constraints on four side boundaries and a bottom boundary of the finite element model; Activate the relevant elements and nodes of the absorbing layer, and apply fixed constraints with six degrees of freedom to the four side boundaries and the bottom boundary of the absorbing layer.

7. The formation loss simulation method according to claim 1, characterized in that: After obtaining the formation loss, the method further includes: Apply gravity to stations and perform station gravity analysis; The seismic wave time history is converted into nodal forces and regional seismic wave simulation and structural dynamic analysis are performed based on the nodal forces.

8. A formation loss simulation device, characterized in that: include: A first analysis module is used to perform an initial stress state analysis on a finite element model of the underground structure through a start command of a simulation system; The second analysis module is used to perform excavation analysis on the finite element model after the initial stress state analysis by removing commands of the simulation system to obtain the reaction force at the top of the excavated hole soil; An output module, used to gradually reduce the reaction force to zero and output the deformation of the top of the rigid body inside the excavated soil; The acquisition module is used to obtain the stratum loss when the deformation of the top of the rigid body inside the excavated soil exceeds a preset deformation; the preset deformation is the product of the predetermined stratum loss rate and the station height.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and running on the processor, characterized in that: When the processor executes the computer program, the formation loss simulation method according to any one of claims 1 to 7 is implemented.

10. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the formation loss simulation method according to any one of claims 1 to 7 is implemented.

Images