Математическая модель процесса вытяжки трубчатой заготовки

(50)

То есть в ячейки, соответствующие степеням свободы r и z данного узла, вносятся стоящие при них в уравнении (50) коэффициенты. Таким образом, удовлетворяется уравнение связи (49).

Рис. 3. Ограничения неподвижных границ

Подвижная граница перемещае

тся в заданном направлении (вертикальном или горизонтальном, в зависимости от вида процесса) на величину DH на каждом шаге:

(51)

где H – полный ход инструмента; i – число шагов решения.

Если после очередного перемещения границы какой-либо узел (или несколько) может оказаться “в теле” инструмента, то ему назначается принудительное перемещение на величину Dh:

или(52)

Здесь R, Z – координаты подвижной границы; r, z – координаты узла.

Таким образом, перемещение узла в направлении перпендикулярной границе степени свободы будет обеспечивать его движение вместе с границей, включая момент касания. В этом случае перемещение узла перпендикулярно границе не запрещается, а подчиняется равенствам:

или(53)

Для системы уравнений это означает внесение в соответствующие данному узлу ячейки матрицы жёсткости коэффициентов 0 и 1 в зависимости от направления движения и в правую часть величины Dh.

Рис. 4 Ограничения подвижных границ

Подвижная наклонная граница моделируется так же, как подвижная вертикальная или горизонтальная. Возникающая при этом ошибка, связанная со смещением узла вдоль границы минимизируется последовательными приближениями (см. Рис. 5).

Рис. 5. Процедура уточнения положения узла при смещении его подвижным инструментом

7. Трение

При пластическом формоизменении на границе контакта материала и инструмента возникает сила трения. При этом направление течения материала зависит от величины силы трения, то есть направление силы трения на границе контакта заранее не известно и может изменяться.

Модуль напряжения трения определяется законом Кулона (54):

,(54)

где sN – нормальное напряжение на границе инструмента; m – коэффициент трения скольжения; tK – касательное к границе напряжение .

Рис. 6. Контакт конечного элемента заготовки с инструментом

При пластическом течении касательное напряжение в элементах tK, контактирующих с внешними телами (матрица или пуансон), не должно превышать предел текучести при сдвиге tS:

(55)

Для выполнения этого условия предложен следующий алгоритм расчёта. Для элемента, лежащего на границе, определяется касательное напряжение tK, по которому вычисляется узловая сила трения:

(56)

где SK – площадь контакта элемента с границей инструмента. Если к узлу примыкают два элемента, лежащие на границе, то узловая сила трения получается суммированием сил, вычисленных по формуле (56). Для определения направления силы трения реализуется следующий алгоритм. Сначала выполняется шаг нагружения без учёта трения (FTP=0). Затем по результатам этого шага в каждом узле на границе контакта определяется нормальная узловая сила на границе контакта, приращение перемещения вдоль неё и, описанным выше способом, сила трения.

В процессе деформирования узлы, скользящие по границе инструмента, могут останавливаться и затем менять направление движения. В неподвижном состоянии узловая сила трения может оказаться меньше рассчитанной по формулам (54) – (56). Для этого введено последовательное уточнение силы FTP прибавлением к ней величины невязки DFTP:

(57)

Здесь n – номер итерации при уточнении силы трения.

Начальное значение невязки определяется силой трения:

(58)

Последующие значения невязки могут менять знак:

(59)

Эта формула обеспечивает рост силы трения в направлении противоположном движению. Если сила трения возрастает настолько, что меняет направление движения узла, то эта же формула обеспечит уменьшение силы трения.

Сходимость итерационного процесса обеспечивается уменьшением величины невязки при смене направления движения:

(60)

Кроме того, рост силы трения FTP ограничен её предельным значением , вычисленным по формулам (54) – (56). Описанная итерационная процедура приближает силу трения к значениям, изображённым на графике (Рис. 7).

Рис. 7. Зависимость силы трения от направления движения узла

Основные выводы

1. Приведенная математическая модель, сделанная на базе метода конечных элементов, является наиболее универсальной и адекватной с точки зрения оценки протекающих в ней процессов.

2. Модель подвижной наклонной границы позволяет более точно представить процесс, а также оценить картину деформирования и течения материала на радиусах скругления инструментов.

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

4. Введение в модель учета повреждаемости заготовки позволило не только проводить контроль за разрушением заготовки непосредственно во время протекания процесса, но и прогнозировать состояние изделия, полученного смоделированным методом.

 Скачать реферат