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

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

Thin interphase/imperfect interface in elasticity with application to coated fiber composites

The imperfect interface conditions which are equivalent to the effect of a thin elastic interphase are derived by a Taylor expansion method in terms of interface displacement and traction jumps. Plane and cylindrical interfaces are analyzed as special cases. The effective elastic moduli of a unidirectional coated fiber composite are obtained on the basis of the derived imperfect interface conditions. High accuracy of the method is demonstrated by comparison of solutions of several problems in terms of the imperfect interface conditions or explicit presence of interphase as a third phase. The problems considered are transverse shear of a coated infinite fiber in infinite matrix and effective transverse bulk and shear moduli and effective axial shear modulus of a coated fiber composite. Unlike previous elastic imperfect interface conditions in the literature, the present ones are valid for the entire range of interphase stiffness, from very small to very large.     

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.     

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.     

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.     

Thermoelastic characteristics of a composite solid with initial stresses

Thermoelastic problem for a composite solid with initial stresses is considered on the basis of the asymptotic homogenization method. The homogenized model is constructed by means of the two-scale asymptotic homogenization techniques. The major result of a present paper is that the effective (homogenized) thermoelastic characteristics of the composite material depend not only on local distributions of all types of material characteristics: local elastic properties, local thermoelastic properties, but also on local initial stresses. Therefore it is shown that for the inhomogeneous (composite) material local initial stresses contribute towards values of the effective characteristics of the material. This kind of interaction is not possible for the homogeneous materials. From the mathematical viewpoint, the asymptotic homogenization procedure is equivalent to the computation of G-limit of the corresponding operator. And the above noted phenomenon is based on the fact that in the considering case the G-limit of a sum is not equal to the sum of G-limits. The developed general homogenized model is illustrated in the particular case of the small initial stresses, which is common for the practical mechanical problems. The explicit formulas for the effective thermoelastic characteristics and numerical results are obtained for a laminated composite solid with the initial stresses.     

The asymptotic homogenization elasticity tensor properties for composites with material discontinuities

The classical asymptotic homogenization approach for linear elastic composites with discontinuous material properties is considered as a starting point. The sharp length scale separation between the fine periodic structure and the whole material formally leads to anisotropic elastic-type balance equations on the coarse scale, where the arising fourth rank operator is to be computed solving single periodic cell problems on the fine scale. After revisiting the derivation of the problem, which here explicitly points out how the discontinuity in the individual constituents’ elastic coefficients translates into stress jump interface conditions for the cell problems, we prove that the gradient of the cell problem solution is minor symmetric and that its cell average is zero. This property holds for perfect interfaces only (i.e., when the elastic displacement is continuous across the composite’s interface) and can be used to assess the accuracy of the computed numerical solutions. These facts are further exploited, together with the individual constituents’ elastic coefficients and the specific form of the cell problems, to prove a theorem that characterizes the fourth rank operator appearing in the coarse-scale elastic-type balance equations as a composite material effective elasticity tensor. We both recover known facts, such as minor and major symmetries and positive definiteness, and establish new facts concerning the Voigt and Reuss bounds. The latter are shown for the first time without assuming any equivalence between coarse and fine-scale energies (Hill’s condition), which, in contrast to the case of representative volume elements, does not identically hold in the context of asymptotic homogenization. We conclude with instructive three-dimensional numerical simulations of a soft elastic matrix with an embedded cubic stiffer inclusion to show the profile of the physically relevant elastic moduli (Young’s and shear moduli) and Poisson’s ratio at increasing (up to 100 %) inclusion’s volume fraction, thus providing a proxy for the design of artificial elastic composites.     

Homogenization methods to approximate the effective response of random fibre-reinforced Composites

In this article a fibre-reinforced composite material is modelled via an approach employing a representative volume element with periodic boundary conditions. The effective elastic moduli of the material are thus derived. In particular, the method of asymptotic homogenization is used where a finite number of fibres are randomly distributed within the representative periodic cell. The study focuses on the efficacy of such an approach in representing a macroscopically random (hence transversely isotropic) material. Of particular importance is the sensitivity of the method to cell shape, and how this choice affects the resulting (configurationally averaged) elastic moduli. The averaging method is shown to yield results that lie within the Hashin–Shtrikman variational bounds for fibre-reinforced media and compares well with the multiple scattering and (classical) self-consistent approximations with a deviation from the latter in the larger volume fraction cases. Results also compare favourably with well-known experimental data from the literature.     

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.     

The effect of interfacial properties on the cohesion of highly filled composite materials

It has long been recognized that the cohesion of composite materials, in low confinement, is strongly affected by the properties of the interfacial transition zone (ITZ) between inclusions and matrix. While the effect of the ITZ on the elasticity properties of composites has been studied by many authors in the context of linear homogenization methods, the upscaling of the cohesion strength of highly filled composite materials has not been addressed. This is the focus of the non-linear homogenization procedure developed in this paper, which is based on the separation of the heterogeneous material system in phases of constant strength properties, a non-linear elastic representation of the limit stress state in each phase, and the definition of appropriate effective strain quantities that capture the morphological features of the microstructure. Applied to a three phase composite model composed of rigid inclusion, interface zone and matrix, the model provides a quantitative means of studying the effect of the interface cohesion and the interface volume fraction on the composite cohesion. In particular, we identify a critical interface-to-matrix cohesion ratio, below which the composite cohesion is smaller than the one of the matrix. Furthermore, the model lends itself readily to the study of the degradation of the interfacial properties in composite materials. This is shown for non-degraded and chemically softened cement-based materials, for which we provide conclusive evidence (1) that the interface strength properties of mortar are far more affected by chemical degradation than the one of the cement paste matrix; and (2) that chemical degradation does affect the mechanical strength performance of the cement paste not only through a change of volume proportions (i.e. increase of porosity), but as well through a pure chemical softening of the solid’s cohesion.     

On the stiffness of materials containing a disordered array of microscopic holes or hard inclusions

We study the macroscopic mechanical behavior of materials with microscopic holes or hard inclusions. Specifically, we deal with the effective elastic moduli of composites whose microgeometry consists of either soft or hard isolated inclusions surrounded by an elastic matrix. We approach this problem by taking the stiffness of the inclusion phase to be a complex variable, which we eventually evaluate at the soft or hard limits. Our main result states that there is a certain class of non-physical, negative-definite values of the elastic moduli of the inclusion phase for which the effective tensor does not have infinities or become otherwise singular.We present applications of this result to the estimation of effective moduli and to homogenization theorems. The first application involves using complexanalytic methods to obtain rigorous and accurate bounds on the effective moduli of the high-contrast composites under consideration. We also discuss the variational estimates of Rubenfeld & Keller, which yield a complementary set of bounds on these moduli. The best bounds are given by a combination of the analytical and variational results. As a second application, we show that certain known theorems of homogenization for materials with holes are simple consequences of our main result, and in this connection we establish corresponding new theorems for materials with hard inclusions. While our rederivation of the homogenization theorems for materials with holes can be closely related to other known constructions, it appears that certain elements provided by our main result are essential in the proof of homogenization for the hard-inclusion case.     

Multiscale homogenization of n-component composites with semi-elliptical random interface defects

Effective elastic characteristics of periodic multicomponent composite materials with random interface defects are studied in the paper. The defects are assumed to be semi-elliptical and lying with major semi axes along the interfaces, where minor and major semi-axes as well as the defects number are given as input random variables. The homogenization approach has a multiscale character—some algebraic approximation is used first to calculate effective elastic parameters of the interphase including all defects located at the same interface. Equations for interphase random elastic parameters are obtained using MAPLE symbolic mathematics in conjunction with probabilistic generalized perturbation method. A different homogenization method is applied at the micro scale, where the cell problem is solved numerically using the Finite Element Method (FEM) program. Since the composites considered exhibit random variations of both elastic properties and the interface defects, the overall homogenized characteristics must be obtained as random quantities, which is realized on the micro scale by the Monte-Carlo simulation. The proposed interface defects model obeys the porosity effects resulting from the nature of some matrices in engineering composites as well as the interface cracks appearing as a result of composites ageing during static or fatigue fracture.     

基于迭代法的非线性弹性均质化研究

微观结构对复合材料的宏观力学性能具有至关重要的影响, 通过合理设计复合材料微观结构可以得到期望的宏观性能. 均质化方法作为一种有效的设计方法, 它从微观结构的角度出发, 利用均匀化的概念, 实现了对复合材料宏观力学性能的预测和设计. 而当考虑非线性因素, 均质化的实现就非常困难. 本文利用双渐近展开方法, 将位移按照宏观位移和微观位移展开, 推导了非线性弹性均质化方程. 通过直接迭代法, 对非线性弹性均质化方程进行了求解, 并给出了具体的迭代方法和实现步骤. 本文基于迭代步骤和非线性弹性均质化方程编写MATLAB 程序, 对3种典型本构关系的周期性多孔材料平面问题进行了计算, 对比细致模型的应变能、最大位移和等效泊松比, 对程序及迭代方法的准确性进行了验证. 之后对一种三元橡胶基复合材料进行多尺度均质化, 将其分为芯丝尺度和层间尺度. 用线弹性的均质化方法得到了芯丝尺度的等效弹性参数, 并将其作为层间尺度的材料参数. 在层间尺度应用非线性弹性均质化方法对结构进行计算, 得到材料的宏观等效性能, 并以实验结果为基准进行评价.      

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.     

Differential geometrical method in elastic composite with imperfect interfaces

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

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.     

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.     

HOMOGENIZATION—BASED TOPOLOGY DESIGN FOR PURE TORSION OF COMPOSITE SHAFTS

In conjunction with the homogenization theory and the finite element method, the mathematical models for designing the corss-section of composite shafts by maximizing the torsion rigidity are developed in this paper. To obtain the extremal torsion rigidity, both the cross-section of the macro scale shaft and the representative microstructure of the composite material are optimized using the new models. The micro scale computational model addresses the problem of finding the periodic microstructures with extreme shear moduli. The optimal microstructure obtained with the new model and the homogenization method can be used to improve and optimize natural or artificial materials. In order to be more practical for engineering applications, cellular materials rather than ranked materials are used in the optimal process in the existence of optimal bounds for the elastic properties. Moreover, the macro scale model is proposed to optimize the cross-section of the torsional shaft based on the tailared composites. The validating optimal results show that the models are very effective in obtaining composites with extreme elastic properties, and the cross-section of the composite shaft with the extremal torsion rigidity. The project supported by the National Natural Science Foundation of China (10172078 and 10102018)     

Homogenization of two-phase elasto-plastic composite materials and structures: Study of tangent operators,cyclic plasticity and numerical algorithms

We develop homogenization schemes and numerical algorithms for two-phase elasto-plastic composite materials and structures. A Hill-type incremental formulation enables the simulation of unloading and cyclic loadings. It also allows to handle any rate-independent model for each phase. We study the crucial issue of tangent operators: elasto-plastic (or “continuum”) versus algorithmic (or “consistent”), and anisotropic versus isotropic. We apply two methods of extraction of isotropic tangent moduli. We compare mathematically the stiffnesses of various tangent operators. All rate equations are discretized in time using implicit integration. We implemented two homogenization schemes: Mori–Tanaka and a double inclusion model, and two plasticity models: classical J₂ plasticity and Chaboche’s model with non-linear kinematic and isotropic hardenings. We consider composites with different properties and present several discriminating numerical simulations. In many cases, the results are validated against finite element (FE) or experimental data. We integrated our homogenization code into the FE program ABAQUS using a user material interface UMAT. A two-scale procedure allows to compute realistic structures made of non-linear composite materials within reasonable CPU time and memory usage; examples are shown.     

Interfacial debonding of two-phase composite with periodic structure

A mixed analytical-numerical (boundary element method) procedure is presented for estimating the effective elastic moduli of a two-phase periodic composite by application of a unit cell. The two-phase composite consists of a metal/polymer matrix and one/three circular ceramic inclusions with adhesive and partial debonding of the interface. The results are displayed numerically with special attention given to development of plastic zones as debonding occurs. Dependence of load-time history is exhibited.