Stationarity of the strain energy density for some classes of anisotropic solids

Homogeneous, anisotropic and linearly elastic solids, subjected to a given state of strain (or stress), are considered. The problem dealt with consists in finding the mutual orientations of the principal directions of strain to the material symmetry axes in order to make the strain energy density stationary. Such relative orientations are described through three Euler’s angles. When the stationarity problem is formulated for the generally anisotropic solid, it is shown that the necessary condition for stationarity demands for coaxiality of the stress and the strain tensors. From this feature, a procedure which leads to closed form solutions is proposed. To this end, tetragonal and cubic symmetry classes, together with transverse isotropy, are carefully dealt with, and for each case all the sets of Euler’s angles corresponding to critical points of the energy density are found and discussed. For these symmetries, three material parameters are then defined, which play a crucial role in ordering the energy values corresponding to each solution.     

Nonlinear Eulerian thermoelasticity for anisotropic crystals

A complete continuum thermoelastic theory for large deformation of crystals of arbitrary symmetry is developed. The theory incorporates as a fundamental state variable in the thermodynamic potentials what is termed an Eulerian strain tensor (in material coordinates) constructed from the inverse of the deformation gradient. Thermodynamic identities and relationships among Eulerian and the usual Lagrangian material coefficients are derived, significantly extending previous literature that focused on materials with cubic or hexagonal symmetry and hydrostatic loading conditions. Analytical solutions for homogeneous deformations of ideal cubic crystals are studied over a prescribed range of elastic coefficients; stress states and intrinsic stability measures are compared. For realistic coefficients, Eulerian theory is shown to predict more physically realistic behavior than Lagrangian theory under large compression and shear. Analytical solutions for shock compression of anisotropic single crystals are derived for internal energy functions quartic in Lagrangian or Eulerian strain and linear in entropy; results are analyzed for quartz, sapphire, and diamond. When elastic constants of up to order four are included, both Lagrangian and Eulerian theories are capable of matching Hugoniot data. When only the second-order elastic constant is known, an alternative theory incorporating a mixed Eulerian–Lagrangian strain tensor provides a reasonable approximation of experimental data.     

Wrinkling of Orthotropic Membranes: An Analysis by the Polar Method

In the present paper, a simple membrane model based on the wrinkle strain approach is revisited with the aim of examining how the material elastic constants affect the static response of anisotropic membranes when wrinkling is taken into account. Employing the polar method, we analyze the role played by the polar moduli, which enable expressing the elasticity matrix components of an anisotropic material in terms of its invariant quantities. With reference to orthotropic materials, we first address the issue of membrane susceptibility to wrinkling by investigating the influence of the three polar parameters characterizing the anisotropic part of the constitutive law. The stress and strain states at any given point in a wrinkled membrane are analyzed by searching for explicit expressions for the principal values of stress and wrinkle strain. Finally, a comparison between our results and those obtained by a numerical solution available in the literature is made in the basic case of a membrane subjected to a pure shear strain state.     

Optimal material orientation of linear and non-linear elastic 3D anisotropic materials

Optimal material orientation problems of linear and non-linear elastic three-dimensional anisotropic materials are studied. Most commonly, the energy based formulation is applied for solving orientational design problems of anisotropic materials, considering elastic energy density as a measure of the stress strain state. The same approach is used in the current study, but the strength criteria based approaches are also discussed. A simple relation between the stationary conditions in terms of Euler angles and the optimality conditions in terms of strains is pointed out. The complexity analysis of the different existing optimality conditions has been performed. The solution of the posed optimization problem is decomposed into the strain level solution, search for global extremes and evaluation of Euler angles (parameters). The results obtained are extended to some nonlinear elastic material models.     

Remarks on Coaxiality of Strain and Stress in Anisotropic Elasticity

Two different consequences of the problem of the extremization of the strain energy by varying the orientation of the material symmetry axes relative to the principal axes of stress are discussed. These two different consequences depend upon whether the stress state is considered as arbitrary and general or as fixed and specific. This revised version was published online in August 2006 with corrections to the Cover Date.     

Extreme conditions of elastic constants and principal axes of anisotropy

This paper describes the derivation of extreme conditions of each elasticity coefficient (Young’s modulus, shear modulus, et al.,) for the general case of linear-elastic anisotropic materials. The stationarity conditions are obtained, and they determine the orthogonal coordinate systems being the principal axes of anisotropy, where the number of independent elasticity constants decreases from 21 to 18 and, in some cases of anisotropy, to 15 or lower. The example of a material with cubic symmetry is given.     

On the modelling of anisotropic material behaviour in viscoplasticity

A uniaxial viscoplastic deformation is motivated as a discrete sequence of stable and unstable equilibrium states and approximated by a smooth family of stable states of equilibrium depending on the history of the mechanical process. Three-dimensional crystal viscoplasticity starts from the assumption that inelastic shearings take place on slip systems, which are known from the particular geometric structure of the crystal. A constitutive model for the behaviour of a single crystal is developed, based on a free energy, which decomposes into an elastic and an inelastic part. The elastic part, the isothermal strain energy, depends on the elastic Green strain and allows for the initial anisotropy, known from the special type of the crystal lattice. Additionally, the strain energy function contains an orthogonal tensor-valued internal variable representing the orientation of the anisotropy axes. This orientation develops according to an evolution equation, which satisfies the postulate of full invariance in the sense that it is an observer-invariant relation. The inelastic part of the free energy is a quadratic function of the integrated shear rates and corresponding internal variables being equivalent to backstresses in order to consider kinematic hardening phenomena on the slip system level. The evolution equations for the shears, backstresses and crystallographic orientations are thermomechanically consistent in the sense that they are compatible with the entropy inequality. While the general theory applies to all types of lattices, specific test calculations refer to cubic symmetry (fcc) and small elastic strains. The simulations of simple tension and compression processes of a single crystal illustrates the development of the crystallographic axes according to the proposed evolution equation. In order to simulate the behaviour of a polycrystal the initial orientations of the anisotropy axes are assumed to be space-dependent but piecewise constant, where each region of a constant orientation corresponds to a grain. The results of the calculation show that the initially isotropic distribution of the orientation changes in a physically reasonable manner and that the intensity of this process-induced texture depends on the specific choice of the material constants.     

A constitutive model for orthotropic elasto-plasticity at large strains

Summary The paper presents a thermodynamically consistent constitutive model for elasto-plastic analysis of orthotropic materials at large strain. The elastic and plastic anisotropies are assumed to be persistent in the material but the anisotropy axes can undergo a rigid rotation due to large plastic deformations. The orthotropic yield function is formulated in terms of the generally nonsymmetric Mandel stress tensor such that its skew-symmetric part is additionally taken into account. Special attention is focused on the convexity of the yield surface resulting in the nine-dimensional stress space. Of particular interest are new convexity conditions which do not appear in the classical theory of anisotropic plasticity. They impose additional constraints on the material constants governing the plastic spin. The role of the plastic spin is further studied in simple shear accompanied by large elastic and large plastic deformations. If the plastic spin is neglected, the shear stress response is characterized by oscillations with an amplitude strictly dependent on the degree of the plastic anisotropy.accepted for publication 2 March 2004     

On Compressible Materials Capable of Sustaining Axisymmetric Shear Deformations. Part 4: Helical Shear of Anisotropic Hyperelastic Materials

Conditions on the form of the strain energy function in order that homogeneous, compressible and isotropic hyperelastic materials may sustain controllable static, axisymmetric anti-plane shear, azimuthal shear, and helical shear deformations of a hollow, circular cylinder have been explored in several recent papers. Here we study conditions on the strain energy function for homogeneous and compressible, anisotropic hyperelastic materials necessary and sufficient to sustain controllable, axisymmetric helical shear deformations of the tube. Similar results for separate axisymmetric anti-plane shear deformations and rotational shear deformations are then obtained from the principal theorem for helical shear deformations. The three theorems are illustrated for general compressible transversely isotropic materials for which the isotropy axis coincides with the cylinder axis. Previously known necessary and sufficient conditions on the strain energy for compressible and isotropic hyperelastic materials in order that the three classes of axisymmetric shear deformations may be possible follow by specialization of the anisotropic case. It is shown that the required monotonicity condition for the isotropic case is much simpler and less restrictive. Restrictions necessary and sufficient for anti-plane and rotational shear deformations to be possible in compressible hyperelastic materials having a helical axis of transverse isotropy that winds at a constant angle around the tube axis are derived. Results for the previous case and for a circular axis of transverse isotropy are included as degenerate helices. All of the conditions derived here have essentially algebraic structure and are easy to apply. The general rules are applied in several examples for specific strain energy functions of compressible and homogeneous transversely isotropic materials having straight, circular, and helical axes of material symmetry.     

Anisotropic Hyperelasticity Based Upon General Strain Measures

This paper presents stress-strain constitutive equations for anisotropic elastic materials. A special attention is given to the logarithmic strain. Assuming a constitutive equation for the specific internal energy the equation governing the Cauchy stress is derived. Mathematical relations presented take a relatively simple form and concern a very wide class of elastic materials. The dependence of third-order elastic constants on the choice of strain measure is shown. The third-order elastic constants measured experimentally in relation to the Green strain are recalculated here for the logarithmic strain. This revised version was published online in July 2006 with corrections to the Cover Date.     

Strain energy density bounds for linear anisotropic elastic materials

Upper and lower bounds are presented for the magnitude of the strain energy density in linear anisotropic elastic materials. One set of bounds is given in terms of the magnitude of the stress field, another in terms of the magnitude of the strain field. Explicit algebraic formulas are given for the bounds in the case of cubic, transversely isotropic, hexagonal and tetragonal symmetry. In the case of orthotropic symmetry the explicit bounds depend upon the solution of a cubic equation, and in the case of the monoclinic and triclinic symmetries, on the solution of sixth order equations.     

Physical invariants for nonlinear orthotropic solids

Seven invariants, with immediate physical interpretation, are proposed for the strain energy function of nonlinear orthotropic elastic solids. Three of the seven invariants are the principal stretch ratios and the other four are squares of the dot product between the two preferred directions and two principal directions of the right stretch tensor. A strain energy function, expressed in terms of these invariants, has a symmetrical property almost similar to that of an isotropic elastic solid written in terms of principal stretches. Ground state and stress–strain relations are given. Using principal axes techniques, the formulation is applied, with mathematical simplicity, to several types of deformations. In simple shear, a necessary and sufficient condition is given for Poynting relation and two novel deformation-dependent universal relations are formulated. Using series expansions and the symmetrical property, the proposed general strain energy function is refined to a particular general form. A type of strain energy function, where the ground state constants are written explicitly, is proposed. Some advantages of this type of function are indicated. An experimental advantage is demonstrated by showing a simple triaxial test can vary a single invariant while keeping the remaining invariants fixed.     

Acoustoelasticity in a weakly anisotropic monoclinic elastic material

The acoustoelasticity in a stressed monoclinic elastic material is analyzed theoretically. It is assumed that the material has weak anisotropy, such that the second-order elastic constants differ slightly from those of an isotropic material and the third-order elastic constants retain general monoclinic anisotropy. The propagation velocities, the polarization directions and the acoustoelastic effects for principal longitudinal and transverse waves are obtained and presented as functions of the elastic constants, principal stresses and directions of principal axes of stress. The coefficients appearing in the formulas are tabulated for Laue groups.     

The elliptic paraboloid failure surface for transversely isotropic materials off-axis loaded

The elliptic paraboloid failure surface has been well established as a potential criterion for yielding and failure of transversely isotropic materials, presenting also the strength differential effect [1]. This was done by extending well established criteria for isotropic materials presenting the strength differential effect (SDE), through an introduction process which maintained basic physical principles for the anisotropic materials. All previous literature concerned the special case where the principal axes of the external loading coincided with the principal strength axes of the material. In this paper the most general case where the two systems of frames are arbitrarily oriented relatively to each other is considered. In this situation the simplifications derived from the coincidence of the external principal stress and material principal strength axes are lost and the material should be considered as a general orthotropic one. The general properties for such types of loading of the transversely isotropic material are established by maintaining the general features of the failure locus invariant. Then, this study completes the investigation of yielding and failure mode of such materials considering the most general case of their loading.     

Crack trajectories influenced by mechanical and thermal disturbance in anisotropic material

This work is concerned with the application of the volume energy density criterion for predicting the crack trajectories as influenced by mechanical and thermal disturbance in an anisotropic material. Two-dimensional linear thermoelasticity is employed in conjunction with the well-known complex potentials such that a linear relationship is obtained for the boundary conditions across the crack or line of discontinuity. Boundary collocation is then used to determine the unknown coefficients from which the contours of the volume energy density in the cracked plate can be obtained. The crack path is assumed to coincide with the loci where dilatation would dominate. This corresponds to the locations of relative minimum energy density which can be found by visual inspection. An equal and opposite mechanical stress and thermal gradient are applied on the cracked plate. The former and latter enhance symmetric and asymmetric crack growth, respectively. They would complete depending on the magnitude of the mechanical and thermal load. Numerical results are presented for three (3) different cases of a plate whose principal axes of material symmetry are tilted to the crack plane. The influence of anisotropy on crack path is found to be secondary.     

Acoustic axes in elasticity

New results are presented for the degeneracy condition of elastic waves in anisotropic materials. The condition for the existence of acoustic axes involves a traceless symmetric third order tensor that must vanish identically. It is shown that all previous representations of the degeneracy condition follow from this acoustic axis tensor. The conditions for existence of acoustic axes in elastic crystals of orthorhombic, tetragonal, hexagonal and cubic (RTHC) symmetry are reinterpreted using the geometrical methods developed here. Application to weakly anisotropic solids is discussed, and it is shown that the satisfaction of the acoustic axes conditions to first order in anisotropy does not in general coincide with true acoustic axes.     

Dynamic damage constitutive relation of mesoscopic heterogenous brittle rock under rotation of principal stress axes

A micro-mechanics-based model is developed to investigate microcrack damage mechanism of four stages of brittle rock under rotation of the principal stress axes. They consist of linear elastic, non-linear hardening, rapid stress drop and strain softening. The frictional sliding crack model is applied to analyze microcracks nucleation, propagation and coalescence. The strain energy density factor approach is applied to determine the critical condition of microcrack nucleation, propagation and coalescence. The inelastic strain increments are formulated within the framework of thermodynamics with internal variables. Rotation of principal stress axes affect the dynamic damage constitutive relationship and the failure strength of brittle rock.     

Effect of texture and microstructure on strain hardening anisotropy for aluminum deformed in uniaxial tension and simple shear

Uniaxial and simple shear stress–strain curves were obtained for a 1050-O aluminum alloy sheet sample in different specimen orientations with respect to the material symmetry axes. For uniaxial tension, a strong anisotropy of strain hardening was observed leading to about 30% difference in uniform tensile elongation between the extreme conditions. For simple shear, the hardening was also significantly different. These results were rationalized with an analysis that accounts for dislocation substructure observations, crystallographic texture measurements and polycrystal modeling of texture-induced strength evolution.     

Incorporating the effects of fabric in the dilatant double shearing model for planar deformation of granular materials

The objective of this paper is to incorporate the effects of fabric and its evolution into the Dilatant Double Shearing Model [Mehrabadi, M.M., Cowin, S.C., 1978. Initial planar deformation of dilatant granular materials. J. Mech. Phys. Solids 26, 269–284] for granular materials in order to capture the anisotropic behavior and the complex response of granular materials in cyclic shear loading. An important consequence of considering the fabric is that one can have unequal shearing rates along the two slip directions. This property leads to the non-coaxiality of the principal axes of stress and strain rate, which is more appropriate for a material that exhibits initial and induced anisotropy. In addition, we employ a fabric-dependent elasticity tensor with orthotropic symmetry. The model developed in this paper also predicts one of the experimentally observed characteristics of granular materials: the gradual concentration of the contact normals towards the maximum principal stress direction.We implement the constitutive equations into ABAQUS/Explicit by writing a user material subroutine in order to predict the strength anisotropy of granular materials in a plane strain biaxial compression test and investigate the mechanical behavior of granular materials under the cyclic shear loading conditions. The predictions from this model show good quantitative agreement with the experiments of [Park, C.S., 1990. Anisotropy in deformation and strength properties of sands in plane strain compression, Masters Thesis, University of Tokyo; Park, C.S., Tatsuoka, F., 1994. Anisotropic strength and deformation of sands in plane strain compression. In: XIII ICSMFE, New Delhi, India; Okada, N., 1992. Energy dissipation in inelastic flow of cohesionless granular media. Ph.D. Thesis, University of California, San Diego].