Overall moduli of solids with microcracks: Load-induced anisotropy

For a linearly elastic brittle solid containing microcracks that may be closed or may undergo frictional sliding, a general method is developed for estimating the overall instantaneous moduli which depend on the loading conditions. When the cracks are all open and when they are randomly distributed, then the overall response is isotropic. The moduli for this case have been obtained by Budiansky and O'C onnell (1976). On the other hand, when some cracks close, and when some closed cracks undergo frictional sliding, then the overall response becomes anisotropic and dependent on the loading conditions, as well as on the loading path. The self-consistent method is used to estimate the overall moduli. The effects of crack closure and loadinduced anisotropy are included. Several illustrative examples are worked out, showing the important influence of the load path on the overall response when crack closure and frictional sliding are involved.     

Overall moduli and constitutive relations of bodies containing multiple bridged microcracks

The overall moduli and non-linear constitutive relations of solids containing randomly oriented and parallel multiple bridged microcracks are studied in this paper. For both linear and non-linear bridging laws, two non-dimensional mixed parameters which play an important role in determining the overall moduli and constitutive properties are identified. It is shown that the overall moduli of a body containing bridged cracks can be derived from those of the body containing unbridged cracks of the same configuration with the crack density parameter being discounted by these two non-dimensional parameters. For non-linear bridging, attention is paid to the power-law bridging which represents the existing bridging laws in the literature. Particular emphasis is placed on the widely used square-root type bridging law. In this case, closed-form expressions for the overall non-linear constitutive relations are obtained. For a general power-law bridging, it is found that the overall constitutive behaviour of the cracked body is governed by Ramberg–Osgood type relations under strong bridging.     

Electromechanical response of (3–0, 3–1) particulate,fibrous, and porous piezoelectric composites with anisotropic constituents: A model based on the homogenization method

An analytical framework based on the homogenization method has been developed to predict the effective electromechanical properties of periodic, particulate and porous, piezoelectric composites with anisotropic constituents. Expressions are provided for the effective moduli tensors of n-phase composites based on the respective strain and electric field concentration tensors. By taking into account the shape and distribution of the inclusion and by invoking a simple numerical procedure, solutions for the electromechanical properties of a general anisotropic inclusion in an anisotropic matrix are obtained. While analytical forms are provided for predicting the electroelastic moduli of composites with spherical and cylindrical inclusions, numerical evaluation of integrals over the composite microstructure is required in order to obtain the corresponding expressions for a general ellipsoidal particle in a piezoelectric matrix. The electroelastic moduli of piezoelectric composites predicted by the analytical model developed in the present study demonstrate excellent agreement with results obtained from three-dimensional finite-element models for several piezoelectric systems that exhibit varying degrees of elastic anisotropy.     

Exact Relations for Effective Tensors of Polycrystals. I. Necessary Conditions

The set of all effective moduli of a polycrystal usually has a nonempty interior. When it does not, we say that there is an exact relation for effective moduli. This can indeed happen as evidenced by recent results [4, 10, 12] on polycrystals. In this paper we describe a general method for finding such relations for effective moduli of laminates. The method is applicable to any physical setting that can be put into the Hilbert space framework developed by Milton[13]. The idea is to use the W-function of Milton[13] that transforms a lamination formula into a convex combination. The method reduces the problem of finding exact relations to a problem from representation theory of SO(d)(d= 2 or 3) corresponding to a particular physical setting. When this last problem is solved, there is a finite amount of calculation required to be done in order to answer the question completely. At present, each candidate relation has to be examined separately in order to confirm the stability under homogenization. We apply our general theory to the settings of conductivity and two‐dimensional elasticity. (Accepted April 4, 1997)     

Dynamic stability analysis of an elastic composite material having a negative-stiffness phase

The rigorous classical bounds of elastic composite materials theory provide limits on the achievable composite stiffnesses in terms of the properties and arrangements of the composite's constituents. These bounds result from the assumption, presumably made for stability reasons, that each constituent material must have positive-definite elastic moduli. If this assumption is relaxed, recently published elasticity analyses and experimental measurements show these bounds can be greatly exceeded, resulting in new materials of enormous potential.The key question is whether a composite material having a non-positive-definite constituent can be stable overall in the practically useful situation of applied traction boundary conditions. Drugan [2007. Elastic composite materials having a negative-stiffness phase can be stable. Phys. Rev. Lett. 98 (5), article no. 055502] first proved the answer is yes, by applying the energy criterion of elastic stability to the basic two- and three-dimensional composites consisting of a cylinder or sphere having non-positive-definite (but strongly elliptic) moduli with a thin positive-definite coating and proving overall stability provided the coating is sufficiently stiff.Here, we perform a complete and direct dynamic stability analysis of the plane strain fundamental elastic composite consisting of a circular cylinder of non-positive-definite material firmly bonded to a positive-definite concentric coating, for the full range of coating thicknesses (i.e., volume fractions). We determine quantitatively the full permissible range of inclusion and coating moduli, as a function of coating thickness, for which the overall composite is stable under dead traction boundary conditions. Among the results, we show that in the thin-coating case, the present dynamic stability analysis leads to precisely the same analytical stability requirements as those derived via the energy criterion by Drugan [2007. Elastic composite materials having a negative-stiffness phase can be stable. Phys. Rev. Lett. 98 (5), article no. 055502], and we derive new analytical stability requirements that are valid for a wider range of coating thickness. At the other extreme, we show that in the case of very thick coatings (corresponding to the dilute case of a matrix-inclusion composite), even an inclusion with merely strongly elliptic moduli can be stabilized by a positive-definite matrix satisfying weak requirements, for which we derive analytical expressions. Overall, our results show that surprisingly weak restrictions on the moduli and thickness of the positive-definite coating are sufficient to stabilize a non-positive-definite inclusion, even one whose moduli are merely strongly elliptic. These results legitimize expanding the search for novel materials with extreme properties to those incorporating a non-positive-definite constituent, and they provide quantitative restrictions on the constituent materials’ moduli and volume fractions, for the geometry examined here, that ensure overall stability of such composite materials.     

Tailored heterogeneity increases overall stability regime of composites having a negative-stiffness inclusion

Recent theoretical and experimental results have shown the possibility of enormous increases in composite material overall elastic stiffness, damping, thermal expansion, piezoelectricity, etc., when the composite contains a tuned non-positive-definite (i.e., negative stiffness) constituent. For such composite materials to have practical utility, they must be stable. Recent research has shown they can be, for a limited range of constituent negative stiffness. This research has treated linear elastic composite materials with homogeneous phases, via the energy method and full dynamic stability analyses.In the present work, we first show how to analyze the composites previously treated by the comprehensive but simpler static stability approach, obtaining closed-form results. We then employ this approach to show that permitting heterogeneity of the positive-definite phase can substantially increase the range of constituent negative stiffness while maintaining overall composite stability. We first treat the positive-definite phase heterogeneity as piecewise homogeneous, and then treat it as continuously-varying. In the continuously-varying heterogeneity case, we seek the radially optimal distribution of the elastic moduli in the coatings, under constant coating average moduli constraint, to permit the most negative possible inclusion stiffness while maintaining overall composite stability. This is accomplished for three coating cases: constant bulk modulus but arbitrarily radially-varying shear modulus; constant shear modulus but arbitrarily radially-varying bulk modulus; and both moduli arbitrarily radially varying. We find the optimal coatings to be: a heterogeneous one with shear modulus being a specific continuously decreasing function of radius for the first case; a homogeneous one for the second case; and a heterogeneous one with both moduli being either Dirac-delta or Heaviside-step decreasing functions of radius for the last case (if the coating moduli are unrestricted in magnitude or have upper limits, respectively). The results show a substantial increase in the permissible inclusion negative stiffness range is provided by coating heterogeneity, while maintaining overall composite stability. Such an increased range of constituent negative stiffness provides an enlarged tuning parameter range for the development of novel, high-performance composite materials.     

Mathematical Models and Methods of the Mechanics of Stochastic Composites

The basic approaches used in mathematical models and general methods for solution of the equations of the mechanics of stochastic composites are generalized. They can be reduced to the stochastic equations of the theory of elasticity of a structurally inhomogeneous medium, to the equations of the theory of effective elastic moduli, to the equations of the theory of elastic mixtures, or to more general equations of the fourth order. The solution of the stochastic equations of the elastic theory for an arbitrary domain involves substantial mathematical difficulties and may be implemented only rather approximately. The construction of the equations of the theory of effective moduli is associated with the problem on the effective moduli of a stochastically inhomogeneous medium, which can be solved by the perturbation method, by the method of moments, or by the method of conditional moments. The latter method is most appropriate. It permits one to determine the effective moduli in a two-point approximation and nonlinear deformation properties. In the structure of equations, the theory of elastic mixtures is more general than the theory of effective moduli; however, since the state equations have not been strictly substantiated and the constants have not been correctly determined, theoretically or experimentally, this theory cannot be used for systematic designing composite structures. A new model of the nonuniform deformation of composites is more promising. It is constructed by performing strict mathematical transformations and averaging the output stochastic equations, all the constants being determined. In the zero approximation, the equations of the theory of effective moduli follow from this model, and, in the first approximation, fourth-order equations, which are more general than those of the theory of mixtures, follow from it     

Mesoscopic natural element analysis of elastic moduli, yield stress and fracture of solids containing a number of voids

A two-dimensional mesoscale simulation method based on the natural element method, which is a kind of meshless method, is developed and applied to the analysis of overall elastic moduli, macro yield stress and void-linking fracture. The calculated results are compared with the theoretical solutions for overall elastic moduli, the experimental macro yield stress for aluminum and brass as well as the improved Gurson’s yield function, and the experimental void linking fracture to show the validity of the proposed method.     

Influence of fiber’s shape and size on overall elastoplastic property for micropolar composites

A general micromechanical method is developed for a micropolar composite with ellipsoidal fibers, where the matrix material is idealized as a micropolar material model. The method is based on a special micro–macro transition method, and the classical effective moduli for micropolar composites can be determined in an analytical way. The influence of both fiber’s shape and size can be analyzed by the proposed method. The effective moduli, initial yield surface and effective nonlinear stress and strain relation for a micropolar composite reinforced by ellipsoidal fibers are examined, it is found that the prediction on the effective moduli and effective nonlinear stress and strain curves are always higher than those based on classical Cauchy material model, especially for the case where the size of fiber approaches to the characteristic length of matrix material. As expected, when the size of fiber is sufficiently large, the classical results (size-independence) can be recovered.     

Rate form of the Eshelby and Hill tensors

Expressions are derived for the rates of change of the S and P tensors for transformed homogeneous inclusions in an anisotropic comparison medium undergoing prescribed changes of its elastic moduli. General results are obtained for ellipsoids and then reduced to yield explicit expressions in terms of the Stroh eigenvalues for cylindrical and disk-shaped inclusions in anisotropic solids and for spherical inclusions in isotropic solids. Applications are illustrated by solving the rate problem for an inhomogeneity in a large volume of a comparison medium, which is shown to be readily adaptable to standard averaging techniques for predictions of rates of change of overall moduli of composite materials experiencing evolution of phase moduli.     

Finite deformations of fibre-reinforced elastic solids with fibre bending stiffness

In the conventional theory of finite deformations of fibre-reinforced elastic solids it is assumed that the strain-energy is an isotropic invariant function of the deformation and a unit vector A that defines the fibre direction and is convected with the material. This leads to a constitutive equation that involves no natural length. To incorporate fibre bending stiffness into a continuum theory, we make the more general assumption that the strain-energy depends on deformation, fibre direction, and the gradients of the fibre direction in the deformed configuration. The resulting extended theory requires, in general, a non-symmetric stress and the couple-stress. The constitutive equations for stress and couple-stress are formulated in a general way, and specialized to the case in which dependence on the fibre direction gradients is restricted to dependence on their directional derivatives in the fibre direction. This is further specialized to the case of plane strain, and finite pure bending of a thick plate is solved as an example. We also formulate and develop the linearized theory in which the stress and couple-stress are linear functions of the first and second spacial derivatives of the displacement. In this case for the symmetric part of the stress we recover the standard equations of transversely isotropic linear elasticity, with five elastic moduli, and find that, in the most general case, a further seven moduli are required to characterize the couple-stress.     

Rigorous bounds on the effective moduli of composites and inhomogeneous bodies with negative-stiffness phases

We review the theoretical bounds on the effective properties of linear elastic inhomogeneous solids (including composite materials) in the presence of constituents having non-positive-definite elastic moduli (so-called negative-stiffness phases). Using arguments of Hill and Koiter, we show that for statically stable bodies the classical displacement-based variational principles for Dirichlet and Neumann boundary problems hold but that the dual variational principle for traction boundary problems does not apply. We illustrate our findings by the example of a coated spherical inclusion whose stability conditions are obtained from the variational principles. We further show that the classical Voigt upper bound on the linear elastic moduli in multi-phase inhomogeneous bodies and composites applies and that it imposes a stability condition: overall stability requires that the effective moduli do not surpass the Voigt upper bound. This particularly implies that, while the geometric constraints among constituents in a composite can stabilize negative-stiffness phases, the stabilization is insufficient to allow for extreme overall static elastic moduli (exceeding those of the constituents). Stronger bounds on the effective elastic moduli of isotropic composites can be obtained from the Hashin–Shtrikman variational inequalities, which are also shown to hold in the presence of negative stiffness.     

The total creep of viscoelastic composites under hydrostatic or antiplane loading

The problem of bounding the total creep (or total stress relaxation) of a composite made of two linear viscoelastic materials and subjected to a constant hydrostatic or antiplane loading is considered. It is done by coupling the immediate and the relaxed responses of the composite, which are pure elastic. The coupled bounds provide the possible range of the total deformation at infinite time as a function of the initial deformation of the composite. For antiplane shear existing bounds for coupled two-dimensional conductivity yield the required coupled bounds, and these are attained by doubly coated cylinder assemblages. The translation method is used to couple the effective bulk moduli of a viscoelastic composite at zero and infinite time. A number of microgeometries are found to attain the bulk modulus bounds. It is shown that the Hashin's composite sphere assemblage does not necessarily correspond to the maximum or minimum overall creep, although it necessarily attains the bounds for effective bulk moduli. For instance, there are cases when the doubly coated sphere microstructure or some special polycrystal arrangements attain the bounds on the total creep.     

On estimation methods for effective moduli of microcracked solids

Summary  Most of the conventional methods for estimating the overall elastic moduli of microcracked solids are defined based on the concept of effective medium or effective field. The formal similarity of these methods is examined in this paper. A one-to-one correspondence relation exists between the effective medium methods and the effective field methods in the sense that they yield identical results. In addition to the conventional estimation techniques, any other number of such approaches may be constructed by appropriately specifying the effective matrix compliance (or stiffness) tensor and the effective stress (or strain) field which a microcrack is assumed to be subjected to. To generate continuous spectra of new methods for estimating the effective elastic moduli, two simple and straightforward approaches are proposed, which contain one or two adjustable parameters in order to yield results of good accuracy. The discussion in this paper can be extended to other kinds of heterogeneous materials. Received 4 October 2000; accepted for publication 30 January 2001     

ANALYTICAL SOLUTIONS FOR ELASTOSTATIC PROBLEMS OF PARTICLE-AND FIBER-REINFORCED COMPOSITES WITH INHOMOGENEOUS INTERPHASES

IntroductionRecently ,theinterphaseincompositeshasattractedtheattentionofmanyresearchers[1- 7].Someinterphasesareproducedbydesignandothersbychemicalreactionsinfabricatingcomposites.Fromtheviewpointofmechanicsofcompositematerials,oneofthefundamentalproblemsofcompositeswithinterphasesistopredicttheeffectivemoduli.Therearemanymethodsforpredictingeffectivemoduliofcomposites,suchasgeneralizedself_consistentmethod(GSCM) ,compositesphereandcompositecylindermodels (CCAandCSA)andIDDestimate ,etc …     

A continuum micromechanical theory of overall plasticity for particulate composites including particle size effect

Classical continuum micromechanics cannot predict the particle size dependence of the overall plasticity for composite materials, a simple analytical micromechanical method is proposed in this paper to investigate this size dependence. The matrix material is idealized as a micropolar continuum, an average equivalent inclusion method is advanced and the Mori–Tanaka's method is extended to a micropolar medium to evaluate the effective elastic modulus tensor. The overall plasticity of composites is predicted by a new secant moduli method based on the second order moment of strain and torsion of the matrix in a framework of micropolar theory. The computed results show that the size dependence is more pronounced when the particle's size approaches to the matrix characteristic length, and for large particle sizes, the prediction coincides with that predicted by classical micromechanical models. The method is analytical in nature, and it can capture the particle size dependence on the overall plastic behavior for particulate composites, and the prediction agrees well with the experimental results presented in literature. The proposed model can be considered as a natural extension of the widely used secant moduli method from a heterogeneous Cauchy medium to a micropolar composite.     

On the isotropic moduli of 2D strain-gradient elasticity

In the present paper, the simplest model of strain-gradient elasticity will be considered, that is, the isotropy in a bidimensional space. Paralleling the definition of the classic elastic moduli, our aim is to introduce second-order isotropic moduli having a mechanical interpretation. A general construction process of these moduli will be proposed. As a result, it appears that many sets can be defined, each of them constituted of 4 moduli: 3 associated with 2 distinct mechanisms and the last one coupling these mechanisms. We hope that these moduli (and the construction process) might be useful for forthcoming investigations on generalized continuum mechanics.     

A numerical simulation for effective elastic moduli of plates with various distributions and sizes of cracks

A fast convergent numerical model is developed to calculate the effective moduli of plates with various distributions and sizes of cracks, in which the crack line is divided into M parts to obtain the unknown traction on the crack line. When M=1, the model reduces to Kachanov's approximation method [Int. J. Solids Struct. 23 (1987) 23]. Six types of crack distributions and three kinds of crack sizes are considered, which are four regular (equilateral triangle, equilateral hexagon, rectangle, and diamond) and two random distributions (random location and orientation, and parallel orientation and random location), and one, two and random crack sizes. Some typical examples are also analyzed using the finite element method (FEM) to validate the present model. Then, the effective moduli associated with the crack distributions and sizes are calculated in detail. The present results for the regular distributions show some very interesting phenomena that have not been revealed before. And for the two random distributions, as the effective moduli depend on samples due to the randomness, the effect of the sample size and number are analyzed first. Then, effective moduli for plates with the three sizes of cracks are calculated. It is found that the effect of crack sizes on the effective moduli is significant for high crack densities, and small for low crack densities, and the random crack size leads to the lowest effective moduli. The present numerical results are compared with several popular micromechanics models to determine which one can provide the optimum estimation of the effective moduli of cracked plates with general crack densities. Furthermore, some existing numerical results are analyzed and discussed.     

Relations Between Relaxation Modulus and Creep Compliance in Anisotropic Linear Viscoelasticity

The results obtained previously for scalar and class P completely monotone relaxation moduli are extended to arbitrary anisotropy. It is shown for general anisotropic viscoelastic media that, if the relaxation modulus is a locally integrable completely monotone function, then the creep compliance is a Bernstein function and conversely. The elastic and equilibrium limits of the two material functions are related to each other. The relaxation modulus or its derivative can be singular at 0. A rigorous general formulation of the relaxation spectrum in an anisotropic viscoelastic medium is given. The effect of Newtonian viscosity on creep compliance is examined. Put some makeup on him and lay him to rest. Anonymous     

Iteration method in linear elasticity of random structure composites containing heterogeneities of noncanonical shape

We consider a linearly elastic composite medium, which consists of a homogeneous matrix containing a statistically inhomogeneous random set of heterogeneities of arbitrary shape. The general integral equations connecting the stress and strain fields in the point being considered with the stress and strain fields in the surrounding points are obtained for the random fields of heterogeneities. The method is based on a recently developed centering procedure where the notion of a perturbator is introduced and statistical averages are obtained without any auxiliary assumptions such as, e.g., effective field hypothesis implicitly exploited in the known centering methods. Effective elastic moduli and the first statistical moments of stresses in the heterogeneities are estimated for statistically homogeneous composites with the general case of both the shape and inhomogeneity of the heterogeneities moduli. The explicit new representations of the effective moduli and stress concentration factors are built by the iteration method in the framework of the quasicristallite approximation but without basic hypotheses of classical micromechanics such as both the EFH and “ellipsoidal symmetry” assumption. Numerical results are obtained for some model statistically homogeneous composites reinforced by aligned identical homogeneous heterogeneities of noncanonical shape. Some new effects are detected that are impossible in the framework of a classical background of micromechanics.