Использован криволинейный конечный элемент срединной линии осесимметрично нагруженной оболочки вращения с матрицей жесткости размером $8\times 8$ при выборе узловых неизвестных в виде перемещений и их первых производных. Определяющие уравнения на шаге нагружения реализованы в двух вариантах. В первом варианте использованы соотношения деформационной теории пластичности, состоящие из выражений упругих и пластических частей. Соотношения между приращениями деформаций и приращениями напряжений определялись дифференцированием используемых уравнений.
