Numerical simulation of the axisymmetric hydrolipoclasy of thin metal sheet, fixed on the contour

Abstract

The results of application the axisymmetric rigid-plastic membrane finite element model to the study of forming a thin contour fixed metal plate under hydrostatic pressure are outlined. A comparison with the known numerical solutions and experimental data is provided.

Full Text

Достаточно распространенной в процессах холодной листовой штамповки является операция, при которой формоизменение закрепленной по контуру круглой заготовки из листового металла осуществляется под действием гидростатического давления (рисунок 1). В подобном процессе гидроформовки давление рабочей жидкости в течение некоторого промежутка времени растет, достигая своего предельного значения, а затем падает. При этом разрыв формуемой оболочки (в случае высокопластичного материала) обычно реализуется на запредельной стадии деформирования. Расчетным путем предельное и запредельное поведение формуемой указанным образом оболочки еще недостаточно изучено. Основная проблема здесь в обеспечении получения надежных числовых результатов при моделировании процесса формоизменения оболочки в запредельной области. Представленная в работах [1, 2] осесимметричная жесткопластическая безмоментная конечноэлементная модель обладает отмеченной вычислительной надежностью, о чем можно судить по результатам ее применения, изложенным в работах [1-5]. В настоящей статье вычислительная модель [1, 2] применяется к исследованию обсуждаемого процесса гидроформовки, включая и его запредельную стадию. Рисунок 1. Схема процесса свободного гидровыпучивания Краткое описание этой модели применительно к рассматриваемому случаю выглядит следующим образом. Предполагаем, что формуемая из листового металла под действием гидростатического давления осесимметричная оболочка является настолько тонкой, что ее допустимо считать безмоментной. Упругими деформациями на фоне больших пластических деформаций пренебрегаем, считая материал оболочки жесткопластическим. При этом используем предложенный Р. Хиллом [6] вариант теории течения (квадратичный критерий текучести) для трансверсально изотропного материала с изотропным упрочнением. Образующую рассматриваемой оболочки в ее исходном недеформированном состоянии разбиваем на достаточно большое количество участков малых размеров с тем, чтобы в течение всего процесса формоизменения можно было бы считать эти участки прямолинейными. Процесс деформирования подобной безмоментной оболочечной модели, состоящей из указанных элементов с прямолинейными образующими, рассматриваем как пошаговый, при котором переход из известного состояния в момент времени в новое состояние (в момент времени ), осуществляется с малыми приращениями деформаций. На данном малом временном интервале (шаге нагружения) формулировку задачи для принятой дискретной модели оболочки выполняем в терминах узловых перемещений с учетом изменения конфигурации оболочки за время . При этом используем цилиндрическую систему координат (, , ). Решение сформулированной физически и геометрически нелинейной задачи для дискретной модели оболочки на шаге нагружения сводится посредством итерационной процедуры к решению последовательности линейных задач. Линеаризация исходной нелинейной системы уравнений на шаге нагружения в рамках такой процедуры осуществляется с использованием методов Ньютона и переменных параметров упругости. Итерационные уточнения выполняются до достижения заданной относительной точности () по перемещениям. Решение соответствующей системы линейных алгебраических уравнений проводится по методу Гаусса. Следует отметить, что в качестве параметра нагружения при реализации описанной вычислительной модели принималось перемещение полюса формуемой оболочки. Давление при этом являлось вычисляемым параметром. Подобный подход обеспечивает устойчивый счет при численном моделировании рассматриваемого процесса гидроформовки с заходом в упомянутую выше запредельную область. Перейдем теперь к результатам исследований с применением описанной конечноэлементной модели. Сначала остановимся на результатах, подтверждающих надежность данной модели при решении рассматриваемого класса задач. Будем при этом использовать экспериментальные и расчетные данные работ [7, 8]. В работе [7] приведены данные по гидровыпучиванию круглых заготовок из листовой меди толщиной (0,927 мм) с радиусом (127 мм). Коэффициент нормальной анизотропии указанного материала равен единице (изотропный материал), а его кривая упрочнения представима в виде: , (1) где: , , . Рисунок 2. Зависимость гидростатического давления от толщиной деформации в полюсе оболочки в случае , Расчеты по предлагаемой модели проводились с числом элементов , с шагом по параметру нагружения и при выборе . О хорошем согласовании результатов такого конечноэлементного расчета с экспериментом позволяют судить данные, представленные на рисунке 2. Здесь - толщинная логарифмическая деформация в полюсе оболочки. Вместе с тем, расчетная модель [7], построенная как и предлагаемая конечноэлементная модель на основе теории течения для материала с изотропным упрочнением, демонстрирует неустойчивый счет по мере приближения к предельному значению параметра . В работе [8] рассмотрено решение задачи о гидровыпучивании круглой мембраны на основе теории течения и деформационной теории. Показано, что расчетная модель, построенная с использованием соотношений теории течения, лучше согласуется с экспериментом, особенно в ситуациях, близких к предельным. Приведено сравнение результатов расчета с экспериментальными данными, полученными Брамлеем и Меллором в 1965 году. В указанных экспериментах использовались заготовки из листовой стали толщиной (1,016 мм) с радиусом (127 мм). Коэффициент нормальной анизотропии данного материала определен в виде , а его кривая упрочнения в виде: , (2) где: , . Рисунок 3. Зависимость безразмерного параметра нагрузки от относительного перемещения полюса оболочки Рисунок 4. Зависимости безразмерного параметра нагрузки от относительного перемещения полюса для разных значений показателя кривой упрочнения На рисунке 3 представлено сравнение экспериментальных и расчетных данных, приведенных в работе [8], с результатами расчета с использованием разработанной конечноэлементной модели. Здесь безразмерный параметр нагрузки определяется формулой: , (3) где: - исходная толщина листового металла. Как видно, результаты расчета на основе предлагаемой конечноэлементной модели лучше согласуются с экспериментальными данными, чем аналогичные результаты, полученные на основе расчетной модели [8]. Отмечаем также, что разработанная модель демонстрирует устойчивый счет на любой стадии процесса формоизменения, в том числе и в запредельной области. Перейдем теперь к изложению результатов параметрических исследований, выполненных с применением предлагаемой конечноэлементной модели. Листовой материал при проведении таких исследований предполагался изотропным (). Показатель n кривой упрочнения, задаваемой в виде (2), варьировался в пределах от 0,15 до 0,5 с шагом 0,05. На рисунке 4, 5, 6 приведены зависимости безразмерных параметров , , (соответственно) от относительного перемещения полюса формуемой оболочки для упомянутых (восьми) значений показателя кривой упрочнения листового металла. Здесь - безразмерный параметр нагрузки, определяемый формулой (3); и - значения площади поверхности оболочки в ее исходном (плоском) и текущем (деформированном) состояниях; - толщинная логарифмическая деформация в полюсе оболочки. Рисунок 5. Зависимости относительного увеличения площади поверхности от относительного перемещения полюса для разных значений показателя кривой упрочнения Рисунок 6. Зависимости толщиной деформации в полюсе от относительного перемещения полюса для разных значений показателя кривой упрочнения Здесь следует учесть то обстоятельство, что в рассматриваемом процессе гидроформовки наиболее интенсивно деформирующимся участком оболочки является окрестность ее полюса (где имеет место наибольшее ее утонение). В этой зоне реализуется состояние однородного двухосного равностороннего растяжения, при котором и . Разрыв формуемой оболочки (происходящий в случае высопластичных листовых металов вследствие шейкообразования - локализации деформации) следует ожидать именно в этой зоне. В качестве критического (предельно допустимого - опасного на разрыв) значения толщинной деформации в обсуждаемом состоянии однородного двухосного равностороннего растяжения данного листового металла можно принять значение , непосредственно предшествующее моменту локализации деформации. Укажем теперь на то, каким образом с использованием диаграмм, представленных на рисунках 4, 5, 6, можно проводить оценки предельно допустимых для рассматриваемого процесса гидроформовки значений параметров , , применительно к заданному металлическому листу (с заданными характеристиками: , , , ). Полагаем сначала, что известно опасное на разрыв значение толщинной деформации. Тогда с использованием графиков на рисунке 6 определяем соответствующее значение . Это позволяет с использованием графиков на рисунках 4 и 5 определить искомые значения (а затем ) и . Если в соответствующем эксперименте непосредственно перед моментом разрыва формуемой оболочки зафиксировано значение перемещения ее полюса, то по графикам на рисунках 4, 5, 6 можно сразу определить искомые значения (а затем ), и . Если, наконец, разрыв оболочки в эксперименте фиксируется при давлении на стадии возрастания зависимости , то, зная значение , по графикам на рисунке 4 определяем соответствующее значение . Затем с использованием графиков на рисунках 5, 6 определяем искомые значения и . Выводы 1. С использованием осесимметричной жесткопластической безмоментной конечноэлементной модели выполнено исследование предельного поведения оболочки, формуемой из закрепленного по контуру круглого металлического листа под действием гидростатического давления. Достоверность получаемых результатов подтверждена сравнением с известными экспериментальными данными. 2. Построены диаграммы, позволяющие определять предельные параметры процессов гидроформовки круглых листовых заготовок с различными геометрическими и физико-механическими характеристиками.