Shape functions for concave quadrilaterals