CN115659743A - Method and device for constructing elastic-plastic damage model of stratified rock mass - Google Patents
Method and device for constructing elastic-plastic damage model of stratified rock mass Download PDFInfo
- Publication number
- CN115659743A CN115659743A CN202211314180.1A CN202211314180A CN115659743A CN 115659743 A CN115659743 A CN 115659743A CN 202211314180 A CN202211314180 A CN 202211314180A CN 115659743 A CN115659743 A CN 115659743A
- Authority
- CN
- China
- Prior art keywords
- model
- damage
- elastic
- plastic
- rock mass
- 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.)
- Pending
Links
Images
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The application discloses a method and a device for constructing a layered rock elastic-plastic damage model, wherein the method comprises the following steps: based on FLAC 3D An embedded transverse isotropic elastic model; combining a preset yield criterion, introducing an independent variable taking equivalent plastic strain as a damage variable, and performing dynamic plastic damage correction on the cohesive force and the internal friction angle of the rock mass through a loss variable; substituting the corrected cohesive force and the internal friction angle into the updated yield criterion; and calculating the effective stress expression after the plasticity correction, and compiling the elastoplasticity damage constitutive model on a preset platform according to the corrected effective stress expression. The embodiment of the application can be used for judging the plastic zone and calculating the plastic stage, and can truly reflect the transverse isotropic deformation characteristic of the stratified rock mass by considering the directionality of rock mass parameters, so that the method is simple, practical and easy to operate.
Description
Technical Field
The application relates to the technical field of rock and soil, in particular to a method and a device for constructing a layered rock elastic-plastic damage model.
Background
In the related technology, a transverse isotropic elastic model and a pervasive joint model can be embedded into a stratified rock body, wherein the transverse isotropic elastic model can be used for simple linear elastic calculation, and the pervasive joint model can adopt a Mohr-Coulomb model to establish respective strength yield criterion for a rock block and a structural surface, so that transverse isotropy and structural surface characteristics can be considered.
However, in the related art, only simple linear elasticity calculation can be performed, the plastic region determination and the calculation of the plastic stage cannot be performed, and the calculation is isotropic, so that the directionality of rock parameters is not considered, and the transverse isotropic deformation characteristic of the stratified rock mass cannot be truly reflected, and a solution is needed.
Disclosure of Invention
The present application is based on the inventors' recognition and problem that:
the stratified rock mass widely exists in tunnel engineering, the physical and mechanical properties of the stratified rock mass tend to be complicated due to long-term geological historical evolution process and tectonic motion, and any constitutive model cannot perfectly fit the stress-strain relation of the rock mass.
Model simplification is an important mode for researching constitutive models of rock masses, the stratified rock masses are generally regarded as transverse isotropic models, and common finite difference numerical software FLAC 3D A transverse isotropic elastic model and a pervasive joint model are embedded in software aiming at a stratified rock mass, transverse isotropy and structural surface characteristics can be considered, but the software has slight defects in application, the transverse isotropic elastic model only can carry out simple linear elastic calculation and cannot carry out plastic region judgment and plastic stage calculation, and the pervasive joint model adopts a Mohr-Coulomb model to respectively aim at rock masses and knotsThe respective strength yield criterion of the tectonic surface is established, but the calculation is isotropic, the directionality of rock mass parameters is not considered, and the transverse isotropic deformation characteristics of the stratified rock mass cannot be truly reflected.
The stratified rock body can be damaged in the plastic deformation stage, the dynamic change of the cohesive force and the internal friction angle of the rock body can be influenced, and FLAC 3D The theoretical stress-strain curve of the Mohr-Coulomb model in the software is straight in the later stage and cannot reflect the strain softening stage.
The application provides a method and a device for constructing an elastic-plastic damage model of a layered rock mass, which are used for solving the problems that in the related technology, only simple linear elastic calculation can be performed, plastic zone discrimination and plastic stage calculation cannot be performed, isotropy is performed in the calculation, the directionality of rock mass parameters is not considered, and the transverse isotropy deformation characteristics of the layered rock mass cannot be truly reflected.
The embodiment of the first aspect of the application provides a method for constructing a lamellar rock elastic-plastic damage model, which comprises the following steps: based on FLAC 3D Determining input parameters of a constitutive model of the stratified rock mass by using the embedded transverse isotropic elastic model; determining a stress conversion matrix of a local coordinate system and a whole coordinate system in the stratified rock mass constitutive model, and defining an elastic stiffness matrix under the local coordinate system; calculating an elastic stiffness matrix under the overall coordinate system through the stress conversion matrix; calculating an imaginary elastic stress increment according to the total strain increment and the elastic stiffness matrix; correcting the stress which accords with the elastic-plastic model in the next step according to the plastic strain; judging whether the material meets a preset yield condition or not based on a preset yield criterion, if so, introducing an independent variable taking equivalent plastic strain as a damage variable, and performing dynamic plastic damage correction on the cohesive force and the internal friction angle of the rock mass through the damage variable; substituting the corrected cohesive force and the internal friction angle into the updated yield criterion; and calculating the effective stress expression after the plasticity correction, and compiling an elastoplasticity damage constitutive model on a preset platform according to the corrected effective stress expression.
Optionally, in one embodiment of the present application, the present application implementsThe method of example further comprises: exporting a dll dynamic link file from the written elastoplasticity damage constitutive model to obtain a data for the FLAC 3D And calling the improved constitutive model file.
Optionally, in an embodiment of the present application, the input parameters include an elastic modulus and a poisson's ratio parallel to the isotropic plane, an elastic modulus perpendicular to the isotropic plane, a poisson's ratio, and a shear modulus.
Optionally, in an embodiment of the present application, the effective stress expression is:
wherein λ is t Is a stress correction factor in the tensile failure criterion.
Optionally, in an embodiment of the present application, the preset yield criterion is a Mohr-Coulomb strength yield criterion.
The embodiment of the second aspect of the present application provides a layered rock mass elastic-plastic damage model building device, including: a first determination module to determine based on FLAC 3D Determining input parameters of a constitutive model of the stratified rock mass by using the embedded transverse isotropic elastic model; the second determination module is used for determining a stress conversion matrix of a local coordinate system and a whole coordinate system in the layered rock mass constitutive model and defining an elastic stiffness matrix under the local coordinate system; the first calculation module is used for calculating an elastic stiffness matrix under the integral coordinate system through the stress conversion matrix; the second calculation module is used for calculating an imaginary elastic stress increment according to the total strain increment and the elastic stiffness matrix; the correcting module is used for correcting the stress which accords with the elastic-plastic model in the next step according to the plastic strain; the judging module is used for judging whether the material meets a preset yield condition or not based on a preset yield criterion, if the material yields, introducing an independent variable taking equivalent plastic strain as a damage variable, and performing dynamic plastic damage correction on the cohesive force and the internal friction angle of the rock mass through the damage variable; a processing module for correcting the cohesion and the internal friction angleBringing in updated yield criteria; and the construction module is used for calculating the effective stress expression after the plasticity correction so as to compile an elastoplasticity damage constitutive model on a preset platform according to the corrected effective stress expression.
Optionally, in an embodiment of the present application, the apparatus in the embodiment of the present application further includes: a exporting module for exporting the written elastoplasticity damage constitutive model to a dll dynamic link file to obtain a data for the FLAC 3D And calling the improved constitutive model file.
Optionally, in an embodiment of the present application, the input parameters include an elastic modulus and a poisson's ratio parallel to the isotropic plane, an elastic modulus perpendicular to the isotropic plane, a poisson's ratio, and a shear modulus.
Optionally, in an embodiment of the present application, the effective stress expression is:
wherein λ is t Is a stress correction factor in the tensile failure criterion.
Optionally, in an embodiment of the present application, the preset yield criterion is a Mohr-Coulomb strength yield criterion.
An embodiment of a third aspect of the present application provides an electronic device, including: the device comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor executes the program to realize the stratified rock mass elastic-plastic damage model building method in the embodiment.
A fourth aspect of the present application provides a computer-readable storage medium, which stores a computer program, and when the program is executed by a processor, the method for constructing the elastic-plastic damage model of the stratified rock mass is implemented.
The embodiment of the application can be based on FLAC 3D An embedded transverse isotropic elastic model is combined with a preset yield criterion and introduces equivalent plastic strain asThe method comprises the steps of damaging independent variables of variables, dynamically correcting plastic damage of a rock mass cohesive force and an internal friction angle through the damaged variables, bringing the corrected cohesive force and the internal friction angle into an updated yield criterion, calculating an effective stress expression after plastic correction, and compiling an elastoplastic damage constitutive model on a platform according to the corrected effective stress expression, so that the calculation of plastic region discrimination and plastic stage can be carried out, the directionality of rock mass parameters is considered, the transverse isotropic deformation characteristics of the stratified rock mass can be truly reflected, and the method is simple, practical and easy to operate. Therefore, the problems that in the related technology, only simple linear elasticity calculation can be carried out, the plastic region judgment and the calculation of the plastic stage cannot be carried out, the calculation is isotropic, the directionality of rock mass parameters is not considered, and the transverse isotropic deformation characteristics of the stratified rock mass cannot be truly reflected are solved.
Additional aspects and advantages of the present application will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the present application.
Drawings
The foregoing and/or additional aspects and advantages of the present application will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:
FIG. 1 is a flow chart of a method for constructing a layered rock elastic-plastic damage model according to an embodiment of the present application;
FIG. 2 is a flowchart of a method for constructing an elastoplasticity damage model of a stratified rock mass according to an embodiment of the present application;
FIG. 3 is a schematic structural diagram of the method and the device for constructing the elastic-plastic damage model of the stratified rock mass according to the embodiment of the application;
fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present application.
Detailed Description
Reference will now be made in detail to embodiments of the present application, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the drawings are exemplary and intended to be used for explaining the present application and should not be construed as limiting the present application.
The method and the device for constructing the elastic-plastic damage model of the stratified rock mass according to the embodiment of the application are described below with reference to the accompanying drawings. Aiming at the problems that only simple linear elastic calculation can be carried out, the plastic region discrimination and the calculation of the plastic stage cannot be carried out, the calculation is isotropic, the directionality of rock parameters is not considered, and the transverse isotropic deformation characteristic of the stratified rock mass cannot be truly reflected in the related technology mentioned in the background technology center, the application provides a stratified rock mass elastic-plastic damage model construction method based on FLAC (FLAC-based computational modeling) 3D The embedded transverse isotropic elastic model is combined with a preset yield criterion, equivalent plastic strain is introduced to serve as an independent variable of a damage variable, dynamic plastic damage correction is carried out on the cohesive force and the internal friction angle of a rock body through the damage variable, the corrected cohesive force and the corrected internal friction angle are brought into the updated yield criterion, an effective stress expression formula after plastic correction is calculated, and an elastic-plastic damage constitutive model is compiled on a platform according to the corrected effective stress expression formula, so that the calculation of plastic region discrimination and plastic stage can be carried out, the directionality of rock body parameters is considered, the transverse isotropic deformation characteristics of the stratified rock body can be truly reflected, and the embedded transverse isotropic elastic model is simple, practical and easy to operate. Therefore, the problems that in the related technology, only simple linear elasticity calculation can be carried out, the plastic region judgment and the calculation of the plastic stage cannot be carried out, the calculation is isotropic, the directionality of rock mass parameters is not considered, and the transverse isotropic deformation characteristics of the stratified rock mass cannot be truly reflected are solved.
Specifically, fig. 1 is a schematic flow chart of a layered rock elastic-plastic damage model building method provided in an embodiment of the present application.
As shown in figure 1, the method for constructing the elastic-plastic damage model of the stratified rock mass comprises the following steps:
in step S101, based on FLAC 3D And determining the input parameters of the constitutive model of the stratified rock mass by using the embedded transverse isotropic elastic model.
It can be understood that the stratified rock mass can be approximately regarded as a transverse isotropic body based on the FLAC in the embodiment of the application 3D The embedded transverse isotropic elastic model determines the input parameters of the constitutive model of the stratified rock mass in the following steps, so that the feasibility of constructing the elastoplasticity damage model of the stratified rock mass is effectively improved.
Optionally, in an embodiment of the application, the input parameters comprise an elastic modulus and a poisson's ratio parallel to the isotropic plane, an elastic modulus perpendicular to the isotropic plane, a poisson's ratio and a shear modulus.
In practical implementation, the elastic modulus E parallel to the isotropic surface can be selected according to the embodiment of the present application 1 And poisson ratio mu 1 Elastic modulus E perpendicular to the isotropic plane 2 Poisson ratio mu 2 And shear modulus G 2 The input parameters can be obtained only through a common rock mechanics test, so that the comprehensiveness and the accuracy of the input parameters of the constitutive model of the stratified rock are improved, and meanwhile, the performability of the elastic-plastic damage model construction of the stratified rock is improved.
In step S102, a stress transformation matrix of the local coordinate system and the global coordinate system in the layered rock mass constitutive model is determined, and an elastic stiffness matrix in the local coordinate system is defined.
It can be understood that, the embodiment of the application can determine the stress conversion matrix L of the local coordinate system and the global coordinate system in the constitutive model of the stratified rock mass, and define the elastic stiffness matrix D under the local coordinate system ′ Through the conversion relation between the local coordinate system and the whole coordinate system, the intuitive setting habit of a user for the inclination angle is met, the convenience of coordinate conversion in the calculation process can be improved, and the accuracy of a stress-strain result is improved.
When the elastic equation of each isobody is transversely observed, the stress conversion matrix L converted from the overall coordinate system to the local coordinate system is as follows:
wherein, { l i ,m i ,n i I =1,2,3 is the cosine of the angle between the ith local coordinate axis and the global coordinate axis.
Further, for the convenience of model calculation and the application habit of local coordinate system in geology, the determined elastic stiffness matrix of the transverse isotropic body under the local coordinate system is as follows:
in step S103, an elastic stiffness matrix in the global coordinate system is calculated from the stress conversion matrix.
In an actual implementation process, the elastic stiffness matrix D in the global coordinate system may be calculated by using the stress transformation matrix L in the embodiment of the present application, that is, the transverse isotropic elastic stiffness matrix in the global coordinate system obtained according to the relationship between the stress and the strain in the local coordinate system and the stress and the strain in the global coordinate system is:
[D]＝[L] T [D′][L]
in step S104, a virtual elastic stress increment is calculated from the total strain increment and the elastic stiffness matrix.
It can be understood that the embodiment of the present application can be based on the total strain increment Δ ∈ i And the assumed elastic stress increment is calculated with the elastic stiffness matrix D, so that the intelligence and convenience of model development are effectively improved.
In some embodiments, when the stress state reaches the initial yield surface, the total strain increase is expressed as:
wherein, the first and the second end of the pipe are connected with each other,
in order to increase the elastic strain of the material,
is the plastic strain increment.
Wherein the increase in plastic strain is according to the plastic flow law
Can be expressed as:
wherein λ represents a stress correction factor, g represents a plastic potential function, σ i Is a stress.
In step S105, a correction is made based on the plastic strain and the stress for the next step conforming to the elastoplastic model.
It is understood that the embodiments of the present application may be based on plastic strain
And stress for next step conforming to elastoplastic model
And correcting to effectively improve the accuracy of the stress-strain result.
Wherein, according to the plastic strain pair
And stress for next step conforming to elastoplastic model
And (3) correcting, namely:
wherein the content of the first and second substances,
the hypothetical elastic test stress calculated from the total strain is shown.
In step S106, whether the material meets the preset yield condition is determined based on the preset yield criterion, and if the material yields, the equivalent plastic strain is introduced as an independent variable of the damage variable, and dynamic plastic damage correction is performed on the cohesive force and the internal friction angle of the rock mass through the damage variable.
It is understood that the embodiments of the present application can determine whether the material satisfies the yield condition based on the yield criterion in the following steps, and when the yield condition is satisfied, the equivalent plastic strain is introduced
The independent variable of the damage variable H is considered that the stratified rock body is damaged along with the expansion and penetration of internal micro-cracks in the deformation process, and the cohesive force c and the internal friction angle of the rock body are considered through the damage variable H
And (4) performing dynamic plastic damage correction to realize plastic zone discrimination and plastic stage calculation.
Further, according to the damage variable H, the internal friction angle c and the cohesive force
After dynamic correction, due to c and
the stress-strain curve can generate a strain softening phenomenon due to the reduction of the value, and the stress-strain curve is more in line with the actual situation compared with the stress-strain curve which is flat after yielding and of an ideal molar coulomb elastoplastic model.
Wherein the damage variable H used in the examples of the present application follows the equivalent plastic strain
Is shown byThe trend toward slower and slower growth, finally, toward 1.
Wherein, in one embodiment of the present application, the predetermined yield criterion is the Mohr-Coulomb strength yield criterion.
As a possible implementation manner, the embodiment of the present application may determine whether to enter the plastic stage according to the Mohr-Coulomb strength yield criterion, when plastic deformation occurs, the internal damage of the sample is continuously accumulated, and the cohesive force between the internal friction angle c and the damage variable is continuously accumulated
Collapse, cohesion c and internal friction angle as plastic deformation accumulates
The damage correction is carried out by adopting a formula (8) for describing the damage evolution process of the stratified rock mass:
wherein, c 0 、
Respectively, the initial cohesion and the internal friction angle, and H represents the damage variable.
Wherein the damage variable H can be expressed as equivalent plastic strain
Dependent variables as independent variables, namely:
wherein a is a dimensionless parameter of the model.
In some embodiments, damage does not actually occur until the plastic phase is entered, since the stratified rock mass is not damaged or can be ignored at the initial compaction phase during deformation, whereas FLAC does not actually occur 3D Middle school theoryThe stress-strain curve cannot reflect the compaction process, so that the elastoplastic damage model of the embodiment of the application is more suitable for FLAC 3D Carrying out secondary development of the constitutive model.
In step S107, the corrected cohesion and internal friction angles are substituted into the updated yield criterion.
It can be understood that the embodiment of the present application can modify the cohesion c and the internal friction angle
And introducing the updated yield criterion, thereby improving the accuracy of model development and better meeting the actual situation.
In step S108, the effective stress expression after the plastic correction is calculated, so as to compile an elastoplastic damage constitutive model on a preset platform according to the corrected effective stress expression.
It can be understood that, in the embodiment of the present application, the effective stress expression after plastic correction in the following steps may be calculated, so as to compile an elastoplastic damage constitutive model on the platform according to the corrected effective stress expression, for example, the elastoplastic damage constitutive model may be compiled on a Visual Studio 2010 platform, which effectively improves convenience and personalization level of model development.
For example, a Visual Studio 2010 platform can be selected to realize secondary development of the constitutive model of the elastic-plastic damage of the stratified rock mass, and the platform provides a C + + module which is convenient to understand and compile to compile model codes, so that secondary development of scientific researchers is facilitated, and the constitutive model with personalized functions is realized.
Wherein the effective stress expression after correcting plastic damage according to the shearing failure criterion is as follows:
wherein, the first and the second end of the pipe are connected with each other,
where θ is the shear expansion angle, λ s Is a stress correction factor in the shear failure criterion.
Wherein the effective stress expression after correcting plastic damage according to the tensile failure criterion is as follows:
wherein λ is t Is a stress correction factor in the tensile failure criterion.
Optionally, in an embodiment of the present application, the method of the embodiment of the present application further includes: exporting a dll dynamic link file from the written elastoplastic damage constitutive model to obtain a FLAC for the FLAC 3D And calling the improved constitutive model file.
In some embodiments, according to the modified effective stress expression, the embodiment of the application can implement compiling and deployment of a transverse isotropic elastoplastic damage model on a C + + plate in a Visual Studio 2010 environment, compile to generate a dll dynamic link file, and place a file source code in an FLAC 3D In the installation directory "\\ plugfiles \ cmodels", if there is no folder, it is in FLAC 3D The folder is newly established in the installation directory, and the dll dynamic link file is loaded for use through a 'model configure plug' command during numerical simulation analysis, so that the folder is provided for the FLAC 3D The called improved constitutive model file is suitable for the stress-strain process of the stratified rock mass and has a good effect.
The working principle of the embodiment of the present application is explained in detail with a specific embodiment as shown in fig. 2.
Step S201: inputting model parameters, i.e., embodiments of the present application may be based on FLAC 3D And determining input parameters of the constitutive model by using the embedded transverse isotropic elastic model.
Step S202: defining an elastic rigidity matrix and a stress conversion matrix under a local coordinate system.
Step S203: and calculating an elastic stiffness matrix under the whole coordinate system, namely calculating the elastic stiffness matrix under the whole coordinate system through a stress conversion matrix.
Step S204: and calculating the virtual elastic stress, namely calculating the virtual elastic stress increment according to the total strain increment and the elastic rigidity matrix.
Step S205: and judging whether the material is yielded or not, namely judging whether the material is yielded or not according to a yielding criterion, executing the step S206 when the material is yielded, and otherwise, executing the step S203.
Step S206: the equivalent plastic strain is calculated and the damage variable is calculated.
Step S207: and judging whether the calculation is stopped or not, and executing the step S210 when the calculation is stopped, otherwise, executing the step S208.
Step S208: correction of cohesion c and internal friction angle
Namely the cohesive force c and the internal friction angle of the rock mass through the damage variable
And (4) performing dynamic plastic damage correction.
Step S209: correcting yield criterion, i.e. viscosity c and internal friction angle after correction
Bringing in updated yield criteria.
Step S210: and calculating an effective stress expression and finishing.
Step S211: and (3) secondary development of the Visual Studio 2010, namely writing an elastoplastic damage constitutive model on a Visual Studio 2010 platform according to the effective stress expression.
In summary, the embodiment of the present application may be based on FLAC 3D Transverse isotropic elastic model codes embedded in software are combined with the Mohr-Coulomb strength yield criterion, damage variables which are equivalent plastic strain to independent variables are introduced to dynamically correct the cohesive force and the internal friction angle of the rock mass, secondary development is realized in a Visual Studio 2010 platform, the construction method of the elastic plastic damage model is simple and practical, and the elastic plastic damage model is suitable for the stress-strain process of the stratified rock mass and has the effectsBetter, simple and practical, and easy operation.
The stratified rock mass elastoplasticity damage model construction method provided by the embodiment of the application can be based on FLAC 3D The embedded transverse isotropic elastic model is combined with a preset yield criterion, equivalent plastic strain is introduced as an independent variable of a damage variable, dynamic plastic damage correction is carried out on the cohesive force and the internal friction angle of the rock body through the damage variable, the corrected cohesive force and the corrected internal friction angle are brought into the updated yield criterion, an effective stress expression after the plastic correction is calculated, and an elastic-plastic damage constitutive model is compiled on a platform according to the corrected effective stress expression, so that the plastic zone judgment and the plastic stage calculation can be carried out, the directionality of rock body parameters is considered, the transverse isotropic deformation characteristic of the stratified rock body can be truly reflected, and the transverse isotropic elastic model is simple, practical and easy to operate. Therefore, the problems that in the related technology, only simple linear elasticity calculation can be carried out, the plastic region judgment and the calculation of the plastic stage cannot be carried out, the calculation is isotropic, the directionality of rock mass parameters is not considered, and the transverse isotropic deformation characteristics of the stratified rock mass cannot be truly reflected are solved.
Next, referring to the attached drawings, the stratified rock mass elastoplasticity damage model construction device provided by the embodiment of the application is described.
Fig. 3 is a block diagram of the stratified rock mass elastoplasticity damage model building device according to the embodiment of the application.
As shown in fig. 3, the stratified rock mass elastoplastic damage model building device 10 includes: a first determination module 100, a second determination module 200, a first calculation module 300, a second calculation module 400, a correction module 500, a judgment module 600, a processing module 700, and a construction module 800.
In particular, a first determination module 100 for FLAC-based 3D And determining the input parameters of the constitutive model of the stratified rock mass by using the embedded transverse isotropic elastic model.
And a second determining module 200, configured to determine a stress transformation matrix of the local coordinate system and the global coordinate system in the layered rock constitutive model, and define an elastic stiffness matrix in the local coordinate system.
The first calculating module 300 is configured to calculate an elastic stiffness matrix in the global coordinate system through the stress transformation matrix.
And a second calculating module 400, configured to calculate the hypothetical elastic stress increment according to the total strain increment and the elastic stiffness matrix.
And the correcting module 500 is used for correcting the stress which accords with the elastic-plastic model in the next step according to the plastic strain.
The determining module 600 is configured to determine whether the material meets a preset yield condition based on a preset yield criterion, if the material yields, introduce an equivalent plastic strain as an independent variable of a damage variable, and perform dynamic plastic damage correction on the cohesive force and the internal friction angle of the rock mass through the damage variable.
A processing module 700 for substituting the corrected cohesion and internal friction angle into the updated yield criterion.
And a building module 800, configured to calculate the effective stress expression after the plasticity correction, and compile an elastoplasticity damage constitutive model on a preset platform according to the corrected effective stress expression.
Optionally, in an embodiment of the present application, the apparatus 10 of the embodiment of the present application further includes: and a derivation module.
Wherein, the export module is used for exporting dll dynamic link files from the written elastoplasticity damage constitutive model to obtain a FLAC 3D And calling the improved constitutive model file.
Optionally, in an embodiment of the application, the input parameters comprise an elastic modulus and a poisson's ratio parallel to the isotropic plane, an elastic modulus perpendicular to the isotropic plane, a poisson's ratio and a shear modulus.
Optionally, in an embodiment of the present application, the effective stress expression is:
wherein λ is t Is a stress correction factor in the tensile failure criterion.
Optionally, in one embodiment of the present application, the preset yield criterion is a Mohr-Coulomb strength yield criterion.
It should be noted that the explanation of the embodiment of the method for constructing the lamellar rock elastoplasticity damage model is also applicable to the device for constructing the lamellar rock elastoplasticity damage model of this embodiment, and is not described herein again.
According to the stratified rock mass elastoplasticity damage model construction device provided by the embodiment of the application, the device can be based on FLAC 3D The embedded transverse isotropic elastic model is combined with a preset yield criterion, equivalent plastic strain is introduced as an independent variable of a damage variable, dynamic plastic damage correction is carried out on the cohesive force and the internal friction angle of the rock body through the damage variable, the corrected cohesive force and the corrected internal friction angle are brought into the updated yield criterion, an effective stress expression after the plastic correction is calculated, and an elastic-plastic damage constitutive model is compiled on a platform according to the corrected effective stress expression, so that the plastic zone judgment and the plastic stage calculation can be carried out, the directionality of rock body parameters is considered, the transverse isotropic deformation characteristic of the stratified rock body can be truly reflected, and the transverse isotropic elastic model is simple, practical and easy to operate. Therefore, the problems that in the related technology, only simple linear elasticity calculation can be performed, the plastic zone discrimination and the calculation of the plastic stage cannot be performed, the calculation is isotropic, the directionality of rock mass parameters is not considered, and the transverse isotropic deformation characteristics of the stratified rock mass cannot be truly reflected are solved.
Fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present application. The electronic device may include:
The processor 402 executes the program to implement the method for constructing the elastic-plastic damage model of the stratified rock mass provided in the above embodiment.
Further, the electronic device further includes:
a communication interface 403 for communication between the memory 401 and the processor 402.
A memory 401 for storing computer programs executable on the processor 402.
If the memory 401, the processor 402 and the communication interface 403 are implemented independently, the communication interface 403, the memory 401 and the processor 402 may be connected to each other through a bus and perform communication with each other. The bus may be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, an Extended ISA (EISA) bus, or the like. The bus may be divided into an address bus, a data bus, a control bus, etc. For ease of illustration, only one thick line is shown in FIG. 4, but this does not indicate only one bus or one type of bus.
Alternatively, in practical implementation, if the memory 401, the processor 402, and the communication interface 403 are integrated on one chip, the memory 401, the processor 402, and the communication interface 403 may complete mutual communication through an internal interface.
The processor 402 may be a Central Processing Unit (CPU), an Application Specific Integrated Circuit (ASIC), or one or more Integrated circuits configured to implement embodiments of the present Application.
The embodiment of the application also provides a computer-readable storage medium, on which a computer program is stored, and when the program is executed by a processor, the method for constructing the elastic-plastic damage model of the stratified rock mass is realized.
In the description herein, reference to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the application. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or N embodiments or examples. Furthermore, various embodiments or examples and features of different embodiments or examples described in this specification can be combined and combined by one skilled in the art without contradiction.
Furthermore, the terms "first", "second" and "first" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one such feature. In the description of the present application, "N" means at least two, e.g., two, three, etc., unless specifically limited otherwise.
Any process or method descriptions in flow charts or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or N executable instructions for implementing steps of a custom logic function or process, and alternate implementations are included within the scope of the preferred embodiment of the present application in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the embodiments of the present application.
The logic and/or steps represented in the flowcharts or otherwise described herein, such as an ordered listing of executable instructions that can be considered to implement logical functions, can be embodied in any computer-readable medium for use by or in connection with an instruction execution system, apparatus, or device, such as a computer-based system, processor-containing system, or other system that can fetch the instructions from the instruction execution system, apparatus, or device and execute the instructions. For the purposes of this description, a "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of the computer-readable medium would include the following: an electrical connection (electronic device) having one or N wires, a portable computer diskette (magnetic device), a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). Further, the computer readable medium could even be paper or another suitable medium upon which the program is printed, as the program can be electronically captured, via optical scanning of the paper or other medium, then compiled, interpreted or otherwise processed in a suitable manner if necessary, and then stored in a computer memory.
It should be understood that portions of the present application may be implemented in hardware, software, firmware, or a combination thereof. In the above embodiments, the N steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. If implemented in hardware, as in another embodiment, any one or combination of the following technologies, which are well known in the art, may be used: a discrete logic circuit having a logic gate circuit for implementing a logic function on a data signal, an application specific integrated circuit having an appropriate combinational logic gate circuit, a Programmable Gate Array (PGA), a Field Programmable Gate Array (FPGA), or the like.
It will be understood by those skilled in the art that all or part of the steps carried by the method for implementing the above embodiments may be implemented by hardware related to instructions of a program, which may be stored in a computer readable storage medium, and when the program is executed, the program includes one or a combination of the steps of the method embodiments.
In addition, functional units in the embodiments of the present application may be integrated into one processing module, or each unit may exist alone physically, or two or more units are integrated into one module. The integrated module can be realized in a hardware mode, and can also be realized in a software functional module mode. The integrated module, if implemented in the form of a software functional module and sold or used as a stand-alone product, may also be stored in a computer readable storage medium.
The storage medium mentioned above may be a read-only memory, a magnetic or optical disk, etc. While embodiments of the present application have been shown and described above, it will be understood that the above embodiments are exemplary and should not be construed as limiting the present application and that changes, modifications, substitutions and alterations in the above embodiments may be made by those of ordinary skill in the art within the scope of the present application.
Claims (10)
1. A method for constructing a lamellar rock elastic-plastic damage model is characterized by comprising the following steps:
based on FLAC 3D Determining input parameters of a constitutive model of the stratified rock mass by using the embedded transverse isotropic elastic model;
determining a stress conversion matrix of a local coordinate system and a whole coordinate system in the stratified rock mass constitutive model, and defining an elastic stiffness matrix under the local coordinate system;
calculating an elastic stiffness matrix under the overall coordinate system through the stress conversion matrix;
calculating an imaginary elastic stress increment according to the total strain increment and the elastic stiffness matrix;
correcting the stress which accords with the elastic-plastic model in the next step according to the plastic strain;
judging whether the material meets a preset yield condition or not based on a preset yield criterion, if so, introducing an independent variable taking equivalent plastic strain as a damage variable, and performing dynamic plastic damage correction on the cohesive force and the internal friction angle of the rock mass through the damage variable;
substituting the corrected cohesive force and the internal friction angle into the updated yield criterion; and
and calculating the effective stress expression after the plasticity correction, and compiling an elastoplasticity damage constitutive model on a preset platform according to the corrected effective stress expression.
2. The method of claim 1, further comprising:
exporting a dll dynamic link file from the written elastoplasticity damage constitutive model to obtain a data for the FLAC 3D And calling the improved constitutive model file.
3. The method of claim 1, wherein the input parameters include elastic modulus and poisson's ratio parallel to the isotropic face, elastic modulus perpendicular to the isotropic face, poisson's ratio, and shear modulus.
4. The method of claim 1, wherein the effective stress expression is:
wherein λ is t Is a stress correction factor in the tensile failure criterion.
5. The method according to any one of claims 1 to 4, wherein the preset yield criterion is the Mohr-Coulomb strength yield criterion.
6. The utility model provides a stratiform rock mass elastoplasticity damage model building device which characterized in that includes:
a first determining module to determine based on FLAC 3D Determining input parameters of a constitutive model of the stratified rock mass by using the embedded transverse isotropic elastic model;
the second determination module is used for determining a stress conversion matrix of a local coordinate system and a whole coordinate system in the layered rock mass constitutive model and defining an elastic stiffness matrix under the local coordinate system;
the first calculation module is used for calculating an elastic stiffness matrix under the integral coordinate system through the stress conversion matrix;
the second calculation module is used for calculating an imaginary elastic stress increment according to the total strain increment and the elastic stiffness matrix;
the correcting module is used for correcting the stress which accords with the elastic-plastic model in the next step according to the plastic strain;
the judging module is used for judging whether the material meets a preset yield condition or not based on a preset yield criterion, if the material yields, introducing an independent variable taking equivalent plastic strain as a damage variable, and performing dynamic plastic damage correction on the cohesive force and the internal friction angle of the rock mass through the damage variable;
the processing module is used for substituting the corrected cohesive force and the internal friction angle into the updated yield criterion; and
and the construction module is used for calculating the effective stress expression after the plasticity correction so as to compile an elastoplasticity damage constitutive model on a preset platform according to the corrected effective stress expression.
7. The apparatus of claim 6, further comprising:
an export module for exporting the written elastoplastic damage constitutive model to a dll dynamic link file to obtain a data for the FLAC 3D And calling the improved constitutive model file.
8. The apparatus of claim 6, wherein the input parameters include elastic modulus and Poisson's ratio parallel to the isotropic plane, elastic modulus perpendicular to the isotropic plane, poisson's ratio, and shear modulus.
9. An electronic device, comprising: a memory, a processor and a computer program stored on the memory and operable on the processor, the processor executing the program to implement the method of constructing a stratified rock mass elastoplasticity damage model as claimed in any one of claims 1 to 5.
10. A computer-readable storage medium on which a computer program is stored, the program being executable by a processor for implementing a stratified rock mass elasto-plastic damage model building method as claimed in any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211314180.1A CN115659743A (en) | 2022-10-25 | 2022-10-25 | Method and device for constructing elastic-plastic damage model of stratified rock mass |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211314180.1A CN115659743A (en) | 2022-10-25 | 2022-10-25 | Method and device for constructing elastic-plastic damage model of stratified rock mass |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115659743A true CN115659743A (en) | 2023-01-31 |
Family
ID=84991443
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211314180.1A Pending CN115659743A (en) | 2022-10-25 | 2022-10-25 | Method and device for constructing elastic-plastic damage model of stratified rock mass |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115659743A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116822330A (en) * | 2023-05-19 | 2023-09-29 | 四川大学 | Proppant elastoplastic embedding process analysis method, device, equipment and storage medium |
-
2022
- 2022-10-25 CN CN202211314180.1A patent/CN115659743A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116822330A (en) * | 2023-05-19 | 2023-09-29 | 四川大学 | Proppant elastoplastic embedding process analysis method, device, equipment and storage medium |
| CN116822330B (en) * | 2023-05-19 | 2024-02-20 | 四川大学 | Proppant elastoplastic embedding process analysis method, device, equipment and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Sivasithamparam et al. | Undrained expansion of a cylindrical cavity in clays with fabric anisotropy: theoretical solution | |
| Monforte et al. | Numerical simulation of undrained insertion problems in geotechnical engineering with the Particle Finite Element Method (PFEM) | |
| Nazem et al. | Arbitrary Lagrangian–Eulerian method for large‐strain consolidation problems | |
| Yin et al. | An anisotropic elastic–viscoplastic model for soft clays | |
| Asferg et al. | A consistent partly cracked XFEM element for cohesive crack growth | |
| Chen et al. | A family of improved yield surfaces and their application in modeling of isotropically over-consolidated clays | |
| CN111695285A (en) | Anisotropic rock mass stress-damage-seepage coupling numerical simulation method | |
| CN115659743A (en) | Method and device for constructing elastic-plastic damage model of stratified rock mass | |
| Wu et al. | A generalized nonlinear failure criterion for frictional materials | |
| Yang et al. | Modelling the plastic anisotropy of Lower Cromer Till | |
| Tian et al. | Application of the orthogonal design method in geotechnical parameter back analysis for underground structures | |
| CN115358099B (en) | Numerical simulation method and system for full-length bonding anchor cable | |
| Iraji et al. | A modification to dense sand dynamic simulation capability of Pastor–Zienkiewicz–Chan model | |
| Yin et al. | Estimation of robustness of time integration algorithms for elasto-viscoplastic modeling of soils | |
| Chen et al. | A two-surface plasticity model for cyclic behavior of saturated clay | |
| Wang et al. | Geotechnical localization analysis based on Cosserat continuum theory and second-order cone programming optimized finite element method | |
| CN113221371B (en) | Method and device for determining critical sliding surface of side slope and terminal equipment | |
| Dai et al. | Removal of singularities in Hoek-Brown criterion and its numerical implementation and applications | |
| CN116542011A (en) | Rock water absorption softening effect particle discrete element simulation method, device and storage medium | |
| Cui et al. | A unified critical state parameter model for sand and overconsolidated clay in the framework of subloading surface theory | |
| Mousavian et al. | Torsion–Shear Behaviour at Interlocking Joints: Calibration of Discrete Element-Deformable Models Using Experimental and Numerical Analyses | |
| CN115659758A (en) | Shield tunnel rock-soil parameter inversion and tunneling parameter optimization method based on approximate model | |
| Weng et al. | An elastoplastic constitutive model for unsaturated loess with internal structure | |
| CN114329701A (en) | Large-deformation tunnel buffer layer support design method | |
| Permanoon et al. | Failure of existing structures with semi-brittle mechanical properties on meso scale and reduction of computational cost using non-linear topology optimization |
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 |