A stochastic model for elasticity tensors with uncertain material symmetries

In this paper, we consider the probabilistic modeling of media exhibiting uncertainties on material symmetries. More specifically, we address both the construction of a stochastic model and the definition of a methodology allowing the numerical simulation (and consequently, the inverse experimental identification) of random elasticity tensors whose mean distance (in a sense to be defined) to a given class of material symmetry is specified. Following the eigensystem characterization of the material symmetries, the proposed approach relies on the probabilistic model derived in Mignolet and Soize (2008), allowing the variance of selected eigenvalues of the elasticity tensor to be partially prescribed. In this context, a new methodology (regarding in particular the parametrization of the model) is defined and illustrated in the case of transversely isotropic materials. The efficiency of the approach is demonstrated by computing the mean distance of the random elasticity tensor to a given material symmetry class, the distance and projection onto the space of transversely isotropic tensors being defined by considering the Riemmanian metric and the Euclidean projection, respectively. It is shown that the methodology allows the above distance to be (partially) reduced as the overall level of statistical fluctuations increases, no matter the initial distance of the mean model used in the simulations. A comparison between this approach and the initial nonparametric approach introduced in Soize (2008) is finally provided.     

Applications of tensor functions to the formulation of yield criteria for anisotropic materials

Yielding of anisotropic materials can be characterized by yield criteria which are scalar-valued functions of the stress tensor and of material tensors, for instance, of rank two or four, characterizing the anisotropic properties of the material. Because of the requirement of invariance, a yield criterion can be expressed as a single-valued function of the integrity basis. In finding an integrity basis involving the stress tensor and material tensors, the constitutive equations are first formulated based on the tensor function theory. Since the plastic work characterizes the yield process, we read from this scalar expression the essential invariants to formulate a yield criterion. Some examples for practical use are discussed in detail.     

Representation Theorems for Hemitropic and Tranversely Isotropic Tensor Functions

A representation theorem for transversely isotropic tensor-valued functions of a symmetric tensor variable is proved. The theorem holds in any finite dimension. The proof is based on the decomposition of a symmetric tensor of dimension N into a scalar, a vector, and a symmetric tensor of dimension N-1, and on the fact that the transverse isotropy of the original function is equivalent to the hemitropy of three functions, one scalar-valued, one vector-valued, and one tensor-valued, of the last two terms in the decomposition. Representation theorems for the three functions are obtained as generalizations of two theorems of W. Noll on isotropic functions. The proofs make use of an appropriate algebraic structure based on alternating forms. The three-dimensional case, as well as those of linear and of hyperelastic functions, are treated as special cases. This revised version was published online in August 2006 with corrections to the Cover Date.     

Membrane Curvatures and Stress-strain Full Fields of Axisymmetric Bulge Tests from 3D-DIC Measurements. Theory and Validation on Virtual and Experimental results

The bulge test is mostly used to analyze equibiaxial tensile stress state at the pole of inflated isotropic membranes. Three-dimensional digital image correlation (3D-DIC) technique allows the determination of three-dimensional surface displacements and strain fields. In this paper, a method is proposed to determine also the membrane stress tensor fields for in-plane isotropic materials, independently of any constitutive equation. Stress-strain state is then known at any surface point which enriches greatly experimental data deduced from the axisymmetric bulge tests. Our method consists, first in calculating from the 3D-DIC experimental data the membrane curvature tensor at each surface point of the bulge specimen. Then, curvature tensor fields are used to investigate axisymmetry of the test. Finally in the axisymmetric case, membrane stress tensor fields are determined from meridional and circumferential curvatures combined with the measurement of the inflating pressure. Our method is first validated for virtual 3D-DIC data, obtained by numerical simulation of a bulge test using a hyperelastic material model. Afterward, the method is applied to an experimental bulge test performed using as material a silicone elastomer. The stress-strain fields which are obtained using the proposed method are compared with results of the finite element simulation of this overall bulge test using a neo-Hookean model fitted on uniaxial and equibiaxial tensile tests.     

Stochastic framework for modeling the linear apparent behavior of complex materials: Application to random porous materials with interphases

This paper is concerned with the modeling of randomness in multiscale analysis of heterogeneous materials. More specifically, a framework dedicated to the stochastic modeling of random properties is first introduced. A probabilistic model for matrix-valued second-order random fields with symmetry propertries, recently proposed in the literature, is further reviewed. Algorithms adapted to the Monte Carlo simulation of the proposed representation are also provided. The derivations and calibration procedure are finally exemplified through the modeling of the apparent properties associated with an elastic porous microstructure containing stochastic interphases.     

基于微面有效应力矢量的各向异性屈服准则

基于微面模型,定义损伤变量为微面上有效承载面积的减少. 将Kachanov的一维有效 应力概念推广到三维,提出微面有效应力矢量的概念. 根据微面的有效应力矢量,将无损材 料的宏观应力张量及不变量与微面应力矢量的积分关系拓展到有损材料,得到了有损材料的 宏观有效应力张量及其不变量与宏观名义应力张量、微面面积损伤组构张量之间的关系. 将 无损材料的以应力张量不变量表示的Drucker-Prager准则推广到有损材料,建立了含缺陷 材料的各向异性屈服准则. 对有损材料,宏观有效应力张量与Murakami的有效应力张量具 有相同的形式,各向异性强度准则与Liu等提出的扩展Hill准则有相同的形式,当不考虑 静水应力对屈服的影响时,它与Hill准则具有相同的形式.     

Invariant formulation of hyperelastic transverse isotropy based on polyconvex free energy functions

In this paper we propose a formulation of polyconvex anisotropic hyperelasticity at finite strains. The main goal is the representation of the governing constitutive equations within the framework of the invariant theory which automatically fulfill the polyconvexity condition in the sense of Ball in order to guarantee the existence of minimizers. Based on the introduction of additional argument tensors, the so-called structural tensors, the free energies and the anisotropic stress response functions are represented by scalar-valued and tensor-valued isotropic tensor functions, respectively. In order to obtain various free energies to model specific problems which permit the matching of data stemming from experiments, we assume an additive structure. A variety of isotropic and anisotropic functions for transversely isotropic material behaviour are derived, where each individual term fulfills a priori the polyconvexity condition. The tensor generators for the stresses and moduli are evaluated in detail and some representative numerical examples are presented. Furthermore, we propose an extension to orthotropic symmetry.     

Non‐linear flux‐splitting schemes with imposed discrete maximum principle for elliptic equations with highly anisotropic coefficients

Families of flux‐continuous, locally conservative, finite‐volume schemes have been developed for solving the general tensor pressure equation of petroleum reservoir simulation on structured and unstructured grids. The schemes are applicable to diagonal and full tensor pressure equation with generally discontinuous coefficients and remove the O(1) errors introduced by standard reservoir simulation schemes when applied to full tensor flow approximation. The family of flux‐continuous schemes is quantified by a quadrature parameterization. Improved convergence using the quadrature parameterization has been established for the family of flux‐continuous schemes. When applied to strongly anisotropic full‐tensor permeability fields the schemes can fail to satisfy a maximum principle (as with other FEM and finite‐volume methods) and result in spurious oscillations in the numerical pressure solution. This paper presents new non‐linear flux‐splitting techniques that are designed to compute solutions that are free of spurious oscillations. Results are presented for a series of test‐cases with strong full‐tensor anisotropy ratios. In all cases the non‐linear flux‐splitting methods yield pressure solutions that are free of spurious oscillations. Copyright © 2010 John Wiley & Sons, Ltd.     

Stochastic modeling and identification of an uncertain computational dynamical model with random fields properties and model uncertainties

This paper is devoted to the construction and to the identification of a probabilistic model of random fields in the presence of modeling errors, in high stochastic dimension and presented in the context of computational structural dynamics. Due to the high stochastic dimension of the random quantities which have to be identified using statistical inverse methods (challenging problem), a complete methodology is proposed and validated. The parametric–nonparametric (generalized) probabilistic approach of uncertainties is used to perform the prior stochastic models: (1) system-parameters uncertainties induced by the variabilities of the material properties are described by random fields for which their statistical reductions are still in high stochastic dimension and (2) model uncertainties induced by the modeling errors are taken into account with the nonparametric probabilistic approach in high stochastic dimension. For these two sources of uncertainties, the methodology consists in introducing prior stochastic models described with a small number of parameters which are simultaneously identified using the maximum likelihood method and experimental responses. The steps of the methodology are explained and illustrated through an application.     

Footing Settlement on a Consolidating Soil Layer with Stochastic Properties

Geotechnical engineering applications are characterized by various sources of uncertainties, most of them attributed to the stochastic nature of soil parameters and their properties. In particular, soil’s inherent random heterogeneity, inexact measurements and insufficient data necessitate numerical methods that incorporate the stochastic soil properties for a realistic representation of the soil behavior. In this paper, the process of consolidation of saturated soils is examined on the basis of the coupled u–p finite element formulation. A generalized Newmark implicit time integration scheme is implemented to treat the time integration of the coupled consolidation equations. A benchmark geotechnical engineering problem of a strip footing resting on a saturated soil layer is analyzed. The soil permeability coefficient k, as well as the elastic modulus E, are treated as lognormal random fields in two dimensions. The investigation of the effect of the spatial variability of the soil properties on the response of a footing–soil system is examined by means of the direct Monte Carlo simulation. The influence of the coefficient of variation and correlation length of the stochastic fields is quantified in terms of footing settlements, as well as excess soil water pore pressure. The effects of spatial variability of the permeability coefficient k and the elastic modulus E on the system response are demonstrated. It is shown that the footing differential settlement, along with generated excess pore pressures, is highly affected by the variation of the soil properties considered, as well as the correlation length of the underlying random fields.     

Continuum modeling of twinning,amorphization, and fracture: theory and numerical simulations

A continuum mechanical theory is used to model physical mechanisms of twinning, solid-solid phase transformations, and failure by cavitation and shear fracture. Such a sequence of mechanisms has been observed in atomic simulations and/or experiments on the ceramic boron carbide. In the present modeling approach, geometric quantities such as the metric tensor and connection coefficients can depend on one or more director vectors, also called internal state vectors. After development of the general nonlinear theory, a first problem class considers simple shear deformation of a single crystal of this material. For homogeneous fields or stress-free states, algebraic systems or ordinary differential equations are obtained that can be solved by numerical iteration. Results are in general agreement with atomic simulation, without introduction of fitted parameters. The second class of problems addresses the more complex mechanics of heterogeneous deformation and stress states involved in deformation and failure of polycrystals. Finite element calculations, in which individual grains in a three-dimensional polycrystal are fully resolved, invoke a partially linearized version of the theory. Results provide new insight into effects of crystal morphology, activity or inactivity of different inelasticity mechanisms, and imposed deformation histories on strength and failure of the aggregate under compression and shear. The importance of incorporation of inelastic shear deformation in realistic models of amorphization of boron carbide is noted, as is a greater reduction in overall strength of polycrystals containing one or a few dominant flaws rather than many diffusely distributed microcracks.     

On the thermodynamics of thermoelastic materials with additional scalar degrees of freedom

This work deals with the thermodynamic formulation of a model for a class of materials containing microstructure which evolves or changes relative to\/ the (global) bulk material. The approach taken here is based on a generalization of the total energy, total energy flux, and total energy supply to take into account the corresponding additional degrees of freedom involved. Restricting attention for simplicity to thermoelastic materials with scalar-valued such degrees of freedom, the thermodynamically-consistent forms of the remaining balance and corresponding constitutive relations for this material class are obtained in the context of the Müller-Liu entropy principle. In particular, the thermodynamically-consistent form of the evolution relation for the additional scalar-valued degrees of freedom obtained in this fashion contains in part\/ the well-known generalized Euler-Lagrange or Ginzburg-Landau relations established in many previous work, with the remaining terms accounting for effects associated with microinertia, or with non-equilibrium processes. Received September 14, 1998     

A Generalized Theory of Stress and Strain Measures in the Classical Continuum Mechanics

A generalized theory of stress and strain tensor measures in the classical continuum mechanics is discussed: the main axioms of the theory are proposed, the general formulas for new tensor measures are derived, arid an energy conjugate theorem is formulated to distinguish the complete Lagrangian class of measures. As a subclass, a simple Lagrangian class of energy conjugate measures of stresses and finite strains is constructed in which the families of holonomic and corotational measures are distinguished. The characteristics of holonomic and corotational measures are studied by comparing the tensor measures of the simple Lagrangian class with one another and with logarithmic measures. For the simple Lagrangian class and its families, their completeness and closure are shown with respect to the choice of a generating pair of energetically conjugate measures. The applications of the new tensor measures in modeling the properties of plasticity, viscoelasticity, and shape memory are mentioned.     

小样本数据下圆柱薄壳初始缺陷不确定性量化的极大熵方法

具有不确定性特征的初始缺陷被认为是导致薄壳结构实际临界载荷值与理论解不相符并呈现分散特征的主要原因.对实际薄壳结构初始缺陷的建模至少需要考虑两个方面的不确定性量化,一是对缺陷分布形式和幅值等固有随机性的量化,二是对小样本量和不准确测量所导致缺陷统计量的不确定性的量化.本文在利用随机场的Karhunen-Loeve展开法对薄壳初始几何缺陷建模的基础上,提出一种基于极大熵原理的缺陷建模方法.首先,采用极大熵分布来估计Karhunen-Loeve随机变量的概率密度函数,以适应不能使用高斯随机场进行缺陷随机场建模的情况.随后,通过将经典的等式约束极大熵模型扩展为区间约束极大熵模型,实现对实际工程中仅能获得少量薄壳结构几何缺陷样本数据所导致的认知不确定性的量化.最后,将所提方法用于对国际缺陷数据库的A-Shell进行缺陷建模和临界载荷预测.研究表明,所提基于区间约束极大熵原理的随机场建模方法在能够有效表征实测数据高阶矩信息的同时,还具备量化小样本数据导致的认知不确定性的能力,并且高斯随机场模型和基于等式约束极大熵原理的随机场模型是本文所提建模方法的两种特殊情况.     

On Non-Symmetric Deformations of an Incompressible Nonlinearly Elastic Isotropic Sphere

A class of non-symmetric deformations of a neo-Hookean incompressible nonlinearly elastic sphere are investigated. It is found via the semi-inverse method that, to satisfy the governing three-dimensional equations of equilibrium and the incompressibility constraint, only three special cases of the class of deformation fields are possible. One of these is the trivial solution, one the solution describing radially symmetric deformation, and the other a (non-symmetric, non-homogeneous) deformation describing inflation and stretching. The implications of these results for cavitation phenomena are also discussed. In the course of this work, we also present explicitly the spherical polar coordinate form of the equilibrium equations for the nominal stress tensor for a general hyperelastic solid. These are more complicated than their counterparts for Cauchy stresses due to the mixed bases (both reference and deformed) associated with the nominal (or Piola-Kirchhoff) stress tensor, but more useful in considering general deformation fields. This revised version was published online in August 2006 with corrections to the Cover Date.     

Dielectric elastomer composites: A general closed-form solution in the small-deformation limit

A solution for the overall electromechanical response of two-phase dielectric elastomer composites with (random or periodic) particulate microstructures is derived in the classical limit of small deformations and moderate electric fields. In this limit, the overall electromechanical response is characterized by three effective tensors: a fourth-order tensor describing the elasticity of the material, a second-order tensor describing its permittivity, and a fourth-order tensor describing its electrostrictive response. Closed-form formulas are derived for these effective tensors directly in terms of the corresponding tensors describing the electromechanical response of the underlying matrix and the particles, and the one- and two-point correlation functions describing the microstructure. This is accomplished by specializing a new iterative homogenization theory in finite electroelastostatics (Lopez-Pamies, 2014) to the case of elastic dielectrics with even coupling between the mechanical and electric fields and, subsequently, carrying out the pertinent asymptotic analysis.Additionally, with the aim of gaining physical insight into the proposed solution and shedding light on recently reported experiments, specific results are examined and compared with an available analytical solution and with new full-field simulations for the special case of dielectric elastomers filled with isotropic distributions of spherical particles with various elastic dielectric properties, including stiff high-permittivity particles, liquid-like high-permittivity particles, and vacuous pores.     

Controllability of viscoelastic stresses for nonlinear Maxwell models

We investigate a class of models for viscoelastic fluids, in which the elastic stress is determined by a conformation tensor, and the conformation tensor is linked to the velocity field by a system of ordinary differential equations. We study the question which values of the conformation tensor can be reached in a homogeneous flow, subject to a given initial condition and arbitrary velocity fields. This problem is a special “easy” case for the question of controllability of viscoelastic flows. For a class of models, we show that constraints on the values of the conformation tensor are given by lower and/or upper bounds on its determinant. The behavior of seemingly similar models, e.g. the PTT, Giesekus and Peterlin dumbbell models, turns out to be surprisingly different.     

基于球形颗粒几何排列的离散元试样高效生成方法

在球体离散元数值模拟中,颗粒的初始排列状态是影响计算效率和计算结果的重要环节。本文采用前进面几何构造算法,提出了一种基于网格搜索的球形颗粒随机排列高效算法。通过求解空间三边方程,满足了粒径设置的任意大小的颗粒依次置入前进面的外侧,并与构成前进面的三个颗粒相互接触。为获得高体积分数的颗粒簇,该算法允许颗粒改变其粒径大小。采用颗粒网格化方法可以简化前进面的搜索,并由此提高排列效率。通过计算平均配位数、体积分数和二阶结构张量的特征值,对不同粒径比下得到的立方体试样进行了分析,得到试样配位数及体积分数均随着粒径比的增大而增大,且得到的试样为各向同性。此外,空间网格的大小和初始颗粒的生成点对随机排列的效率均会产生显著的影响。最后,对非规则铁路道砟进行了精细构造及压碎模拟,发现DEM模拟得到的应力-应变曲线与试验结果基本吻合,验证了该算法得到的颗粒试样在模拟道砟裂纹起裂、扩展等过程的有效性。