Multilayered structures, such as sandwich panels, composite laminated structures are extensively used to build large part of next generation aircraft, spacecraft and advanced components of ship and automotive vehicles.Analysis of failure mechanisms in layered structures demands an accurate evaluation of strain and stress fields in each lamina.Advanced shell models with variable kinematic as well as zig-zag theories have been proposed in last three decades literature.However, the solutions of real structures with complex geometries and boundary conditions requires the use of computational methods.Among these numerical methods, a relevant role is played by finite element method, FEM [1].On the other hand, numerical models could introduce numerical mechanisms that put severe limitations on the use of well established continuum mechanic theories.This is the case of the shell formulation which demands specific efforts to overcome shear and membrane locking [2].Efficient shell elements are known for linear and nonlinear analysis of traditional shell structures made by isotropic one layered materials.