On an efficient octic order sub-parametric finite element method on curved domains

Abstract

An efficient finite element method using the sub-parametric transformations of the curved boundaries for the 45-node octic order triangular element is presented in this paper. The one-sided curved and two-sided straight triangle in the universal coordinate system (x,y) is mapped onto an isosceles triangle (ξ,η)∈[0,1] and ξ+η≤1 in the (ξ,η)-coordinate system by using an isoparametric coordinate transformation. Using such a transformation, the curved edge of the triangular element is altered implicitly by a newly defined octic curve. Such an octic curve is obtained by using the topology of the curved edge. With an apparent condition of arc being a parabola, a new transformation is developed to find a unique curve for each edge that passes through the points on the original curved arc. Overall, the new sub-parametric transformations are developed, these include the transformation for the curved edge on the original boundary, in addition to the parameters from the curved triangle's interior. Several numerical experiments are presented to show the efficacy of the proposed method. The computational results obtained for a boundary value problem using octic order basis functions outperform the results obtained using the existing higher order basis functions such as quadratic (6-node), cubic (10-node), quartic (16-node), quintic (21-node), sextic (28-node), and septic (36-node) order triangular elements (TEs) with curves. The method implemented in this paper can be easily extendable to quasi-static and dynamic crack evolution problems within linear and nonlinear elasticity.