Finite element formulation for anisotropic coupled piezoelectro‐hygro‐thermo‐viscoelasto‐dynamic problems

A finite element algorithm has been developed for the efficient analysis of smart composite structures with piezoelectric polymer sensors or/and actuators based on piezoelectro‐hygro‐thermo‐viscoelasticity. Variational principles for anisotropic coupled piezoelectro‐hygro‐thermo‐viscoelasto‐dynamic problems have also been proposed in this study. As illustrative studies, dynamic responses in laminated composite beams and plates with PVDF sensors and actuators are obtained as functions of time using the present finite element procedures. The voltage feedback control scheme is utilized. The proposed numerical method can be used for analysing problems in the design of smart structures as well as smart sensors and actuators. Copyright © 1999 John Wiley & Sons, Ltd.     

An X‐FEM and level set computational approach for image‐based modelling: Application to homogenization

The advances in material characterization by means of imaging techniques require powerful computational methods for numerical analysis. The present contribution focuses on highlighting the advantages of coupling the extended finite elements method and the level sets method, applied to solve microstructures with complex geometries. The process of obtaining the level set data starting from a digital image of a material structure and its input into an extended finite element framework is presented. The coupled method is validated using reference examples and applied to obtain homogenized properties for heterogeneous structures. Although the computational applications presented here are mainly two‐dimensional, the method is equally applicable for three‐dimensional problems. Copyright © 2010 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 stabilized finite element formulation for monolithic thermo‐hydro‐mechanical simulations at finite strain

An adaptively stabilized monolithic finite element model is proposed to simulate the fully coupled thermo‐hydro‐mechanical behavior of porous media undergoing large deformation. We first formulate a finite‐deformation thermo‐hydro‐mechanics field theory for non‐isothermal porous media. Projection‐based stabilization procedure is derived to eliminate spurious pore pressure and temperature modes due to the lack of the two‐fold inf‐sup condition of the equal‐order finite element. To avoid volumetric locking due to the incompressibility of solid skeleton, we introduce a modified assumed deformation gradient in the formulation for non‐isothermal porous solids. Finally, numerical examples are given to demonstrate the versatility and efficiency of this thermo‐hydro‐mechanical model. Copyright © 2015 John Wiley & Sons, Ltd.     

An efficient finite element procedure for analyzing three‐phase porous media based on the relaxed Picard method

Effective simulation of the solid‐liquid‐gas coupling effect in unsaturated porous media is of great significance in many diverse areas. Because of the strongly nonlinear characteristics of the fully coupled formulations for the three‐phase porous media, an effective numerical solution scheme, such as the finite element method with an efficient iterative algorithm, has to be employed. In this paper, an efficient finite element procedure based on the adaptive relaxed Picard method is developed for analyzing the coupled solid‐liquid‐gas interactions in porous media. The coupled model and the finite element analysis procedure are implemented into a computer code PorousH2M, and the proposed procedure is validated through comparing the numerical simulations with the experimental benchmarks. It is shown that the adaptive relaxed Picard method has salient advantage over the traditional one with respect to both the efficiency and the robustness, especially for the case of relatively large time step sizes. Compared with the Newton‐Raphson scheme, the Picard method successfully avoids the unphysical ‘spurious unloading’ phenomenon under the plastic deformation condition, although the latter shows a better convergence rate. The proposed procedure provides an important reference for analyzing the fully coupled problems related to the multi‐phase, multi‐field coupling in porous media. Copyright © 2014 John Wiley & Sons, Ltd.     

Time‐step constraints in transient coupled finite element analysis

In transient finite element analysis, reducing the time‐step size improves the accuracy of the solution. However, a lower bound to the time‐step size exists, below which the solution may exhibit spatial oscillations at the initial stages of the analysis. This numerical ‘shock’ problem may lead to accumulated errors in coupled analyses. To satisfy the non‐oscillatory criterion, a novel analytical approach is presented in this paper to obtain the time‐step constraints using the θ‐method for the transient coupled analysis, including both heat conduction–convection and coupled consolidation analyses. The expressions of the minimum time‐step size for heat conduction–convection problems with both linear and quadratic elements reduce to those applicable to heat conduction problems if the effect of heat convection is not taken into account. For coupled consolidation analysis, time‐step constraints are obtained for three different types of elements, and the one for composite elements matches that in the literature. Finally, recommendations on how to handle the numerical ‘shock’ issues are suggested. Copyright © 2015 John Wiley & Sons, Ltd.     

Two‐scale diffusion–deformation coupling model for material deterioration involving micro‐crack propagation

In this paper, we introduce a two‐scale diffusion–deformation coupled model that represents the aging material deterioration of two‐phase materials involving micro‐crack propagations. The mathematical homogenization method is applied to relate the micro‐ and macroscopic field variables, and a weak coupling solution method is employed to solve the two‐way coupling phenomena between the diffusion of scalar fields and the deformation of quasi‐brittle solids. The macroscopic mechanical behavior represented by the derived two‐scale two‐way coupled model reveals material nonlinearity due to micro‐scale cracking induced by the scalar‐field‐induced deformation, which can be simulated by the finite cover method. After verifying the fundamental validity of the proposed model and the analysis method, we perform a simple numerical example to demonstrate their ability to predict aging material deterioration. Copyright © 2010 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.     

Fatigue life prediction for butt‐welded joints considering weld‐induced residual stresses and initial damage,relaxation of residual stress,and elasto‐plastic fatigue damage

Fatigue damage of butt‐welded joints is investigated by a damage mechanics method. First, the weld‐induced residual stresses are determined by using a sequentially coupled thermo‐mechanical finite element analysis. The plastic damage of material is then calculated with the use of Lemaitre's plastic damage model. Second, during the subsequent fatigue damage analysis, the residual stresses are superimposed on the fatigue loading, and the weld‐induced plastic damage is considered as the initial damage via an elasto‐plastic fatigue damage model. Finally, the fatigue damage evolution, the relaxation of residual stress, and the fatigue lives of the joints are evaluated using a numerical implementation. The predicted results agree well with the experimental data.     

Parallel computing of high‐speed compressible flows using a node‐based finite‐element method

An efficient parallel computing method for high‐speed compressible flows is presented. The numerical analysis of flows with shocks requires very fine computational grids and grid generation requires a great deal of time. In the proposed method, all computational procedures, from the mesh generation to the solution of a system of equations, can be performed seamlessly in parallel in terms of nodes. Local finite‐element mesh is generated robustly around each node, even for severe boundary shapes such as cracks. The algorithm and the data structure of finite‐element calculation are based on nodes, and parallel computing is realized by dividing a system of equations by the row of the global coefficient matrix. The inter‐processor communication is minimized by renumbering the nodal identification number using ParMETIS. The numerical scheme for high‐speed compressible flows is based on the two‐step Taylor–Galerkin method. The proposed method is implemented on distributed memory systems, such as an Alpha PC cluster, and a parallel supercomputer, Hitachi SR8000. The performance of the method is illustrated by the computation of supersonic flows over a forward facing step. The numerical examples show that crisp shocks are effectively computed on multiprocessors at high efficiency. Copyright © 2003 John Wiley & Sons, Ltd.     

Time‐domain hybrid formulations for wave simulations in three‐dimensional PML‐truncated heterogeneous media

We are concerned with the numerical simulation of wave motion in arbitrarily heterogeneous, elastic, perfectly‐matched‐layer‐(PML)‐truncated media. We extend in three dimensions a recently developed two‐dimensional formulation, by treating the PML via an unsplit‐field, but mixed‐field, displacement‐stress formulation, which is then coupled to a standard displacement‐only formulation for the interior domain, thus leading to a computationally cost‐efficient hybrid scheme. The hybrid treatment leads to, at most, third‐order in time semi‐discrete forms. The formulation is flexible enough to accommodate the standard PML, as well as the multi‐axial PML. We discuss several time‐marching schemes, which can be used à la carte, depending on the application: (a) an extended Newmark scheme for third‐order in time, either unsymmetric or fully symmetric semi‐discrete forms; (b) a standard implicit Newmark for the second‐order, unsymmetric semi‐discrete forms; and (c) an explicit Runge–Kutta scheme for a first‐order in time unsymmetric system. The latter is well‐suited for large‐scale problems on parallel architectures, while the second‐order treatment is particularly attractive for ready incorporation in existing codes written originally for finite domains. We compare the schemes and report numerical results demonstrating stability and efficacy of the proposed formulations. Copyright © 2014 John Wiley & Sons, Ltd.     

Two‐scale computational modelling of water flow in unsaturated soils containing irregular‐shaped inclusions

The focus of this paper is two‐dimensional computational modelling of water flow in unsaturated soils consisting of weakly conductive disconnected inclusions embedded in a highly conductive connected matrix. When the inclusions are small, a two‐scale Richards’ equation‐based model has been proposed in the literature taking the form of an equation with effective parameters governing the macroscopic flow coupled with a microscopic equation, defined at each point in the macroscopic domain, governing the flow in the inclusions. This paper is devoted to a number of advances in the numerical implementation of this model. Namely, by treating the micro‐scale as a two‐dimensional problem, our solution approach based on a control volume finite element method can be applied to irregular inclusion geometries, and, if necessary, modified to account for additional phenomena (e.g. imposing the macroscopic gradient on the micro‐scale via a linear approximation of the macroscopic variable along the microscopic boundary). This is achieved with the help of an exponential integrator for advancing the solution in time. This time integration method completely avoids generation of the Jacobian matrix of the system and hence eases the computation when solving the two‐scale model in a completely coupled manner. Numerical simulations are presented for a two‐dimensional infiltration problem. Copyright © 2014 John Wiley & Sons, Ltd.     

An optimal high‐order non‐reflecting finite element scheme for wave scattering problems

A new finite element scheme is proposed for the numerical solution of time‐harmonic wave scattering problems in unbounded domains. The infinite domain in truncated via an artificial boundary ?? which encloses a finite computational domain Ω. On ?? a local high‐order non‐reflecting boundary condition (NRBC) is applied which is constructed to be optimal in a certain sense. This NRBC is implemented in a special way, by using auxiliary variables along the boundary ??, so that it involves no high‐order derivatives regardless of its order. The order of the scheme is simply an input parameter, and it may be arbitrarily high. This leads to a symmetric finite element formulation where standard C⁰ finite elements are used in Ω. The performance of the method is demonstrated via numerical examples, and it is compared to other NRBC‐based schemes. The method is shown to be highly accurate and stable, and to lead to a well‐conditioned matrix problem. Copyright © 2002 John Wiley & Sons, Ltd.     

Variational principles in dissipative electro‐magneto‐mechanics: A framework for the macro‐modeling of functional materials

This paper presents a general framework for the macroscopic, continuum‐based formulation and numerical implementation of dissipative functional materials with electro‐magneto‐mechanical couplings based on incremental variational principles. We focus on quasi‐static problems, where mechanical inertia effects and time‐dependent electro‐magnetic couplings are a priori neglected and a time‐dependence enters the formulation only through a possible rate‐dependent dissipative material response. The underlying variational structure of non‐reversible coupled processes is related to a canonical constitutive modeling approach, often addressed to so‐called standard dissipative materials. It is shown to have enormous consequences with respect to all aspects of the continuum‐based modeling in macroscopic electro‐magneto‐mechanics. At first, the local constitutive modeling of the coupled dissipative response, i.e. stress, electric and magnetic fields versus strain, electric displacement and magnetic induction, is shown to be variational based, governed by incremental minimization and saddle‐point principles. Next, the implications on the formulation of boundary‐value problems are addressed, which appear in energy‐based formulations as minimization principles and in enthalpy‐based formulations in the form of saddle‐point principles. Furthermore, the material stability of dissipative electro‐magneto‐mechanics on the macroscopic level is defined based on the convexity/concavity of incremental potentials. We provide a comprehensive outline of alternative variational structures and discuss details of their computational implementation, such as formulation of constitutive update algorithms and finite element solvers. From the viewpoint of constitutive modeling, including the understanding of the stability in coupled electro‐magneto‐mechanics, an energy‐based formulation is shown to be the canonical setting. From the viewpoint of the computational convenience, an enthalpy‐based formulation is the most convenient setting. A numerical investigation of a multiferroic composite demonstrates perspectives of the proposed framework with regard to the future design of new functional materials. Copyright © 2011 John Wiley & Sons, Ltd.     

Local and non‐local thermomechanical modeling of elastic‐plastic materials undergoing large strains

This paper deals with the numerical analysis of instabilities for elastic‐plastic materials undergoing large deformations in non‐isothermal conditions. The considered isotropic model is fully thermomechanically coupled and includes temperature‐induced softening, which is another source of strain localization next to geometrical effects. Due to complexity of the model, a symbolic‐numerical tool Ace is used for the preparation of user‐supplied subroutines for the finite element method. The computational verification is performed using two benchmark tests: necking of circular bar in tension and shear banding of elongated rectangular plate in plain strain conditions. The attention is focused on mesh dependence of the numerical results and the regularizing effect of heat conduction. The research reveals that the conductivity influences the shear band width and ductility of the material response; however, for the adiabatic case, the results are discretization sensitive, and another regularization is needed. A new gradient‐enhanced thermomechanical model is developed that introduces an internal length parameter governing the size of the shear band caused by thermal softening. The numerical verification of the non‐local model is performed for the adiabatic case. Subsequently, the simultaneous application of the gradient enhancement and heat conduction in the model is analyzed, which reproduces an evolving shear band. Copyright © 2016 John Wiley & Sons, Ltd.     

A novel hybrid FS‐FEM/SEA for the analysis of vibro‐acoustic problems

Predicting the frequency response of a complex vibro‐acoustic system becomes extremely difficult in the mid‐frequency regime. In this work, a novel hybrid face‐based smoothed finite element method/statistical energy analysis (FS‐FEM/SEA) method is proposed, aiming to further improve the accuracy of ‘mid‐frequency’ predictions. According to this approach, the whole vibro‐acoustic system is divided into a combination of a plate subsystem with statistical behaviour and an acoustic cavity subsystem with deterministic behaviour. The plate subsystem is treated using the recently developed FS‐FEM, and the cavity subsystem is dealt with using the SEA. These two different types of subsystems can be coupled and interacted through the so‐called diffuse field reciprocity relation. The ensemble average response of the system is calculated, and the uncertainty is confined and treated in the SEA subsystems. The use of FS‐FEM ‘softens’ the well‐known ‘overly stiff’ behaviour in the standard FEM and reduces the inherent numerical dispersion error. The proposed FS‐FEM/SEA approach is verified and its features are examined by various numerical examples. Copyright © 2015 John Wiley & Sons, Ltd.     

An adaptive non‐conforming finite‐element method for Reissner–Mindlin plates

Adaptive algorithms are important tools for efficient finite‐element mesh design. In this paper, an error controlled adaptive mesh‐refining algorithm is proposed for a non‐conforming low‐order finite‐element method for the Reissner–Mindlin plate model. The algorithm is controlled by a reliable and efficient residual‐based a posteriori error estimate, which is robust with respect to the plate's thickness. Numerical evidence for this and the efficiency of the new algorithm is provided in the sense that non‐optimal convergence rates are optimally improved in our numerical experiments. Copyright © 2003 John Wiley & Sons, Ltd.     

Coupled meshfree and fractal finite element method for mixed mode two‐dimensional crack problems

This paper presents a coupling technique for integrating the element‐free Galerkin method (EFGM) with the fractal finite element method (FFEM) for analyzing homogeneous, isotropic, and two‐dimensional linear‐elastic cracked structures subjected to mixed‐mode (modes I and II) loading conditions. FFEM is adopted for discretization of the domain close to the crack tip and EFGM is adopted in the rest of the domain. In the transition region interface elements are employed. The shape functions within interface elements which comprise both the EFG and the finite element (FE) shape functions, satisfies the consistency condition thus ensuring convergence of the proposed coupled EFGM–FFEM. The proposed method combines the best features of EFGM and FFEM, in the sense that no special enriched basis functions or no structured mesh with special FEs are necessary and no post‐processing (employing any path independent integrals) is needed to determine fracture parameters, such as stress‐intensity factors (SIFs) and T‐stress. The numerical results show that SIFs and T‐stress obtained using the proposed method are in excellent agreement with the reference solutions for the structural and crack geometries considered in the present study. Also, a parametric study is carried out to examine the effects of the integration order, the similarity ratio, the number of transformation terms, and the crack length to width ratio on the quality of the numerical solutions. A numerical example on mixed‐mode condition is presented to simulate crack propagation. As in the proposed coupled EFGM–FFEM at each increment during the crack propagation, the FFEM mesh (around the crack tip) is shifted as it is to the new updated position of the crack tip (such that FFEM mesh center coincides with the crack tip) and few meshless nodes are sprinkled in the location where the FFEM mesh was lying previously, crack‐propagation analysis can be dramatically simplified. Copyright © 2010 John Wiley & Sons, Ltd.