CN111062162B - Numerical modeling and application method of rock and soil material accurate constitutive model - Google Patents
Numerical modeling and application method of rock and soil material accurate constitutive model Download PDFInfo
- Publication number
- CN111062162B CN111062162B CN201911273901.7A CN201911273901A CN111062162B CN 111062162 B CN111062162 B CN 111062162B CN 201911273901 A CN201911273901 A CN 201911273901A CN 111062162 B CN111062162 B CN 111062162B
- Authority
- CN
- China
- Prior art keywords
- stress
- stress path
- geotechnical
- compared
- path
- 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.)
- Active
Links
- 239000000463 material Substances 0.000 title claims abstract description 99
- 238000000034 method Methods 0.000 title claims abstract description 66
- 239000002689 soil Substances 0.000 title claims abstract description 23
- 239000011435 rock Substances 0.000 title description 9
- 238000012360 testing method Methods 0.000 claims abstract description 34
- 230000014509 gene expression Effects 0.000 claims abstract description 15
- 230000004069 differentiation Effects 0.000 claims abstract description 6
- 230000035945 sensitivity Effects 0.000 claims abstract description 6
- 238000004458 analytical method Methods 0.000 claims description 15
- 238000010276 construction Methods 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 19
- 239000004927 clay Substances 0.000 description 12
- 230000006835 compression Effects 0.000 description 6
- 238000007906 compression Methods 0.000 description 6
- 230000003993 interaction Effects 0.000 description 4
- 230000006399 behavior Effects 0.000 description 3
- 238000005094 computer simulation Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 description 1
- 238000009864 tensile test Methods 0.000 description 1
Images
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention discloses a numerical modeling method of a geotechnical material accurate constitutive model, which comprises the following steps: a, performing a triaxial test on a rock-soil sample under a given stress path to obtain a series of plastic body strain values and plastic shear strain values of the rock-soil sample on a p-q stress plane; respectively carrying out interpolation fitting on the strain value and the plastic shear strain value of the plastic body along the directions p and q by adopting a spline function to obtain a polynomial expression f for describing the strain field and the shear strain field of the plastic body 1 (p, q) and f 2 (p, q); to f 1 (p, q) and f 2 (p, q) respectively carrying out partial differentiation on p, q to obtain 4 partial derivatives; and establishing an equation expression of the constitutive model for determining the geotechnical material under the specified given stress path based on the four partial derivatives. The constitutive equation of the rock-soil material can quantitatively describe the basic characteristics of plastic deformation of the rock-soil material: the correlation of pressure sensitivity, shear expansion and stress path can describe the deformation process of the rock-soil material as truly as possible.
Description
Technical Field
The invention belongs to the technical field of computer simulation and modeling of geotechnical engineering, and particularly relates to a method and a device for determining a constitutive model of a geotechnical material, computer equipment and a storage medium.
Background
Prior to the last 60 s of the century, only relatively simple constitutive equations, such as linear elastic models, were used in stress analysis of geotechnical structures due to the lack of high-speed and large-memory-capacity computers. In the early stage of the last 60 th century, as large-scale high-speed computers begin to develop rapidly and new numerical techniques such as a finite element method, a boundary element method and the like appear, conditions are created for using complex and accurate constitutive equations, and therefore the development of an elastic-plastic constitutive model of the rock-soil material is promoted. Currently, super electronic computers have appeared, providing a wider development space for geotechnical engineering computer simulation and emulation.
Unlike metals, body strain occurs in plastic deformation of rock and soil, so that plastic deformation behaviors of the rock and soil are quite complex, at present, the true mechanism of the plastic deformation behaviors of rock and soil materials is not satisfactorily explained, and the constitutive model of the existing rock and soil materials cannot reflect some important characteristics of the plastic deformation of the rock and soil materials, such as stress path correlation.
The constitutive model of the geotechnical material generally comprises two parts, namely a mathematical expression of a constitutive equation and a test determination method of parameters in the constitutive equation. The principle of interaction between mass distribution and deformation of geotechnical materials is proposed in 2010 by the WangJingtao of the applicant, and the principle indicates that the mass distribution-deformation interaction is a root cause of basic characteristics of plastic deformation behavior of the geotechnical materials, and is also represented by interaction between plastic strain and shear strain. Based on the principle, the isoline of the plastic body strain and the generalized shear strain in the stress space is adopted
Respectively, the trajectories of the volume and shear yield, wherein
And
plastic volume strain and generalized shear strain, respectively, p and q mean normal stress and generalized shear stress, respectively, f 1 (p, q) and f 2 (p, q) are two nonlinear functions in a stress space, the two groups of yield tracks respectively describe strain and shear strain fields of the plastic body, and the constitutive equation of the geotechnical material is established by utilizing the symmetry of the two strain fields:
or
Wherein,
and
are respectively f 1 And f 2 Is the stress increment vector (d) along a given stress path p ，d q ) The constitutive equation quantitatively reflects the interaction of mass distribution and deformation, so that the basic characteristics of the plastic deformation of the geotechnical material can be described: pressure sensitivity, shear expansion, and stress path dependence, where a stress path is the process by which the stress state of a material element changes during structural deformation, and stress path dependence is the effect of a stress path on the plastic deformation of a material. Experiments have shown that the stress path dependence is a comprehensive manifestation of pressure sensitivity and shear expansion. However, the method for determining the corresponding parameters, namely the four coefficients of the constitutive equation, is also lacked.
Disclosure of Invention
The invention aims to overcome the defects and provides a numerical modeling and application method of a geotechnical material accurate constitutive model so as to establish the accurate constitutive model capable of quantitatively describing the basic characteristics of the geotechnical material plastic deformation.
In order to solve the problems, the specific technical scheme of the invention is as follows:
a numerical modeling method of a geotechnical material accurate constitutive model comprises the following steps:
a, performing a triaxial test on a rock-soil sample under a given stress path to obtain a series of plastic body strain values and plastic shear strain values of the rock-soil sample on a p-q stress plane;
b, respectively carrying out interpolation fitting on the strain value and the plastic shear strain value of the plastic body along the p and q directions by adopting a spline function to obtain a polynomial expression f for describing the strain and shear strain field of the plastic body 1 (p, q) and f 2 (p，q)；
C for f 1 (p, q) and f 2 (p, q) respectively carrying out partial differentiation on p, q to obtain 4 partial derivatives;
d, establishing a constitutive model equation expression for determining the geotechnical material under the specified given stress path based on the four partial derivatives:
the constitutive equation of the geotechnical material can quantitatively describe the basic characteristics of plastic deformation of the geotechnical material: pressure sensitivity, shear swell and stress path dependence.
Preferably, in the step a, the given stress path is a stress path occurring during plastic deformation of the geotechnical engineering structure.
Preferably, in step B, the spline function is a B-spline curve.
The invention also discloses an application method of the numerical modeling of the rock-soil material accurate constitutive model, which comprises the following steps:
s1, performing stress analysis on a geotechnical engineering structure to be simulated by adopting the elastic model of the geotechnical material based on a finite element program to obtain a primary stress path of a geotechnical material unit;
s2, constructing a to-be-optimized constitutive model of the geotechnical material by applying the numerical modeling method of the accurate constitutive model of the geotechnical material, wherein the given stress path is the preliminary stress path, and performing stress analysis on the geotechnical engineering structure based on a finite element program to obtain to-be-compared stress paths of the material units;
s3, judging whether the types of the preliminary stress path and the stress path to be compared are consistent or not:
s31, if the types of the preliminary stress path and the to-be-compared stress path are judged to be all consistent, taking the to-be-optimized constitutive model constructed in the current round as the accurate constitutive model of the geotechnical material;
s32, if the type of the preliminary stress path is not consistent with that of the stress path to be compared:
s321, constructing a to-be-optimized constitutive model of the geotechnical material by applying the numerical modeling method of the accurate constitutive model of the geotechnical material, wherein the given stress path is the to-be-compared stress path, and performing stress analysis on the to-be-optimized constitutive model based on the finite element method to obtain the to-be-compared stress path of the material unit constitutive wheel;
s322, judging whether the type of the stress path to be compared in the previous round is consistent with that of the stress path to be compared in the current round;
s3221, if the type of the stress path to be compared in the previous round is judged to be consistent with that of the stress path to be compared in the current round, taking the constitutive model to be optimized constructed in the current round as the accurate constitutive model of the geotechnical material;
s3222, if it is determined that the type of the stress path to be compared in the previous round is not consistent with the type of the stress path to be compared in the current round, returning to step S321 until it is determined that the type of the stress path to be compared in the previous round is consistent with the type of the stress path to be compared in the current round.
The numerical modeling method of the accurate constitutive model of the geotechnical material disclosed by the invention is characterized in that a series of triaxial tests are carried out on various types of stress paths under the determined stress paths (which are stress paths possibly generated in all parts of the structure in the deformation process of the actual geotechnical engineering structure), and based on the results of the triaxial tests, 4 coefficients of corresponding constitutive equations are determined by adopting a numerical technology. Meanwhile, based on the constitutive equations determined by the method, the constitutive equations are determined under different types of stress paths, and each set of established constitutive equations is only suitable for a specific type of stress path, so that multiple sets of constitutive equations suitable for different stress paths in the deformation process of the actual geotechnical engineering structure can be obtained, and the deformation process can be described as truly as possible.
The invention also provides an application method of the numerical modeling of the accurate constitutive model of the geotechnical material, which comprises the steps of firstly determining a primary stress path (particularly a stress path type) of a geotechnical engineering structure based on an elastic model, then applying the numerical modeling method of the accurate constitutive model of the geotechnical material as described in the embodiment I, carrying out a three-axis test under the specified primary stress path, constructing the to-be-optimized constitutive model of the geotechnical material of the wheel, comparing the types of the primary stress path and the stress path to be compared, if the type of the primary stress path is the same as that of the stress path to be compared, selecting the to-be-optimized constitutive model of the current time as the accurate constitutive model of the geotechnical material, if the types of the primary stress path and the stress path to be compared are different, or a new stress path type appears, and the like, applying the numerical modeling method of the accurate constitutive model of the geotechnical material as described in the embodiment I, carrying out the three-axis test under the specified stress path to-be-compared, comparing the stress path and comparing the stress path with the current stress path until the current stress path is consistent.
Firstly, in the method for applying the numerical modeling of the accurate constitutive model of the geotechnical material, the four coefficients of the constitutive equation are not only functions of stress points but also related to stress path changes because the numerical modeling method is applied to the constitutive model to be optimized, and the constitutive model to be optimized, which is established by the method, has the capability of reflecting the stress path correlation. Meanwhile, constitutive equations are constructed based on the numerical modeling method, the constitutive equations are determined under different types of stress paths (the types of the stress paths are determined based on the actual geotechnical engineering structure), each set of constitutive equations is only suitable for a specific type of stress paths, and therefore multiple sets of constitutive equations suitable for different stress paths generated in the deformation process of the actual geotechnical engineering structure can be obtained. Therefore, in the stress analysis of the geotechnical engineering structure, for different points in the geotechnical engineering structure, the corresponding constitutive equation is selected according to the actually experienced stress path, and on the basis, the finally constructed constitutive equation can describe the deformation process of the geotechnical material as truly as possible by further continuously optimizing the constitutive equation.
Drawings
FIG. 1 is a flow chart of an embodiment of a numerical modeling method of a geotechnical material accurate constitutive model of the invention;
fig. 2 (a) is a volumetric yield trace of clay based on a triaxial test taken under CTC stress path;
FIG. 2 (b) is a volumetric yield trace of clay based on a triaxial test taken under a TC stress path;
FIG. 2 (c) is a plot of the volumetric yield trace of clay obtained based on a triaxial test under the RTC stress path;
fig. 3 (a) is a shear yield trace of clay obtained based on a triaxial test under CTC stress path;
FIG. 3 (b) is a shear yield trace of clay obtained based on a triaxial test under a TC stress path;
FIG. 3 (c) is a shear yield trace of clay obtained based on a triaxial test under the RTC stress path;
FIG. 4 is a flowchart of an embodiment of a numerical modeling application method of a geotechnical material accurate constitutive model of the present invention.
Detailed Description
Example one
In the deformation process of the geotechnical engineering structure, material units of all parts of the structure can experience different stress paths, mainly comprising Conventional Triaxial Compression (CTC), triaxial Compression (TC), simplified triaxial compression (RTC), corresponding tensile tests and the like, and a constitutive equation which cannot reflect the influence of the stress paths generates a large error in stress analysis, as shown in figure 1, the invention provides a numerical modeling method of a geotechnical material accurate constitutive model, which can obtain a plurality of groups of constitutive equations suitable for different stress paths in the deformation process of the actual geotechnical engineering structure, so that the deformation process can be described as truly as possible, and the method comprises the following steps:
a, performing a triaxial test on a rock-soil sample under a given stress path to obtain a series of plastic body strain values and plastic shear strain values of the rock-soil sample on a p-q stress plane;
b, respectively carrying out interpolation fitting on the strain value and the plastic shear strain value of the plastic body along the p and q directions by adopting a spline function to obtain a polynomial expression f for describing the strain and shear strain field of the plastic body 1 (p, q) and f 2 (p，q)；
C for the f 1 (p, q) and f 2 (p, q) respectively carrying out partial differentiation on p, q to obtain 4 partial derivatives;
d, establishing a constitutive model equation expression for determining the geotechnical material under the specified given stress path based on the four partial derivatives:
the constitutive equation of the geotechnical material can quantitatively describe the basic characteristics of plastic deformation of the geotechnical material: pressure sensitivity, shear swell and stress path dependence.
In the deformation process of a complex geotechnical structure, stress paths experienced by material units of various parts of the structure are different, mainly comprising Conventional Triaxial Compression (CTC), triaxial Compression (TC), simplified triaxial compression (RTC), and the like, wherein a volume yield trajectory of the clay obtained based on a triaxial test under a CTC stress path is shown in fig. 2 (a), a volume yield trajectory of the clay obtained based on a triaxial test under a TC stress path is shown in fig. 2 (b), a shear yield trajectory of the clay obtained based on a triaxial test under an RTC stress path is shown in fig. 3 (c), a shear yield trajectory of the clay obtained based on a triaxial test under a CTC stress path is shown in fig. 3 (a), a shear yield trajectory of the clay obtained based on a triaxial test under a TC stress path is shown in fig. 3 (b), a shear trajectory of the clay obtained based on a triaxial test under an RTC stress path is shown in fig. 3 (c), volume and shear trajectories are both curves observed from the yield trajectories, and for different stress paths, their shapes, orientations, and the like are all of the curvature of the paths sufficiently affect the yield trajectory, thereby completely determining yield stress path.
To truly describe this deformation process, constitutive equations must be employed that reflect the effects of different stress paths. The numerical modeling method of the accurate constitutive model of the geotechnical material disclosed by the invention is characterized in that a series of triaxial tests are carried out on each type of stress path under a given stress path, based on the results of the triaxial tests, 4 coefficients of each constitutive equation are determined by adopting a numerical technique, and different from the existing geotechnical material constitutive equations, the 4 coefficients of each constitutive equation are not determined by some parameters but are nonlinear functions in a stress space, namely, four coefficients are not only functions of stress points but also related to stress path changes, so that the established constitutive equation has the capability of reflecting stress path correlation. Meanwhile, based on the constitutive equations determined by the method, the constitutive equations are determined under different types of stress paths, and each set of established constitutive equations is only suitable for a specific type of stress path, so that multiple sets of constitutive equations suitable for different stress paths in the deformation process of the actual geotechnical engineering structure can be obtained, and the deformation process can be described as truly as possible.
In the step A, under a given stress path, a series of triaxial tests are carried out aiming at each type of stress path, and plastic body strain values and plastic shear strain values of rock-soil test samples on a p-q stress plane under different stress paths are obtained; as a preferred solution, in the present embodiment, the given path is a stress path occurring during deformation of the actual geotechnical structure. Since the established constitutive equation coefficient functions are related to the stress paths, determining them must be based on the results of triaxial experiments at a given stress path.
In the step B, interpolation fitting is respectively carried out on the strain value of the plastic body strain value and the shear strain value of the plastic shear strain value along the directions p and q by adopting a spline function, and a polynomial expression f for describing the strain and shear strain fields of the plastic body is obtained 1 (p, q) and f 2 (p, q), in this step, since the coefficients of the constitutive equation are all nonlinear functions in the stress space as described above, they are determined only by numerical methods.
As a preferable solution, in this embodiment, in the step B, the spline function is a B-spline curve, that is, interpolation fitting is performed by using the B-spline function, and polynomial expressions of the spline function and the B-spline curve are obtained. The parametric curve represented by B-spline is called parametric B-spline curve, or B-spline curve, for example, an n +1 th (n) th order B-spline curve may be represented as:
meanwhile, as a preferable scheme, in the present embodiment, in the step B, interpolation fitting is performed along p and q directions by using the same pitch to obtain the polynomial expression describing the strain and shear strain field of the plastic body, in this step, for convenience of subsequent calculation, the same pitch is selected when performing interpolation fitting calculation, and further, when performing interpolation fitting along p and q directions, the same pitch may be 50kPa or 100kPa.
Finally, for the above two polynomial expressions f 1 (p, q) and f 2 (p, q) partial differentiation of p, q is performed separately, and 4 partial derivatives are obtained as 4 coefficient functions of the constitutive equation in the stress space, i.e. four coefficients in the following equation:
example two
As shown in fig. 4, the invention also provides an application method of the numerical modeling of the geotechnical material accurate constitutive model, which comprises the following steps:
s1, performing stress analysis on a geotechnical engineering structure to be simulated by adopting the elastic model of the geotechnical material based on a finite element program to obtain a primary stress path of a geotechnical material unit;
s2, constructing a to-be-optimized constitutive model of the geotechnical material by applying the numerical modeling method of the accurate constitutive model of the geotechnical material, wherein the given stress path is the initial stress path, and performing stress analysis on the geotechnical engineering structure based on a finite element program to obtain to-be-compared stress paths of the material units;
s3, judging whether the types of the preliminary stress path and the stress path to be compared are consistent or not:
s31, if the types of the preliminary stress path and the stress path to be compared are judged to be all consistent, taking the constitutive model to be optimized constructed in the current round as the accurate constitutive model of the geotechnical material;
s32, if the type of the preliminary stress path is not consistent with that of the stress path to be compared:
s321, constructing a to-be-optimized constitutive model of the geotechnical material by applying the numerical modeling method of the accurate constitutive model of the geotechnical material as described in the first embodiment, wherein the given stress path is the to-be-compared stress path, and performing stress analysis on the geotechnical engineering structure by using the to-be-optimized constitutive model by using the finite element method to obtain the to-be-compared stress path of the material unit constitutive wheel,
s322, judging whether the type of the stress path to be compared in the previous round is consistent with that of the stress path to be compared in the current round;
s3221, if the type of the stress path to be compared in the previous round is judged to be consistent with that of the stress path to be compared in the current round, taking the constitutive model to be optimized constructed in the current round as the accurate constitutive model of the geotechnical material;
s3222, if it is determined that the type of the stress path to be compared in the previous round is not consistent with the type of the stress path to be compared in the current round, returning to step S321 until it is determined that the type of the stress path to be compared in the previous round is consistent with the type of the stress path to be compared in the current round.
The invention provides an application method of numerical modeling of a geotechnical material accurate constitutive model, which comprises the steps of firstly determining a preliminary stress path of a geotechnical engineering structure based on an elastic model (, then applying the numerical modeling method of the geotechnical material accurate constitutive model as described in the embodiment I, carrying out a triaxial test under the specified preliminary stress path, constructing a to-be-optimized constitutive model of a geotechnical material of the wheel, comparing the preliminary stress path with the type of the to-be-compared stress path, if the preliminary stress path is the same as the type of the to-be-compared stress path, selecting the to-be-optimized constitutive model as the accurate constitutive model of the geotechnical material, if the preliminary stress path is different from the type of the to-be-compared stress path, or if a new stress path type appears, applying the numerical modeling method of the to-be-optimized constitutive model of the geotechnical material as described in the embodiment I, carrying out the triaxial test under the specified to-be-compared stress path, and comparing the to-be-compared stress path with the accurate constitutive model of the geotechnical material if the type is the to-be-compared, selecting the to-optimized constitutive model as the accurate constitutive model of the geotechnical material if the type is the same, comparing stress path, and comparing the cyclic stress path until the to-stress path is continued.
Firstly, the constitutive equation is constructed by applying the numerical modeling method shown in the first embodiment, and four coefficients of the constitutive equation are not only functions of stress points, but also are related to stress path change, so that the constructed constitutive equation has the capability of reflecting stress path correlation. Meanwhile, the constitutive equations constructed based on the numerical modeling method shown in the first embodiment are determined under different types of stress paths (the types of the stress paths are determined based on the actual geotechnical engineering structure), each set of constitutive equation established is only suitable for a specific type of stress path, and therefore multiple sets of constitutive equations suitable for different stress paths in the deformation process of the actual geotechnical engineering structure can be obtained, and corresponding constitutive equations are selected for different points in the geotechnical engineering structure according to the stress paths actually experienced by the points. On the basis, further through continuous iterative optimization of the constitutive equation, the finally constructed constitutive equation can describe the deformation process of the geotechnical material as truly as possible.
In the deformation process of an actual geotechnical engineering structure, stress paths experienced by material units in all parts of the structure are different. The stress path is a control factor for the development of plastic deformation of geotechnical materials, and therefore, the correlation of the stress path must be considered, and for this reason, it is first necessary to determine the type of stress path that may occur during the deformation of the actual geotechnical engineering structure. In the step S1, stress analysis is carried out on the geotechnical engineering structure to be simulated by adopting the elastic model of the geotechnical material based on a finite element program, and a primary stress path of the geotechnical material unit, namely a stress path type possibly generated in the deformation process of the actual geotechnical engineering structure, is obtained.
After the preliminary stress path is determined, performing a series of triaxial tests according to a modeling method shown in the first embodiment under the determined stress path to determine a constitutive equation, in the first embodiment, assuming that the preliminary stress path includes stress paths of three types, namely CTC, TC and RTC, performing a triaxial test based on the CTC stress path, a triaxial test based on the TC stress path, and a triaxial test based on the RTC stress path, for the rock and soil sample, respectively, and then acquiring a series of plastic body strain values and plastic shear strain values of the rock and soil sample on a p-q stress plane based on respective triaxial test results; and respectively carrying out interpolation fitting on the strain value and the plastic shear strain value of the plastic body along the directions p and q by adopting a spline function to obtain a polynomial expression f for describing the strain field and the shear strain field of the plastic body 1 (p, q) and f 2 After (p, q), the f is treated again 1 (p, q) and f 2 (p, q) respectively carrying out partial differentiation on p, q to obtain 4 partial derivatives; and establishing constitutive model equation expressions of the geotechnical materials in the three groups of stress paths, namely establishing three groups of constitutive equations of the geotechnical materials respectively. After determining the three sets of constitutive equations, stress analysis is performed on the geotechnical engineering structure based on a finite element program to obtain stress paths to be compared of the material units.
In step S3, it is determined whether the types of the preliminary stress path and the stress path to be compared are consistent: in step S2, if the constitutive equation to be optimized and the stress path to be compared determined by the finite element program are also CTC, TC, and RTC, the correctness of the constitutive equation is verified, the constitutive model to be optimized constructed in the current round is used as the accurate constitutive model of the geotechnical material, and in this embodiment, ABAQUS software is also selected by the finite element program.
If the two are different, if the stress paths to be compared determined based on the constitutive equation to be optimized and the finite element program comprise three stress paths including CTC, TC and RTC and other types of stress paths such as CRPS (constant principal stress ratio), the missing stress paths exist in the constitutive model to be optimized constructed in the current round, S321 applies the numerical modeling method of the geotechnical material accurate constitutive model described in the first embodiment, respective triaxial tests are respectively carried out under the newly determined stress paths, based on the results of the triaxial tests, the numerical technology is sampled to determine 4 coefficients of the constitutive equation under each stress path, then the geotechnical engineering structure is subjected to stress analysis based on the finite element program to obtain the stress paths to be compared of the current round of the material units; s322, continuously comparing whether the types of the stress path to be compared obtained in the first round are consistent with the types of the stress path to be compared obtained in the current round, if the types of the stress path to be compared obtained in the first round are consistent with the types of the stress path to be compared obtained in the current round, S3221 shows that the constitutive model to be optimized constructed in the current round has no missing stress path types, and taking the constitutive model to be optimized constructed in the second round as the accurate constitutive model of the geotechnical material; s3222, if the two types are not the same, returning to step S321 until it is determined that the types of the stress paths to be compared are the same as the type of the next stress path to be compared, and ending the loop.
The above description is only a preferred embodiment of the present invention, and not intended to limit the scope of the present invention, and all modifications of equivalent structures and equivalent processes, which are made by using the contents of the present specification and the accompanying drawings, or directly or indirectly applied to other related technical fields, are included in the scope of the present invention.
Claims (3)
1. An application method of numerical modeling of a geotechnical material accurate constitutive model is characterized by comprising the following steps:
s1, performing stress analysis on a to-be-simulated geotechnical engineering structure by adopting an elastic model of the geotechnical material based on a finite element program to obtain a primary stress path of a geotechnical material unit;
s2, constructing a constitutive model to be optimized of the geotechnical material, and carrying out stress analysis on the geotechnical engineering structure based on a finite element program to obtain a stress path to be compared of the material units;
the construction of the constitutive model to be optimized of the geotechnical material comprises the following steps:
a, performing a series of triaxial tests on a rock-soil sample according to stress paths of various types under a given stress path to obtain plastic body strain values and plastic shear strain values of the rock-soil sample on a p-q stress plane under different stress paths, wherein the given stress path is the primary stress path;
b, respectively carrying out interpolation fitting on the strain value and the plastic shear strain value of the plastic body along the p and q directions by adopting a spline function to obtain a polynomial expression f for describing the strain and shear strain field of the plastic body 1 (p, q) and f 2 (p，q)；
C for f 1 (p, q) and f 2 (p, q) respectively carrying out partial differentiation on p, q to obtain 4 partial derivatives;
d, establishing a constitutive model equation expression for determining the geotechnical material under the given stress path based on the 4 partial derivatives:
the constitutive equation of the geotechnical material can quantitatively describe the basic characteristics of plastic deformation of the geotechnical material: the dependence of pressure sensitivity, shear swell and stress path;
s3, judging whether the types of the preliminary stress path and the stress path to be compared are consistent or not:
s31, if the types of the preliminary stress path and the to-be-compared stress path are judged to be all consistent, taking the to-be-optimized constitutive model constructed in the current round as the accurate constitutive model of the geotechnical material;
s32, if the type of the preliminary stress path is not consistent with that of the stress path to be compared:
s321, constructing a to-be-optimized constitutive model of the geotechnical material based on the to-be-compared stress path, and performing stress analysis on the to-be-optimized constitutive model based on a finite element method to obtain the to-be-compared stress path of the material unit in the corresponding wheel;
s322, judging whether the type of the stress path to be compared in the previous round is consistent with that of the stress path to be compared in the current round;
s3221, if the type of the stress path to be compared in the previous round is judged to be consistent with that of the stress path to be compared in the current round, taking the constitutive model to be optimized constructed in the current round as the accurate constitutive model of the geotechnical material;
s3222, if it is determined that the type of the stress path to be compared in the previous round is not consistent with the type of the stress path to be compared in the current round, returning to step S321 until it is determined that the type of the stress path to be compared in the previous round is consistent with the type of the stress path to be compared in the current round.
2. The method for applying numerical modeling of an accurate constitutive model of geotechnical materials according to claim 1, wherein in the step A, the given stress path is a stress path occurring during plastic deformation of geotechnical engineering structure.
3. The method for applying numerical modeling of an accurate constitutive model of geotechnical materials according to claim 1, wherein in the step B, the spline function is a B-spline curve.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911273901.7A CN111062162B (en) | 2019-12-12 | 2019-12-12 | Numerical modeling and application method of rock and soil material accurate constitutive model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911273901.7A CN111062162B (en) | 2019-12-12 | 2019-12-12 | Numerical modeling and application method of rock and soil material accurate constitutive model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111062162A CN111062162A (en) | 2020-04-24 |
| CN111062162B true CN111062162B (en) | 2022-11-08 |
Family
ID=70298948
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911273901.7A Active CN111062162B (en) | 2019-12-12 | 2019-12-12 | Numerical modeling and application method of rock and soil material accurate constitutive model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111062162B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111783332B (en) * | 2020-06-12 | 2022-10-21 | 广西大学 | Finite element simulation method for vibration accumulation deformation of compression non-linear change and shear sudden change characteristic material |
| CN113155612B (en) * | 2021-04-16 | 2022-09-30 | 浙江科技学院 | Deformation prediction method for microfiber mixed silica sol solidified calcareous sand |
| CN113673050B (en) * | 2021-08-11 | 2022-04-12 | 中国科学院力学研究所 | Excavated object identification method based on strain response of movable arm of excavating head |
| CN115597664B (en) * | 2022-10-18 | 2024-02-06 | 中煤科工开采研究院有限公司 | Method and system for constructing digital fully-mechanized coal mining face of coal mine |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7356449B2 (en) * | 2002-04-16 | 2008-04-08 | The Board Of Trustees Of The University Of Illinois | Method and program product for solid mechanics modelling workbench and dynamic display |
| US20130289953A1 (en) * | 2012-01-24 | 2013-10-31 | The University Of Akron | Self-optimizing, inverse analysis method for parameter identification of nonlinear material constitutive models |
| CN106568645B (en) * | 2016-10-12 | 2019-03-22 | 河海大学 | A kind of rock-soil material complexity unloading stress path test method |
| CN108375505B (en) * | 2018-02-08 | 2020-04-14 | 内蒙古大学 | High-precision linear stress path test method for frozen soil |
| CN109064560B (en) * | 2018-06-27 | 2023-05-12 | 青岛理工大学 | Geological condition modeling method in subway tunnel construction process |
| CN109284523A (en) * | 2018-07-19 | 2019-01-29 | 同济大学 | A kind of rock soil medium Progressive failure, class solid-liquid phase change behavior analogy method |
| CN109446731B (en) * | 2018-12-11 | 2024-06-21 | 合肥哈工热气球数字科技有限公司 | Geotechnical engineering numerical simulation method based on ABAQUS |
| CN110008599B (en) * | 2019-04-09 | 2023-06-06 | 江西理工大学 | Water-soil coupling landslide simulation method based on high-order double-sleeve double-phase object particle method |
| CN110057675A (en) * | 2019-05-13 | 2019-07-26 | 湖南大学 | A kind of soil body disturbed after e-logp ' curve determination method |
-
2019
- 2019-12-12 CN CN201911273901.7A patent/CN111062162B/en active Active
Non-Patent Citations (1)
| Title |
|---|
| 建立岩土本构模型的数值方法;王靖涛;《华中科技大学学报》;20020331;第19卷(第1期);第44-47页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111062162A (en) | 2020-04-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111062162B (en) | Numerical modeling and application method of rock and soil material accurate constitutive model | |
| WO2021169598A1 (en) | Digital twinning-based physical simulation method for production line and system thereof | |
| CN114492211B (en) | Residual oil distribution prediction method based on autoregressive network model | |
| CN110631792B (en) | Seismic hybrid test model updating method based on convolutional neural network | |
| CN111666699B (en) | Rock mass engineering trans-scale simulation calculation method based on REV full-area coverage | |
| CN115630558B (en) | Method for predicting assembly deformation of composite material component | |
| CN104298857B (en) | Mechanism reliability calculating method based on multi-factor coupling | |
| CN109446731B (en) | Geotechnical engineering numerical simulation method based on ABAQUS | |
| CN111339594B (en) | DIC technology-based near-field dynamics parameter experiment inversion system and use method | |
| CN104809273A (en) | Creep deformation describing method | |
| CN110926973A (en) | Concrete damage accumulation analysis method under cyclic load | |
| CN104615840A (en) | Digital simulation model correcting method and system | |
| CN114638048A (en) | Three-dimensional spray pipe flow field rapid prediction and sensitivity parameter analysis method and device | |
| Shen et al. | Enhancing deep neural networks for multivariate uncertainty analysis of cracked structures by POD-RBF | |
| CN106980735B (en) | Numerical simulation method for thermal fracture of brittle material | |
| CN116227287A (en) | Fracture fluid flow fluid-solid coupling simulation method based on linear complementation method | |
| CN105844025A (en) | Non-probability aerothermoelasticity reliability designing method for hypersonic control surface | |
| CN110110406B (en) | Slope stability prediction method for achieving LS-SVM model based on Excel computing platform | |
| CN109388833B (en) | Elastic element structure optimization design method based on fatigue life | |
| CN104834795B (en) | Band connection structure nonlinear contact with friction simulated behavior method and system | |
| CN113611377A (en) | Method for simulating hybrid control creep fatigue deformation by using crystal plastic model | |
| CN111474905B (en) | Parameter drift fault diagnosis method in manufacturing process of electromechanical product | |
| CN107145665B (en) | Roadway surrounding rock stress modeling and prediction method | |
| CN118036416B (en) | Contact fatigue damage analysis method, device, equipment and storage medium | |
| CN111506998B (en) | Method for constructing parameter drift fault feature sample library in manufacturing process of electromechanical product |
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 |