Разработаны расчетные методы и созданы программные средства для численного определения НДС и усилий противоизгиба холоднолистового четырехвалкового прокатного стана с учетом остаточных термонапряжений. Расчет термонапряжений осуществлен посредством решения задачи термоупруговязкопластичности для среды с нестабильной структурой. В основу решения объемной контактной задачи для определения рабочих напряжений от прокатки положен конечно-элементный комплекс ANSYS. Представлены результаты численной реализации применительно к реальному прокатному стану КВАРТО 500/1300х1700. Установлено, что приближенное аналитическое решение, основанное на балочной схеме, дает заниженные значения усилий противоизгиба, а также что остаточные напряжения от термообработки валков значительно больше эксплуатационных напряжений от прокатки.