# Multipoint Approximation of Statistical Descriptors of Local Strain and Stress Fields in Heterogeneous Media Using Integral Equation Method

Authors: Mikhail Tashkinov
Journal: Advances in Mathematical Physics
 (2018)

---

## Abstract

This paper is devoted to derivation of analytic expressions for statistical descriptors of stress and strain fields in heterogeneous media. Multipoint approximations of solutions of stochastic elastic boundary value problems for representative volume elements are investigated. The stress and strain fields are represented in the form of random coordinate functions, for which analytical expressions for the first- and second-order statistical central moments are obtained. Such moments characterize distribution of fields under prescribed loading of a representative volume element and depend on the geometry features and location of components within a volume. The information of the internal geometrical structure is taken into account by means of multipoint correlation functions. Within the framework of the second approximation of the boundary value problem, the correlation functions up to the fifth order are required to calculate the statistical characteristics. Using the method of Green’s functions, analytical expressions for the moments in distinct phases of the microstructure are obtained explicitly in a form of integral equations. Their analysis and comparison with previously obtained results are performed.


---

## Body

## 1. Introduction

Heterogeneous materials nowadays are the essential link of advanced engineering solutions. The advantages of such materials are underpinned by an ability to obtain tailored properties by combining constituents and their characteristics on microscale level. In the context of solving problems connected with evaluation of microstructural behavior of such heterogeneous materials and media, it is necessary to take into account mutual influence of inhomogeneities. As a rule, modeling of the effective response of materials is performed within the concept of representative volume element (RVE). In these frameworks, wide range of mechanical approaches can be used to establish connection between different scales of materials with complex microstructure. Randomly reinforced media can be distinguished in a separate class of materials where local mechanical characteristics can be estimated using statistical instruments and the theory of random functions. Thus, moments of the stress and strain fields in components can be used as the statistical measures of local mechanical behavior [1]. The morphological properties of the microstructure can be formalized using correlation functions that describe interaction of the inhomogeneities with each other and can sense geometrical effects, such as shape, distribution, clustering, and percolation [2].Statistical assessment of mechanical behavior can be done analytically for each phase of media. Kroener [3] and Beran [4] developed statistical mathematical formulations in order to connect geometrical correlation functions and properties of heterogeneous materials. Explicit relations for second moments of stresses were obtained in [5] and utilized second- and third-order interactions between inclusions. Several exact solutions for second-order moment of stress fields were offered for some cases of deterministic structures. One of such models is media with regular structure [6–8]. The method of integral equations can be implemented for obtaining analytical expressions for local statistics. It presumes that mechanical properties of microstructural components are defined with conventional phenomenological equations and criteria while microscopic strain and stress fields are computed using the solutions of stochastic boundary value problems (SBVPs) with rapidly oscillating coefficients. The numerous ways of resolving and closing the integral equations were offered [9–21].The commonly applied simplification when constructing statistical metrics for heterogeneous RVEs is that only lower order interaction between the microstructural elements is being considered. This work contains explicit derivation of the expressions for local statistical descriptors of stress and strain in heterogeneous media using multipoint high order approximation of SBVP solution and second derivative expansion of Green’s function approach.One of the ground concepts used in statistical analysis of random structures and also applied in this work is based on perturbation approach [4, 16, 22, 23] and presumes that local stress and strain fields as well as local stiffness tensor depend on position of the radius vector and can be decomposed into an averaged value and random addition:(1)σijr→=σijr→+σ′ijr→εijr→=εijr→+ε′ijr→Cijklr→=Cijklr→+C′ijklr→,where σ′ijr→=ε′ijr→=C′ijklr→=0, r→ is radius vector with components (x1,x2,x3). The averaging procedure for any coordinate-dependent value fr→ is determined as(2)fr→=1V∫Vfr→dV,where integration is performed over the representative volume V.


## 2. Approximation of Stochastic Boundary Value Problem Solution

Stochastic boundary value problem (SBVP) of elasticity theory for RVEs of heterogeneous materials with random microstructural elements can be defined with the following equations:(3)σij,jr→=0(4)εijr→=12ui,jr→+uj,ir→(5)σijr→=Cijklr→εklr→(6)uir→r→∈Гu=eijrj,where Cijklr→ is stiffness tensor, eij are components of tensor of small strains which define applied load, and rj are coordinates of the surface points. In this formulation it is assumed that conditions of ideal contact are satisfied on the interphase surface. The boundary value problem is statistically nonlinear since its physical equations (5) contain multiplication of random fields.Due to ergodicity hypothesis volume averaging is equal to statistical averaging:εij∗=εij(r→, σij∗=σij(r→.For two-phase heterogeneous media, containing matrix and one type of inclusions, the stiffness tensorCijklr→ is defined as(7)Cijklr→=λr→CijklI+1-λr→CijklM,where CijklI and CijklM are constant stiffness tensors of inclusions and matrix, respectively. The indicator function λ(r→) is defined as(8)λr→=1,r→∈VI0,r→∈VM,where VI is spatial volume of inclusions and VM is volume of matrix.After averaging (7), the constant tensor, depending on volume fraction, is obtained:(9)Cijklr→=pCijklI+1-pCijklM,where p=λ(r→ is volume fraction of inclusions.Going back to solution of the boundary value problem, (3)–(6) can be transformed into the following representation:(10)Cijklr→12uk,lr→+ul,kr→,j=0.Due to symmetrical properties of the stiffness tensorCijkl(r→)=Cijlk(r→ it can be reduced to(11)Cijklr→uk,lr→,j=0.Applying the perturbation approach, (11) can be expressed through mean and fluctuation:(12)Cijkl+C′ijklr→ukr→,l+u′k,lr→,j=0.Here, averaged valuesCijkl and uk(r→ are expressed via (9) and boundary conditions (6).Expanding (12) gives the following:(13)Cijklu′k,lr→,j=-C′ijklr→ekl+C′ijklr→u′k,lr→,j.Denoting the right part asPij,j(r→):(14)Cijklu′k,ljr→=-Pij,jr→.Solution of this differential equation in displacementsu′k(r→) can be found using Green’s function method. According to definition of Green’s functions [24], for the constant tensor Cijkl Green’s function Gknr→,r→1 and its derivatives are equal 0 in infinity and satisfy(15)CijklGkn,jlr→,r→1=-δinδr→-r→1.Here,δ(r→-r→1 is Dirac’s delta function and δin is Kronecker symbol.Using Green’s function, fluctuations of displacementsu′i(r→ are obtained as the following integral equation:(16)u′ir→=∫V1Gijr→,r→1∂Pjnr→1∂x1ndV1.Such representation of fluctuations of the local displacement field can be used in expressions for statistical descriptors of random fields-multipoint moments of various orders.


## 3. Statistical Descriptors of Local Stress and Strain Fields

In the general case of a multicomponent material, an indicator functionλCr→ takes the value 1 if the radius vector is in component C and is equal to 0 in all other cases. Then the local statistical moments of the first and second order in the component C of heterogeneous media are determined with the following relations [20]:(17)εijC=eij+1λCλ′Cr→ε′ijr→σijC=σij+1λCλ′Cr→σ′ijr→ε′ijr→ε′αβr→C=ε′ijr→ε′αβr→+eijeαβ-εijCεαβC+1λCλ′Cr→ε′ijr→ε′αβr→+eijλ′Cr→ε′αβr→+eαβλ′Cr→ε′ijr→σ′ijr→σ′αβr→C=σ′ijr→σ′αβr→+σijσαβ-σijCσαβC+1λCλ′Cr→σ′ijr→σ′αβr→+σijλ′Cr→σ′αβr→+σαβλ′Cr→σ′ijr→.Random fields of fluctuations of stresses and strains with regard to the statement of the boundary value problem are defined as(18)ε′ijr→=12u′i,jr→+u′j,ir→σ′ijr→=σijr→-σij=C′ijklr→ekl-C′ijklr→ε′klr→+Cijklr→ε′klr→.Thus, the moments for the unknown statistical characteristics depend on the fluctuations of the field of the indicator functionλ′C(r→), the displacements u′i(r→, and strains ε′ij(r→).Returning to two-phase materials, taking into account the latest transformations, the moments containing stress fluctuations can be expressed in terms of the moments containing strain fluctuations:(19)σ′ijr→σ′αβr→=ekleϕhC-ijklC-αβϕhDλ2+eklC-ijklCαβϕhλ′r→ε′ϕhr→+eϕhC-αβϕhCijklλ′r→ε′klr→+CijklCαβϕhε′ϕhr→ε′klr→-eklC-ijklC-αβϕh1-2pλ′r→ε′ϕhr→-eϕhC-ijklC-αβϕh1-2pλ′r→ε′klr→-C-αβϕhCijklλ′r→ε′ϕhr→ε′klr→-C-ijklCαβϕhλ′r→ε′ϕhr→ε′klr→+C-ijklC-αβϕh1-2pλ′r→ε′klr→ε′ϕhr→+Dλ2ε′klr→ε′ϕhr→λ′r→σ′ijr→=eklC-ijklDλ2-C-ijkl1-2pλ′r→ε′klr→+Cijklλ′r→ε′klr→λ′r→σ′ijr→σ′αβr→=ekleϕhC-ijklC-αβϕhDλ3+eklC-ijkl1-2pCαβϕhλ′r→ε′ϕhr→+eϕhC-αβϕh1-2pCijklλ′r→ε′klr→+CijklCαβϕhλ′r→ε′ϕhr→ε′klr→-eklC-ijklC-αβϕh1-3Dλ2λ′r→ε′ϕhr→-eϕhC-ijklC-αβϕh1-3Dλ2λ′r→ε′klr→-C-αβϕhCijkl+C-ijklCαβϕh1-2pλ′r→ε′klr→ε′ϕhr→+Dλ2ε′klr→ε′ϕhr→+C-ijklC-αβϕh1-3Dλ2λ′r→ε′klr→ε′ϕhr→+Dλ3ε′klr→ε′ϕhr→.Here, the tensorC-ijkl=CijklI-CijklM is the difference of the stiffness tensors of the inclusions and the matrix.In order to calculate the local moments of the strain and stress fields, it is necessary to know the analytical form of the momentε′kl(r→)ε′ϕh(r→) and also mixed unconditional moments λ′(r→)ε′kl(r→) and λ′(r→)ε′ϕh(r→)ε′kl(r→). These moments are not related to certain phase and characterize the RVE as homogeneous medium. All moments necessary for calculating the statistical characteristics depend on the geometry of the microstructure, its physical properties, and loading conditions. The physical properties of the components of the structure are determined by the tensors of the moduli of elasticity of inclusions and the matrix CijklI, CijklM. The geometry of the structure is characterized by the volume fraction p, the central moments of the indicator function Dλn, and multipoint high order correlation functions.


## 4. Integral Equations for Statistical Descriptors

After differentiating and expanding, expression (16) can be rewritten as(20)u′i,jr→=∫V1Gim,jr→,r→1C′mnklr→1ekl+C′mnklr→1u′k,lr→1,1ndV1.The displacement fluctuationsu′(r→, with respect to which the equation is being solved, are both on the left and on the right side of the expression, so the successive approximations are used for the solution. The recurrence sequence has the following form:(21)u′i,jχr→=∫V1Gim,jr→,r→1C′mnklr→1ekl+C′mnklr→1u′k,lχ-1r→1,1ndV1,where χ is order of approximation.This paper utilizes the second approximation, which contributes to the accuracy of the solution by taking into account multipoint structural correlation functions:(22)u′i,j2r→=ekl∫V1Gim,jr→,r→1C′mnklr→1,1ndV1+eoq∫V1∫V11Gim,jr→,r→1C′mnklr→1Gkf,lr→1,r→11C′fsoqr→11,11s,1ndV11dV1.The first term in expression (22) is the solution of the problem in the first approximation.Integrals in (21) and (22) as well as expressions for the moments ε′kl(r→)ε′ϕh(r→), λ′(r→)ε′kl(r→), λ′(r→)ε′ϕh(r→)ε′kl(r→) can be converted to the second derivative of Green’s function using Stieltjes convolution. The Stieltjes convolution of the two functions f(x) and g(x) over the interval (a,b) is defined as function Ψ(t):(23)Ψt≡∫abft-xdgx≡∫abgt-xdfxΨt≡∫abft-xg′xdx≡∫abgt-xf′xdx.Application of this definition to the integrals containing the first derivative of Green’s function and the derivative of the indicator function in the solutions of the boundary value problem (22) gives the next form of the integrals:(24)∫V1Gim,jr→,r→1C′mnklr→1,1ndV1=∫V1Gim,jnr→,r→1C′mnklr→1dV1∫V1∫V11Gim,jr→,r→1C′mnklr→1Gkf,lr→1,r→11C′mnklr→1,11s,1ndV11dV1=∫V1∫V11Gim,jnr→,r→1Gkf,lsr→1,r→11C′mnklr→1C′mnklr→1dV11dV1.For a two-phase inhomogeneous material, the solution using the second derivatives of the Green’s functions takes the form(25)u′i,j2r→=eklC-mnkl∫V1Gim,jnr→,r→1λ′r→1dV1+eoqC-mnklC-fsoq∫V1∫V11Gim,jnr→,r→1Gkf,lsr→1,r→11λ′r→1λ′r→11dV11dV1.The product of the fluctuations of the indicator functionλ′(r→1)λ′(r→11) is the correlation function Kλ2r→,r→1, which characterizes the relative position of the components in microstructure.By definition, Green’s functionGkf(r→,r→1 satisfying the expression (15) has a singularity at points r→-r→1=0. To apply numerical integration methods, the second derivative of Green’s function can be represented as the sum of the singular and formal components Gkf,ls(r→,r→1)=Gkf,lss(r→,r→1)+Gkf,lsf(r→,r→1).Then the integrals containing Green’s function can be decomposed into two terms [16]:(26)∫VGkf,lsr→,r→1Kλ2r→,r→1dV=∫V0Gkf,lssr→,r→1Kλ2r→,r→1dV+∫V-V0Gkf,lsfr→,r→1Kλ2r→,r→1dV,where V0 is the singularity region of the second derivative of Green’s function.For an isotropic tensorCijkl and in case when the elements of a representative volume are much smaller than the characteristic dimensions of the structure, Green’s function for an infinite medium can be used [1, 20]:(27)Gknr→,r→1=Aδknr→-r→1+Brk-r1krn-r1nr→-r→13,where rk denotes the kth component of the radius vector r→ and the coefficients A and B are linked with the effective Lame constants λ and μ by the relations(28)A=λ+3μ8πμλ+2μ,B=λ+μ8πμλ+2μ.Here,λ=λIp+λM(1-p, μ=μIp+μM(1-p.Expressions representing the singular and formal components for the second derivative of the tensor (27) were obtained in [16](29)Gkf,lssr→,r→1=-δr→-r→13μδkfδls-15κδkfls=δr→-r→1gkflsGkf,lsfr→,r→1=18πμr→-r→132-κδkf3nfns-δfs-κ3nkfnls+nklnfs+nksnfl+nfsnkl+nflnks-2Ikfls-15nknfnlns,where κ=λ+μ/λ+2μ, ni=ri-r1i/r→-r→1, nij=ninj, Ikfls=1/2δklδfs+δksδfl.Taking into account the definition of the generalized Dirac function∫δ(r→-r→1f(r→1)dr→1=fr→, the factor δ(r→-r→1 in the singular component of the second derivative of Green’s tensor leads to a transition from an integral to a constant value; that is, the integral operator Gkf,lsf(r→,r→1) becomes a constant tensor:(30)∫VGkf,lsr→,r→1Kλ2r→,r→1dV=gkflsDλ2+∫V-V0Gkf,lsfr→,r→1Kλ2r→,r→1dV.The momentsε′kl(r→)ε′ϕh(r→), λ′(r→)ε′kl(r→), λ′(r→)ε′ϕh(r→)ε′kl(r→) can be written in the second derivatives of Green’s function, taking into account the expansion (30). The integrals containing the Green functions and the moment functions here and below will be written in abbreviated form as follows:(31)I1,2,11,21r→,r→1,r→2,r→11,r→21=Dλ5∫V1∫V2∫V11∫V21Gim,jnfr→,r→1+Gjm,infr→,r→1Gαγ,βηfr→,r→2+Gβγ,αηfr→,r→2=Dλ5∫V1∫V2∫V11∫V21×Gkf,lsfr→x,r→11Gϕv,hwfr→x,r→21fλ5r→,r→1,r→2,r→11,r→21dV21dV11dV2dV1.Here the subscripts in parentheses denote the presence of the corresponding factors containing Green’s function in the integral and the upper indices determine the order and argument of the correlation function. The indexx of the radius vector is uniquely determined by the integrands and the Dirac convolution rules. The multiplicity of the integral is determined by the number of independent variables r→i. The absence of variables r→i in brackets means that there is no correlation function in the integrand.Then, in the second approximation of the problem solution, the mixed moments are written in the following form:(32)ε′ijr→ε′αβr→=14A1+2A2+A4λ′r→ε′ijr→=12B1+B2λ′r→ε′ijr→ε′αβr→2=14C1+2C2+C4.The termsA2, A4, B2, C2, C4in these expressions are an additive introduced by the second approximation that refines the solution in the first approximation.Taking into account the second derivative of Green’s function and the expansion into the singular and formal components, the terms of the sought moments have the following form:(33)A1=C-mnklC-γηϕhekleϕhI1,2r→1,r→2+gijmn+gjminI2r→,r→2+gαγβη+gβγαηI1r→,r→1+gijmn+gjmingαγβη+gβγαηDλ2A2=eklebdC-mnklC-γηϕhC-vwbdI1,2,21r→1,r→2,r→21+gimjn+gjminI2,21r→,r→2,r→21+gαγβη+gβγαηI1,21r→,r→1,r→21+gϕvhwI1,2r→,r→1,r→2+gαγβη+gβγαηgϕvhw1-2pI1r→,r→1+gimjn+gjmingϕvhw1-2pI2r→,r→2+gimjn+gjmingαγβη+gβγαη1-2pI21r→,r→21+gimjn+gjmingαγβη+gβγαηggvhwDλ3A4=eoqebdC-mnklC-γηϕhC-fsoqC-vwbdgimjn+gjmingαγβη+gβγαηgkflsgϕvhwDλ4+gimjn+gjmingαγβη+gβγαηgkfls1-3Dλ2I21r→,r→21+gimjn+gjmingkflsgϕvhw1-3Dλ2I2r→,r→2+gαγβη+gβγαηgkflsgϕvhw1-3Dλ2I1r→,r→1+gimjn+gjmingαγβη+gβγαη1-2pI11,21r→,r→11,r→21+Dλ2I11,21r→11,r→21+gimjn+gjmingkfls1-2pI2,21r→,r→2,r→21+Dλ2I2,21r→2,r→21+gimjn+gjmingϕvhw1-2pI2,11r→,r→2,r→11+Dλ2I2,11r→2,r→11+gαγβη+gβγαηgkfls1-2pI1,21r→,r→1,r→21+Dλ2I1,21r→1,r→21+gαγβη+gβγαηgϕvhw1-2pI1,11r→,r→1,r→11+Dλ2I1,11r→1,r→11+gkflsgϕvhw1-2pI1,2r→,r→1,r→2+Dλ2I1,2r→1,r→2+gimjn+gjminI2,11,21r→,r→2,r→11,r→21+gαγβη+gβγαηI1,11,21r→,r→1,r→11,r→21+gkflsI1,2,21r→,H→1,r→2,r→21+gϕvhwI1,2,11r→,r→1,r→2,r→11+I1,2,11,21r→1,r→2,r→11,r→21B1=eklC-mnklgimjn+gjminDλ2+I1r→,r→1B2=eoqC-mnklC-fsoqgimjn+gjmingkflsDλ3+gimjn+gjmin1-2pI11r→,r→11+gkfls1-2pI1r→,r→1+I1,11r→,r→1,r→11C1=ekleϕhC-mnklC-γηϕhI1,2r→,r→1,r→2+gimjn+gjmin1-2pI2r→,r→2+gαγβη+gβγαη1-2pI1r→,r→1+gimjn+gjmingαγβη+gβγαηDλ2C2=eklebdC-mnklC-γηϕhC-vwbdgimjn+gjmingαγβη+gβγαηgϕvhwDλ4+gimjn+gjmingαγβη+gβγαη1-3Dλ2I21r→,r→21+gimjn+gjmingϕvhw1-3Dλ2I2r→,r→2+gαγβη+gβγαηgϕvhw1-3Dλ2I1r→,r→1+gimjn+gjmin1-2pI2,21r→,r→2,r→21+Dλ2I2,21r→2,r→21+gαγβη+gβγαη1-2pI1,21r→,r→1,r→21+Dλ2I1,21r→1,r→21+gϕvhw1-2eI1,2r→,r→1,r→2+Dλ2I1,2r→1,r→2+I1,2,21r→,r→1,r→2,r→21C4=eoqebdC-mnklC-γηϕhC-fsoqC-vwbdgimjn+gjmingαγβη+gβγαηgkflsgϕvhwDλ5+gimjn+gjmingαγβη+gβγαηgkflsp-1Dλ4-p3I21r→,r→21+gimjn+gjmingαγβη+gβγαηgϕvhwp-1Dλ4-p3I11r→,r→11+gimjn+gjmingkflsgϕvhwp-1Dλ4-p3I2r→,r→2+gαγβη+gβγαηgkflsgϕvhwp-1Dλ4-p3I1r→,r→1+gimjn+gjmingαγβη+gβγαη1-3Dλ2I11,21r→,r→11,r→21+Dλ3I11,21r→11,r→21+gimjn+gjmingkfls1-3Dλ2I2,21r→,r→2,r→21+Dλ3I2,21r→2,r→21+gimjn+gjmingϕvhw1-3Dλ2I2,11r→,r→2,r→11+Dλ3I2,11r→2,r→11+gαγβη+gβγαηgkfls1-3Dλ2I1,21r→,r→1,r→21+Dλ3I1,21r→1,r→21+gαγβη+gβγαηgϕvhw1-3Dλ2I1,11r→,r→1,r→11+Dλ3I1,11r→1,r→11+gkflsgϕvhw1-3Dλ2I1,2r→,r→1,r→2+Dλ3I1,2r→1,r→2+gimjn+gjmin1-2pI2,11,21r→,r→2,r→11,r→21+Dλ2I2,11,21r→2,r→11,r→21+p3Dλ2I2,11,21+gαγβη+gβγαη1-2pI1,11,21r→,r→1,r→11,r→21+Dλ2I1,11,21r→1,r→11,r→21+p3Dλ2I1,11,21+gkfls1-2pI1,2,21r→,r→1,r→2,r→21+Dλ2I1,2,21r→1,r→2,r→21+p3Dλ2I1,2,21+gϕvhw1-2pI1,2,11r→,r→1,r→2,r→11+Dλ2I1,2,11r→1,r→2,r→11+p3Dλ2I1,2,11+I1,2,11,21r→,r→1,r→2,r→11,r→21.These expressions explicitly contain the structural correlation functions up to the 5th order and represent a generalized case of expressions for the statistical characteristics of fields. In previously obtained results, the correlation functions of second- and third-order maximum were taken into account based on first-order approximations of the solution of a boundary value problem. For instance, in the two-point approximation only pair interactions are taken into account; therefore correlation functions of higher order than binary ones are excluded [16, 20].Eliminating all the terms containing the formal component from the last expressions, a singular approximation for statistical characteristics is derived [16]. Because of negligibility of terms containing correlation functions, the singular approximation does not take into account the real distribution of components in a representative volume but depends only on the volume fraction:(34)ε′ijr→ε′αβr→=14C-mnklC-γηϕhekleϕhgijmn+gjmingαγβη+gβγαηDλ2+2eklebdC-mnklC-γηϕhC-vwbdgimjn+gjmingαγβη+gβγαηggvhwDλ3+eoqebdC-mnklC-γηϕhC-fsoqC-vwbdgimjn+gjmingαγβη+gβγαηgkflsgϕvhwDλ4λ′r→ε′ijr→=12eklC-mnklgimjn+gjminDλ2+eoqC-mnklC-fsoqgimjn+gjmingkflsDλ3λ′r→ε′ijr→ε′αβr→=14ekleϕhC-mnklC-γηϕhgimjn+gjmingαγβη+gβγαηDλ2+2eklebdC-mnklC-γηϕhC-vwbdgimjn+gjmingαγβη+gβγαηgϕvhwDλ4+eoqebdC-mnklC-γηϕhC-fsoqC-vwbdgimjn+gjmingαγβη+gβγαηgkflsgϕvhwDλ5.Unlike the singular approximation, multipoint methods allow to take into account the morphological features of the microstructure of inhomogeneous media by means of a set of correlation functions.It is obvious that the contribution made by the second approximation of the solution of the SBVP is significant in terms of the complication of analytic expressions; however, it contains terms that do not depend on the correlation functions and which can be neglected during numerical integration due to the insignificance of their contribution to the values of the moments.Numerical calculation of the obtained expressions can be simplified by specifying types of loading due to the symmetry properties of the tensors in the expressions. Numerical calculation of the values of the multidimensional integrals can be implemented in the commercial packages, which offer a large variation of integration schemes. Singularity isolation allows integration without additional analysis of integrands. That significantly increases the computational speed. The solution of the obtained integral equations can be effectively parallelized.Thus, by specifying the elastic characteristics of the components of the composite, the stress-strain state (using the macroscopic strain tensor), knowing the geometry of the structure (the volume fraction, the size of the inclusions) and the analytical expressions for the microstructural correlation functions, the moments of the first- and second-orders of the fields of local strain and stress can be calculated using methods of numerical integration.Some examples of the numerical implementation of the derived stochastic micromechanical model are presented in [19, 25], where calculations for a case study of elastic and elastoplastic composites with different microstructural properties were obtained for various loading conditions of the RVE. It is worth mentioning that contribution of the second approximation of the integral equation (21) is more significant when the correlation functions of 4th and 5th order demonstrate more oscillatory behavior. The latter happens when inclusions are unevenly distributed inside the RVE and due to the fact that correlation functions of each subsequent order appeared to be more sensitive to the features of microstructural morphology. Empirically it was established that the third-order functions reflect changes in the form of inclusions, while fourth-order functions describe how inclusions are grouped and how their sizes are distributed. Thus, it was shown that the second approximation term in (22) almost does not make any addition for RVEs with uniformly distributed inclusions (the difference between results of local fields modeling between the first and second solution approximation was around 3%). However, for more complex structures, when inclusions form clusters, the effect of the microstructure becomes more valuable and is explicitly accounted in the model with help of the high order correlation functions.


## 5. Conclusions

The paper presents the derivation of analytical expressions for the moments of the first and second order of local stress and strain in the components of a two-phase structurally inhomogeneous medium using multipoint correlation functions and second derivatives of Green’s functions. The advantage of this method is that it allows to take into account the information about the internal structure of the medium, formalized by multipoint correlation functions, with high precision.The obtained statistical moments have a number of practical applications. The statistical description of local microstructural fields of stresses and strains is one of the tools for assessing the initiation of damage growth and understanding of nonlinear mechanical phenomena. Based on the statistical metrics, a model for the multiscale failure probability calculation can be implemented. Representation of the failure criteria in a statistical form is required for that, followed by statistical evaluation of the values of the critical quantities. Homogenization methods based on multipoint statistical instruments allow obtaining a fast numerical estimate for effective properties with allowance for microscale interactions.



---

*Source: 5785193-2018-05-06.xml*