NODAL ACCURACY IMPROVEMENT AND SUPER-CONVERGENT COMPUTATION IN FEM ANALYSIS OF FEMOL SECOND ORDER ODES

HUANG Ze-min,YUAN Si

DOI: https://doi.org/10.6052/j.issn.1000-4750.2021.06.S002

2022-01-01

Abstract:Elements with degree \begin{document}$ m $\end{document} is used in finite element method (FEM) to solve the second order ordinary differential equations (ODEs) derived from the FEM of lines (FEMOL). The interior displacement of elements generally has a convergence order of \begin{document}$ m + 1 $\end{document}, while the nodal displacements can achieve a convergence order of \begin{document}$ 2m $\end{document}. The super-convergence computation using the element energy projection (EEP) method usually has a convergence order of \begin{document}$ \min (m + 2,2m) $\end{document}, which benefits from the nodal displacements of a higher convergence order but also limits its accuracy by the nodal displacements of elements with lower degrees. In this paper, a modified EEP (M-EEP) method is proposed. With the EEP solution, the nodal displacement accuracy is improved first, and then the interior displacement of elements is recovered, which leads to a modified EEP solution. Numerical experiments show that improved nodal displacements can achieve a convergence order of \begin{document}$ 2m + 2 $\end{document}, and the interior displacements of elements always have a convergence order of \begin{document}$ m + 2 $\end{document} without the constraint of order \begin{document}$ 2m $\end{document}. For linear elements, the interior displacement of M-EEP solution does not have the limitation of second-order convergence from the traditional FEM solution and can achieve the remarkable third-order convergence, equivalent to the convergence order of quadratic elements.