УДК 621.98
Численное моделирование процесса калибровки менисковых деталей жидкой технологической средой
А.О. Мазуркевич, А.В. Петюков, В.Д. Баскаков
Методом конечно-элементного анализа проведено исследование напряженно-деформированного состояния детали менисковой формы на стадии нагрузки жидкой технологической средой и последующей разгрузки. Получена зависимость пружинения от геометрических параметров заготовки и внешнего прикладываемого давления. Установлено доминирующее влияние радиуса кривизны мениска и его высоты на величину пружинения.
Ключевые слова: численное моделирование, ANSYS Mechanical APDL, калибровка менисковых деталей, регрессионный анализ.
Numerical simulation of meniscus shape parts calibration by liquid process medium
A.O. Mazurkevich, A.V. Petyukov, V.D. Baskakov
The research of the stress-strain state of the meniscus shape part has been carried out at the stage of the liquid process medium loading and subsequent discharge by the finite element analysis. The relation of the springing and the workpiece geometric parameters and the external pressure applied is obtained. The dominant influence of the meniscus curvature radius and of its height on the springing value is established.
Keywords: numerical simulation, ANSYS Mechanical APDL, meniscus parts calibration, regression analysis.
/'Лдной из важнейших деталей, определяющих эффективность действия снарядоформирующих зарядов, является кумулятивная облицовка менисковой формы, изготовляемая из стали или меди. Она позволяет сформировать при взрыве удлиненный поражающий элемент, пробивающий танковую броню с больших расстояний [1].
В деталях подобного класса предъявляют высокие требования к геометрической точности. В частности, отклонение от сферичности не должно превышать 0,03...0,05 мм, а разностенность* — 0,03...0,05 мм.
Одна из наиболее перспективных технологий получения облицовок включает операцию калибровки внутренней поверхности и позволяет обеспечить точность ее формы и в дальнейшем, при механической обработке, обеспечить снижение разностенности детали [2]. Основным способом проведения этой операции является деформирование сферической заготовки за счет внешнего гидростатического давления
МАЗУРКЕВИЧ Алексей Олегович
студент
ПЕТЮКОВ Андрей Вячеславович
студент кафедры «Высокоточные летательные аппараты»
БАСКАКОВ Владимир Дмитриевич
кандидат технических наук, доцент кафедры
«Технологии ракетно-космического машиностроения» (МГТУ им. Н.Э. Баумана)
e-mail: lex19900310@yandex.ru.ru
* Разностенность — разность между наибольшей и наименьшей толщиной детали в плоскости поперечного сечения.
ВШЗЕХЭШШ] выкшшех ©аведжшй
Рис. 1. Схема проведения калибровки:
1 — сферическая сегментная заготовка;
2 — сферическая матрица; Н — высота сегментной заготовки; S — толщина заготовки; R3 — радиус заготовки; R — радиус матрицы;
Р — прикладываемое давление
P (давление жидкой технологической среды, либо полимера) до полного соприкосновения внутренней поверхности заготовки, имеющей в начале операции радиус R, с поверхностью матрицы радиусом Лм, при этом Лз > Лм (рис. 1).
В данной работе проводилось численное моделирование калибровки на обжим сферической заготовки при помощи программного комплекса ANSYS Mechanical APDL с целью выявления основных закономерностей ее поведения.
C помощью программного комплекса ANSYS построена геометрическая модель процесса калибровки в лагранжевой постановке в трехмерной геометрии. В силу симметрии задачи рассматривалась 1/4 часть модели.
На рисунке 2 представлена конечно-элементная модель системы тел, участвующих в процессе калибровки. В качестве материала оправки использовалась модель абсолютно жесткого тела. Для материала заготовки была выбрана модель упругопластического материала с кинематическим упрочнением. Параметры материалов приведены в табл. 1. Заготовка разбита 20-узловыми объемными конечными элементами. Число элементов, на которые разбита заготовка, составило 4 800, оправка — 1 850.
Таблица 1
Механические свойства материала заготовки
Свойство Медь М1 Сталь Ст3
Плотность р, кг/м3 8940 7850
Модуль упругости Е, ГПа 110 210
Коэффициент Пуассона V 0,35 0,28
Предел текучести стт, МПа 120 235
Предел прочности ств, МПа 220 410
Модуль упрочнения Е, МПа 166 764
Рис. 2. Конечно-элементная модель системы тел, участвующих в процессе калибровки:
1 — заготовка; 2 — матрица
В начальный и последующие моменты времени к заготовке прикладывалось давление, закон изменения которого приведен на рис. 3. В плоскости симметрии ХОУ запрещалось перемещение вдоль оси Z и поворот вокруг осей У и X всех узлов заготовки, лежащих в этой плоскости. Аналогично для плоскости YOZ запрещалось перемещение вдоль оси X и поворот вокруг осей У и Z узлов, лежащих в этой плоскости. Для задания контактных условий между матрицей и заготовкой использовался, так называемый штрафной метод (penalty formulation). При этом сила отталкивания соприкасающихся ячеек материалов прямо пропорциональна глубине их взаимопроникнове-
Р, M Па
О [ 2 3 4 (/С
Рис. 3. Изменение во времени t прикладываемого к заготовке давления Р: Рт — максимальный уровень давления
ния. Трение между матрицей и заготовкой не учитывалось.
В данной работе пружинение заготовки П в процессе /-го моделирования калибровки определялось как модуль вектора перемещения крайней точки грани, принадлежащей внутренней поверхности заготовки, при разгрузке (максимальное перемещение заготовки вследствие упругого пружинения).
В качестве величин, имеющих решающее влияние на результат калибровки, были выбраны следующие (табл. 2): начальный радиус заготовки Лз; начальная толщина облицовки 5; начальная высота облицовки Н; максимальное давление нагружения на этапе калибровки Рт.
Предполагая, что зависимость между факторами и величиной пружинения является линейной, для построения искомого полинома был использован метод множественного регрессионного анализа [3]. В соответствии с выбранным подходом численное моделирование проводилось по плану полного факторного эксперимента, в котором каждый фактор принимал два значения. Для построения одного полинома проводилось п = 24 =16 расчетов.
Таблица 2
Диапазон изменения факторов, влияющих на пружинение
Уровни факторов Яз, мм 5, мм Н, мм Рт, МПа (Ст3) Р, т МПа (М1)
Основной уровень 62,4 2,86 11,7 216,25 110
Интервал варьирования 5,2 0,78 3,9 136,25 70
Верхний уровень 67,6 3,64 15,6 352,5 180
Нижний уровень 57,2 2,08 7,8 80 40
В соответствии с матрицей планирования численных расчетов был проведен ряд вычислений, в результате которых получены величины пружинения для заготовок, выполненных из меди М1 и стали Ст3.
Рассмотрим подробнее результаты моделирования одной из расчетных задач. Подтверждением корректности использования симметрии задачи с целью минимизации расчетов, а также правильности выбора граничных усло-
вий, служит осесимметричное распределение модуля вектора перемещения N узлов заготовки (рис. 4). В целом аналогичную картину демонстрирует и распределение интенсивностей напряжений. Неравномерность распределения величины интенсивности напряжений по толщине говорит о правильности выбора типа конечного элемента для заготовки (рис. 5).
Рис. 4. Распределение модуля вектора перемещения N узлов заготовки в конечный момент времени
Рис. 5. Распределение интенсивности напряжений ст; в узлах заготовки в конечный момент времени
На основе полученных результатов попытаемся установить зависимость между начальными параметрами задач (факторами) и величинами упругого пружинения заготовки Пг, при этом предположим, что искомая связь описывается полиномом вида
П = а 0 + ,
(1)
} =1
где П — оцененное значение пружинения, мм; а0 — постоянная составляющая пружинения, мм; а! — коэффициент влияния у-го фак-
ВШзехэшш] выкшшпх ©аведжшй
тора на величину пружинения; х. — значение у-го фактора; / = 4 — число факторов.
Определение коэффициентов линейного полинома (1) для заготовок, выполненных из указанных материалов, было проведено с использованием линейного множественного регрессионного анализа (табл. 3).
Таблица 3
Значения коэффициентов уравнения регрессии
Коэффициент влияния Размерность Ст3 М1
а0 (постоянная составляющая) мм 0,157 0,161
а1 (радиус заготовки Яз) — -0,456 -0,928
а2 (толщина заготовки 5) — -0,019 -0,014
а3 (высота заготовки Н) — 5,518 5,175
а4 (максимальное давление Ря) м/МПа -2,5-10-8 —1,2-10-7
Близость результатов наблюдений и оценки с использованием полинома (1) была установлена путем вычисления коэффициента множественной корреляции Я. Величина Я для заготовок, выполненных из сплава М1, составила 0,95, а из стали Ст3 — 0,93. Значимость коэффициента множественной корреляции оценивали с помощью критерия Фишера. Связь считают статистически значимой в том случае, если расчетное значение РЯ (коэффициент Фишера) превышает табличное ¥т. Расчетное значение Рк для
стали составляет 19,10, для меди — 26,25; ¥т = = 14,45 для уровня значимости 0,01 и числа степеней свободы V! = 11, v2 = 4 [3].
Аналогично по ^критерию Стьюдента (для уровня значимости 0,01): расчетное значение ?Я для Ст3 составляет 24,64, для М1 — 33,28Г равно 3,11 для уровня значимости 0,01 и числа степеней свободы V! = 11 [3].
Близкое к единице значение коэффициента множественной корреляции Я подтверждает гипотезу о линейном характере зависимости величины пружинения от выбранных факторов в заданном диапазоне варьирования факторов.
Запишем полином в безразмерном виде:
П = а 0 + а1Я3 + а 2 5 + а 3 Н + а4 Ря
(2)
где П =П/Ям; а0 = а0 /Ям; 5 = 5/Ям; Н = Н / Ям; аА = а4 а Т / Ям; Ря = Рт / а т.
Полином в форме (2) работает в следующем диапазоне значений факторов:
Я3 = [1,1...1,3]; 5 = [0,04...0,07]; Н = [0,15...0,30]; Ря =[Рт1П;1,5].
Здесь Рт1п — безразмерное давление, минимально необходимое для обеспечения полного прилегания заготовки к матрице.
Полученная полиномиальная зависимость позволяет провести оценку пружинения сферической заготовки (см. рис. 1), величины прикладываемого внешнего воздействия, а также размеров, которые лежат в диапазоне варьирования выбранных факторов калибровки.
Характер воздействия факторов на величину пружинения определяется знаком коэффициентов влияния а. Отрицательная величина коэффициентов влияния радиуса заготовки, толщины заготовки, а также прикладываемого давления говорит о том, что данные факторы стремятся снизить пружинение. Сходность знаков коэффициентов полиномов, полученных для заготовок, выполненных из стали Ст3 и сплава М1, подтверждает индифферентность характера их влияния к материалу заготовки.
При помощи полученных зависимостей были установлены абсолютные значения вкладов, которые вносят факторы в переменную часть полинома пружинения (рис. 6).
Диаграмма демонстрирует качественное совпадение распределения вкладов факторов для медной и стальной заготовок.
На основе приведенной диаграммы вкладов, а также знаков коэффициентов влияния, можно сделать вывод о том, что для уменьшения величины пружинения необходимо калибровать заготовки с большим припуском на механообработку (с большей толщиной). Повышение максимального давления калибровки приводит к незначительному уменьшению пружинения. Исходя из этого, значение данного параметра рекомендуется принимать минимально необходимым, при котором обеспечи-
Рис. 6. Диаграмма вкладов факторов в переменную составляющую пружинения для различных материалов заготовки
вается закрытие зазора между матрицей и заготовкой, что согласуется с общепринятой практикой и позволяет повысить срок службы штамповой оснастки.
С помощью полинома (2) можно решать как прямые задачи оценки пружинения при заданных значениях конструктивных параметров заготовки, так и обратные — по допустимой величине пружинения определить параметры заготовки. Однако, в условиях реального производства факторы R3, 5, H, Pm являются случайными значениями. Поэтому также представляет интерес оценка диапазона разброса случайной составляющей пружинения ЛП:
Л П = ax ЛR3 + a2Л5 + a3 ЛH + a4 ЛPm,
где ЛП = ЛП/RM; AR3 = ЛRз /RM; Л5 = Л5/_Я„; ЛИ = ЛИ /Rм; ЛPm = ЛPm /ат ; ЛRз, Л5, ЛН, ЛРт — нестабильность радиуса заготовки, ее толщины, ее высоты и прикладываемого давления соответственно.
Таким образом, для оценки рассеяния величины пружинения необходимо определить диапазоны колебания факторов. Эта информация может быть получена по результатам измерений партии заготовок перед калибровкой, а также из паспортных данных прессового оборудования.
Выводы
1. Предложен подход к построению аналитических соотношений для оценки номинальной и случайной составляющих величины пружинения однотипных заготовок после калибровки жидкими или эластичными средами путем обработки результатов численного эксперимента и построения на этой основе регрессионной зависимости.
2. На основе изложенного подхода получены значения коэффициентов уравнения регрессии для оценки пружинения, а также нестабильности пружинения, заготовок типа сферический сегмент, выполненных из Ст3 и М1.
3. Анализ полученных соотношений показал, что для заготовки типа сферический сегмент наибольшее влияние на величину пружи-нения оказывают ее геометрические параметры: радиус R3 и высота H.
Литература
1. Средства поражения и боеприпасы / А.В. Бабкин, В.А. Велданов, Е.Ф. Грязнов и др.; под общ. ред. В.В. Селиванова. М.: Изд-во МГТУ им. Н.Э. Баумана, 2008. 984 с.
2. Тарасов В.А., Баскаков В.Д., Круглое П.В. Методика проектирования технологии изготовления высокоточных деталей боеприпасов // Оборонная техника. 2000. № 1-2. С. 89—92.
3. Львовский Е.Н. Статистические методы построения эмпирических формул. М.: Высш. шк., 1988. 239 с.
Статья поступила в редакцию 02.05.2012