Unified formulae of variational bounds for multiphase anisotropic elastic composites

Using the spherical and deviator decomposition of the polarization and strain tensors, we present a general algorithm for the calculation of variational bounds of dimension d for any type of anisotropic linear elastic composite as a function of the properties of the comparison body. This procedure is applied in order to obtain analytical expressions of bounds for multiphase, linear elastic composites with cubic symmetry where the geometric shapes of the inclusions are arbitrary. For the validation, it can be proved that for the isotropic particular case, the bounds coincide with those recently reported by Gibiansky and Sigmund. On the other hand, based on this general procedure some, classical bounds reported by Hashin for transversely isotropic composites, are reproduced. Numerical calculations and some comparisons with other models and experimental data are shown.     

Effective elastic shear stiffness of a periodic fibrous composite with non-uniform imperfect contact between the matrix and the fibers

In this contribution, effective elastic moduli are obtained by means of the asymptotic homogenization method, for oblique two-phase fibrous periodic composites with non-uniform imperfect contact conditions at the interface. This work is an extension of previous reported results, where only the perfect contact for elastic or piezoelectric composites under imperfect spring model was considered. The constituents of the composites exhibit transversely isotropic properties. A doubly periodic parallelogram array of cylindrical inclusions under longitudinal shear is considered. The behavior of the shear elastic coefficient for different geometry arrays related to the angle of the cell is studied. As validation of the present method, some numerical examples and comparisons with theoretical results verified that the present model is efficient for the analysis of composites with presence of imperfect interface and parallelogram cell. The effect of the non uniform imperfection on the shear effective property is observed. The present method can provide benchmark results for other numerical and approximate methods.     

Variational bounds for anisotropic elastic multiphase composites with different shapes of inclusions

In the present work, unified formulae for the overall elastic bounds for multiphase transversely isotropic composites with different geometrical types of inclusions embedded in a matrix are calculated, including the spherical and long or short continuous cylindrical fiber cases. The influence of the different geometrical configurations of the inclusions on the composites is studied. The transversely isotropic effective bounds are obtained by applying the variational formulation for anisotropic composites developed by Willis, which relies on expressions for the static transversely isotropic Green’s function. Some numerical calculations and comparisons with the effective coefficients derived from the self-consistent approach, asymptotic homogenization method, and finite element method (FEM) are shown for different aspect ratio values, exhibiting good agreement.     

A pattern-based method for bounding the effective response of a nonlinear composite

This paper deals with the prediction of the effective properties of nonlinear composites. Rather than bounding the effective energy, this work aims at bounding directly the effective stress-strain response, by extending a method originally introduced by Milton and Serkov (J. Mech. Phys. Solids 48 (2000) 1295) and recently refined by Talbot and Willis (Proc. Roy. Soc. 460 (2004) 2705). In this paper, bounding the effective response is achieved by introducing a linear comparison composite with the same micro-geometry as the given nonlinear composite, as Ponte Castañeda (J. Mech. Phys. Solids 39 (1991) 45) did for the energy. It is found that any lower bound for the energy of the linear comparison composite generates a corresponding bound for the stress-strain response of the nonlinear composite. A selection of examples is presented to illustrate the method and compare the bounds obtained with existing results.     

Transversely isotropic sequentially laminated composites in finite elasticity

A general expression for the energy-density function of sequentially laminated composites is derived. For the class of neo-Hookean composites in the limit of small deformations well-known results for linear transversely isotropic composites are recovered. However, it is shown that under large deformations these composites are not isotropic. Transversely isotropic composites are obtained with sequentially-coated composites in which the next rank composite is constructed by lamination of the previous composite with thin layers of the matrix phase. The transverse behavior of this sequentially-coated composite is neo-Hookean with shear modulus in the form of the Hashin-Shtrikman bounds for the corresponding class of linear composites. Comparison of the behaviors of these composites with recent estimates for transversely isotropic composites reveals good agreement up to relatively large deformations and volume fractions of the inclusion phase.     

On the effective behavior of nonlinear inelastic composites: I. Incremental variational principles

A new method for determining the overall behavior of composite materials comprised of nonlinear inelastic constituents is presented. Upon use of an implicit time-discretization scheme, the evolution equations describing the constitutive behavior of the phases can be reduced to the minimization of an incremental energy function. This minimization problem is rigorously equivalent to a nonlinear thermoelastic problem with a transformation strain which is a nonuniform field (not even uniform within the phases). In this first part of the study the variational technique of Ponte Castañeda is used to approximate the nonuniform eigenstrains by piecewise uniform eigenstrains and to linearize the nonlinear thermoelastic problem. The resulting problem is amenable to simpler calculations and analytical results for appropriate microstructures can be obtained. The accuracy of the proposed scheme is assessed by comparison of the method with exact results.     

Overall dynamic constitutive relations of layered elastic composites

A method for the homogenization of a layered elastic composite is presented. It allows direct, consistent, and accurate evaluation of the averaged overall frequency-dependent dynamic material constitutive relations without the need for a point-wise solution of the field equations. When the spatial variation of the field variables is restricted by Bloch-form (Floquet-form) periodicity, then these relations together with the overall conservation and kinematical equations accurately yield the displacement or stress mode-shapes and, necessarily, the dispersion relations. The method can also give the point-wise solution of the elastodynamic field equations (to any desired degree of accuracy), which, however, is not required for the calculation of the average overall properties. The resulting overall dynamic constitutive relations are general and need not be restricted by the Bloch-form periodicity.The formulation is based on micromechanical modeling of a representative unit cell of the composite. For waves in periodic layered composites, the overall effective mass-density and compliance (stiffness) are always real-valued whether or not the corresponding unit cell (representative volume element used as a unit cell) is geometrically and/or materially symmetric. The average strain and linear momentum are coupled and the coupling constitutive parameters are always each others' complex conjugates. We separate the overall constitutive relations, which depend only on the composition and structure of the unit cell, from the overall field equations which hold for any elastic composite; i.e., we use only the local field equations and material properties to deduce the overall constitutive relations. Finally, we present solved numerical examples to further clarify the structure of the averaged constitutive relations and to bring out the correspondence of the current method with recently published results.     

Modeling the macroscopic behavior of two-phase nonlinear composites by infinite-rank laminates

A new approach is proposed for estimating the macroscopic behavior of two-phase nonlinear composites with random, particulate microstructures. The central idea is to model composites by sequentially laminated constructions of infinite rank whose macroscopic behavior can be determined exactly. The resulting estimates incorporate microstructural information up to the two-point correlation functions, and require the solution to a Hamilton–Jacobi equation with the inclusion concentration and the macroscopic fields playing the role of ‘time’ and ‘spatial’ variables, respectively. Because they are realizable, by construction, these estimates are guaranteed to be convex, to satisfy all pertinent bounds, to exhibit no duality gap, and to be exact to second order in the heterogeneity contrast. Sample results are provided for two- and three-dimensional power-law composites, and are compared with other homogenization estimates, as well as with numerical simulations available from the literature. The estimates are found to give physically sensible predictions for all the cases considered, even for extreme values of the nonlinearity and heterogeneity contrast. Interestingly, in the case of isotropic porous materials under hydrostatic loadings, the estimates agree exactly with standard Gurson-type models for viscoplastic porous media.     

Neo-Hookean fiber-reinforced composites in finite elasticity

The response of a transversely isotropic fiber-reinforced composite made out of two incompressible neo-Hookean phases undergoing finite deformations is considered. An expression for the effective energy-density function of the composite in terms of the properties of the phases and their spatial distribution is developed. For the out-of-plane shear and extension modes this expression is based on an exact solution for the class of composite cylinder assemblages. To account for the in-plane shear mode we incorporate an exact result that was recently obtained for a special class of transversely isotropic composites. In the limit of small deformation elasticity the expression for the effective behavior agrees with the well-known Hashin-Shtrikman bounds. The predictions of the proposed constitutive model are compared with corresponding numerical simulation of a composite with a hexagonal unit cell. It is demonstrated that the proposed model accurately captures the overall response of the periodic composite under any general loading modes.     

Analytical procedure for determining the linear and nonlinear effective properties of the elastic composite cylinder

In this work we consider a cylindrical structure composed of a nonlinear core (inhomogeneity) surrounded by a different nonlinear shell (matrix). We elaborate a technique for determining its linear elastic moduli (second order elastic constants) and the nonlinear elastic moduli, which are called Landau coefficients (third order elastic constants). Firstly, we develop a nonlinear perturbation method which is able to turn the initial nonlinear elastic problem into a couple of linear problems. Then, we prove that only the solution of the first linear problem is necessary to calculate the linear and nonlinear effective properties of the heterogeneous structure. The following step consists in the exact solution of such a linear problem by means of the complex elastic potentials. As result we obtain the exact closed forms for the linear and nonlinear effective elastic moduli, which are valid for any volume fraction of the core embedded in the external shell.