On the basis of a cubic approximation of tangential displacements through the thickness of stack, we derive the nonlinear equations for multilayer shells. This method allows us to make more accurate estimation of tangential stress distribution. The main relationships are derived using the Hamilton variation principle.