Method, system and storage medium for simulating seismic force response of prefabricated subway stations

Abstract

The present invention provides a method, system and storage medium for simulating the anti-seismic force response of an assembled subway station. The method comprises the following steps: S1, constructing a geometric model of the assembled subway station in a finite element software platform, wherein the geometric model comprises a large-volume concrete part and a joint part, wherein the large-volume concrete part is simulated using a polymer rod model, and the joint part is simulated using an axial shear coupling unit model; S2, establishing a soil unit around the assembled subway station; and connecting it to the structure of the assembled subway station through a common node method to simulate the rigid connection between the soil unit and the assembled subway station; then using the MYS model to simulate the nonlinear behavior of the soil unit under earthquake action, and finally capturing the cumulative plastic deformation of the soil unit through multiple yield surfaces.

Description

Technical Field

[0001] The invention relates to a method, system and storage medium for simulating anti-vibration force response of an assembled subway station. Background Art

[0002] With the rapid development of urban rail transit systems, prefabricated subway stations, as an efficient method for constructing underground infrastructure, have attracted more and more attention. Compared with traditional cast-in-place structures, prefabricated subway stations use the method of splicing prefabricated components, which has the advantages of short construction period, small impact on urban traffic, and high resource utilization. However, the seismic performance of prefabricated structures is significantly different from that of traditional cast-in-place structures, and its seismic design and performance evaluation have always been one of the difficulties in the field of structural engineering. Under the action of earthquakes, the structural response of prefabricated subway stations will be affected by multiple factors, mainly including the nonlinear behavior of mass concrete, the mechanical properties of assembly joints, and soil-structure interaction. These factors are difficult to accurately simulate in traditional analysis methods, so it is particularly important to propose a new simulation method for the seismic force response of prefabricated subway stations. Under the action of earthquakes, the structural response of prefabricated subway stations will be affected by multiple factors, mainly including the nonlinear behavior of mass concrete, the mechanical properties of assembly joints, and soil-structure interaction. These factors are difficult to accurately simulate in traditional analysis methods. Summary of the invention

[0003] The present invention provides a method, system and storage medium for simulating the anti-vibration force response of an assembled subway station, which can effectively solve the above-mentioned problems.

[0004] The present invention is achieved in that:

[0005] The present invention provides a method for simulating the anti-vibration force response of an assembled subway station, comprising the following steps:

[0006] S1, constructing a geometric model of the prefabricated subway station in a finite element software platform, wherein the geometric model includes a large-volume concrete part and a joint part, wherein the large-volume concrete part is simulated using a polymer rod model, and the joint part is simulated using an axial-shear coupling unit model; the simulation of the large-volume concrete part using the polymer rod model specifically includes: in a two-dimensional model, defining that there are four mutually intersecting but independent single-axis rods in the structure of the prefabricated subway station, and in a three-dimensional model, defining that there are nine single-axis rods in the structure of the prefabricated subway station, and the stress contribution of each single-axis rod is summarized into an overall stress state through coordinate transformation;

[0007] S2, establish a soil unit around the prefabricated subway station; and connect it with the structure of the prefabricated subway station through a common node method to simulate the rigid connection between the soil unit and the prefabricated subway station; then use the MYS model to simulate the nonlinear behavior of the soil unit under the action of an earthquake, and finally obtain the seismic performance of the prefabricated subway station.

[0008] The present invention further provides a seismic force response simulation system for a prefabricated subway station, comprising:

[0009] processor;

[0010] A memory for storing processor executable instructions; wherein the processor is configured to call the instructions stored in the memory to execute the above method.

[0011] The present invention further provides a computer-readable storage medium storing a computer program, wherein the computer program can be executed by a processor of a device where the computer-readable storage medium is located to implement the above method.

[0012] The beneficial effects of the present invention are:

[0013] 1. Good real-time calculation: The present invention realizes efficient simulation of the nonlinear behavior of large-volume concrete and complex joints in prefabricated subway stations through the Aggregated Line Model (ALM) and the axial shear coupling unit model. Compared with the traditional three-dimensional concrete model, the ALM model has higher calculation efficiency and is simple to use. It greatly reduces the calculation time and improves the stability of the calculation without reducing the calculation accuracy.

[0014] 2. High calculation accuracy: The aggregate rod model (ALM) in the present invention can accurately simulate the nonlinear behavior of large-volume prefabricated structures under load, including nonlinear hysteresis behavior, stiffness degradation, strength degradation, etc. The axial-shear coupling unit model can effectively consider the mechanical behavior of the joints in the prefabricated subway station structure under the combined action of axial force and shear force, making the bearing capacity of the joints closer to reality and improving the accuracy of seismic analysis of prefabricated structures.

[0015] 3. Wide applicability: The seismic simulation method of the present invention is both applicable and flexible, and is applicable to both prefabricated subway stations and static and dynamic analysis of other large-volume concrete structures. In addition, the parameters of the polymer rod model are easy to adjust, and can be flexibly adapted to other prefabricated structures according to specific engineering requirements. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for use in the embodiments will be briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present invention and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without paying creative work.

[0017] Figure 1 It is a flow chart of a method for simulating the anti-seismic force response of a prefabricated subway station provided in an embodiment of the present invention.

[0018] Figure 2 It is a schematic diagram of a polymer rod model used in a method for simulating the anti-seismic force response of a prefabricated subway station provided in an embodiment of the present invention.

[0019] Figure 3 It is a schematic diagram of calculating the shear capacity of the joint part in the shear direction in the anti-seismic force response simulation method for a prefabricated subway station provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0020] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, rather than all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention. Therefore, the following detailed description of the embodiments of the present invention provided in the drawings is not intended to limit the scope of the invention claimed for protection, but merely represents the selected embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0021] In the description of the present invention, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of the features. In the description of the present invention, the meaning of "plurality" is two or more, unless otherwise clearly and specifically defined.

[0022] Reference Figure 1 As shown, an embodiment of the present invention provides a method for simulating the anti-vibration force response of a prefabricated subway station, comprising the following steps:

[0023] S1, constructing a geometric model of the prefabricated subway station in a finite element software platform, wherein the geometric model includes a mass concrete part and a joint part, wherein the mass concrete part is simulated using a polymer rod model, and the joint part is simulated using an axial shear coupling unit model;

[0024] S2, establish a soil unit around the prefabricated subway station; and connect it with the structure of the prefabricated subway station through a common node method to simulate the rigid connection between the soil unit and the prefabricated subway station; then use the MYS model (Multi-yield Surface Plasticity Model) to simulate the nonlinear behavior of the soil unit under earthquake action, and finally obtain the seismic performance of the prefabricated subway station.

[0025] In step S1, the finite element software platform may be a platform such as OpenSees.

[0026] The mass concrete part is simulated using the polymer rod model specifically including:

[0027] In the two-dimensional model, the structure of the prefabricated subway station is defined as having four mutually intersecting but independent single-axis rods, while in the three-dimensional model, the structure of the prefabricated subway station is defined as consisting of nine single-axis rods, and the stress contribution of each single-axis rod is summarized into an overall stress state through coordinate transformation.

[0028] For details, see Figure 2 As shown, in the two-dimensional model, it is defined that there are four mutually intersecting but independent single-axis rods in the structure of the prefabricated subway station, while in the three-dimensional model, it is defined that there are nine single-axis rods in the structure of the prefabricated subway station, and the stress contribution of each single-axis rod is summarized into the overall stress state through coordinate transformation. The specific calculation process may include:

[0029] In the global coordinate system (the global coordinate system is the default coordinate system), the four intersecting but independent single-axis rods can be split into two first rods with a 45-degree angle to the coordinate system and two second rods that coincide with the coordinate system. For the two first rods with a 45-degree angle to the global coordinate system, the transmitted strain vector is converted to the local coordinate system where the two first rods are located (the local coordinate system is the coordinate of each rod) through coordinate transformation. The specific method is to multiply a transformation matrix T(θ) on the left, where θ is the rotation angle from the global coordinate system to the local coordinate system, and its direction is defined as left rotation as positive:

[0030]

[0031] Among them, ε l ,ε gare the strain vectors in the local coordinate system and the global coordinate system, respectively.

[0032] In this local coordinate system, the normal strains in the two directions calculated by the above formula are transferred into the constitutive relation of the uniaxial rod, and the normal stresses in these two directions can be obtained and form a stress vector. The stress vector can be transformed from the local coordinate system to the global coordinate system by multiplying the stress vector on the left by the transformation matrix in the opposite direction to that in the above formula, that is:

[0033]

[0034] Among them, σ l ,σ g are the normal stresses in the local coordinate system and the global coordinate system (N / mm 2 ).

[0035] At this point, we can get the stress vectors that the two first rods that are 45 degrees to the overall coordinate system contribute to the overall stress state. In addition, the calculation idea of ​​the contribution of the two second rods that coincide with the overall coordinate system to the overall stress state is the same as the calculation idea of ​​the first rod that is 45 degrees to the coordinate system. However, since they coincide with the coordinate system, no coordinate conversion is involved. The transmitted strain is directly input into the constitutive structure of the uniaxial rod to obtain the stress vector that contributes to the overall stress state in this state. By adding the stress vectors obtained in the above two cases, we can get the stress state formula generated by these four rods in two dimensions:

[0036]

[0037] Among them, σ x , σ y , τ xy Respectively represent the normal stress along the x-axis, the normal stress along the y-axis, and the shear stress in the xy direction; σ 1 x , σ 1 y , τ 2 xy and σ 2 x , σ 2 y , τ 2 xy Similarly, the superscripts 1 and 2 represent the stress provided by a pair of 45-degree first rods and the stress provided by a pair of 90-degree second rods, respectively.

[0038] In the three-dimensional case, the stress vector is composed of the third cross rod with a 45-degree angle with the coordinate axis and three mutually perpendicular fourth cross rods on three sets of coordinate planes. The stress in each direction is aggregated by the rods that contribute to the component force in that direction. The stress state can be expressed as:

[0039]

[0040] Among them, σ x , σ y , σ z , τ xy , τ yz , τ xz They represent the normal stress along the x-axis, the normal stress along the y-axis, the normal stress along the z-axis, and the shear stress in the xy, xz, and xz directions. 1 x , σ 1 y , τ 1 xy represents the normal stress and shear stress calculated by the 45-degree crossbar on the xoy surface, σ 2 y , σ 2 z , τ 2 yz and σ 3 x , σ 3 z , τ 3 xz and σ 4 x , σ 4 y , σ 4 z represents the normal stress and shear stress calculated for three rods coincident with the coordinate axes.

[0041] See also Figure 3 As shown, the joint part is simulated by an axial shear coupling unit model, and the specific calculation process may include:

[0042] The mechanical relationship between shear stress and material tensile strength can be expressed by the Moore stress circle, from which we can get:

[0043] .

[0044] Furthermore, the shear stress of the joint part (grouting type) when it is sheared is:

[0045]

[0046] In addition, the shear resistance of the joint part also needs to consider the reduction factor of the joint shear failure , therefore, it can be obtained by using the formula:

[0047]

[0048] Obtain the shear resistance of the joint part in the shear direction ,in, is the reduction factor, which is determined by experiment; is the cross-sectional area of ​​the axial shear coupling unit (mm 2 ); is the axial tensile strength of the joint part (N / mm 2 ); is the shear stress of the axis-shear coupling unit in the tangential direction (N / mm 2 ); It is the angle of principal tensile stress when the material fails.

[0049] In step S2, the step of using the MYS model to simulate the nonlinear behavior of the soil unit under earthquake action and finally obtaining the seismic performance of the prefabricated subway station specifically includes:

[0050] Apply a horizontal seismic wave time history curve at the bottom of the soil unit and set shear boundary conditions on the left and right boundaries of the soil unit to simulate the simple shear motion caused by the earthquake;

[0051] The horizontal displacement, inter-layer displacement difference and stress distribution of the large-volume concrete parts and joint parts in the prefabricated subway station are recorded, and finally the seismic performance of the prefabricated subway station is verified through displacement difference, stiffness degradation and hysteresis behavior indicators.

[0052] The embodiment of the present invention further provides a seismic force response simulation system for a prefabricated subway station, comprising:

[0053] processor;

[0054] A memory for storing processor executable instructions; wherein the processor is configured to call the instructions stored in the memory to execute the above method.

[0055] The embodiment of the present invention further provides a computer-readable storage medium storing a computer program. The computer program can be executed by a processor of a device where the computer-readable storage medium is located to implement the above method.

[0056] It should be understood by those skilled in the art that embodiments of the present invention may be provided as methods, devices, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present invention 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, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0057] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program 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 produce 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.

[0058] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising 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.

[0059] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing 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.

[0060] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for simulating the anti-seismic force response of a prefabricated subway station, characterized in that: The following steps are included: S1, constructing a geometric model of the prefabricated subway station in a finite element software platform, wherein the geometric model includes a large-volume concrete part and a joint part, wherein the large-volume concrete part is simulated using a polymer rod model, and the joint part is simulated using an axial-shear coupling unit model; the simulation of the large-volume concrete part using the polymer rod model specifically includes: in a two-dimensional model, defining that there are four mutually intersecting but independent single-axis rods in the structure of the prefabricated subway station, and in a three-dimensional model, defining that there are nine single-axis rods in the structure of the prefabricated subway station, and the stress contribution of each single-axis rod is summarized into an overall stress state through coordinate transformation; In the two-dimensional model, it is defined that there are four mutually intersecting but independent single-axis rods in the structure of the prefabricated subway station, and in the three-dimensional model, it is defined that there are nine single-axis rods in the structure of the prefabricated subway station, and the stress contribution of each single-axis rod is summarized into the overall stress state through coordinate transformation, specifically including: By formula: Obtain the stress state generated by the four intersecting but independent uniaxial rods in the two-dimensional model, where σ x , σ y , τ xy Represents the normal stress along the x-axis (N / mm 2 ), normal stress along the y-axis and shear stress along the xy-axis (N / mm 2 );σ 1 x , σ 1 y , τ 2 xy and σ 2 x , σ 2 y , τ 2 xy , superscripts 1 and 2 represent the stress provided by a pair of 45-degree cross bars and a pair of 90-degree cross bars, respectively; By formula: Obtain the stress state generated by the nine uniaxial rods in the three-dimensional model; where σ x , σ y , σ z , τ xy , τ yz , τ xz Respectively represent the normal stress along the x-axis, the normal stress along the y-axis, the normal stress along the z-axis, and the shear stress in the xy, yz, and xz directions; σ 1 x , σ 1 y , τ 1 xy represents the normal stress and shear stress calculated by the 45-degree crossbar on the xoy surface, σ 2 y , σ 2 z , τ 2 yz ; σ 3 x , σ 3 z , τ 3 xz ; σ 4 x , σ 4 y , σ 4 z represents the normal stress and shear stress calculated for three bars coinciding with the coordinate axes; The joint part is simulated by the axial shear coupling unit model, which specifically includes: through the formula: Obtain the shear resistance of the joint part in the shear direction ,in, is the reduction factor, which is determined by experiment; is the cross-sectional area of ​​the axial shear coupling unit (mm 2 ); is the axial tensile strength of the joint part (N / mm 2 ); is the shear stress of the axis-shear coupling unit in the tangential direction (N / mm 2 ); is the angle of the principal tensile stress when the material fails; S2, establish a soil unit around the prefabricated subway station; and connect it with the structure of the prefabricated subway station through a common node method to simulate the rigid connection between the soil unit and the prefabricated subway station; then use the MYS model to simulate the nonlinear behavior of the soil unit under the action of an earthquake, and finally obtain the seismic performance of the prefabricated subway station.

2. The anti-seismic force response simulation method for a prefabricated subway station according to claim 1, characterized in that: In step S2, the step of using the MYS model to simulate the nonlinear behavior of the soil unit under earthquake action and finally obtaining the seismic performance of the prefabricated subway station specifically includes: Apply a horizontal seismic wave time history curve at the bottom of the soil unit and set shear boundary conditions on the left and right boundaries of the soil unit to simulate the simple shear motion caused by the earthquake; The horizontal displacement, inter-layer displacement difference and stress distribution of the large-volume concrete parts and joint parts in the prefabricated subway station are recorded, and finally the seismic performance of the prefabricated subway station is verified through displacement difference, stiffness degradation and hysteresis behavior indicators.

3. A seismic force response simulation system for a prefabricated subway station, characterized in that: include: processor; A memory for storing processor-executable instructions; wherein the processor is configured to call the instructions stored in the memory to execute the method according to any one of claims 1-2.

4. A computer-readable storage medium, characterized in that: A computer program is stored, and the computer program can be executed by a processor of the device where the computer-readable storage medium is located to implement the method according to any one of claims 1 to 2.

Images