УДК 621.983.044.7.001.573
ВЛИЯНИЕ МАГНИТНО-ИМПУЛЬСНОЙ ОБРАБОТКИ НА ОБЪЕМ НЕСПЛОШНОСТЕЙ ДЕТАЛЕЙ ЛЕТАТЕЛЬНЫХ АППАРАТОВ
© 2017 Н.В. Курлаев, Н.А Рынгач, К.Н. Бобин
Новосибирский государственный технический университет
Статья поступила в редакцию 30.05.2017
В данной статье рассмотрено влияние характера импульсного нагружения на дефекты сплошности деталей летательных аппаратов из алюминиевого сплава.
Ключевые слова: математическая модель, магнитно-импульсная обработка, дефекты.
В технологических процессах изготовления деталей, связанных с пластическим деформированием, в структуре материалов возникает деформационная пористость. При определённых условиях обработки объем пор, пустот и микротрещин, не достигших критического размера, можно уменьшить импульсной обработкой давлением, в том числе, импульсным магнитным полем (ИМП) [1]. Однако параметры импульсных воздействий нуждаются в исследовании. При этом важно знать влияние основных факторов нагружения и способы их оптимизации для уменьшения пористости при импульсной обработке давлением.
Большое разнообразие методов и режимов обработки, в том числе магнитно-импульсной обработки (МИО), факторов воздействия, разброс исходных параметров материалов, разнообразие и сложность методов измерения полученных результатов затрудняют исследование механизма импульсного воздействия и, зачастую, приводят к противоречивым и даже отрицательным результатам. В процессах импульсной обработки давлением, в частности, при МИО необходим правильный выбор режимов для улучшения эксплуатационных показателей деталей. В противном случае может произойти снижение исходных эксплуатационных показателей. Для более точного исследования процессов импульсной обработки возникает необходимость проведения математического моделирования, в основном, с помощью численных методов с соответствующей оптимизацией технологических режимов по основным факторам: амплитуда, форма и время действия импульсной нагрузки, а также коэффициент пористости и динамический предел текучести материала.
Курлаев Николай Васильевич, доктор технических наук, профессор, заместитель заведующего кафедрой самолето- и вертолетостроения. E-mail: kurlaev@corp.nstu.ru Рынгач Николай Анатольевич, кандидат технических наук, доцент. E-mail: ryngach@corp.nstu.ru Бобин Константин Николаевич, кандидат технических наук, доцент. E-mail: bobin@corp.nstu.ru
Анализ процессов импульсной обработки производится на основе численного моделирования в рамках механики сплошной среды [1]. Используются уравнения, являющиеся следствием законов сохранения массы, импульса и энергии. Замыкается система соотношениями упругопластического поведения материала Прандтля-Рейса. В качестве критерия пластичности используется критерий Мизеса. При моделировании МИО электромагнитная сила ^ будет входить в уравнения движения.
vp=Ур; р = Gj j + f; рё = ау£у
p = p(р,ё); £j = оДуу + Vjt); stj +Asj = i^Ej
z-SjP+sj; e„S,
1=0
упругость; Я>0 - пластичность; stisti <— Y0
3
(1)
В системе уравнений (1): х. - координаты материальных частиц, V. - компоненты вектора скорости, г - плотность, V - объем, е - удельная внутренняя энергия ст.., 5.., (-.. - компоненты тензоров напряжения, девиатора напряжений и тензора скоростей деформаций соответственно, р - давление, 5.. - символ Кронекера, т -модуль сдвига, У0 - динамический предел текучести.
Численная реализация алгоритма для решения этих уравнений выполняется в рамках явной разностной схемы модифицированных методом типа Уилкинса [2]. В расчетной области строится разностная сетка из треугольных ячеек, в узлах которых определяются значения координат и скоростей, а в ячейках плотность, компоненты тензоров напряжений и скоростей деформаций, внутренняя энергия и т.д. Решение реализуется с помощью программного комплекса ККШ24 ИТПМ СО РАН.
В настоящей работе проводится численное моделирование импульсного воздействия на материал с порами, в частности, ударное воздействие и МИО давлением ИМП (рис. 1).
При ударе расчетной схемой является осе-симметричный объем (цилиндр) с круглой порой на оси симметрии (см. рис. 1,а). При моделировании импульсной нагрузки с помощью удара цилиндру задается начальная скорость v0 и он контактирует с жесткой преградой. Бо-
о
V
Уг
%
г
й Гг,
у-— 1
г
/ £
1 Ни
для полного смыкания поры. При этом отмечалось, что время смыкания поры (от начала до окончания ее пластического деформирования) для различных v0 требуется различное.
35 а. Упор, Ю'8 м3
1НТР
^_______
ь
— 10м/с 70м/с 110м/с 190м/с -230м/с
мкс
а) б)
Рис. 1. Расчетные схемы: а - при ударном воздействии; б - при МИО
ковая поверхность также является жесткой. От контактной поверхности цилиндра с преградой распространяется упругопластическая волна сжатия. При подходе к поре волна взаимодействует с ней, и происходит изменение объема поры.
При магнитно-импульсной обработке (см. рис. 1.б) импульсная нагрузка, моделирующая воздействие ИМП, в виде удельной объемной силы (.) сконцентрирована в скин-слое, от которого аналогично воздействию поршня передается импульсное воздействие на материал цилиндра. Также возникает упругопластическая волна сжатия, воздействующая на пору. При этом для МИО характерно бесконтактное воздействие и возможность регулирования амплитуды и времени импульсной нагрузки.
Задачи импульсной обработки (см. рис.1) рассматриваются для одной поры. Но можно ввести понятие пористости и использовать коэффициент пористости а, как характеристику материала:
а= ^ОР / ^ША , (2)
где VПОp - объем пор в расчетном элементе, VmIA - объем расчетного элемента. При этом используются: а0 - исходная пористость и аК - конечная пористость.
При решении задачи удара цилиндра из А/сплава (У0=300 МПа) с порой (а0« 0,015) о жесткую преграду определялось изменение объема поры (рис. 2). При малых скоростях воздействие упругое и объем поры не меняется. При увеличении v0 появляется изменение объема поры, причем остаточное. При достижении определенной скорости удара v0=vm пора смыкалась полностью. Скорость V является минимально необходимой
0 2 4 6 8 10
Рис. 2. Изменение объёма поры в А/-сплаве (а0=0,015, У0=300 МПа) при различных скоростях удара
При исследовании зависимости vm(а0) для абсолютного размера пор установлено, что для различных радиусов пор (гП) в телах с одинаковым коэффициентом пористости (а0 ) примерно одинакова и минимально необходимая скорость удара ^т). Результаты расчетов представлены на рис. 3, из которого видно, что осредненная зависимость vm от а0 практически линейная. 300 Ут,м/с
250 200 150 100
50 0
"гп=0,2 мм "гп=0,15 мм гп=0,125 мм Тп=0,05 мм
ап
0,004 0,008 0,012 0,016 0,02 0,024 0,028
Рис. 3. Зависимость V от а
т о
От момента удара тела с порой о жёсткую преграду до того как пора (поры) полностью или частично сомкнется и перестанет пластически деформироваться проходит определённый промежуток времени (см. рис. 2). Непосредственное время деформирования поры (всех пор), отсчитываемое от начала пластического деформирования до полного смыкания или окончания пластического деформирования, назовем временем смыкания и обозначим как V*.
Численные расчеты показывают, что минимально необходимая для полного смыкания скорость удара vm (см. рис. 3) и время смыкания
пор Ь* (рис. 4) зависят линейно от исходной пористости а0 и не зависят от размеров пор. Следует отметить, что это наблюдается как для одной поры, так и для множества пор (3, 6, 9 пор).
полного смыкания поры импульс давления в зоне поры достигает максимального значения (рис. 8). При дальнейшем увеличении у0 величина 1П уменьшается из-за большего влияния инерционных сил. Таким образом, максимальное значение 1П соответствует минимальной нагрузке необходимой для полного смыкания пор. Отметим, что это свойство можно использовать для постановки задач оптимизации.
При моделировании воздействия ИМП (см. рис. 1б) принимаем схему МИО с помощью плоского индуктора (рис. 9).
Электромагнитная сила
/г = 1 X В,
(4)
Рис. 4. Зависимость ¿*от а0
Как отмечалось в [1], при интенсивных нагрузках в импульсно нагружаемом элементе (см. рис. 1) шаровая часть (р) тензора деформаций (5.) будет много больше девиаторной части (хг?). Поэтому, часто при описании процесса в качестве характеристики напряженного состояния рассматривается только давление р=р(г,е).
Изменение давления в зоне поры (рп) при ударе (для v0 = 170м/с, У0=200МПа, а0« 0,01) показано на рис. 5, где более темные изобары соот-
как сила, действующая на единицу объема в скин-слое (см. рис1.а), будет входить в уравнение изменения импульса (р = О},} + ), где £ в общем случае имеет составляющие £х, £, £ .
Принимается, что сила ИМП действует по оси х, которая совпадает с осью симметрии цилиндрического элемента. Поэтому все составляющие £ , кроме одной (£), будут равны нулю.
В соответствии с известными соотношениями для индукции ИМП Б(х.уЛ) и плотности тока }(х,у,1) [26], -
В(х, у, г) = Вт ■ ехр^- д^ • si^OJt (5)
Рис. 5. Типовая динамика развития давления (рп) в зоне поры при ударе
ветствуют более высоким величинам давления. Фиксируя среднее значение давления в слоях на расстоянии < 0,1 г0 от границ поры в различные моменты времени Ь можем построить зависимость рП(Ь) (рис. 6). На основании зависимостей рП(Ь) для различных амплитуд нагрузки (0) и соответствующего времени смыкания пор (Ь*) (см. рис. 6) в качестве характеристики процесса смыкания пор можно ввести импульс давления в зоне поры (1П), который имеет вид:
и ](х,уЛ) =
Мод
■ехр
д,
х П
О--+ —
до 4
(6)
- уравнение электромагнитной силы £ имеет
вид:
Iп = I Рп (!) Ж
(3)
Импульс давления в зоне поры (1П) зависит от исходной пористости, практически линейно возрастая при увеличении а0 (рис. 7). Величина 1П также зависит от амплитуды нагрузки, в частности, от величины у0 при ударе. Первоначально при увеличении у0 величина 1П возрастает, при достижении у0=у - минимальной скорости для
Рис. 6. Зависимость развития давления
в зоне поры рп(Ь): Ь* и Ь2*- деформационное время смыкания (при скоростях удара v0 = 140м/с и v0 = 170 м/с, соответственно)
2500
2000 -
1500 -
1000
500
а
0 0,01 0,02 0,03 0,04
Рис. 7. Зависимость импульса давления смыкания пор ^ от а0
1600 п
1200
800
400
МПа х мкс
м/с
50 70 90 110 130 150 170 Рис. 8. Зависимость 1П(у0)
К
190
Индуктор
Обр абат ываемая деталь
Г 2J
шш
Пора
Ж
S?
S*
S3,И 2/Mo Гз,И ;
. Mo Im ПИ .
Bm =-
1m
\2WH CU2 LI2 —-; WH =-=-
L " 2 2
L = Lc + LH-3
, (8)
с геометрическими параметрами:
1Р = ЯИ - гИ; х <ё3, (9)
где у ЗИ - удельная электропроводимость материала заготовки или индуктора, ^ - амплитуда тока; WН - накапливаемая энергия заряда магнитно-импульсной установки (МИУ); L - общая индуктивность - индуктивность МИУ, LИ З -индуктивность системы «индуктор-заготовка»); пи - число витков индуктора; ^ - ширина токовой полосы индуктора.
Индуктивность системы « индуктор-заготовка» вычисляется в зависимости от магнитного объема ^^ по упрощенной формуле, предложенной в работе [3]:
М(Умпи (10)
L
и - з
■Элементарный объем
Рис. 9. Схема взаимного расположения детали, содержащей расчетный элементарный объем с порой, и шины плоского индуктора ^и и г11 - внешний и внутренний радиусы шины) при обработке ИМП
/х (х,уЛ) = ^ф-ехрР\~^гХ\ - + П) п^0* -, (7)
где Вп1 - максимум магнитной индукции на поверхности детали; dЗ - толщина скин-слоя заготовки; ш0 - магнитная проницаемость вакуума (шд = 1,257 мкГн/м); ю - частота разрядного контура.
Толщина скин-слоя заготовки или индуктора ^ЗИ), другие величины, входящие в уравнение (7), определяются известными электрическими соотношениями для МИО [3]:
Магнитный объем ^^ определяется из следующего выражения:
¥м = (Ив + 8И + 83 )БИ, (11)
где hВ - толщина воздушного зазора между индуктором и заготовкой, а SИ - площадь обработки индуктором
В качестве характеристики магнитно-импульсного нагружения используем удельную энергию МИУ:
= . (12)
Аналогично ударному нагружению при МИО существуют минимально необходимые нагрузки ИМП, при которых достигается полное смыкание пор. Типовой пример изменения динамики объема поры (а0«0,01) для Al-сплава (У0=200МПа) представлен на рис. 10, отражающем динамику сокращения объема поры во времени при различной энергии МИО. Расчеты производились для индуктора с числом витков пи =1, с шириной токовой полосы ^=20 мм и с площадью обработки SИ « 22,4 см2.
Изменение давления в зоне поры при ударе и при МИО аналогичны друг другу. Зависимости рП(0 для различных амплитуд, времени действия нагрузки при различных методах импульсной обработки, - ударе и МИО, - показаны в виде сравнительных графиков на рис. 11.
На рис.11. представлено: во-первых, развитие давления в зоне поры во времени ф^)) на тело с пористостью а0=0,01 при различных скоростях удара У0=140м/с - ряд 1, и v0=170м/с
l
I, МКС
Рис. 10. Динамика изменения объема поры
(а0 «0,01) при МИО: 1 - №*=3,21Дж/мм2; 2 - w*=3,57,Дж/мм2; 3 -w*=3,97Дж/мм2; 4 - w* «4,53Дж/мм2
3500 п 3000 -2500 2000 -1500 -1000 -500 -0
■<л>
МПа х мкс
* МКС
о
10
Рис. 11. Развитие давления в зоне поры при ударе: ряд 1 - 140 м/с; ряд 2 -170 м/с; и при воздействии ИМП: ряд 3 - Ь =20 мкс (при ^н=10,2кДж); ряд 4 - Ь =20мкс (при 1^я=13,4кДж); ряд 5 - Ь =36мкс; ряд 6 - Ь =50мкс
- ряд 2. Во-вторых, сравнивается МИО при различных энергиях заряда МИУ: ^Н=10,2 кДж
- ряд 3 и Р^Н=13,4 кДж - ряд 4, при одинаковом времени воздействия ИМП (1=20 мкс). В-третьих, для МИО показаны графики рП(Ь) при различном времени действия ИМП (полное смыкание пор во всех случаях достигается): ряды 3 и 4 - при Ь=20 мкс; ряд 5 - при Ь=36 мкс, ряд 6 - при Ь=50 мкс.
Пользуясь данными р(), представленными на рис. 11, и пользуясь соотношением (3) построим зависимость импульса давления в зоне поры от времени смыкания Ь* (рис. 12) для величин нагрузок минимально необходимых для полного смыкания пор.
При построении зависимости импульса давления в зоне поры 1П во время её смыкания Ь* (см. рис. 12) определено, что для более растянутого времени импульсного воздействия величина импульса давления в зоне поры 1П(Ь*) требуется больше. Это говорит об уменьшении влияния инерционных сил. Это характерно как для ударного воздействия, так и для МИО.
Очевидно, что импульс давления в зоне поры (1П) является характеристикой, объеди-
Рис. 12. Зависимость импульса давления в зоне поры 1П от времени её смыкания Ь*
няющей разные процессы импульсного на-гружения. Зависимость 1П(Ь*) (см. рис. 12) для импульсных нагрузок, минимально необходимых для полного смыкания пор, является единой для удара и МИО. Таким образом, по импульсу давления в зоне поры можно определить характеристику процесса независимо от вида нагружения.
При построении зависимости IП^*) (рис. 13) определено, что аналогично зависимости IП(у0) (см. рис. 5), импульс давления в зоне поры сначала растет при росте амплитуды импульсной нагрузки, а затем начинает убывать.
При этом максимальное значение 1П соответствует минимальной нагрузке необходимой для полного смыкания пор при МИО. Это свойство импульса давления в зоне поры (1П), как характеристики процесса импульсной обработки, можно использовать для постановки задач оптимизации требуемой величины энергии заряда накапливаемой МИУ.
3000 -,
2500
2000
1500
1000
500 -
1П'
МПа X мкс
w*, Дж/мм2
Рис. 13. Зависимость 1П для МИО А/-сплавов при 1 « 20 мкс:
1 - У0=100 МПа; 2 - У0=200 МПа; 3 - У0= 300 МПа
Для определения энергетических режимов МИО, требуемых для уменьшения пористости, предлагается инженерная методика сведения результатов численных расчетов и оптимизации амплитуды импульсной нагрузки в номограмму (рис. 14). По номо-
Рис. 14. Номограмма определения энергетических режимов обработки ИМП типовой МИУ (с собственной частотой £с~20(25кГц) для различных сплавов сплавов (гЭЛ, У0) с различной исходной пористостью (а0) и заданным коэффициентом смыкания (ц*= 1 - аК/а0) при обработке индукторами с числом витков пИ и шириной токовой полосы I
свойств обраба-Уп и удельного
грамме в зависимости от тываемого материала (а0, электрического сопротивления - гЭЛ), параметров индуктора (пИ и 1р), требуемой степени смыкания пор (ц* = 1 - аК/а0) определяется удельная энергия МИУ, приведенная к площади обработки индуктором
Таким образом:
- при численном моделировании импульсной обработки давлением, применяемой для уменьшения объема пор в материалах, определено, что результаты расчетов не зависят от размера и количества пор. Решая задачи об одной поре можно обобщать результаты через коэффициент пористости а;
- введен параметр - импульс давления в зоне поры (1П), который позволяет получать характеристики процесса смыкания пор независимо от способа нагружения;
- зависимость изменения импульса давления в зоне поры от амплитуды импульсной нагрузки имеет экстремум. Максимальная величина 1П соответствует минимальной величине нагрузки, необходимой для полного смыкания
пор, что является основой для постановки задач оптимизации;
- результаты численного моделирования и оптимизации параметров импульсной нагрузки можно обобщить в виде номограмм для инженерных расчетов по определению требуемых технологических режимов.
СПИСОК ЛИТЕРАТУРЫ
1. Моделирование воздействия ИМП на несплошности в структуре материалов деталей / Н.В. Курлаев, А.И. Гулидов, В.Б. Юдаев, Н.А. Рынгач // Проблемы машиностроения и надежности машин. 2003. № 2. С.80-86.
2. Влияние импульсной обработки на дефекты сплошности и долговечность материалов / Н.В. Курлаев, А.И. Гулидов, Л.А Мержиевский, Н.А Рынгач // Научный вестник НГТУ. 2005. № 1(18). С. 97-110.
3. Белоусов В. С. Физические основы лазерной и магнитно-импульсной обработки: Учебное пособие. Новосибирск: НЭТИ, 1991. 63 с.
THE INFLUENCE OF ELECTOMAGNETIC PULSE TREATMENT ON THE VOLUME DISCONTINUITY IN AIRCRAFT PARTS
© 2017 N.V. Kurlaev, N.A. Ryngach., K.N. Bobin
Novosibirsk State Technical University
In this article considers the influence of the nature of impulse loading on the defects of continuity of aircraft parts made of aluminum alloy.
Keywords: mathematical model, magnetic-pulse processing defects.
Nikolai Kurlaev, Doctor of Technics, Professor, Deputy Head at the Aircraft and Helicopter Manufacturing Department. E-mail: kurlaev@corp.nstu.ru
Nicholai Ryngach, Candidate of Technics, Associate Professor. E-mail: ryngach@corp.nstu.ru Konstantin Bobin, Candidate of Technics, Associate Professor. E-mail: bobin@corp.nstu.ru