Finite strain inelasticity for isotropy,a simple and efficient finite element formulation

We consider a formulation of associative isotropic J₂‐elastoplasticity at finite inelastic strains and aspects of its numerical implementation. The essential ingredients include the multiplicative decomposition of the deformation gradient in elastic and inelastic parts, the definition of a convex elastic domain in stress space and a material representation of the constitutive equations for general non‐Cartesian coordinate charts. On the numerical side we propose a stress update algorithm for elasto‐plastic response, including isotropic hardening. The finite element formulation is based on assumed strain and enhanced strain variational principles, for a complete outline see [3]. Remarkably the formulation is very similar to the case of infinitesimal plasticity: (i) The scheme of linear return mapping algorithm takes the form of standard return mapping of the infinitesimal theory for the case of isotropic elastic response. (ii) The algorithmic elastoplastic moduli have a similar structure as in the linear case. Together with an exact fulfillment of plastic incompressibility by means of a simple correction one achieves an advantageously efficient finite element formulation. Its performance is documented by a numerical example.     

Simulation of micro-cutting considering finite deformation crystal plasticity

Micro-machining processes on metalic microstructures are influenced by the crystal structure, i. e. the grain orientation. Furthermore, the chip formation underlies large deformations. To perform finite element simulations of micro-cutting processes, a large deformation material model is necessary in order to model the hyperelastic and finite plastic material behaviour. In the case of cp-titanium material with hcp-crystal structure the anisotropic behaviour must be considered by an appropriate set of slip planes and slip directions. In the present work the impact of the grain orientation on the plastic deformation is demonstrated by means of finite element simulations of a finite deformation single slip crystal plasticity model. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Self-consistent finite element method: A new method of predicting effective properties of inclusion media

The self-consistent method is a microchemical model for predicting the effective elastic properties of an inclusion medium. A numerical method based on self-consistent theory, namely the self-consistent finite element method, is developed. This new method can be applied to finding the determination of the effective properties of multiphase media with arbitrarily shaped and anisotropic inclusions. Applications to fibre composites demonstrate the implementation and accuracy of the method. This method can be extended to the elastoplastic and finite deformation case.     

Application of Material Forces to Hyperelastodynamics

In the framework of numerical simulation, the analysis of fracture and damage mechanisms play an increasing role. In this context, the configurational force method is an useful tool for investigating structures with the finite element method [1] to predict the material behaviour more precisely. In this paper, we discuss the configurational forces in the context of high strain rate loaded hyperelastic structures [6]. Balance of momentum in the material space and the weak formulation are shown with focus on the explicite finite element method. Finally, an application will be presented based on a detailed treatment of structures subjected to short time loading. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Finite element modelling of anisotropic elastic–viscoplastic behaviour of metals

An implementation of the unified theory of visco-plasticity of Bodner in a three-dimensional finite element program for the analysis of anisotropic inelastic behaviour of selected metals is presented in this paper. A derivation of an effective hardening parameter for the anisotropic (directional) deformation state is also given in this paper using some basic assumptions introduced by Bodner. The effect of the imposed strain rate on the level of the stress–strain curve is also investigated. A comparison of the results of the present finite element model with some published theoretical and experimental results for pure titanium and 2024-T4 aluminium alloy is also made.     

Finite Element Models of Hyperelastic Materials Based on a New Strain Measure

To construct constitutive equations for hyperelastic materials, one increasingly often proposes new strain measures, which result in significant simplifications and error reduction in experimental data processing. One such strain measure is based on the upper triangular (QR) decomposition of the deformation gradient. We describe a finite element method for solving nonlinear elasticity problems in the framework of finite strains for the case in which the constitutive equations are written with the use of the QR-decomposition of the deformation gradient. The method permits developing an efficient, easy-to-implement tool for modeling the stress–strain state of any hyperelastic material.     

A generalized polyconvex hyperelastic model for anisotropic solids

A generalized polyconvex hyperelastic model for anisotropic solids is presented. The strain energy function is formulated in terms of convex functions of generalized invariants and is given by a series with an arbitrary number of terms. The model addresses solids with orthotropic or transversely isotropic material symmetry as well as fiber-reinforced materials. Special cases of the strain energy function suitable for anisotropic elastomers and soft biological tissues are proposed. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Extension of the Arruda-Boyce model to the modelling of the curing process of polymers

A phenomenologically motivated finite strain general framework to simulate the curing of polymer have been developed and discussed in our recently published papers [2,4]. The Arruda-Boyce model is a classical hyperelastic model for polymeric materials. This contribution presents an extension of the Arruda-Boyce model towards modelling the curing process of polymers following our previous framework. In this paper, we will show how to model the elastic behaviour and shrinkage effects of the polymer curing process in the isothermal case using the Arruda-Boyce model. Several numerical examples have been demonstrated to verify our newly proposed modified approach in case of curing process. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Automatic energy–momentum conserving time integrators for hyperelastic waves

An energy–momentum conserving time integrator coupled with an automatic finite element algorithm is developed to study longitudinal wave propagation in hyperelastic layers. The Murnaghan strain energy function is used to model material nonlinearity and full geometric nonlinearity is considered. An automatic assembly algorithm using algorithmic differentiation is developed within a discrete Hamiltonian framework to directly formulate the finite element matrices without recourse to an explicit derivation of their algebraic form or the governing equations. The algorithm is illustrated with applications to longitudinal wave propagation in a thin hyperelastic layer modeled with a two-mode kinematic model. Solution obtained using a standard nonlinear finite element model with Newmark time stepping is provided for comparison.     

A comparison between inverse form finding and shape optimization methods for anisotropic hyperelasticity in logarithmic strain space

In this paper an inverse mechanical formulation and a Limited-Broyden-Flechter-Goldfarb-Shanno method for shape optimization are compared. Both methods deal with the determination of the undeformed shape of an hyperelastic part knowing its deformed configuration and the applied loads. We consider anisotropic hyperelastic materials that are formulated in the logarithmic strain space. Beside the theoretical aspects, we present a numerical example. We established that no difference could be found between the node coordinates on the undeformed sheets computed with both methods. However the convergence to the solution is faster for the inverse mechanical formulation compared to the L-BFGS algorithm. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Anisotropic mesh generation methods based on ACVT and natural metric for anisotropic elliptic equation

Anisotropic meshes are known to be well-suited for problems which exhibit anisotropic solution features.Defning an appropriate metric tensor and designing an efcient algorithm for anisotropic mesh generation are two important aspects of the anisotropic mesh methodology.In this paper,we are concerned with the natural metric tensor for use in anisotropic mesh generation for anisotropic elliptic problems.We provide an algorithm to generate anisotropic meshes under the given metric tensor.We show that the inverse of the anisotropic difusion matrix of the anisotropic elliptic problem is a natural metric tensor for the anisotropic mesh generation in three aspects:better discrete algebraic systems,more accurate fnite element solution and superconvergence on the mesh nodes.Various numerical examples demonstrating the efectiveness are presented.     

Calculation of Stresses and Consistent Tangent Moduli from Automatic Differentiation of Hyerelastic Strain Energy Functions Through the Use of Hyper Dual Numbers

Many materials as e.g. engineering rubbers, polymers and soft biological tissues are often described by hyperelastic strain energy functions. For their finite element implementation the stresses and consistent tangent moduli are required and obtained mainly in terms of the first and second derivative of the strain energy function. Depending on its mathematical complexity in particular for anisotropic media the analytic derivatives may be expensive to be calculated or implemented. Then numerical approaches may be a useful alternative reducing the development time. Often-used classical finite difference schemes are however quite sensitive with respect to perturbation values and they result in a poor accuracy. The complex-step derivative approximation does never suffer from round-off errors, cf. [1], [2], but it can only provide first derivatives. A method which also provides higher order derivatives is based on hyper dual numbers [3]. This method is independent on the choice of perturbation values and does thus neither suffer from round-off errors nor from approximation errors. Therefore, here we make use of hyper dual numbers and propose a numerical scheme for the calculation of stresses and tangent moduli which are almost identical to the analytic ones. Its uncomplicated implementation and accuracy is illustrated by some representative numerical examples. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The integration of analysis and testing for the simulation of the response of hyperelastic materials

A family of hyperelastic finite elements capable of modeling arbitrarily large strains for axisymmetric and plane strain analyses has been developed. Constitutive behavior is determined by the selection of a strain energy density function for which user-supplied coefficients are required. Selective reduced integration for the volumetric strain energy terms allows for successful modeling of nearly incompressible materials. Available strain energy density functions are as follows: Mooney-Rivlin, Blatz-Ko, power law, and a nine-term Mooney expansion. The Ogden Strain Energy (OSE) law has also been implemented. The OSE law defines the strain energy relationship entirely in terms of the three principal components of stretch. This differs from the approach of other strain energy formulations, such as the Mooney law in which the strain energy is written as a function of strain invariants. The OSE law as implemented in this formulation is designed to facilitate the user's task of converting physical test data to the numerical (algebraic) form required for input. The family of hyperelastic finite elements has been integrated into ANSYS Revision 4.2 via the user element interface. Numerous verification solutions have been performed. As a representative example, a comparison with a closed-form solution for a Mooney-Rivlin type material is presented. Finally, the difficulties of obtaining test data in the form of user-supplied constants is discussed in the context of the comparison of experimental measurements and analytical simulation of an elastomeric test specimen.     

Anisotropic mesh refinement for finite element methods based on error reduction

In this paper, we propose an anisotropic adaptive refinement algorithm based on the finite element methods for the numerical solution of partial differential equations. In 2-D, for a given triangular grid and finite element approximating space V, we obtain information on location and direction of refinement by estimating the reduction of the error if a single degree of freedom is added to V. For our model problem the algorithm fits highly stretched triangles along an interior layer, reducing the number of degrees of freedom that a standard h-type isotropic refinement algorithm would use.     

On two different inverse form finding methods for hyperelastic and elastoplastic materials

This contribution focuses on two different developments of mechanical-computational methods for the optimal determination of the initial shape of formed functional components knowing the deformed configuration, the applied loads and the boundary conditions. The first method uses an inverse mechanical formulation and can be applied to materials with hyperelastic behaviors. For materials with elastoplastic properties this method is not advocated, without knowing the final plastic strains, due to the non uniqueness of the solution. The second method uses a shape optimization formulation in the sense of an inverse problem via successive iterations of the direct problem. For hyperelastic materials the inverse mechanical formulation is preferred for its velocity and the non exhibition of possible mesh distortions. In the shape optimization formulation mesh distortions can be avoided by an update of the reference configuration of the functional part. Both methods are using a formulation in the logarithmic strain space. A numerical example for materials with isotropic elastoplastic behaviors illustrates the shape optimization formulation. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Formulation of hybrid Trefftz finite element method for elastoplasticity

The present investigation provides a hybrid Trefftz finite element approach for analysing elastoplastic problems. A dual variational functional is constructed and used to derive hybrid Trefftz finite element formulation for elastoplasticity of bulky solids. The formulation is applicable to either strain hardening or elastic-perfectly plastic materials. A solution algorithm based on initial stress formulation is introduced into the new element model. The performance of the proposed element model is assessed by three examples and comparison is made with results obtained by other approaches. The hybrid Trefftz finite element approach is demonstrated to be particularly suited for nonlinear analysis of two-dimensional elastoplastic problems.     

An Energy-Stable Parametric Finite Element Method for Simulating Solid-State Dewetting Problems in Three Dimensions

We propose an accurate and energy-stable parametric finite element method for solving the sharp-interface continuum model of solid-state dewetting in three-dimensional space. The model describes the motion of the film\slash vapor interface with contact line migration and is governed by the surface diffusion equation with proper boundary conditions at the contact line. We present a weak formulation for the problem, in which the contact angle condition is weakly enforced. By using piecewise linear elements in space and backward Euler method in time, we then discretize the formulation to obtain a parametric finite element approximation, where the interface and its contact line are evolved simultaneously. The resulting numerical method is shown to be well-posed and unconditionally energy-stable. Furthermore, the numerical method is generalized to the case of anisotropic surface energies in the Riemannian metric form. Numerical results are reported to show the convergence and efficiency of the proposed numerical method as well as the anisotropic effects on the morphological evolution of thin films in solid-state dewetting.     

非协调有限元的全塑性分析

本文采用满足相容条件的非协调有限元模型以解决全塑性分析中有限元解的数值精度问题.文中讨论了该模型适用于全塑性分析的机理和判据,还设计了一个确定塑性极限载荷的算法.     

A positive cell vertex Godunov scheme for a Beeler–Reuter based model of cardiac electrical activity

The monodomain model is a widely used model in electrocardiology to simulate the propagation of electrical potential in the myocardium. In this paper, we investigate a positive nonlinear control volume finite element scheme, based on Godunov's flux approximation of the diffusion term, for the monodomain model coupled to a physiological ionic model (the Beeler–Reuter model) and using an anisotropic diffusion tensor. In this scheme, degrees of freedom are assigned to vertices of a primal triangular mesh, as in conforming finite element methods. The diffusion term which involves an anisotropic tensor is discretized on a dual mesh using the diffusion fluxes provided by the conforming finite element reconstruction on the primal mesh and the other terms are discretized by means of an upwind finite volume method on the dual mesh. The scheme ensures the validity of the discrete maximum principle without any restriction on the transmissibility coefficients. By using a compactness argument, we obtain the convergence of the discrete solution and as a consequence, we get the existence of a weak solution of the original model. Finally, we illustrate by numerical simulations that the proposed scheme successfully removes nonphysical oscillations in the propagation of the wavefront and maintains conduction velocity close to physiological values.     

The mathematical modeling of the electric field in the media with anisotropic objects

We present a numerical scheme for modeling the electric field in the media with tensor conductivity. This scheme is based on vector finite element method in frequency domain. The numerical computations of the electric field in the anisotropic medium are done. The conductivity of the anisotropic medium is positive defined dense tensor in general case. We consider the electric field from anisotropic layer, inclined anisotropic layer and some anisotropic objects in isotropic half-space.