Stochastic Modelling of Symmetric Positive Definite Material Tensors

Sharana Kumar Shivanand,Bojana Rosić,Hermann G. Matthies
2024-02-09
Abstract:Spatial symmetries and invariances play an important role in the behaviour of materials and should be respected in the description and modelling of material properties. The focus here is the class of physically symmetric and positive definite tensors, as they appear often in the description of materials, and one wants to be able to prescribe certain classes of spatial symmetries and invariances for each member of the whole ensemble, while at the same time demanding that the mean or expected value of the ensemble be subject to a possibly 'higher' spatial invariance class. We formulate a modelling framework which not only respects these two requirements$-$positive definiteness and invariance$-$but also allows a fine control over orientation on one hand, and strength/size on the other. As the set of positive definite tensors is not a linear space, but rather an open convex cone in the linear space of physically symmetric tensors, we consider it advantageous to widen the notion of mean to the so-called Fréchet mean on a metric space, which is based on distance measures or metrics between positive definite tensors other than the usual Euclidean one. It is shown how the random ensemble can be modelled and generated, independently in its scaling and orientational or directional aspects, with a Lie algebra representation via a memoryless transformation. The parameters which describe the elements in this Lie algebra are then to be considered as random fields on the domain of interest. As an example, a 2D and a 3D model of steady-state heat conduction in a human proximal femur, a bone with high material anisotropy, is modelled with a random thermal conductivity tensor, and the numerical results show the distinct impact of incorporating into the constitutive model different material uncertainties$-$scaling, orientation, and prescribed material symmetry$-$on the desired quantities of interest.
Numerical Analysis,Computational Engineering, Finance, and Science
What problem does this paper attempt to address?
The paper is primarily dedicated to addressing the issue of how to reasonably represent a set of tensors with spatial symmetry and positive definiteness when describing material properties. Specifically, the research focuses on physically symmetric and positive definite (SPD) second-order tensors, which frequently appear in the description of material behavior. The goal of the paper is to precisely control the mean or expected value of the entire set while maintaining positive definiteness and specific spatial symmetry or invariance, allowing it to conform to possibly higher-order symmetry categories. To achieve this goal, the authors propose a modeling framework that not only respects the requirements of positive definiteness and invariance but also allows for fine control over orientation and intensity/magnitude. Since the set of positive definite tensors is not a linear space but an open convex cone in the linear space of physically symmetric tensors, the paper considers extending the concept of the mean to the so-called Fréchet mean, which is based on the distance metric or metric between positive definite tensors rather than the usual Euclidean metric. The paper discusses in detail how to model and generate random sets independently of scaling and directional aspects, using Lie algebra representations through memoryless transformations. Additionally, the paper discusses how to treat the parameters describing these Lie algebra elements as random fields over the domain of interest. As an example, the paper presents 2D and 3D models of steady-state heat conduction in the proximal femur of the human body, where the femur exhibits high material anisotropy. Its random thermal conductivity tensor is used for simulation, and numerical results show the significant impact of incorporating different material uncertainties—scaling, orientation, and prescribed material symmetry—into the configuration model on the quantities of interest. This provides a new perspective for understanding material uncertainty and quantifying uncertainty.