On Continuum Damage Modeling of Fiber Reinforced Viscoelastic Composites with Microcracks in terms of Invariants

A continuum damage model is developed for the linear viscoelastic behavior of composites with microcracks consisting of an isotropic matrix reinforced by two arbitrarily independent and inextensible fiber families. Despite the fact that the matrix material is isotropic, the model in consideration bears the characteristic of directed media included in the transverse isotropy symmetry group solely due to its fibers distributions and the existence of microcracks. Using the basic laws of continuum damage mechanics and equations belonging to kinematics and deformation geometries of fibers, the constitutive functions have been obtained. It has been detected as a result of the thermodynamic constraints that the stress potential function is dependent on two symmetric tensors and two vectors, whereas the dissipative stress function is dependent on four symmetric tensors and two vectors. To determine arguments of the constitutive functionals, findings relating to the theory of invariants have been used as a method because of the fact that isotropy constraint is imposed on the material. As a result the linear constitutive equations of elastic stress, dissipative stress, and strain energy density release rate have been written in terms of material coordinate description. Using these expressions, total stress has been found.


Introduction
In various fields of industry, fiber reinforced composite materials are commonly used for load carrying components.In particular, these materials have become indispensable in advanced aerospace and aviation structures [1].Fiber reinforced composites can display time dependent behavior (stress relaxation or creep) as a manifestation of one or both constituents deforming in a viscous manner [2].Therefore, it is significant to develop a damage model that also includes the rate dependent nature of the viscoelastic material subjected to external loading.Damage modeling of composites is not an easy task, given the evolution of various damage mechanisms (matrix cracking, fiber breakage, interfacial debonding, transverse ply cracking, and ply delamination) in composite materials [3].The effective stress theory assumes that the damage is uniformly distributed, where in fiber reinforced composites the complex damage occurs in a distributed manner with the damage mechanisms occurring in experimental work [1].Some researches have conducted studies on the formulation of damage modeling on various composites [4][5][6][7][8][9][10].
Continuum damage mechanics (CDM), developed in the last few decades, provides a framework for incorporating the effects of damage induced stiffness softening, anisotropy, and so forth in constitutive equations [11].Phenomenological continuum damage mechanics (CDM) models have been used extensively for modeling creep and fatigue damage by several authors as described in Lemaitre and Chaboche [12], Krajicinovic [13], and Voyiadjis and Kattan [14].Chaboche has proposed a thermodynamic framework for continuum damage mechanics (CDM) for isotropic damage [15].Furthermore, an updated continuum damage model was proposed, in which the fracture energy density, a function of the stress triaxiality, temperature, and strain rate in the transition region were taken as the critical damage factor [16].Li et al. presented a damage model and its numerical solution by means of fast Fourier transforms (FFT) [17].Much progress has been [18] made in the development of a consistent continuum framework for damage mechanics.Darabi  phenomenon in materials that tend to self-heal.This framework extends the well-known Kachanov's [19] effective configuration and the concept of the effective stress space to selfhealing materials by introducing the healing natural configuration in order to incorporate the microdamage healing effects [18].Continuum damage mechanics (CDM) has been proven to be a valuable tool for the determination of properties deterioration due to microcracks of composites.In continuum damage mechanics (CDM) models, families of damage are represented as continuous internal state variables, and concepts of irreversible thermodynamics are used in deriving constitutive equations [20].A number of models of damage behavior of fiber reinforced composites are based on the methods of continuum damage mechanics (CDM).For the case of damage in viscoelastic materials two sets of internal variables are needed.Weitsman [21] presented a formulation to treat this case.In this formulation, damage in the form of microcracks is represented by two symmetric second rank tensors related to the total area of active and passive microcracks in a representative volume element (RVE) of the material and viscoelasticity is introduced through  scalar valued internal state variables (ISVs).Since the scale of molecular rearrangement is much smaller than the scale of microcracks, the effect of damage on viscoelastic retardation time is ignored.The formulation uses Taylor series expansion of Gibbs free energy (Figure 1) [2].
Voyiadjis and Kattan [22] have presented a new formulation to link continuum damage mechanics with the concept of fabric tensors within the framework of classical elasticity assuming small strains.The new theory was called damage mechanics with fabric tensors.A fourth-rank damage tensor was used and its exact relationship to the fabric tensors was illustrated.A [23] model of damage mechanics for directional data was formulated using fabric tensors [22].Qiang et al. have described accurately by a fourth-order fabric tensor of the cracks within the framework of noninteraction approximation the compliance increment due to the presence of generally nonuniformly distributed cracks.They have expressed that the compliance increment can be approximated by the corresponding second-order fabric tensor [24].Spencer has formulated constitutive equations for flow of fabric-reinforced composite materials which show viscous and viscoplastic response at the forming temperature [24,25].
The fabric tensor is a symmetric second rank tensor which characterizes the geometric arrangement of the porous material microstructure.In developing this result it is assumed that the matrix material of the porous elastic solid is isotropic and, thus, that the anisotropy of the porous elastic solid is determined by the fabric tensor [26].Fabric tensor has proved to be an effective tool statistically characterizing directional data in a smooth and frame-indifferent form.Directional data arising from microscopic physics and mechanics can be summed up as tensor-valued orientation distribution functions [27].The subject of fabric tensor has been studied by a number of different researchers [28][29][30][31].There are studies that establish a relationship between the fabric tensor with damage tensor.However, such studies are based on experimental data.In these studies, statistical distribution along a certain plane of the microcracks in the material has been considered and an average density of microcracks has been used.Our study is not based on experimental data and is phenomenological.
In this study in the scope of continuum damage mechanics (CDM), considering representative volume element (RVE), the mechanical representation of damage was expressed with two interior conditional variables showing the properties of a second degree symmetrical tensor.It is assumed that an element from two different continuous fiber families is placed on each point of the composite material.In the scope of continuum mechanics, balance equations have been summarized and the combined form was stated for the energy equation and entropy inequality.Considering necessary constitutive axioms, after determining the arguments affecting the stress potential it was further proceeded to the formulation of the constitutive theory and a model has been formed for the damaged viscoelastic composite with isotropic matrix material.Considering fully orthogonal transformations for the material coordinate system, after determining the common invariants affecting both the stress potential and the dissipative stress, constitutive equations have been obtained related with the "elastic stress, " the "dissipative stress, " and the "strain energy density release rate." On the other hand, inextensibility of the fibers and incompressibility of the composite are acceptable as broad recognition in practice in terms of formulation.Due to some technological requirements, it is aspired that specific construction elements have rather elastic properties, provided that they have high durability in certain directions.Fiberreinforced composite materials are produced sticking fibers in a polymeric matrix which is elastic but with low strength.These fibers are manufactured from high strength graphite or bor.They can be easily bent due to the very small size of their cross section.Assume that inextensibility of the fibers is a reasonable approach since the rigidity of the fibers is very high according to the rigidity of the matrix [32].
The constitutive model developed in this study is applicable for natural or biological fiber composites, such as living tissues, or synthetic composites for implant applications, because the incompressible (nearly incompressible) behavior of composites is often used for natural or biological fiber composites, such as living tissues, or synthetic composites for implant applications.In our study [33], we are concerned with developing the continuum damage mechanics model for elastic behavior of composites having microcracks consisting of an isotropic matrix reinforced by independent and inextensible two families of arbitrarily fibers.

The Representation of Damage
In some researches, in order to be able to define the damage variable, a representative volume element (RVE) has been considered that has a  number of microcracks.While the open or active part of any th microcrack has been shown by A () , its closed or passive surface has been shown by A * () .Active or passive surfaces of a crack can switch positions among each other depending on stress, temperature, and humidity percentage.Despite that, Weitsman states that these open and closed surfaces can be selected as independent variables characterizing the state of a material at a certain time range [21,34].
Stress and strain at the macrolevel are average values over the representative volume element (RVE) volume.Infinitesimal deformations can also be considered among these macrovalues.To fully consider the behaviors of representative volume element (RVE), it is necessary to deal with a  number of crack parameters representing A () and A * () (no sum on ,  = 1, . . ., ) surfaces.Because the real shape of these surfaces is unknown on the mesoscale, assuming them to be equivalent plane surfaces, Weitsman represented them by vectors A () =  () n () and A * () =  * () n () .Here, n ()  stands for a unit normal vector of a microcrack surface [34].However, since each microcrack possesses two equal and opposite surfaces, the description should not discriminate between (+n  ) and (−n  ), which is accomplished by choosing the dyadic product as follows to represent the th flaw.This choice follows an analogous approach by Spencer [35]: In view of the paucity of information regarding the positions and sizes of the individual microflaws, it is proposed to represent damage by the dyadic sums of all H () , ( = 1, 2, . . ., ) and H * () , ( = 1, 2, . . ., ) contained in the representative volume element (RVE), namely, by Obviously, H and H * are symmetric, second rank tensors.As the constitutive variable, in this study we are going to deal with only one damage tensor taking into consideration only the effect of open microsurfaces.In this study, due to the existence of fiber distributions and microvoids in the material it is assumed that the material has gained directed medium characteristics, that is, an anisotropic structure has appeared due to the damage and the fibers.We assume that initially the material was isotropic and that the anisotropy is only caused by the dispersion of microvoids and fibers.For a medium like that, the role of material description vectors will be played by the vector A(X, ) representing the mean values in the representative volume element (RVE) and the vector Ȧ(X, ) representing the change in time of the preceding vector.We believe that, by dividing these vectors by the area of any characteristic surface pertaining to representative volume element (RVE), we render them dimensionless.
On the other hand, because the material will not be able to detect the positive and negative sides of microvoid surfaces, we had previously specified that the dependence on vectors A(X, ) and Ȧ(X, ) can be expressed by a product of tensors: (3)

Kinematics of Fibers Deformation and Thermomechanic Balance Equations
Fiber families are represented by continuous vector fields B(X) and Z(X) before deformation and by continuous vector fields b(x) and z(x) after deformation.Relationships given below are true for  and -fiber families [35,36]: where  and  are, respectively, arc length of fiber before and after deformation,   and   are fiber vector components before deformation,   and   are fiber vector components after deformation,  , =   /  is deformation gradient,   and   are rates of extension of fiber families, and   =  ,  , is Green deformation tensor.
The local thermomechanical balance equations can be summarized as follows [32,37]: conservation of mass in material representation: balance of linear momentum: balance of moment of momentum: conservation of energy: second law of thermodynamics: Here, k stands for the velocity field in a continuous medium,  0 for mass density before deformation,  for mass density after deformation,  ≡ det[ , ] =  0 /(x, ) for Jacobian, k for acceleration,   for stress tensor,   for the mechanical volumetric force per unit of mass,  for internal energy density per unit of mass,   for heat flux vector, ℎ for heat source per unit of mass,  for entropy density per unit of mass, (X, ) for the absolute temperature of a material point  at a moment , and   for permutation tensor.

Thermodynamic Constraints and Modeling Constitutive Equations
Local energy equation ( 9) is then suitably combined with the entropy inequality (10) and, using a Legendre transformation such as  ≡  −  for free energy, entropy inequality is obtained as follows in the material form: Terms relating to the new values have been provided as follows: Here, Σ stands for thermodynamic stress potential,  for generalized free energy density,   for deformation (strain) rate tensor,  , =   /  for the deformation gradient of the reverse motion,   for the stress tensor on material coordinates,   for the heat flux vector on material coordinates, and  , for the temperature gradient on material coordinates.
For the use of inequality (11), which is a general expression of entropy generation, we need to know on which independent variables the thermodynamic potential Σ depends as well as the mechanism of such dependence.Based on the material selected, the arguments of Σ and variables it depends on have been found using constitutive axioms.Using the results of the axioms of causality, determinism, objectivity, smooth neighborhood, and admissibility axioms [32,37], the arguments for Σ allowing two fiber families and viscoelastic structures, bearing cracks, and where these cracks are believed to be changing with time, exposed to thermomechanical loading without heat transfer can be summarized as follows: (  , 0,   , 0,   ,   , ) = 0, In expression (18),   is called as the strain energy density release.The definition   ≡ −  is used in order to deal with a positive value.In expression (22), dissipative stress tensor    was defined with expression (19).   is considered as an energetic stress.All energetic stresses are not elastic stresses, because a part of the energetic stress occurring in the material is spent for the formation of microcracks.But in this study,    is called the elastic stress and it was defined with expression (17).
It is expressly understood from ( 17) and ( 18) that the elastic stress and the strain energy density release rate are derived from the stress potential Σ. From ( 19) it is understood that the dissipative stress is obtained in form of a matrix whose arguments are known.In this case, open forms of Σ and    , which are known as constitutive functions with clear arguments, should be found.
However, firstly we should consider the constraints imposed by the material symmetry axiom onto the material under consideration.Because the symmetry group of the material under consideration is the fully orthogonal group, property [  ] −1 = [  ]  , det S = ±1 is true for the symmetry operation [  ].Therefore, each material point conversion matches an orientation of the material medium.Such conversion should, for every [  ], be in the following form and leave constitutive functionals form invariant.Mathematically this means the validity of the following conversions: On the other hand, both incompressibility of the medium and inextensibility of the fiber are broadly recognized in terms of formulation.Assuming the medium to be incompressible and the fiber to be inextensible, the following requirements should be satisfied [32]: Thus, the constitutive equation for the elastic stress is obtained as follows in material coordinates: In this expression, , Γ  , and Γ  are Lagrange coefficients and are defined by field equations and boundary conditions. −1  ≡  ,  , is the Piola deformation tensor.

Determination of Elastic Stress and Strain Energy Density Release Rate Constitutive Equations
Since the matrix has been assumed to be isotropic, relation ( 24) is expressed as follows; Here, the orthogonal matrix indicating the symmetry group {  } will be expressed for [  ] ∈ (3) and property On the other hand, since Σ has been assumed to be the analytical function of its arguments, such arguments, which are expected to remain invariant under orthogonal transformations belonging to the symmetry group, should depend on a finite number of invariants: Instead of the first three invariants of the Green deformation tensor C, we can use the principal invariants as follows: Given the incompressibility of the composite, inextensibility of the fiber families, and the fact that B and Z are unit vectors, the invariants ,  8 ,  10 ,  5 , and  6 in expressions ( 29) and ( 30) are equal to 1 thus eliminating the dependence of Σ on these invariants.As a result, the invariants on which Σ depends are expressed as follows: Taking the derivative of expression (31) according to   and   and substituting it into ( 27) and ( 18), the following expressions are obtained: It is understood that, as always, repeated indices will undergo summation.If derivatives of invariants appearing in these equations according to   and   are taken from expressions ( 29) and (30) and substituted afterwards, constitutive equation of the elastic stress in nonlinear form is obtained as follows: More concrete form of the constitutive equations given by ( 33) and ( 34) can be obtained provided that Lagrange coefficients −, Γ  , and Γ  and the derivatives of Σ based on its invariants are known.It has been already stated that −, Γ  , and Γ  can be obtained from field equations and boundary conditions.To obtain the derivatives of Σ according to its invariants, how Σ depends on the invariants it is shown to depend on in expression (31) should be estimated.In this study, the matrix material has been considered as an isotropic medium.According to the fact that Σ is an analytical function of those invariants, assuming that this function is analytic, the stress potential is expanded in the power series around natural condition.To obtain a quadratic theory, the terms in this series expanding should be kept to second order; therefore, the stress potential can be represented by a polynomial [32,37].However, the grade and number of terms of the polynomial representing Σ depend on the size of its deformation invariant and their shares of interaction in the case [38][39][40].
In this study, mechanical interactions and effect of damage have been assumed to be linear.Furthermore, considering that the material remains insensitive to directional changes along fibers, double components of fiber vectors have been included in the operation.Because mechanical interactions and effect of damage are assumed to be linear, the stress should remain linear according to the deformation tensor and the damage tensor.Therefore, function Σ could be represented by a second degree polynomial according to the invariants it depends on.On the other hand, because internal energy is defined as a positive definite form, for a polynomial to be positively defined and for the order of invariants not to affect Σ, the polynomial must have symmetric coefficients; that is, it must be in a quadratic form.Accordingly, if polynomial approximation is selected, the following expression can be recorded for the stress potential Σ in terms of the existing invariants: Σ = Σ .       , (,  = 1, 2, 4, 7, 9, 11, . . ., 21) ,   =   . ( The derivatives of Σ based on its invariants in ( 33) and ( 34) are obtained from expression (35) as follows: At this stage, expressions ( 29) and ( 30) have shown on what the invariants in expression (36) depend.Due to the existence of the relationship   =   + 2  between the Green deformation tensor and the strain tensor and assuming that mechanic interactions are linear (  ≅ Ẽ = (1/2)( , + , )), those invariants that depend on the Green deformation tensor (  ) have been expressed in terms of strain tensor (  ), which is a more useful parameter.
Terms after the third term on the right side of ( 33) and all terms of the right side of (34) have been calculated using the partial derivatives given in expression (36) and invariants that depend on the strain tensor (  ).Due to the assumptions made in this study, of the first grade components of the strain tensor   and the damage tensor   and of the external multiplication components of vectors   and   , those whose number is even have been taken into consideration.Thus, in the beginning, the elastic stress is expressed for the condition without stress and without load (with the term  1   assumed to be zero) by taking common coefficients into common parenthesis: The strain energy density release rate is expressed as follows by taking common coefficients into parenthesis in the beginning, without microcracks (the term  1   is taken here as zero) being obtained as follows: Coefficients [  ( = 1, 2, 3, . . ., 22) and   ( = 1, 2, 3, . . ., 20)] in (37) and (38) have been depending on the medium temperature  and   .
In a composite material that consists of an isotropic matrix reinforced by two arbitrary independent and inextensible fiber families, the medium is assumed to be incompressible and homogeneous, has microcracks, and shows linear viscoelastic behavior.Equation (37) is the linear constitutive equation of elastic stress.First, second and third terms of (37) are hydrostatic pressure and contributions of fiber tensions to the elastic stress, respectively; fourth and eighth terms combined are the contribution of the elastic deformation; fifth and twenty sixth terms combined are the contribution of the damage tensor; sixth, tenth, and twenty-third terms are the stress arising of the interaction between the fiber distribution B and the damage tensor; seventh, fifteenth, and twenty-fourth terms are the stress arising of the interaction between the fiber distribution Z and the damage tensor; ninth and thirteenth terms are the stress arising of the interaction between the fiber distribution B and the elastic deformation; eleventh, sixteenth, twenty-first, and twenty-fifth terms are the contribution produced by the triple interaction between the fiber fields B and Z and the damage tensor; twelfth term is the contribution of fiber distribution B; fourteenth and eighteenth terms are the stress arising of the interaction between the fiber distribution Z and the elastic deformation; seventeenth term is the contribution of fiber distribution Z; nineteenth term is the stress produced by the interaction between the fiber field B and the fiber field Z; and twentieth and twenty second terms are the contribution produced by the triple interaction between the fiber fields B and Z and the elastic deformation field.
Equation (38) is the linear constitutive equation of strain energy density release rate.First and nineteenth terms combined are the contribution of the elastic deformation; second term is the contribution of the damage tensor; third and seventh terms combined are the strain energy density release arising of the interaction between the fiber distribution B and the damage tensor; fourth and twelfth terms are the strain energy density release arising of the interaction between the fiber distribution Z and the damage tensor; fifth term is the contribution of fiber distribution B; sixth and ninth terms are the strain energy density release arising of the interaction between the fiber distribution B and the elastic deformation; eighth, thirteenth, and seventeenth terms are the contribution produced by the triple interaction between the fiber fields B and Z and the damage tensor; tenth term is the contribution of fiber distribution Z; eleventh and fourteenth terms are the strain energy density release arising of the interaction between the fiber distribution Z and the elastic deformation; fifteenth term is the strain energy density release produced by the interaction between the fiber field B and the fiber field Z; and sixteenth and eighteenth terms are the contribution produced by the triple interaction between the fiber fields B and Z and the elastic deformation field.

Determination of Dissipative Stress Constitutive Equation
It is assumed that the viscoelastic behavior of the medium in consideration is in conformity with Kelvin-Voigt model.It has been determined that the dissipative stress depends on deformation, deformation rate, damage, damage rate, and distributions of fibers, yielding in expressions ( 19)- (21).Additional constraints imposed on the dissipative stress by constitutive functions originate from the material symmetry of the medium.The structure of the dissipative stress should be in compliance with the following transformation for each orthogonal matrix [  ] ∈ (3) belonging to the symmetry group of the material: where the matrix is isotropic; relation (39) is valid for each orthogonal matrix of the fully orthogonal group.The dissipative stress is an isotropic function of the symmetric matrices   , Ċ  ,   , Ḣ  and polar vectors   and   .For simplicity of notation, dependence of the dissipative stress tensor on  has not been denoted.To obtain the explicit expression of the tensor component    in terms of its invariant arguments, the following way has been followed according to the theory of invariants [41]:   is an arbitrary vector,    and the vector   are multiplied on the right and on the left by scalar multiplication, and the product is defined by a scalar function R.
Here the R scalar is an isotropic function of the symmetric matrices   , Ċ  ,   , Ḣ  and absolute vectors   ,   , and   .Taking the partial derivative of expression (40) according to   and   , the following expression can be recorded [32]: Because the left side of this expression is independent from the vector V, V=0 should be true for inequality (41).In this situation, isotropic tensor function    is expressed as follows: In this situation, to find the tensor    from relationship (42), structure of the scalar R should be determined that depends on the arguments   , Ċ  ,   , Ḣ  ,   ,   ,   and second degree partial derivative of this function according to the vector V should be calculated at V = 0. Let us first remove the arbitrary vector V from the arguments of the scalar function R and define a scalar function with arguments   , Ċ  ,   , Ḣ  ,   ,   : For Ξ, which is an isotropic function, to keep the invariant under orthogonal coordinate transformations, its arguments must depend on a finite number of invariants.Using the methods in the theory of invariants [41], 98 invariants of the four symmetric tensors   , Ċ  ,   , Ḣ  and the two polar vectors   ,   , which are independent of each other, have been expressed in the following list: However, the main function that needs to be obtained is in the form of arguments of the scalar isotropic function R, that is,   , Ċ  ,   , Ḣ  ,   ,   ,   .Because the scalar R is a bilinear function of the vector V, it depends on the following invariants in addition to those provided in (44): Because mechanical interactions and effect of damage have been assumed to be linear, coefficients in (47) are each a scalar function of invariants that do not contain square or higher grade terms of tensors C and Ċ or terms in the forms (C Ċ), (CH), (C Ḣ), ( ĊH), ( Ċ Ḣ), and (H Ḣ) in (44).Besides, taking into account that values of the invariants  36 and  57 are equal to 1 due to the inextensibility of fiber families and values of the invariants  33 and  34 are equal to 1 because they are unit vectors pertaining to the distribution of fibers B and Z before deformation, invariants on which the said coefficients depend have been recorded as follows in terms of (  ) and ( Ė  ).
Coefficients in (47) can be recorded as a scalar function of the invariants in (48) can be expressed as follows.
Here the following definitions have been used.Using expressions (49) in (47), expressing the tensors   and Ċ  in terms of the tensors   and Ė  , taking mechanical interactions, effect of damage interactions, and assumptions made concerning the fibers into account, the following has been expressed: Kelvin-Voigt body composing a viscoelastic model is defined as a continuum model where stress depends on the rate of deformation together with the deformation.The stress constitutive equation of the linear Kelvin-Voigt body has been expressed as   = Σ  (, x)∈  + Γ  (, x) ∈  .Σ  coefficient bears the symmetry characteristics given as Σ  = Σ  = Σ  = Σ  , whereas Γ  coefficient bears the symmetry characteristics given as Γ  = Γ  = Γ  .This symmetry characteristic decreases the number of components of material module Σ  to 21 and the number of components of material module Γ  to 36.Besides, as a result of Onsager principle in the linear reversible thermodynamics, it has been assumed that relaxation tensor Γ  has a symmetry characteristic as Γ  = Γ  .According to this symmetry condition, the number of components pertaining to material module Γ  is thus reduced to 21.

Concluding Remarks
This paper presents a continuum damage model based on fundamental concepts of continuum mechanics for the linear viscoelastic behavior of incompressible composites with microcracks that consist of an isotropic matrix reinforced by inextensible two families of fibers having an arbitrary distribution."Damage" is expressed by two symmetric, second rank tensors which are related to the total areas of "active" and "passive" microcracks within a representative volume element of the multifractured material.The matrix material has been assumed to be an isotropic medium; however, due to the distributions of fibers and the existence of microcracks, it has gained the property of a directed object, thus gaining the appearance of an anisotropic structure.Inextensibility of the fibers is an assumption which is widely used in practice.Thus, fiber families are assumed to be inextensible and composite medium is assumed to be incompressible.In this context, the composite expresses itself behaviorally in terms of the elastic stress, the dissipative stress, and the strain energy density release rate.To obtain a more concrete expression of nonlinear constitutive equations of the elastic stress and the strain energy density release rate given by expressions (33) and (34), derivatives of Σ must be known according to the arguments it depends on.Thus, stress potential Σ has been represented by a second degree polynomial and its derivatives according to its invariants have been calculated.In this study, mechanical interactions and effect of damage have been assumed to be linear.Furthermore, since the matrix material has to remain insensitive to directional changes along fibers, even-numbered exterior products of vector fields representing fiber distributions have been considered.The linear constitutive equations of the elastic stress and the strain energy density release rate have been given by expressions (37) and (38).From these equations, it can be seen that the deformation field, the damage tensor, distribution of fibers, and interactions of them all contribute to the creation of elastic stress and strain energy density release rate.
Equation (47) obtained by using the invariants of the arguments depends on the dissipative stress, since the dissipative stress is a vector-valued isotropic function.Coefficients in this equation have been expressed in terms of the invariants on which they depend and each term has been calculated.When these calculations are made, the mechanical interaction and the effect of damage have been assumed to be linear and even number vector components of fiber vectors have been included in the operations, since the composite remains indifferent to changes of direction along the fibers.The linear constitutive equation of the dissipative stress has been expressed by (52).As understood from this equation, separate and common interactions of the fiber distribution fields with the rate of deformation and the change in time of the damage tensor contribute to the formation of dissipative stress.Substituting (37) and (52) into the expression (22), the total stress has been obtained in (53).
This paper is concerned with developing the continuum damage mechanics model for viscoelastic behavior of composites having microcracks consisting of an isotropic matrix reinforced by independent and inextensible two families of arbitrarily fibers.The elastic stress and strain energy density release rate are expressed in terms of the thermodynamic stress potential, a function of the left Cauchy-Green tensor, the damage tensor, and the fiber distributions.The dissipative stress depends on these quantities and the rate of deformation and the change in time of the damage tensor.The elastic stress and strain energy density release rate and dissipative stress are treated in separate sections.The appropriate invariants used as arguments for thermodynamic stress potential function and dissipative stress function are introduced and the elastic stress, strain energy density release rate, and dissipative stress constitutive equations are worked out.This is followed by specializations for incompressibility, the special case when thermodynamic stress potential is a quadratic polynomial in the invariants and a linearization based on small strains.
The resulting expressions show that the first is attributable to the "inherent" viscoelastic behavior of the undamaged material, while the second is due to the change in time of the damage tensor.As can be noted from constitutive equations for elastic stress, strain energy density release rate, and dissipative stress, the coupling of damage and viscoelasticity introduces very significant complexities in the composite material's response and substantial work remains to be done, both experimentally and analytically, to attain a quantitative and practical understanding of the phenomenon.
In this paper, as a result the constitutive equations relating to the elastic stress, strain energy density release rate, and dissipative stress have been written in terms of material coordinate descriptions.This paper develops a mathematical model based on methods of continuum damage mechanics.After this paper, practical problems will be solved by forming B(X) and Z(X) vector fields for various fiber distributions whose parametric equations in the material medium are in the form of X=X() and necessary interpretations will be made in a more concrete way.Also, in a future work we will study the development of numerical methods for this model.