INFLUENCE OF ANISOTROPY OF SHEETS FROM ALUMINUM-LITHIUM ALLOY 1441 ON THE SPRINGBACK AFTER STRETCH-WRAP FORMING
- Authors: Surudin S.1, Erisov Y.1,2, Petrov I.1, Razzhivin V.1
-
Affiliations:
- Samara National Research University named after Academician S.P. Korolyov
- Samara Federal Research Center of Russian Academy of Science
- Issue: Vol 21, No 5 (2019)
- Pages: 46-51
- Section: Articles
- URL: https://journals.eco-vector.com/1990-5378/article/view/88446
- ID: 88446
Cite item
Full Text
Abstract
Full Text
ВВЕДЕНИЕ Освоение новых процессов формообразования обтяжкой обшивок летательных аппаратов из алюминиевых сплавов требует надежных знаний о их механических свойствах, в том числе их анизотропии. Кроме того, необходимы сведения о влиянии реологических свойств деформируемых материалов на напряженно-деформированное состояние в очаге деформации, определяющее характер развития деформации и формирующее качество готового изделия. Одним из наиболее важных факторов процесса обтяжки является пружинение, которое влияет как на правильность геометрической формы оболочки, что важно для аэродинамических показателей, так и на процесс установки обшивки на каркас самолета (трудоемкость сборки планера самолета). Пружинение обшивки после обтяжки, в основном, является следствием таких факторов, как неодинаковое внешнее трение, неодинаковое реологическое поведение деформируемого материала, задаваемое анизотропией свойств листовых полуфабрикатов [1-3]. Однако, в работах последних лет, посвященных решению частных технологических задач по совершенствованию процессов формообразования обтяжкой вопрос о влиянии анизотропии механических свойств на пружинение затрагивается слабо, либо не затрагивается совсем. Так в работах [4-6] описываются особенности применения упругих элементов при формообразовании обтяжкой. Дифференциальной обтяжке с контролем зоны наибольшего напряжения и деформации посвящены работы [7-8]. Учет анизотропии свойств при построении диаграммы предельного деформирования применительно к процессу обтяжки показан в работе [9]. В публикациях [10-12] акцент направлен на описании особенностей конструкции обтяжного пресса и его зажимных устройств. Моделированию операций формообразования обтяжкой с учетом перемещения зажимных устройств и положения стола пресса посвящены работы [13-14]. Кроме того в работах [15-16] рассматриваются вопросы учета влияния пружинения при применении многосекционного пуансона. В работе [17] изучено влияние анизотропии свойств на разнотолщиность обшивки после обтяжки. Более того, на производстве в качестве входного контроля, как правило, учитываются только механические свойства материалов, и практически не учитывается анизотропия этих свойств. В связи с этим в данной работе проведено исследование влияния анизотропии свойств листовой заготовки на пружинение после обтяжки. МЕТОДИКА МОДЕЛИРОВАНИЯ Моделируемая схема обтяжки осуществляется в следующей последовательности. Заготовка предварительно обертывается вокруг обтяжного пуансона, с последующем поднутрением, а затем под действием подъема пунсона подвергается деформированию. Данная схема обтяжки может быть реализована на прессе, оснащенном системой автоматизации, например на прессе типа FEKD, за счет синхронного движения левого и правого балансиров с зажимами и стола пресса с установленным на нем обтяжным пуансоном [18-19]. Данная кинематическая схема была смоделирована в программе PAM-STAMP 2G следующим образом. Для элементов заготовки в области зажимных губок задавалось перемещение, соответствующее кинематике пресса. После оборачивания заготовки вокруг пуансона моделировалось его перемещение, при этом зажимные губки фиксировались в неизменном положении. Геометрия обшивки приведена на рисунке 1. Радиусы линий, пересекающихся в точке О соответственно равны: продольный по формообразующему контуру =766 мм, поперечный по второй кривизне =9187,6 мм. Прогиб поперечного контура оболочки в точке О равен f0 = 71 мм. Длина формообразующего контура обтяжного пуансона 2L равна 2616 мм. В конечно-элементной модели использовались 4-х узловые оболочечные конечные элементы с 5 точками интегрирования по толщине заготовки. Размеры листовой заготовки равны 2620х2280 мм, толщина 2 мм. Для описания поведения материала заготовки - алюминий-литиевого сплава 1441 - использовалась модель ортотропного упругопластического материала (Hill48 model) [20]. Направление прокатки совпадало с направлением обтяжки. Для описания анизотропии использовались коэффициенты поперечной деформации mij, представляющие собой отношение логарифмической деформации по ширине образца eb к деформации по его длине el: , (1) где индекс i соответствует направлению поперечной деформации при линейном растяжении образца вдоль направления j [3]. Упрочнение материала в ходе пластической деформации подчинялось закону Свифта [21]: , (2) где σeq - интенсивность напряжений; εeq - интенсивность деформаций; ε0 - деформация начала пластического течения; n - степень упрочнения; k - коэффициент упрочнения. Данные показатели определялись через механические свойства: ; ; , (3) где - предел прочности, - предел текучести, - равномерное относительное удлинение. ПЛАНИРОВАНИЕ ЭКСПЕРИМЕНТА Для исследования влияния анизотропии свойств на пружинение применялось центральное композиционное планирование, которое включает полный и дробный факторный эксперименты и некоторое число дополнительных опытов, зависящее от числа факторов. В качестве переменных факторов модели использовались механические свойства материала: предел прочности, предел текучести, равномерное относительное удлинение, коэффициенты поперечной деформации. Уровни факторов и интервалы варьирования приведены в таблице 1. План эксперимента состоял из 46 сочетаний факторов. В качестве отклика использовалось величина пружинения, которая оценивалась, как среднее расстояние между точками кромок обшивки до и после пружинения (рисунок 2). Рассматривалось только центральное сечение вдоль направления обтяжки. Необходимо отметить, что перемещение участков оболочки, находящихся в зажимах пресса (часть оболочки, лежащая ниже горизонтальной плоскости, рисунок 3) не анализировалось. РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ После моделирования всех вариантов обтяжки был проведен статистический анализ результатов и определены коэффициенты математической модели. Разработанная регрессионная модель (в кодированных значениях) зависимости перемещения кромок обшивки S от анизотропии механических свойств имеет следующий вид: (4) где - кодированные значения переменных факторов (таблица 1). По данной модели были построены эпюры распределения величины пружинения при различных значениях механических свойств (рисунок 4-5, значения пружинения указаны в мм). Анализируя полученные зависимости, можно отметить, что механические свойства влияют на пружинение не значительно (рисунок 4) в отличие от показателей анизотропии (рисунок 5). Это также подтверждается статистической значимостью коэффициентов модели (4), исходя из которой в модели можно пренебречь пределами прочности и текучести. В целом, при увеличении относительного удлинения величина пружинения уменьшается. Уменьшение пружинения будет наблюдаться при обтяжке листов, показатели анизотропии которых больше 0,5 и равны между собой (рисунок 5). При низких значениях показателей величина пружинения увеличивается, несмотря на высокую пластичность заготовки. При этом наименьшее пружинение будет в том случае, если деформирующее усилие приложено вдоль направления максимального показателя анизотропии (1 и 21> 0,5), а показатель анизотропии в поперечном направлении обтяжки - наименьший (12≈0,5). Используя известные методы [22] поиска глобального минимума функции (4), была определена оптимальная анизотропия механических свойств, обеспечивающая минимальное пружинение мм для данной схемы обтяжки обшивок из алюминий-литиевого сплава 1441: = 18%, 21= = 0,65, 12= 0,45 (величина sВ и s0,2 значения не имеют. Стоит также отметить, что величина пружинения и разнотолщинность обшивки после обтяжки связаны между собой: чем меньше разнотолщинность, тем меньше и величина пружинения. Данный вывод подтверждается результатами работы [17], в которой было установлено влияние анизотропии свойства на разнотолщинность. ЗАКЛЮЧЕНИЕ В результате проведенной работы были найдены оптимальные значения механических свойств листов из сплава 1441, минимизирующие пружинение. Определено, что для минимизации пружинения листовую заготовку необходимо так ориентировать относительно стола пресса, чтобы направление обтяжки совпадало с направлением максимального показателя анизотропии, а поперечное направление обтяжки совпадало с минимальным показателем анизотропии. На предприятиях необходимо обеспечивать выходной контроль не только по механическим свойствам, но и по коэффициентам поперечной деформации, поскольку именно они оказывают наибольшее влияние на получение необходимой формы изделия. Рис. 1. Оболочка двояковыпуклой формы: а - поверхность оболочки двойной кривизны; б - формообразующий продольный контур; в - максимальный поперечный контур Таблица 1. Уровни и интервалы варьирования факторов Рис. 2. Схема к определению величины пружинения: 1 - заготовка после пружинения; 2 - заготовка до пружинения Рис. 3. Расстояние между обшивкой до и после пружинения, мм (= 405 МПа, = 280 МПа, = 14%,= = 0,65, = 0,45) Рис. 4. Распределение величины пружинения в зависимости от механических свойств при постоянных показателях анизотропии (= 0,5): а - = 14%; б - = 16%; в - = 18% Рис. 5. Распределение величины пружинения в зависимости от уровня показателей анизотропии при постоянных механических свойствах (= 417,5 МПа, = 302,5 МПа, = 16%): а - = 0,35; б - = 0,5; в - = 0,65About the authors
Sergei Surudin
Samara National Research University named after Academician S.P. Korolyov
Email: innosam63@gmail.com
Candidate Degree in Engineering, Associate Professor of Department of Metal Forming of Samara University
Yaroslav Erisov
Samara National Research University named after Academician S.P. Korolyov; Samara Federal Research Center of Russian Academy of Science
Email: erisov@ssau.ru
Candidate Degree in Engineering, Associate Professor of Department of Metal Forming, Engineer, Samara Federal Research Center of Russian Academy of Science
Ilya Petrov
Samara National Research University named after Academician S.P. Korolyov
Email: ilpetrof110895@yandex.ru
PhD Student of Metal Forming Department of Samara University
Vasiliy Razzhivin
Samara National Research University named after Academician S.P. Korolyov
Email: vasia.razzhivin@yandex.ru
PhD Student of Metal Forming Department of Samara University
References
- Dem’yanenko E.G. A technique of shaping the barrel-type parts // Russian Aeronautics, 2014. 57 (2), pp. 204-211.
- A manufacturability improvement of glass-fiber reinforced aluminum laminate by forming an effective crystallographic texture in V95 alloy sheets / F.V. Grechnikov, V.V. Antipov, Y.A. Erisov, A.F. Grechnikova // Russian Journal of Non-Ferrous Metals. 2015. 56 (1), pp. 39-43.
- Erisov Y.A., Grechnikov F.V., Surudin S.V. Yield function of the orthotropic material considering the crystallographic texture // Structural Engineering and Mechanics. 2016. Т. 58. № 4. С. 677-687.
- Park J.-W. , Kim J., Kang B.-S. Study on multiple die stretch forming for curved surface of sheet metal // International Journal of Precision Engineering and Manufacturing. 2014. V. 15, Issue 11. pp. 2429-2436.
- Seo Y.-H., Kang B.-S., Kim J. Study on relationship between design parameters and formability in flexible stretch forming process // International Journal of Precision Engineering and Manufacturing. 2012. 13(10). pp. 1797-1804.
- Numerical simulation on the local stress and local deformation in multi-point stretch forming process / S. Wang, Z. Cai, M. Li, Y. Lan // International Journal of Advanced Manufacturing Technology. 2012. 60 (9-12). pp. 901-911.
- Numerical simulation of different clamping modes on stretch forming parts / X. Chen, M.Z. Li, W.Z. Fu, Z.Y. Cai // Advanced Materials Research. 2011. №189-193. pp. 1922-1925.
- Minimum deformation path sheet metal stretch-forming based on loading at discrete points / Z.-Y. Cai, Z. Yang, C.-J. Che, M.-Z. Li // International Journal of Advanced Manufacturing Technology, 1 April 2016, pp. 1-10.
- Sheet metal forming limits under stretch-bending with anisotropic hardening / J. He, Z.C. Xia, X. Zhu, D. Zeng, S. Li // International Journal of Mechanical Sciences. 2013. №75. pp. 244-256.
- Формообразование профильных заготовок с помощью листового обтяжного пресса / А.А. Кривенок, А.В. Станкевич, С.И. Феоктистов, Р.Ф. Крупский, С.В. Белых // Ученые записки Комсомольского-на-Амуре государственного технического университета. Науки о природе и технике. 2013. №2. 1(14). С. 4-8.
- Мироненко В.В., Чеславская А.А.,. Белых С.В. Моделирование обтяжки летательных аппаратов с учетом эффектов, возникающих в зонах зажатия губками // Ученые записки Комсомольского-на-Амуре государственного технического университета. 2014. №2(18), С. 13-18.
- Малащенко А.Ю. Конечно-элементное моделирование процесса изготовления гибридных листовых деталей // Вестник Иркутского государственного технического университета. 2013. №4. С. 40-43.
- Моделирование кинематики движения рабочих элементов обтяжного пресса FET / Р.Ф. Крупский, А.А. Кривенок, А.В. Станкевич, С.В. Белых, В.В. Мироненко // Вестник Иркутского государственного технического университета. 2014. №9(92). С. 40-45.
- Определение положения пуансона в рабочем пространстве обтяжного пресса FET в процессе технологической подготовки производства / С.В. Белых, А.А. Кривенок, В.В. Мироненко, В.А. Мишагин // Вестник иркутского государственного технического университета. 2013. №12(83). С. 36-41.
- Numerical simulation of multi-point stretch forming and controlling on accuracy of formed workpiece / W. Liu, Y.-Y. Yang, M.-Z. Li // International Journal of Advanced Manufacturing Technology. 2010. 50 (1-4). pp. 61-66.
- Wang S., Cai Z., Li M. Numerical investigation of the influence of punch element in multi-point stretch forming process // International Journal of Advanced Manufacturing Technology. 2012. 49 (5-8). pp. 475-483.
- Сурудин С.В., Ерисов Я.А., Петров И.Н. Оптимизация анизотропии механических свойств обшивочных листов с целью минимизации разнотолщинности при обтяжке // Вектор науки Тольяттинского государственного университета. - 2017. - №2(40) - С.49-55.
- Statistical analysis of stretch shaping process of biconvex skin / V.A. Mikheev, S.D. Smol’nikov, S.V. Surudin, D.V. Savin // Russian Aeronautics. 2016. No.V. 59, Iss. 1.P. 145-150.
- Моделирование кинематической схемы последовательной обтяжки оболочек двояковыпуклой формы на обтяжном прессе FEKD / В.А Михеев. Ф.В. Гречников, С.Г. Дементьев, В.П. Самохвалов, Д.В. Савин, С.В. Сурудин // Известия Самарского научного центра Российской академии наук. 2014. Т. 16. № 6-1. С. 172-179.
- Hill R. The Mathematical Theory of Plasticity // Oxford, Clarendon Press, 1950.
- Gronostajski Z. The Constitutive Equations for FEM Analysis // Journal of Materials Processing Technology. 2000, V.106, pp. 40-44.
- Самарский А.А., Гулин А.В. Численные методы. М.: Наука, 1989. - 432 c.