On non-associative anisotropic finite plasticity fully coupled with isotropic ductile damage for metal forming

In this work, non-associative finite strain anisotropic elastoplasticity fully coupled with ductile damage is considered using a thermodynamically consistent framework. First, the kinematics of large strain based on multiplicative decomposition of the total transformation gradient using the rotating frame formulation, is recalled and different objective derivatives defined. By using different anisotropic equivalent stresses (quadratic and non-quadratic) in yield function and in plastic potential, the evolution equations for all the dissipative phenomena are deduced from the generalized normality rule applied to the plastic potential while the consistency condition is still applied to the yield function. The effect of the objective derivatives and the equivalent stresses (quadratic or non-quadratic) on the plastic flow anisotropy and the hardening evolution with damage is considered. Numerical aspects mainly related to the time integration of the fully coupled constitutive equations are discussed. Applications are made to the AISI 304 sheet metal by considering different loading paths as tensile, shear, plane tensile and bulge tests. For each loading path the effect of the rotating frame, the equivalent stress (quadratic or non-quadratic) and the normality rule (with respect to yield function or to the plastic potential) are discussed on the light of some available experimental results.     

A formulation of anisotropic continuum elastoplasticity at finite strains. Part I: Modelling

A constitutive model for anisotropic elastoplasticity at finite strains is developed together with its numerical implementation. An anisotropic elastic constitutive law is described in an invariant setting by use of structural tensors and the elastic strain measure C_e. The elastic strain tensor as well as the structural tensors are assumed to be invariant in relation to superimposed rigid body rotations. An anisotropic Hill-type yield criterion, described by a non-symmetric Eshelby-like stress tensor and further structural tensors, is developed, where use is made of representation theorems for functions with non-symmetric arguments. The model also considers non-linear isotropic hardening. Explicit results for the specific case of orthotropic anisotropy are given. The associative flow rule is employed and the features of the inelastic flow rule are discussed in full. It is shown that the classical definition of the plastic material spin is meaningless in conjunction with the present formulation. Instead, the study motivates an alternative definition, which is based on the demand that such a quantity must be dissipation-free, as the plastic material spin is in the case of isotropy. Equivalent spatial formulations are presented too. The full numerical treatment is considered in Part II.     

A simple framework for full-network hyperelasticity and anisotropic damage

A formulation of a constitutive behaviour law is proposed for hyperelastic materials, such that damage induced anisotropy can be accounted for continuously. The full-network approach with directional damage is adopted as a starting point. The full-network law with elementary strain energy density based on the inverse Langevin is chosen as a reference law which is cast into the proposed framework. This continuum formalism is then rewritten using spherical harmonics to capture damage directionality. The proposed formalism allows for an efficient (and systematic) expansion of complex non-linear anisotropic constitutive laws. A low order truncated expression of the resulting behaviour is shown to reproduce accurately the stress-strain curves of the exact behaviour laws.     

A novel approach for anisotropic hardening modeling. Part II: Anisotropic hardening in proportional and non-proportional loadings,application to initially isotropic material

The modeling of anisotropic hardening, in particular for non-proportional loading paths, is a challenging task for advanced macroscopic models. The complex distortion of the yield locus is related to the activation and cross-hardening of different slip systems, depending on crystallographic orientations. These physical mechanisms can be taken into account in polycrystalline models but the computation times are enormous. The novel approach detailed in Part I (Rousselier et al., 2009) consists in: (i) drastically reducing the number of crystallographic orientations to save the computation cost, (ii) applying a parameter calibration procedure to obtain a good agreement with the experimental database. This methodology is first applied here to the anisotropic hardening in the proportional loadings of the strongly anisotropic aluminum alloy of Part I. Very good modeling is achieved with only eight crystallographic orientations. Different levels of additional hardening in biaxial proportional loading as compared to uniaxial loading can be modeled with the same polycrystalline model. For this, only the parameter calibration has to be performed with different databases. The same methodology has also been applied for the modeling of isotropic behavior. The best compromise between model accuracy and numerical cost is obtained with fourteen orientations. The deviations from isotropy are acceptable in all loading directions. Different levels of hardening in orthogonal loading: simple shear followed by simple tension, are achieved without any modification of the model equations. Only the parameter calibration has to be performed with different hardening levels in the database. FE calculations of a deep drawing test have been performed. The CPU time of the polycrystalline model is only five times larger than that with the simple von Mises model. The CPU time with texture evolution is further increased by a factor of two. The effects of texture evolution in rolling of the initially isotropic fcc material have been investigated. The resulting texture and hardening are qualitatively good.     

Finite element analysis of Volterra dislocations in anisotropic crystals: A thermal analogue

The present work gives a systematic and rigorous implementation of Volterra dislocations in ordinary two-dimensional finite elements using the thermal analogue and the integral representation of dislocations through the stresses. The full fields are given for edge dislocations in anisotropic crystals, and the Peach–Koehler forces are found for some important examples.     

Elements of a finite strain-gradient thermomechanical theory for material growth and remodeling

A second-gradient theory in finite strains is proposed to deal with the phenomena of material growth and remodeling, as happens in biomechanics, on account of mass transport and morphogenetic species. It involves first-order and second-order transplants (local structural rearrangements) and two material connections on the material manifold. It is shown that the evolution of these structural changes or “material inhomogeneities” is governed by Eshelby-like stress and hyperstress tensors. A thermodynamically admissible set of constitutive equations is proposed. The complexity due to the finite-strain gradient theory is a necessity in order to accommodate mass exchanges and diffusion of species.     

Plane strain consolidation of soil layer with anisotropic permeability

This paper presents an alternative analytical technique to study a plane strain consolidation of a poroelastic soil by taking into account the anisotropy of permeability. From the governing equations of a saturated poroelastic soil, the relationship of basic variables for a point of a soil layer is established between the ground surface （z=0） and the depth z in the Laplace-Fourier transform domain. Combined with the boundary conditions, an exact solution is derived for plane strain Biot＇s consolidation of a finite soil layer with anisotropic permeability in the transform domain. Numerical inversions of the Laplace transform and the Fourier transform are adopted to obtain the actual solution in the physical domain. Numerical results of plane strain Biot＇s consolidation for a single soil layer show that the anisotropic of permeability has a great influence on the consolidation behavior of the soils.     

Modeling of anisotropic wound healing

Biological soft tissues exhibit non-linear complex properties, the quantification of which presents a challenge. Nevertheless, these properties, such as skin anisotropy, highly influence different processes that occur in soft tissues, for instance wound healing, and thus its correct identification and quantification is crucial to understand them. Experimental and computational works are required in order to find the most precise model to replicate the tissues' properties. In this work, we present a wound healing model focused on the proliferative stage that includes angiogenesis and wound contraction in three dimensions and which relies on the accurate representation of the mechanical behavior of the skin. Thus, an anisotropic hyperelastic model has been considered to analyze the effect of collagen fibers on the healing evolution of an ellipsoidal wound. The implemented model accounts for the contribution of the ground matrix and two mechanically equivalent families of fibers. Simulation results show the evolution of the cellular and chemical species in the wound and the wound volume evolution. Moreover, the local strain directions depend on the relative wound orientation with respect to the fibers.     

The scattering of harmonic elastic anti-plane shear waves by two collinear cracks in anisotropic material plane by using the non-local theory

In this paper, the dynamic behavior of two collinear cracks in the anisotropic elasticity material plane subjected to the harmonic anti-plane shear waves is investigated by use of the nonlocal theory. To overcome the mathematical difficulties, a one-dimensional nonlocal kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress field near the crack tips. By use of the Fourier transform, the problem can be solved with the help of a pair of triple integral equations, in which the unknown variable is the displacement on the crack surfaces. To solve the triple integral equations, the displacement on the crack surfaces is expanded in a series of Jacobi polynomials. Unlike the classical elasticity solutions, it is found that no stress singularity is present near crack tips. The nonlocal elasticity solutions yield a finite hoop stress at the crack tips, thus allowing us to using the maximum stress as a fracture criterion. The magnitude of the finite stress field not only depends on the crack length but also on the frequency of the incident waves and the lattice parameter of the materials.     

Strain localization analysis using a large deformation anisotropic elastic–plastic model coupled with damage

Sheet metal forming processes generally involve large deformations together with complex loading sequences. In order to improve numerical simulation predictions of sheet part forming, physically-based constitutive models are often required. The main objective of this paper is to analyze the strain localization phenomenon during the plastic deformation of sheet metals in the context of such advanced constitutive models. Most often, an accurate prediction of localization requires damage to be considered in the finite element simulation. For this purpose, an advanced, anisotropic elastic–plastic model, formulated within the large strain framework and taking strain-path changes into account, has been coupled with an isotropic damage model. This coupling is carried out within the framework of continuum damage mechanics. In order to detect the strain localization during sheet metal forming, Rice’s localization criterion has been considered, thus predicting the limit strains at the occurrence of shear bands as well as their orientation. The coupled elastic–plastic-damage model has been implemented in Abaqus/implicit. The application of the model to the prediction of Forming Limit Diagrams (FLDs) provided results that are consistent with the literature and emphasized the impact of the hardening model on the strain-path dependency of the FLD. The fully three-dimensional formulation adopted in the numerical development allowed for some new results – e.g. the out-of-plane orientation of the normal to the localization band, as well as more realistic values for its in-plane orientation.     

On a finite-strain viscoplastic law coupled with anisotropic damage: theoretical formulations and numerical applications

Based on a dissipation inequality at finite strains and the effective stress concept, a Chaboche-type infinitesimal viscoplastic theory is extended to finite-strain cases coupled with anisotropic damage. The anisotropic damage is described by a rank-two symmetric tensor. The constitutive law is formulated in the corotational material coordinate system. Thus, the evolution equations of all internal variables can be expressed in terms of their material time derivatives. The numerical algorithm for implementing the material model in a finite element programme is also formulated, and several numerical examples are shown. Comparing the numerical simulations with experimental observations indicates that the present material model can describe well the primary, secondary and tertiary creep. It can also predict the anisotropic damage modes observed in experiments correctly.     

The extension of Beltrami's ellipsoid to anisotropic bodies

Summary The spectral decomposition of the compliance, stiffness, and failure tensors for transversely isotropic materials was studied and their characteristic values were calculated using the components of these fourth-rank tensors in a Cartesian frame defining the principal material directions. The spectrally decomposed compliance and stiffness or failure tensors for a transversely isotropic body (fiber-reinforced composite), and the eigenvalues derived from them define in a simple and efficient way the respective elastic eigenstates of the loading of the material. It has been shown that, for the general orthotropic or transversely isotropic body, these eigenstates consist of two double components, ₁ and ₂ which are shears (₂ being a simple shear and ₁, a superposition of simple and pure shears), and that they are associated with distortional components of energy. The remaining two eigenstates, with stress components ₃, and ₄, are the orthogonal supplements to the shear subspace of ₁ and ₂ and consist of an equilateral stress in the plane of isotropy, on which is superimposed a prescribed tension or compression along the symmetry axis of the material. The relationship between these superimposed loading modes is governed by another eigenquantity, the eigenangle .The spectral type of decomposition of the elastic stiffness or compliance tensors in elementary fourth-rank tensors thus serves as a means for the energy-orthogonal decomposition of the energy function. The advantage of this type of decomposition is that the elementary idempotent tensors to which the fourth-rank tensors are decomposed have the interesting property of defining energy-orthogonal stress states. That is, the stress-idempotent tensors are mutually orthogonal and at the same time collinear with their respective strain tensors, and therefore correspond to energy-orthogonal stress states, which are therefore independent of each other. Since the failure tensor is the limiting case for the respective _x, which are eigenstates of the compliance tensor S, this tensor also possesses the same remarkable property.An interesting geometric interpretation arises for the energy-orthogonal stress states if we consider the projections of _x in the principal₃D stress space. Then, the characteristic state ₂ vanishes, whereas stress states ₁, ₃ and ₄ are represented by three mutually orthogonal vectors, oriented as follows: The ₃ and ₄ lie on the principal diagonal plane (₃₁₂) with subtending angles equaling (–/2) and (-), respectively. On the positive principal ₃-axis, is the eigenangle of the orthotropic material, whereas the ₁-vector is normal to the (₃₁₂)-plane and lies on the deviatoric -plane. Vector ₂ is equal to zero.It was additionally conclusively proved that the four eigenvalues of the compliance, stiffness, and failure tensors for a transversely isotropic body, together with value of the eigenangle , constitute the five necessary and simplest parameters with which invariantly to describe either the elastic or the failure behavior of the body. The expressions for the _x-vector thus established represent an ellipsoid centered at the origin of the Cartesian frame, whose principal axes are the directions of the ₁-, ₃- and ₄-vectors. This ellipsoid is a generalization of the Beltrami ellipsoid for isotropic materials.Furthermore, in combination with extensive experimental evidence, this theory indicates that the eigenangle alone monoparametrically characterizes the degree of anisotropy for each transversely isotropic material. Thus, while the angle for isotropic materials is always equal to _i = 125.26° and constitutes a minimum, the angle || progressively increases within the interval 90–180° as the anisotropy of the material is increased. The anisotropy of the various materials, exemplified by their ratiosE _L/2G_L of the longitudinal elastic modulus to the double of the longitudinal shear modulus, increases rapidly tending asymptotically to very high values as the angle approaches its limits of 90 or 180°.     

Modelling the slight compressibility of anisotropic soft tissue

In order to avoid the numerical difficulties in locally enforcing the incompressibility constraint using the displacement formulation of the Finite Element Method, slight compressibility is typically assumed when simulating transversely isotropic, soft tissue. The current standard method of accounting for slight compressibility of hyperelastic soft tissue assumes an additive decomposition of the strain-energy function into a volumetric and a deviatoric part. This has been shown, however, to be inconsistent with the linear theory for anisotropic materials. It is further shown here that, under hydrostatic tension or compression, a transversely isotropic cube modelled using this additive split is simply deformed into another cube regardless of the size of the deformation, in contravention of the physics of the problem. A remedy for these defects is proposed here: the trace of the Cauchy stress is assumed linear in both volume change and fibre stretch. The general form of the strain-energy function consistent with this model is obtained and is shown to be a generalisation of the current standard model. A specific example is used to clearly demonstrate the differences in behaviour between the two models in hydrostatic tension and compression.     

A computational study of a model of single-crystal strain-gradient viscoplasticity with an interactive hardening relation

The behavior of a model of single-crystal strain-gradient viscoplasticity is investigated. The model is an extension of a rate-independent version, and includes a new hardening relation that has recently been proposed in the small-deformation context (Gurtin and Reddy, 2014), and which accounts for slip-system interactions due to self and latent hardening. Energetic and dissipative effects, each with its corresponding length scale, are included. Numerical results are presented for a single crystal with single and multiple slip systems, as well as an ensemble of grains. These results provide a clear illustration of the effects of accounting for slip-system interactions.     

Phase field modeling of fracture in anisotropic brittle solids

A phase field model of fracture that accounts for anisotropic material behavior at small and large deformations is outlined within this work. Most existing fracture phase field models assume crack evolution within isotropic solids, which is not a meaningful assumption for many natural as well as engineered materials that exhibit orientation-dependent behavior. The incorporation of anisotropy into fracture phase field models is for example necessary to properly describe the typical sawtooth crack patterns in strongly anisotropic materials. In the present contribution, anisotropy is incorporated in fracture phase field models in several ways: (i) Within a pure geometrical approach, the crack surface density function is adopted by a rigorous application of the theory of tensor invariants leading to the definition of structural tensors of second and fourth order. In this work we employ structural tensors to describe transverse isotropy, orthotropy and cubic anisotropy. Latter makes the incorporation of second gradients of the crack phase field necessary, which is treated within the finite element context by a nonconforming Morley triangle. Practically, such a geometric approach manifests itself in the definition of anisotropic effective fracture length scales. (ii) By use of structural tensors, energetic and stress-like failure criteria are modified to account for inherent anisotropies. These failure criteria influence the crack driving force, which enters the crack phase field evolution equation and allows to set up a modular structure. We demonstrate the performance of the proposed anisotropic fracture phase field model by means of representative numerical examples at small and large deformations.     

Nonlocal gradient-dependent modeling of plasticity with anisotropic hardening

This work addresses the formulation of the thermodynamics of nonlocal plasticity using the gradient theory. The formulation is based on the nonlocality energy residual introduced by Eringen and Edelen (1972). Gradients are introduced for those variables associated with isotropic and kinematic hardening. The formulation applies to small strain gradient plasticity and makes use of the evanescent memory model for kinematic hardening. This is accomplished using the kinematic flux evolution as developed by Zbib and Aifantis (1988). Therefore, the present theory is a four nonlocal parameter-based theory that accounts for the influence of large variations in the plastic strain, accumulated plastic strain, accumulated plastic strain gradients, and the micromechanical evolution of the kinematic flux. Using the principle of virtual power and the laws of thermodynamics, thermodynamically-consistent equations are derived for the nonlocal plasticity yield criterion and associated flow rule. The presence of higher-order gradients in the plastic strain is shown to enhance a corresponding history variable which arises from the accumulation of the plastic strain gradients. Furthermore, anisotropy is introduced by plastic strain gradients in the form of kinematic hardening. Plastic strain gradients can be attributed to the net Burgers vector, while gradients in the accumulation of plastic strain are responsible for the introduction of isotropic hardening. The equilibrium between internal Cauchy stress and the microstresses conjugate to the higher-order gradients frames the yield criterion, which is obtained from the principle of virtual power. Microscopic boundary conditions, associated with plastic flow, are introduced to supplement the macroscopic boundary conditions of classical plasticity. The nonlocal formulation developed here preserves the classical assumption of local plasticity, wherein plastic flow direction is governed by the deviatoric Cauchy stress. The theory is applied to the problems of thin films on both soft and hard substrates. Numerical solutions are presented for bi-axial tension and simple shear loading of thin films on substrates.