The in-plane element formulation accounts for Dx and Dy of the local coordinate system. The in-plane stiffness matrix of the element is based on the standard isoparametric formulation [Ref 1, 2, 3]. However, when the element is rectangular in shape, incompatible modes may be optionally added to the formulation [Ref. 3]. An incompatible element, when applied, yields results of high quality especially when used to model in-plane bending. Full two by two numerical integration is used to calculate the in-plane stiffness matrix of the element.