Open Access
 Issue Mechanics & Industry Volume 18, Number 6, 2017 607 5 https://doi.org/10.1051/meca/2017032 06 March 2018

## 1 Introduction

Quasi-brittle materials are characterized by an asymmetric behavior in tension and compression, and a dependence of the damage evolution to the confining pressure (influence of microcracks closure). This damage effectiveness has been sometimes modeled by using two different damage variables to separate the opening and closing states [1]. Others authors had introduced a set of variables or dependent parameters, often inspired by the mechanisms at the microstructure level such as the blocked energy by closed defects [2,3], or a 4th order tensor parameter [4], but also open microcracks ratio [5], and the opening mode cracks concept [6,7], etc. We choose here a microplane formalism introduced in 1938 (Taylor [8]) and first applied to metal plasticity (Batdorf Budianski [9]). Since, they have been applied to brittle materials (concrete, ceramics, rocks, ice, etc.). For example, we can cite the work of Bazant et al. [1014] Carol et al. [15,16] and Ramm et al. [17,18]. A detailed comparison of the main classes of those models is presented by Benelfellah et al. [19].

Microplane models consider that each direction of the solid angle chosen fixed, orients a microplane with its associated behavior. The global deformations are projected on each microplane. There are several types of formulations and we will focus on the VDT one (Volumetric, Deviatoric and Tangential). A constitutive law in the plane defines the state of stress associated from the projected deformations. The global stress is built respecting the principle of virtual work, from the stress on each microplane, through an integral over the whole solid angle. In the following simulations, the solid angle is discretized into 42 equally spaced angular sectors, ensuring the ability to describe the elastic isotropy.

## 2 VDT model with effectiveness of damage

Bazant and Gambarova [8] presented this approach to simulate the response of concrete submitted to compressive loading. This formulation is based on the decomposition of the strain tensor for a given direction n in a tangential part , and a normal one which is composed of a volumetric part, ϵV and a deviatoric one ϵD (Fig. 1). (1) (2)

The consequences of opening/closing of microcracks on the effect of damage are introduced into the free energy of each microplane using effectiveness functions denoted αV(ϵV), αD(ϵD) and αT(ϵT), and [20].

The free energy for a given microplane is thus expressed: (3) where kV, μD and μT are defined as the volumetric, deviatoric and tangential moduli, respectively. Applied to microplane, the second principle leads to the equation (4) and must be verified for any loading path. It results a first part giving the state variables laws at the microplane level and a second part giving the dissipation due to damage growth whose the positivity is usually guaranteed by the damage irreversibility. Three dissipative terms, dependent on the effectiveness evolution, appear (second line of Eq. (4)), and raise two questions. Should these terms be zero? If they are not null, what their physical senses? One approach is to assume non-dissipative effectiveness mechanisms (no hysteresis phenomena for example), and these terms should then always remain zero. A second approach is to introduce a dissipation, due to the effectiveness, whose the positivity and the physical interpretation must be given. (4)

For simplicity's sake, the function of volumetric effectiveness αV(ϵV) is a Heaviside function of the observable variable ϵV which guarantees the nullity of the term , that results in no spurious dissipation, and in a continuously differentiable free energy. We accept a partial effectiveness of 50% on the deviatoric damage when microcracks are “closed” (ϵD ≤ 0). The tangential damage is considered always effective (αT(ϵT) = 1). Thus: (5) (6)

 Fig. 1VDT decomposition [14].

## 3 Identification and damage evolution law

Tests called “H10” (respectively “H05”) consist in two main steps; a hydrostatic confinement one until 10 MPa (respectively 5 MPa) followed by a cyclic longitudinal compression one. Tests called “TS” referred to simple tensile test and “CS” to the simple compressive test [21].

It is difficult to identify microplane models because the experimental overall response results from the sum of the elementary responses on each plane, and thus is different from the “local” response of the microplane. It is therefore necessary to apply at first basic behaviors on microplane to identify afterwards the consequences once the global stress tensor is obtained.

Although a coupled damage parameters modeling would have been preferable, for pragmatic reasons, we propose a first approach with 3 independent damage variables. These damage variables are expected to evolve according to the associated thermodynamic forces Fdi, calculated from the free energy. Like the model proposed by Bazant et al., A dependency of the damage evolution at the volume change using the p parameter is added. The evolution laws are given by the following relationships: (7) (8) where is the initial damage of the material (i replaces indices V, D or T). a3, a4 and a5 are parameters to identify. is the negative volumetric strain attained after the isostatic confinement of 10 MPa. Thus, the variable p is equal to 1 at the start of the compression phase of the test H10 () and it is higher than 1 at the beginning of the compressive test CS without confinement. Figure 2 illustrates the effect of the parameters on the damage evolution (Tab. 1).

Dependence on the pressure of the response of the material and its initial microcracking led to consider this material pre-damaged.One can see the difference of the bulk and deviatoric moduli in tension and compression because of the damage effectiveness parameters αV, αD and the consolidating effect of the pressure on the deviatoric response (Fig. 3).

 Fig. 2Effect of the damage parameters on the damage variable evolution.
Table 1

Model parameters identified for an energetic material.

 Fig. 3Microplane volumetric (A) and deviatoric (B) strain–stress responses. Damage and effectiveness effects of a tensile loading on compression one.

## 4 Simulation and comparison to experimental results

To introduce residual strain, a parabolic isotropic yield criteria in the plane (effective pressure − effective octahedral stress ) is “added” to the microplane model. The determination of the yield surface and the flow law, and the identification of the strain hardening are carried out from Le [20]. (9)

The effective stress is determined using the 4th order damage operator (I − D) −1 calculated by the microplane part of the model: .

The simulations of tensile test, compressive test and compressive test with confining pressure at 5 and 10 MPa for the complete model are compared with the experimental responses in Figure 4. The simulations are interrupted just after the last experimental point. The experimental responses in tensile test and compression test are quite well reproduced. Zooming in on Figure 4 could show a slight difference in the tension and compression moduli. In a confined compressive test, a higher modulus at the origin is found and the damage evolve more slowly because of the confining pressure term in the evolution law. Figure 5 compares the simulation to the response of a confined multi cyclic test (10 MPa). The reported points along the axes for a stress of −10 MPa correspond to residual strain. In this figure the residual strains of the model are similar to the experimental ones. The balance between damage and plasticity observed during discharges is quite good.

Given the lack of experimental data for non-proportional multiaxial loading, the influence of the effectiveness and the damage anisotropy can be analyzed only by the model. The damage distribution is presented for a given value of the longitudinal strain (1%) in a plane containing the loading direction e1 (Fig. 6). Firstly, there is a strong anisotropy of the tangential damage (“□”), which is maximum in the longitudinal and transverse directions, and the deviatoric damage (“+”), which is maximum in a 45° direction relative to the loading one e1. Note that in the loading direction, the effective deviatoric damage (“O”) is half that calculated (“+”) because of the effectiveness parameter. For the same level of 1% strain, we note that the damage for compressive test with confining pressure is less developed than in compressive test without pressure, finding thus the limiting effect of the confining pressure on the damage evolution.

 Fig. 4Simulation of tensile test, compression test and compression with confining pressure (5–10 MPa) for an aggregate energetic material.
 Fig. 5Simulation of cyclic compression with confinement (10 MPa) for an aggregate energetic material.
 Fig. 6Rosette of damage for 1% of longitudinal strain in compressive test (A) and triaxial compressive test with 10 MPa of confinement (B). e1 is the loading direction.

## 5 Conclusion

In this article two phenomena of effectiveness have been introduced into the VDT model respecting the thermodynamic framework. To ensure zero dissipation, these mechanisms are governed by the associated deformations. Simple Heaviside functions have been used to manage the stiffness changes on each microplane. With regard to deviatoric damage, the loss or the recovery of the modulus is partial when changing from tension to compression. No dissipative mechanism being introduced, the stiffness discontinuities occur for a null elastic strain to cancel the strain jump that could results from the variation of stiffness. For compression loads, pressure effect has been take into account in the evolution of the damage. This allows to identify the model from the triaxial test “H10” and then translate the behavior observed for compressive test. As commonly seen in some concrete models, we have integrated two damage laws to control tension-compression asymmetry. The plastic behavior has been “added” to the damageable elastic model to predict permanent deformations experimentally observed. Anisotropic effective stress formalism has permitted a decoupling of plasticity and damage. The treatment of results has been adapted and the parameters redetermined. Comparisons between calculations and experiments show a quite good agreement even if a final optimization of the parameters is necessary to further increase thereof.

## References

1. P. Ladeveze, A. Gasser, O. Allix, Damage mechanisms modeling for ceramic composites, J. Eng. Mater. Technol. 116 (3) (1994) 331–336 [CrossRef] [Google Scholar]
2. A. Boursin, J.L. Chaboche, F. Roudolf, Mécanique de l'endommagement avec conditions unilatérales et stockage d'énergie élastique, Comptes Rendus de l'Académie des Sciences Paris, Série II (323) 369–376, 1996. [Google Scholar]
3. J.F. Maire, J.L. Chaboche, A new formulation of continuum damage mechanics (CDM) for composite materials, Aerosp. Sci. Technol. 1 (4) (1997) 247–257 [CrossRef] [Google Scholar]
4. D. Halm, A. Dragon, Y. Charles, A modular damage model for quasi-brittle solids interaction between initial and induced anisotropy, Arch. Appl. Mech. 72 (2002) 498–510 [CrossRef] [Google Scholar]
5. S. Baste, R. El Guerjouma, G. Alain, Mesure de l’endommagement anisotrope d’un composite céramique-céramique par une méthode ultrasonore, Revue de Physique Appliquée (Paris) 7 (1989) 721–731 [CrossRef] [EDP Sciences] [Google Scholar]
6. A. Thionnet, From fracture to damage mechanics: a behavior law for microcracked composites using the concept of crack opening mode, Compos. Struct. 92 (3) (2010) 780–794 [CrossRef] [Google Scholar]
7. A. Thionnet, J. Renard. Modelling unilateral damage effect in strongly anisotropic materials by the introduction of the loading mode in damage mechanics, Int. J. Solids Struct. 36 (28) (1999) 4269–4287 [CrossRef] [Google Scholar]
8. G.I. Taylor, Plastic strain in metals, J. Inst. Metals 62 (1938) 307–324 [Google Scholar]
9. S. Batdorf, B. Budiansky, A mathematical theory of plasticity based on the concept of slip, Technical Note No. 1871, National Advisory Committee for Aeronautics, 1949 [Google Scholar]
10. Z.P. Bazant, P.G. Gambarova, Crack shear in concrete: crack band microplane model, J. Struct. Eng. 110 (9) (1984) 2015–2035 [Google Scholar]
11. Z.P. Bazant, B.H. Oh, Deformation of cracked net-reinforced concrete walls, J. Struct Eng. 109 (1) (1983) 93–108 [CrossRef] [Google Scholar]
12. Z.P. Bazant, B.H. Oh, Microplane model for progressive fracture of concrete and rock, J. Eng. Mech. 111 (4) (1985) 559–582 [CrossRef] [Google Scholar]
13. A. Beghini, Z.P. Bazant, Y. Zhou, O. Gouirand, F.C. Caner, Microplane model m5f for multiaxial behavior and fracture of fiber-reinforced concrete, J. Eng. Mech. 133 (1) (2007) 66–75 [CrossRef] [Google Scholar]
14. I. Carol, M. Jiràsek, Z. Bazant, A thermodynamically consistent approach to microplane theory. part I. free energy and consistent microplane stresses, Int. J. Solids Struct. 38 (17) (2001) 2921–2931 [CrossRef] [Google Scholar]
15. I. Carol, K. Willam, Spurious energy dissipation/generation in stiffness recovery models for elastic degradation and damage, Int. J. Solids Struct. 33 (20–22) (1996) 2939–2957 [CrossRef] [Google Scholar]
16. E. Kuhl, P. Steinmann, I. Carol, A thermodynamically consistent approach to microplane theory. Part II. Dissipation and inelastic constitutive modeling, Int. J. Solids. Struct. 38 (17) (2001) 2933–2952 [CrossRef] [Google Scholar]
17. E. Kuhl, E. Ramm, R. Borst, An anisotropic gradient damage model for quasi-brittle materials, Comput. Methods in Appl. Mech. Eng. 183 (1–2) (2000) 87–103 [CrossRef] [Google Scholar]
18. M. Leukart, E. Ramm, A comparison of damage models formulated on different material scales, in: Twelfth Int. Workshop on Comput. Mech. Mat. Compt. Mat. Sci. 28, 2003, pp. 749–762 [Google Scholar]
19. A. Benelfellah, A. Frachon, M. Gratton, M. Caliez, D. Picart, Analytical and numerical comparison of discrete damage models with induced anisotropy, Eng. Fracture Mech. 121–122 (2014) 28–39 [CrossRef] [Google Scholar]
20. V.D. Le, Modélisation et identification du comportement plastique visco-élastique endommageable d' un matériau agrégataire, Thesis, Université François Rabelais de Tours, 2007 [Google Scholar]
21. D. Picart, A. Benelfellah, J.L. Brigolle, A. Frachon, M. Gratton, M. Caliez, Characterization and modeling of the anisotropic damage of a high explosive composition, Eng. Fract. Mech. 131 (2014) 525–537 [Google Scholar]

Cite this article as: A. Benelfellah, M. Gratton, M. Caliez, A. Frachon, D. Picart, VDT microplane model with anisotropic effectiveness and plasticity, Mechanics & Industry 18, 607 (2017)

## All Tables

Table 1

Model parameters identified for an energetic material.

## All Figures

 Fig. 1VDT decomposition [14]. In the text
 Fig. 2Effect of the damage parameters on the damage variable evolution. In the text
 Fig. 3Microplane volumetric (A) and deviatoric (B) strain–stress responses. Damage and effectiveness effects of a tensile loading on compression one. In the text
 Fig. 4Simulation of tensile test, compression test and compression with confining pressure (5–10 MPa) for an aggregate energetic material. In the text
 Fig. 5Simulation of cyclic compression with confinement (10 MPa) for an aggregate energetic material. In the text
 Fig. 6Rosette of damage for 1% of longitudinal strain in compressive test (A) and triaxial compressive test with 10 MPa of confinement (B). e1 is the loading direction. In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.