.. _r7.01.44-cas_tests: Verification and validation ========================== Test cases --------- The references of the test cases and associated documentation are given in :numref:`r7.01.44-table_cas_tests`. .. _r7.01.44-table_cas_tests: .. list-table:: Cas-tests et documentations associées. *-**Test case reference** - **Documentation reference** - **Description** *-*ssnv160f* - [v6.04.160] - Isotropic compression test *-*ssnv205c* - [v6.04.205] - Cyclic shear test *-*ssnv207c* - [v6.04.207] - Cyclic shear test with micro-discharge *-*comp012i* - [v6.07.112] - Compatibility test with command CALC_GEO_MECA *-*wtnv122e* - [v7.31.122] - Undrained triaxial compression test Examples of answers to the material point -------------------------------------- The :numref:`r7.01.44-compression_isotrope`, :numref:`r7.01.44-compression_triaxiale` and :numref:`r7.01.44-cisaillement_cyclique` present some answers obtained using the MTest executable (answers to the hardware point): * Isotropic compression test. * Triaxial compression tests for several confinement stresses. * Cyclic shear tests for several deformation amplitudes. The parameters used in these simulations are grouped together in :numref:`r7.01.44-table_parametres_simules`. .. _r7.01.44-table_parametres_simules: .. list-table:: Paramètres du modèle simulé. *-**Intervention** - **Appellation** - **Definition** - **Symbol** - **Value** * - Elasticity - | *BulkModulus* | *ShearModulus* | *ShearModulusRatio* - | Total compressibility module | Total shear modulus | Ratio of the shear modulus of component 1 to the total shear modulus - |:math: `K` |:math: `\ mu` |:math: `\ rho` - | 516 MPa | 238 MPa | 0.1 * - Component 1 - | *critstateSlope* | *initCritPress* | *IncoplastIndex* | *IsoHardRatio* | *isoHardIndex* - | Critical state slope | Initial critical pressure | Plastic incompressibility index | Homothetic reduction ratio of the initial elasticity domain | Hardening index by homothetic enlargement of the initial elasticity domain - |:math: `M` |:math: `p_ {c0} ` |:math: `\ beta` |:math: `\ eta` |:math: `\ omega` - |:math: `1.38` |:math: `100` kPa |:math: `30` |:math: `0.99` |:math: `32` * - Component 2 - | *HypDistortion* | *HyperExponent* | *miNcritPress* - | Reference distortion of the "modified hyperbolic" relationship:eq: `calibration_shear_3` | Curvature parameter for the "modified hyperbolic" relationship:eq: `calibration_shear_3` | Minimum pressure at which the critical state:eq: `systeme_critical_status` is achievable - |:math: `\ gamma_ {\ mathrm {hyp}}` |:math: `n_ {\ mathrm {hyp}} ` |:math: `C` - |:math: `2.10^ {-4} ` |:math: `0.78` |:math: `448` kPa The :numref:`r7.01.44-compression_isotrope` shows online the response of the isotropic compression test for the values of :numref:`r7.01.44-table_parametres_simules`. In the charging phase, the response is elastic until reaching the pressure :math:`-\sigma_m=2p_{c0}(1-\eta)=0.2` kPa (:math:`I_p=0`) and then becomes irreversible, with activation of only the plasticity of the first component (:math:`I_p=1`). The discharge is elastic. In comparison, the broken line shows the answer in case :math:`\eta=0`, with the other parameters unchanged. In this situation, the response remains elastic up to a pressure of :math:`-\sigma_m=200` kPa. Thus, using the :math:`\eta` parameter allows to introduce stronger initial compaction for low pressure levels, where volume deformation :math:`\varepsilon_v` is higher under compression. This possibility is interesting for modeling behavior strongly. contracting soils with high porosity. .. _r7.01.44-compression_isotrope: .. figure:: images/compression_isotrope.svg :align: center :width: 800 Isotropic pressure-controlled compression test :math:`-\sigma_m`. :numref:`r7.01.44-compression_triaxiale` shows triaxial compression test responses for three :math:`p_{\mathrm{conf}}` confinement levels. We observe that the higher the lockdown, the greater the pressure :math:`-\sigma_m` reached asympotically. In addition, when this pressure exceeds :math:`C=448` kPa, the volume deformation :math:`\varepsilon_v` stabilizes. On the other hand, when this pressure is lower, as can be clearly seen for :math:`p_{\mathrm{conf}}=100` kPa, the volume deformation continues to increase, indicating dilating behavior. This observation is consistent with the critical state equations referenced in :eq:`systeme_etat_critique` -4. .. _r7.01.44-compression_triaxiale: .. figure:: images/compression_triaxiale.svg :align: center :width: 800 Triaxial compression tests with three confinement stresses :math:`-\sigma_{xx}=-\sigma_{yy}=p_{\mathrm{conf}}`. The :numref:`r7.01.44-cisaillement_cyclique` presents, in a continuous line, the response of shear tests at constant pressure :math:`-\sigma_{m}=200` kPa for :numref:`r7.01.44-table_parametres_simules` values, with multiple :math:`\gamma=2\varepsilon_{xy}\in[10^{-6};10^{-2}]` distortion amplitudes. For each distortion value, one cycle is made according to chronology :math:`\gamma\rightarrow-\gamma\rightarrow\gamma`. Shear stress :math:`\sigma_{xy}` has hysteresis loops, which are not closed due to the fact that a single cycle does not make it possible to stabilize the work hardening associated with the volume components of the internal variables. The responses are compared to those predicted with :math:`\rho=0` and :math:`\rho=1`, the other parameters remaining unchanged. The predictions of model CSSM then join those of the second component alone (case :math:`\rho=0`), and deviate more significantly than those of the first component alone (case :math:`\rho=1`) on the evolution of the standardized secant shear modulus :math:`\mu_{\mathrm{secant}}/\mu` and reduced depreciation. .. _r7.01.44-cisaillement_cyclique: .. figure:: images/cisaillement_cyclique.svg :align: center :width: 1200 Cyclic shear tests at constant pressure :math:`-\sigma_{m}=200` kPa. Example of experimental comparison ------------------------------------ The :numref:`r7.01.44-comparaison_rho_0`, :numref:`r7.01.44-comparaison_rho_1` and :numref:`r7.01.44-comparaison_rho_01` compare the predictions of the model CSSM with experimental data on loose rock materials from two tests triaxial compression at :math:`p_{\mathrm{conf}}=100` kPa and :math:`p_{\mathrm{conf}}=200` kPa confinements, as well as cyclic shear tests. On triaxial compression tests, model CSSM was calibrated "at best" by adjusting the five parameters :math:`p_{c0},M,\beta,\omega` and :math:`C`, assuming the other parameters of :numref:`r7.01.44-table_parametres_simules` and considering three distinct cases that have already been the subject of previous comparisons: * :math:`\rho=0\leftrightarrow` effect of model component 2 only on shear behavior. * :math:`\rho=1\leftrightarrow` effect of component 1 of the model only. * :math:`\rho=0.1\leftrightarrow` effect of the two components of the model. The last case corresponds to considering that the shear stiffness in component 1 :math:`\rho\mu` ("static" stiffness) is ten times lower than the total shear stiffness :math:`\mu` ("dynamic" stiffness). This ratio of ten is not unreasonable, compared to the relationships observed between the dynamic and static modules for ground materials. The :numref:`r7.01.44-comparaison_rho_0` compares the predictions for the case where :math:`\rho=0` with the experimental data. It is clear that component 2 of the model alone fails to reproduce realistically the evolutions of the equivalent von Mises stress :math:`\sigma_{eq}` and volume deformation :math:`\varepsilon_v` observed during compression tests triaxial. On the other hand, component 2 offers a fairly accurate representation of the standardized secant shear modulus :math:`\mu_{\mathrm{secant}}/\mu` during shear tests. As for reduced damping, the prediction progressively overstates the experimental evolution with the increase in distortion level :math:`\gamma`. This prediction is very similar to that of the Iwan model but also of Hujeux [:ref:`r7.01.23 `]. .. _r7.01.44-comparaison_rho_0: .. figure:: images/comparaison_exp_CSSM_rho_0.svg :align: center :width: 1600 Comparison of the predictions in case :math:`\rho=0` to the experimental data. The :numref:`r7.01.44-comparaison_rho_1` compares the predictions for the case where :math:`\rho=1` with the experimental data. In this case, component 1 of the model significantly improves the evolution of the equivalent von Mises stress during compression tests. triaxial. However, only the volume deformation for confinement :math:`p_{\mathrm{conf}}=200` kPa is correctly reproduced. For lockdown :math:`p_{\mathrm{conf}}=100` kPa, the prediction lacks dilatance. With regard to shear tests, the predictions of the standardized secant shear modulus :math:`\mu_{\mathrm{secant}}/\mu` and of the reduced damping are greatly deteriorated compared to to those obtained for case :math:`\rho=0`. .. _r7.01.44-comparaison_rho_1: .. figure:: images/comparaison_exp_CSSM_rho_1.svg :align: center :width: 1600 Comparison of the predictions in case :math:`\rho=1` to the experimental data. The :numref:`r7.01.44-comparaison_rho_01` finally presents the predictions in case :math:`\rho=0.1`. This combination of components 1 and 2 of the model allows: * To improve the observed evolution of dilatance for :math:`p_{\mathrm{conf}}=100` kPa confinement during triaxial compression tests, while maintaining a correct representation of the equivalent von Mises stress. This representation is similar to that obtained with component 1 alone in the case where :math:`\rho=0`. * To obtain predictions close to those of component 2 for shear tests when :math:`\rho=1`. In summary, combining the two components of the model makes it possible to take advantage of the best characteristics of each individual component, thus offering an improved representation for both types of monotonic and cyclic tests. .. _r7.01.44-comparaison_rho_01: .. figure:: images/comparaison_exp_CSSM.svg :align: center :width: 1600 Comparison of the predictions in case :math:`\rho=0.1` to the experimental data. .. How to get the gray box here-under .. **Note:** The previous comparisons validate the relevance of model CSSM on shear loads, with calibration parameters :math:`p_{c0},M,\beta,\omega` and :math:`C` on triaxial compression tests. However, in order to assess the capabilities of the model more comprehensively and rigorously, it will be necessary to compare it to isotropic compression tests. This step will make it possible to estimate more accurately the parameter :math:`\eta`, here taken close to one in :numref:`r7.01.44-table_parametres_simules`, and refine the value of :math:`\omega`, which significantly influence the initial compaction, as we saw on :numref:`r7.01.44-compression_isotrope` with :math:`\eta\in\{0,0.99\}`.