Computational structural and material stability analysis in finite electro‐elasto‐statics of electro‐active materials

Dielectric materials like electro‐active polymers (EAPs) exhibit coupled electro‐mechanical behavior at large strains. They respond by a deformation to an applied electrical field and are used in advanced industrial environments as sensors and actuators, for example, in robotics, biomimetics and smart structures. In field‐activated or electronic EAPs, the electric activation is driven by Coulomb‐type electrostatic forces, resulting in Maxwell stresses. These materials are able to provide finite actuation strains, which can even be improved by optimizing their composite microstructure. However, EAPs suffer from different types of instabilities. This concerns global structural instabilities, such as buckling and wrinkling of EAP devices, as well as local material instabilities, such as limit‐points and bifurcation‐points in the constitutive response, which induce snap‐through and fine scale localization of local states. In this work, we outline variational‐based definitions for structural and material stability, and design algorithms for accompanying stability checks in typical finite element computations. The formulation starts from stability criteria for a canonical energy minimization principle of electro‐elasto‐statics, and then shifts them over to representations related to an enthalpy‐based saddle point principle that is considered as the most convenient setting for numerical implementation. Here, global structural stability is analyzed based on a perturbation of the total electro‐mechanical energy, and related to statements of positive definiteness of incremental finite element tangent arrays. We base the local material stability on an incremental quasi‐convexity condition of the electro‐mechanical energy, inducing the positive definiteness of both the incremental electro‐mechanical moduli as well as a generalized acoustic tensor. It is shown that the incremental arrays to be analyzed in the stability criteria appear within the enthalpy‐based setting in a distinct diagonal form, with pure mechanical and pure electrical partitions. Applications of accompanying stability analyses in finite element computations are demonstrated by means of representative model problems. Copyright © 2015 John Wiley & Sons, Ltd.     

Variational‐based modeling of micro‐electro‐elasticity with electric field‐driven and stress‐driven domain evolutions

Recently, increasing interest in so‐called functional or smart materials with electromechanical coupling has been shown such as ferroelectric piezoceramics. These materials are characterized by microstructural properties, which can be changed by external stress and electric field stimuli, and hence find use as the active components in sensors and actuators. The electromechanical coupling effects result from the existence and rearrangement of microstructural domains with uniformly oriented electric polarization. The understanding and efficient simulation of these highly nonlinear and dissipative mechanisms, which occur on the microscale of ferroelectric piezoceramics, are a key challenge of the current research. This paper does not offer a substantially new physical model of these phenomena but a new mathematical modeling approach based on a rigorous exploitation of rate‐type variational principles. This provides a new insight in the structure of the coupled problem, where the governing field equations appear as the Euler equations of a variational statement. We outline a variational‐based micro‐electro‐elastic model for the microstructural evolution of both electrically and mechanically driven electric domains in ferroelectric ceramics, which also incorporates the surrounding free space. To this end, we extend recently developed multifield incremental variational principles of electromechanics from local to gradient‐extended dissipative response and specialize it by a Ginzburg–Landau‐type phase field model, where the thickness of the domain walls enters the formulation as a length scale. This serves as a natural starting point for a canonical compact, symmetric finite element implementation, considering the mechanical displacement, the microscopic polarization, and the electric potential induced by the polarization as the primary fields. The latter is defined on both the solid domain and a surrounding free space. Numerical simulations treat domain wall motions for electric field‐driven and stress‐driven loading processes, including the expansion of the electric potential into the free space. Copyright © 2012 John Wiley & Sons, Ltd.     

Computational electro‐elasticity and magneto‐elasticity for quasi‐incompressible media immersed in free space

In this work, a mixed variational formulation to simulate quasi‐incompressible electro‐active or magneto‐active polymers immersed in the surrounding free space is presented. A novel domain decomposition is used to disconnect the primary coupled problem and the arbitrary free‐space mesh update problem. Exploiting this decomposition, we describe a block‐iterative approach to solving the linearised multiphysics problem, and a physically and geometrically based, three‐parameter method to update the free space mesh. Several application‐driven example problems are implemented to demonstrate the robustness of the mixed formulation for both electro‐elastic and magneto‐elastic problems involving both finite deformations and quasi‐incompressible media. Copyright © 2016 John Wiley & Sons, Ltd.     

Phase‐field modeling of ductile fracture at finite strains: A robust variational‐based numerical implementation of a gradient‐extended theory by micromorphic regularization

This work provides a robust variational‐based numerical implementation of a phase field model of ductile fracture in elastic–plastic solids undergoing large strains. This covers a computationally efficient micromorphic regularization of the coupled gradient plasticity‐damage formulation. The phase field approach regularizes sharp crack surfaces within a pure continuum setting by a specific gradient damage modeling with geometric features rooted in fracture mechanics. It has proven immensely successful with regard to the analysis of complex crack topologies without the need for fracture‐specific computational structures such as finite element design of crack discontinuities or intricate crack‐tracking algorithms. The proposed gradient‐extended plasticity‐damage formulation includes two independent length scales that regularize both the plastic response as well as the crack discontinuities. This ensures that the damage zones of ductile fracture are inside of plastic zones or vice versa and guarantees on the computational side a mesh objectivity in post‐critical ranges. The proposed setting is rooted in a canonical variational principle. The coupling of gradient plasticity to gradient damage is realized by a constitutive work density function that includes the stored elastic energy and the dissipated work due to plasticity and fracture. The latter represents a coupled resistance to plasticity and damage, depending on the gradient‐extended internal variables that enter plastic yield functions and fracture threshold functions. With this viewpoint on the generalized internal variables at hand, the thermodynamic formulation is outlined for gradient‐extended dissipative solids with generalized internal variables that are passive in nature. It is specified for a conceptual model of von Mises‐type elasto‐plasticity at finite strains coupled with fracture. The canonical theory proposed is shown to be governed by a rate‐type minimization principle, which fully determines the coupled multi‐field evolution problem. This is exploited on the numerical side by a fully symmetric monolithic finite element implementation. An important aspect of this work is the regularization towards a micromorphic gradient plasticity‐damage setting by taking into account additional internal variable fields linked to the original ones by penalty terms. This enhances the robustness of the finite element implementation, in particular, on the side of gradient plasticity. The performance of the formulation is demonstrated by means of some representative examples. Copyright © 2016 John Wiley & Sons, Ltd.     

A linearised hp–finite element framework for acousto‐magneto‐mechanical coupling in axisymmetric MRI scanners

We propose a new computational framework for the treatment of acousto‐magneto‐mechanical coupling that arises in low‐frequency electro‐magneto‐mechanical systems such as magnetic resonance imaging scanners. Our transient Newton–Raphson strategy involves the solution of a monolithic system obtained from the linearisation of the coupled system of equations. Moreover, this framework, in the case of excitation from static and harmonic current sources, allows us to propose a simple linearised system and rigorously motivate a single‐step strategy for understanding the response of systems under different frequencies of excitation. Motivated by the need to solve industrial problems rapidly, we restrict ourselves to solving problems consisting of axisymmetric geometries and current sources. Our treatment also discusses in detail the computational requirements for the solution of these coupled problems on unbounded domains and the accurate discretisation of the fields using hp–finite elements. We include a set of academic and industrially relevant examples to benchmark and illustrate our approach. Copyright © 2017 The Authors. International Journal for Numerical Methods in Engineering Published by John Wiley & Sons, Ltd.     

A novel contact interaction formulation for voxel‐based micro‐finite‐element models of bone

Voxel‐based micro‐finite‐element (μFE) models are used extensively in bone mechanics research. A major disadvantage of voxel‐based μFE models is that voxel surface jaggedness causes distortion of contact‐induced stresses. Past efforts in resolving this problem have only been partially successful, ie, mesh smoothing failed to preserve uniformity of the stiffness matrix, resulting in (excessively) larger solution times, whereas reducing contact to a bonded interface introduced spurious tensile stresses at the contact surface. This paper introduces a novel “smooth” contact formulation that defines gap distances based on an artificial smooth surface representation while using the conventional penalty contact framework. Detailed analyses of a sphere under compression demonstrated that the smooth formulation predicts contact‐induced stresses more accurately than the bonded contact formulation. When applied to a realistic bone contact problem, errors in the smooth contact result were under 2%, whereas errors in the bonded contact result were up to 42.2%. We conclude that the novel smooth contact formulation presents a memory‐efficient method for contact problems in voxel‐based μFE models. It presents the first method that allows modeling finite slip in large‐scale voxel meshes common to high‐resolution image‐based models of bone while keeping the benefits of a fast and efficient voxel‐based solution scheme.     

A coupled hp‐finite element scheme for the solution of two‐dimensional electrostrictive materials

As part of the ongoing research within the field of computational analysis for the coupled electro‐magneto‐mechanical response of smart materials, the problem of linearised electrostriction is revisited and analysed for the first time using the framework of hp‐finite elements. The governing equations modelling the physics of the dielectric are suitably modified by introducing a new total Cauchy stress tensor (A. Dorfmann and R.W. Ogden. Nonlinear electroelasticity. Acta Mechanica, 174:167–183, 2005), which includes the electrostrictive effect and a staggered partitioned scheme for the numerical solution of the coupling phenomena. With the purpose of benchmarking numerical results, the problem of an infinite electrostrictive plate with a circular/elliptical dielectric insert is revisited. The presented analytical solution is based on the theoretical framework for two‐dimensional electrostriction proposed by Knops (R.J. Knops. Two‐dimensional electrostriction. Quarterly Journal of Mechanics and Applied Mathematics, 16:377–388, 1963) and uses classical techniques of complex variable analysis. Our presentation, to the best of our knowledge, provides the first correct closed form expression for the solution to the infinite electrostrictive plate with a circular/elliptical dielectric insert, correcting the errors made in previous presentations of this problem. We use this analytical solution to assess the accuracy, efficiency and robustness of the hp‐formulation in the case of nearly incompressible electrostrictive materials. Copyright © 2012 John Wiley & Sons, Ltd.     

An error‐estimate‐free and remapping‐free variational mesh refinement and coarsening method for dissipative solids at finite strains

A variational h‐adaptive finite element formulation is proposed. The distinguishing feature of this method is that mesh refinement and coarsening are governed by the same minimization principle characterizing the underlying physical problem. Hence, no error estimates are invoked at any stage of the adaption procedure. As a consequence, linearity of the problem and a corresponding Hilbert‐space functional framework are not required and the proposed formulation can be applied to highly non‐linear phenomena. The basic strategy is to refine (respectively, unrefine) the spatial discretization locally if such refinement (respectively, unrefinement) results in a sufficiently large reduction (respectively, sufficiently small increase) in the energy. This strategy leads to an adaption algorithm having O(N) complexity. Local refinement is effected by edge‐bisection and local unrefinement by the deletion of terminal vertices. Dissipation is accounted for within a time‐discretized variational framework resulting in an incremental potential energy. In addition, the entire hierarchy of successive refinements is stored and the internal state of parent elements is updated so that no mesh‐transfer operator is required upon unrefinement. The versatility and robustness of the resulting variational adaptive finite element formulation is illustrated by means of selected numerical examples. Copyright © 2008 John Wiley & Sons, Ltd.     

Formulation and numerical exploitation of mixed variational principles for coupled problems of Cahn–Hilliard‐type and standard diffusion in elastic solids

This work develops variational principles for the coupled problem of standard and extended Cahn–Hilliard‐type species diffusion in solids undergoing finite elastic deformations. It shows that the coupled problem of diffusion in deforming solids, accounting for phenomena like swelling, diffusion‐induced stress generation and possible phase segregation caused by the diffusing species, is related to an intrinsic mixed variational principle. It determines the rates of deformation and concentration along with the chemical potential, where the latter plays the role of a mixed variable. The principle characterizes a canonically compact model structure, where the three governing equations involved, that is, the mechanical equilibrium condition, the mass balance for the species content and a microforce balance that determines the chemical potential, appear as the Euler equation of a variational statement. The existence of the variational principle underlines an inherent symmetry of the coupled deformation–diffusion problem. This can be exploited in the numerical implementation by the construction of time‐discrete and space‐discrete incremental potentials, which fully determine the update problems of typical time stepping procedures. The mixed variational principles provide the most fundamental approach to the monolithic finite element solution of the coupled deformation–diffusion problem based on low‐order basis functions. They induce in a natural format the choice of symmetric solvers for Newton‐type iterative updates, providing a speedup and reduction of data storage when compared with non‐symmetric implementations. This is a strong argument for the use of the developed variational principles in the computational design of deformation–diffusion problems. Copyright © 2014 John Wiley & Sons, Ltd.     

Mixed finite element method for coupled thermo‐hydro‐mechanical process in poro‐elasto‐plastic media at large strains

A mixed finite element for coupled thermo‐hydro‐mechanical (THM) analysis in unsaturated porous media is proposed. Displacements, strains, the net stresses for the solid phase; pressures, pressure gradients, Darcy velocities for pore water and pore air phases; temperature, temperature gradients, the total heat flux are interpolated as independent variables. The weak form of the governing equations of coupled THM problems in porous media within the element is given on the basis of the Hu–Washizu three‐filed variational principle. The proposed mixed finite element formulation is derived. The non‐linear version of the element formulation is further derived with particular consideration of the THM constitutive model for unsaturated porous media based on the CAP model. The return mapping algorithm for the integration of the rate constitutive equation, the consistent elasto‐plastic tangent modulus matrix and the element tangent stiffness matrix are developed. For geometrical non‐linearity, the co‐rotational formulation approach is utilized. Numerical results demonstrate the capability and the performance of the proposed element in modelling progressive failure characterized by strain localization and the softening behaviours caused by thermal and chemical effects. Copyright © 2005 John Wiley & Sons, Ltd.     

Mixed plate bending elements based on least‐squares formulation

A finite element formulation for the bending of thin and thick plates based on least‐squares variational principles is presented. Finite element models for both the classical plate theory and the first‐order shear deformation plate theory (also known as the Kirchhoff and Mindlin plate theories, respectively) are considered. High‐order nodal expansions are used to construct the discrete finite element model based on the least‐squares formulation. Exponentially fast decay of the least‐squares functional, which is constructed using the L₂ norms of the equations residuals, is verified for increasing order of the nodal expansions. Numerical examples for the bending of circular, rectangular and skew plates with various boundary conditions and plate thickness are presented to demonstrate the predictive capability and robustness of the new plate bending elements. Plate bending elements based on this formulation are shown to be insensitive to both shear‐locking and geometric distortions. Copyright © 2004 John Wiley & Sons, Ltd.     

Cyclic steady states of nonlinear electro‐mechanical devices excited at resonance

We present an efficient numerical method to solve for cyclic steady states of nonlinear electro‐mechanical devices excited at resonance. Many electro‐mechanical systems are designed to operate at resonance, where the ramp‐up simulation to steady state is computationally very expensive – especially when low damping is present. The proposed method relies on a Newton–Krylov shooting scheme for the direct calculation of the cyclic steady state, as opposed to a naïve transient time‐stepping from zero initial conditions. We use a recently developed high‐order Eulerian–Lagrangian finite element method in combination with an energy‐preserving dynamic contact algorithm in order to solve the coupled electro‐mechanical boundary value problem. The nonlinear coupled equations are evolved by means of an operator split of the mechanical and electrical problem with an explicit as well as implicit approach. The presented benchmark examples include the first three fundamental modes of a vibrating nanotube, as well as a micro‐electro‐mechanical disk resonator in dynamic steady contact. For the examples discussed, we observe power law computational speed‐ups of the form S  = 0.6·ξ  ^? 0.8, where ξ is the linear damping ratio of the corresponding resonance frequency. Copyright © 2016 John Wiley & Sons, Ltd.     

A two‐dimensional ordinary,state‐based peridynamic model for linearly elastic solids

Peridynamics is a non‐local mechanics theory that uses integral equations to include discontinuities directly in the constitutive equations. A three‐dimensional, state‐based peridynamics model has been developed previously for linearly elastic solids with a customizable Poisson's ratio. For plane stress and plane strain conditions, however, a two‐dimensional model is more efficient computationally. Here, such a two‐dimensional state‐based peridynamics model is presented. For verification, a 2D rectangular plate with a round hole in the middle is simulated under constant tensile stress. Dynamic relaxation and energy minimization methods are used to find the steady‐state solution. The model shows m‐convergence and δ‐convergence behaviors when m increases and δ decreases. Simulation results show a close quantitative matching of the displacement and stress obtained from the 2D peridynamics and a finite element model used for comparison. Copyright © 2014 John Wiley & Sons, Ltd.     

Monolithic modelling of electro‐mechanical coupling in micro‐structures

The purpose of the present work is to model and to simulate the coupling between the electric and mechanical fields. A new finite element approach is proposed to model strong electro‐mechanical coupling in micro‐structures with capacitive effect. The proposed approach is based on a monolithic formulation: the electric and the mechanical fields are solved simultaneously in the same formulation. This method provides a tangent stiffness matrix for the total coupled problem which allows to determine accurately the pull‐in voltage and the natural frequency of electro‐mechanical systems such as MEMs. To illustrate the methodology results are shown for the analysis of a micro‐bridge. Copyright © 2005 John Wiley & Sons, Ltd.     

A three‐scale compressible microsphere model for hyperelastic materials

This paper presents a seminumerical homogenization framework for porous hyperelastic materials that is open for any hyperelastic microresponse. The conventional analytical homogenization schemes do apply to a limited number of elementary hyperelastic constitutive models. Within this context, we propose a general numerical scheme based on the homogenization of a spherical cavity in an incompressible unit hyperelastic solid sphere, which is denoted as the mesoscopic representative volume element (mRVE). The approach is applicable to any hyperelastic micromechanical response. The deformation field in the sphere is approximated via nonaffine kinematics proposed by Hou and Abeyaratne (JMPS 40:571‐592,1992). Symmetric displacement boundary conditions driven by the principal stretches of the deformation gradient are applied on the outer boundary of the mRVE. The macroscopic quantities, eg, stress and moduli expressions, are obtained by analytically derived pointwise geometric transformations. The macroscopic expressions are then computed numerically through quadrature rules applied in the radial and surface directions of the sphere. A three‐scale compressible microsphere model is derived from the developed seminumerical homogenization framework where the micro‐meso transition is based on the nonaffine microsphere model at every point of the mRVE. The numerical scheme developed for the derivation of macroscopic homogenized stresses and moduli terms as well as the modeling capability of the three‐scale microsphere model is investigated through representative boundary value problems.     

Extension of the node‐to‐segment contact element for surface‐expansion‐dependent contact laws

A class of friction laws depending on the measure of contact surface expansion is defined in the paper within the continuum contact mechanics framework. The nominal and spatial forms of constitutive relations are discussed, including incremental penalty relations. Further, an extended node‐to‐segment element is derived which is capable of treating surface‐expansion‐dependent contact laws in a consistent way. The approach is suitable for any kind of node‐to‐segment contact elements. Finally, the computational efficiency of the extended element as well as other possible approaches are illustrated by numerical examples relevant to metal forming applications. Copyright © 2001 John Wiley & Sons, Ltd.     

Analysis of three‐dimensional crack initiation and propagation using the extended finite element method

An Erratum has been published for this article in International Journal for Numerical Methods in Engineering 2005, 63(8): 1228. We present a new formulation and a numerical procedure for the quasi‐static analysis of three‐dimensional crack propagation in brittle and quasi‐brittle solids. The extended finite element method (XFEM) is combined with linear tetrahedral elements. A viscosity‐regularized continuum damage constitutive model is used and coupled with the XFEM formulation resulting in a regularized ‘crack‐band’ version of XFEM. The evolving discontinuity surface is discretized through a C⁰ surface formed by the union of the triangles and quadrilaterals that separate each cracked element in two. The element's properties allow a closed form integration and a particularly efficient implementation allowing large‐scale 3D problems to be studied. Several examples of crack propagation are shown, illustrating the good results that can be achieved. Copyright © 2005 John Wiley & Sons, Ltd.