A method for numerical analysis of the elastoplastic stress-strain state of thin layered shells of revolution under axisymmetric loading is proposed. Constitutive equations describing the elastoplastic deformation of isotropic materials with allowance for the stress mode are used. Numerical results are presented Keywords: constitutive equations, elastoplastic deformation, isotropic material, stress-strain state Introduction. The paper [1] outlines a method of successive approximations for solving boundary-value problems of plasticity using constitutive equations [12][13][14] that describe the elastoplastic deformation of isotropic materials along paths of small curvature with allowance for the stress mode. The practical convergence of the method was tested by way of an example. The stress state is characterized in [12][13][14] by the stress mode angle [3], which is defined as the angle between the direction of the octahedral shear stress and the negative direction of the projection of the principal axis (along which the principal stresses are minimum) of the stress deviator onto the octahedral plane. In [2], a method of solving the axisymmetric problem of the theory of thin shells is outlined and numerical results obtained with the same constitutive equations are presented. These equations relating the components of the stress tensors and the linear strain components (to be called strains, for short) were experimentally validated in [8,9,[12][13][14]. It was assumed that the strains can be represented as the sum of elastic and plastic components and that the stress deviators and the deviators of plastic strain differential are coaxial. The equations include two nonlinear functions found experimentally. One of these functions relates the mean stress, the mean strain, and the stress mode angle. The other function relates the shear-stress intensity, the shear-strain intensity, and the stress mode angle. If we replace the former nonlinear function by a linear relation between the first invariants of the stress and strain tensors and assume that the latter function is independent of the stress mode and determined from simple-tension tests, then the constitutive equations go over into the widely used equations [5-7, etc.] describing deformation along paths of small curvature [5], which in the case of active loading coincide with the equations of incremental plasticity [3, 10, etc.] associated with the von Mises yield criterion (Prandtl-Reuss equations).Simplified constitutive equations [12][13][14] were considered in [9], where the former nonlinear function was replaced by a linear relation between the mean stress and mean strain and the possibility of such a replacement was justified. The simplified constitutive equations [9] were used in [11] to solve a spatial problem of plasticity for a body of revolution based on a method of successive approximations.In contrast to [11] and support of [1, 2], the present paper outlines a method to solve axisymmetric problems of plasticity for thin layered shells allowing for the ...