Elastoplastic modeling of circular fiber-reinforced ductile matrix composites considering a finite RVE

A micromechanical elastoplastic damage model considering a finite RVE is proposed to predict the overall elastoplastic damage behavior of circular fiber-reinforced ductile (matrix) composites. The constitutive damage model proposed in our preceding work (Kim and Lee, 2009) considering a finite Eshelby’s tensor (Li et al., 2005, Wang et al., 2005) is extended to accommodate the elastoplastic behavior of the composites. On the basis of the exterior-point Eshelby’s tensor for circular inclusions and the ensemble-averaged effective yield criterion, a micromechanical framework for predicting the effective elastoplastic damage behavior of ductile composites is derived. A series of numerical simulations are carried out to illustrate stress–strain response of the proposed micromechanical framework and to examine the influence of a Weibull parameter on the elastoplastic behavior of the composites. Furthermore, comparisons between the present predictions and experimental data available in the literature are made to further assess the predictive capability of the proposed model.     

An elastoplastic multi-level damage model for ductile matrix composites considering evolutionary weakened interface

An elastoplastic multi-level damage model considering evolutionary weakened interface is developed in this work to predict the effective elastoplastic behavior and multi-level damage evolution in particle reinforced ductile matrix composites (PRDMCs). The elastoplastic multi-level damage model is micromechanically derived on the basis of the ensemble-volume averaging procedure and the first-order effects of eigenstrains. The Eshelby’s tensor for an ellipsoidal inclusion with slightly weakened interface [Qu, J., 1993a. Eshelby tensor for an elastic inclusion with slightly weakened interfaces. Journal of Applied Mechanics 60 (4), 1048–1050; Qu, J., 1993b. The effect of slightly weakened interfaces on the overall elastic properties of composite materials. Mechanics of Materials, 14, 269–281] is adopted to model particles having mildly or severely weakened interface, and a multi-level damage model [Lee, H.K., Pyo, S.H., in press. Multi-level modeling of effective elastic behavior and progressive weakened interface in particulate composites. Composites Science and Technology] in accordance with the Weibull’s probabilistic function is employed to describe the sequential, progressive weakened interface in the composites. Numerical examples corresponding to uniaxial, biaxial and triaxial tension loadings are solved to illustrate the potential of the proposed micromechanical framework. A series of parametric analysis are carried out to investigate the influence of model parameters on the progression of weakened interface in the composites. Furthermore, the present prediction is compared with available experimental data in the literature to verify the proposed elastoplastic multi-level damage model.     

Micromechanics-based elastic damage modeling of particulate composites with weakened interfaces

A micromechanical framework is proposed to predict the effective elastic behavior and weakened interface evolution of particulate composites. The Eshelby’s tensor for an ellipsoidal inclusion with slightly weakened interface [Qu, J., 1993a. Eshelby tensor for an elastic inclusion with slightly weakened interfaces. Journal of Applied Mechanics 60 (4), 1048–1050; Qu, J., 1993b. The effect of slightly weakened interfaces on the overall elastic properties of composite materials. Mechanics of Materials 14, 269–281] is adopted to model spherical particles having imperfect interfaces in the composites and is incorporated into the micromechanical framework. Based on the Eshelby’s micromechanics, the effective elastic moduli of three-phase particulate composites are derived. A damage model is subsequently considered in accordance with the Weibull’s probabilistic function to characterize the varying probability of evolution of weakened interface between the inclusion and the matrix. The proposed micromechanical elastic damage model is applied to the uniaxial, biaxial and triaxial tensile loadings to predict the various stress–strain responses. Comparisons between the present predictions with other numerical and analytical predictions and available experimental data are conducted to assess the potential of the present framework.     

Strain gradient solution for a finite-domain Eshelby-type anti-plane strain inclusion problem

A solution for the finite-domain Eshelby-type inclusion problem of a finite elastic body containing an anti-plane strain inclusion of arbitrary cross-sectional shape prescribed with a uniform eigenstrain and a uniform eigenstrain gradient is derived in a general form using a simplified strain gradient elasticity theory (SSGET). The formulation is facilitated by an extended Betti’s reciprocal theorem and an extended Somigliana’s identity based on the SSGET and suitable for anti-plane strain problems. The disturbed displacement field is obtained in terms of the SSGET-based Green’s function for an infinite anti-plane strain elastic body. The solution reduces to that of the infinite-domain anti-plane strain inclusion problem when the boundary effect is not considered. The problem of a circular cylindrical inclusion embedded concentrically in a finite cylindrical elastic matrix undergoing anti-plane strain deformations is analytically solved by applying the general solution, with the Eshelby tensor and its average over the circular cross section of the inclusion obtained in closed forms. This Eshelby tensor, being dependent on the position, inclusion size, matrix size, and a material length scale parameter, captures the inclusion size and boundary effects, unlike existing ones. It reduces to the classical linear elasticity-based Eshelby tensor for the circular cylindrical inclusion in an infinite matrix if both the strain gradient and boundary effects are suppressed. Numerical results quantitatively show that the inclusion size effect can be quite large when the inclusion is small and that the boundary effect can dominate when the inclusion volume fraction is high. However, the inclusion size effect is diminishing with the increase of the inclusion size, and the boundary effect is vanishing as the inclusion volume fraction becomes sufficiently low.     

Eshelby tensors for an ellipsoidal inclusion in a microstretch material

Eshelby tensors for an ellipsoidal inclusion in a microstretch material are derived in analytical form, involving only one-dimensional integral. As micropolar Eshelby tensor, the microstretch Eshelby tensors are not uniform inside of the ellipsoidal inclusion. However, different from micropolar Eshelby tensor, it is found that when the size of inclusion is large compared to the characteristic length of microstretch material, the microstretch Eshelby tensor cannot be reduced to the corresponding classical one. The reason for this is analyzed in details. It is found that under a pure hydrostatic loading, the bulk modulus of a microstretch material is not the same as the one in the corresponding classical material. A modified bulk modulus for the microstretch material is proposed, the microstretch Eshelby tensor is shown to be reduced to the modified classical Eshelby tensor at large size limit of inclusion. The fully analytical expressions of microstretch Eshelby tensors for a cylindrical inclusion are also derived.     

The Arithmetic Mean Theorem of Eshelby Tensor for a Rotational Symmetrical Inclusion

In 1997, H. Nozaki and M. Taya found numerically that for any regular polygonal inclusion except for a square, both the Eshelby tensor at the center and the average Eshelby tensor over the inclusion domain are equal to the Eshelby tensor for a circular inclusion and independent of the orientation of the inclusion. Then in 2001, these remarkable properties were mathematically justified by Kawashita and Nozaki. In this paper, a more radical property is presented for a rotational symmetrical inclusion: For any N-fold (N is an integer greater than 2 and unequal to 4) rotational symmetrical inclusion, the arithmetic mean of the Eshelby tensors at N rotational symmetrical points in the inclusion is the same as the Eshelby tensor for a circular inclusion and independent of the orientation of the inclusion. It follows that the Eshelby tensor at the center and the average Eshelby tensor over the rotational symmetrical inclusion domain are identical to the Eshelby tensor for a circular inclusion and independent of the orientation of the inclusion as well. This paper shows that although the Eshelby property does not hold for non-ellipsoidal inclusions, the Eshelby tensor for a rotational symmetrical inclusion satisfies the arithmetic mean property. Mathematics Subject Classifications (2000) 73C02.     

An elastoplastic damage model for metal matrix composites considering progressive imperfect interface under transverse loading

An elastoplastic damage model considering progressive imperfect interface is proposed to predict the effective elastoplastic behavior and multi-level damage progression in fiber-reinforced metal matrix composites (FRMMCs) under transverse loading. The modified Eshelby’s tensor for a cylindrical inclusion with slightly weakened interface is adopted to model fibers having mild or severe imperfect interfaces [Lee, H.K., Pyo, S.H., 2009. A 3D-damage model for fiber-reinforced brittle composites with microcracks and imperfect interfaces. J. Eng. Mech. ASCE. doi:10.1061/(ASCE)EM.1943-7889.0000039]. An elastoplastic model is derived micromechanically on the basis of the ensemble-volume averaging procedure and the first-order effects of eigenstrains. A multi-level damage model [Lee, H.K., Pyo, S.H., 2008a. Multi-level modeling of effective elastic behavior and progressive weakened interface in particulate composites. Compos. Sci. Technol. 68, 387–397] in accordance with the Weibull’s probabilistic function is then incorporated into the elastoplastic multi-level damage model to describe the sequential, progressive imperfect interface in the composites. Numerical examples corresponding to uniaxial and biaxial transverse tensile loadings are solved to illustrate the potential of the proposed micromechanical framework. A series of parametric analysis are carried out to investigate the influence of model parameters on the progression of imperfect interface in the composites. Furthermore, a comparison between the present prediction and experimental data in the literature is made to assess the capability of the proposed micromechanical framework.     

Special properties of Eshelby tensor for a regular polygonal inclusion

When studying the regular polygonal inclusion in 1997, Nozaki and Taya discovered numerically some remarkable properties of Eshelby tensor: Eshelby tensor at the center and the averaged Eshelby tensor over the inclusion domain are equal to that of a circular inclusion and independent of the orientation of the inclusion. Then Kawashita and Nozaki justified the properties mathematically. In the present paper, some other properties of a regular polygonal inclusion are discovered. We find that for an N-fold regular polygonal inclusion except for a square, the arithmetic mean of Eshelby tensors at N rotational symmetrical points in the inclusion is also equal to the Eshelby tensor for a circular inclusion and independent of the orientation of the inclusion. Furthermore, in two corollaries, we point out that Eshelby tensor at the center, the averaged Eshelby tensor over the inclusion domain, and the line integral average of Eshelby tensors along any concentric circle of the inclusion are all identical with the arithmetic mean.The project supported by the National Natural Science Foundation of China (10172003 and 10372003) The English text was polished by Keren Wang.     

Strain gradient solution for a finite-domain Eshelby-type plane strain inclusion problem and Eshelby’s tensor for a cylindrical inclusion in a finite elastic matrix

A solution for the finite-domain Eshelby-type inclusion problem of a finite elastic body containing a plane strain inclusion prescribed with a uniform eigenstrain and a uniform eigenstrain gradient is derived in a general form using a simplified strain gradient elasticity theory (SSGET). The formulation is facilitated by an extended Betti’s reciprocal theorem and an extended Somigliana’s identity based on the SSGET and suitable for plane strain problems. The disturbed displacement field is obtained in terms of the SSGET-based Green’s function for an infinite plane strain elastic body, which differs from that in earlier studies using the three-dimensional Green’s function. The solution reduces to that of the infinite-domain inclusion problem when the boundary effect is suppressed. The problem of a cylindrical inclusion embedded concentrically in a finite plane strain cylindrical elastic matrix of an enhanced continuum is analytically solved for the first time by applying the general solution, with the Eshelby tensor and its average over the circular cross section of the inclusion obtained in closed forms. This Eshelby tensor, being dependent on the position, inclusion size, matrix size, and a material length scale parameter, captures the inclusion size and boundary effects, unlike existing ones. It reduces to the classical elasticity-based Eshelby tensor for the cylindrical inclusion in an infinite matrix if both the strain gradient and boundary effects are not considered. Numerical results quantitatively show that the inclusion size effect can be quite large when the inclusion is very small and that the boundary effect can dominate when the inclusion volume fraction is very high. However, the inclusion size effect is diminishing with the increase of the inclusion size, and the boundary effect is vanishing as the inclusion volume fraction becomes sufficiently low.     

Micromechanics-based elastic-damage analysis of laminated composite structures

Based on the micromechanics-based constitutive model, derived in our preceding work [Lee, H.K., Pyo, S.H., 2009. A 3D-damage model for fiber-reinforced brittle composites with microcracks and imperfect interfaces. Journal of Engineering Mechanics-ASCE, in press, doi:10.1061/(ASCE)EM.1943.7889.0000039.], incorporating a multi-level damage model and a continuum damage model, the overall elastic behavior and damage evolution of laminated composite structures are studied in detail. The constitutive model is implemented into the finite element program ABAQUS using a user-subroutine UMAT to solve boundary value problems of the composite structures. The validity of the implemented constitutive model is assured by comparing the predicted stress–strain curves with experimental data available in literature under uniaxial tension with various fiber orientations. The results show that the proposed micromechanics-based constitutive model accurately predict the overall elastic-damage behavior of laminated composite structures having different material compositions, thicknesses and boundary conditions.     

Some general properties of Eshelby’s tensor fields in transport phenomena and anti-plane elasticity

Consider an infinite thermally conductive medium characterized by Fourier’s law, in which a subdomain, called an inclusion, is subjected to a prescribed uniform heat flux-free temperature gradient. The second-order tensor field relating the gradient of the resulting temperature field over the medium to the uniform heat flux-free temperature gradient is referred to as Eshelby’s tensor field for conduction. The present work aims at deriving the general properties of Eshelby’s tensor field for conduction. It is found that: (i) the trace of Eshelby’s tensor field is equal to the characteristic function of the inclusion, independently of the latter’s shape; (ii) the isotropic part of Eshelby’s tensor field over the inclusion of arbitrary shape is identical to Eshelby’s tensor field over a 2D circular or 3D spherical inclusion; (iii) when the medium is made of an isotropic material and when the inclusion has some specific rotational symmetries, the value of the Eshelby’s tensor field evaluated at the inclusion gravity center and the symmetric average of Eshelby’s tensor fields are both equal to Eshelby’s tensor field for a 2D circular or 3D spherical inclusion. These results are then extended, with the help of a linear transformation, to the general case where the medium consists of an anisotropic conductive material. The method elaborated and results obtained by the present work are directly transposable to the physically analogous transport phenomena of electric conduction, dielectrics, magnetism, diffusion and flow in porous media and to the mathematically identical phenomenon of anti-plane elasticity.     

Electroelastic behavior modeling of piezoelectric composite materials containing spatially oriented reinforcements

In this work, a modeling of electroelastic composite materials is proposed. The extension of the heterogeneous inclusion problem of Eshelby for elastic to electroelastic behavior is formulated in terms of four interaction tensors related to Eshelby’s electroelastic tensors. Analytical formulations of interaction tensors are presented for ellipsoidal inclusions. These tensors are basically used to derive the self-consistent model, Mori–Tanaka and dilute approaches. Numerical solutions are based on numerical computations of these tensors for various types of inclusions. Using the obtained results, effective electroelastic moduli of piezoelectric multiphase composites are investigated by an iterative procedure in the context of self-consistent scheme. Generalised Mori–Tanaka’s model and dilute approach are re-formulated and the three models are deeply analysed. Concentration tensors corresponding to each model are presented and relationships of effective coefficients are given. Numerical results of effective electroelastic moduli are presented for various types of piezoelectric inclusions and for various orientations and compared to existing experimental and theoretical ones.     

Plastic constitutive behavior of short-fiber/particle reinforced composites

A cylindrical cell model based on continuum theory for plastic constitutive behavior of short-fiber/particle reinforced composites is proposed. The composite is idealized as uniformly distributed periodic arrays of aligned cells, and each cell consists of a cylindrical inclusion surrounded by a plastically deforming matrix. In the analysis, the non-uniform deformation field of the cell is decomposed into the sum of the first order approximate field and the trial additional deformation field. The precise deformation field are determined based on the minimum strain energy principle. Systematic calculation results are presented for the influence of reinforcement volume fraction and shape on the overall mechanical behavior of the composites. The results are in good agreement with the existing finite element analyses and the experimental results. This paper attempts to stimulate the work to get the analytical constitutive relation of short-fiber/particle reinforced composites.     

Micro–macro analysis of shape memory alloy composites

The aim of the paper is to develop a micro–macro approach for the analysis of the mechanical behavior of composites obtained embedding long fibers of Shape Memory Alloys (SMA) into an elastic matrix. In order to determine the overall constitutive response of the SMA composites, two homogenization techniques are proposed: one is based on the self-consistent method while the other on the analysis of a periodic composite. The overall response of the SMA composites is strongly influenced by the pseudo-elastic and shape memory effects occurring in the SMA material. In particular, it is assumed that the phase transformations in the SMA are governed by the wire temperature and by the average stress tensor acting in the fiber. A possible prestrain of the fibers is taken into account in the model.Numerical applications are developed in order to analyze the thermo-mechanical behavior of the SMA composite. The results obtained by the proposed procedures are compared with the ones determined through a micromechanical analysis of a periodic composite performed using suitable finite elements.Then, in order to study the macromechanical response of structural elements made of SMA composites, a three-dimensional finite element is developed implementing at each Gauss point the overall constitutive laws of the SMA composite obtained by the proposed homogenization procedures. Some numerical applications are developed in order to assess the efficiency of the proposed micro–macro model.     

Strain gradient solution for the Eshelby-type polyhedral inclusion problem

The Eshelby-type problem of an arbitrary-shape polyhedral inclusion embedded in an infinite homogeneous isotropic elastic material is analytically solved using a simplified strain gradient elasticity theory (SSGET) that contains a material length scale parameter. The Eshelby tensor for a polyhedral inclusion of arbitrary shape is obtained in a general analytical form in terms of three potential functions, two of which are the same as the ones involved in the counterpart Eshelby tensor based on classical elasticity. These potential functions, as volume integrals over the polyhedral inclusion, are evaluated by dividing the polyhedral inclusion domain into tetrahedral duplexes, with each duplex and the associated local coordinate system constructed using a procedure similar to that employed by Rodin (1996. J. Mech. Phys. Solids 44, 1977–1995). Each of the three volume integrals is first transformed to a surface integral by applying the divergence theorem, which is then transformed to a contour (line) integral based on Stokes' theorem and using an inverse approach different from those adopted in the existing studies based on classical elasticity. The newly derived SSGET-based Eshelby tensor is separated into a classical part and a gradient part. The former contains Poisson's ratio only, while the latter includes the material length scale parameter additionally, thereby enabling the interpretation of the inclusion size effect. This SSGET-based Eshelby tensor reduces to that based on classical elasticity when the strain gradient effect is not considered. For homogenization applications, the volume average of the new Eshelby tensor over the polyhedral inclusion is also provided in a general form. To illustrate the newly obtained Eshelby tensor and its volume average, three types of polyhedral inclusions – cubic, octahedral and tetrakaidecahedral – are quantitatively studied by directly using the general formulas derived. The numerical results show that the components of the SSGET-based Eshelby tensor for each of the three inclusion shapes vary with both the position and the inclusion size, while their counterparts based on classical elasticity only change with the position. It is found that when the inclusion is small, the contribution of the gradient part is significantly large and should not be neglected. It is also observed that the components of the averaged Eshelby tensor based on the SSGET change with the inclusion size: the smaller the inclusion, the smaller the components. When the inclusion size becomes sufficiently large, these components are seen to approach (from below) the values of their classical elasticity-based counterparts, which are constants independent of the inclusion size.     

Inclusions and inhomogeneities in strain gradient elasticity with couple stresses and related problems

At small length scales and/or in presence of large field gradients, the implicit long wavelength assumption of classical elasticity breaks down. Postulating a form of second gradient elasticity with couple stresses as a suitable phenomenological model for small-scale elastic phenomena, we herein extend Eshelby’s classical formulation for inclusions and inhomogeneities. While the modified size-dependent Eshelby’s tensor and hence the complete elastic state of inclusions containing transformation strains or eigenstrains is explicitly derived, the corresponding inhomogeneity problem leads to integrals equations which do not appear to have closed-form solutions. To that end, Eshelby’s equivalent inclusion method is extended to the present framework in form of a perturbation series that then can be used to approximate the elastic state of inhomogeneities. The approximate scheme for inhomogeneities also serves as the basis for establishing expressions for the effective properties of composites in second gradient elasticity with couple stresses. The present work is expected to find application towards nano-inclusions and certain types of composites in addition to being the basis for subsequent non-linear homogenization schemes.     

Modeling for piezoelectric-shape memory alloy composites

A new concept of a piezoelectric ceramic/shape memory alloy (SMA) composite is proposed with aim of using this as a new actuator material with fast actuation speed and large strain. To prove the new concept, a new model is constructed based on Eshelby formulation where linear piezoelectric constitutive equations and bi-linear superelastic equations of SMA are used. The predictions of the strain induced by applied stress and electric field are made for two simple designs of piezo–SMA composites, 1-D series and 1-D parallel laminated composites. The proposed model indicated that 1-D parallel laminate provides the highest strain induced under bias stress and applied electric field among other composite geometries.     

Micro-mechanical behaviors of SMA composite materials under hygro-thermo-elastic strain fields

An analytical procedure to evaluate the behavior of shape memory alloy (SMA) composite under hygrothermal environment is presented. The SMA wires are considered as inclusions embedded in a homogeneous matrix medium of the composite. The inhomogeneity associated with the phase transformation and thermal strains in the SMA wire as well as the hygrothermal strain in the matrix is homogenized using Eshelby’s equivalent inclusion method. In the present work, a similar approach adopted for SMA composites by Marfia and Sacco [Marfia, S., Sacco, E., 2005. Micromechanics and homogenisation of SMA-wire-reinforced materials. J. Appl. Mech. 72 (2), 259–268.] is considered in order to validate the response of SMA composite subjected to thermo-elastic strain field. However, in the present approach, certain modifications and new derivations for the inelastic strain tensors is carried out. First, the constitutive laws for the SMA wire and matrix are expressed in terms of the average strain in the composite. The evolutionary equations used to characterize the pseudoelastic (PE) behavior of the SMA wire are redefined in terms of the eigen strains (phase transformation and thermal strains) occurring in the SMA wire, which are then expressed in terms of the average strain in the composite. Further, the SMA composite constitutive law under coupled hygro-thermo-elastic strain fields is proposed. The generic homogenized hygric and thermal inelastic composite tensors required for the proposed hygro-thermo-elastic constitutive law are derived. Finally, the SMA composite lamina is characterized using Eshelby’s equivalent inclusion method. Using the proposed modifications and derivations, the analytical results are validated for the case of thermo-elastic strain fields and the procedure is then extended to evaluate the SMA composite behavior under hygro-thermo-elastic strain fields. The results include the effect of thermo-elastic and hygro-thermo-elastic strains on the transformation stresses and the nature of hysteresis due to hygric and thermo-elastic strains.     

SMA短纤维复合材料的热胀系数和相变应变系数

基于Eshelby的等效夹杂模型、Mori和Tanaka的场平均法,考虑到形状记忆合金(SMA)的强物理非线性,发展了增量型的等效夹杂模型(Incremental Equivalent Inclusion Model).讨论了SMA短纤维增强的铝基复合材料的热胀系数和相变应变系数.特别研究了SMA短纤维复合材料纤维几何尺寸和体积分数等参数对SMA复合材料的热胀系数和相变应变系数的影响.这些工作对于指导材料设计和了解SMA复合材料热机械特性是非常有意义的.