Stochastic multi-scale analysis of homogenised proper ties considering uncertainties in cellular solid micros- tructures using a first-order perturbation (original) (raw)
Related papers
Micromechanically based stochastic finite elements: length scales and anisotropy
Probabilistic Engineering Mechanics, 1996
The present stochastic finite element (SFE) study amplifies a recently developed micromechanically based approach in which two estimates (upper and lower) of the finite element stiffness matrix and of the global response need first to be calculated. These two estimates correspond, respectively, to the principles of stationary potential and complementary energy on which the SFE is based. Both estimates of the stiffness matrix are anisotropic and tend to converge towards one another only in the infinite scale limit; this points to the fact that an approximating meso-scale continuum random field is neither unique nor isotropic. The SFE methodology based on this approach is implemented in a Monte Carlo sense for a conductivity (equivalently, out-of-plane elasticity) problem of a matrix-inclusion composite under mixed boundary conditions. Two versions are developed: in one an exact calculation of all the elements' stiffness matrices from the microstructure over the entire finite element mesh is carried out, while in the second one a second-order statistical characterization of the mesoscale continuum random field is used to generate these matrices.
Computer Methods in Applied Mechanics and Engineering, 2019
An inverse Mean-Field Homogenization (MFH) process is developed to improve the computational efficiency of non-linear stochastic multiscale analyzes by relying on a micro-mechanics model. First full-field simulations of composite Stochastic Volume Element (SVE) realizations are performed to characterize the homogenized stochastic behavior. The uncertainties observed in the non-linear homogenized response, which result from the uncertainties of their micro-structures, are then translated to an incrementalsecant MFH formulation by defining the MFH input parameters as random effective properties. These effective input parameters, which correspond to the micro-structure geometrical information and to the material phases model parameters, are identified by conducting an inverse analysis from the full-field homogenized responses. Compared to the direct finite element analyzes on SVEs, the resulting stochastic MFH process reduces not only the computational cost, but also the order of uncertain parameters in the composite micro-structures, leading to a stochastic Mean-Field Reduced Order Model (MF-ROM). A data-driven stochastic model is then built in order to generate the random effective properties under the form of a random field used
Stochastic finite element analysis of composite structures based on material microstructure
Composite Structures, 2015
The linking of microstructure uncertainty with the random variation of material properties at the macroscale is particularly needed in the framework of the stochastic finite element method (SFEM) where arbitrary assumptions are usually made regarding the probability distribution and correlation structure of the macroscopic mechanical properties. This linking can be accomplished in an efficient manner by exploiting the excellent synergy of the extended finite element method (XFEM) and Monte Carlo simulation (MCS) for the computation of the effective properties of random two-phase composites. The homogenization is based on Hill's energy condition and involves the generation of a large number of random realizations of the microstructure geometry based on a given volume fraction of the inclusions and other parameters (shape, number, spatial distribution and orientation). In this paper, the mean value, coefficient of variation and probability distribution of the effective elastic modulus and Poisson ratio are computed taking into account the material microstructure. The effective properties are used in the framework of SFEM to obtain the response of a composite structure and it is shown that the response variability can be significantly affected by the random microstructure.
Perturbation stochastic finite element-based homogenization of polycrystalline materials
Journal of Mechanics of Materials and Structures, 2011
This work presents a study of the influence on the macroscopic (homogenized) elastic properties of polycrystalline materials induced by uncertainties in the material texture and microstructure geometry. Since many microelectromechanical systems are made of materials deposited as thin films with < < 111> > fiber texture, we study the variance of the homogenized elastic properties of the material around its nominal < < 111> > texture. To perform this analysis, the perturbation stochastic finite element method (PSFEM) is coupled to the mathematical theory of homogenization leading to a second-order perturbation-based homogenization method. This method is able to evaluate the mean and variance of a given homogenized property as a function of the grain property uncertainty. The multiscale formulation is implemented in a plane-stress linear elastic finite element framework based on a multigrain periodic unit cell generated by Voronoi tessellation. This perturbation-based homogenization method is verified against Monte Carlo simulations, showing its effectiveness and limitations. Then, through applications, it is evaluated how different levels of uncertainty in grains induce uncertainty in the macroscopic elastic properties of the polycrystalline material. In particular, the influence of the unit cell is studied. Finally, by coupling the PSFEM with the Monte Carlo method, the effects on the macroscopic properties of uncertainty of both the geometry and orientation of the grains is estimated.
Computational Mechanics, 2014
The computationally random homogenization analysis of a two-phase heterogeneous materials is addressed in the context of linear elasticity where the randomness of constituents' moduli and microstructural morphology together with the correlation among random moduli are fully considered, and random effective quantities such as effective elastic tensor and effective stress as well as effective strain energy together with their numerical characteristics are then sought for different boundary conditions. Based on the finite element method and Monte-carlo method, the RVE with randomly distributing particles determined by a numerical convergence scheme is firstly generated and meshed, and two types of boundary conditions controlled by average strain are then applied to the RVE where the uncertainty existing in the microstructure is accounted for simultaneously. The numerical characteristics of random effective quantities such as coefficients of variation and correlation coefficients are then evaluated, and impacts of different factors on random effective quantities are finally investigated and revealed as well.
Micromechanically based stochastic finite elements
Finite Elements in Analysis and Design, 1993
A stochastic finite element method for analysis of effects of spatial variability of material properties is developed with the help of a micromechanics approach. The method is illustrated by evaluating the first and second moments of the global response of a membrane with microstructure of a spatially random inclusionmatrix composite under a deterministic uniformly distributed load. It is shown that two mesoscale random continuum fields have to be introduced to bound the material properties and, in turn, the global response from above and from below. The intrinsic scale dependence of these two random fields is dictated by the choice of the finite element mesh.
ASME J. Risk Uncertainty Part B, 2019
This paper addresses the effect of microstructure uncertainties on elastic properties of nanocomposites using finite element analysis (FEA) simulations. Computer-simulated microstructures were generated to reflect the variability observed in nanocomposite microstructures. The effect of waviness, agglomeration, and orientation of carbon nanotubes (CNTs) were investigated. Generated microstructures were converted to image-based 2D FEA models. Two hundred different realizations of microstructures were generated for each microstructure type to capture the stochastic response. The results confirm previously reported findings and experimental results. The results show that for a given fiber volume fraction, CNTs orientation, waviness, and agglomeration result in different elastic properties. It was shown that while a given microstructural feature will improve the elastic property, it will increase the variability in the elastic properties.
Scale-dependent homogenization of random composites as micropolar continua
European Journal of Mechanics - A/Solids, 2015
A multitude of composite materials ranging from polycrystals to rocks, concrete, and masonry overwhelmingly display random morphologies. While it is known that a Cosserat (micropolar) medium model of such materials is superior to a Cauchy model, the size of the Representative Volume Element (RVE) of the effective homogeneous Cosserat continuum has so far been unknown. Moreover, the determination of RVE properties has always been based on the periodic cell concept. This study presents a homogenization procedure for disordered Cosserat-type materials without assuming any spatial periodicity of the microstructures. The setting is one of linear elasticity of statistically homogeneous and ergodic two-phase (matrix-inclusion) random microstructures. The homogenization is carried out according to a generalized Hill-Mandel type condition applied on mesoscales, accounting for non-symmetric strain and stress as well as couple-stress and curvature tensors. In the setting of a two-dimensional elastic medium made of a base matrix and a random distribution of disk-shaped inclusions of given density, using Dirichlet-type and Neumann-type loadings, two hierarchies of scale-dependent bounds on classical and micropolar elastic moduli are obtained. The characteristic length scales of approximating micropolar continua are then determined. Two material cases of inclusions, either stiffer or softer than the matrix, are studied and it is found that, independent of the contrast in moduli, the RVE size for the bending micropolar moduli is smaller than that obtained for the classical moduli. The results point to the need of accounting for: spatial randomness of the medium, the presence of inclusions intersecting the edges of test windows, and the importance of additional degrees of freedom of the Cosserat continuum.
Computer Methods in Applied Mechanics and Engineering, 2013
Recognizing that modern materials contain multiple phases with inherent random microstructure and in situ constituent material properties that are oft uncharacterizable with exactness, this research uses benchmark computational studies to unveil scenarios where uncertainties significantly affect macroscopic behavior. The benchmark studies, which serve as numerical experiments capturing the main features of a wide class of problems in solid mechanics, suggest a generalized uncertainty propagation criterion D whose assessment may be used to understand whether uncertainties may non-negligibly propagate to apparent system properties. The new D-criterion combines four features of a microstructured material system: the microstructure size (micro), material property correlation length (micro), structure size (macro), and global length scale of loading (macro). We compare the D-criterion with the familiar one from classical homogenization that relates a statistical volume element size to the size of the microstructure it contains. This work is motivated by the knowledge paucity of uncertainty's place in multiscale analysis and the great computational expense, either in runtime or implementation, of uncertainty propagation methods. The uncertainties considered are material property uncertainty of the separate phases comprising each problem's microstructure and the statistical description of the microstructure itself. The problems involve quasi-static elasto-plasticity, wave propagation in viscoelastic composite, and pattern transition in periodic voided elastomer. Besides the criterion, observed physical behavior for the benchmark problems is also presented in the context of input uncertainty. A close examination of the boundary condition effect on uncertainty propagation suggests that localization increases the characteristic microstructure length, which we represent with a simple deformationdependent microstructure evolution function and that increases uncertainty in apparent system properties. This work provides a systematic framework for modelers to understand when they may reliably waive the consideration of uncertainty for their application or when their application mandates it.