Adaptive finite element analysis of free vibration of elastic membranes via element energy projection technique

Purpose A general strategy is developed for adaptive finite element (FE) analysis of free vibration of elastic membranes based on the element energy projection (EEP) technique. Design/methodology/approach By linearizing the free vibration problem of elastic membranes into a series of linear equivalent problems, reliable a posteriori point-wise error estimator is constructed via EEP super-convergent technique. Hierarchical local mesh refinement is incorporated to better deal with tough problems. Findings Several classical examples were analyzed, confirming the effectiveness of the EEP-based error estimation and overall adaptive procedure equipped with a local mesh refinement scheme. The computational results show that the adaptively-generated meshes reasonably catch the difficulties inherent in the problems and the procedure yields both eigenvalues with required accuracy and mode functions satisfying user-preset error tolerance in maximum norm. Originality/value By reasonable linearization, the linear-problem-based EEP technique is successfully transferred to two-dimensional eigenproblems with local mesh refinement incorporated to effectively and flexibly deal with singularity problems. The corresponding adaptive strategy can produce both eigenvalues with required accuracy and mode functions satisfying user-preset error tolerance in maximum norm and thus can be expected to apply to other types of eigenproblems.

An adaptive nonlinear finite element analysis of minimal surface problem based on element energy projection technique

Purpose This paper aims to propose a new adaptive strategy for two-dimensional (2D) nonlinear finite element (FE) analysis of the minimal surface problem (MSP) based on the element energy projection (EEP) technique. Design/methodology/approach By linearizing nonlinear problems into a series of linear problems via the Newton method, the EEP technique, which is an effective and reliable point-wise super-convergent displacement recovery strategy for linear FE analysis, can be directly incorporated into the solution procedure. Accordingly, a posteriori error estimate in maximum norm was established and an adaptive 2D nonlinear FE strategy of h-version mesh refinement was developed. Findings Three classical known surfaces, including a singularity problem, were analysed. Moreover, an example whose analytic solution is unavailable was considered and a comparison was made between present results and those computed by the MATLAB PDE toolbox. The results show that the adaptively-generated meshes reflect the difficulties inherent in the problems and the proposed adaptive analysis can produce FE solutions satisfying the user-preset error tolerance in maximum norm with a fair adaptive convergence rate. Originality/value The EEP technique for linear FE analysis was extended to the nonlinear procedure of MSP and can be expected to apply to other 2D nonlinear problems. The employment of the maximum norm makes point-wisely error control on the sought surfaces possible and makes the proposed method distinguished from other adaptive FE analyses.

Adaptive Finite Element Method of Lines with Local Mesh Refinement in Maximum Norm Based on Element Energy Projection Method

The reliable and efficient self-adaptive analysis is a modern goal of various numerical computations. Most adaptivity methods, however, adopt energy norm to measure errors, which may not be the most natural and convenient means, e.g., for problems with locally singular gradient of displacement. Based on the Element Energy Projection (EEP) super-convergent technique in the Finite Element Method of Lines (FEMOL) which is a general and powerful semi-discrete method, reliable error estimates of displacements in maximum norm can be obtained anywhere on the FEMOL mesh and hence adaptive FEMOL by maximum norm becomes feasible. However, to tackle singularity problems effectively and efficiently, an automatic and flexible local mesh refinement strategy is required to generate meshes of high quality for more efficient adaptive FEMOL analysis. Taking the two-dimensional Poisson equation as the model problem, the paper firstly introduces the FEMOL and EEP methods with interface sides resulting from local mesh refinement. Then a local mesh refinement strategy and corresponding adaptive algorithm are presented. The numerical results given show that the proposed adaptive FEMOL with local mesh refinement can produce displacement solutions satisfying the specified tolerances in maximum norm and the adaptively-generated meshes reasonably reflect the local difficulties inherent in the physical problems without much redundant accuracy.

Adaptive finite element method for eigensolutions of regular second and fourth order Sturm-Liouville problems via the element energy projection technique

Purpose The purpose of this paper is to present a numerically adaptive finite element (FE) method for accurate, efficient and reliable eigensolutions of regular second- and fourth-order Sturm–Liouville (SL) problems with variable coefficients. Design/methodology/approach After the conventional FE solution for an eigenpair (i.e. eigenvalue and eigenfunction) of a particular order has been obtained on a given mesh, a novel strategy is introduced, in which the FE solution of the eigenproblem is equivalently viewed as the FE solution of an associated linear problem. This strategy allows the element energy projection (EEP) technique for linear problems to calculate the super-convergent FE solutions for eigenfunctions anywhere on any element. These EEP super-convergent solutions are used to estimate the FE solution errors and to guide mesh refinements, until the accuracy matches user-preset error tolerance on both eigenvalues and eigenfunctions. Findings Numerical results for a number of representative and challenging SL problems are presented to demonstrate the effectiveness, efficiency, accuracy and reliability of the proposed method. Research limitations/implications The method is limited to regular SL problems, but it can also solve some singular SL problems in an indirect way. Originality/value Comprehensive utilization of the EEP technique yields a simple, efficient and reliable adaptive FE procedure that finds sufficiently fine meshes for preset error tolerances on eigenvalues and eigenfunctions to be achieved, even on problems which proved troublesome to competing methods. The method can readily be extended to vector SL problems.

An Adaptive FEM for Buckling Analysis of Nonuniform Bernoulli-Euler Members via the Element Energy Projection Technique

This paper presents a new adaptive finite element (FE) procedure for accurate, efficient, and reliable computation of the critical buckling loads and the associated modes of nonuniform Bernoulli-Euler members. After the conventional FE solution on a given mesh has been obtained, a novel conceptual and practical strategy is introduced, in which the FE solution of the eigenproblem is equivalently viewed as the FE solution of an associated linear problem. This strategy allows the recently developed element energy projection (EEP) technique to be readily used to calculate super-convergent FE solutions for buckling modes, which are subsequently used to estimate the FE solution errors and to further guide mesh refinements, until the accuracy of the FE solution matches the user-preset error tolerance. Numerical examples are given to show the effectiveness, efficiency, accuracy, and reliability of the proposed method, which also paves the way for development of an adaptive FE solver for more general and complicated eigenproblems.