CN111046502A - Method and device for calculating stiffness of soil spring of pipeline crossing fault - Google Patents
Method and device for calculating stiffness of soil spring of pipeline crossing fault Download PDFInfo
- Publication number
- CN111046502A CN111046502A CN201911109371.2A CN201911109371A CN111046502A CN 111046502 A CN111046502 A CN 111046502A CN 201911109371 A CN201911109371 A CN 201911109371A CN 111046502 A CN111046502 A CN 111046502A
- Authority
- CN
- China
- Prior art keywords
- soil
- finite element
- pipeline
- model
- dimensional finite
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
The invention relates to a method and a device for calculating the stiffness of a soil spring of a pipeline crossing fault, wherein the method comprises the steps of establishing a three-dimensional finite element model of the pipeline crossing fault; and acquiring a three-way soil spring stiffness model based on the three-dimensional finite element model and the pipe trench parameters. The method and the device for calculating the stiffness of the soil spring of the pipeline crossing the fault provided by the embodiment of the invention determine a method for calculating the three-way soil spring on the basis of considering the pipe ditch parameters. A three-dimensional finite element model of the pipeline crossing fault is established by using ABAQUS finite element software, and the three-dimensional finite element model is verified by an ASCE guideline method. On the basis of the three-dimensional finite element model, a three-way soil spring stiffness model of the pipeline crossing fault is constructed by using 1stOpt fitting software. And a reference function is provided for pipeline design and safety evaluation under the action of faults.
Description
Technical Field
The invention relates to the field of pipeline design, in particular to a method and a device for calculating the stiffness of a soil spring of a pipeline crossing fault.
Background
With the accelerated construction and development of natural gas pipe networks in China, high-grade steel gas pipelines inevitably pass through fault zones, and the safe operation of the pipelines can be effectively guaranteed by performing anti-seismic analysis on the high-grade steel gas pipelines. The establishment of the soil spring stiffness model is the basis for researching the pipeline reaction rule.
The earthquake is one of the potential damage causes of the long-distance pipeline, and the analysis of the reaction rule of the pipeline under the action of the fault can effectively prevent and reduce the damage in the future earthquake. In the currently common research methods, Liu Xiao Ben, Zhang hong and the like provide a calculation method of X80 pipeline design strain under the action of a three-dimensional oblique reverse fault based on finite element analysis results; the Cheng Xu Dong and Pangming have proposed the regression formula of the maximum tensile strain and the maximum compressive strain of the buried pipeline crossing the inclined slip fault.
However, the above method does not take into account the trench parameters. In fact, the deformation of the pipeline can be effectively reduced by changing the parameters of the pipe ditch so as to solve the problem of overhigh strain when the natural gas pipeline passes through the fault. Currently, the determination of the earth spring stiffness generally adopts a calculation method in an ASCE guide. However, the determination of the stiffness of the soil spring in the ASCE guide is established based on the same soil characteristics in an infinite large range outside the pipeline, and cannot reflect the situation that the characteristics of backfill soil in the pipe ditch are different from the original soil outside the pipe ditch, so that the influence of the size and the shape of the pipe ditch on the deformation of the pipeline cannot be reflected, and the pipe ditch parameter is one of main control parameters and cannot be ignored.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a method and a device for calculating the stiffness of a soil spring of a pipeline crossing fault.
The technical scheme for solving the technical problems is as follows:
in a first aspect, the invention provides a method for calculating the stiffness of a soil spring of a pipeline crossing fault, which comprises the following steps:
establishing a three-dimensional finite element model of the pipeline crossing fault;
and acquiring a three-way soil spring stiffness model based on the three-dimensional finite element model and the pipe trench parameters.
Further, the establishing of the three-dimensional finite element model of the pipeline under the fault action specifically includes:
defining a pipeline crossing fault as a dual nonlinear problem of material nonlinearity and geometric nonlinearity, establishing a three-dimensional finite element model of the pipeline crossing fault, and determining a load boundary condition and a pipe-soil characteristic.
Further, after building the three-dimensional finite element model of the conduit-traversing fault, the method further comprises:
and carrying out validity verification on the three-dimensional finite element model based on ASCE guidelines.
Further, the performing validity verification on the three-dimensional finite element model based on the ASCE guideline specifically includes: and obtaining a three-way soil spring stiffness calculation result of the three-dimensional finite element model, and comparing and verifying the three-way soil spring stiffness calculation result with a three-way soil spring stiffness calculation result obtained by adopting an ASCE (automatic sampling and accounting) guideline.
Further, the three-way soil spring stiffness model comprises a horizontal and transverse soil spring stiffness calculation model, a pipe shaft direction soil spring stiffness calculation model and a vertical direction soil spring stiffness calculation model.
Further, the calculation model of the horizontal transverse soil spring stiffness is as follows:
in the formula, ps: the pressure applied to the outer surface of the horizontal and transverse pipeline in unit length is kN/m;
c: cohesive force of clay, MPa;
d: pipe diameter, m;
h: depth of pipeline burial, m;
ρ0: sand density, kg/m 3;
b: widening margin m;
h: thickness of the cushion layer, m;
β slope, rad;
the calculation model of the stiffness of the soil spring in the pipe shaft direction is as follows:
fs＝(0.313+1.194H+0.052H2)(-0.616+6.969D-0.623D2)[2.165sin(β)+2.033cos(β)](2.852-0.087b-0.005b2)
in the formula (f)sThe friction force of unit length between the soil and the outer surface of the pipeline in the pipe axis direction is expressed, kN/m;
the vertical direction soil spring stiffness calculation model comprises a vertical upward soil spring stiffness calculation model and a vertical downward soil spring stiffness calculation model;
wherein, the calculation model of the stiffness of the vertical upward soil spring is as follows:
(94.243H-34.792H2+4.372H3-80.145)(1.654D-6.194D2+1.291D3-0.232)
in the formula, qs1Representing the pressure applied to the pipeline in unit length vertically upwards, kN/m;
the calculation model of the vertical downward soil spring stiffness is as follows:
in the formula, qs2Indicating the pressure experienced by a vertically downward unit length pipe; ω represents the clay friction angle.
In a second aspect, the present invention provides a device for calculating the soil spring rate of a pipeline crossing fault, comprising:
the finite element model establishing module is used for establishing a three-dimensional finite element model of the pipeline crossing fault;
and the soil spring stiffness calculation module is used for acquiring a three-way soil spring stiffness model based on the three-dimensional finite element model and the pipe trench parameters.
Further, the finite element model building module is specifically configured to:
defining a pipeline crossing fault as a dual nonlinear problem of material nonlinearity and geometric nonlinearity, establishing a three-dimensional finite element model of the pipeline crossing fault, and determining a load boundary condition and a pipe-soil characteristic.
In a third aspect, an embodiment of the present invention provides an electronic device, including a processor, a communication interface, a memory, and a bus, where the processor and the communication interface, the memory complete communication with each other through the bus, and the processor may call a logic instruction in the memory to perform the steps of the method provided in the first aspect.
In a fourth aspect, the invention provides a non-transitory computer readable storage medium having stored thereon a computer program which, when executed by a processor, performs the steps of the method as provided in the first aspect.
The method and the device for calculating the stiffness of the soil spring of the pipeline crossing the fault provided by the embodiment of the invention determine a method for calculating the three-way soil spring on the basis of considering the pipe ditch parameters. A three-dimensional finite element model of the pipeline crossing fault is established by using ABAQUS finite element software, and the three-dimensional finite element model is verified by an ASCE guideline method. On the basis of the three-dimensional finite element model, a three-way soil spring stiffness model of the pipeline crossing fault is constructed by using 1stOpt fitting software. And a reference function is provided for pipeline design and safety evaluation under the action of faults.
Drawings
FIG. 1 is a schematic flow chart of a method for calculating the stiffness of a soil spring of a pipeline crossing fault according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a soil-in-pipe contact analysis model provided by an embodiment of the present invention;
fig. 3 is a transverse cross-sectional view of a buried pipeline provided by an embodiment of the present invention;
FIG. 4(a) is a schematic diagram of a three-dimensional finite element model of a pipeline crossing fault according to an embodiment of the present invention;
FIG. 4(b) is a schematic diagram of meshing a three-dimensional finite element model according to an embodiment of the present invention;
FIG. 5 is a schematic diagram showing a relationship between a calculation result of a triaxial soil spring fitting formula and actual working condition data;
FIG. 6 is a schematic structural diagram of a device for calculating the earth spring stiffness of a pipeline crossing fault according to an embodiment of the present invention;
fig. 7 is a schematic physical structure diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Aiming at the problem that the reaction rule research of the pipeline under the action of the current fault does not fully consider the parameters of the pipe ditch, the embodiment of the invention provides a method for calculating the stiffness of a soil spring of the pipeline crossing the fault. Fig. 1 is a schematic flow chart of a method for calculating the stiffness of a soil spring of a pipeline crossing fault according to an embodiment of the present invention, and with reference to fig. 1, the method includes:
In the finite element model, the soil-in-pipe interaction was simulated using soil springs, including pipe-axis direction soil springs, horizontal lateral soil springs, and vertical direction soil springs. Currently, the determination of the earth spring stiffness generally adopts a calculation method in an ASCE guide. However, the determination of the spring stiffness of the soil in the ASCE guide is based on the same soil characteristics in an infinite range outside the pipeline, and cannot reflect the situation when the characteristics of the backfill soil in the pipe trench are different from the original soil outside the pipe trench, so that the influence of the size and shape of the pipe trench on the deformation of the pipeline cannot be reflected.
In order to study the law of action between pipe and soil, the present embodiment establishes a three-dimensional pipe and soil contact analysis model before step 101 is executed. Fig. 2 is a schematic diagram of a pipe-soil contact analysis model provided in an embodiment of the present invention, and fig. 3 is a transverse cross-sectional view of a buried pipeline provided in an embodiment of the present invention. As shown in fig. 2 and 3, the embodiment of the invention can obtain the buckling displacement of the field soil in all directions and the corresponding yield force by applying the displacement condition.
Furthermore, a three-dimensional finite element model of the pipeline crossing fault is established by using ABAQUS finite element software, and the buried pipeline reacts to the double nonlinear problem of material nonlinearity and geometric nonlinearity under the action of the fault. Therefore, the pipeline in the embodiment adopts the shell unit model to better simulate the actual situation. The pipeline model adopts a four-node curved surface thin shell reduction integral unit (S4R) and limited membrane strain. The soil model uses eight-node linear hexahedron cells (C3D 8). The three-dimensional finite element model and meshing are shown in fig. 4(a) and 4 (b). FIG. 4(a) is a schematic diagram of a three-dimensional finite element model of a pipeline crossing fault according to an embodiment of the present invention; fig. 4(b) is a schematic diagram of meshing a three-dimensional finite element model according to an embodiment of the present invention.
And 102, acquiring a three-way soil spring stiffness model based on the three-dimensional finite element model and the pipe trench parameters.
Specifically, on the basis of the three-dimensional finite element model, the pipe ditch parameters are considered, and multiple groups of buckling displacement-yield force data are obtained by changing single variables such as field conditions, backfill soil properties, pipe ditch sizes, pipeline sizes and pipeline burial depths. And then, fitting a three-way soil spring stiffness model considering the pipe ditch parameters by using a global optimization algorithm of 1stOpt software. The pipe trench parameters in this embodiment may include the soil covering margin, the soil thickness, the over-digging depth, the widening margin, the backfill fine soil, the side slope ratio, and the like of the pipe trench. By setting the pipe trench parameters, the reliability of the three-way soil spring stiffness model is improved.
The method for calculating the stiffness of the soil spring of the pipeline crossing the fault determines a method for calculating a three-way soil spring on the basis of considering pipe ditch parameters. A three-dimensional finite element model of the pipeline crossing fault is established by using ABAQUS finite element software, and the three-dimensional finite element model is verified by an ASCE guideline method. On the basis of the three-dimensional finite element model, a three-way soil spring stiffness model of the pipeline crossing fault is constructed by using 1stOpt fitting software. And a reference function is provided for pipeline design and safety evaluation under the action of faults.
Based on the content of the foregoing embodiment, in step 101, the establishing a three-dimensional finite element model of a pipeline under a fault specifically includes:
defining a pipeline crossing fault as a dual nonlinear problem of material nonlinearity and geometric nonlinearity, establishing a three-dimensional finite element model of the pipeline crossing fault, and determining a load boundary condition and a pipe-soil characteristic.
Specifically, in the present embodiment, a west-gas-east pipeline is taken as a research object, and a process of establishing a three-dimensional finite element model of a pipeline under a fault effect is described. The pipeline crossing fault is defined as a dual nonlinear problem of material nonlinearity and geometric nonlinearity, and the pipeline adopts a shell unit model in the embodiment. Wherein, the shell element model pipeline adopts four nodes, a curved surface thin shell reduction integral unit (S4R) and limited membrane strain; an eight-node linear hexahedron cell (C3D8) of the soil model. A three-dimensional finite element model is established, and the model and meshing are shown in fig. 4(a) and 4 (b). FIG. 4(a) is a schematic diagram of a three-dimensional finite element model of a pipeline crossing fault according to an embodiment of the present invention; fig. 4(b) is a schematic diagram of meshing a three-dimensional finite element model according to an embodiment of the present invention. A master-slave contact algorithm is employed in the pipe-soil contact setup: the main surface is the outer surface of the pipeline, and the secondary surface is the surface of the soil body at the pipe-soil contact position. The interaction between the pipe and the soil interface is described in terms of both normal and tangential action. The tangential action defines the friction coefficient between the soil and the pipe through a penalty function, and the normal action is set as hard contact, namely when the contact pressure between the contact surfaces is considered to be zero or negative, the two contact surfaces are separated, and the mutual constraint relationship is removed.
Further, load boundary conditions are determined. The load boundary condition includes a load condition and a boundary condition, wherein:
(1) load condition
Analyzing displacement components according to different fault types, and applying the displacement components to the pipe-soil unit as displacement load conditions, wherein for the walk-slip fault, the motion occurs in a horizontal plane, the intersection angle of a pipeline crossing fault is β, the fault inclination angle is psi, and the pipeline crossing horizontal displacement is deltasWhen the fault is slipped, the axial displacement component delta x and the lateral displacement component delta y along the pipeline are respectively as follows:
pipeline crossing vertical displacement is deltapThe horizontal component Δ x, the horizontal lateral component Δ y, and the vertical component Δ z along the pipeline are:
for having horizontal displacement deltasAnd a vertical displacement deltapThe displacement component of the fault activity is:
δptaking a positive value for the positive fault and a negative value for the negative fault. DeltasThe right slip fault is positive and the left slip fault is negative the intersection angle β between the pipe and fault is defined as the angle between the direction of motion of the right ground of the fault and the axis of the pipe to the right.
(2) Boundary condition
The boundary conditions are set such that the upper surface is a free surface, the lower surface constrains all degrees of freedom of displacement, and the other surfaces constrain their normal degrees of freedom. The load-displacement curve can be obtained in two ways, namely, the displacement is solved by the given load, or the corresponding load is solved by the given displacement. In order to avoid the convergence problem caused by the reduction of rigidity (even approaching 0) when the limit condition is reached, the model adopts the method of solving load by given displacement, and the resistance of the soil body in each direction is obtained by applying the displacement.
Further, determining pipe-soil characteristics, including pipeline material characteristics and soil parameter characteristics, wherein:
(1) pipe material properties
The pipe material adopts a Ramberg-Osgood constitutive model, and the parameters of the pipe material are determined by fitting a real stress-strain curve of the material measured by a tensile experiment.
In the formula:
epsilon: true strain;
σ: axial tensile stress, MPa;
e: modulus of elasticity, N/mm;
σ0: yield stress, MPa;
α, r Ramberg-Osgood parameters.
Table 1 shows the sigma values of the steel pipes of the usual grades0α and r values.
TABLE 1 Ramberg-Osgood parameters for conventional pipe steels
Note: longitudinal submerged arc welded pipe technology based on strain design according to western gas and east gas transmission line natural gas pipeline engineering
The steel pipes produced under the operation condition are divided into two grades of X80HD1 and X80HD2, which are respectively suitable for bending the steel pipes
Areas with a bending strain under load of 1.0% and less than 1.5%.
(2) Characteristic of soil parameters
Soil parameters were obtained based on the existing literature and are shown in table 2 below:
TABLE 2 soil Material basis parameters
Based on the content of the foregoing embodiment, as an alternative embodiment, after the step 101 of establishing a three-dimensional finite element model of a pipeline crossing fault, the method further includes:
and carrying out validity verification on the three-dimensional finite element model based on ASCE guidelines.
Specifically, the calculation result of the stiffness of the three-way soil spring of the three-dimensional finite element model is obtained, and is compared with the calculation result of the stiffness of the three-way soil spring by adopting an ASCE guideline for verification, and the verification result is shown in table 3:
TABLE 3 model verification results
As can be seen from Table 3, the error can be controlled within 10%, and the three-dimensional finite element model is feasible.
Based on the content of the foregoing embodiment, as an optional embodiment, in step 102, the three-way soil spring stiffness model includes a horizontal transverse soil spring stiffness calculation model, a pipe axis direction soil spring stiffness calculation model, and a vertical direction soil spring stiffness calculation model.
Wherein, the horizontal soil spring stiffness calculation model is as follows:
in the formula, ps: the pressure applied to the outer surface of the horizontal and transverse pipeline in unit length is kN/m;
c: cohesive force of clay, MPa;
d: pipe diameter, m;
h: depth of pipeline burial, m;
ρ0: sand density, kg/m 3;
b: widening margin m;
h: thickness of the cushion layer, m;
β slope, rad;
the calculation model of the stiffness of the soil spring in the pipe shaft direction is as follows:
fs＝(0.313+1.194H+0.052H2)(-0.616+6.969D-0.623D2)[2.165sin(β)+2.033cos(β)](2.852-0.087b-0.005b2)
in the formula (f)sThe friction force of unit length between the soil and the outer surface of the pipeline in the pipe axis direction is expressed, kN/m;
the vertical direction soil spring stiffness calculation model comprises a vertical upward soil spring stiffness calculation model and a vertical downward soil spring stiffness calculation model;
the calculation model of the stiffness of the vertical upward soil spring is as follows:
(94.243H-34.792H2+4.372H3-80.145)(1.654D-6.194D2+1.291D3-0.232)
in the formula, qs1Representing the pressure applied to the pipeline in unit length vertically upwards, kN/m;
the calculation model of the vertical downward soil spring stiffness is as follows:
in the formula, qs2Indicating the pressure experienced by a vertically downward unit length pipe; ω represents the clay friction angle.
Based on the content of the above embodiment, in order to verify the accuracy of the three-way soil spring stiffness model, in this embodiment, according to geological survey data of the second-line engineering of the west gas and east gas transmission, 40 groups of actual strain data of the second-line pipelines of the west gas and east gas transmission are extracted and compared with the calculation result of the three-way soil spring stiffness fitting formula for analysis. The fitting formula of the three-way soil spring stiffness refers to a three-way soil spring stiffness model. Fig. 5 is a schematic diagram showing a relationship between a calculation result of a triaxial soil spring fitting formula and actual condition data, and an abscissa and an ordinate in fig. 5 respectively indicate a strain value of a pipeline calculated by using the triaxial soil spring fitting formula provided by the present patent and a strain value of the pipeline under the actual condition. As shown in fig. 5, the maximum strain error is 3.986%, the minimum strain error is 0.011915%, and comparative analysis shows that the fitting calculation result and the actual working condition data have the same variation trend, and the fitting formula has high fitting degree and certain accuracy.
Fig. 6 is a schematic structural diagram of an apparatus for calculating a stiffness of a soil spring of a pipeline crossing fault according to an embodiment of the present invention, and referring to fig. 6, the present invention provides an apparatus for calculating a stiffness of a soil spring of a pipeline crossing fault, including:
a finite element model establishing module 601, configured to establish a three-dimensional finite element model of a pipeline crossing fault;
and the soil spring stiffness calculation module 602 is configured to obtain a three-way soil spring stiffness model based on the three-dimensional finite element model and the pipe trench parameters.
Specifically, the device for calculating the stiffness of the earth spring of the pipeline crossing the fault according to the embodiment of the present invention is specifically configured to execute the steps of the method for calculating the stiffness of the earth spring of the pipeline crossing the fault according to the embodiment of the present invention, and since the method for calculating the stiffness of the earth spring of the pipeline crossing the fault has been described in detail in the embodiment, no further description is given to functional modules of the device for calculating the stiffness of the earth spring of the pipeline crossing the fault.
The method for calculating the stiffness of the soil spring of the pipeline crossing the fault determines a method for calculating a three-way soil spring on the basis of considering pipe ditch parameters. A three-dimensional finite element model of the pipeline crossing fault is established by using ABAQUS finite element software, and the three-dimensional finite element model is verified by an ASCE guideline method. On the basis of the three-dimensional finite element model, a three-way soil spring stiffness model of the pipeline crossing fault is constructed by using 1stOpt fitting software. And a reference function is provided for pipeline design and safety evaluation under the action of faults.
Based on the content of the foregoing embodiment, the finite element model building module is specifically configured to:
defining a pipeline crossing fault as a dual nonlinear problem of material nonlinearity and geometric nonlinearity, establishing a three-dimensional finite element model of the pipeline crossing fault, and determining a load boundary condition and a pipe-soil characteristic.
Fig. 7 is a schematic entity structure diagram of an electronic device according to an embodiment of the present invention, and as shown in fig. 7, the electronic device may include: a processor (processor)701, a communication Interface (Communications Interface)702, a memory (memory)703 and a communication bus 704, wherein the processor 701, the communication Interface 702 and the memory 703 complete communication with each other through the communication bus 704. The processor 701 may call a computer program stored on the memory 703 and executable on the processor 701 to perform the method for calculating the earth spring stiffness of the pipe crossing fault provided by the above embodiments, for example, the method includes: establishing a three-dimensional finite element model of the pipeline crossing fault; and acquiring a three-way soil spring stiffness model based on the three-dimensional finite element model and the pipe trench parameters.
Embodiments of the present invention further provide a non-transitory computer-readable storage medium, on which a computer program is stored, where the computer program is implemented to perform the method for calculating the earth spring stiffness of a pipeline crossing a fault provided in the foregoing embodiments, for example, the method includes: establishing a three-dimensional finite element model of the pipeline crossing fault; and acquiring a three-way soil spring stiffness model based on the three-dimensional finite element model and the pipe trench parameters.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. A method for calculating the stiffness of a soil spring of a pipeline crossing fault is characterized by comprising the following steps of;
establishing a three-dimensional finite element model of the pipeline crossing fault;
and acquiring a three-way soil spring stiffness model based on the three-dimensional finite element model and the pipe trench parameters.
2. The method according to claim 1, wherein the establishing of the three-dimensional finite element model of the faulted pipe comprises in particular:
defining a pipeline crossing fault as a dual nonlinear problem of material nonlinearity and geometric nonlinearity, establishing a three-dimensional finite element model of the pipeline crossing fault, and determining a load boundary condition and a pipe-soil characteristic.
3. The method of claim 1, wherein after building the three-dimensional finite element model of the cross-fault pipeline, the method further comprises:
and carrying out validity verification on the three-dimensional finite element model based on ASCE guidelines.
4. The method according to claim 3, wherein the validation of the three-dimensional finite element model based on ASCE guidelines comprises:
and obtaining a three-way soil spring stiffness calculation result of the three-dimensional finite element model, and comparing and verifying the three-way soil spring stiffness calculation result with a three-way soil spring stiffness calculation result obtained by adopting an ASCE (automatic sampling and accounting) guideline.
5. The method of claim 1, wherein the three-way soil spring rate model comprises a horizontal transverse soil spring rate calculation model, a pipe axis direction soil spring rate calculation model, and a vertical direction soil spring rate calculation model.
6. The method of claim 5, wherein the horizontal lateral earth spring rate calculation model is:
in the formula, ps: the pressure applied to the outer surface of the horizontal and transverse pipeline in unit length is kN/m;
c: cohesive force of clay, MPa;
d: pipe diameter, m;
h: depth of pipeline burial, m;
ρ0: sand density, kg/m 3;
b: widening margin m;
h: thickness of the cushion layer, m;
β slope, rad;
the calculation model of the stiffness of the soil spring in the pipe shaft direction is as follows:
fs＝(0.313+1.194H+0.052H2)(-0.616+6.969D-0.623D2)[2.165sin(β)+2.033cos(β)](2.852-0.087b-0.005b2)
in the formula (f)sThe friction force of unit length between the soil and the outer surface of the pipeline in the pipe axis direction is expressed, kN/m;
the vertical direction soil spring stiffness calculation model comprises a vertical upward soil spring stiffness calculation model and a vertical downward soil spring stiffness calculation model;
wherein, the calculation model of the stiffness of the vertical upward soil spring is as follows:
in the formula, qs1Representing the pressure applied to the pipeline in unit length vertically upwards, kN/m;
the calculation model of the vertical downward soil spring stiffness is as follows:
in the formula, qs2Indicating the pressure experienced by a vertically downward unit length pipe; ω represents the clay friction angle.
7. An earth spring rate calculation device for a pipeline crossing fault, comprising:
the finite element model establishing module is used for establishing a three-dimensional finite element model of the pipeline crossing fault;
and the soil spring stiffness calculation module is used for acquiring a three-way soil spring stiffness model based on the three-dimensional finite element model and the pipe trench parameters.
8. The apparatus of claim 7, wherein the finite element modeling module is specifically configured to:
defining a pipeline crossing fault as a dual nonlinear problem of material nonlinearity and geometric nonlinearity, establishing a three-dimensional finite element model of the pipeline crossing fault, and determining a load boundary condition and a pipe-soil characteristic.
9. An electronic device, comprising a processor, a communication interface, a memory and a bus, wherein the processor, the communication interface and the memory communicate with each other via the bus, and the processor can call logic instructions in the memory to execute the method according to claim 7 or 8.
10. A non-transitory computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the method according to claim 7 or 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911109371.2A CN111046502B (en) | 2019-11-13 | 2019-11-13 | Soil spring stiffness calculation method and device for pipeline crossing fault |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911109371.2A CN111046502B (en) | 2019-11-13 | 2019-11-13 | Soil spring stiffness calculation method and device for pipeline crossing fault |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111046502A true CN111046502A (en) | 2020-04-21 |
| CN111046502B CN111046502B (en) | 2023-08-25 |
Family
ID=70231885
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911109371.2A Active CN111046502B (en) | 2019-11-13 | 2019-11-13 | Soil spring stiffness calculation method and device for pipeline crossing fault |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111046502B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112651151A (en) * | 2020-12-01 | 2021-04-13 | 中国石油大学(北京) | Test device for simulating pipeline to pass through slip fault and determination method |
| CN112949117A (en) * | 2021-02-02 | 2021-06-11 | 天津城建大学 | Three-dimensional strain analysis method for asphalt pavement based on multi-dimensional parameters |
| CN113722841A (en) * | 2020-05-25 | 2021-11-30 | 中国石油天然气股份有限公司 | Method and device for establishing three-dimensional model of long-distance pipeline |
| CN114440002A (en) * | 2020-11-02 | 2022-05-06 | 中国石油天然气集团有限公司 | Laying structure and method for relieving strain of oil-gas pipeline under action of slip fault |
| CN115130259A (en) * | 2022-06-24 | 2022-09-30 | 武汉大学 | Method and device for predicting axial displacement of cross-fault pipeline interface and building model |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102943945A (en) * | 2012-12-07 | 2013-02-27 | 中国石油集团工程设计有限责任公司 | Earthquake-proof protection method for pipeline projects passing through fault zones |
| CN104933269A (en) * | 2015-07-13 | 2015-09-23 | 中国石油天然气集团公司 | Design method for oil and gas pipeline crossing earthquake fault |
| CN105155376A (en) * | 2015-08-17 | 2015-12-16 | 苏州同尚工程设计咨询有限公司 | Non-pressed pipeline reinforcing structure for road surface and construction method for non-pressed pipeline reinforcing structure |
| CN105760622A (en) * | 2016-03-16 | 2016-07-13 | 上海市城市建设设计研究总院 | Soil spring design method of shield tunnel beam-spring model and calculation model |
| CN106777454A (en) * | 2016-11-09 | 2017-05-31 | 中国石油天然气集团公司 | A kind of method for designing of Pipeline Crossing Program strike-slip fault |
| CN107563014A (en) * | 2017-08-11 | 2018-01-09 | 西南石油大学 | The computational methods of the strain of pipeline flexing and critical fault displcement under a kind of faulting |
| CN107741324A (en) * | 2017-10-13 | 2018-02-27 | 北京工业大学 | A kind of housing washer fault section diagnosis method |
| CN108457311A (en) * | 2017-12-13 | 2018-08-28 | 上海交通大学 | A kind of deep base pit enclosure wall stress deformation quick calculation method considering the coupling of wall soil |
| CN207961907U (en) * | 2018-01-31 | 2018-10-12 | 中国石油大学(北京) | Buried oil-gas pipeline antishock device and buried oil-gas pipeline system |
-
2019
- 2019-11-13 CN CN201911109371.2A patent/CN111046502B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102943945A (en) * | 2012-12-07 | 2013-02-27 | 中国石油集团工程设计有限责任公司 | Earthquake-proof protection method for pipeline projects passing through fault zones |
| CN104933269A (en) * | 2015-07-13 | 2015-09-23 | 中国石油天然气集团公司 | Design method for oil and gas pipeline crossing earthquake fault |
| CN105155376A (en) * | 2015-08-17 | 2015-12-16 | 苏州同尚工程设计咨询有限公司 | Non-pressed pipeline reinforcing structure for road surface and construction method for non-pressed pipeline reinforcing structure |
| CN105760622A (en) * | 2016-03-16 | 2016-07-13 | 上海市城市建设设计研究总院 | Soil spring design method of shield tunnel beam-spring model and calculation model |
| CN106777454A (en) * | 2016-11-09 | 2017-05-31 | 中国石油天然气集团公司 | A kind of method for designing of Pipeline Crossing Program strike-slip fault |
| CN107563014A (en) * | 2017-08-11 | 2018-01-09 | 西南石油大学 | The computational methods of the strain of pipeline flexing and critical fault displcement under a kind of faulting |
| CN107741324A (en) * | 2017-10-13 | 2018-02-27 | 北京工业大学 | A kind of housing washer fault section diagnosis method |
| CN108457311A (en) * | 2017-12-13 | 2018-08-28 | 上海交通大学 | A kind of deep base pit enclosure wall stress deformation quick calculation method considering the coupling of wall soil |
| CN207961907U (en) * | 2018-01-31 | 2018-10-12 | 中国石油大学(北京) | Buried oil-gas pipeline antishock device and buried oil-gas pipeline system |
Non-Patent Citations (9)
| Title |
|---|
| DA HA: "Centrifuge Modeling of Earthquake Effects on Buried High-Density Polyethylene(HDPE) Pipelines Crossing Fault Zones", JOURNAL OF GEOTECHNICAL AND GEOENVIRONMENTAL ENGINEERING, pages 1501 - 1515 * |
| XIAOTING GU: "RESEARCH ON ASEISMATIC MEASURES OF GAS PIPELINE CROSSING A FAULT FOR STRAIN-BASED DESIGN", ASME, pages 1 - 10 * |
| 刘学杰: "地下管道穿越断层的应变设计方法", pages 81 - 85 * |
| 刘爱文: "地震断层作用下埋地管线的反应分析", 地震工程与工程振动, pages 22 - 27 * |
| 夏梦莹: "基于壳单元的连续型采空区埋地管道应变分析", 完整性与可靠性, pages 256 - 262 * |
| 李渊博: "岩土崩塌冲击作用下埋地管道应力与变形分析", 后勤工程学院学报, pages 31 - 35 * |
| 杨娜: "场地土对跨断层压力管道变形反应的影响研究", pages 029 - 673 * |
| 缪文振: "含腐蚀缺陷高强度管线钢在活动断层 中的反应分析", 中国优秀硕士论文电子期刊网, pages 019 - 72 * |
| 顾晓婷;张宏;王国丽;赵乐晋;: "大口径天然气管线穿越断层的管沟设计研究", no. 08, pages 113 - 115 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113722841A (en) * | 2020-05-25 | 2021-11-30 | 中国石油天然气股份有限公司 | Method and device for establishing three-dimensional model of long-distance pipeline |
| CN113722841B (en) * | 2020-05-25 | 2022-11-04 | 中国石油天然气股份有限公司 | Method and device for establishing three-dimensional model of long-distance pipeline |
| CN114440002A (en) * | 2020-11-02 | 2022-05-06 | 中国石油天然气集团有限公司 | Laying structure and method for relieving strain of oil-gas pipeline under action of slip fault |
| CN112651151A (en) * | 2020-12-01 | 2021-04-13 | 中国石油大学(北京) | Test device for simulating pipeline to pass through slip fault and determination method |
| CN112651151B (en) * | 2020-12-01 | 2023-06-20 | 中国石油大学(北京) | Test device and determination method for simulating pipeline to pass through walk-slip fault |
| CN112949117A (en) * | 2021-02-02 | 2021-06-11 | 天津城建大学 | Three-dimensional strain analysis method for asphalt pavement based on multi-dimensional parameters |
| CN115130259A (en) * | 2022-06-24 | 2022-09-30 | 武汉大学 | Method and device for predicting axial displacement of cross-fault pipeline interface and building model |
| CN115130259B (en) * | 2022-06-24 | 2023-10-03 | 武汉大学 | Cross-fault pipeline interface axial displacement prediction and model construction method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111046502B (en) | 2023-08-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111046502B (en) | Soil spring stiffness calculation method and device for pipeline crossing fault | |
| Ng et al. | Three-dimensional centrifuge and numerical modeling of the interaction between perpendicularly crossing tunnels | |
| CN106777454B (en) | Design method for pipeline crossing slip fault | |
| Zhang et al. | A new finite element model of buried steel pipelines crossing strike-slip faults considering equivalent boundary springs | |
| CN108385727B (en) | Method for calculating reasonable buried depth of submarine tunnel constructed by immersed tube method | |
| Kunert et al. | Nonlinear FEM strategies for modeling pipe–soil interaction | |
| Gantes et al. | Evaluation of seismic protection methods for buried fuel pipelines subjected to fault rupture | |
| Cheong et al. | 3D FE analyses of buried pipeline with elbows subjected to lateral loading | |
| Shi et al. | Buried pipeline responses to ground displacements induced by adjacent static pipe bursting | |
| Talebi et al. | A refined nonlinear analytical method for buried pipelines crossing strike‐slip faults | |
| Xu et al. | Review of dynamic response of buried pipelines | |
| He et al. | Numerical modelling of transverse cracking in embankment dams | |
| Zhang et al. | Effects of karst cave shape on the stability and minimum safety thickness of tunnel surrounding rock | |
| Triantafyllaki et al. | Offshore pipeline performance under strike-slip fault movements | |
| Qian et al. | Comparative study on interface elements, thin-layer elements, and contact analysis methods in the analysis of high concrete-faced rockfill dams | |
| Yang et al. | Theoretical analysis and finite element simulation of pipeline structure in liquefied soil | |
| Hassani et al. | A 3D numerical modeling of polyethylene buried pipes affected by fault movement | |
| CN116796396A (en) | Analysis method for deformation of underlying tunnel caused by foundation pit excavation and precipitation | |
| CN103835278B (en) | Buried pipe is subject to the safety monitoring appraisal procedure that contiguous Construction Vibration affects under water | |
| Xu et al. | Integrate pipe-soil interaction model with the vector form intrinsic finite element method-nonlinear analysis of free-span | |
| Pan et al. | Probabilistic seismic responses and failure analyses of free-spanning subsea pipelines under offshore spatial earthquake motions | |
| Triantafyllaki et al. | Design of offshore gas pipelines against active tectonic fault movement | |
| Zhang et al. | Dynamic back analysis of soil deformation during the construction of deep cantilever foundation pits | |
| Melissianos et al. | Protection measures for buried steel pipelines subjected to fault rupture | |
| CN115544702B (en) | Oil-gas pipeline immersed tube construction stress checking method, device, equipment and storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |