A higher-order zig-zag theory for laminated composite and sandwich structures is proposed. The proposed theory satisfies the interlaminar continuity conditions and free surface conditions of transverse shear stresses. Moreover, the number of unknown variables involved in present model is independent of the number of layers. Compared to the zig-zag theory available in literature, the merit of present theory is that the first derivatives of transverse displacement have been taken out from the in-plane displacement fields, so that the C 0 interpolation functions is only required during its finite element implementation. To obtain accurately transverse shear stresses by integrating three-dimensional equilibrium equations within one element, a six-node triangular element is employed to model the present zig-zag theory. Numerical results show that the present zig-zag theory can predict more accurate in-plane displacements and stresses in comparison with other zig-zag theories. Moreover, it is convenient to obtain transverse shear stresses by integration of equilibrium equations, as the C 0 shape functions is only used when implemented in a finite element.