Translate this page:
Please select your language to translate the article


You can just close the window to don't translate
Library
Your profile

Back to contents

Software systems and computational methods
Reference:

A method for constructing a multidimensional model for finding values of molding process parameters

Bulychev Roman Nikolaevich

Postgraduate Student, Department of Engineering and Computer Graphics, East Siberian State University of Technology and Management

670013, Russia, respublika Buryatiya, g. Ulan-Ude, ul. Klyuchevskaya, 40v, of. stroenie 1

Buly4ev.roman@yandex.ru

DOI:

10.7256/2454-0714.2019.3.30474

Received:

07-08-2019


Published:

08-09-2019


Abstract: The object of study is the geometric model of the incremental molding process of parts from sheet material.The subject of the study is a graphical model for finding the optimal values of the parameters of the molding process based on multidimensional descriptive geometry. The author in the article examined the main optimizing factors and process parameters. Particular attention is paid to the problems of constructing geometric models for determining optimizing factors (contact zone, processing time and surface quality of the resulting part) and parameter optimization (profile and trajectory of the shaping tool) for incremental shaping.The research method is a method of constructing a graphical optimization model of the process using the projection drawing of Radishchev for multidimensional space. Mathematical modeling was also applied to verify the correctness of the obtained optimal parameters. The novelty of the study lies in the application of multidimensional descriptive geometry methods for a multifactor, multicomponent system in the study of the incremental molding process in order to improve the quality of manufactured parts. Main conclusions. Using a graphical optimization model based on multidimensional descriptive geometry, a range of values of incremental molding parameters was obtained for a product of predicted quality. The process was simulated with the obtained parameters. An incremental molding of a conical-shaped sheet metal part was carried out. The results obtained correspond to the specified quality criteria.


Keywords:

incremental sheet forming, layer-by-layer deformation, sheet stamping, Radishchev's drawing, optimization, parameter optimization, optimizing factors, process parameters, multidimensional geometry, modeling


Введение

Инкрементальная листовая формовка — это новый процесс обработки материала давлением. Основным преимуществом инкрементальной формовки, в сравнении с традиционными процессами листовой штамповки, являются снижение затрат и времени на подготовку производства детали. Это преимущество делает данную технологию привлекательной для единичного и мелкосерийного производства [6].

Несмотря на преимущество, инкрементальная формовка имеет ряд недостатков, которые сдерживают применение данного процесса.

Первый недостаток процесса накладывает на форму конечной детали и называется — «закон синуса» [25], согласно которому невозможно получить вертикальные стенки, так как толщина готовой детали будет стремиться к нулю. Для данного процесса максимально допустимый угол стенки, который можно получить находиться в пределах 60 градусов. Согласно закону синуса при таком значении угла, толщина стенки уменьшается наполовину. Предельная величина зависит от пластичных свойств листового материала.

Второй недостаток процесса инкрементальной формовки связан с формообразующим инструментом. Рабочий профиль может формообразующего инструмента может иметь различную геометрию, но лучше всех зарекомендовала себя сферическая форма [12]. За счет локального приложения нагрузки формообразующим инструментом на заготовку, остается след в виде канавки на обрабатываемой поверхности [18]. Чем больше пятно контакта, тем отчетливее видны дефекты в виде гребней на детали. Радиус профиля рабочей части сферического формообразующего инструмента необходимо учитывать, что бы обрабатываемая часть поверхности готовой детали соответствовала стандартам качества штампуемых изделий из листа.

Третий недостаток процесса инкрементального формообразования — это время необходимое для получения из листовой заготовки готовой детали. Время подготовки производства и создания формообразующей оснастки для инкрементального формообразования намного меньше традиционных процессов листовой штамповки. Но машинное время инкрементальной формовки значительно больше остальных процессов обработки металлов давлением. Машинное время процесса инкрементальной формовки сопоставимо со временем необходимым при фрезеровании детали на станке с ЧПУ [6].

Все вышеперечисленные недостатки инкрементальной формовки некритичны при запуске процесса получения детали на станке ЧПУ с оптимальными параметрами.

Авторы [17, 20, 24] выделяют два основных параметра оказывающих наибольшее влияние на процесс — форма профиля рабочей части инструмента Rt и расстояние между витками траектории Δz движения инструмента.

В зарубежных и отечественных работах исследовался вопрос оптимизации основных параметров процесса формовки. Так в работе [23] были получены оптимальные значения параметров процесса с заданной толщиной готовой детали. В работах [15],[22] поиск оптимальных параметров проводился, чтобы получить поверхность детали с заданными предельными отклонениями от исходной геометрии. В работе [16] оптимизация проводилась по нагрузкам, приходящимся на приводы станка с ЧПУ, чтобы можно было использовать менее мощное оборудование. Так же в зарубежных исследованиях вопросом оптимизации параметров процесса инкрементальной формовки занимались авторы работ [19],[14]. В отечественных исследованиях авторами работ [5],[7] проводилась оптимизация параметров для процессов локального деформирования, схожих по принципу с инкрементальной формовкой, а именно для ротационной вытяжки.

Все рассмотренные работы объединяет то, что при такой оптимизации не учитывается многогранность факторов данного процесса, а значит, оптимизация по одному фактору не дает объективной картины процесса, т.к. оптимальные значения параметров для одного фактора могут не соответствовать оптимальным значениям для другого фактора. Чтобы учесть все моменты процесса инкрементальной формовки нужно оптимизировать параметры процесса сразу по нескольким факторам. Применение традиционного математического моделирования для оптимизации данного процесса осложняется в связи с необходимостью выявления закономерностей, учитывающее большое число факторов. Математическое моделирование многофакторных процессов связано с выполнением большого объема вычислений и отсутствием наглядного представления об объекте исследования.

Данную задачу, учитывающую множество независимых параметров и факторов, можно решить, с помощью многомерной начертательной геометрии, которая даст наглядную интерпретацию полученных результатов. При построении графической оптимизационной модели для решения многокомпонентной системы часто используют проекционный чертеж Радищева [9]. Данный метод можно применить и для оптимизации процесса инкрементального формообразования. Решению этой задачи посвящена настоящая статья.

Описание метода построения многомерной модели

В многомерной начертательной геометрии исходные данные и полученные результаты представляются в виде графической модели в форме линейных и нелинейных подпространств многомерного пространства и отношений между ними.

Возможности методов геометрического моделирования на основе применения чертежа Радищева с целью оптимизации рассмотрены в работах [2, 3, 10]. В работах сформулированы теоретические основы и рассмотрены задачи построения оптимизационных моделей и алгоритмы их решения.

Рассмотрим алгоритм построения области пересечения гиперповерхности с гиперплоскостью уровня, который реализуется следующим образом: гиперповерхность описывает зависимость оптимизирующих факторов от компонентов системы, а гиперплоскость задает требуемое значение оптимизирующего фактора.

В данной работе при построении графической оптимизационной модели использовались два параметров процесса формовки (x1 – радиус профиля рабочей части формообразующего инструмента Rt, x2 – расстояние между витками траектории Δz движения инструмента), в зависимости от значений трех оптимизирующих факторов (y1 - качество поверхности полученной инкрементальной формовкой детали, y2 - машинное время процесса формовки, y3 – усилие формования), которые геометрически представляют собой 3-поверхности, заданные 1-поверхностями уровня, описывающие зависимость между параметрами процесса формовки и оптимизирующими факторами.

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

Постановка задачи

В общем виде математическую постановку задачи можно сформулировать следующим образом. Компонентами системы являются параметры процесса формовки X={x1,…,xn} и критерии качества Y={y1,…,ym}. Каждому набору параметров X соответствуют числовые значения критериев Y. Каждому из Y можно задать значения y1* и возможные величины отклонения Δi от этих значений: |yi- yi*|≤ Δi.

Для заданных значений критериев Y*={y1*,…,ym*} требуется определить общую область значений параметров в многомерном пространстве En+m.

В частности, найдем оптимальные значения параметров инкрементальной формовки, форма профиля рабочей части инструмента Rt и расстояние между витками траектории Δz движения инструмента (рис. 1а), методом построения многомерной модели. Ограничим диапазон поиска значения оптимальных параметров для Rt от 6 до 15мм. (x1) и для Δz от 0,1 до 10 мм. (x2). Процесс формовки будем рассматривать для детали формой конуса (рис. 1б). Угол наклона стенки конуса составляет 60 градусов. Материал заготовки АМцМ [4], толщиной 1 мм.

Щелкните на изображении чтобы закрыть окно

Рисунок 1 – Постановка задачи

а - параметры для оптимизации процесса, б – параметры детали

Для поиска оптимальных значений выберем три оптимизирующих фактора. Нужно чтобы параметры процесса были подобны таким образом, чтобы высота гребня от формообразующего инструмента была не более 0,2 мм, чтобы поверхность детали удовлетворяла требованиям [8], при этом время формовки должно превышать 15 минут. При формовании детали будет использоваться фрезерный станок с ЧПУ, которой сможет развить усилие деформирования не больше 3x10 кН.

Решение задачи оптимизации

а. оптимизация по высоте гребня от следа формообразующего инструмента.

После формообразования на поверхности детали от локального приложения нагрузки формообразующим инструментом остаются следы в виде гребней (рис. 2а).

Щелкните на изображении чтобы закрыть окно

Рисунок 2 – Высота гребня

а – гребни от инструмента, б – расчет высоты гребня t

Высота гребней t зависит от радиуса профиля формообразующего инструмента Rt и шага между витками траектории Δz (рис. 2б) и ее можно определить из выражения:

По полученным значениям t, для диапазона значений x1 и x2, построимгиперповерхность y1(x1,x2) и гиперплоскость оптимизирующего фактора y1=y1оптим. (рис. 3).

Щелкните на изображении чтобы закрыть окно

Рисунок 3 – оптимизация по высоте гребня от следа формообразующего инструмента

На пересечении гиперповерхности y1(x1,x2) и гиперплоскости оптимизирующего фактора y1=y1оптим, получим кривую оптимальных значений x1 и x2 по высоте гребня.

б. оптимизация по времени инкрементальной формовки.

Инкрементальная формовка проходит с постоянной скоростью перемещение инструмента. Оптимальная скорость перемещения инструмента для формовки заготовки из алюминиевых сплавов равна vобр.=1000 мм/мин [21]. Определяем время формовки по расстоянию, которое проходит формообразующий инструмент. Стратегией формовки конуса выбрана послойная обработка. Это значит, что траектория состоит из горизонтального и вертикального перемещения. Инструмент делает один полный оборот в горизонтальной плоскости, а затем внедряется в заготовку на величину Δz. Для нахождения длины траектории составлена программа в Mathcad, в которой определяется количество витков траектории, а затем определяется горизонтальное и вертикальное перемещение в зависимости от значений x1 и x2 (рис. 4). По полученным данным определяем время формовки tобр. в зависимости от значений x1 и x2.

По полученным значениям tобр, для диапазона значений x1 и x2, построимгиперповерхность y2(x1,x2) и гиперплоскость оптимизирующего фактора y2=y2оптим. (рис. 5).

Щелкните на изображении чтобы закрыть окно

Рисунок 4 - Определение количества витков и длины траектории формовки в MathCad

Щелкните на изображении чтобы закрыть окно

Рисунок 5 – Оптимизация по времени формовки

На пересечении гиперповерхности y2(x1,x2) и гиперплоскости оптимизирующего фактора y2=y2оптим, получим кривую оптимальных значений x1 и x2 по времени формовки.

в. оптимизация по деформирующему усилию.

Наибольшее усилие формования приходиться на ось Z [21], оно на порядок выше, чем по направлениям X и Y. Для расчета деформирующего усилия по направлению Z воспользуемся выражениями, приведенными в работе [11]:

, (1)

где,

S – толщина стенки;

σs – напряжение текучести(σs = 170 МПа для листа АМцМ 1мм.);

Fz- проекция площади контакта на ось Z (рис. 6)

Для определения геометрической площади контакта заготовки с формообразующем инструментом и проекции на ось Z (для различных значений Rt и Δz) используем моделирование с помощью CAD пакета Siemens NX, методика определения площади пятна контакта описана статье [13].

Щелкните на изображении чтобы закрыть окно

Рисунок 6 – Определение зоны контакта

На основе полученных проекций площадей для диапазона значений x1 и x2 и формулы (1) получим значениям деформирующего усилия по оси Z. Построимгиперповерхность y3(x1,x2) и гиперплоскость оптимизирующего фактора y3=y3оптим. (рис. 7).

Щелкните на изображении чтобы закрыть окно

Рисунок 7 – оптимизация по деформирующему усилию

На пересечении гиперповерхности y3(x1,x2) и гиперплоскости оптимизирующего фактора y3=y3оптим, получим кривую оптимальных значений x1 и x2 по времени формовки.

г. результат оптимизации

При пересечении кривых с оптимальными значениями x1,x2 для факторов y1,y2,y3 получим зону, в которой выполняются заданные условия y1=y1оптим, y2=y2оптим и y3=y3оптим (рис. 8).

Щелкните на изображении чтобы закрыть окно

Рисунок 8 – Результат оптимизации

Так для заданной детали в виде конуса оптимальными параметрами формовки получились радиус инструмента Rt=10 мм. и шаг между витками Δz=3 мм. С этими параметрами данный процесс укладывается в ограничения заданными оптимальными факторами ограничения.

Для подтверждения оптимальных параметров проведем математическое моделирование процесса с применением ПО PamStamp.

Математическое моделирование инкрементальной формовки

С помощью CAD пакета Siemens NX смоделирован конус, заготовка, формообразующий инструмент и опора (рис. 9а). С помощью CAM пакета Siemens NX получена управляющая программа (рис. 9б), которая понадобилась для математического моделирования в среде Pam-Stamp. Методика моделирования процесса инкрементальной формовки описана в статье [1].

Щелкните на изображении чтобы закрыть окно

Рисунок 9 – Результат оптимизации

а - управляющая программа, б – расчетная модель

Результаты математического моделирования процесса представлены на рисунке 10. Отклонения полученной детали от первоначально заданной панели (рис.9.в), при данном способе формовки получились в пределах ±0,2мм (рис.10.а). На рисунке 10б показано распределение толщины детали, которое соответствует закону синуса. На рисунке 10в показано распределение деформаций готовой детали. Усилия на приводы станка с ЧПУ сопоставимы с расчетными значениями, порядка 3кН, и находятся в пределах рабочего диапазона оборудования.

Щелкните на изображении чтобы закрыть окно

Рисунок 10 – Результат математического моделирования в PamStamp

а – отклонения, б – распределение толщины, в – деформации

Для ее изготовления была спроектирована и изготовлена оснастка. Инкрементальная формовка конуса производилась из заготовки материала АМцМ лист толщиной 1мм, а результаты представлены на рисунке 10. Результат формовки сходиться с результатом математического моделирования. Процесс проходил с полученными оптимальными параметрами.

Щелкните на изображении чтобы закрыть окно

Рисунок 11 – Результат инкрементальной формовки конуса

Заключение

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

Исследование процесса инкрементальной формовки как многофакторной, многокомпонентной системы методами многомерной начертательной геометрии помогут повысить качество и расширить номенклатуры выпускаемых изделий.

References
1. Ayusheev T.V. Modelirovanie protsessa izgotovleniya detalei metodom inkremental'noi formovki [Tekst] / T.V. Ayusheev, R.N. Bulychev // Problemy mekhaniki sovremennykh mashin. Materialy VII mezhdunarodnoi nauchnoi konferentsii. Tom
2. Izd-vo VSGUTU, 2018.-C. 104-108 2. Volkov, V.Ya. Mnogomernaya ischislitel'naya geometriya: monografiya [Tekst] / V.Ya. Volkov, V.Yu. Yurkov.-Omsk: Izd-vo OmGPU, 2008.-244 s
3. Volkov, V.Ya. Graficheskie optimizatsionnye modeli mnogofaktornykh protsessov: monografiya [Tekst] / V.Ya. Volkov, M.A. Chizhik.-Omsk: Izd-vo OGIS, 2009.-101 s
4. GOST 21631-76. Listy iz alyuminiya i alyuminievykh splavov. Tekhnicheskie usloviya [Tekst].-Moskva: Izd-vo standartov, 1977.-29 s
5. Gusev A.I. Optimizatsiya parametrov protsessa rotatsionnoi vytyazhki [Tekst] / A.I. Gusev // Issledovaniya v oblasti plastichnosti i obrabotki metallov davleniem, Tul'skii Politekhnicheskii Institut, 1981,-Vyp.8,-S. 101-105
6. Krivoshein V. A. Perspektivy ispol'zovaniya tekhnologii inkremental'noi formovki v sovremennom proizvodstve [Tekst] / V.A. Krivoshein, A.A. Antsifirov, Yu.V. Maistrov // Izvestiya vysshikh uchebnykh zavedenii. Mashinostroenie, 2014.-Vyp. 11 (656).-S. 84-89
7. Novik F.S. Optimizatsiya protsessov tekhnologii metallov metodami planirovaniya eksperimentov [Tekst] / F.S. Novik, Ya.B. Arsov // M.: mashinostroenie.-1980.-T. 304.-S. 42-45
8. OST 1 02507-92. Samolety dozvukovye. Obshchie trebovaniya k kachestvu vneshnei poverkhnosti [Tekst].-Vzamen OST 1 02507-84; vved. 1993-01-01-Moskva: Izd-vo standartov, 1993.-10 s
9. Radishchev V.P. O primenenii geometrii chetyrekh izmerenii k postroeniyu ravnovesnykh fiziko-khimicheskikh diagramm [Tekst] / V.P. Radishchev // Izv. SFKhA.-M., 1947.-T. 15-S. 129-134
10. Chizhik M.A. Graficheskie optimizatsionnye modeli mnogoparametricheskikh tekhnologicheskikh protsessov legkoi promyshlennosti [Tekst] / M.A. Chizhik, V.Ya. Volkov // Inzhenernyi vestnik Dona, 2012.-T. 21.-№. 3.-S. 87-94
11. Chumadin A.S. Issledovanie protsessov posloinogo deformirovaniya listovykh zagotovok s ispol'zovanie frezernogo stanka s ChPU [Tekst] / A.S. Chumadin, D.A. Baturin // Kuznechno-shtampovochnoe proizvodstvo. Obrabotka materialov davleniem, 2014.-Vyp.7.-C. 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 [Tekst] // Bull., Series D. – 2016.-T. 78.-Vyp. 1.-S. 173-184
13. Adams D. A new model for contact geometry in single-point incremental forming [Tekst] / D. Adams, J. Jeswiet // Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture, 2015.-T. 229. №. 6.-S. 982-989.-DOI: 10.1177/0954405414535580
14. Attanasio A. Optimization of tool path in two points incremental forming [Tekst] / A. Attanasio, E. Ceretti, C. Giardini // Journal of Materials Processing Technology, 2006.-T. 177.-Vyp. 1-3.-S. 409-412.-DOI: 10.1016/j.jmatprotec.2006.04.047
15. Azaouzi M. Tool path optimization for single point incremental sheet forming using response surface method [Tekst] / 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 [Tekst] / R. Bahloul, H. Arfa, H. Belhadjsalah // Key Engineering Materials, 2013.-T. 554-557.-S. 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 [Tekst] / S. Gatea, H. Ou, G. McCartney // The International Journal of Advanced Manufacturing Technology, 2016.-T. 87.-Vyp. 1-4.-S. 479-499. – DOI: 10.1007/s00170-016-8426-6
18. Guangcheng Z. Forming process of automotive body panel based on incremental forming technology [Tekst] / Z. Guangcheng, X. Jinbo, S. Xiaofan, Z. Xun, L. Chuankai // Metallurgical and Mining Industry, 2015.-T. 1. – Vyp. 12.-S. 350-357
19. Gulati V. Process parameters optimization in single point incremental forming [Tekst] / V. Gulati [i dr.] // Journal of The Institution of Engineers, 2016.-T. 97.-Vyp. 2.-S. 185-193.-DOI: 10.1007/s40032-015-0203-z
20. Hagan E. A review of conventional and modern single-point sheet metal forming methods [Tekst] / E. Hagan, J. Jeswiet // Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture, 2003.-T. 217.-Vyp. 2.-S. 213-225.-DOI: 10.1243/095440503321148858
21. Jeswiet J. Rapid prototyping of a headlight with sheet metal [Tekst] / J. Jeswiet, E. Hagan // Canadian Institute of Mining, Metallurgy and Petroleum, 2001.-S. 109-114
22. Li Y. Investigation and optimization of deformation energy and geometric accuracy in the incremental sheet forming process using response surface methodology [Tekst] / Y. Li [i dr.]. The International Journal of Advanced Manufacturing Technology, 2015.-T. 79.-Vyp. 9-12.-S. 2041-2055.-DOI 10.1007/s00170-015-6986-5
23. Moayedfar M. Parameter optimization and stretch enhancement of AISI 316 sheet using rapid prototyping technique [Tekst] / Moayedfar M. [i dr.] // IOP Conference Series: Materials Science and Engineering. – IOP Publishing, 2017.-T. 257.-Vyp. 1.-S. 1-13.-DOI: 10.1088/1757-899X/257/1/012006
24. Nimbalkar D.H. Review of incremental forming of sheet metal components [Tekst] / D.H. Nimbalkar, V.M. Nandedkar // International Journal of Engineering Research and Applications.-2013.-vol. 3,-issue 5.-S. 39-51.
25. Young D. Wall thickness variations in single-point incremental forming [Tekst] / D. Young, J. Jeswiet // Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture, 2004.-T. 218.-№. 11.-S. 1453-1459.-DOI: 10.1243/0954405042418400