“…The integrand of the weak form of the boundary value problem is no longer polynomial; hence, the Gauss integration rule generally used in FEM is unsuitable. Several alternatives have been proposed to surpass such limitation (Fries and Belytschko, 2010), and one frequently used, and herein studied, is the subdivision of elements (Mo€ es et al, 1999;Sukumar et al, 2000;Park et al, 2009;Pereira et al, 2009;Torres et al, 2015;Chen et al, 2018;Zhou and Chen, 2019).…”