CN112668229A - Data simulation method for crack generation and expansion of initial dam of tailing dam - Google Patents
Data simulation method for crack generation and expansion of initial dam of tailing dam Download PDFInfo
- Publication number
- CN112668229A CN112668229A CN202110068826.1A CN202110068826A CN112668229A CN 112668229 A CN112668229 A CN 112668229A CN 202110068826 A CN202110068826 A CN 202110068826A CN 112668229 A CN112668229 A CN 112668229A
- Authority
- CN
- China
- Prior art keywords
- stress
- dam
- tailing dam
- hydraulic fracture
- state
- 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
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention provides a data simulation method for the crack generation and expansion of a tailing dam initial dam, which is used for pushing and establishing a finite element numerical simulation model for describing the hydraulic fracture generation and expansion process by combining the theory of the Biao consolidation so as to comprehensively and efficiently predict and judge the hydraulic fracture generation and expansion process, thereby ensuring the operation safety of the tailing dam.
Description
Technical Field
The invention relates to the technical field of hydraulic engineering of a tailing dam, in particular to a data simulation method for crack generation and expansion of an initial dam of the tailing dam.
Background
Hydraulic fracturing is the most interesting and controversial problem in tailing dam engineering. Hydraulic fracturing can lead to the destruction of the dam impervious body with catastrophic consequences. The reservoir water storage pressure can cause cracks on the upstream surface of the core wall, so that a channel for centralized water seepage is formed and the dam is damaged, and the problem and the debate are still questioned and debated in engineering and academia until now. In order to accurately determine the occurrence state of hydraulic fracture in the hydraulic tailing dam, effective and accurate prediction and judgment on the occurrence state and the evolution state of the hydraulic fracture are needed, but the effective and accurate prediction and judgment on the occurrence state and the evolution state of the hydraulic fracture are usually established based on a dispersion fracture theory at present, and the method cannot comprehensively and efficiently perform corresponding prediction and judgment.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a data simulation method for the occurrence and expansion of the initial dam crack of a tailing dam, and the numerical simulation method for the occurrence and expansion of hydraulic fracture comprises the following steps: step S1, constructing a finite element equation about the tailing dam, and calculating to obtain displacement, stress and pore pressure field corresponding to hydraulic fracture about the tailing dam; step S2, calculating the main stress parameter related to the hydraulic fracture according to the displacement, stress and pore pressure field of the hydraulic fracture; step S3, according to the main stress parameter of the hydraulic fracture, the tensile strength of the tailing dam is judged; step S4, according to the pulling rope strength judgment processing result, processing the tailing dam in different modes of cracking state or non-cracking state; therefore, the numerical simulation method for hydraulic fracture occurrence and expansion is combined with the Biao consolidation theory to deduce and establish a finite element numerical simulation model for describing the hydraulic fracture occurrence and expansion process, so that the hydraulic fracture occurrence and expansion process is comprehensively and efficiently predicted and judged, and the operation safety of the tailing dam is ensured.
The invention provides a data simulation method for crack generation and expansion of a tailing dam initial dam, which is characterized by comprising the following steps of:
step S1, constructing a finite element equation about the tailing dam, and calculating to obtain displacement, stress and pore pressure field about the corresponding hydraulic fracture of the tailing dam;
step S2, calculating to obtain a main stress parameter related to the hydraulic fracture according to the displacement, the stress and the pore pressure field of the hydraulic fracture;
step S3, according to the main stress parameter of the hydraulic fracture, the tensile strength of the tailing dam is judged;
step S4, according to the pulling rope strength judgment processing result, carrying out different mode processing aiming at the condition that the tailing dam is in a cracking state or in an uncracked state;
further, in the step S1, constructing finite element equations about the tailings dam, and calculating the displacement, stress and pore pressure field corresponding to the hydraulic fracture of the tailings dam specifically includes,
step S101, determining a control equation set of a proportional-integral theory about the tailing dam according to a calculation model about fluid-solid coupling effective stress;
step S102, carrying out transformation processing on the control equation set of the Biaoh consolidation theory to obtain a finite element format equation set of a two-dimensional Biaoh consolidation theory;
step S103, solving a finite element format equation set of the two-dimensional Biao consolidation theory to calculate and obtain the displacement, the stress and the pore pressure field of the hydraulic fracture;
further, in the step S101, according to the calculation model of the fluid-solid coupling effective stress, the specific control equation set of the proportion consolidation theory about the tailing dam is determined to comprise,
determining the system of the control equation of the Biao consolidation theory, which is composed of the following formulas (1) to (5) together, according to a calculation model about the effective stress of fluid-solid coupling
σij'＝σij-qδij＝Hijklεkl (2)
εkl＝1/2(uk,l+ul,k) (3)
In the above formulas (1) to (5), i, j are directions of the cauchy coordinate system; sigmaij,σij' Total stress and effective stress, respectively; fiA volumetric force component in the i direction; q is pore water pressure; deltaijKnown as the Kronecker (Kronecker) notation; hijklIs the constitutive relation tensor of the soil; epsilonijIs the strain tensor; u. ofijIs the displacement component; v. ofiIs the pore water flow velocity component; theta is the volume strain of soil; t is time; k is a radical ofijIs the permeability coefficient; gamma raywIs the volume weight of water;
alternatively, the first and second electrodes may be,
in step S102, transforming the proportional control equation set to obtain a finite element format equation of the two-dimensional proportional control theory specifically includes,
transforming the control equation set of the Biaoo consolidation theory to obtain a finite element format equation set of the two-dimensional Biaoo consolidation theory, which is formed by the following formulas (6) to (10)
In the above formulas (6) to (10), { δ }eRepresenting unit node displacement; { p }eRepresenting the pore water pressure of the unit node;
representing cell edge interface forces; { F } represents the unit volume force; v isnRepresenting cell boundary flow rate; [ K ]]Representing a permeability coefficient matrix; { K }eIs a cell stiffness matrix; [ K ]c]eCoupling a matrix for the cells; [ K ]s]eIs a unit seepage matrix; n is a radical of1...N4Is a unit shape function;
further, in the step S2, the calculating of the main stress parameter related to the hydraulic fracture according to the displacement, the stress and the pore pressure field of the hydraulic fracture specifically includes,
step S201, constructing a stress-strain relation of the compacted clay in the tailing dam in a stretching state according to the displacement, the stress and the pore pressure field of the hydraulic fracture;
step S202, determining a relation equation between stress and strain of the compacted clay in a stretching state according to the stress-strain relation of the compacted clay in the stretching state;
step S203, solving a relation equation between stress and strain of the compacted clay in a tensile state to calculate the minimum effective principal stress of the compacted clay as the principal stress parameter;
further, in the step S202, determining the relation equation between the stress and the strain of the compacted clay under the tensile state according to the relation between the stress and the strain of the compacted clay under the tensile state specifically comprises,
determining the relation equation between the stress and the strain of the compacted clay in the tensile state according to the stress-strain relation of the compacted clay in the tensile state in the following formula (11)
In the above formula (11), ε represents a stress; σ represents strain; e represents the initial elastic modulus of the stretching section; mu.s0The initial tangent elastic modulus reduction coefficient is used for describing the initial damage degree of the soil body before the soil body bears the tensile stress; epsilonfIs the peak strain; a. the1And B1Is a material constant, and
and
further, in the step S3, the determining the tensile strength of the tailing dam according to the main stress parameter of the hydraulic fracture specifically includes,
step S301A, acquiring the minimum effective main stress of the tailing dam about the compacted clay from the main stress parameters of the hydraulic fracture;
step S302A, comparing the minimum effective principal stress with the tensile strength of the compacted clay;
alternatively, the first and second electrodes may be,
in the step S3, the determining the tensile strength of the tailing dam according to the main stress parameter of the hydraulic fracture specifically includes,
step S301B, obtaining the maximum main stress F of the hydraulic fracture according to the main stress parameter of the hydraulic fractureMAnd minimum principal stress Fm；
Step S302B, calculating the hydraulic fracture pressure u according to the following formula (12)
In the formula (12), u is hydraulic fracture pressure, rho is compacted dry density of the tailing dam, delta is tensile strength of soil, and F isMIs the maximum principal stress of the hydraulic fracture, FmA minimum principal stress for the hydraulic fracture;
step S303B, calculating the tensile strength k of the tailing dam according to the following formula (13)
In the formula (13), k is tensile strength of the tailing dam, u is hydraulic fracture pressure, c is cohesion,
the internal friction angle is set, S is the contact area of the tailing dam and water, and S is the cross-sectional area of the tailing dam;
further, in the step S4, the performing different modes of processing on the tailings dam in a cracked state or an uncracked state according to the result of the pulling rope strength determination processing specifically includes,
step S401, if the result of the stay cord strength judgment processing indicates that the minimum effective principal stress is smaller than the tensile strength of the compacted clay, performing corresponding mode processing on the tailings dam in an uncracked state;
step S402, if the result of the stay cord strength judgment processing indicates that the minimum effective principal stress is greater than or equal to the tensile strength of the compacted clay, performing corresponding mode processing on the tailing dam in a cracking state;
further, in the step S401, if the result of the judgment processing of the rope pulling strength indicates that the minimum effective principal stress is smaller than the tensile strength of the compacted clay, the corresponding mode processing for the tailings dam in the non-cracked state specifically includes,
if the result of the stay cord strength judgment processing indicates that the minimum effective principal stress is less than the tensile strength of the compacted clay, repeating the steps S1-S4 aiming at the current tailing dam to realize the next load step numerical simulation;
further, in the step S402, if the result of the judgment processing of the rope pulling strength indicates that the minimum effective principal stress is greater than or equal to the tensile strength of the compacted clay, the corresponding mode processing for the tailing dam in the cracking state specifically includes,
s4021, calculating and obtaining a crack direction corresponding to a crack state of the tailing dam aiming at the current crack state;
step S4022, according to the crack direction, unit stiffness correction processing is carried out on the current crack state of the tailing dam so as to determine a corresponding step length correction mode;
s4023, setting an initial pore pressure field and an external load corresponding to the current cracking state of the tailing dam according to the step length correction mode;
further, in the step S4021, specifically calculating and obtaining a crack direction corresponding to the crack state of the tailing dam according to the current crack state includes,
s40211, determining the anisotropic characteristics of a cracked soil body in terms of stress-strain relationship and permeability characteristics according to the cracking state of the tailing dam at present;
step S40212, determining a relation curve between a permeability coefficient and normal effective stress according to the anisotropic characteristics of the cracked soil body in terms of stress-strain relation and permeability characteristic, and calculating to obtain a crack direction corresponding to the cracked state.
Compared with the prior art, the numerical simulation method for hydraulic fracture generation and expansion comprises the following steps: step S1, constructing a finite element equation about the tailing dam, and calculating to obtain displacement, stress and pore pressure field corresponding to hydraulic fracture about the tailing dam; step S2, calculating the main stress parameter related to the hydraulic fracture according to the displacement, stress and pore pressure field of the hydraulic fracture; step S3, according to the main stress parameter of the hydraulic fracture, the tensile strength of the tailing dam is judged; step S4, according to the pulling rope strength judgment processing result, processing the tailing dam in different modes of cracking state or non-cracking state; therefore, the numerical simulation method for hydraulic fracture generation and expansion is combined with the Biao consolidation theory to push and establish a finite element numerical simulation model for describing the hydraulic fracture generation and expansion process, so that the hydraulic fracture generation and expansion process is comprehensively and efficiently predicted and judged, and the operation safety of the tailing dam is ensured.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
The technical solution of the present invention is further described in detail by the accompanying drawings and embodiments.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a schematic flow chart of a data simulation method for crack generation and expansion of an initial dam of a tailing dam provided by the invention.
Detailed Description
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 only a part of the embodiments of the present invention, and not all of the embodiments. 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.
Fig. 1 is a schematic flow chart of a data simulation method for crack occurrence and propagation in an initial dam of a tailing dam according to an embodiment of the present invention. The numerical simulation method for hydraulic fracture generation and expansion comprises the following steps:
step S1, constructing a finite element equation about the tailing dam, and calculating to obtain displacement, stress and pore pressure field corresponding to hydraulic fracture about the tailing dam;
step S2, calculating the main stress parameter related to the hydraulic fracture according to the displacement, stress and pore pressure field of the hydraulic fracture;
step S3, according to the main stress parameter of the hydraulic fracture, the tensile strength of the tailing dam is judged;
and step S4, according to the pulling rope strength judgment processing result, carrying out different mode processing aiming at the condition that the tailing dam is in a cracking state or in an uncracked state.
Preferably, in the step S1, constructing finite element equations about the tailings dam, and calculating the displacement, stress and pore pressure field corresponding to the hydraulic fracture of the tailings dam from the finite element equations,
s101, determining a control equation set of a proportional-integral theory about the tailing dam according to a calculation model about fluid-solid coupling effective stress;
step S102, carrying out transformation processing on the control equation set of the Biaoh consolidation theory to obtain a finite element format equation set of a two-dimensional Biaoh consolidation theory;
and S103, solving a finite element format equation set of the two-dimensional Biao consolidation theory to calculate and obtain the displacement, the stress and the pore pressure field of the hydraulic fracture.
Preferably, in step S101, the determination of the theoretical control equation set of the doot consolidation for the tailing dam specifically includes,
determining the system of the control equation of the Biao consolidation theory, which is composed of the following formulas (1) to (5) together, according to a calculation model of the fluid-solid coupling effective stress
σij'＝σij-qδij＝Hijklεkl (2)
εkl＝1/2(uk,l+ul,k) (3)
In the above formulas (1) to (5), i, j are directions of the cauchy coordinate system; sigmaij,σij' Total stress and effective stress, respectively; fiA volumetric force component in the i direction; q is pore water pressure; deltaijKnown as the Kronecker (Kronecker) notation; hijklIs the constitutive relation tensor of the soil; epsilonijIs the strain tensor; u. ofijIs the displacement component; v. ofiIs the pore water flow velocity component; theta is the volume strain of soil; t is time; k is a radical ofijIs the permeability coefficient; gamma raywIs the volume weight of water.
Preferably, in step S102, the transformation processing is performed on the system of equations for controlling the proportional-integral theory to obtain the finite element format equation of the two-dimensional proportional-integral theory specifically includes,
the control equation set of the Biaoo consolidation theory is transformed to obtain a finite element format equation set of the two-dimensional Biaoo consolidation theory which is formed by the following formulas (6) to (10) together
In the above formulas (6) to (10), { δ }eRepresenting unit node displacement; { p }eRepresenting the pore water pressure of the unit node;
representing cell edge interface forces; { F } represents the unit volume force; v isnRepresenting cell boundary flow rate; [ K ]]Representing a permeability coefficient matrix; { K }eIs a cell stiffness matrix; [ K ]c]eCoupling a matrix for the cells; [ K ]s]eIs a unit seepage matrix; n is a radical of1...N4Is a unit shape function.
Preferably, in step S2, calculating a principal stress parameter for the hydraulic fracture based on the displacement, stress and pore pressure field of the hydraulic fracture specifically includes,
step S201, constructing a stress-strain relation of the compacted clay in the tailing dam in a stretching state according to the displacement, the stress and the pore pressure field of the hydraulic fracture;
step S202, determining a relation equation between stress and strain of the compacted clay in a stretching state according to the stress-strain relation of the compacted clay in the stretching state;
step S203, solving a relation equation between the stress and the strain of the compacted clay in the tensile state to calculate the minimum effective principal stress of the compacted clay as the principal stress parameter.
Preferably, in this step S202, determining the equation of relationship between the stress and strain of the compacted clay in tension from the relationship of stress and strain of the compacted clay in tension specifically comprises,
from the relationship of stress and strain in the compacted clay in a stretched state, an equation of relationship between stress and strain in a stretched state of the compacted clay of the following formula (11) is determined
In the above formula (11), ε represents a stress; σ represents strain; e represents the initial elastic modulus of the stretching section; mu.s0The initial tangent elastic modulus reduction coefficient is used for describing the initial damage degree of the soil body before the soil body bears the tensile stress; epsilonfIs the peak strain; a. the1And B1Is a material constant, and
and
preferably, in the step S3, the determining the tensile strength of the tailings dam according to the main stress parameter of the hydraulic fracture specifically includes,
step S301A, acquiring the minimum effective main stress of the tailing dam about the compacted clay from the main stress parameter of the hydraulic fracture;
step S302A, the minimum effective principal stress is compared to the tensile strength of the compacted clay.
Preferably, in the step S3, the determining the tensile strength of the tailings dam according to the main stress parameter of the hydraulic fracture specifically includes,
step S301B, obtaining the maximum main stress FM and the minimum main stress F of the hydraulic fracture according to the main stress parameter of the hydraulic fracturem；
Step S302B, calculating the hydraulic fracture pressure u according to the following formula (12)
In the formula (12), u is hydraulic fracture pressure, rho is compacted dry density of the tailing dam, delta is tensile strength of soil, and F isMMaximum principal stress for the hydraulic fracture, FmIs the minimum principal stress of the hydraulic fracture;
step S303B, calculating the tensile strength k of the tailing dam according to the following formula (13)
In the formula (13), k is tensile strength of the tailing dam, u is hydraulic fracture pressure, c is cohesion,
the internal friction angle is set, S is the contact area of the tailing dam and water, and S is the cross-sectional area of the tailing dam;
through the steps S301B-S303B, the tensile strength of the tailing dam is obtained, and further the tensile strength of the tailing dam can be judged to obtain the judgment result.
Preferably, in the step S4, the performing different modes of processing for the tailings dam in a cracked state or in an uncracked state according to the result of the pulling rope strength judging processing specifically includes,
step S401, if the result of the stay cord strength judgment processing indicates that the minimum effective principal stress is smaller than the tensile strength of the compacted clay, performing corresponding mode processing on the tailings dam in an uncracked state;
step S402, if the result of the stay cord strength judgment processing indicates that the minimum effective principal stress is greater than or equal to the tensile strength of the compacted clay, corresponding mode processing is carried out on the tailing dam in a cracking state.
Preferably, in step S401, if the result of the rope pulling strength determining process indicates that the minimum effective principal stress is smaller than the tensile strength of the compacted clay, the performing the corresponding mode process for the tailings dam in the uncracked state specifically includes,
if the result of the stay cord strength judgment processing indicates that the minimum effective principal stress is less than the tensile strength of the compacted clay, repeating the steps S1-S4 for the current tailings dam to realize the next load step numerical simulation.
Preferably, in step S402, if the result of the rope pulling strength determining process indicates that the minimum effective principal stress is greater than or equal to the tensile strength of the compacted clay, the performing of the corresponding mode process for the cracked state of the tailing dam specifically includes,
s4021, calculating and obtaining a crack direction corresponding to the crack state of the tailing dam aiming at the current crack state of the tailing dam;
step S4022, according to the crack direction, unit stiffness correction processing is carried out on the current crack state of the tailing dam so as to determine a corresponding step length correction mode;
and S4023, setting an initial pore pressure field and an external load corresponding to the current cracking state of the tailing dam according to the step length correction mode.
Preferably, in the step S4021, calculating and obtaining a crack direction corresponding to the crack state of the tailing dam according to the current crack state includes,
s40211, determining the anisotropic characteristics of a cracked soil body in terms of stress-strain relationship and permeability characteristics according to the cracking state of the tailing dam at present;
step S40212, determining a relation curve between the permeability coefficient and the normal effective stress according to the anisotropic characteristics of the cracked soil body in terms of stress-strain relation and permeability characteristic, and calculating to obtain the crack direction corresponding to the cracked state.
From the content of the embodiment, the numerical simulation method for hydraulic fracture occurrence and expansion is pushed to and establishes a finite element numerical simulation model for describing the hydraulic fracture occurrence and expansion process by combining the theory of the Biao consolidation, so that the hydraulic fracture occurrence and expansion process is comprehensively and efficiently predicted and judged, and the operation safety of the tailing dam is ensured.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such 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 such modifications and variations.
Claims (10)
1. A data simulation method for crack generation and expansion of a tailing dam initial dam is characterized by comprising the following steps of:
step S1, constructing a finite element equation about the tailing dam, and calculating to obtain displacement, stress and pore pressure field about the corresponding hydraulic fracture of the tailing dam;
step S2, calculating to obtain a main stress parameter related to the hydraulic fracture according to the displacement, the stress and the pore pressure field of the hydraulic fracture;
step S3, according to the main stress parameter of the hydraulic fracture, the tensile strength of the tailing dam is judged;
and step S4, according to the pulling rope strength judgment processing result, carrying out different mode processing aiming at the condition that the tailing dam is in a cracking state or in an uncracked state.
2. The tailing dam initial dam crack occurrence and propagation data simulation method of claim 1, characterized by:
in step S1, constructing finite element equations for the tailings dam, and calculating the displacement, stress and pore pressure field corresponding to the hydraulic fracture of the tailings dam from the finite element equations,
step S101, determining a control equation set of a proportional-integral theory about the tailing dam according to a calculation model about fluid-solid coupling effective stress;
step S102, carrying out transformation processing on the control equation set of the Biaoh consolidation theory to obtain a finite element format equation set of a two-dimensional Biaoh consolidation theory;
and S103, solving a finite element format equation set of the two-dimensional Biao consolidation theory to calculate and obtain the displacement, the stress and the pore pressure field of the hydraulic fracture.
3. The tailing dam initial dam crack occurrence and propagation data simulation method of claim 2, characterized by:
in the step S101, according to the calculation model of the fluid-solid coupling effective stress, the specific control equation set of the proportion consolidation theory about the tailing dam is determined to comprise,
determining the system of the control equation of the Biao consolidation theory, which is composed of the following formulas (1) to (5) together, according to a calculation model about the effective stress of fluid-solid coupling
σij'＝σij-qδij＝Hijklεkl (2)
εkl＝1/2(uk,l+ul,k) (3)
In the above formulas (1) to (5), i, j are directions of the cauchy coordinate system; sigmaij,σij' Total stress and effective stress, respectively; fiA volumetric force component in the i direction; q is pore water pressure; deltaijKnown as the Kronecker (Kronecker) notation; hijklIs the constitutive relation tensor of the soil; epsilonijIs the strain tensor; u. ofijIs the displacement component; v. ofiIs the pore water flow velocity component; theta is the volume strain of soil; t is time; k is a radical ofijIs the permeability coefficient; gamma raywIs the volume weight of water;
alternatively, the first and second electrodes may be,
in step S102, transforming the proportional control equation set to obtain a finite element format equation of the two-dimensional proportional control theory specifically includes,
transforming the control equation set of the Biaoo consolidation theory to obtain a finite element format equation set of the two-dimensional Biaoo consolidation theory, which is formed by the following formulas (6) to (10)
In the above formulas (6) to (10), { δ }eRepresenting unit node displacement; { p }eRepresenting the pore water pressure of the unit node;
representing cell edge interface forces; { F } represents the unit volume force; v isnRepresenting cell boundary flow rate; [ K ]]Representing a permeability coefficient matrix; { K }eIs a cell stiffness matrix; [ K ]c]eCoupling a matrix for the cells; [ K ]s]eIs a unit seepage matrix; n is a radical of1...N4Is in the form of a unitAnd (4) counting.
4. The tailing dam initial dam crack occurrence and propagation data simulation method of claim 1, characterized by:
in step S2, the calculating of the principal stress parameter related to the hydraulic fracture according to the displacement, stress and pore pressure field of the hydraulic fracture specifically includes,
step S201, constructing a stress-strain relation of the compacted clay in the tailing dam in a stretching state according to the displacement, the stress and the pore pressure field of the hydraulic fracture;
step S202, determining a relation equation between stress and strain of the compacted clay in a stretching state according to the stress-strain relation of the compacted clay in the stretching state;
step S203, solving a relation equation between the stress and the strain of the compacted clay in the tensile state to calculate the minimum effective principal stress of the compacted clay, wherein the minimum effective principal stress is used as the principal stress parameter.
5. The tailing dam initial dam crack occurrence and propagation data simulation method of claim 4, characterized by comprising the following steps:
in the step S202, determining the relation equation between the stress and the strain of the compacted clay in the tensile state according to the relation between the stress and the strain of the compacted clay in the tensile state specifically comprises,
determining the relation equation between the stress and the strain of the compacted clay in the tensile state according to the stress-strain relation of the compacted clay in the tensile state in the following formula (11)
In the above formula (11), ε represents a stress; σ represents strain; e represents the initial elastic modulus of the stretching section; mu.s0The initial tangent elastic modulus reduction coefficient is used for describing the initial damage degree of the soil body before the soil body bears the tensile stress; epsilonfIs the peak strain; a. the1And B1Is a material constant, and
and
6. the tailing dam initial dam crack occurrence and propagation data simulation method of claim 1, characterized by:
in the step S3, the determining the tensile strength of the tailing dam according to the main stress parameter of the hydraulic fracture specifically includes,
step S301A, acquiring the minimum effective main stress of the tailing dam about the compacted clay from the main stress parameters of the hydraulic fracture;
step S302A, comparing the minimum effective principal stress with the tensile strength of the compacted clay;
alternatively, the first and second electrodes may be,
in the step S3, the determining the tensile strength of the tailing dam according to the main stress parameter of the hydraulic fracture specifically includes,
step S301B, obtaining the maximum main stress F of the hydraulic fracture according to the main stress parameter of the hydraulic fractureMAnd minimum principal stress Fm；
Step S302B, calculating the hydraulic fracture pressure u according to the following formula (12)
In the formula (12), u is hydraulic fracture pressure, rho is compacted dry density of the tailing dam, delta is tensile strength of soil, and F isMIs the maximum principal stress of the hydraulic fracture, FmA minimum principal stress for the hydraulic fracture;
step S303B, calculating the tensile strength k of the tailing dam according to the following formula (13)
In the formula (13), k is tensile strength of the tailing dam, u is hydraulic fracture pressure, c is cohesion,
the internal friction angle is S, the contact area of the tailing dam and water is S, and the cross-sectional area of the tailing dam is S.
7. The tailing dam initial dam crack occurrence and propagation data simulation method of claim 6, characterized by:
in step S4, the specific steps of performing different modes of processing on the tailings dam in a cracked state or an uncracked state according to the result of the pulling rope strength determination processing include,
step S401, if the result of the stay cord strength judgment processing indicates that the minimum effective principal stress is smaller than the tensile strength of the compacted clay, performing corresponding mode processing on the tailings dam in an uncracked state;
step S402, if the result of the stay cord strength judgment processing indicates that the minimum effective principal stress is greater than or equal to the tensile strength of the compacted clay, performing corresponding mode processing on the tailing dam in a cracking state.
8. The tailing dam initial dam crack occurrence and propagation data simulation method of claim 7, characterized by:
in the step S401, if the result of the judgment processing of the strength of the pulling rope indicates that the minimum effective principal stress is smaller than the tensile strength of the compacted clay, the corresponding mode processing specifically includes that the tailings dam is in an uncracked state,
and if the result of the stay cord strength judgment processing indicates that the minimum effective principal stress is smaller than the tensile strength of the compacted clay, repeating the steps S1-S4 aiming at the current tailing dam so as to realize the next load step numerical simulation.
9. The tailing dam initial dam crack occurrence and propagation data simulation method of claim 7, characterized by:
in the step S402, if the result of the judgment processing of the rope pulling strength indicates that the minimum effective principal stress is greater than or equal to the tensile strength of the compacted clay, performing corresponding mode processing on the tailing dam in a cracking state specifically includes,
s4021, calculating and obtaining a crack direction corresponding to a crack state of the tailing dam aiming at the current crack state;
step S4022, according to the crack direction, unit stiffness correction processing is carried out on the current crack state of the tailing dam so as to determine a corresponding step length correction mode;
and S4023, setting an initial pore pressure field and an external load corresponding to the current cracking state of the tailing dam according to the step length correction mode.
10. The tailing dam initial dam crack occurrence and propagation data simulation method of claim 9, characterized by:
in the step S4021, specifically calculating and obtaining a crack direction corresponding to the crack state of the tailing dam according to the current crack state of the tailing dam,
s40211, determining the anisotropic characteristics of a cracked soil body in terms of stress-strain relationship and permeability characteristics according to the cracking state of the tailing dam at present;
step S40212, determining a relation curve between a permeability coefficient and normal effective stress according to the anisotropic characteristics of the cracked soil body in terms of stress-strain relation and permeability characteristic, and calculating to obtain a crack direction corresponding to the cracked state.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110068826.1A CN112668229B (en) | 2021-01-18 | 2021-01-18 | Data simulation method for crack generation and expansion of initial dam of tailing dam |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110068826.1A CN112668229B (en) | 2021-01-18 | 2021-01-18 | Data simulation method for crack generation and expansion of initial dam of tailing dam |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112668229A true CN112668229A (en) | 2021-04-16 |
| CN112668229B CN112668229B (en) | 2022-02-22 |
Family
ID=75415672
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110068826.1A Active CN112668229B (en) | 2021-01-18 | 2021-01-18 | Data simulation method for crack generation and expansion of initial dam of tailing dam |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112668229B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102799713A (en) * | 2012-06-26 | 2012-11-28 | 武汉大学 | Numerical simulation method for hydraulic fracture of core wall of rock-fill dam |
| US20130090902A1 (en) * | 2010-06-28 | 2013-04-11 | Yao Yao | Method and System for Modeling Fractures in Ductile Rock |
| CN111695285A (en) * | 2020-06-17 | 2020-09-22 | 大连海事大学 | Anisotropic rock mass stress-damage-seepage coupling numerical simulation method |
-
2021
- 2021-01-18 CN CN202110068826.1A patent/CN112668229B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130090902A1 (en) * | 2010-06-28 | 2013-04-11 | Yao Yao | Method and System for Modeling Fractures in Ductile Rock |
| CN102799713A (en) * | 2012-06-26 | 2012-11-28 | 武汉大学 | Numerical simulation method for hydraulic fracture of core wall of rock-fill dam |
| CN111695285A (en) * | 2020-06-17 | 2020-09-22 | 大连海事大学 | Anisotropic rock mass stress-damage-seepage coupling numerical simulation method |
Non-Patent Citations (2)
| Title |
|---|
| 李全明 等: "土石坝水力劈裂发生过程的有限元数值模拟", 《岩土工程学报》 * |
| 李全明 等: "尾矿库在线监测系统中位移数据分析方法探讨", 《中国安全生产科学技术》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112668229B (en) | 2022-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chatterjee et al. | Numerical simulations of pipe–soil interaction during large lateral movements on clay | |
| Paavilainen et al. | Parameter effects on simulated ice rubbling forces on a wide sloping structure | |
| CN111144016B (en) | Multi-axis stress state concrete variable four-parameter damage structure | |
| Zheng et al. | Arch-dam crack deformation monitoring hybrid model based on XFEM | |
| Caselunghe et al. | Structural element approaches for soil-structure interaction | |
| Haiderali et al. | Cyclic lateral loading of monopiles for offshore wind turbines | |
| CN114969952A (en) | Building collapse risk assessment method and device, computer equipment and storage medium | |
| Aron et al. | Structural element approaches for soil-structure interaction | |
| CN112668229B (en) | Data simulation method for crack generation and expansion of initial dam of tailing dam | |
| Gong et al. | The simulation of high compressive stress and extrusion phenomenon for concrete face slabs in CFRDs under strong seismic loads | |
| Lemmen et al. | The influence of foundation stiffness on the behaviour of surface strip foundations on sand | |
| CN111414656A (en) | Segmental casting method and analysis method for large-volume radiation-proof concrete wall | |
| Yao et al. | An SBFEM-Based model for hydraulic fracturing in quasi-brittle materials | |
| Hori et al. | Model test and consolidation analysis of failure of a loose sandy embankment dam during seepage | |
| Schlicke et al. | Computational prediction of restraint-induced macrocrack patterns in concrete walls | |
| Yin | A constitutive model with two yield surfaces for soils | |
| Fink et al. | Push-out test parametric simulation study of a new sheet-type shear connector | |
| Song et al. | Seismic damage analysis of the outlet piers of arch dams using the finite element sub-model method | |
| Zahra et al. | Experimental investigation of dry joint surface and closure characteristics of interlocking blocks under compression | |
| Grimstad et al. | Stability analyses of quick clay using FEM and an anisotropic strain softening model with internal length scale | |
| Adachi et al. | An elasto-viscoplastic constitutive model with strain-softening for soft sedimentary rocks | |
| Trindade | Numerical modeling of the post-cracking behavior of SFRC and its application on design of beams according to fib Model Code 2010. | |
| CN113312697A (en) | High-pressure grouting anti-slip stability prediction method for riparian retaining wall | |
| CN112948993A (en) | Fretting fatigue life prediction method considering abrasion influence under spectral load | |
| CN113627052A (en) | Rock-fill dam rheological numerical simulation method considering hydraulic coupling effect |
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 |