Scalar, vectorial, and tensorial damage parameters from the mesoscopic background/Mesoskoopse taustaga skalaarsed, vektoriaalsed ja tensoriaalsed purunemisparameetrid.1. INTRODUCTION
1.1. Phenomenological definitions of damage parameters
Numerous damage models have incorporated scalar, vectorial, or tensorial damage variables that can be characterized at the macro-scale, for example, by the change in compliance. A scalar damage parameter was introduced [1-4] to account for the decrease in the stiffness of the material with progressing damage. Two scalar damage parameters were proposed  to account independently for the change in hydrostatic energy and the remaining part of the elastic energy with increasing damage. A different reason for introducing two scalar damage parameters was the healing of cracks under compression . In composites and fibre reinforced materials it is reasonable to introduce independent scalar order parameters for the prescribed directions, given by the fibre orientation.
A second order damage tensor was defined , accounting for the reduction of the effective surface area that transmits forces. The resulting effective stress is expressed in terms of the damage tensor . For definitions of a second order damage tensor see also [3,8-11]. For parallel microcracks a second order damage tensor was associated with the dyadic product of crack orientation with itself times a scalar parameter . This definition coincides with our definition from the mesoscopic point of view in the special case of parallel microcracks.
A fourth order damage tensor has been introduced. It can be understood as mapping the elastic tensor (in a linearized strain theory) of the virgin material to the elastic tensor of the damaged material, or as mapping the respective stiffness tensors. For a summary of damage parameters of different orders see also [13,14].
For a constitutive theory of damaged materials with a thermodynamic background see [2,15]. A thermodynamic theory of damage, including the interpretation of failure as loss of thermodynamic stability, can be found in . For a comparison to experimental results see .
An alternative choice of damage variable is one that incorporates salient aspects of damage morphology in its definition. Such 'micro-mechanically inspired' damage models involving scalar, tensor, or 'Fabric tensor' representations of damage were introduced in the study of heterogeneous materials containing voids or various crack-like surface discontinuities [18-23].
Our aim here is to show how damage parameters of different tensor order can be defined from the mesoscopic background. The different damage parameters correspond to different levels of macroscopic approximation of the mesoscopic distribution of crack sizes and orientations. As it was shown in [24,25] in the case of a scalar damage parameter, the mesoscopic theory leads not only to the definition of damage parameters, but also to equations of motion for them. On the example of the free energy density we will show the general form of a constitutive equation for the different choices of a damage parameter. In the case of a rotation symmetry of crack orientations, the different forms of the constitutive equation can be reduced to a form with two scalar parameters: the average crack size and a scalar orientational order parameter.
1.2. Mesoscopic theory of complex materials and application to material damage
The mesoscopic theory was developed in order to deal with complex materials within continuum mechanics . The idea was to enlarge the domain of the field quantities by an additional variable, characterizing the internal degree of freedom connected with the internal structure of the material. In a simple model the microcrack is described as a flat, rotation symmetric surface, a so-called penny-shaped crack. The single crack is characterized then by the crack radius l and the orientation of the crack surface normal, a unit vector n, see Fig. 1. In the mesoscopic theory l and n are the additional variables, i.e. mesoscopic field quantities depend on position x of the continuum element, time t, crack radius l, and crack orientation n. We will abbreviate this set of mesoscopic variables with (*).
In addition we make here the following simplifying assumptions:
1. The diameter of the cracks is much smaller than the linear dimension of the continuum element. Under this assumption the cracks can be treated as an internal structure of the continuum element. The cracks are assumed to be small enough so that there is a whole range of crack sizes and orientations in the volume element.
2. The cracks are fixed to the material. Therefore their motion is coupled to the motion of representative volume elements.
3. The cracks cannot rotate independently of the material, i.e. the rotation velocity is determined by the antisymmetric part of the time derivative of the deformation gradient of the surrounding material.
4. The number of cracks is fixed, there is no production of cracks, but very short cracks preexist in the virgin material.
5. The cracks cannot decrease in area, but can only enlarge, meaning that cracks cannot heal.
[FIGURE 1 OMITTED]
To summarize our model, the microcrack is characterized by a unit vector n representing the orientation of the surface normal and by the radius l of the spherical crack surface. These parameters will be taken as the additional variables in the mesoscopic theory.
Beyond the use of additional variables the mesoscopic concept introduces a statistical element, the so-called mesoscopic distribution function. In our case this is a distribution of crack lengths and orientations in the continuum element at position x and time t, called here crack distribution function (CDF). The distribution function is the probability density of finding a crack of length l and orientation n in the continuum element. The elements are material elements, including the same material and the same cracks for all times. Macroscopic quantities are calculated from mesoscopic ones as averages over crack sizes and crack orientations.
1.3. Mesoscopic balance equation of crack number
Field quantities such as mass density, momentum density, angular momentum density, and energy density are defined on the mesoscopic space. For distinguishing these fields from the macroscopic ones we add the word 'mesoscopic'. In addition to mass density we introduce the crack number density N as the density of an extensive quantity. The mesoscopic crack number density N(l, n, x, t) is the number density, counting only cracks of length l and orientation n.
1.3.1. Balance of crack number
In our model the cracks move together with the material element. Therefore their flux is a convective flux, having a part in position space, a part in orientation space, and a part in the length interval. There is no production and no supply of crack number. Therefore we have for the crack number density N:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1.1)
We have used spherical coordinates for the mesoscopic variables crack length l [member of] [0, [infinity]] and crack orientation n [member of] [S.sup.2], and we represent the divergence with respect to the mesoscopic variables in spherical coordinates. The covariant derivative on the unit sphere is denoted by [[nabla].sub.n] and the material velocity by v. In our model all cracks within the continuum element move with this velocity. We designate by u(x) = [??] the orientation change velocity, which is not the same for all cracks in the continuum element. It is related to the angular velocity [omega](x, t) by the relation
u(x) = [omega] x n. (1.2)
This angular velocity is the same for all cracks in the element. It is determined by the rotation of the surrounding material.
1.4. Definition of the distribution function and equation of motion
Because of its definition as probability density the distribution function is the number fraction
f(l,n,x,t) = N(l,n,x,t) / N(x,t) (1.3)
in volume elements, where the number density N(x, t) is non-zero. Here N(x, t) is the macroscopic number density of cracks of any length and orientation. As the distribution function in equation (1.3) is not well defined if N(x, t) = 0, we define in addition that in this case f (l, n, x, t) = 0. As there is no creation of cracks in our model the distribution function will be zero for all times in these volume elements. In all other volume elements with a nonzero crack number it is normalized
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1.4)
With respect to crack length it is supposed that the distribution function has a compact support, meaning that in a sample there cannot exist cracks larger than the sample size.
We obtain from the mesoscopic balance of crack number density a balance of the CDF f(l,n,x,t), by inserting its definition:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1.5)
The right hand side is equal to zero, as for the co-moving observer the total number of cracks in a volume element does not change in time.
A growth law for the single crack l is needed in equation (1.5). For example the Rice-Griffith dynamics, which is motivated from macroscopic thermodynamic considerations, can be applied.
In  the mesoscopic theory was specialized to damaged material with penny-shaped cracks. The balance equations and the differential equation for the crack size distribution function were derived. Using the Rice-Griffith differential equation for the size of a single crack, the time evolution of the whole distribution of cracks under load, as well as the evolution of the average crack size, was investigated . In  two different growth laws for the single crack under load were considered. Finally, the dynamics of a second order damage tensor was derived in  under the assumption of a simplified single crack growth law under an effective stress.
2. MESOSCOPIC DEFINITIONS OF DAMAGE PARAMETERS OF DIFFERENT ORDERS
Definitions of a scalar damage parameter, a vectorial damage parameter, and a damage tensor are given, based on the mesoscopic distribution function. A possible choice of a scalar damage parameter is the average crack length. However, a scalar parameter is not sufficient, because the crack growth introduces an anisotropy into the material. In order to account for this anisotropy, it is necessary to define a vectorial or a tensorial damage parameter. Starting from the mesoscopic distribution function, the more natural way is to define a damage tensor of second order. This second order tensor is the second moment of the orientation distribution function, depending on the crack orientation vector. A case of special interest is a distribution with a rotation symmetry, the uniaxial case. In this case the damage tensor can be expressed in terms of a scalar quantity and a unit vector, which is the orientation of the rotation symmetry axis. In this case we can define easily a damage vector from the second order damage tensor. This damage vector has the orientation of the rotation symmetry axis. A representation of the second order damage tensor in the general case without rotation symmetry is also given. It is shown how a vectorial damage parameter can be defined in general without rotation symmetry in terms of eigenvectors of the damage tensor.
2.1. Damage parameter of second order
We define the second order damage tensor as the second orientational moment of the distribution function:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2.1)
where [??] denotes the symmetric traceless part of the dyadic product, and D is a second order symmetric traceless tensor. This definition of the second order damage parameter accounts for the crack length distribution as well as for the orientation distribution.
2.2. Vectorial damage parameter defined from the second order tensor
Due to the symmetry the second order tensor damage parameter D has a spectral decomposition with orthonormal eigenvectors d, e, and f and eigenvalues [delta], [epsilon] and [phi]:
D = [delta]dd + [epsilon]ee + [phi]ff. (2.2)
Because D is traceless we have
[delta] + [epsilon] + [phi] = 0. (2.3)
Therefore, not all eigenvalues can have the same sign. The following cases concerning the signs of the eigenvalues are possible:
1. One eigenvalue is positive and two eigenvalues are negative, for instance
[delta] > 0, [epsilon] < 0, [phi] < 0. (2.4)
In this case we chose the eigenvector (here d) corresponding to the single positive eigenvalue as the unit vector defining the orientation of the vector damage parameter.
2. One eigenvalue is negative and two eigenvalues are positive, for instance
In this case we chose the eigenvector (heref) corresponding to the single negative eigenvalue as the unit vector defining the orientation of the vector damage parameter.
3. All eigenvalues are zero:
[delta] = 0, [epsilon] = 0, [phi] = 0. (2.6)
In this case D = 0 and we have an isotropic orientation distribution. In this case no vector damage parameter can be defined.
4. One eigenvalue is zero, and the two others have opposite signs, for instance:
[delta] = 0, [epsilon] > 0, [phi] = -[epsilon] < 0. (2.7)
In this very special case we could define a vector damage parameter, having the orientation of the wedge-product of the two eigenvectors.
The length of the damage vector can be defined as the absolute value of the corresponding eigenvalue. The definition of the damage vector in terms of an eigenvector naturally leads to the symmetry of an orientation, namely the damage vector and the reversed one cannot be distinguished.
2.3. Special case of the uniaxial distribution function
If there exists a rotation symmetry axis of the distribution function, two eigenvalues coincide, either the two positive ones, or the two negative ones. In both cases the tensor damage parameter is of the form:
D = (lS(l)) [??] = (lS(l)) (dd - 1/3 1) = (lS(l)) (dd = 1/3 (dd + ee + ff)) (2.8)
with the unit tensor 1 and a scalar parameter S, denoted as scalar orientational order parameter. The unit vector d is the orientation of the rotation symmetry axis. S(l) is a measure of the degree of parallel order of the cracks. It is zero if the orientations are distributed isotropically and has the value 1 in the case all cracks are oriented parallel. The orientational order can be different for different crack sizes, therefore S is a function of crack radius l. The average () here is the average over all crack lengths:
[FIGURE 2 OMITTED]
(lS(l)) [[integral].sup.[infinity].sub.0] = f(l,x,t)lS(l)[l.sup.2]dl. (2.9)
For the eigenvalues this corresponds to
[delta] = (lS(l)) (1 - 1/3) = 2/3 (lS(l)), (2.10>
[epsilon] = - 1/3 (lS(l)), (2.11)
[phi] = 1 1/3 (lS(l)), (2.12)
For positive values of S we have one positive eigenvalue and two negative ones. For negative values of S two eigenvalues are positive and one is negative. In both cases the definition of the vector damage parameter given in the previous section leads to the eigenvector d as the orientation of the damage vector. In the case of rotation symmetric orientation distributions this is the orientation of the rotation symmetry axis. The case of positive values of S corresponds to a distribution where the crack-normals are more or less parallel to the rotation symmetry axis. For negative values of S crack orientations are concentrated in a plane perpendicular to the rotation symmetry axis (see Fig. 2).
For the damage vector we find in the uniaxial case:
[??] = 2/3 (lS(l))d. (2.13)
It depends on the degree of orientafional order and on the average crack length.
2.4. Case of a small deviation of the distribution function from rotation symmetry
If the deviation of the orientation distribution from rotation symmetry is small, the two eigenvalues of equal sign differ only by a small amount s, and the damage tensor is of the form:
D = (2/3 (S(l)l)dd + (- 1/3(S(l) -s(l))l)) ee + (- 1/3((S(l) - s(l))l)) ff). (2.14)
In this case we can still define the damage vector the sam///e way as in the uniaxial case////:
[??] = 2/3 (lS(l))d. (2.15)
The scalar order parameter S is a measure of the degree of order and the biaxiality parameter s is a measure of the deviation of the orientation distribution from rotation symmetry.
2.5. Scalar damage parameters
One possible definition of a scalar damage parameter is the average crack length:
D = [[integral].sup.[infinity].sub.0] If (l,x, t)[l.sup.2]dl. (2.16)
In the rotation symmetric case another scalar measure of damage is
[D.sub.S] = (lS(l)) = [[integral].sup.[infinity].sub.0] lS(l)f(l,x,t)[l.sup.2]dl. (2.17)
It is the average crack radius, projected onto the plane perpendicular to the rotation symmetry axis. Under the assumption that only a load applied perpendicular to the crack surface causes crack growth and that a stress vector component in the crack plane has no effect, it is reasonable that the damage parameter (2.17) is relevant for describing the progressive damage under external load in the case of an anisotropic crack distribution.
Another possible choice of a scalar damage parameter would be the number fraction of cracks with radius l, exceeding a certain critical length L, [[integral].sup.[infinity].sub.L] f[l.sup.2]dl. A comparison of results obtained with the different choices of scalar damage parameters is left for a future work.
3. EXAMPLES OF CONSTITUTIVE FUNCTIONS FOR THE DIFFERENT DAMAGE DESCRIPTIONS
As an example of a constitutive function we will consider the free energy density. We will start out with a representation theorem for the free energy depending on the second order damage tensor. Then we will show how the constitutive equation simplifies in the case of a rotation symmetric distribution function.
We will assume that constitutive quantities depend on the equilibrium variables strain tensor S and temperature and in addition on the damage parameter. The temperature dependence will not be denoted explicitly, as it is a scalar quantity. All material coefficients may depend on temperature.
3.1. Free energy as a function of strain and damage in the case of a second order damage tensor
The most general polynomial form of the energy density up to second order in each variable is given by a representation theorem [30,31]:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3.1)
This form is the simplest and natural extension of linear elasticity considering the damage. The coefficients still can be arbitrary functions of temperature.
In the case of rotation symmetry we have:
D = [D.sub.S][??] (3.2)
and the scalar products can be calculated as:
D x D = [D.sup.2.sub.S] [??] x [??] = [D.sup.2.sub.S] (dd - 1/3 1) x (dd - 1/3 1) = [D.sup.2.sub.S] (1/3 -dd + 1/9 1) (3.3)
tr(D x D) = 2/3 [D.sup.2.sub.S). (3.4)
The expression for the free energy simplifies to:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]
The free energy density is expressed here in terms of the vector d (the rotation symmetry axis) and the scalar damage parameter [D.sub.s].
3.2. Special case of uniaxial strain in the z-direction and symmetry axis of the distribution in the same direction
This situation occurs (approximately) in a uniaxial tension experiment (see Fig. 3).
The assumption that in all volume elements the CDF is rotation symmetric with the z-direction as symmetry axis is an approximation, valid for small deformations. In the case of large deformations, the rotation of the volume element cannot be neglected, and the cracks rotate with the material element. This leads to a rotation of the symmetry axis of the orientation distribution function in the volume element, which depends on the position of the volume element. The (local) symmetry axis of the CDF does not coincide with the direction of the applied strain anymore.
[FIGURE 3 OMITTED]
In the geometry with the global rotation symmetry around the z-axis the only interesting components of tensors are the z-z-components and traces. For the strain and the damage vector we have:
S = [[epsilon].sub.zz][e.sub.z][e.sub.z] =: [epsilon] [e.sub.z][e.sub.z], (3.6) D = [e.sub.z] (3.7)
and the free energy density reduces to a function of three scalar quantities [epsilon], D, [D.sub.S], where D is the average crack length (in any direction) and [D.sub.S] is the average crack length in a plane orthogonal to the z-direction. Its value depends on the anisotropy of the average crack length distribution.
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3.8)
C. P. thanks the Deutsche Forschungsgemeinschaft (DFG), contract PA 410/5-1, for financial support. This research was supported by OTKA T048489, by the EU-I3HP project, and the Bolyai scholarship to P. V
Received 5 February 2008, in revised form 13 May 2008
[1.] Kachanov, L. M. On the time to failure under creep conditions. Izv. AN SSSR, Otd. Tekhn. Nauk, 1958, 8, 26-31.
[2.] Lemaitre, J. A Course on Damage Mechanics. Springer-Verlag, Berlin, 1996.
[3.] Kachanov, M. Continuum model of medium with microcracks. J. Eng. Mech. Diu, 1980, 106(EMS), 1039-1051.
[4.] Murakami, S. Notion of continuum damage mechanics and its application to anisotropic creep damage theory. J. Eng. Mat. Technol., 1983, 105, 99.
[5.] Ju, J. W. On energy-based coupled elastoplastic damage theories: Constitutive modelling and computational aspects. Int. J. Solids Struct., 1989, 25(7), 803-833.
[6.] Fremond, M. and Nedjar, B. Damage, gradient of damage and principle of virtual power. Int. J. Solids Struct., 1996, 33(8), 1083-1103.
[7.] Murakami, S. and Ohno, N. A continuum theory of creep and creep damage. In 3rd Creep in Structures Symposium, Leicester (Ponter, D. R. and Hayhorst, A. R., eds). NTAM, Springer, Berlin, 1980, 422-443.
[8.] Lecki, F. A. and Onat, E. T. Physical nonlinearities in structural analysis. In Tensorial Nature of Damage Measuring Internal Variables. Springer-Verlag, Berlin, 1981, 140-155.
[9.] Chaboche, J.-L. Development of continuum damage mechanics for elastic solids sustaining anisotropic and unilateral damage. Int. J. Damage Mech, 1993, 2, 311-329.
[10.] Kanatani, K.-I. Distribution of directional data and fabric tensors. Int. J. Eng. Sci., 1984, 22(2), 149-164.
[11.] Rizzi, E. and Carol, I. A formulation of anisotropic elastic damage using compact tensor formalism. J. Elasticity, 2001, 64(23),85-109.
[12.] Betten, J., Sklepus, A. and Zolochevsky, A. A constitutive theory for creep behavior of initially isotropic materials sustaining unilateral damage. Mech. Res. Comm., 2003, 30, 251-256.
[13.] Maire, J. F. and Chaboche, J. L. A new formulation of continuum damage mechanics (cdm) for composite materials. Aerospace Sci. Technol., 1997, 4, 247-257.
[14.] Cauvin, A. and Testa, R. B. Damage mechanics: basic variables in continuum theories. Int. J. Solids Strucz, 1999, 36(5), 747-761.
[15.] Murakami, S. and Kamiya, K. Constitutive and damage evolution equations. Int. J. Mech. Sci., 1997, 39(4), 47386.
[16.] Van, P. Internal thermodynamic variables and the failure of microcracked materials. J. Non-Equil. Thermodyn., 2001, 26(2), 167-189.
[17.] Van, P. and Vasarhelyi, B. Second law of thermodynamics and the failure of rock materials. In Rock Mechanics in the National Interest Vl (Tinucci, J. P., Elsworth, D. and Heasley, K. A., eds). Balkema Publishers, Lisse-Abingdon-Exton(PA)-Tokyo, 2001, 767-773. Proceedings of the 9th North American Rock Mechanics Symposium, Washington, USA, 2001.
[18.] Nemat-Nasser, S. and Hori, M. Micromechanics: Overall Properties of Heterogeneous Materials. North-Holland, Amsterdam, 1993.
[19.] Krajcinovic, D. and Fonseka, G. U. The continuous damage theory of brittle materials. Part 1: General theory. J. Appl. Mech, 1981, 48, 809-815.
[20.] Krajcinovic, D. Damage Mechanics. Elsevier, Amsterdam, 1996. North-Holland Series in Applied Mathematics and Mechanics.
[21.] Krajcinovic, D. and Silva, M. A. G. Statistical aspects of the continuous damage mechanics. Int. J. Solids Strucz, 1982, 18, 551-562.
[22.] Luo, D. M., Takezono, K. and Tao, S. The mechanical behavior analysis of cfcc with overall anisotropic damage by the micromacro scale method. Int. J. Damage Mech, 2003, 12(2), 141-162.
[23.] Rizzi, E. and Carol, I. Dual orthotropic damage-effect tensors with complementary structures. Int. J E'ng. Sci., 2003, 41(13-14), 1445-1495.
[24.] Papenfuss, C. Damage evolution in micro-cracked materials under load. In Trends in Continuum Physics (Maruszewski, B. T., Muschik, W. and Radowicz, A., eds). World Scientific, 2004, 223-233.
[25.] Papenfuss, C., Bohme, T., Herrmann, H., Muschik, W. and Verhas, J. Dynamics of the size and orientation distribution of microcracks and evolution of macroscopic damage parameters. J. Non-E'quilib. Thermodyn., 2007, 32(2), 129-143.
[26.] Papenfuss, C. Theory of liquid crystals as an example of mesoscopic continuum mechanics. Comp. Mater. Sci., 2000, 19, 45-52.
[27.] Van, P., Papenfuss, C. and Muschik, W Mesoscopic dynamics of microcracks. Phys. Rev. E., 2000, 62(5), 6206-6215.
[28.] Van, P., Papenfuss, C. and Muschik, W. Griffith cracks in the mesoscopic microcrack theory. J. Phys. A, 2004, 37(20), 5315-5328. Published online: Condensed Matter, abstract, cond-mat/0211207; 2002.
[29.] Papenfuss, C., Van, E and Muschik, W Mesoscopic theory of microcracks. Arch. Mech, 2003, 55(5-6), 481-499.
[30.] Smith, G. F. On isotropic integrity bases. Arch. Ration. Mech. An., 1964, 17, 282-292.
[31.] Pipkin, A. C. and Rivlin, R. S. The formulation of constitutive equations in continuum physics 1. Arch. Ration. Mech. An., 1959, 4, 129-144.
Christina Papenfuss (a) * and Peter Van (b)
(a) Technische Universtat Berlin, Strasse des 17. Juni 135, 10623 Berlin and Technische Fachhochschule Berlin, Luxemburger Str. 10, 13353 Berlin, Germany
(b) KFKI Research Institute for Particle and Nuclear Physics, H-1525, Budapest, P.O. Box 49, Hungary; email@example.com
* Corresponding author, firstname.lastname@example.org