УДК 621.98.042
АВТОМАТИЗИРОВАННЫЙ РАСЧЕТ ТЕХНОЛОГИЧЕСКИХ ПАРАМЕТРОВ ДРОБЕУДАРНОГО ФОРМООБРАЗОВАНИЯ КРУПНОГАБАРИТНЫХ ПАНЕЛЕЙ
© 2013 Ю.С. Андряшина
Национальный исследовательский Иркутский государственный технический университет
Поступила в редакцию 29.11.2013
Приведена методика расчёта технологических параметров дробеударного формообразования крупногабаритных панелей, и описана работа программного модуля.
Ключевые слова: дробеударное формообразование, приведенная толщина, удельный изгибающий момент, NetBeans (Java) Unigraphics NX
Автоматизация технологических процессов на производственном рынке занимает лидирующую позицию среди методов повышения производительности труда, улучшения качества продукции и оптимизации процессов управления. Одним из методов формообразования крупногабаритных листовых панелей самолетов является дробеударное формообразование (ДУФ). Этот метод заключается в обработке наружной (аэродинамической) поверхности детали потоком дроби. На Иркутском авиационном заводе (ИАЗ) применяется метод комбинированного формообразования, при котором продольная кривизна детали достигается гибкой на универсальном прессовом или валковом оборудовании, а для получения поперечной кривизны используется процесс ДУФ, который проводится в два этапа: предварительное формообразование и доводка. На сегодняшний день доводочная операция проводится в ручном режиме, а процесс предварительного формообразования реализуется с программным управлением, при этом образуется 70-90% от требуемой кривизны детали. Процесс ДУФ реализуется на дробеметной установке УДФ-4 контактного типа путем формирования на обрабатываемой аэродинамической поверхности детали полос обработки, высота которых определяется размером уплотнительного устройства дробемет-ного аппарата (ДА). После ДуФ производится зачистка обработанной поверхности абразивными лепестковыми кругами, ширина которых равна высоте полос обработки, создаваемых ДА [1].
Для выполнения автоматизированного расчета технологических параметров процесса ДУФ был создан программный модуль, который работает следующим образом. На поверхность электронной модели обрабатываемой детали, созданной в системе автоматизированного проектирования Siemens Unigraphics NX, в интерактивном режиме на основе опыта формообразования подобных деталей наносятся полосы обработки. Далее в автоматическом режиме происходит разбиение полос
Андряшина Юлия Сергеевна, научный сотрудник. E-mail: andryashina_ys@mail. ru
на прямоугольные участки, вычисление радиуса кривизны и приведенной толщины участков, на основе которых определяется изгибающий момент, потребный для формообразования. Под приведенной толщиной понимается толщина прямоугольного элемента сечения участка, осевой момент инерции которого эквивалентен толщине криволинейного сечения участка. Для интеграции системы с NX используются библиотеки интерфейса программирования NXOpen (Java). В качестве среды программирования применяется NetBeans (Java). Структура диалогового окна была создана посредством модуля UIBlock Styler системы Unigraphics NX. Диалоговое окно программного модуля представлено на рисунках 1а, 1б.
Цель расчета параметров ДУФ - определение режимов обработки дробью - частоты вращения ДА и времени обработки (расход дроби в ДА контактного типа является постоянным), потребных для предварительного формообразования i-го участка детали. Под участком детали подразумевается некоторая прямоугольная область с постоянными или усреднёнными значениями компонент кривизны и жесткости. Функциональные возможности Siemens Unigraphics NX позволяют вычислить приведенную толщину сечения участка (рис. 2) и радиус кривизны аэродинамической поверхности (рис. 3).
Автоматический расчет технологических параметров процесса ДУФ реализован с использованием зависимостей, приведенных в работе [3]. На основе полученных геометрических параметров участка детали вычисляется удельный изгибающий момент Mxi, потребный для формообразования поперечной кривизны рассматриваемого участка при ДУФ:
Mx = 0,09375kпk0EH3 кт
(1)
где ^ - коэффициент, учитывающий уменьшение требуемой кривизны при предварительном формообразовании, kП = 0,7-0,9; k() - коэффициент, учитывающий влияние ранее обработанных участков,
и определяемый опытным путём при формообразовании первой детали; Н - приведённая тол-
ЯТ - требуемое значение радиуса кривизны аэро-
динамической поверхности детали; г0 - коорди-
щина расчетного Участка детали в продольном на- ната центра тяжести сечения участка (расстояние правлении; кТ - требуемая поперечная кривизна от поверхности).
нейтрального слоя участка, кТ = 1 /(Щр - г0 ),
а)
Рассчение полосы обработки Выберите объекты А
* Аэродинамическая поверхность (0)
* Направление рассечения (0)
Укажите парам етры А
Ширина полосы
Отмена
б)
Рис. 1. Программный модуль «Определение параметров дробеударного формообразования»: а) диалоговое окно блока «Режимы обработки»; б) диалоговое окно блока «Рассечение полосы обработки»
Рис. 2. Расчет эквивалентного прямоугольного элемента с помощью UшgrapЫcs NX
Режимы обработки дробью, обеспечивающие достижение потребного момента Мх определяются на основе уравнения равновесия моментов внешних и внутренних сил при изгибе:
МХ = МДУФх
(2)
1ДУФ,
момента М дуФ^, обусловленного удлинением поверхностного слоя и вызывающего формоизменение детали при ДУФ.
Рис. 3. Расчет радиуса кривизны аэродинамической поверхности участка с помощью Unigraphics
Ж
где М - поперечная компонента изгибающего М
Для определения изгибающего момента
дyФj необходимо найти внутренние силовые факторы процесса - растягивающую силу РДУФ и расстояние 2ДУФ от обработанной поверхности до точки её приложения. Компоненты МдуФ^ связаны
с внутренними силовыми факторами следующей зависимостью
M
ДУФ j = РДУФ (z0k X J,k = x, У
(3)
Здесь координаты z0k - центра тяжести: z0x - поперечного сечения образуемого участка, z0y - продольных сечений участка (усреднённая). При постоянных значениях диаметра и расхода шариков
г, ДУФ
сочетание параметров гдуФ и zc соответствует режиму ДУФ - скорости дроби и времени обработки. Данные параметры определяются следующим образом.
Изготавливаются образцы-полоски из материала обрабатываемых деталей, которые обрабатываются на исследуемом режиме, представляющем сочетание частоты вращения ЫдУФ и скорости 5дуФ подачи ДА. Процесс обработки разбивается на несколько рабочих ходов одинаковой продолжительности. После каждого рабочего хода производится измерение стрелы прогиба fo6, а также удлинения ALo6 образцов при помощи индикаторных приспособлений. Производится аппроксимация экспериментальных данных эмпирическими зависимостями стрелы прогиба fo6 и относительного удлинения £об образцов (еоб =ALo6/Lo6), где Ьоб -длина образца) от времени обработки
fo6 = af [1 - eXP(-^ f Тдуф )]
ГДУФ'
Fs = as[1 - ехР(-ЬЕТДУф)]
(4)
(5)
Р = ЕН f
1oô E11oôFoô •
об Zo =
Н
2 f 3 f
Н
(6)
где Е - модуль упругости материала образца; Ьоб, Ноб - длина и толщина образца.
Для определения компонент силы, действующей на листовую деталь, используем соотношение РдуФх =РдуФу = 1,5 Роб. При одинаковых условиях обработки образцов и деталей и допущении о симметричном характере НДС при обработке дробью можно считать, что zДУФ = . Для определения параметров ДУФ, необходимых для достижения значения изгибающего момента Мх,
потребного для получения поперечной кривизны участка при заданной частоте ШдУФ вращения ДА, решаем задачу нахождения значений внутренних силовых факторов РдуФ и 2Д , при которых
выполняется условие (2). Для этого используем уравнения (4) и (5) с набором коэффициентов, хранящихся в базе данных и соответствующих заданному значению ЫдУФ. Результатом расчета является значение времени обработки ТдУФ, соответствующее частоте вращения ДА по которой определяется подача &дуФ, используемая при обработке
s
ДУФ
Т / T
' ТДУФ ' 1 ДУФ
(7)
где ЬдУФ - длина зоны обработки ДА.
Выбор оптимальных режимов обработки осуществляем, минимизируя суммарное машинное время технологической последовательности ДУФ-зачистка. Для этого методом перебора имеющихся в базе данных значений частоты вращения ДА решаем следующее уравнение
Тф-в = ТДУФ + Твач = ТЕ (Vdp ) = min ТФ
(8)
где Тф_з - суммарное машинное время сочетания «ДУФ-зачистка»; ТдУФ и Тзач - машинное время соответствующих операций; Удр - скорость дроби.
Машинное время Тзач, необходимое для зачистки обработанного при ДУФ участка, при условии одинаковой ширины полос обработки при ДУФ и зачистке, находим из выражения
где а, Ъ, с, а„ Ъ£, се - коэффициенты, определяемые на основе экспериментальных значений и
£об.
Значения коэффициентов а,, Ъ, с, а„ Ъе, се в уравнениях (4) и (5) определяются для всего технологического диапазона значений частоты ШдУФ вращения ДА, используемых при обработке деталей, и хранятся в базе данных вместе с измеренными значениями глубины На отпечатков, образующихся при ДУФ с частотой вращения ДА. Уравнения (4) и (5) дают возможность определить значения стрелы прогиба ,об и удлинения £об образца при заданном времени обработки, на основе которых можно найти удельную растягивающую силу Роб, действующую на образец-полоску, и коор-
об
динату 2с точки ее приложения
Т =
ваш
Т
a„
a
(9)
где Ь - длина обрабатываемого участка; ,?зач - подача при зачистке; азач - общий припуск, удаляемый при зачистке, азач =Н( - (0,02-0,03), мм; На -максимальная глубина отпечатка при ДУФ; а( -припуск, удаляемый при зачистке за один рабочий ход.
Результатом расчета является набор режимных параметров и ^д^, при котором достигается минимальное значение суммарного машинного времени Тфз совмещенной операции ДУФ-зачистка. Полученные данные обрабатываются программой постпроцессинга, генерирующей файл, содержащий управляющую программу для дробеметной установки. Эта управляющая программа представляет собой последовательность директив, выполняемых командным процессором с применением системы ЧПУ Smumeric 840Б 81 (рис.
4).
Разработанное программное обеспечение апробировано в производственных условиях Иркутского авиационного завода при изготовлении опытных образцов панелей крыла самолета БЕ-200. Завершение выполняемых в настоящее время работ позволит полностью отказаться от проведения дорогостоящих опытных работ на конструктивно-подобных образцах при запуске в производство новых изделий.
др
2
Программа постпроцессинга
N1 XO Y0
N2 M3 S600
N3 Gl XO Y399 F0.9587
N4 Gl XI35 Y399 Fl.2177
N5 Gl X271 Y399 Fl.2707
N6 Gl X406 Y399 Fl.2529
N7 Gl X541 Y399 Fl.0748
N8 Gl X677 Y399 F1.2477
N53 Gl X6357 Y397 F0.7003
N54 M30
Рис. 4. Результат работы программного модуля и фрагмент программы постпроцессинга
Представленная в рамках данной статьи работа проводится при финансовой поддержке правительства Российской Федерации (Минобрнауки России) по комплексному проекту 2012-218-03-120 «Автоматизация и повышение эффективности процессов изготовления и подготовки производства изделий авиатехники нового поколения на базе Научно-производственной корпорации «Иркут» с научным сопровождением Иркутского государственного технического университета» согласно постановлению Правительства Российской Федерации от 9 апреля 2010 г. №218.
СПИСОК ЛИТЕРАТУРЫ:
1. Вепрев, А.А. О создании отраслевой технологии дробеударного формообразования панелей / А.А.
Вепрев и др. // Авиационная промышленность. 2009. №2. С. 24-29.
Вепрев, АА. Автоматизация производства длинномерных панелей и обшивок на иркутском авиационном заводе /АА. Вепрев и др. // Наука и технологии в промышленности. 2013. №1-2. С. 49-52. Пашков, А.Е. Технологические связи в процессе изготовления длинномерных листовых деталей. -Иркутск: Изд-во ИрГТУ, 2005. 140 с. Пашков, А.Е. Математическое моделирование процесса формирования факела дроби в дробеметном аппарате контактного типа / А.Е. Пашков, В.В. Герасимов В.В. // Вестник ИрГТУ. 2011. №12. С. 4852.
THE AUTOMATED CALCULATION OF PEEN SHAPING TECHNOLOGICAL PARAMETERS OF LARGE-SIZE PANELS
© 2013 Yu.S. Andryashina National Research Irkutsk State Technical University
The method of calculation the peen shaping technological parameters of large-size panels is given, and operation of program module is described.
Key words: peen shaping, specified thickness, specific bending moment, NetBeans (Java) Unigraphics NX
Yuliya Andryashina, Research Fellow. E-mail: [email protected]