... a higher order of approximations, much better accuracy and continuous nodal stress. Moreover, the linear dependence problem associated with many PU-based methods with high-order approximations is eliminated in the present element. A number of numerical examples indicate the high accuracy and robustness...