We propose an iterationless approach to studying subcritical and supercritical deformations of shells of complex shape having anisotropic layers. We obtain the relationships which are necessary for derivation of general matrices and vectors when using the finite element method and describe the procedure of recalculation of the stress state when organizing a step-by-step process. We present the results of solution of some testing problems.