Purpose This study aims to propose a general and efficient adaptive strategy with local mesh refinement for two-dimensional (2D) finite element (FE) analysis based on the element energy projection (EEP) technique. Design/methodology/approach In view of the inflexibility of the existing global dimension-by-dimension (D-by-D) recovery method via EEP technique, in which displacements are recovered through element strips, an improved element D-by-D recovery strategy was proposed, which enables the EEP recovery of super-convergent displacements to be implemented mostly on a single element. Accordingly, a posteriori error estimate in maximum norm was established and an EEP-based adaptive FE strategy of h-version with local mesh refinement was developed. Findings Representative numerical examples, including stress concentration and singularity problems, were analyzed; the results of which show that the adaptively generated meshes reasonably reflect the local difficulties inherent in the physical problems and the proposed adaptive analysis can produce FE displacement solutions satisfying the user-specified tolerances in maximum norm with an almost optimal adaptive convergence rate. Originality/value The proposed element D-by-D recovery method is a more efficient and flexible displacement recovery method, which is implemented mostly on a single element. The EEP-based adaptive FE analysis can produce displacement solutions satisfying the specified tolerances in maximum norm with an almost optimal convergence rate and thus can be expected to apply to other 2D problems.

Abstract:The finite element (FE) method has been widely used in predicting the structural responses of asphalt pavements. However, the three-dimensional (3D) modeling in general-purpose FE software systems such as ABAQUS requires extensive computations and is relatively time-consuming. To address this issue, a specific computational code EasyFEM was developed based on the finite layer method (FLM) for analyzing structural responses of asphalt pavements under a static load. Basically, it is a 3D FE code that requires only a one-dimensional (1D) mesh by incorporating analytical methods and using Fourier series in the other two dimensions, which can significantly reduce the computational time and required resources due to the easy implementation of parallel computing technology. Moreover, a newly-developed Element Energy Projection (EEP) method for super-convergent calculations was implemented in EasyFEM to improve the accuracy of solutions for strains and stresses over the whole pavement model. The accuracy of the program is verified by comparing it with results from BISAR and ABAQUS for a typical asphalt pavement structure. The results show that the predicted responses from ABAQUS and EasyFEM are in good agreement with each other. The EasyFEM with the EEP post-processing technique converges faster compared with the results derived from ordinary EasyFEM applications, which proves that the EEP technique can improve the accuracy of strains and stresses from EasyFEM. In summary, the EasyFEM has a potential to provide a flexible and robust platform for the numerical simulation of asphalt pavements and can easily be post-processed with the EEP technique to enhance its advantages.

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.