Разработан программный комплекс определения геометрических характе- ристик и анализа касательных напряжений призматических стержней про- извольного поперечного сечения на этапе препроцессорной и постпроцессор- ной обработки данных при конечно-элементном анализе. На основе принципа возможных перемещений получены вариационные функционалы для задачи Сен-Венана о кручении призматического стержня и об изгибе поперечной силой, не вызывающей закручивания, которые непосредственно используют- ся для получения разрешающих соотношений метода конечных элементов. На основе теоремы Бетти о взаимности работ определяются координаты центра изгиба. Все соотношения сформулированы относительно функции коробления, что позволяет избежать проблем, связанных с неоднознач- ностью при использовании функции напряжений Прандтля в случае много- связной области.