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.     

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.     

Finite element solution of free‐surface ship‐wave problems

An unstructured finite element solver to evaluate the ship‐wave problem is presented. The scheme uses a non‐structured finite element algorithm for the Euler or Navier–Stokes flow as for the free‐surface boundary problem. The incompressible flow equations are solved via a fractional step method whereas the non‐linear free‐surface equation is solved via a reference surface which allows fixed and moving meshes. A new non‐structured stabilized approximation is used to eliminate spurious numerical oscillations of the free surface. Copyright © 1999 John Wiley & Sons, Ltd.     

Implementation of an exact finite reduction scheme for steady‐state reaction‐diffusion equations

In this paper we propose the numerical solution of a steady‐state reaction‐diffusion problem by means of application of a non‐local Lyapunov–Schmidt type reduction originally devised for field theory. A numerical algorithm is developed on the basis of the discretization of the differential operator by means of simple finite differences. The eigendecomposition of the resulting matrix is used to implement a discrete version of the reduction process. By the new algorithm the problem is decomposed into two coupled subproblems of different dimensions. A large subproblem is solved by means of a fixed point iteration completely controlled by the features of the original equation, and a second problem, with dimensions that can be made much smaller than the former, which inherits most of the non‐linear difficulties of the original system. The advantage of this approach is that sophisticated linearization strategies can be used to solve this small non‐linear system, at the expense of a partial eigendecomposition of the discretized linear differential operator. The proposed scheme is used for the solution of a simple non‐linear one‐dimensional problem. The applicability of the procedure is tested and experimental convergence estimates are consolidated. Numerical results are used to show the performance of the new algorithm. Copyright © 2006 John Wiley & Sons, Ltd.     

An implicit–explicit solution method for electro‐osmotic flow through three‐dimensional micro‐channels

This paper presents a comprehensive finite‐element modelling approach to electro‐osmotic flows on unstructured meshes. The non‐linear equation governing the electric potential is solved using an iterative algorithm. The employed algorithm is based on a preconditioned GMRES scheme. The linear Laplace equation governing the external electric potential is solved using a standard pre‐conditioned conjugate gradient solver. The coupled fluid dynamics equations are solved using a fractional step‐based, fully explicit, artificial compressibility scheme. This combination of an implicit approach to the electric potential equations and an explicit discretization to the Navier–Stokes equations is one of the best ways of solving the coupled equations in a memory‐efficient manner. The local time‐stepping approach used in the solution of the fluid flow equations accelerates the solution to a steady state faster than by using a global time‐stepping approach. The fully explicit form and the fractional stages of the fluid dynamics equations make the system memory efficient and free of pressure instability. In addition to these advantages, the proposed method is suitable for use on both structured and unstructured meshes with a highly non‐uniform distribution of element sizes. The accuracy of the proposed procedure is demonstrated by solving a basic micro‐channel flow problem and comparing the results against an analytical solution. The comparisons show excellent agreement between the numerical and analytical data. In addition to the benchmark solution, we have also presented results for flow through a fully three‐dimensional rectangular channel to further demonstrate the application of the presented method. Copyright © 2007 John Wiley & Sons, Ltd.     

Hybrid‐stress six‐node prismatic elements

This paper presents three novel hybrid‐stress six‐node prismatic elements. Starting from the element displacement interpolation, the equilibrating non‐constant stress modes for the first element are identified and orthogonalized with respect to the constant stress modes for higher computational efficiency. For the second element, the non‐constant stress modes are non‐equilibrating and chosen for the sake of stabilizing the reduced‐integrated element. The first two elements are intended for three‐dimensional continuum analysis with both passing the patch test for three‐dimensional continuum elements. The third element is primarily intended for plate/shell analysis. Shear locking is alleviated by a new assumed strain scheme which preserves the element accuracy with respect to the twisting load. Furthermore, the Poisson's locking along the in‐plane and out‐of‐plane directions is overcome by using the hybrid‐stress modes of the first element. The third element passes the patch test for plate/shell elements. Unless the element assumes the right prismatic geometry, it fails the patch test for three‐dimensional continuum elements. It will be seen that all the proposed elements are markedly more accurate than the conventional fully integrated element. Copyright © 2004 John Wiley & Sons, Ltd.     

Energy and momentum conserving elastodynamics of a non‐linear brick‐type mixed finite shell element

The paper presents aspects of the finite element formulation of momentum and energy conserving algorithms for the non‐linear dynamic analysis of shell‐like structures. The key contribution is a detailed analysis of the implementation of a Simó–Tarnow‐type conservation scheme in a recently developed new mixed finite shell element. This continuum‐based shell element provides a well‐defined interface to strain‐driven constitutive stress updates algorithms. It is based on the classic brick‐type trilinear displacement element and is equipped with specific gradient‐type enhanced strain modes and shell‐typical assumed strain modifications. The excellent performance of the proposed dynamic shell formulation with respect to conservation properties and numerical stability behaviour is demonstrated by means of three representative numerical examples of elastodynamics which exhibit complex free motions of flexible structures undergoing large strains and large rigid‐body motions. Copyright © 2001 John Wiley & Sons, Ltd.     

Efficient generation of pareto‐optimal topologies for compliance optimization

In multi‐objective optimization, a design is defined to beit pareto‐optimal if no other design exists that is better with respect to one objective, and as good with respect to other objectives. In this paper, we first show that if a topology is pareto‐optimal, then it must satisfy certain properties associated with the topological sensitivity field, i.e. no further comparison is necessary. This, in turn, leads to a deterministic, i.e. non‐stochastic, method for efficiently generating pareto‐optimal topologies using the classic fixed‐point iteration scheme. The proposed method is illustrated, and compared against SIMP‐based methods, through numerical examples. In this paper, the proposed method of generating pareto‐optimal topologies is limited to bi‐objective optimization, namely compliance–volume and compliance–compliance. The future work will focus on extending the method to non‐compliance and higher dimensional pareto optimization. Copyright © 2011 John Wiley & Sons, Ltd.     

Magnetic resonance imaging reconstruction via non‐convex total variation regularization

Magnetic resonance imaging (MRI) reconstruction model based on total variation (TV) regularization can deal with problems such as incomplete reconstruction, blurred boundary, and residual noise. In this article, a non‐convex isotropic TV regularization reconstruction model is proposed to overcome the drawback. Moreau envelope and minmax‐concave penalty are firstly used to construct the non‐convex regularization of L₂ norm, then it is applied into the TV regularization to construct the sparse reconstruction model. The proposed model can extract the edge contour of the target effectively since it can avoid the underestimation of larger nonzero elements in convex regularization. In addition, the global convexity of the cost function can be guaranteed under certain conditions. Then, an efficient algorithm such as alternating direction method of multipliers is proposed to solve the new cost function. Experimental results show that, compared with several typical image reconstruction methods, the proposed model performs better. Both the relative error and the peak signal‐to‐noise ratio are significantly improved, and the reconstructed images also show better visual effects. The competitive experimental results indicate that the proposed approach is not limited to MRI reconstruction, but it is general enough to be used in other fields with natural images.     

Projection‐based model reduction for contact problems

To be feasible for computationally intensive applications such as parametric studies, optimization, and control design, large‐scale finite element analysis requires model order reduction. This is particularly true in nonlinear settings that tend to dramatically increase computational complexity. Although significant progress has been achieved in the development of computational approaches for the reduction of nonlinear computational mechanics models, addressing the issue of contact remains a major hurdle. To this effect, this paper introduces a projection‐based model reduction approach for both static and dynamic contact problems. It features the application of a non‐negative matrix factorization scheme to the construction of a positive reduced‐order basis for the contact forces, and a greedy sampling algorithm coupled with an error indicator for achieving robustness with respect to model parameter variations. The proposed approach is successfully demonstrated for the reduction of several two‐dimensional, simple, but representative contact and self contact computational models. Copyright © 2015 John Wiley & Sons, Ltd.     

Shifted FSAI preconditioners for the efficient parallel solution of non‐linear groundwater flow models

The present paper investigates the performance of a shifted factorized sparse approximate inverse as a parallel preconditioner for the iterative solution to the linear systems arising in the finite element discretization of non‐linear groundwater flow models. The shift strategy is based on an inexpensive preconditioner update exploiting the structure of the coefficient matrix. The proposed algorithm is experimented with in the parallel simulation of a large‐scale real multi‐aquifer system characterized by a stochastic distribution of the hydraulic conductivity. The numerical results show that the shifted factorized sparse approximate inverse algorithm may yield an overall computational gain up to 300% with respect to the non‐shifted scheme with an excellent parallel efficiency. Copyright © 2011 John Wiley & Sons, Ltd.     

An extension of Hill's three‐component model to include different fibre types in finite element modelling of muscle

Most of the proposed versions of the Hill's model use a sliding‐element theory, considering a single sarcomere. However, a muscle represents a collection of different fibre types with a large range in contractile properties among them. An extension of Hill's three‐component model is proposed here to take into account different fibre types. We present a model consisting of a number of sarcomeras of different types coupled in parallel with the connective tissue. Each sarcomere is modelled by one non‐linear elastic element connected in series with one non‐linear contractile element. Using the finite element method, in an incremental‐iterative scheme of calculating equilibrium configurations of a muscle, the key step is the determination of stresses corresponding to strain increments. The stress calculation procedure for the extended Hill's model is reduced to the solution of a number of independent non‐linear equations with respect to the stretch increments of the serial elastic elements in each sarcomere. Since the distribution of the specific fibre type is non‐uniform over the muscle volume, we have material heterogeneity which we modelled by using the so‐called ‘Generalized Isoparametric Element Formulation’ for functionally graded materials (FGMs). The proposed computational scheme is built in our FE package PAK, so that muscles of complex three‐dimensional shapes can be modelled. In numerical examples, we illustrate the main characteristics of the developed numerical model and some possibilities of realistic modelling of muscle functioning. Copyright © 2006 John Wiley & Sons, Ltd.     

Trust‐region based return mapping algorithm for implicit integration of elastic‐plastic constitutive models

A new method for the solution of the non‐linear equations forming the core of constitutive model integration is proposed. Specifically, the trust‐region method that has been developed in the numerical optimization community is successfully modified for use in implicit integration of elastic‐plastic models. Although attention here is restricted to these rate‐independent formulations, the proposed approach holds substantial promise for adoption with models incorporating complex physics, multiple inelastic mechanisms, and/or multiphysics. As a first step, the non‐quadratic Hosford yield surface is used as a representative case to investigate computationally challenging constitutive models. The theory and implementation are presented, discussed, and compared with other common integration schemes. Multiple boundary value problems are studied and used to verify the proposed algorithm and demonstrate the capabilities of this approach over more common methodologies. Robustness and speed are then investigated and compared with existing algorithms. Through these efforts, it is shown that the utilization of a trust‐region approach leads to superior performance versus a traditional closest‐point projection Newton–Raphson method and comparable speed and robustness to a line search augmented scheme. Copyright © 2017 John Wiley & Sons, Ltd.     

Improved many‐to‐one sweeping

This paper presents PMOST, the polymorphic many‐to‐one sweep tool. Through its polymorphic design it is able to project nodes quickly and with high quality through ‘modular’ node projection schemes. PMOST currently employs two different node projection schemes: Faceted and BoundaryError node projectors. Implementation details are presented for both of these methods. The proposed many‐to‐one sweep tool was compared to existing software with respect to speed and element quality. Performance testing using PMOST showed that the Faceted node projection scheme was the most efficient while the BoundaryError method produced less element distortion while approximating the efficiency of the Faceted scheme. We also demonstrate improved performance with respect to both element quality and algorithm efficiency when compared to other existing node projection methods. Published in 2005 by John Wiley & Sons, Ltd.     

A two‐grid method for expanded mixed finite‐element solution of semilinear reaction–diffusion equations

We present a scheme for solving two‐dimensional semilinear reaction–diffusion equations using an expanded mixed finite element method. To linearize the mixed‐method equations, we use a two‐grid algorithm based on the Newton iteration method. The solution of a non‐linear system on the fine space is reduced to the solution of two small (one linear and one non‐linear) systems on the coarse space and a linear system on the fine space. It is shown that the coarse grid can be much coarser than the fine grid and achieve asymptotically optimal approximation as long as the mesh sizes satisfy H=O(h^1/3). As a result, solving such a large class of non‐linear equation will not be much more difficult than solving one single linearized equation. Copyright © 2003 John Wiley & Sons, Ltd.     

A regularized force‐based beam element with a damage–plastic section constitutive law

A new beam finite element is presented, with a generalized section constitutive law based on damage mechanics and plasticity, to analyse the cyclic structural response of plane frames. Both displacement‐based and force‐based (FB) approaches are used and compared, to demonstrate the significant advantages of the FB formulation in the presence of material non‐linearity. In order to overcome the analytical problems and the pathological mesh dependency of the numerical response in the presence of strain‐softening post‐peak behaviour, a classical non‐local regularization procedure is adopted first, based on the integral definition of the associated variable governing the damaging evolution process. Subsequently, for the FB element a new simple regularization technique is proposed based on a selected integration procedure along the element length, which predefines the location of the Gauss points in the beam region, where the localization phenomena take place. As for the other computational aspects, an iterative element state determination is adopted for the FB formulation and a local predictor–corrector algorithm is used to solve the incremental evolution problems of the damage and plastic internal variables. Finally, some examples are shown on simple beams and frames, subjected to monotonically increasing and cyclic loading conditions. Copyright © 2006 John Wiley & Sons, Ltd.     

A fast implementation of the FETI‐DP method: FETI‐DP‐RBS‐LNA and applications on large scale problems with localized non‐linearities

As parallel and distributed computing gradually becomes the computing standard for large scale problems, the domain decomposition method (DD) has received growing attention since it provides a natural basis for splitting a large problem into many small problems, which can be submitted to individual computing nodes and processed in a parallel fashion. This approach not only provides a method to solve large scale problems that are not solvable on a single computer by using direct sparse solvers but also gives a flexible solution to deal with large scale problems with localized non‐linearities. When some parts of the structure are modified, only the corresponding subdomains and the interface equation that connects all the subdomains need to be recomputed. In this paper, the dual–primal finite element tearing and interconnecting method (FETI‐DP) is carefully investigated, and a reduced back‐substitution (RBS) algorithm is proposed to accelerate the time‐consuming preconditioned conjugate gradient (PCG) iterations involved in the interface problems. Linear–non‐linear analysis (LNA) is also adopted for large scale problems with localized non‐linearities based on subdomain linear–non‐linear identification criteria. This combined approach is named as the FETI‐DP‐RBS‐LNA algorithm and demonstrated on the mechanical analyses of a welding problem. Serial CPU costs of this algorithm are measured at each solution stage and compared with that from the IBM Watson direct sparse solver and the FETI‐DP method. The results demonstrate the effectiveness of the proposed computational approach for simulating welding problems, which is representative of a large class of three‐dimensional large scale problems with localized non‐linearities. Copyright © 2005 John Wiley & Sons, Ltd.     

A finite difference solution of non‐linear systems of radiative–conductive heat transfer equations

A finite difference solution for a system of non‐linear integro–differential equations modelling the steady‐state combined radiative–conductive heat transfer is proposed. A new backward–forward finite difference scheme is formulated for the Radiative Transfer Equation. The non‐linear heat conduction equation is solved using the Kirchhoff transformation associated with a centred finite difference scheme. The coupled system of equations is solved using a fixed‐point method, which relates to the temperature field. An application on a real insulator composed of silica fibres is illustrated. The results show that the method is very efficient. Copyright © 2002 John Wiley & Sons, Ltd.     

An essentially non‐oscillatory semi‐Lagrangian method for tidal flow simulations

We develop an essentially non‐oscillatory semi‐Lagrangian method for solving two‐dimensional tidal flows. The governing equations are derived from the incompressible Navier–Stokes equations with assumptions of shallow water flows including bed frictions, eddy viscosity, wind shear stresses and Coriolis forces. The method employs the modified method of characteristics to discretize the convective term in a finite element framework. Limiters are incorporated in the method to reconstruct an essentially non‐oscillatory algorithm at minor additional cost. The central idea consists in combining linear and quadratic interpolation procedures using nodes of the finite element where departure points are localized. The resulting semi‐discretized system is then solved by an explicit Runge–Kutta Chebyshev scheme with extended stages. This scheme adds in a natural way a stabilizing stage to the conventional Runge–Kutta method using the Chebyshev polynomials. The proposed method is verified for the recirculation tidal flow in a channel with forward‐facing step. We also apply the method for simulation of tidal flows in the Strait of Gibraltar. In both test problems, the proposed method demonstrates its ability to handle the interaction between water free‐surface and bed frictions. Copyright © 2009 John Wiley & Sons, Ltd.     

Fast BEM–FEM mortar coupling for acoustic–structure interaction

A coupling algorithm based on Lagrange multipliers is proposed for the simulation of structure–acoustic field interaction. Finite plate elements are coupled to a Galerkin boundary element formulation of the acoustic domain. The interface pressure is interpolated as a Lagrange multiplier, thus, allowing the coupling of non‐matching grids. The resulting saddle‐point problem is solved by an approximate Uzawa‐type scheme in which the matrix–vector products of the boundary element operators are evaluated efficiently by the fast multipole boundary element method. The algorithm is demonstrated on the example of a cavity‐backed elastic panel. Copyright © 2005 John Wiley & Sons, Ltd.