Issue 
Mechanics & Industry
Volume 18, Number 6, 2017



Article Number  607  
Number of page(s)  5  
DOI  https://doi.org/10.1051/meca/2017032  
Published online  06 March 2018 
Regular Article
VDT microplane model with anisotropic effectiveness and plasticity
^{1}
Institut Pprime (UPR CNRS 3346) − ENSMA, Université de Poitiers,
CNRS 1, avenue Clément Ader – BP
40109–86961
Futuroscope Cedex, France
^{2}
Laboratoire de Mécanique et Rhéologie, E.A. 2640, B.P. 3410, INSACVL, Université Franҫois Rabelais de Tours,
41034
Blois Cedex, France
^{3}
Laboratoire de PhysicoChimie des Matériaux Appliqués, ENSAM CASABLANCA, HASSAN II University of Casablanca,
Casablanca, Morocco
^{4}
CEA, DAM, Le Ripault,
37260
Monts, France
^{*} email: a.benelfellah@gmail.com
Received:
1
December
2015
Accepted:
22
September
2017
The openingclosing state of the microcracks is a kinematic phenomenon usually modeled using a set of damage effectiveness variables, which results in different elastic responses for the same damage level. In this work, the microplane model with volumetric, deviatoric and tangential decomposition denoted VDT is modified. The influence of the confining pressure is taken into account in the damage variables evolution laws. For a better understanding of the mechanisms introduced into the model, the damage rosettes are presented for a strain given level. The model is confirmed through comparisons of the simulations with the experimental results of monotonic, and cyclic tensile and compressive testing with different levels of confining pressure.
Key words: Induced anisotropy / damage effect / microcracks / constitutive modeling / microplane / plasticity
© AFM, EDP Sciences 2017
1 Introduction
Quasibrittle 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. [10–14] 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 k_{V}, μ_{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 nondissipative 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)
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 F^{di}, 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). a_{3}, a_{4} and a_{5} 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 predamaged.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. 2
Effect of the damage parameters on the damage variable evolution. 
Model parameters identified for an energetic material.
Fig. 3
Microplane 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 nonproportional 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 e_{1} (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 e_{1}. 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. 4
Simulation of tensile test, compression test and compression with confining pressure (5–10 MPa) for an aggregate energetic material. 
Fig. 5
Simulation of cyclic compression with confinement (10 MPa) for an aggregate energetic material. 
Fig. 6
Rosette of damage for 1% of longitudinal strain in compressive test (A) and triaxial compressive test with 10 MPa of confinement (B). e_{1} 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 tensioncompression 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
 P. Ladeveze, A. Gasser, O. Allix, Damage mechanisms modeling for ceramic composites, J. Eng. Mater. Technol. 116 (3) (1994) 331–336 [CrossRef] [Google Scholar]
 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]
 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]
 D. Halm, A. Dragon, Y. Charles, A modular damage model for quasibrittle solids interaction between initial and induced anisotropy, Arch. Appl. Mech. 72 (2002) 498–510 [CrossRef] [Google Scholar]
 S. Baste, R. El Guerjouma, G. Alain, Mesure de l’endommagement anisotrope d’un composite céramiquecéramique par une méthode ultrasonore, Revue de Physique Appliquée (Paris) 7 (1989) 721–731 [CrossRef] [EDP Sciences] [Google Scholar]
 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]
 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]
 G.I. Taylor, Plastic strain in metals, J. Inst. Metals 62 (1938) 307–324 [Google Scholar]
 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]
 Z.P. Bazant, P.G. Gambarova, Crack shear in concrete: crack band microplane model, J. Struct. Eng. 110 (9) (1984) 2015–2035 [CrossRef] [Google Scholar]
 Z.P. Bazant, B.H. Oh, Deformation of cracked netreinforced concrete walls, J. Struct Eng. 109 (1) (1983) 93–108 [CrossRef] [Google Scholar]
 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]
 A. Beghini, Z.P. Bazant, Y. Zhou, O. Gouirand, F.C. Caner, Microplane model m5f for multiaxial behavior and fracture of fiberreinforced concrete, J. Eng. Mech. 133 (1) (2007) 66–75 [CrossRef] [Google Scholar]
 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]
 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]
 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]
 E. Kuhl, E. Ramm, R. Borst, An anisotropic gradient damage model for quasibrittle materials, Comput. Methods in Appl. Mech. Eng. 183 (1–2) (2000) 87–103 [CrossRef] [Google Scholar]
 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]
 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]
 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]
 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 [CrossRef] [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
All Figures
Fig. 1
VDT decomposition [14]. 

In the text 
Fig. 2
Effect of the damage parameters on the damage variable evolution. 

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

In the text 
Fig. 4
Simulation of tensile test, compression test and compression with confining pressure (5–10 MPa) for an aggregate energetic material. 

In the text 
Fig. 5
Simulation of cyclic compression with confinement (10 MPa) for an aggregate energetic material. 

In the text 
Fig. 6
Rosette of damage for 1% of longitudinal strain in compressive test (A) and triaxial compressive test with 10 MPa of confinement (B). e_{1} is the loading direction. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.