CC BY-NC
27Программные системы и вычислительные методы
Правильная ссылка на статью:
Булычев Р.Н. — Метод построения многомерной модели для поиска значений параметров процесса формовки // Программные системы и вычислительные методы. - 2019. - № 3. DOI: 10.7256/2454-0714.2019.3.30474 URL: https//nbpublish.com'library_read_article.php?id=30474
Метод построения многомерной модели для поиска значений параметров процесса формовки
Булычев Роман Николаевич
аспирант, кафедра Инженерная и Компьютерная Графика, Восточно-Сибирский Государственный
Университет Технологий и Управления
670013, Россия, республика Бурятия, г. Улан-Удэ, ул. Ключевская, 40в, оф строение 1
И Buly4ev.roman@yandex.ru
Статья из рубрики "Математическое моделирование и вычислительный эксперимент"
Аннотация.
Объектом исследования являются геометрические модели процесса инкрементальной формовки деталей из листового материала. Предметом исследования является графическая модель для поиска оптимальных значений параметров процесса формовки на базе многомерной начертательной геометрии. Автором в статье рассмотрены основные оптимизирующие факторы и параметры процесса. Особое внимание уделяется задачам построения геометрических моделей определения оптимизирующих факторов (зона контакта, время обработки и качество поверхности получаемой детали) и оптимизации параметров (профиль и траектория движения формообразующего инструмента) для инкрементального формообразования. Методом исследования является способ построении графической оптимизационной модели процесса с использованием проекционного чертеж Радищева для многомерного пространства. Так же было применено математическое моделирование для проверки правильности полученных оптимальных параметров. Новизна исследования заключается в применении для многофакторной, многокомпонентной системы методов многомерной начертательной геометрии в исследовании процесса инкрементальной формовки с целью повышения качества изготавливаемых деталей. Основные выводы. С помощью графической оптимизационной модели на базе многомерной начертательной геометрии был получен диапазон значений параметров инкрементальной формовки для изделия прогнозируемого качества. Проведено моделирование процесса с полученными параметрами. Проведена инкрементальная формовка детали из листового металла конической формы. Полученные результаты соответствует заданным критериям качества.
Ключевые слова: инкрементальная формовка, послойное деформирование, листовая штамповка, чертеж Радищева, оптимизация, оптимизация параметров, оптимизирующие факторы, параметры процесса, многомерная геометрия, моделирование
DOI:
10.7256/2454-0714.2019.3.30474
Дата направления в редакцию:
08-08-2019
Дата рецензирования:
08-08-2019
Введение
Инкрементальная листовая формовка — это новый процесс обработки материала давлением. Основным преимуществом инкрементальной формовки, в сравнении с традиционными процессами листовой штамповки, являются снижение затрат и времени на подготовку производства детали. Это преимущество делает данную технологию привлекательной для единичного и мелкосерийного производства
Несмотря на преимущество, инкрементальная формовка имеет ряд недостатков, которые сдерживают применение данного процесса.
Первый недостаток процесса накладывает на форму конечной детали и называется —
«закон синуса» согласно которому невозможно получить вертикальные стенки, так как толщина готовой детали будет стремиться к нулю. Для данного процесса максимально допустимый угол стенки, который можно получить находиться в пределах 60 градусов. Согласно закону синуса при таком значении угла, толщина стенки уменьшается наполовину. Предельная величина зависит от пластичных свойств листового материала.
Второй недостаток процесса инкрементальной формовки связан с формообразующим инструментом. Рабочий профиль может формообразующего инструмента может иметь
различную геометрию, но лучше всех зарекомендовала себя сферическая форма За счет локального приложения нагрузки формообразующим инструментом на заготовку,
остается след в виде канавки на обрабатываемой поверхности I18!. Чем больше пятно контакта, тем отчетливее видны дефекты в виде гребней на детали. Радиус профиля рабочей части сферического формообразующего инструмента необходимо учитывать, что бы обрабатываемая часть поверхности готовой детали соответствовала стандартам качества штампуемых изделий из листа.
Третий недостаток процесса инкрементального формообразования — это время необходимое для получения из листовой заготовки готовой детали. Время подготовки производства и создания формообразующей оснастки для инкрементального формообразования намного меньше традиционных процессов листовой штамповки. Но машинное время инкрементальной формовки значительно больше остальных процессов обработки металлов давлением. Машинное время процесса инкрементальной формовки
сопоставимо со временем необходимым при фрезеровании детали на станке с ЧПУ [6].
Все вышеперечисленные недостатки инкрементальной формовки некритичны при запуске процесса получения детали на станке ЧПУ с оптимальными параметрами.
Авторы [17, 20, 24] выделяют два основных параметра оказывающих наибольшее влияние на процесс — форма профиля рабочей части инструмента Rt и расстояние между витками траектории Ах движения инструмента.
В зарубежных и отечественных работах исследовался вопрос оптимизации основных параметров процесса формовки. Так в работе ^^^ были получены оптимальные значения параметров процесса с заданной толщиной готовой детали. В работах ИДЫ22! поиск оптимальных параметров проводился, чтобы получить поверхность детали с заданными предельными отклонениями от исходной геометрии. В работе [16! оптимизация проводилась по нагрузкам, приходящимся на приводы станка с ЧПУ, чтобы можно было использовать менее мощное оборудование. Так же в зарубежных исследованиях вопросом оптимизации параметров процесса инкрементальной формовки занимались
авторы работ ЦЯЫ!4!. в отечественных исследованиях авторами работ проводилась оптимизация параметров для процессов локального деформирования, схожих по принципу с инкрементальной формовкой, а именно для ротационной вытяжки.
Все рассмотренные работы объединяет то, что при такой оптимизации не учитывается многогранность факторов данного процесса, а значит, оптимизация по одному фактору не дает объективной картины процесса, т.к. оптимальные значения параметров для одного фактора могут не соответствовать оптимальным значениям для другого фактора. Чтобы учесть все моменты процесса инкрементальной формовки нужно оптимизировать параметры процесса сразу по нескольким факторам. Применение традиционного математического моделирования для оптимизации данного процесса осложняется в связи с необходимостью выявления закономерностей, учитывающее большое число факторов. Математическое моделирование многофакторных процессов связано с выполнением большого объема вычислений и отсутствием наглядного представления об объекте исследования.
Данную задачу, учитывающую множество независимых параметров и факторов, можно решить, с помощью многомерной начертательной геометрии, которая даст наглядную интерпретацию полученных результатов. При построении графической оптимизационной модели для решения многокомпонентной системы часто используют проекционный
чертеж Радищева Данный метод можно применить и для оптимизации процесса инкрементального формообразования. Решению этой задачи посвящена настоящая с та тья .
Описание метода построения многомерной модели
В многомерной начертательной геометрии исходные данные и полученные результаты представляются в виде графической модели в форме линейных и нелинейных подпространств многомерного пространства и отношений между ними.
Возможности методов геометрического моделирования на основе применения чертежа
Радищева с целью оптимизации рассмотрены в работах 12—3—В работах сформулированы теоретические основы и рассмотрены задачи построения оптимизационных моделей и алгоритмы их решения.
Рассмотрим алгоритм построения области пересечения гиперповерхности с гиперплоскостью уровня, который реализуется следующим образом: гиперповерхность описывает зависимость оптимизирующих факторов от компонентов системы, а гиперплоскость задает требуемое значение оптимизирующего фактора.
В данной работе при построении графической оптимизационной модели использовались два параметров процесса формовки (х 1 - радиус профиля рабочей части формообразующего инструмента Rt , х 2 - расстояние между витками траектории Ах движения инструмента), в зависимости от значений трех оптимизирующих факторов (у 1 -качество поверхности полученной инкрементальной формовкой детали, у 2 - машинное время процесса формовки, у 3 - усилие формования), которые геометрически представляют собой 3-поверхности, заданные 1-поверхностями уровня, описывающие зависимость между параметрами процесса формовки и оптимизирующими факторами.
Пересечения гиперповерхности с гиперплоскостью уровня для каждого фактора процесса даст три кривые, которые образуют область оптимальных значений параметров. Значения компонентов, лежащие в полученной области, будут решением задачи по нахождению оптимальных значений параметров многофакторного процесса инкрементальной формовки.
Постановка задачи
В общем виде математическую постановку задачи можно сформулировать следующим образом. Компонентами системы являются параметры процесса формовки Х={ х1,..., хп} и критерии качества У={ у1,..., ут} . Каждому набору параметров X соответствуют числовые значения критериев У . Каждому из У можно задать значения у1* и возможные величины отклонения А i от этих значений: | у'1- у/'*1< А i .
Для заданных значений критериев У*={ у1*,...г ут*} требуется определить общую область значений параметров в многомерном пространстве Еп+т .
В частности, найдем оптимальные значения параметров инкрементальной формовки, форма профиля рабочей части инструмента Rt и расстояние между витками траектории Ах движения инструмента (рис. 1а), методом построения многомерной модели. Ограничим диапазон поиска значения оптимальных параметров для № от 6 до 15мм. (х1 ) и для Ах от 0,1 до 10 мм. (х2 ). Процесс формовки будем рассматривать для детали формой конуса (рис. 1б). Угол наклона стенки конуса составляет 60 градусов. Материал
заготовки АМцМ толщиной 1 мм.
Рисунок 1 - Постановка задачи
а - параметры для оптимизации процесса, б - параметры детали
Для поиска оптимальных значений выберем три оптимизирующих фактора. Нужно чтобы параметры процесса были подобны таким образом, чтобы высота гребня от формообразующего инструмента была не более 0,2 мм, чтобы поверхность детали
удовлетворяла требованиям при этом время формовки должно превышать 15 минут. При формовании детали будет использоваться фрезерный станок с ЧПУ, которой сможет развить усилие деформирования не больше 3x10 кН.
Решение задачи оптимизации
а. оптимизация по высоте гребня от следа формообразующего инструмента.
После формообразования на поверхности детали от локального приложения нагрузки формообразующим инструментом остаются следы в виде гребней (рис. 2а).
Рисунок 2 - Высота гребня
а - гребни от инструмента, б - расчет высоты гребня t
Высота гребней t зависит от радиуса профиля формообразующего инструмента Rt и шага между витками траектории Ах (рис. 2б) и ее можно определить из выражения:
По полученным значениям t , для диапазона значений х 1 и х 2 , построимгиперповерхность у 1( х 1, х 2) и гиперплоскость оптимизирующего фактора у 1= у 1оптим . (рис. 3).
Рисунок 3 - оптимизация по высоте гребня от следа формообразующего инструмента
На пересечении гиперповерхности у 1( х 1, х 2) и гиперплоскости оптимизирующего фактора у 1= у 1оптим, получим кривую оптимальных значений х 1 и х 2 по высоте гребня.
б. оптимизация по времени инкрементальной формовки.
Инкрементальная формовка проходит с постоянной скоростью перемещение инструмента.
Оптимальная скорость перемещения инструмента для формовки заготовки из
алюминиевых сплавов равна v обр. =1000 мм/мин l^ü. Определяем время формовки по расстоянию, которое проходит формообразующий инструмент. Стратегией формовки конуса выбрана послойная обработка. Это значит, что траектория состоит из горизонтального и вертикального перемещения. Инструмент делает один полный оборот в горизонтальной плоскости, а затем внедряется в заготовку на величину Az . Для нахождения длины траектории составлена программа в Mathcad, в которой определяется количество витков траектории, а затем определяется горизонтальное и вертикальное перемещение в зависимости от значений x 1 и x 2 (рис. 4). По полученным данным определяем время формовки t обр. в зависимости от значений x1 и x 2 .
По полученным значениям t обр , для диапазона значений x 1 и x 2 , построимгиперповерхность y 2( x 1, x 2) и гиперплоскость оптимизирующего фактора y 2= y 2оптим . (рис. 5).
Рисунок 4 - Определение количества витков и длины траектории формовки в MathCad
Рисунок 5 - Оптимизация по времени формовки
На пересечении гиперповерхности у 2( х 1, х 2) и гиперплоскости оптимизирующего фактора у 2= у 2оптим, получим кривую оптимальных значений х 1 и х 2 по времени формовки.
в. оптимизация по деформирующему усилию .
Наибольшее усилие формования приходиться на ось Z оно на порядок выше, чем по
направлениям X и Y . Для расчета деформирующего усилия по направлению Z воспользуемся выражениями, приведенными в работе -ШЬ
Pz = -—Fz
2 я* , (1)
где,
S - толщина стенки;
as - напряжение текучести(о^ = 170 МПа для листа АМцМ 1мм.); F z- проекция площади контакта на ось Z (рис. 6)
Для определения геометрической площади контакта заготовки с формообразующем инструментом и проекции на ось Z (для различных значений Rt и Az) используем моделирование с помощью CAD пакета Siemens NX, методика определения площади
пятна контакта описана статье -Ш1.
Рисунок 6 - Определение зоны контакта
На основе полученных проекций площадей для диапазона значений х 1 и х 2 и формулы (1) получим значениям деформирующего усилия по оси 2 . Построимгиперповерхность у 3( х 1, х 2) и гиперплоскость оптимизирующего фактора у 3= у Зоптим . (рис. 7).
Рисунок 7 - оптимизация по деформирующему усилию
На пересечении гиперповерхности у 3( х 1, х 2) и гиперплоскости оптимизирующего фактора у 3= у Зоптим, получим кривую оптимальных значений х 1 и х 2 по времени формовки.
г. результат оптимизации
При пересечении кривых с оптимальными значениями х1,х2 для факторов у1 , у2 , у3 получим зону, в которой выполняются заданные условия у1 = у1 оптим , у2 = у2 оптим и у3 = у3 оптим (рис. 8).
Рисунок 8 - Результат оптимизации
Так для заданной детали в виде конуса оптимальными параметрами формовки получились радиус инструмента Rt =10 мм. и шаг между витками Az = 3 мм. С этими параметрами данный процесс укладывается в ограничения заданными оптимальными факторами ограничения.
Для подтверждения оптимальных параметров проведем математическое моделирование процесса с применением ПО PamStamp.
Математическое моделирование инкрементальной формовки
С помощью CAD пакета Siemens NX смоделирован конус, заготовка, формообразующий инструмент и опора (рис. 9а). С помощью CAM пакета Siemens NX получена управляющая программа (рис. 9б), которая понадобилась для математического моделирования в среде Pam-Stamp. Методика моделирования процесса инкрементальной формовки описана в
статье -Ш.
Рисунок 9 - Результат оптимизации
а - управляющая программа, б - расчетная модель
Результаты математического моделирования процесса представлены на рисунке 10. Отклонения полученной детали от первоначально заданной панели (рис.9.в), при данном способе формовки получились в пределах ±0,2мм (рис.10.а). На рисунке 10б показано распределение толщины детали, которое соответствует закону синуса. На рисунке 10в показано распределение деформаций готовой детали. Усилия на приводы станка с ЧПУ сопоставимы с расчетными значениями, порядка 3кН, и находятся в пределах рабочего диапазона оборудования.
Рисунок 10 - Результат математического моделирования в PamStamp а - отклонения, б - распределение толщины, в - деформации
Для ее изготовления была спроектирована и изготовлена оснастка. Инкрементальная формовка конуса производилась из заготовки материала АМцМ лист толщиной 1мм, а результаты представлены на рисунке 10. Результат формовки сходиться с результатом математического моделирования. Процесс проходил с полученными оптимальными параметрами.
Рисунок 11 - Результат инкрементальной формовки конуса
Заключение
Применение графической многомерной модели и разработанных на ее основе алгоритмов решения задачи позволяет получать значения параметров процесса инкрементальной формовки детали из листового материала, при этом наглядно оценивать процесс формовки, оперативно устанавливать оптимальные режимы, параметры, а также прогнозировать характеристики исследуемого процесса.
Исследование процесса инкрементальной формовки как многофакторной, многокомпонентной системы методами многомерной начертательной геометрии помогут повысить качество и расширить номенклатуры выпускаемых изделий.
Библиография
1. Аюшеев Т.В. Моделирование процесса изготовления деталей методом инкрементальной формовки [Текст] / Т.В. Аюшеев, Р.Н. Булычев // Проблемы механики современных машин. Материалы VII международной научной конференции. Том
2. Изд-во ВСГУТУ, 2018.-C. 104-108 2. Волков, В.Я. Многомерная исчислительная геометрия: монография [Текст] / В.Я. Волков, В.Ю. Юрков.-Омск: Изд-во ОмГПУ, 2008.-244 с
3. Волков, В.Я. Графические оптимизационные модели многофакторных процессов: монография [Текст] / В.Я. Волков, М.А. Чижик.-Омск: Изд-во ОГИС, 2009.-101 с
4. ГОСТ 21631-76. Листы из алюминия и алюминиевых сплавов. Технические условия [Текст].-Москва: Изд-во стандартов, 1977.-29 с
5. Гусев А.И. Оптимизация параметров процесса ротационной вытяжки [Текст] / А.И. Гусев // Исследования в области пластичности и обработки металлов давлением, Тульский Политехнический Институт, 1981,-Вып.8,-С. 101-105
6. Кривошеин В. А. Перспективы использования технологий инкрементальной формовки в современном производстве [Текст] / В.А. Кривошеин, А.А. Анцифиров, Ю.В. Майстров // Известия высших учебных заведений. Машиностроение, 2014.-Вып. 11 (656).-С. 84-89
7. Новик Ф.С. Оптимизация процессов технологии металлов методами планирования экспериментов [Текст] / Ф.С. Новик, Я.Б. Арсов // М.: машиностроение.-1980.-Т.
304.-С. 42-45
8. ОСТ 1 02507-92. Самолеты дозвуковые. Общие требования к качеству внешней поверхности [Текст].-Взамен ОСТ 1 02507-84; введ. 1993-01-01-Москва: Изд-во стандартов, 1993.-10 с
9. Радищев В.П. О применении геометрии четырех измерений к построению равновесных физико-химических диаграмм [Текст] / В.П. Радищев // Изв. СФХА.-М., 1947.-Т. 15-С. 129-134
10. Чижик М.А. Графические оптимизационные модели многопараметрических технологических процессов легкой промышленности [Текст] / М.А. Чижик, В.Я. Волков // Инженерный вестник Дона, 2012.-Т. 21.-№. 3.-С. 87-94
11. Чумадин А.С. Исследование процессов послойного деформирования листовых заготовок с использование фрезерного станка с ЧПУ [Текст] / А.С. Чумадин, Д.А. Батурин // Кузнечно-штамповочное производство. Обработка материалов давлением, 2014.-Вып.7.-С. 29-32
12. Abass K. A study to comparing spherical, ellipse and flat forming tool profile effect in single point incremental forming by finite element analysis [Текст] // Bull., Series D. -2016.-Т. 78.-Вып. 1.-С. 173-184
13. Adams D. A new model for contact geometry in single-point incremental forming [Текст] / D. Adams, J. Jeswiet // Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture, 2015.-Т. 229. №. 6.-С. 982-989.-DOI: 10.1177/0954405414535580
14. Attanasio A. Optimization of tool path in two points incremental forming [Текст] / A. Attanasio, E. Ceretti, C. Giardini // Journal of Materials Processing Technology, 2006.Т. 177.-Вып. 1-3.-С. 409-412.-D0I: 10.1016/j.jmatprotec.2006.04.047
15. Azaouzi M. Tool path optimization for single point incremental sheet forming using response surface method [Текст] / M. Azaouzi, N. Lebaal // Simulation Modeling Practice and Theory, 2012.-T. 24.-C. 49-58.-DOI: 10.1016/j.simpat.2012.01.008
16. Bahloul R. Application of response surface analysis and genetic algorithm for the optimization of single point incremental forming process [Текст] / R. Bahloul, H. Arfa, H. Belhadjsalah // Key Engineering Materials, 2013.-T. 554-557.-С. 1265-1272.-DOI: 10.4028/www.scientific.net/KEM. 554-557.1265
17. Gatea S. Review on the influence of process parameters in incremental sheet forming [Текст] / S. Gatea, H. Ou, G. McCartney // The International Journal of Advanced Manufacturing Technology, 2016.-Т. 87.-Вып. 1-4.-С. 479-499. - DOI: 10.1007/s00170-016-8426-6
18. Guangcheng Z. Forming process of automotive body panel based on incremental forming technology [Текст] / Z. Guangcheng, X. Jinbo, S. Xiaofan, Z. Xun, L. Chuankai // Metallurgical and Mining Industry, 2015.-T. 1. - Вып. 12.-С. 350-357
19. Gulati V. Process parameters optimization in single point incremental forming [Текст] / V. Gulati [и др.] // Journal of The Institution of Engineers, 2016.-Т. 97.-Вып. 2.-С. 185-193.-DOI: 10.1007/s40032-015-0203-z
20. Hagan E. A review of conventional and modern single-point sheet metal forming methods [Текст] / E. Hagan, J. Jeswiet // Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture, 2003.-Т. 217.-Вып. 2.-С. 213-225.-DOI: 10.1243/095440503321148858
21. Jeswiet J. Rapid prototyping of a headlight with sheet metal [Текст] / J. Jeswiet, E. Hagan // Canadian Institute of Mining, Metallurgy and Petroleum, 2001.-С. 109-114
22. Li Y. Investigation and optimization of deformation energy and geometric accuracy in
the incremental sheet forming process using response surface methodology [Текст] / Y. Li [и др.]. The International Journal of Advanced Manufacturing Technology, 2015.-Т. 79.-Вып. 9-12.-С. 2041-2055.-D01 10.1007/s00170-015-6986-5
23. Moayedfar M. Parameter optimization and stretch enhancement of AISI 316 sheet using rapid prototyping technique [Текст] / Moayedfar M. [и др.] // IOP Conference Series: Materials Science and Engineering. - IOP Publishing, 2017.-Т. 257.-Вып. 1.-С. 1-13.-DOI: 10.1088/1757-899X/257/1/012006
24. Nimbalkar D.H. Review of incremental forming of sheet metal components [Текст] / D.H. Nimbalkar, V.M. Nandedkar // International Journal of Engineering Research and Applications.-2013.-vol. 3,-issue 5.-С. 39-51.
25. Young D. Wall thickness variations in single-point incremental forming [Текст] / D. Young, J. Jeswiet // Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture, 2004.-Т. 218.-№. 11.-С. 1453-1459.-DOI: 10.1243/0954405042418400
