How quadrilateral elements represented in most FEA solvers for structural mechanics problems?
I know that Nastran-95 represent quadrilaterals for membrane and bending elements as a set of triangles:

At the same time there are a lot of literature where stiffness matrix for quadrilaterals is computed using ξ-η coordinates:
The question is which one of both methods is used in modern FEA packages?