Three-dimensional stochastic analysis using a perturbation-based homogenization method for elastic properties of composite material considering microscopic uncertainty

This paper discusses evaluation of influence of microscopic uncertainty on a homogenized macroscopic elastic property of an inhomogeneous material. In order to analyze the influence, the perturbation-based homogenization method is used. A higher order perturbation-based analysis method for investigating stochastic characteristics of a homogenized elastic tensor and an equivalent elastic property of a composite material is formulated.As a numerical example, macroscopic stochastic characteristics such as an expected value or variance, which is caused by microscopic uncertainty in material properties, of a homogenized elastic tensor and homogenized equivalent elastic property of unidirectional fiber reinforced plastic are investigated. The macroscopic stochastic variation caused by microscopic uncertainty in component materials such as Young’s modulus or Poisson’s ratio variation is evaluated using the perturbation-based homogenization method. The numerical results are compared with the results of the Monte-Carlo simulation, validity, effectiveness and a limitation of the perturbation-based homogenization method is investigated. With comparing the results using the first-order perturbation-based method, effectiveness of a higher order perturbation is also investigated.     

Asymptotic solutions by the successive disordering method

We develop the periodic componentmethod [1] and represent the solution of a stochastic boundary value elasticity problem for a random quasiperiodic structure with a given disordering degree of inclusions in the matrix via the deviations from the corresponding solution for a random structure with a smaller disordering degree. An example in which the tensor of elastic properties of a composite is calculated is used to illustrate the asymptotic and differential approaches of the successive disordering method. The asymptotic approach permits representing the quasiperiodic structure with a given chaos coefficient and the desired tensor of effective elastic properties as a result of small successive disordering of an originally ideally periodic structure of a composite with known tensor of elastic properties. In the differential approach, a differential equation is obtained for the tensor of effective elastic properties as a function of the chaos coefficient. Its solution coincides with the solution provided by the asymptotic approach. The solution is generalized to the case of piezoactive composites, and a numerical analysis of the effective properties is performed for a PVF (polyvinylidene fluoride) piezoelectric with various quasiperiodic structures on the basis of the cubic structure with spherical inclusions of a high-module elastic material.     

Stochastic homogenization analysis on elastic properties of fiber reinforced composites using the equivalent inclusion method and perturbation method

This paper describes a methodology for evaluation of influence of microscopic uncertainty in material properties and geometry of a microstructure on a homogenized macroscopic elastic property of an inhomogeneous material. For the analysis of the stochastic characteristics of a homogenized elastic property, the first-order perturbation method is used. In order to analyze the influence of microscopic geometrical uncertainty, the perturbation-based equivalent inclusion method is formulated. In this paper, an analytical form of the perturbation term using the equivalent inclusion method is provided.As a numerical example, macroscopic stochastic characteristics such as an expected value or variance of the homogenized elastic tensor of a unidirectional fiber reinforced plastic, which is caused by microscopic uncertainty in material properties or geometry of a microstructure, are estimated with computing the first order perturbation term of the homogenized elastic tensor. Compared the results of the proposed method with the results of the Monte-Carlo simulation, validity, effectiveness and a limitation of the perturbation-based homogenization method is investigated.     

Sensitivity of the macroscopic elasticity tensor to topological microstructural changes

A remarkably simple analytical expression for the sensitivity of the two-dimensional macroscopic elasticity tensor to topological microstructural changes of the underlying material is proposed. The derivation of the proposed formula relies on the concept of topological derivative, applied within a variational multi-scale constitutive framework where the macroscopic strain and stress at each point of the macroscopic continuum are volume averages of their microscopic counterparts over a representative volume element (RVE) of material associated with that point. The derived sensitivity—a symmetric fourth order tensor field over the RVE domain—measures how the estimated two-dimensional macroscopic elasticity tensor changes when a small circular hole is introduced at the microscale level. This information has potential use in the design and optimisation of microstructures.     

A MULTISCALE MECHANICAL MODEL FOR MATERIALS BASED ON VIRTUAL INTERNAL BOND THEORY

Only two macroscopic parameters are needed to describe the mechanical properties of linear elastic solids, i.e. the Poisson's ratio and Young's modulus. Correspondingly, there should be two microscopic parameters to determine the mechanical properties of material if the macroscopic mechanical properties of linear elastic solids are derived from the microscopic level. Enlightened by this idea, a multiscale mechanical model for material, the virtual multi-dimensional internal bonds (VMIB) model, is proposed by incorporating a shear bond into the virtual internal bond (VIB) model. By this modification, the VMIB model associates the macro mechanical properties of material with the microscopic mechanical properties of discrete structure and the corresponding relationship between micro and macro parameters is derived. The tensor quality of the energy density function, which contains coordinate vector, is mathematically proved. Prom the point of view of VMIB, the macroscopic nonlinear behaviors of material could be attributed to the evolution of virtual bond distribution density induced by the imposed deformation. With this theoretical hypothesis, as an application example, a uniaxial compressive failure of brittle material is simulated. Good agreement between the experimental results and the simulated ones is found.     

基于扭曲Kagome点阵结构的一模材料设计

零能模式超材料指弹性矩阵的特征值中有若干为零的弹性材料,根据零特征值的个数可将其分类为一模至五模材料。当前,针对五模材料已有较深入研究,并在水声和弹性波调控方面获得重要应用,而对其他类型零能模式材料的研究尚未展开。本文对扭曲Kagome周期桁架这样一类欠约束点阵材料的有效弹性性质进行了研究,结果表明通过调节点阵材料的微观几何构型和杆件刚度,该类结构能够涵盖一系列一模材料谱系。针对给定一模弹性张量,发展了软-硬模式分离的微结构逆向优化设计策略。通过特定一模材料中的波传播现象对有效性质预测和微结构设计进行了数值验证。     

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.     

Anisotropy of plane complex elastic bodies

The problem of the search of the invariants of an anisotropic elastic tensor representing the mechanical response of a complex elastic body in a two dimensional space is addressed, in particular for a tensor that does not possess all the tensor symmetries typical of classical elasticity. The invariants of the stiffness tensor are found and all the possible types of orthotropy are discussed.     

PMN-PT夹层弹性半空间结构中SH波的传播Propagation of SH waves in sandwich structures consisting of PMN-PT single crystal layer and elastic half-spaces

分析了弹性上下半空间和PMN‐PT单晶层组成的夹层结构中SH波的传播性质,PMN‐PT单晶沿[011]c方向极化,宏观上呈mm2对称,且晶体沿角度θ方向切割。基于正交各向异性压电材料和各向同性弹性材料的基本方程,得到了夹层结构中SH波传播时行列式形式的频散方程。通过对数值算例进行分析可以看出,PMN‐PT单晶的切割角度和弹性材料属性对结构中的相速度有很大影响,因此波的某些传播性能可以通过材料的设计以及晶体切割的方向来实现,这些结论为声表面波器件的开发和应用提供了理论依据。     

基于虚内键理论的材料多尺度力学模型

宏观上线弹性材料的力学属性只需杨氏模量和泊松比两个相互独立的参量来控制;相应地,微观上也需要两个相互独立的参量来控制.基于这个思想,在原VIB模型中引入了切向键,并提出了VMIB模型.该模型在材料的宏观力学属性与微观虚拟键力学属性之间建立起了一座桥梁.考虑到模型中能量密度函数含有坐标轴方向向量一项,该文对能量密度函数的张量性进行了严格的数学证明,并将VMIB模型初步应用到脆性材料的单轴受压破坏.     

Mindlin second-gradient elastic properties from dilute two-phase Cauchy-elastic composites Part II: Higher-order constitutive properties and application cases

Starting from a Cauchy elastic composite with a dilute suspension of randomly distributed inclusions and characterized at first-order by a certain discrepancy tensor (see part I of the present article), it is shown that the equivalent second-gradient Mindlin elastic solid: (i) is positive definite only when the discrepancy tensor is negative defined; (ii) the non-local material symmetries are the same of the discrepancy tensor, and (iii) the non-local effective behaviour is affected by the shape of the RVE, which does not influence the first-order homogenized response. Furthermore, explicit derivations of non-local parameters from heterogeneous Cauchy elastic composites are obtained in the particular cases of: (a) circular cylindrical and spherical isotropic inclusions embedded in an isotropic matrix, (b) n-polygonal cylindrical voids in an isotropic matrix, and (c) circular cylindrical voids in an orthotropic matrix.     

Optimizing multifunctional materials: Design of microstructures for maximized stiffness and fluid permeability

Topology optimization is used to systematically design periodic materials that are optimized for multiple properties and prescribed symmetries. In particular, mechanical stiffness and fluid transport are considered. The base cell of the periodic material serves as the design domain and the goal is to determine the optimal distribution of material phases within this domain. Effective properties of the material are computed from finite element analyses of the base cell using numerical homogenization techniques. The elasticity and fluid flow inverse homogenization design problems are formulated and existing techniques for overcoming associated numerical instabilities and difficulties are discussed. These modules are then combined and solved to maximize bulk modulus and permeability in periodic materials with cubic elastic and isotropic flow symmetries. The multiphysics problem is formulated such that the final design is dependent on the relative importance, or weights, assigned by the designer to the competing stiffness and flow terms in the objective function. This allows the designer to tailor the microstructure according to the materials’ future application, a feature clearly demonstrated by the presented results. The methodology can be extended to incorporate other material properties of interest as well as the design of composite materials.     

Second-order homogenization of periodic materials based on asymptotic approximation of the strain energy: formulation and validity limits

In this paper a second-order homogenization approach for periodic material is derived from an appropriate representation of the down-scaling that correlates the micro-displacement field to the macro-displacement field and the macro-strain tensors involving unknown perturbation functions. These functions take into account of the effects of the heterogeneities and are obtained by the solution of properly defined recursive cell problems. Moreover, the perturbation functions and therefore the micro-displacement fields result to be sufficiently regular to guarantee the anti-periodicity of the traction on the periodic unit cell. A generalization of the macro-homogeneity condition is obtained through an asymptotic expansion of the mean strain energy at the micro-scale in terms of the microstructural characteristic size ?; the obtained overall elastic moduli result to be not affected by the choice of periodic cell. The coupling between the macro- and micro-stress tensor in the periodic cell is deduced from an application of the generalised macro-homogeneity condition applied to a representative portion of the heterogeneous material (cluster of periodic cell). The correlation between the proposed asymptotic homogenization approach and the computational second-order homogenization methods (which are based on the so called quadratic ansätze) is obtained through an approximation of the macro-displacement field based on a second-order Taylor expansion. The form of the overall elastic moduli obtained through the two homogenization approaches, here proposed, is analyzed and the differences are highlighted. An evaluation of the developed method in comparison with other recently proposed in literature is carried out in the example where a three-phase orthotropic material is considered. The characteristic lengths of the second-order equivalent continuum are obtained by both the asymptotic and the computational procedures here analyzed. The reliability of the proposed approach is evaluated for the case of shear and extensional deformation of the considered two-dimensional infinite elastic medium subjected to periodic body forces; the results from the second-order model are compared with those of the heterogeneous continuum.     

Mindlin second-gradient elastic properties from dilute two-phase Cauchy-elastic composites. Part I: Closed form expression for the effective higher-order constitutive tensor

It is shown that second-order homogenization of a Cauchy-elastic dilute suspension of randomly distributed inclusions yields an equivalent second gradient (Mindlin) elastic material. This result is valid for both plane and three-dimensional problems and extends earlier findings by Bigoni and Drugan [Bigoni, D., Drugan, W.J., 2007. Analytical derivation of Cosserat moduli via homogenization of heterogeneous elastic materials. J. Appl. Mech. 74, 741–753] from several points of view: (i) the result holds for anisotropic phases with spherical or circular ellipsoid of inertia; (ii) the displacement boundary conditions considered in the homogenization procedure is independent of the characteristics of the material; (iii) a perfect energy match is found between heterogeneous and equivalent materials (instead of an optimal bound). The constitutive higher-order tensor defining the equivalent Mindlin solid is given in a surprisingly simple formula. Applications, treatment of material symmetries and positive definiteness of the effective higher-order constitutive tensor are deferred to Part II of the present article.     

Averaging Anisotropic Elastic Constant Data

A method of averaging the data on the anisotropic elastic constants of a material is presented. The anisotropic elastic constants are represented by the elasticity tensor which is expressed as a second rank tensor in a space of six dimensions. The method consists of averaging eigenbases of different measurements of the elasticity tensor, then averaging the eigenvalues referred to the average eigenbasis. The eigenvalues and eigenvectors are obtained by using a representation of the stress-strain relations due, in principle, to Kelvin [17, 18]. The formulas for the representation of the averaged elasticity tensor are simple and concise. The applications of these formulas are illustrated using previously reported data, and are contrasted with the traditional analysis of the same data by Hearmon [9]. An interesting result that emerges from this analysis is a method dealing with variable composition anisotropic elastic materials whose elastic constants depend upon the particular composition. In the case of porous isotropic materials, for example, it is customary to regress the Young's modulus against porosity. The results of this paper suggest a structure or paradigm for extending to anisotropic materials this empirical method of regressing elastic constant data against composition or porosity.     

Décomposition de Kelvin et concept de contraintes effectives multiples pour les matériaux anisotropes

The effective stress concept, now classical in continuum damage mechanics, is generalized to the case of an initial anisotropy. In order to be used for both damage–elasticity and damage–(visco-)plasticity coupling, the effective stress should not depend on the elastic properties. Kelvin decomposition of the elasticity tensor allows to define such a stress for isotropic and cubic symmetries. For other material symmetries, the concept of multiple effective stresses is proposed. To cite this article: R. Desmorat, C. R. Mecanique 337 (2009).     

A finite-deformation theory of plasticity

Based on a multiplicative decomposition of local deformation into elastic and plastic deformations general constitutive equations of elastic-plastic materials are proposed. Two alternative approaches are discussed: one in which the elastic deformation is used as an independent variable, and the other in which the stress is one of the independent variables. The appropriate material symmetries are defined, and it is shown that the plastic spin is absent in the theory of isotropic materials. Analysis of a simple extension is given as an example.     

Effective elastic properties of periodic composite medium

A new method is presented for calculating the bulk effective elastic stiffness tensor of a two-component composite with a periodic microstructure. The basic features of this method are similar to the one introduced by Bergman and Dunn (1992) for the dielectric problem. It is based on a Fourier representation of an integro-differential equation for the displacement field, which is used to produce a continued-fraction expansion for the elastic moduli. The method enabled us to include a much larger number of Fourier components than some previously proposed Fourier methods. Consequently our method provides the possibility of performing reliable calculations of the effective elastic tensor of periodic composites that are neither dilute nor low contrast, and are not restricted to arrays of nonoverlapping inclusions. We present results for a cubic array of nonoverlapping spheres, intended to serve as a test of quality, as well as results for a cubic array of overlapping spheres and a two dimensional hexagonal array of circles (a model for a fiber reinforced material) for comparison with previous work.     

Corotational rates in constitutive modeling of elastic-plastic deformation

The principal axes technique is used to develop a new hypoelastic constitutive model for an isotropic elastic solid in finite deformation. The new model is shown to produce solutions that are independent of the choice of objective stress rate. In addition, the new model is found to be equivalent to the isotropic finite elastic model; this is essential if both models describe the same material.

The new hypoelastic model is combined with an isotropic flow rule to form an elastic-plastic rate constitutive equation. Use of the principal axes technique ensures that the stress tensor is coaxial with the elastic stretch tensor and that solutions do not depend on the choice of objective stress rate. The flow rule of von Mises and a parabolic hardening law are used to provide an example of application of the new theory. A solution is obtained for the prescribed deformation of simple rectilinear shear of an isotropic elastic and isotropic elastic-plastic material.