Исследования в области повышения надежности технологических трубопроводов, работающих в нестационарных термомеханических условиях, трудоемкие и дорогостоящие, при условии проведения натурных испытаний. Использование достоверных и достаточно точных математических методов и алгоритмов позволяет создавать программные продукты, позволяющие осуществлять такие исследования с намного меньшими материальными и физическими затратами. Также применение ЭВМ для решения такого рода задач существенно экономит время, затрачиваемое на исследование.
Рассматривается тонкостенный цилиндрический элемент какой-либо конструкции под воздействием малоцикловых термомеханических нагрузок. Данный элемент может нагружаться внутренним давлением q и осевым усилием P при температурном воздействии T. Внешние нагрузки и температурное поле изменяются во времени независимо друг от друга и вызывают циклическое упругопластическое деформирование оболочки. Рассматривается работа конструкции в упругопластической стадии с учетом истории нагружения. При расчете учитывается возможность изменения знаков напряжений и пластических деформаций, а также физико-механических свойств конструкционного материала в процессе нагружения. Решение задачи строится в главных осях x, y, z рассматриваемой конструкции.
Полная система уравнений при расчете конструкций состоит из трех групп уравнений. Первую группу уравнений для статических задач составляют уравнения равновесия, которым должны удовлетворять компоненты напряженного состояния. На границе тела напряжения должны удовлетворять статическим граничным условиям.
Уравнения равновесия не зависят от физико-механических свойств материала и условий нагрева и являются оригинальными для каждой конструкции. Для расчета элементов конструкций шаговым методом эти уравнения и граничные условия записывают в приращениях.
,
Вторую группу уравнений составляет уравнение Лапласа:
,
где , , , , - приращения меридионального и кольцевого напряжений, внутреннего давления и осевого усилия на этапе нагружения; D, h - диаметр и толщина стенки трубопровода; и - главные радиусы кривизны.
Для цилиндрической элемента , а .
Третью группу уравнений составляют уравнения пластического течения, связывающие напряжения и деформации. Их конкретный вид определяется применяемой теорией пластичности и принятой для расчета математической моделью кинетики процесса упругопластического деформирования конструкционного материала. В данном случае рассматриваем деформирование линейно-анизотропного материала.
Уравнения упругопластического деформирования материала, связывающие приращения напряжений и деформаций, записываются в матричной форме:
= + , (1)
где матрицы и определяются выражениями:
= ;
=
= ,
= ,
;
H = .
Рассматривая деформации пластического течения, полагаем, что в пространстве девиаторов напряжений существует область, в пределах которой поведение материала упругое. Границы этой области определяют с заданным допуском поверхность текучести, конфигурация и положение которой являются функционалами процесса нагружения. Начальная поверхность текучести является поверхностью Мизеса. Уравнение поверхности текучести принимаем в форме [2]:
,
где:
, - активные напряжения;
, - компоненты девиатора напряжений;
.
= - приращение накопленной пластической деформации
- интенсивность активных напряжений.
Задача решается для конструкционного материала с линейным анизотропным упрочнением.
Предполагается, что программа нагружения конструкции задана. Определены также физико-механические характеристики конструкционного материала.
Программа нагружения конструкции разбивается на ряд малых этапов, расчет которых выполняется последовательно. Считается, что в общем случае программа нагружения имеет предысторию, и в момент начала численного анализа имеется полная информация о напряженно-деформированном состоянии конструкции.
Расчет конструкции шаговым методом включает две основных процедуры. Первая процедура связана с нахождением параметров напряженно-деформированного состояния для рассматриваемой конструкции.
В рассмотрение вводится вектор состояния
= ,
который полностью характеризует напряженно-деформированное состояние исследуемой конструкции.
Пусть в некоторой точке процесса, соответствующей началу очередного n-го этапа нагружения, все компоненты вектора состояния известны. Располагая этой информацией, вычисляется значения функций ( ), а также значения упругих постоянных материала в исследуемой точке процесса нагружения. Полагая, что все эти параметры сохраняют свои значения на протяжении n-го этапа нагружения, определяются коэффициенты разрешающей системы уравнений и находятся параметры напряженно-деформированного состояния.
В результате находятся приращения напряжений , , , деформаций , , на n-м этапе нагружения. Числовые значения искомых величин определяются на каждом этапе вычисления. На каждом этапе точке ставится в соответствие параметр plast (признак пластичности), который принимает значение 0, если в рассматриваемой точке материал деформируется упруго, или 1, если имеет место пластическое течение. В процессе нахождения параметров НДС конструкции формирование матрицы в уравнении (1) выполняется с учетом значения этого параметра:
= , если plast = 0;
= + , если plast = 1.
Найдя параметры напряженно-деформированного состояния конструкции, переходят к следующей процедуре расчета, связанной с определением векторов состояния в узловых точках конструкции в конце n-го этапа нагружения.
Вычисляется напряжения в конце n-го этапа по формулам:
, (j = x, y, z),
где - напряжения в начале n-го этапа нагружения.
Затем вычисляется значения активных напряжений в конце этапа нагружения:
, (j = x, y, z),
где , (j = x, y, z),
.
Дальнейший анализ параметров состояния различается для упругих и пластических узловых точек.
Для упругой точки параметр Plast=0. Для упругой точки можно предложить следующий алгоритм ее анализа и расчета:
1. Вычислить значение A:
A = < ,
где - радиус поверхности текучести;
2. Выяснить возможное состояние точки:
В данном случае узловая точка остается упругой. Признак пластичности сохраняется, продолжение расчета. {plast = 0; Down =0 ; Exit};
В этом случае рассматриваемая узловая точка переходит в пластическое состояние. Назначаем этой точке параметр пластичности Plast = 1 и производим перерасчет этапа нагружения. {plast = 1; Exit};
Рассматриваемая узловая точка вышла за пределы поверхности текучести. Необходимо уменьшение этапа нагружения, назначаем этой точке параметр пластичности Plast = 1 и повторно решаем краевую задачу. Коэффициент необходимого уменьшения нагрузки на этапе равен Down = 0,5.
Для пластической точки параметр Plast=1. Для пластической узловой точки можно предложить следующий алгоритм ее анализа и расчета:
1. Вычислить значение параметр по формуле
2. Вычислить скорость изменения параметра .
- < 0 (2)
Выполнение условия (2) означает упругую разгрузку. В этом случае необходимо назначить этой точке параметр пластичности Plast = 0 и произвести перерасчет этапа нагружения.
- ≥ 0 (3)
Если выполняется условие (3), то необходимо вычислить приращения:
,
= ,
= ,
= .
Для упругих точек приращение радиуса поверхности текучести определяется по формуле:
= .
- Определить вектор состояния в конце рассматриваемого этапа нагружения конструкции и перейти к расчету следующего этапа:
= {Z0} + {ΔZ}.
Схема рассмотренного алгоритма приведена на рисунке 1. Рассмотренный метод использован при разработке программного обеспечения компьютерного анализа упругопластического деформирования трубчатых элементов конструкций. Использование разработанного программного продукта для изучения особенностей поведения конструкции при сложных термомеханических воздействиях хорошо себя зарекомендовало, т.к. позволяет получать точные результаты, которые невозможно получить при расчетах ручным способом. При этом повышается надежность эксплуатации оборудования в целом, а значит, вероятность возникновения аварийных ситуаций сводится к минимуму, что, несомненно, актуально в современном мире. Когда гораздо дешевле и проще предотвратить катастрофу, чем бороться с ее последствиями. Постоянный мониторинг призван решить эту проблему, обеспечив своевременное информирование о возникающих неполадках.
Рис. 1. Схема алгоритма программы
Рецензенты:
- Луганцев Леонид Дмитриевич, доктор технических наук, профессор, заведующий кафедрой «Системы автоматизированного проектирования» Московского государственного университета инженерной экологии (МГУИЭ), г. Москва.
- Мокрова Наталья Владиславовна, доктор технических наук, профессор, заведующий кафедрой «Информационные технологии» Московского государственного университета инженерной экологии (МГУИЭ), г. Москва.