Periodization of random media and representative volume element size for linear composites

Several existing numerical studies show that the effective linear properties of random composites can be accurately estimated using small volumes subjected to periodic boundary conditions – more suitable than homogeneous strain or stress boundary conditions – providing that a sufficient number of realizations are considered. Introducing the concept of periodization of random media, this Note gives a new definition of representative volume element which leads to estimates of its minimum size in agreement with existing theoretical results. A qualitative convergence criterion for the numerical simulations is proposed and illustrated with finite element computations. To cite this article: K. Sab, B. Nedjar, C. R. Mecanique 333 (2005).     

Determination of the size of the representative volume element for random composites: statistical and numerical approach

The representative volume element (RVE) plays a central role in the mechanics and physics of random heterogeneous materials with a view to predicting their effective properties. A quantitative definition of its size is proposed in this work. A RVE size can be associated with a given precision of the estimation of the wanted overall property and the number of realizations of a given volume V of microstructure that one is able to consider. It is shown to depend on the investigated morphological or physical property, the contrast in the properties of the constituents, and their volume fractions. The methodology is applied to a specific random microstructure, namely a two-phase three-dimensional Voronoı̈ mosaic. Finite element simulations of volumes of different sizes are performed in the case of linear elasticity and thermal conductivity. The volumes are subjected to homogeneous strain, stress or periodic boundary conditions. The effective properties can be determined for large volumes and a small number of realizations. Conversely, smaller volumes can be used providing that a sufficient number of realizations are considered. A bias in the estimation of the effective properties is observed for too small volumes for all types of boundary conditions. The variance of computed apparent properties for each volume size is used to define the precision of the estimation. The key-notion of integral range is introduced to relate this error estimation and the definition of the RVE size. For given wanted precision and number of realizations, one is able to provide a minimal volume size for the computation of effective properties. The results can also be used to predict the minimal number of realizations that must be considered for a given volume size in order to estimate the effective property for a given precision. The RVE sizes found for elastic and thermal properties, but also for a geometrical property like volume fraction, are compared.     

Residual stresses in random elastic composites: nonlocal micromechanics-based models and first estimates of the representative volume element size

Random elastic composites with residual stresses are examined in this paper with the aim of understanding how the prestress may influence the overall mechanical properties of the composite. A fully non-local effective response is found in perfect analogy with the un-prestressed case examined in (Drugan and Willis, J. Mech. Phys. Solids 44(4):497–524, 1996). The second gradient approximation is considered and the impact of the residual stresses on the estimate of the RVE size is studied whenever the local response is used to describe the mechanical properties of the heterogeneous medium. To this aim, total and incremental formulations are worked out in this paper and the influence of both uniform and spatially varying prestresses are studied. Among other results, it is shown how rapid oscillations of relatively “small” residual stresses in most cases may result in the impossibility of describing the overall behavior of the composite with a local constitutive equation. On the other hand, prestresses with relatively high amplitudes and slow spatial oscillations may even reduce the RVE size required for approximating the mechanical properties of un-prestressed heterogeneous media with a local constitutive equation.     

On the minimum size of representative volume element: An experimental investigation

In this paper, the minimum size of the representative volume element (RVE) of a heterogeneous material is determined experimentally using the digital image correlation technique. Experiments of uniaxial compression and thermal expansion were conducted on PBS 9501, a high explosive simulant material. The minimum size of the RVE of the PBS 9501 heterogeneous material, where the average crystal diameter of the material is of the order of 100 μm, was determined to be approximately 1.5 mm. This result is consistent with numerical calculations on polycrystalline materials and some other composites.     

Generalized self-consistent estimation of the apparent isotropic elastic moduli and minimum representative volume element size of heterogeneous media

Under investigation is a heterogeneous material consisting of an elastic homogeneous isotropic matrix in which layered elastic isotropic inclusions or pores are embedded. The generalized self-consistent model (GSCM) is extended so as to be capable of estimating the apparent elastic properties of a finite-size specimen smaller than a representative volume element (RVE). The kinematical or static apparent shear modulus is determined as a root of a cubic polynomial equation instead of a quadratic polynomial equation as in the classical GSCM of Christensen and Lo [Christensen, R.M., Lo, K.H., 1979. Solutions for effective shear properties in three phase sphere and cylinder models. J. Mech. Phys. Solids 27, 315–330]. It turns out that the extended GSCM establishes a link between the composite sphere assemblage model (CSAM) of Hashin [Hashin, Z., 1962. The elastic moduli of heterogeneous materials. J. Appl. Mech. 29, 143–150] and the classical GSCM. Demanding that the normalized distance between the kinematical and static apparent moduli of a finite-size specimen be smaller than a certain tolerance, the minimum RVE size is estimated in a closed form.     

Exterior statistics based boundary conditions for representative volume elements of elastic composites

Statistically equivalent representative volume elements or SERVEs are representations of the microstructure that are used for micromechanical simulations to generate homogenized material constitutive responses and properties (Swaminathan et al., 2006a, Ghosh, 2011). Typically, a SERVE is generated from the parent microstructure as a statistically equivalent region, whose size is determined from the requirements of convergence of macroscopic properties. Standard boundary conditions, such as affine transformation-based displacement boundary conditions (ATDBCs), uniform traction boundary conditions (UTBCs) or periodic boundary conditions (PBCs) are conventionally applied on the SERVE boundary for micromechanical simulations. However, when the microstructure is characterized by arbitrary, nonuniform distributions of heterogeneities, these simple boundary conditions do not represent the effect of regions exterior to the SERVE. Improper boundary conditions can result in significantly larger than optimal SERVE domains, needed for converged properties. In an attempt to overcome the limitations of the conventional boundary conditions on the SERVE, this paper explores the effect of boundary conditions that incorporate the statistics of the exterior region on the SERVE of elastic composites. Using Green's function based interaction kernels, coupled with statistical functions of the microstructural characteristics like one-point and two-point correlation functions, a novel exterior statistics-based boundary condition or ESBC is derived for the SERVE. The advantages of the ESBC are established by comparing with results of simulations using conventional boundary conditions. Results of the SERVE simulations subjected to ESBCs are also compared with those from other popular methods like statistical volume element (SVE) and weighted statistical volume element (WSVE). The proposed ESBCs offer significant advantages over other methods in the SERVE-based analysis of heterogeneous materials.     

A unified periodical boundary conditions for representative volume elements of composites and applications

An explicit unified form of boundary conditions for a periodic representative volume element (RVE) is presented which satisfies the periodicity conditions, and is suitable for any combination of multiaxial loads. Starting from a simple 2-D example, we demonstrate that the “homogeneous boundary conditions” are not only over-constrained but they may also violate the boundary traction periodicity conditions. Subsequently, the proposed method is applied to: (a) the simultaneous prediction of nine elastic constants of a unidirectional laminate by applying multiaxial loads to a cubic unit cell model; (b) the prediction of in-plane elastic moduli for [±θ]_n angle-ply laminates. To facilitate the analysis, a meso/micro rhombohedral RVE model has been developed for the [±θ]_n angle-ply laminates. The results obtained are in good agreement with the available theoretical and experimental results.     

Determining a representative element volume for DEM simulations of samples with non-circular particles

Numerical studies on the number of particles or system size required to attain a representative element volume (REV) for discrete element method (DEM) simulations of granular materials have almost always considered samples with spherical or circular particles. This study considers how many particles are needed to attain a REV for 2D samples of 2-disc cluster particles where the particle aspect ratio (AR) was systematically varied. Dense and loose assemblies of particles were simulated. The minimum REV was assessed both by considering the repeatability of static packing characteristics and the shearing behaviour in biaxial compression tests, and by investigating the effect of sample size on the measured characteristics and observed shearing behaviour. The repeatability of the data considered generally improved with increasing sample size. The packing characteristics of the dense samples were more repeatable suggesting that the minimum REV reduces with increasing packing density. The minimum REV was observed to be sensitive to the characteristic measured. Although the overall responses of the samples during shear deformation were similar irrespective of the sample sizes, the smaller the sample size, the higher the fluctuations observed in the responses. Analysis of the coefficient of variation of the fluctuations around the critical state stress ratio can provide insight as to whether a REV is attained. The particle AR influences the effect of sample size on shearing characteristics and thus the minimum number of particles required to attain a REV; this can be explained by the influence of AR on the number of contacts within the samples.     

Identification of the representative crack length evolution in a multi-level interface model for quasi-brittle masonry

It is proposed here to identify the law of crack length evolution with a small number of parameters governing a recently presented model (Rekik and Lebon, submitted for publication) describing the interface behavior in damaged masonry. Studies on non-confined medium- and large-sized masonry structures have shown that it is necessary to obtain a linear increasing crack in the post-peak part of the “stress–strain or –displacement” diagram. In confined masonry structures showing softening and sliding parts, the results obtained with this crack evolution failed to match the experimental data. The crack lengths identified in the post-peak part at several points on the experimental “stress–displacement” diagram show that the representative crack length is a bilinear or trilinear function describing the increase in the crack length with respect to the decrease in the shear stress. Numerical studies on medium- and large-sized masonry structures consisting of the same materials subjected to various loads were performed to determine the ultimate crack length, and the results are relatively insensitive to the size of the masonry and the type of the load applied. The numerical local fields determined in the elementary and full-scale structures investigated were used to test the validity of the present model at the local scale, as well as to obtain an additional unilateral condition in the case of compressed masonry structures in order to prevent overlapping between the masonry components.     

Influence of the particle size and volume fraction on micro-damage of the composites

The bulk and shear modulus of metal matrix composites with various volume fractions of particles are modified based on the Eshelby’s equivalent inclusion method combined with self-consistent scheme. By introducing the modified modulus, a new model, which can predict the particle size effects on the stress–strain relation under interfacial debonding damage between matrix and particles, is established. The results obtained from the present investigation show a better agreement with the experimental data.     

Assessment of existing and introduction of a new and robust efficient definition of the representative volume element

Accurate numerical homogenization necessitates the thorough determination of the Representative Volume Element (RVE). There exists several seminal works on the notion of the RVE in homogenization, its definitions and methods of determination for efficient computation of composite effective properties. The objective of the current work is to assess the ability of numerical RVE determination methods to deliver accurate effective properties of composite materials. This paper demonstrates that common and well-established RVE determination methods, based on studying the convergence rate of the effective properties with respect to the volume element size, are invalid for the case of composites reinforced by randomly oriented fibers and yield erroneous estimates of their effective properties. Following the failure of traditional RVE determination methods, we proposed a new RVE determination criterion that is not based on the average property stability, but its statistical variations. Our new proposed criterion has been shown to be more accurate than other criteria in computing the effective properties of composites for aspect ratios up to 60. Moreover, the proposed criterion does not necessitate a convergence study over the volume element size, hence reducing considerably the RVE determination cost. Finally, our work questions the validity of many published works dealing with composites including heterogeneities of high aspect ratios.     

Quantification of stochastically stable representative volumes for random heterogeneous materials

This paper details a procedure to determine lower bounds on the size of representative volume elements (RVEs) by which the size of the RVE can be quantified objectively for random heterogeneous materials. Here, attention is focused on granular materials with various distributions of inclusion size and volume fraction of inclusions. An extensive analysis of the RVE size dependence on the various parameters is performed. Both deterministic and stochastic parameters are analysed. Also, the effects of loading mode and the parameter of interest are studied. As the RVE size is a function of the material, some material properties such as Young's modulus and Poisson's ratio are analysed as factors that influence the RVE size. The lower bound of RVE size is found as a function of the stochastically distributed volume fraction of inclusions; thus the stochastic stability of the obtained results is assessed. To this end a newly defined concept of stochastic stability (DH-stability) is introduced by which stochastic effects can be included in the stability considerations. DH-stability can be seen as an extension of classical Lyapunov stability. As is shown, DH-stability provides an objective tool to establish the lower bound nature of RVEs for fluctuations in stochastic parameters.     

A representative volume element based on translational symmetries for FE analysis of cracked laminates with two arrays of cracks

A methodology is proposed for the construction of a representative volume element (RVE) for analysis of laminated composites containing two arrays of ply cracks running in different directions. The only requirement is that the cracks in any ply are uniformly spaced, and if more than one ply of a given orientation is cracked, then the crack spacing of individual plies must only be in exact multiples of each other. The spacing of cracks in the two directions can be fully independent. The RVE is constructed through a systematic consideration of translational symmetries present in the cracked laminate. As a result, the boundary conditions on the RVE can be imposed without compromising accuracy. Examples of the application of the RVE methodology are given to illustrate its broad capability and a finite element (FE) stress analysis is performed for these cases to illustrate results such as the crack surface displacements, local stress fields and RVE-averaged elastic properties. For one case, the average properties are compared with experimental results, showing good agreement.     

高导无氧铜杆的冲击拉伸断裂系列实验以及基于样本体积单元的分析

采用Hopkinson装置和一种基于一级气体炮的高速冲击拉伸断裂装置,研究了无刻槽高导无氧铜 (OFHC)杆在一系列冲击拉伸速度下的断裂。当冲击拉伸速度大于40m/s时,断裂位置总在冲击拉伸端附 近,此速度被确定为OFHC的实验临界冲击拉伸速度。一种受单轴冲击拉伸荷载的、中心含椭球空穴的样本 体积单元被用于数值模拟所含空穴的增长与失稳的过程。OFHC的J-C与Z-A 本构关系用于描述基体材料 的动态响应。讨论了空穴失稳条件并提出以空穴形状演化为判据,比较了空穴失稳时的样本体积单元平均径 向应变与无刻槽杆的冲击断裂应变。也用这种样本体积单元模型分析了OFHC的实验临界冲击拉伸速度。     

Monte Carlo simulation of complex cohesive fracture in random heterogeneous quasi-brittle materials

A numerical method is developed to simulate complex two-dimensional crack propagation in quasi-brittle materials considering random heterogeneous fracture properties. Potential cracks are represented by pre-inserted cohesive elements with tension and shear softening constitutive laws modelled by spatially-varying Weibull random fields. Monte Carlo simulations of a concrete specimen under uni-axial tension were carried out with extensive investigation of the effects of important numerical algorithms and material properties on numerical efficiency and stability, crack propagation processes and load-carrying capacities. It was found that the homogeneous model led to incorrect crack patterns and load–displacement curves with strong mesh-dependence, whereas the heterogeneous model predicted realistic, complicated fracture processes and load-carrying capacity of little mesh-dependence. Increasing the variance of the tensile strength random fields with increased heterogeneity led to reduction in the mean peak load and increase in the standard deviation. The developed method provides a simple but effective tool for assessment of structural reliability and calculation of characteristic material strength for structural design.     

Effects of deformed volume, volume fraction and particle size on the deformation behaviour of W/Cu composites

Tungsten/copper (W/Cu) particle reinforced composites were used to investigate the scaling effects on the deformation and fracture behaviour. The effects of the volume fraction and the particle size of the reinforcement (tungsten particles) were studied. W/Cu-80/20, 70/30 and 60/40 wt.% each with tungsten particle size of 10 μm and 30 μm were tested under compression and shear loading. Cylindrical compression specimens with different volumes (D_S = H) were investigated with strain rates between 0.001 s⁻¹ and about 5750 s⁻¹ at temperatures from 20 °C to 800 °C. Axis-symmetric hat-shaped shear specimens with different shear zone widths were examined at different strain rates as well. A clear dependence of the flow stress on the deformed volume and the particle size was found under compression and shear loading. Metallographic investigation was carried out to show a relation between the deformation of the tungsten particles and the global deformation of the specimens. The size of the deformed zone under either compression or shear loading has shown a clear size effect on the fracture of the hat-shaped specimens.The quasi-static flow curves were described with the material law from Swift. The parameters of the material law were presented as a function of the temperature and the specimen size. The mechanical behaviour of the composite materials were numerically computed for an idealized axis-symmetric hat-shaped specimen to verify the determined material law.