The Representation Theorem for Linear, Isotropic Tensor Functions in Even Dimensions

In this note we give a proof of the representation theorem for linear, isotropic, tensor functions, which only assumes invariance under proper orthogonal tensors. This revised version was published online in August 2006 with corrections to the Cover Date.     

Objective tensorial representation of the pressure–strain correlations of turbulence

An explicit algebraic model for the fluctuating pressure–strain rate correlations of turbulence is developed by the use of representation theorems for tensor-valued isotropic tensors, and by invoking the principle of objectivity. The resulting model differs from others by the absence of the vorticity tensor from its formulation. The new model is calibrated by reference to data from homogeneous shear flows, and its potential as a practical tool for the analysis of turbulent flows is demonstrated by numerical simulations of a benchmark two-dimensional shear layer.     

Objective Tensor Rates and Applications in Formulation of Hyperelastic Relations

Following Ogden, a class of objective (Lagrangian and Eulerian) tensors is identified among the second-rank tensors characterizing continuum deformation, but a more general definition of objectivity than that used by Ogden is introduced. Time rates of tensors are determined using convective rates. Sufficient conditions of objectivity are obtained for convective rates of objective tensors. Objective convective rates of strain tensors are used to introduce pairs of symmetric stress and strain tensors conjugate in a generalized sense. The classical definitions of conjugate Lagrangian (after Hill) and Eulerian (after Xiao et al.) stress and strain tensors are particular cases of the definition of conjugacy of stress and strain tensors in the generalized sense used in the present paper. Pairs of objective stress and strain tensors conjugate in the generalized sense are used to formulate constitutive relations for a hyperelastic medium. A family of objective generalized strain tensors is introduced, which is broader than Hill’s family of strain tensors. The basic forms of the hyperelastic constitutive relations are obtained with the aid of pairs of Lagrangian stress and strain tensors conjugate after Hill (the strain tensors in these pairs belong to the family of generalized strain tensors). A method is presented for generating reduced forms of the constitutive relations with the aid of pairs of Lagrangian and Eulerian stress and strain tensors conjugate in the generalized sense which are obtained from pairs of Lagrangian tensors conjugate after Hill by mapping tensor fields on one configuration of a deformable body to tensor fields on another configuration.      

Orientation tensors in simple flows of dilute suspensions of non-Brownian rigid ellipsoids,comparison of analytical and approximate solutions

General analytical solutions are obtained for the planar orientation structure of rigid ellipsoid of revolutions subjected to an arbitrary homogeneous flow in a Newtonian fluid. Both finite and infinite aspect ratio particles are considered. The orientation structure is described in terms of two-dimensional, time-dependent tensors that are commonly employed in constitutive equations for anisotropic fluids such as fiber suspensions. The effect of particle aspect ratio on the evolution of orientation structure is studied in simple shear and planar elongational flows. With the availability of analytical solutions, accuracies of quadratic closure approximations used for nonhomogeneous flows are analyzed, avoiding numerical integration of orientation distribution function. In general, fourth-order orientation evolution equations with sixth-order quadratic closure approximations yield more accurate representations compared to the commonly used second-order evolution equations with fourth-order quadratic closure approximations. However, quadratic closure approximations of any order are found to give correct maximum orientation angle (i.e., preferred direction) results for all particle aspect ratios and flow cases.     

Relationship of crack fabric tensors of different orders

A lower-order crack fabric tensor can be determined accurately by a higher-order one, but in general the reverse does not hold. In this paper, the approximate dependence of a higher-order fabric tensor on a lower-order one is established based on the properties of orientation distribution functions (ODF). As a demonstration of its application, the approximate relationship is used to simplify the fabric-tensor dependent compliance increment due to the presence of the cracks.     

On the Explicit Determination of the Polar Decomposition in n-Dimensional Vector Spaces

A method for the explicit determination of the polar decomposition (and the related problem of finding tensor square roots) when the underlying vector space dimension n is arbitrary (but finite), is proposed. The method uses the spectral resolution, and avoids the determination of eigenvectors when the tensor is invertible. For any given dimension n, an appropriately constructed van der Monde matrix is shown to play a key role in the construction of each of the component matrices (and their inverses) in the polar decomposition.     

Complete Asymptotic Expansions of Solutions of the System of Elastostatics in the Presence of an Inclusion of Small Diameter and Detection of an Inclusion

We consider the system of elastostatics for an elastic medium consisting of an imperfection of small diameter, embedded in a homogeneous reference medium. The Lamé constants of the imperfection are different from those of the background medium. We establish a complete asymptotic formula for the displacement vector in terms of the reference Lamé constants, the location of the imperfection and its geometry. Our derivation is rigorous, and based on layer potential techniques. The asymptotic expansions in this paper are valid for an elastic imperfection with Lipschitz boundaries. In the course of derivation of the asymptotic formula, we introduce the concept of (generalized) elastic moment tensors (Pólya–Szegö tensor) and prove that the first order elastic moment tensor is symmetric and positive (negative)-definite. We also obtain estimation of its eigenvalue. We then apply these asymptotic formulas for the purpose of identifying with high precision the order of magnitude of the diameter of the elastic inclusion, its location, and its elastic moment tensors.     

On the Comparison of Two Strategies to Formulate Orthotropic Hyperelasticity

The main goal of this work is to clarify the relation between two strategies to formulate constitutive equations for orthotropic materials at large strains. On the one hand, the classical approach is based on the incorporation of structural tensors into the free energy function via an enriched set of invariants. On the other hand, a fictitious isotropic configuration is introduced which renders an anisotropic, undeformed reference configuration via an appropriate linear tangent map. This formulation results in a reduced (with respect to the more general setting based on structural tensors) but nevertheless physically motivated set of invariants which are related to the invariants defined by structural tensors. As a main conceptual advantage standard isotropic constitutive equations can be applied and moreover, due to the reduced set of physically motivated invariants, the numerical treatment within a finite element setting becomes manageable.     

Evolution of fabric tensors in damage mechanics of solids with micro-cracks: Part II – Evolution of length and orientation of micro-cracks with an application to uniaxial tension

In Part II of this work, the equations of thermodynamics are employed in order to derive the exact evolution equations of the fabric tensors defined in Part I (companion paper). In this regard, a thermodynamic force that is associated with the fabric tensor is defined and utilized in the derivation of the evolution equations. A special case of uniaxial tension is solved in order to illustrate the theory.We also derive specific uncoupled equations for the evolution of the length and orientation of micro-cracks. In this regard, some interesting results are obtained. It is concluded that the micro-crack length and orientation cannot evolve simultaneously for the same set of micro-cracks. However, two different sets of micro-cracks may be considered in the same RVE where in one set the micro-crack length evolves, while in the second set the micro-crack orientation evolves.     

Evolution of fabric tensors in damage mechanics of solids with micro-cracks: Part I – Theory and fundamental concepts

The evolution of fabric tensors based on micro-crack distributions is formulated based on sound thermodynamic principles. In Part I of this work, the exact definition of fabric tensors based on micro-crack distributions is presented. This definition is seen to incorporate both the orientation and length of a micro-crack. In this regard, the micro-crack distribution is assumed to be radially symmetric, i.e. symmetric about a line through the origin.The equations of thermodynamics are employed in order to derive the exact evolution equations of the fabric tensors defined in the first part. In this regard, a thermodynamic force that is associated with the fabric tensor is defined and utilized in the derivation of the evolution equations. The application of the theory to the case of uniaxial tension is derived in Part II (companion paper) of this work.     

Effective Diffusivity Tensors of Point-Like Molecules in Anisotropic Porous Media by Monte Carlo Simulation

Monte Carlo simulations of random walks in anisotropic structured media are performed to determine the dependence of effective diffusivities on geometrical properties. The anisotropic media used in this study are periodic systems, which are generated by extending primitive, face-centered, and body-centered unit cells indefinitely in all axial directions. Results of simulations compare well with published experimental data and the calculations by the volume averaging method. In addition, these results suggest that if the 2D media with percolation thresholds subtantially differ from those of 3D, 2D approximations of 3D media are not satisfactory. When percolation thresholds are the same, the effective diffusivity tensors depend solely on the porosity. This fact has been suggested for isotropic media and it seems to hold for anisotropic media.     

The effective properties of piezocomposites, part II: The effective electroelastic moduli

Since piezoelectric ceramic/polymer composites have been widely used as smart materials and smart structures, it is more and more important to obtain the closed-from solutions of the effective properties of piezocomposites with piezoelectric ellipsoidal inclusions. Based on the closed-from solutions of the electroelastic Eshelby's tensors obtained in the part I of this paper and the generalized Budiansky's energy-equivalence framework, the closed-form general relations of effective electroelastic moduli of the piezocomposites with piezoelectric ellipsoidal inclusions are given. The relations can be applicable for several micromechanics models, such as the dilute solution and the Mori-Tanaka's method. The difference among the various models is shown to be the way in which the average strain and the average electric field of the inclusion phase are evaluated. Comparison between predicted and experimental results shows that the theoretical values in this paper agree quite well with the experimental results. These expression can be readily utilized in analysis and design of piezocomposites. The project supported by the National Natural Science Foundation of China     

Invariant Properties for Finding Distance in Space of?Elasticity Tensors

Using orthogonal projections, we investigate distance of a given elasticity tensor to classes of elasticity tensors exhibiting particular material symmetries. These projections depend on the orientation of the elasticity tensor; hence the distance is obtained as the minimization of corresponding expressions with respect to the action of the orthogonal group. These expressions are stated in terms of the eigenvalues of both the given tensor and the projected one. The process of minimization is facilitated by the fact that, as we prove, the traces of the corresponding Voigt and dilatation tensors are invariant under these orthogonal projections. For isotropy, cubic symmetry and transverse isotropy, we formulate algorithms to find both the orientation and the eigenvalues of the elasticity tensor endowed with a particular symmetry and closest to the given elasticity tensor.      

Functional relation between two symmetric second-rank tensors

The functional relationship between two symmetric second-rank tensors is considered. A new interpretation of the components of the tensors as projections onto an orthogonal tensor basis is given. It is shown that the constitutive relations can be written in the form of six functions each of which depends on one variable. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 5, pp. 134–137, September–October, 2007.     

Determination of permeability tensors for two-phase flow in homogeneous porous media: Theory

In this paper we continue previous studies of the closure problem for two-phase flow in homogeneous porous media, and we show how the closure problem can be transformed to a pair of Stokes-like boundary-value problems in terms of pressures that have units of length and velocities that have units of length squared. These are essentially geometrical boundary value problems that are used to calculate the four permeability tensors that appear in the volume averaged Stokes' equations. To determine the geometry associated with the closure problem, one needs to solve the physical problem; however, the closure problem can be solved using the same algorithm used to solve the physical problem, thus the entire procedure can be accomplished with a single numerical code.Nomenclature a a vector that maps V onto , m^-1. - A a tensor that maps V onto . - A area of the - interface contained within the macroscopic region, m². - A area of the -phase entrances and exits contained within the macroscopic region, m². - A area of the - interface contained within the averaging volume, m². - A area of the -phase entrances and exits contained within the averaging volume, m². - Bo Bond number (= (=(–)g²/). - Ca capillary number (= v/). - g gravitational acceleration, m/s². - H mean curvature, m^-1. - I unit tensor. - permeability tensor for the -phase, m². - viscous drag tensor that maps V onto V. - ^* dominant permeability tensor that maps onto v , m². - ^* coupling permeability tensor that maps onto v , m². - characteristic length scale for the -phase, m. - l characteristic length scale representing both and , m. - L characteristic length scale for volume averaged quantities, m. - n unit normal vector directed from the -phase toward the -phase. - n unit normal vector representing both n and n . - n unit normal vector representing both n and n . - P pressure in the -phase, N/m². - p superficial average pressure in the -phase, N/m². - p intrinsic average pressure in the -phase, N/m². - p –p , spatial deviation pressure for the -phase, N/m². - r ₀ radius of the averaging volume, m. - r position vector, m. - t time, s. - v fluid velocity in the -phase, m/s. - v superficial average velocity in the -phase, m/s. - v intrinsic average velocity in the -phase, m/s. - v –v , spatial deviation velocity in the -phase, m/s. - V volume of the -phase contained within the averaging volmue, m³. - averaging volume, m³. Greek Symbols V /, volume fraction of the -phase. - viscosity of the -phase, Ns/m². - density of the -phase, kg/m³. - surface tension, N/m. - (v +v ^T ), viscous stress tensor for the -phase, N/m².     

On the modelling of anisotropic elastic and inelastic material behaviour at large deformation

The purpose of this work is the formulation and discussion of an approach to the modelling of anisotropic elastic and inelastic material behaviour at large deformation. This is done in the framework of a thermodynamic, internal-variable-based formulation for such a behaviour. In particular, the formulation pursued here is based on a model for plastic or inelastic deformation as a transformation of local reference configuration for each material element. This represents a slight generalization of its modelling as an elastic material isomorphism pursued in earlier work, allowing one in particular to incorporate the effects of isotropic continuum damage directly into the formulation. As for the remaining deformation- and stress-like internal variables of the formulation, these are modelled in a fashion formally analogous to so-called structure tensors. On this basis, it is shown in particular that, while neither the Mandel nor back stress is generally so, the stress measure thermodynamically conjugate to the plastic “velocity gradient”, containing the difference of these two stress measures, is always symmetric with respect to the Euclidean metric, i.e., even in the case of classical or induced anisotropic elastic or inelastic material behaviour. Further, in the context of the assumption that the intermediate configuration is materially uniform, it is shown that the stress measure thermodynamically conjugate to the plastic velocity gradient is directly related to the Eshelby stress. Finally, the approach is applied to the formulation of metal plasticity with isotropic kinematic hardening.     

Continuum representation of systems of dislocation lines: A general method for deriving closed-form evolution equations

Plasticity is governed by the evolution of, in general anisotropic, systems of dislocations. We seek to faithfully represent this evolution in terms of density-like variables which average over the discrete dislocation microstructure. Starting from T. Hochrainer's continuum theory of dislocations (CDD) (Hochrainer, 2015), we introduce a methodology based on the ‘Maximum Information Entropy Principle’ (MIEP) for deriving closed-form evolution equations for dislocation density measures of different order. These equations provide an optimum representation of the kinematic properties of systems of curved and connected dislocation lines with the information contained in a given set of density measures. The performance of the derived equations is benchmarked against other models proposed in the literature, using discrete dislocation dynamics simulations as a reference. As a benchmark problem we study dislocations moving in a highly heterogeneous, persistent-slip-band like geometry. We demonstrate that excellent agreement with discrete simulations can be obtained in terms of a very small number of averaged dislocation fields containing information about the edge and screw components of the total and excess (geometrically necessary) dislocation densities. From these the full dislocation orientation distribution which emerges as dislocations move through a channel-wall structure can be faithfully reconstructed.     

Separation of screw-sensed particles in a homogeneous shear field

The Stokes motions of three-dimensional screw-sensed slender particles in a homogeneous shear field are investigated, including the effects of buoyancy. Conclusions are drawn about the possibility of achieving a separation of mixtures of right- and left-handed particles. The linearity of the Stokes equations allows complex flows to be solved by adding the effects of the several terms which describe the flow in which the particle is immersed. The homogeneous shear flow considered here consists of three such terms; solutions for a series of 12 unit motions are sufficient to determine the hydrodynamic resistance tensors. The forces and torques experienced by screw-sensed particles are calculated from these 51 resistance tensors, using slender-filament theory. The results allow an estimate of the range of buoyancy parameters for which gravitational sedimentation can be neglected. The fundamental component of the particle motion is a rotation, at approximately the same angular velocity as that of the fluid. Superimposed on this are variations, of large period, in the particle orientation. A phase plane analysis is used to find the terminal orientations. Very long calculation times are required for the phase portrait. An approximate method based on azimuthally-averaged equations is developed to avoid the requirements for long time integration.     

Effective properties of a piezocomposite containing shape memory alloy and inert inclusions

Received July 2, 2001 / Published online February 4, 2002     

Pure Shear – A Footnote

A result on pure shear provides the motivation for the determination of some new general results relating real second order Cartesian tensors.