Size-dependent effective elastic moduli of particulate composites with interfacial displacement and traction discontinuities

This work aims at estimating the size-dependent effective elastic moduli of particulate composites in which both the interfacial displacement and traction discontinuities occur. To this end, the interfacial discontinuity relations derived from the replacement of a thin uniform interphase layer between two dissimilar materials by an imperfect interface are reformulated so as to considerably simplify the characteristic expressions of a general elastic imperfect model which is adopted in the present work and include the widely used Gurtin–Murdoch and spring-layer interface models as particular cases. The elastic fields in an infinite body made of a matrix containing an imperfectly bonded spherical particle and subjected to arbitrary remote uniform strain boundary conditions are then provided in an exact, coordinate-free and compact way. With the aid of these results, the elastic properties of a perfectly bonded spherical particle energetically equivalent to an imperfectly bonded one in an infinite matrix are determined. The estimates for the effective bulk and shear moduli of isotropic particulate composites are finally obtained by using the generalized self-consistent scheme and discussed through numerical examples.     

Effective elastic moduli of composites reinforced by particle or fiber with an inhomogeneous interphase

The solution of the strain energy change of an infinite matrix due to the presence of one spherical particle or cylindrical fiber surrounded by an inhomogeneous interphase is the basis of solving effective elastic moduli of corresponding composites based on various micromechanics models. In order to find out the strain energy change, the composite sphere or cylinder, i.e., the spherical particle or cylindrical fiber together with its interphase, is replaced by an effective homogeneous particle or fiber. Independent governing differential equations for each modulus of the effective particle or fiber are derived by extending the replacement method [J. Mech. Phys. Solids 12 (1964) 199]. As far as the strain energy changes of the infinite matrix subjected to various far-field stress systems are concerned, the present model is simple. Meanwhile, FEM analysis is carried out for a verification, which shows that the model can lead to rather accurate results for most practical interphases. Besides, to check the validity of the model further when the interactions among composite cylinders exist, the two problems of an infinite matrix containing two composite cylinders and the effective moduli of composites with the equilateral triangular distribution of composite cylinders are analyzed using FEM. The FEM results show that the model is still rather accurate, especially for the case of interphase properties varying between those of fiber and matrix. Therefore, composite spheres or cylinders are assumed as the effective homogeneous particles or fibers and simple expressions of the effective moduli of composites containing the composite spheres or cylinders are obtained. Furthermore, the present model is compared with some existing models that are based on very complicated derivations.     

Effective property of piezoelectric composites containing coated nano-elliptical fibers with interfacial debonding

Based on both the spring layer interface model and the Gurtin-Murdoch surface/interface model, the anti-plane shear problem is studied for piezoelectric composites containing coated nano-elliptical fibers with imperfect interfaces. By using the complex function method and the technique of conformal mapping, the exact solutions of the electroelastic fields in fiber, coating, and matrix of piezoelectric nanocomposites are derived under far-field anti-plane mechanical and in-plane electrical loads. Furthermore, the generalized self-consistent method is used to accurately predict the effective electroelastic moduli of the piezoelectric nanocomposites containing coated nano-elliptical fibers with imperfect interfaces. Numerical examples are illustrated to show the effects of the material constants of the imperfect interface layers, the aspect ratio of the fiber section, and the fiber volume fraction on the effective electroelastic moduli of the piezoelectric nanocomposites. The results indicate that the effective electroelastic moduli of the piezoelectric nanocomposites can be significantly reduced by the interfacial debonding, but it can be improved by the surface/interface stresses at the small scale, which provides important theoretical reference for the design and optimization of piezoelectric nanodevices and nanostructures.     

Influence of imperfect elastic contact condition on the antiplane effective properties of piezoelectric fibrous composites

A fiber-reinforced periodic piezoelectric composite, where the constituents exhibit transverse isotropic properties, is considered. The fiber cross-section is circular and the periodicity is the same in two orthogonal directions. Imperfect mechanic contact conditions at the interphase between the matrix and fibers are represented in parametric form. In order to analyze the influence of the imperfect interface effect over the behavior of the composite, the effective axial piezoelectric moduli are obtained by means of the Asymptotic Homogenization Method. Some numerical examples are given.     

非理想界面对三相复合材料应力场的影响

对非理想界面的三相复合材料,提出了计算弹性应力场的微观力学模型,在适当的简化假设下,对带界相的颗粒增强和纤维增强复合材料,得到了应力场的计算公式。以剪切载荷为例给出了数值例子。给出的数值结果表明非理想界面对三相复合材料应力场的影响。     

On two models of arbitrarily curved three-dimensional thin interphases in elasticity

In boundary value problems involving thin interphases, it is often desirable to have a model of an interphase which makes possible to solve for the fields in the adjacent media without having to solve for the fields in it. This is usually achieved in the literature by replacing the interphase by a geometrical surface with appropriately designed “imperfect interface” conditions on it. In the present study, carried out in the setting of elasticity, another option is explored: the geometry of the interphase is left intact, and conditions are devised for the displacements and tractions pertaining to the media adjacent to the interphase and evaluated at both sides of it such that they will simulate the presence of the interphase. Those conditions do not involve the fields within the interphase, yet they depend on its material properties and on those of the adjacent media as well, and make possible to solve for the fields in the adjacent media without having to solve for the fields in the interphase. The formulation is given in a parallel orthogonal curvilinear coordinate system suitable for the modeling arbitrarily curved three-dimensional interphases of constant thickness. Both types of the above described interphase models are tested in the setting of a coated infinite fiber embedded in a matrix which is subjected to an anti-plane shear loading and an in-plane transverse shear loading at infinity, and their predictions are compared with the exact solutions for the fields in the three-phase configuration consisting of the interphase and its adjacent media. The model in which the interphase geometry is left intact is observed to perform generally better than the one in which the interphase is replaced by an interface.     

Differential geometrical method in elastic composite with imperfect interfaces

I.IntroductionWhethertheinterfacesofcompositematerialsareperfectornotwillaffectitsmacromechanicaloreffectivepropertiesimportantly.Butsofar,almostallofthestudiesontheeffectivepropertiesofcompositematerialsarebasedontheassumptionthattheinterfacesareperfectl"2].Infact,thisisnotappropriateforallinterfaces[31.Thusthestudiesonmechanicalpropertyofcompositematerialswithimperfaceintert'acehavebeenconsideredrecentlyinsomeliteratures.Hashin16]hasextendedtheelasticextremumprinciplesofminimumpotentialandm…     

Elastic interaction of partially debonded circular inclusions. II. Application to fibrous composite

A complete analytical solution has been obtained of the elasticity problem for a plane containing periodically distributed, partially debonded circular inclusions, regarded as the representative unit cell model of fibrous composite with interface damage. The displacement solution is written in terms of periodic complex potentials and extends the approach recently developed by Kushch et al. (2010) to the cell type models. By analytical averaging the local strain and stress fields, the exact formulas for the effective transverse elastic moduli have been derived. A series of the test problems have been solved to check an accuracy and numerical efficiency of the method. An effect of interface crack density on the effective elastic moduli of periodic and random structure FRC with interface damage has been evaluated. The developed approach provides a detailed analysis of the progressive debonding phenomenon including the interface cracks cluster formation, overall stiffness reduction and damage-induced anisotropy of the effective elastic moduli of composite.     

无限介质中单个纳米涂层夹杂反平面问题的两个模型

采用零厚度界面模型和界面层模型研究了无限介质中单个纳米涂层圆柱形夹杂的反平面问题,利用复变函数方法获得了两种模型夹杂、涂层和基体内应力场的封闭解析解．研究表明,当界面层模型中的界面相厚度趋于零时,界面层模型可以解析地退化为零厚度界面模型．数值算例分析了界面模量不同取值时应力场的分布和应力的尺度依赖性．本文结果丰富了对纳米夹杂力学行为的认识,并可为直接采用零厚度界面模型有困难的纳米夹杂问题的研究提供有价值的参考．     

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.     

Effective mechanical properties of unidirectional composites in the presence of imperfect interface

In this paper, the equivalent inclusion method is implemented to estimate the effective mechanical properties of unidirectional composites in the presence of an imperfect interface. For this purpose, a representative volume element containing three constituents, a matrix, and interface layer, and a fiber component, is considered. A periodic eigenstrain defined in terms of Fourier series is then employed to homogenize non-dilute multi-phase composites. In order to take into account the interphase imperfection effects on mechanical properties of composites, a stiffness parameter in terms of a matrix and interphase elastic modulus is introduced. Consistency conditions are also modified accordingly in such a way that only the part of the fiber lateral stiffness is to be effective in estimating the equivalent composite mechanical properties. Employing the modified consistency equations together with the energy equivalence relation leads to a set of linear equations that are consequently used to estimate the average values of eigenstrain in non-homogeneous phases. It is shown that for composites with both soft and hard reinforcements, largest stiffness parameter that indicates complete fiber–matrix interfacial debonding causes the same equivalent lateral properties.     

Effective longitudinal shear moduli of periodic fibre-reinforced composites with radially-graded fibres

This paper presents a closed-form expression for the homogenized longitudinal shear moduli of a linear elastic composite material reinforced by long, parallel, radially-graded circular fibres with a periodic arrangement. An imperfect linear elastic fibre-matrix interface is allowed. The asymptotic homogenization method is adopted, and the relevant cell problem is addressed. Periodicity is enforced by resorting to the theory of Weierstrass elliptic functions. The equilibrium equation in the fibre domain is solved in closed form by applying the theory of hypergeometric functions, for new wide classes of grading profiles defined in terms of special functions. The effectiveness of the present analytical procedure is proved by convergence analysis and comparison with finite element solutions. A parametric analysis investigating the influence of microstructural and material features on the effective moduli is presented. The feasibility of mitigating the shear stress concentration in the composite by tuning the fibre grading profile is shown.     

Rigorous Bounds on the Torsional Rigidity of Composite Shafts with Imperfect Interfaces

We derive upper and lower bounds for the torsional rigidity of cylindrical shafts with arbitrary cross-section containing a number of fibers with circular cross-section. Each fiber may have different constituent materials with different radius. At the interfaces between the fibers and the host matrix two kinds of imperfect interfaces are considered: one which models a thin interphase of low shear modulus and one which models a thin interphase of high shear modulus. Both types of interface will be characterized by an interface parameter which measures the stiffness of the interface. The exact expressions for the upper and lower bounds of the composite shaft depend on the constituent shear moduli, the absolute sizes and locations of the fibers, interface parameters, and the cross-sectional shape of the host shaft. Simplified expressions are also deduced for shafts with perfect bonding interfaces and for shafts with circular cross-section. The effects of the imperfect bonding are illustrated for a circular shaft containing a non-centered fiber. We find that when an additional constraint between the constituent properties of the phases is fulfilled for circular shafts, the upper and lower bounds will coincide. In the latter situation, the fibers are neutral inclusions under torsion and the bounds recover the previously known exact torsional rigidity.      

Constitutive equations for elastoplastic composites with imperfect bonding

Closed-form constitutive relations are given for the prediction of the overall response of unidirectional fiber-reinforced composites having constituents that are elastoplastic materials. In these equations the damage mode of imperfect bonding between the fiber and matrix phases is incorporated. The interface decohesion is represented by two parameters that completely determine the degree of adhession at the interfaces in thenormal and tangential directions. Perfect contact, perfectly lubricated contact, and complete debonding are obtained as special cases. In the elastic region, the average stress-strain relations are given in terms of the effective elastic moduli of the damaged composite, all of which are given by closed-form expressions. The derived constitutive equations can be readily implemented for the analysis of metal matrix composites.     

含柔性涂层的颗粒增强复合材料弹性模量估计

采用线弹簧型弱界面模型来模拟柔性涂层,研究柔性涂层对复合材料宏观弹性模量的影响.首先利用Mori-Tanaka方法和弱界面球形夹杂问题的弹性解,获得单夹杂内部的平均应力和平均应变,进而求得具有柔性涂层的复合材料的宏观弹性模量,并研究界面柔度对复合材料弹性模量的影响.     

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.     

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.     

Influence of interfaces on effective properties of nanomaterials with stochastically distributed spherical inclusions

The method of conditional moments is generalized to include evaluation of the effective elastic properties of particulate nanomaterials and to investigate the size effect in those materials. Determining the effective constants necessitates finding a stochastically averaged solution to the fundamental equations of linear elasticity coupled with surface/interface conditions (Gurtin–Murdoch model). To obtain such a solution the system of governing stochastic differential equations is first transformed to an equivalent system of stochastic integral equations. Using statistical averaging, the boundary-value problem is then converted to an infinite system of linear algebraic equations. A two-point approximation is considered and the stress fluctuations within the inclusions are neglected in order to obtain a finite system of algebraic equations in terms of component-average strains. Closed-form expressions are derived for the effective moduli of a composite consisting of a matrix and randomly distributed spherical inhomogeneities. As a numerical example a nanoporous material is investigated assuming a model in which the interface effects influence only the bulk modulus of the material. In that model the resulting shear modulus is the same as for the material without surface effects. Dependence of the bulk moduli on the radius of nanopores and on the pore volume fraction is analyzed. The results are compared to, and discussed in the context of other theoretical predictions.     

A continuum approach to the analysis of the stress field in a fiber reinforced composite with a transverse crack

The stress field in a periodically layered composite with an embedded crack oriented in the normal direction to the layering and subjected to a tensile far-field loading is obtained based on the continuum equations of elasticity. This geometry models the 2D problem of fiber reinforced materials with a transverse crack. The analysis is based on the combination of the representative cell method and the higher-order theory. The representative cell method is employed for the construction of Green’s functions for the displacements jumps along the crack line. The problem of the infinite domain is reduced, in conjunction with the discrete Fourier transform, to a finite domain (representative cell) on which the Born–von Karman type boundary conditions are applied. In the framework of the higher-order theory, the transformed elastic field is determined by a second-order expansion of the displacement vector in terms of local coordinates, in conjunction with the equilibrium equations and these boundary conditions. The accuracy of the proposed approach is verified by a comparison with the analytical solution for a crack embedded in a homogeneous plane.Results show the effects of crack lengths, fiber volume fractions, ratios of fiber to matrix Young’s moduli and matrix Poisson’s ratio on the resulting elastic field at various locations of interest. Comparisons with the predictions obtained from the shear lag theory are presented.     

功能梯度界面相对周期分布纤维增强复合材料反平面剪切的影响

基于复变函数理论和边界配点法,探索了功能梯度界面相在周期均匀分布纤维增强复合材料反平面剪切问题中所起的作用.由于纤维在复合材料基体中的周期分布是均匀的,将其简化成含一功能梯度界面相夹杂的方形单胞.采用分层均匀化方法,将功能梯度界面相离散成K层界面层.当K足够大时,每个界面层可视为匀质材料,同时计算得到的复合材料宏观性能趋于定值.根据单胞内的基体、界面相和夹杂的几何外形特点,分别给出复势函数的级数形式,这些复势函数在各组分的相邻界面应满足连续性条件,在单胞的外边界应满足周期性边界条件和远场加载条件,从而确定复势函数中的待定系数,进而根据平均场理论确定复合材料有效模量.主要探讨了夹杂体积分数、各组分模量、功能梯度界面相的模量渐变形式等因素对纤维增强复合材料性能的影响.结果表明:不管基体模量相对于夹杂模量是大还是小,都有对应的界面相模量渐变形式可使夹杂周围的等效应力集中系数减小;另外还发现仅当夹杂模量较大时,功能梯度界面相模量的变化方式对复合材料有效模量产生不可忽视的影响.