A p-adaptive scaled boundary finite element method based on maximization of the error decrease rate

This study enhances the classical energy norm based adaptive procedure by introducing new refinement criteria, based on the projection-based interpolation technique and the steepest descent method, to drive mesh refinement for the scaled boundary finite element method. The technique is applied to p-adaptivity in this paper, but extension to h- and hp-adaptivity is straightforward. The reference solution, which is the solution of the fine mesh formed by uniformly refining the current mesh, is used to represent the unknown exact solution. In the new adaptive approach, a projection-based interpolation technique is developed for the 2D scaled boundary finite element method. New refinement criteria are proposed. The optimum mesh is assumed to be obtained by maximizing the decrease rate of the projection-based interpolation error appearing in the current solution. This refinement strategy can be interpreted as applying the minimisation steepest descent method. Numerical studies show the new approach out-performs the conventional approach.     

An a posteriori error estimator for the p‐ and hp‐versions of the finite element method

An a posteriori error estimator is proposed in this paper for the p‐ and hp‐versions of the finite element method in two‐dimensional linear elastostatic problems. The local error estimator consists in an enhancement of an error indicator proposed by Bertóti and Szabó (Int. J. Numer. Meth. Engng. 1998; 42 :561–587), which is based on the minimum complementary energy principle. In order to obtain the local error estimate, this error indicator is corrected by a factor which depends only on the polynomial degree of the element. The proposed error estimator shows a good effectivity index in meshes with uniform and non‐uniform polynomial distributions, especially when the global error is estimated. Furthermore, the local error estimator is reliable enough to guide p‐ and hp‐adaptive refinement strategies. Copyright © 2004 John Wiley & Sons, Ltd.     

Adaptive selection of polynomial degrees on a finite element mesh

The problem of finding a nearly optimal distribution of polynomial degrees on a fixed finite element mesh is discussed. An a posteriori error estimator based on the minimum complementary energy principle is proposed which utilizes the displacement vector field computed from the finite element solution. This estimator, designed for p- and hp-extensions, is conceptually different from estimators based on residuals or patch recovery which are designed for h-extension procedures. The quality of the error estimator is demonstrated by examples. The results show that the effectivity index is reasonably close to unity and the sequences of p-distributions obtained with the error indicators closely follow the optimal trajectory. © 1998 John Wiley & Sons, Ltd.     

An h‐hierarchical adaptive procedure for the scaled boundary finite‐element method

The scaled boundary finite‐element method (a novel semi‐analytical method for solving linear partial differential equations) involves the solution of a quadratic eigenproblem, the computational expense of which rises rapidly as the number of degrees of freedom increases. Consequently, it is desirable to use the minimum number of degrees of freedom necessary to achieve the accuracy desired. Stress recovery and error estimation techniques for the method have recently been developed. This paper describes an h‐hierarchical adaptive procedure for the scaled boundary finite‐element method. To allow full advantage to be taken of the ability of the scaled boundary finite‐element method to model stress singularities at the scaling centre, and to avoid discretization of certain adjacent segments of the boundary, a sub‐structuring technique is used. The effectiveness of the procedure is demonstrated through a set of examples. The procedure is compared with a similar h‐hierarchical finite element procedure. Since the error estimators in both cases evaluate the energy norm of the stress error, the computational cost of solutions of similar overall accuracy can be compared directly. The examples include the first reported direct comparison of the computational efficiency of the scaled boundary finite‐element method and the finite element method. The scaled boundary finite‐element method is found to reduce the computational effort considerably. Copyright © 2002 John Wiley & Sons, Ltd.     

Guaranteed energy error bounds for the Poisson equation using a flux‐free approach: Solving the local problems in subdomains

A method to compute guaranteed upper bounds for the energy norm of the exact error in the finite element solution of the Poisson equation is presented. The bounds are guaranteed for any finite element mesh however coarse it may be, not just in the asymptotic regime. The bounds are constructed by employing a subdomain‐based a posteriori error estimate which yields self‐equilibrated residual loads in stars (patches of elements). The proposed approach is an alternative to standard equilibrated residual methods providing sharper bounds. The use of a flux‐free error estimator improves the effectivities of the upper bounds for the energy while retaining the certainty of the bounds. Copyright © 2009 John Wiley & Sons, Ltd.     

An automatic adaptive refinement finite element procedure for 2D elastostatic analysis

An automatic adaptive refinement procedure for finite element analysis is presented. The procedure is applied to two-dimensional elastostatic problems to obtain solutions of prescribed accuracy. Through the combined use of new mesh generator using contour developed by Lo¹ and the concept of strain energy concentration, high-quality graded finite element meshes are generated. The whole process is fully automatic and no user intervention is required during the successive cycles of the mesh refinements. The Zienkiewicz and Zhu² error estimator is found to be effective and has been adopted for the present implementation. In the numerical examples tested, the error estimator gives an accurate error norm estimation and the effectivity index of the estimator converges to a value close to unity.     

EXTENSION OF THE ZIENKIEWICZ–ZHU ERROR ESTIMATOR TO SHAPE SENSITIVITY ANALYSIS

The application of the Zienkiewicz–Zhu estimator was extended to the estimation of the discretization error arising from shape sensitivity analysis using the finite element method. The sensitivity error was quantified from the sensitivity of the energy norm by using an estimator specially developed for this purpose. Sensitivity analyses were carried out using the discrete analytical approach, which introduced no additional errors other than the discretization error. In this work, direct nodal averaging was used for linear triangular elements and the SPR technique for quadratic elements in order to obtain the smoothed stress and sensitivities fields. Two examples with an exact solution are used to analyse the effectivity of the proposed estimator and its convergence with the h-adaptive refinement. © 1997 by John Wiley & Sons, Ltd.     

Computation of Single Eigenfrequencies and Eigenfunctions of Plate and Shell Structures Using an H-adaptive FE-method

The accuracy of the computation of eigenfrequencies and eigenfunctions with FE-methods can be substantially improved with efficient adaptive procedures. For such an adaptive analysis of plate and shell structures a simple a-posteriori error estimator or indicator for the error in the energy norm and L ₂-norm of the eigenfunction u _h for shell and plate structures is proposed.Both indicators hite shall represent the correct convergence of the estimated error. The estimator for the error in the energy norm is used to enlarge adaptively the dimension of the finite element subspace. On the basis of numerical examples the efficiency and the quality of the improved solution is discussed. In order to validate the quantity of the estimated error Aitken’s extrapolation technique is applied.     

Mesh adaptation using a four-noded quadrilateral plate bending element

This paper uses a four-noded quadrilateral Reissner-Mindlin plate bending element for mesh adaptation. To overcome the problems of “locking” in the thin plate limit and zero energy modes the transverse shear strains are not evaluated from the displacements and rotations but are separately interpolated. The element is used in a hierarchical mesh adaptation which uses a node-averaging-based energy norm error estimator. Several examples illustrate the use of the adaptive algorithm.     

An hp‐adaptive finite element method for scattering problems in computational electromagnetics

An hp‐adaptive finite element (FE) approach is presented for a reliable, efficient and accurate solution of 3D electromagnetic scattering problems. The radiation condition in the far field is satisfied automatically by approximation with infinite elements (IE). Near optimal discretizations that can effectively resolve local rapid variations in the scattered field are sought adaptively by mesh refinements blended with graded polynomial enrichments. The p‐enrichments need not be spatially isotropic. The discretization error can be controlled by a self‐adaptive process, which is driven by implicit or explicit a posteriori error estimates. The error may be estimated in the energy norm or in a quantity of interest. A radar cross section (RCS) related linear functional is used in the latter case. Adaptively constructed solutions are compared to pure uniform p approximations. Numerical, highly accurate, and fairly converged solutions for a number of generic problems are given and compared to previously published results. Copyright © 2004 John Wiley & Sons, Ltd.     

Adaptive Poly-FEM for the analysis of plane elasticity problems

In this work, we present an adaptive polygonal finite element method (Poly-FEM) for the analysis of two-dimensional plane elasticity problems. The generation of meshes consisting of n ? sided polygonal finite elements is based on the generation of a centroidal Voronoi tessellation (CVT). An unstructured tessellation of a scattered point set, that minimally covers the proximal space around each point in the point set, is generated whereby the method also includes tessellation of nonconvex domains. In this work, we propose a region by region adaptive polygonal element mesh generation. A patch recovery type of stress smoothing technique that utilizes polygonal element patches for obtaining smooth stresses is proposed for obtaining the smoothed finite element stresses. A recovery type a ? posteriori error estimator that estimates the energy norm of the error from the recovered solution is then adopted for the Poly-FEM. The refinement of the polygonal elements is then made on an region by region basis through a refinement index. For the numerical integration of the Galerkin weak form over polygonal finite element domains, we resort to classical Gaussian quadrature applied to triangular subdomains of each polygonal element. Numerical examples of two-dimensional plane elasticity problems are presented to demonstrate the efficiency of the proposed adaptive Poly-FEM.     

Stress recovery and error estimation for the scaled boundary finite‐element method

The scaled boundary finite‐element method is a novel semi‐analytical technique, combining the advantages of the finite element and the boundary element methods with unique properties of its own. This paper develops a stress recovery procedure based on a modal interpretation of the scaled boundary finite‐element method solution process, using the superconvergent patch recovery technique. The recovered stresses are superconvergent, and are used to calculate a recovery‐type error estimator. A key feature of the procedure is the compatibility of the error estimator with the standard recovery‐type finite element estimator, allowing the scaled boundary finite‐element method to be compared directly with the finite element method for the first time. A plane strain problem for which an exact solution is available is presented, both to establish the accuracy of the proposed procedures, and to demonstrate the effectiveness of the scaled boundary finite‐element method. The scaled boundary finite‐element estimator is shown to predict the true error more closely than the equivalent finite element error estimator. Unlike their finite element counterparts, the stress recovery and error estimation techniques work well with unbounded domains and stress singularities. Copyright © 2002 John Wiley & Sons, Ltd.     

A posteriori error estimation for extended finite elements by an extended global recovery

This contribution presents an extended global derivative recovery for enriched finite element methods (FEMs), such as the extended FEM along with an associated error indicator. Owing to its simplicity, the proposed scheme is ideally suited to industrial applications. The procedure is based on global minimization of the L₂ norm of the difference between the raw strain field (C_?1) and the recovered (C₀) strain field. The methodology engineered in this paper extends the ideas of Oden and Brauchli (Int. J. Numer. Meth. Engng 1971; 3 ) and Hinton and Campbell (Int. J. Numer. Meth. Engng 1974; 8 ) by enriching the approximation used for the construction of the recovered derivatives (strains) with the gradients of the functions employed to enrich the approximation employed for the primal unknown (displacements). We show linear elastic fracture mechanics examples, both in simple two‐dimensional settings, and for a three‐dimensional structure. Numerically, we show that the effectivity index of the proposed indicator converges to unity upon mesh refinement. Consequently, the approximate error converges to the exact error, indicating that the error indicator is valid. Additionally, the numerical examples suggest a novel adaptive strategy for enriched approximations in which the dimensions of the enrichment zone are first increased, before standard h‐ and p‐adaptivities are applied; we suggest to coin this methodology e‐adaptivity. Copyright © 2008 John Wiley & Sons, Ltd.     

Error estimation and mesh adaptivity for the virtual element method based on recovery by compatibility in patches

In this article, a recovery by compatibility in patches (RCP)-based a posteriori error estimator is proposed for the virtual element method (VEM), and it is utilized to drive adaptive mesh refinement processes in two-dimensional elasticity problems. In RCP, recovered stresses are obtained by minimizing the complementary energy of patches of elements over a set of assumed equilibrated stress modes. To this aim, the explicit knowledge of displacements is only needed along the patch boundaries and no knowledge of superconvergent points is required, so making the RCP naturally suitable for the VEM. The a posteriori error estimation is conducted by comparing the stress field of a standard displacement-based VEM solution and the stress field obtained through RCP. The procedure is simple, and it does not require ad hoc modifications for small patches. The capability of this RCP-based error estimator to drive adaptive mesh refinements is successfully demonstrated through various numerical examples.     

An L∞–Lp mesh‐adaptive method for computing unsteady bi‐fluid flows

This paper discusses the contribution of mesh adaptation to high‐order convergence of unsteady multi‐fluid flow simulations on complex geometries. The mesh adaptation relies on a metric‐based method controlling the L ^p‐norm of the interpolation error and on a mesh generation algorithm based on an anisotropic Delaunay kernel. The mesh‐adaptive time advancing is achieved, thanks to a transient fixed‐point algorithm to predict the solution evolution coupled with a metric intersection in the time procedure. In the time direction, we enforce the equidistribution of the error, i.e. the error minimization in L ^∞ norm. This adaptive approach is applied to an incompressible Navier–Stokes model combined with a level set formulation discretized on triangular and tetrahedral meshes. Applications to interface flows under gravity are performed to evaluate the performance of this method for this class of discontinuous flows. Copyright © 2010 John Wiley & Sons, Ltd.     

Error estimation and adaptivity for discontinuous failure

Discontinuous failure is simulated via the introduction of a geometrical discontinuity. The cohesive zone is modelled via the use of the partition‐of‐unity property of the finite element interpolation. By this approach, a crack can pass through the elements without any restriction to the underlying mesh. Despite such a feature, it has been confirmed that a sufficiently fine mesh discretization still needs to be ensured in order to obtain a correct crack path and mechanical response. The p‐adaptive scheme, driven by error in an energy norm measure or a goal‐oriented measure, has been examined due to its implementational simplicity. The results have shown that, if considering only increasing the polynomial degree, the p‐approach can greatly improve the results. Copyright © 2008 John Wiley & Sons, Ltd.     

An adaptive finite element approach for frictionless contact problems

The derivation of an a posteriori error estimator for frictionless contact problems under the hypotheses of linear elastic behaviour and infinitesimal deformation is presented. The approximated solution of this problem is obtained by using the finite element method. A penalization or augmented‐Lagrangian technique is used to deal with the unilateral boundary condition over the contact boundary. An a posteriori error estimator suitable for adaptive mesh refinement in this problem is proposed, together with its mathematical justification. Up to the present time, this mathematical proof is restricted to the penalization approach. Several numerical results are reported in order to corroborate the applicability of this estimator and to compare it with other a posteriori error estimators. Copyright © 2001 John Wiley & Sons, Ltd.     

An adaptive remeshing procedure for discontinuous finite element limit analysis

An adaptive remeshing procedure is proposed for discontinuous finite element limit analysis. The procedure proceeds by iteratively adjusting the element sizes in the mesh to distribute local errors uniformly over the domain. To facilitate the redefinition of element sizes in the new mesh, the interelements discontinuous field of elemental bound gaps is converted into a continuous field, ie, the intensity of bound gap, using a patch‐based approximation technique. An analogous technique is subsequently used for the approximation of element sizes in the old mesh. With these information, an optimized distribution of element sizes in the new mesh is defined and then scaled to match the total number of elements specified for each iteration in the adaptive remeshing process. Finally, a new mesh is generated using the advancing front technique. This adaptive remeshing procedure is repeated several times until an optimal mesh is found. Additionally, for problems involving discontinuous boundary loads, a novel algorithm for the generation of fan‐type meshes around singular points is proposed explicitly and incorporated into the main adaptive remeshing procedure. To demonstrate the feasibility of our proposed method, some classical examples extracted from the existing literary works are studied in detail.     

Adaptive isogeometric analysis based on a combined r-h strategy

In the present work, an r-h adaptive isogeometric analysis is proposed for plane elasticity problems. For performing the r-adaption, the control net is considered to be a network of springs with the individual spring stiffness values being proportional to the error estimated at the control points. While preserving the boundary control points, relocation of only the interior control points is made by adopting a successive relaxation approach to achieve the equilibrium of spring system. To suit the noninterpolatory nature of the isogeometric approximation, a new point-wise error estimate for the h-refinement is proposed. To evaluate the point-wise error, hierarchical B-spline functions in Sobolev spaces are considered. The proposed adaptive h-refinement strategy is based on using De-Casteljau’s algorithm for obtaining the new control points. The subsequent control meshes are thus obtained by using a recursive subdivision of reference control mesh. Such a strategy ensures that the control points lie in the physical domain in subsequent refinements, thus making the physical mesh to exactly interpolate the control mesh and thereby allowing the exact imposition of essential boundary conditions in the classical isogeometric analysis (IGA). The combined r-h adaptive refinement strategy results in better convergence characteristics with reduced errors than r- or h-refinement. Several numerical examples are presented to illustrate the efficiency of the proposed approach.     

A posteriori error estimation in non‐linear FE analyses of shell structures

During the last decade, significant scientific efforts were made in the area of quality assurance of numerical results obtained by means of the finite element method (FEM). These efforts were based on adaptive remeshing controlled by an estimated error. This paper reports on the extension of error estimation to non‐linear shell analysis involving strain‐hardening and softening plasticity. In the context of incremental‐iterative analyses, an incremental error estimator is introduced. It is based on the rate of work. The stress recovery technique proposed by Zienkiewicz and Zhu (Int. J. Numer. Meth. Engng 1992; 33 :1331) is modified to allow for discontinuities of certain stress components in case of localization arising from, e.g. cracking of concrete. The developed error estimator is part of a calculation scheme for adaptive non‐linear FE analysis. If the estimated error exceeds a prespecified threshold value in the course of an adaptive analysis, a new mesh is generated. After mesh refinement the state variables are transferred from the old to the new mesh and the calculation is restarted at the load level which was attained with the old mesh. The performance of the proposed error estimator is demonstrated by means of adaptive calculations of a reinforced concrete (RC) cooling tower. The influence of the user‐prescribed error threshold on the numerical results is investigated. Copyright © 2002 John Wiley & Sons, Ltd.