The effective elastic moduli of columnar composites made of cylindrically anisotropic phases with rough interfaces

The present work aims to determine the effective elastic moduli of a composite having a columnar microstructure and made of two cylindrically anisotropic phases perfectly bonded at their interface oscillating quickly and periodically along the circular circumferential direction. To achieve this objective, a two-scale homogenization method is elaborated. First, the micro-to-meso upscaling is carried out by applying an asymptotic analysis, and the zone in which the interface oscillates is correspondingly homogenized as an equivalent interphase whose elastic properties are analytically and exactly determined. Second, the meso-to-macro upscaling is accomplished by using the composite cylinder assemblage model, and closed-form solutions are derived for the effective elastic moduli of the composite. Two important cases in which rough interfaces exhibit comb and saw-tooth profiles are studied in detail. The analytical results given by the two-scale homogenization procedure are shown to agree well with the numerical ones provided by the finite element method and to verify the universal relations existing between the effective elastic moduli of a two-phase columnar composite.     

Size-dependent effective thermoelastic properties of nanocomposites with spherically anisotropic phases

Composites made of semi-crystalline polymers and nanoparticles have a spherulitic microstructure which can be reasonably represented by a spherically anisotropic volume element. Due to the high surface-to-volume ratio of a nanoparticle, the particle-matrix interface stress, usually neglected in determining the effective elastic moduli of particle-reinforced composites, may have a non-negligible effect. To account for the latter in estimating the effective thermoelastic properties of a composite consisting of nanoparticles embedded in a semi-crystalline polymeric matrix, this work adopts a coherent interface model for the nanoparticle-matrix interface and proposes an extended version of the classical generalized-self consistent method. In particular, Eshelby's formulae widely used to calculate the elastic energy change of a homogeneous medium due to the introduction of an inhomogeneity are extended to the thermoelastic case. The nanoparticle size effect on the effective thermoelastic moduli of the composite are theoretically shown and numerically illustrated.     

Stress singularities near the apex of a cylindrically polarized piezoelectric wedge

Summary In this paper, the stress singularities for a cylindrically polarized piezoelectric wedge are investigated. The characteristic equations are derived analytically by using the extended Lekhnitskii formulation. The piezoelectric material (PZT-4) is polarized in the radial, circular or axial direction, respectively. Similar to the rectilinearly polarized piezoelectric problem, the inplane and antiplane stress fields are uncoupled. The results show the variations of the singularity orders with the changes of the wedge angle, material constants, polarized direction, and the boundary conditions.     

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.     

Interaction of an anti-plane singularity with interfacial anti-cracks in cylindrically anisotropic composites

Anti-plane problem for a singularity interacting with interfacial anti-cracks (rigid lines) under uniform shear stress at infinity in cylindrically anisotropic composites is investigated by utilizing a complex potential technique in this paper. After obtaining the general solution for this problem, the closed solution for the interface containing one anti-crack is presented analytically. In addition, the complex potentials for a screw dislocation dipole inside matrix are obtained by the superimposing method. Expressions of stress singularities around the anti-crack tips, image forces and torques acting on the dislocation or the center of dipole are given explicitly. The results indicate that the anisotropy properties of materials may weaken the stress singularity near the anti-crack tip for the singularity being a concentrated force but enhance the one for the singularity being a screw dislocation and change the equilibrium position of screw dislocation. The presented solutions are valid for anisotropic, orthotropic or isotropic composites and can be reduced to some new or previously known results.     

Asymptotic determination of effective elastic properties of composite materials with fibrous square-shaped inclusions

We propose an asymptotic approach for evaluating effective elastic properties of two-components periodic composite materials with fibrous inclusions. We start with a nontrivial expansion of the input elastic boundary value problem by ratios of elastic constants. This allows to simplify the governing equations to forms analogous to the transport problem. Then we apply an asymptotic homogenization method, coming from the original problem on a multi-connected domain to a so called cell problem, defined on a characterizing unit cell of the composite. If the inclusions' volume fraction tends to zero, the cell problem is solved by means of a boundary perturbation approach. When on the contrary the inclusions tend to touch each other we use an asymptotic expansion by non-dimensional distance between two neighbouring inclusions. Finally, the obtained “limiting” solutions are matched via two-point Padé approximants. As the results, we derive uniform analytical representations for effective elastic properties. Also local distributions of physical fields may be calculated. In some partial cases the proposed approach gives a possibility to establish a direct analogy between evaluations of effective elastic moduli and transport coefficients. As illustrative examples we consider transversally-orthotropic composite materials with fibres of square cross section and with square checkerboard structure. The obtained results are in good agreement with data of other authors.     

Size-dependent effective elastic constants of solids containing nano-inhomogeneities with interface stress

The fundamental framework of micromechanical procedure is generalized to take into account the surface/interface stress effect at the nano-scale. This framework is applied to the derivation of the effective moduli of solids containing nano-inhomogeneities in conjunction with the composite spheres assemblage model, the Mori-Tanaka method and the generalized self-consistent method. Closed-form expressions are given for the bulk and shear moduli, which are shown to be functions of the interface properties and the size of the inhomogeneities. The dependence of the elastic moduli on the size of the inhomogeneities highlights the importance of the surface/interface in analysing the deformation of nano-scale structures. The present results are applicable to analysis of the properties of nano-composites and foam structures.     

On the overall elastic moduli of composites with spherical coated fillers

A micromechanical approach is presented to estimate the overall linear elastic moduli of three phase composites consisting of two phase coated spherical particles randomly dispersed in a homogeneous isotropic matrix. The theoretical method is based on Eshelby’s equivalent inclusion method and its recent extension by Shodja and Sarvestani [J. Appl. Mech. 68 (2001) 3] to evaluate the local field variables in case of double (multi) inhomogeneities. Using Tanaka–Mori theorem [J. Elasticity 2 (1972) 199] and a decomposition of Green’s function integral equation, the pair-wise average phase values of stress and strain in two interacting coated particles are estimated. Following Ju and Chen [Acta Mech. 103 (1994) 103; Acta Mech. 103 (1994) 123] the ensemble phase volume average of stress and strain fields can be evaluated within a representative volume element containing a finite number of coated particles. Comparisons with classical bounds are presented to illustrate the accuracy of the proposed method.     

Stochastic analysis of a thermoelastic problem in functionally graded plates with uncertain material properties

The statistics (i.e., mean and variance) of temperature and thermal stress are analytically obtained in functionally graded material (FGM) plates with uncertainties in the thermal conductivity and coefficient of linear thermal expansion. These FGM plates are assumed to have arbitrary nonhomogeneous thermal and mechanical properties through the entire thickness of plate and are subjected to deterministic convective heating. The stochastic temperature and thermal stress fields are analysed by assuming the FGM plate is multilayered with distinct, random thermal conductivity and coefficient of linear thermal expansion in each layer. Vodicka’s method, which is a type of integral transform method, and a perturbation method are employed to obtain the analytical solutions for the statistics. The autocorrelation coefficients of each random property and cross-correlation coefficients between different random properties are expressed in exponential function forms as a non-homogeneous Markov random field of discrete space. Numerical calculations are performed for FGM plates composed of partially stabilised zirconia (PSZ) and austenitic stainless steel (SUS304), which have the largest dispersion of the random properties at the place where the volume fractions of the two constituent materials are both 0.5. The effects of the spatial change in material composition, thermal boundary condition and correlation coefficients on the standard deviations of the temperature and thermal stress are discussed.     

An explicit expression of the effective moduli for composite materials filled with coated inclusions

The obvious shortcoming of the generalized self-consistent method (GSCM) is that the effective shear modulus of composite materials estimated by the method can not be expressed in an explicit form. This is inconvenient in engineering applications. In order to overcome that shortcoming of GSCM, a reformation of GSCM is made and a new micromechanical scheme is suggested in this paper. By means of this new scheme, both the effective bulk and shear moduli of an inclusion-matrix composite material can be obtained and be expressed in simple explicit forms. A comparison with the existing models and the rigorous Hashin-Shtrikman bounds demonstrates that the present scheme is accurate. By a two-step homogenization technique from the present new scheme, the effective moduli of the composite materials with coated spherical inclusions are obtained and can also be expressed in an explicit form. The comparison with the existing theoretical and experimental results shows that the present solutions are satisfactory. Moreover, a quantitative comparison of GSCM and the Mori-Tanaka method (MTM) is made based on a unified scheme. The project supported by the National Natural Science Foundation of China under the Contract NO. 19632030 and 19572008, and China Postdoctoral Science Foundation     

Retardation of a crack with connections between the faces using an induced thermoelastic stress field

This paper considers local temperature variations near the tip of a crack in the presence of regions in which the crack faces interact. It is assumed that these regions are adjacent to the crack tip and are comparable in size to the crack size. The problem of local temperature variations consists of delay or retardation of crack growth. For a crack with connections between the crack faces subjected to external tensile loads, an induced thermoelastic stress field, and the stresses at the connections preventing crack opening, the boundary-value problem of the equilibrium of the crack reduces to a system of nonlinear singular integrodifferential equations with a Cauchy kernel. The normal and tangential stresses at the connections are found by solving this system of equations. The stress intensity factors are calculated. The energy characteristics of cracks with tip regions are considered. The limiting equilibrium condition for cracks with tip regions is formulated using the criterion of limiting stretching of the connections.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 1, pp. 133–143, January–February, 2005     

Prediction of the effective elastic moduli of materials with irregularly-shaped pores based on the pore projected areas

Statistical modeling is used to correlate geometric parameters of pores with their contributions to the overall Young’s moduli of linearly elastic solids. The statistical model is based on individual pore contribution parameters evaluated by finite element simulations for a small pore subset selected using the design of experiments approach, so there is no need to solve the elasticity problem for all pores in the material. A polynomial relating pore geometric parameters to the contribution parameters is then fitted to the results of the simulations. We found a good correlation between normalized projected areas of the pores on three coordinate planes and their contributions to the corresponding effective Young’s moduli. The model is applied and validated for two large sets of pore geometries obtained by X-ray microcomputed tomography of a carbon/carbon and a 3D woven carbon/epoxy composite specimens.     

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.     

The influence of normal stress anisotropy in predicting scalar dispersion with the v–f model

A numerical study of scalar dispersion is presented to investigate the effectiveness of pairing the v²–f turbulence model with algebraic models for the scalar flux. This approach is contrasted with utilizing a full Second Moment Closure (SMC) as the flow field input to the scalar model. Predictions of scalar transport in a turbulent channel and over a wavy wall are compared to available DNS databases. The latter case includes a scalar release from a point source and therefore detailed comparisons of the three-component turbulent scalar flux are reported. It is found that the transported variable v², representing the near wall turbulent velocity fluctuation scale, can be used to increase the level of normal stress anisotropy provided to algebraic scalar models and thereby improve mean scalar prediction over that of the Standard Gradient Diffusion Hypothesis (SGDH). Improvement is most significant in the near wall region. Three specifications of the normal stresses, derived from v², are considered to provide the link from the v²–f model to the algebraic flux models used to close the scalar transport equation. Barycentric maps are used to examine the state of turbulence anisotropy in each case. As the anisotropy in the normal stress specification becomes more accurate, improvements are realized in the prediction of the spanwise flux as well as the mean concentration.