Задача расчета кинематических и силовых параметров процессов обработки металлов давлением является физически нелинейной и может быть записана в виде суммы линейного L(u) и N(u) нелинейного операторов:
Рассмотренный подход соответствует методу упругих решений А. А. Ильюшина, при котором неупругая задача сводится к последовательности упругих. Технологические задачи о больших деформациях при высокой температуре должны рассматриваться как задачи о течении нелинейно-вязкой жидкости. Поэтому разрешающие уравнения для задач обработки металлов давлением записываются аналогично тому, как это принято в гидродинамике:
В качестве характеристики вязкости деформируемого металла использован известный в теории обработки металлов давлением коэффициент жесткости:
Для решения упругой задачи, соответствующей рассматриваемой, применяется гранично—элементное уравнение:
В настоящем решении используется известная аналогия между течением вязкой жидкости и упруго-деформируемым несжимаемым телом (коэффициент Пуассона равен 0,5). Благодаря такому подходу величины Ui равны скоростям пластического течения в искомых точках на границе области. Для практического решения интегралы уравнения (5) заменяются суммами:
Входящие в уравнения (5) и (6) величины со звездочками полагаются искомыми, определяющиеся решением Кельвина. Для решения линейных задач используются только зависимости (5) и (6). В данном же случае, для решения нелинейной задачи пластической деформации используется итерационный цикл, подчиняющейся зависимости:
Один из способов реализации метода упругих решений, расчет на основе алгоритма с дополнительными массовыми силами и поправками к граничным условиям. В этом случае соответствующие зависимости имеют вид:
Отличие связано с тем, что в задачах пластического течения необходимо рассчитывать параметры не только на границе расчетной области но и внутри ее. Компоненты массовой силы в каждой внутренней ячейки области определяются как:
В работе расчеты выполнены с использованием линеаризующих алгоритмов с дополнительными массовыми силами и переменной вязкости. В этом случае на каждом последующем шаге вновь решается задача для того же упругого тела, но уже находящегося под действием новой нагрузки:
Укрупненный алгоритм решения задачи показан на рис. 1.
Расчет компонентов массовой силы производим в последовательности соответствующей блок-схеме показанной на рис.2
Последовательность вычислений по алгоритму переменной вязкости можно геометрически проиллюстрировать рис.3.
Рис. 3
Действительно, для каждого значения деформации находится значение отрезка CB = ω(εn)εn и в следующем приближении полагается, что
Как видно из рис. 3, процесс последовательных приближений (итерационный процесс) сходится к точному значению деформации εточн, если функция ω(ε) непрерывна и удовлетворяет условиям:
Однако можно заметить, что скорость сходимости существенно зависит от вида функции ω(ε). Если материал обладает большим упрочнением, т. е кривая мало отклоняется от прямой σ = Eε или, что то же самое, функция ω(ε) мала, то уже третье, четвертое приближение дает достаточно точное значение деформации.
И наоборот, если материал обладает малым упрочнением, то может потребоваться значительное число итераций (приближений), чтобы получить значение деформаций с требуемой точностью / 4 /.
Следует отметить. что рассмотренная методика является универсальной и пригодна для расчета поковок любой конфигурации.
Изучить методику расчетов процесса пластической деформации с применением метода граничных элементов и рассчитать перемещения при штамповке заданной поковки.
В качестве первого тестового примера используем задачу стяжении плоского кольца внутренним усилием рис.4
Размеры кольца: r = 20 мм и R = 40 мм, со следующими граничными условиями — нормальное напряжение на внутреннем контуре:
Касательные напряжения на поверхности контакта отсутствуют. На наружном контуре σn = σs = 0. Вследствие симметрии достаточно рассмотреть только четверть кольца. Граничные элементы, аппроксимируются четвертью кольца показанного на рис.4. а данные по ним приведены в табл. 1.
Таблица 1.
Как видно из граничным условий, эта задача в напряжениях. Кроме того, была выбрана одна линия полевых точек, находящаяся на оси симметрии х. В двадцати точках, лежащих на этой линии, были вычислимы смещения и напряжения, по которым, в свою очередь, была найдена массовая сила и поправки.
Полученные результаты хорошо согласуются с теоретическим расчетов выполненные методом линий скольжения / 2 /.
Во втором тестовом примере (рис. 4.1) наружный контур был описан 50 граничными элементами, а внутренняя область была разбита на 20 ячеек в которых вычислялись смещения, напряжения и массовые силы.
Рис 4.1
Трение на контакте было максимальным. Нормальные и касательные напряжения на свободной поверхности отсутствуют. Были наложены условия симметрии, что позволило упростить ввод данных и рассмотреть четверть заготовки. На третьем шаге расчетов была достигнута требуемая точность и сходимость итерационного цикла. Определяющим параметром было образование «бочки» на боковой поверхности (по 2мм) на сторону (рис. 4.2). Полученные результаты хорошо согласуются с результатами экспериментов Унксова /3/.
Рис 4.2
Фрагмент файла исходных данных для расчета задачи о полосе представлены ниже.
Сделать выводы по работе.