Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Основы электронно-лучевой обработки материалов

..pdf
Скачиваний:
8
Добавлен:
12.11.2023
Размер:
18.54 Mб
Скачать

кое поле всегда имеет максимум, расчетное температурное рас­ пределение монотонно, но с увеличением диаметра электронного луча кривые обнаруживают тенденцию к формированию макси­ мума.

Это объясняется ростом крутизны распределения тепловых источников вдоль оси (от поверхности к максимуму) и углубле­ нием максимума с увеличением диаметра электронного луча. В модели Арчарда соответствующие параметры не зависят от диаметра, но всегда достаточны для формирования максимума

стационарного

 

распределения

 

 

 

температур

(<7тах/<7пов ^

3,2;

 

 

Таблица 23

^тах

0,44/*о).

 

 

 

 

Сравнение температуры

 

 

 

 

на поверхности [29] с рассчитанной

В

гауссовой

модели

нали­

 

по модели Арчардан

чие

максимума

стационарного

 

 

 

распределения

 

температур

на

D/ r0

Т [29]

т

глубине

обусловлено

отсут­

 

 

 

ствием разрыва на поверхности

 

 

 

(бесконечная среда).

стаци­

0,25

270

279

Сравнение

значений

0,5

221

223

онарных

температур на

по­

0,7

195

196

верхности

[29 ] и рассчитанных

 

 

удовлетвори­

по модели Арчарда (табл. 23) свидетельствует об их

тельной сходимости. С уменьшением диаметра электронного луча их расхождение увеличивается.

Приведенные результаты расчетов температурных распределе­ ний в обрабатываемом материале показывают, что в некоторых случаях, например при малых диаметрах электронного луча {Dirо < 1) или когда необходима информация о максимуме рас­ пределения, расчеты следует осуществлять с использованием более точного пространственного распределения энергетических потерь электронов.

Распределение температур в материале при воздействии элек­ тронного луча можно также получить экспериментально, например по данным металлографического анализа структур при сварке [201 ]• На рис. 67 представлена зависимость Т (х, /) —■Т = f (X2) для ряда металлов.

Как видно, совпадение расчетов по модели мгновенного плос­ кого источника на поверхности полубесконечного тела с экспери­ ментальными данными удовлетворительное.

СВЯЗЬ ПАРАМЕТРОВ ЭЛЕКТРОННОГО ЛУЧА С ХАРАКТЕРИСТИКАМИ ГЛУБОКОГО ПРОПЛАВЛЕНИЯ В МЕТАЛЛАХ

Если тепловые процессы в зоне обработки при воздействии электронного луча малой интенсивности хорошо описываются тео­ рией теплопроводности, то расчет температурного поля при глу-

121

 

боком

проплавлении или

свер­

 

лении

отверстия

(при

плотно­

 

сти потока <72 > 105-ь 106 Вт/см2)

 

в рамках

теории

теплопровод­

]Z!

ности

затруднителен,

так

как

V.

возникает необходимость учета

потерь

энергии

на

фазовые

 

переходы (в основном на ис­

Рис. 68. Модель проплавления для рас­

парение).

из первых

попыток

чета оптимальной погонной энергии

Одна

U91]

получить

связь

глубины

про­

 

плавления с параметрами

элек­

тронного луча и теплофизическими характеристиками материала

сделана в работе

японских исследователей [192].

[192], ориенти­

Глубина проплавления Я, согласно работам

ровочно определяется

выражением

 

н -

1

IU

[(СТПЛ+ ^пл)Р+

>

2,1

vd

где dT/dR — температурный градиент в направлении R границы плавления, град/см (рис. 68). Отсюда погонная энергия, необ­ ходимая для получения глубины проплавления Я,

 

~

= 2,1Hd [ (сТпл 4- Ln>

+

- L ] .

(179)

 

Приближенное значение величины градиента

 

 

 

я

Н

т +

i s - ) 7' -

 

<180>

где

а — коэффициент

температуропроводности.

 

 

Подставляя

выражение

(180)

в

(179), получим

 

я /

2,1Я dp(c7'njI -f- Lnjl) [ 1

( 4 - + - J - )

_____ 7 пл____

v

Р {°ТЛ Л +

7 пл)

 

 

 

 

 

 

 

 

(181)

Из формулы (181) следует очевидный факт: для расплавления металла кроме тепловложення Hd (срТпл + рЬПЛ) необходимо вво­ дить дополнительную теплоту, возмещающую потери на тепло­ проводность. Следовательно, термический КПД

%

1______

 

(182)

5МГПл (1 jdv - f - 1 /2 а )

|

 

Р ( ^ 7 Пл + 7 П Л )

J

 

С увеличением v н а н уменьшением X термический КПД уве­ личивается.

122

 

 

 

 

 

Рис. 69. Модель цилиндри­

 

 

 

 

~7Г

Т"

(

 

л

ческой полости

в мате­

 

 

 

 

 

г

Ь

1

риале

при

глубоком

про­

 

 

 

 

ч;

 

 

 

 

плавлении

(а) и схема от­

 

 

 

 

 

 

 

 

 

вода

тепла

от

передней

 

 

 

 

 

 

 

 

 

стенки полости (б):

 

 

 

 

 

 

 

 

 

 

 

1 — полость;

2

— стенка

 

 

 

 

 

 

 

 

 

жидкого металла;

3 — объ­

 

 

 

 

 

 

 

 

 

ем

расплавленного металла

 

 

 

 

 

 

 

 

 

при

 

образовании

полости

 

 

 

 

 

 

 

 

 

[214]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис.

70.

Система

координат

при

перемещении

 

 

 

 

 

полости в

материале [214]

 

 

 

 

 

 

 

 

 

 

При

 

диаметре

электронного

луча

 

 

 

 

 

2 мм и скорости

сварки v =

0,5 ч-4,0

см/с

 

 

 

 

 

для

алюминия т)7- == 3 + 20% ,

для углеро­

 

 

 

 

 

дистой

и

нержавеющей

сталей

t]T=

 

 

 

 

 

=

14 + 25% ,

для

титана

 

и

циркония

 

 

 

 

 

г\т=

10 + 25%

[192].

 

 

 

 

С учетом

величины эффективного

К П Д электронно-лучевого

нагрева т]и и процесса испарения получим формулу для

глу­

бины проплавления

[219]

 

 

 

 

 

 

 

 

 

 

 

Н.

0,24/1/ (1 — 0,7ft -

 

1,09 - 10 -9Z l/) —

 

( Т е +

273,15)*

(183)

0 , 5 d v

[р (с Т е -|- L n j l ) -f- k L wспр + 5 Х (0,5а -j- 1 / v d )

Т пл]

 

 

 

где [3 — отношение числа отраженных электронов к общему числу излучаемых; Z — порядковый номер элемента; температура ванны жидкого металла; k — отношение количества испаривше­ гося материала к количеству расплавленного (k «=* 0,01); е — коэф­ фициент черноты-, о 0— постоянная Стефана—Ъолъцмана.

Величина эффективного КПД в формуле (183)

Ли = 1 -------------

1,09-10-4 U ,

(184)

где Е — отношение средней энергии отраженных электронов к энергии первичных электронов (Е ^ 0,7).

Формулу (183) можно упростить [219]:

t j _

r\a ( J U / v d ) + A / v

 

 

В +

С /v

 

где обобщающие коэффициенты

 

 

4 = — 4,18cdeo0 ЦТЫ +

273,15)4;

В — 2,1 (cpfTпл ■

рЕпл

kpL„cn-{- 2,5срГпл);

С === 10,57/Гпл;

/ = Г

е/Гпл=

1,1 - 1,6.

Недостатком формулы (183) является необходимость экспери­ ментального определения температуры Те ванны жидкого металла и коэффициента k-

123

Hj2
I 2 D V (A ,8 )d 8
Рис. 71. Значения интеграла в выра­ жении (191) в зависимости от вели­ чины параметра А = vb/2a [214]

Использование модели, в кото­ рой среднюю температуру жид­ кого металла определяют по ве­ личине давления пара в канале, позволило более точно оценить величину потерь энергии на теплопроводность и получить более точную зависимость пара­ метров электронного луча от ха­ рактеристик глубокого проплав-

Рассмотрим сущность этой мо­ дели. На рис. 69 показана модель

цилиндрической полости в мате­ риале при глубоком проплавлении. Предполагается, что полость су­

ществует стационарно и на образование канала, т. е. передней части полости (рис. 69, б), уходит ровно половина введенной энергии электронного луча. Другая половина тратится на по­ догрев расплавленного металла и его перемещение в сторону, противоположную движению луча.

Количество теплоты, затраченной на плавление половинного объема полости и теряемой за счет теплопроводности, в соответ­ ствии со схемой, приведенной на рис. 62, составляет

Евв ~ <2ПЛ-f- QT.

Скорость перемещения расплавленного металла из передней части полости в заднюю (в см3/с)

V* - 2vHb.

Теплосодержание расплавленного металла

+ Pl -

где Тп — температура поверхности жидкого металла. Тогда ба­ ланс затрат введенной энергии

0,24-0,5Ю =

vbH [^ - ( Г п -

Тпл) -f Ln„] + QT.

(185)

Потери теплоты за счет теплопроводности (рис. 70)

 

 

я / 2

 

 

 

 

Я /2

 

 

QT = H

|

g(Q)(bdQ)

 

- 2

ЪН j

x [ ^ ( 0 ) ] r=6d0.

(186)

 

 

 

 

 

 

 

0

 

 

Полагая

V =

(Тпл Т)

и V0 =

(Тпл— Т х ), получим

диф-

ференциальное уравнение теплопроводности

 

 

 

д 2У

д 2У .

у

д у

_

1

д у

(187)

 

 

д у 2

' д г 2 '

а

д у

 

a

dt

 

124

Для

стационарного случая (при

t —* ос)

решение

(187)

имеет

вид [76]

 

 

 

 

 

 

 

 

 

 

 

 

 

V (0, А)

 

__

0 — A r cos

S

e

7» О4)

Кп (Ar) cos «0,

(188)

 

 

 

 

 

 

 

 

 

К„(А)

где А =

vbl2a\

г — г/Ь;

1п,

Кп — функции

Бесселя

первого и

второго

рода соответственно.

 

 

 

 

 

 

Градиент температуры

при г — 1

 

 

 

 

 

 

 

^ • ^ Ь И т т ) J

=

 

 

 

=

A icos 0 -f- е~л cos0

^

£rtcos п0 х

 

 

 

 

 

 

 

 

 

 

- Л=з0

 

 

 

 

 

XI

 

/»-1 (Л) + 7д+, И)

 

ДЛд (А)

 

(189)

 

 

 

 

 

 

 

 

 

 

 

 

где £„ =

1 при /t = 0 и Ел =

2 при п >

1.

 

 

 

Выражение (186)

с

учетом уравнения (189):

 

 

 

 

 

 

я/2

 

 

 

 

 

 

 

 

 

Qr =

26Я

|

-

Х(Тпл7-т.а>)

DV{Q, A)dQ =

 

 

 

 

 

 

О

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

л/2

 

 

 

 

 

 

 

= m

(Гпл -

Гео) }

2DV (0,

Л) (Ю.

 

(190)

 

 

 

 

 

 

 

 

о

 

 

 

 

 

Подставив выражение

(190) в уравнение (185), получим

 

 

 

Н

 

 

 

А[2К(Тп -

Т„А +

4рйГпл] +

 

 

 

0,24-0,5 IU

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Гя/2

 

 

_ 1

 

(191)

 

+

4 ^ пл

 

Т„)

}

2DV (0, A)dQ

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

Значения определенного интеграла в выражении (191) пред­

ставлены на рис.

71.

 

 

 

 

 

 

 

 

 

Для случая, когда температура поверхности металла равна

температуре его

плавления

(Тп =

Тпл)9 выражение

(191)

упро­

щается:

 

 

 

 

 

 

Гя/2

 

 

 

-1

 

 

Я

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

2DV (0, Л)сШ +(4раГпл)Л

.

0,24-0, ЬШ

Ч

Г п л - Г с э )

 

 

 

 

 

 

 

(192)

Среднее значение температуры поверхности жидкого металла в выражении (191) следует определять исходя из того, что давле­

125

ние пара в канале уравновешивает давление, обусловленное поверхностным натяжением жидкого металла и силой тяжести [214]:

Р > (Ро + Ра) = 4 ~ + ('5я -

Для случая полости в

железе давление пара р ^ 7,1 х

X 10“2 кгс/см2, Тп ^

2360° С, для алюминия р & 3,3 • 10"2 кгс/см2,

Тп & 1890° С (при

г = b =

0,05 см).

Сравним теперь между собой результаты расчетов по уравне­ нию (191) и выражению (181).

Перепишем выражение (181) в виде

1

Ш

_ _ _I™_ _ _1 Г

1

] .

(193)

2,1

i/d

Р (рТ пл ~\~^пл) J L

1+2,5V a+ БК!(Щ

 

 

Для случая, когда d =

4b, уравнение (193) представится в виде

-ЯГ = {2.1 [Р (сТпл +

1пл)] [SK +

8 + 2 , 5 % ) - ^ ] У

(194)

Сравнение

результатов

расчета по уравнениям (194)

и

(193)

показывает, что разница (в 1,5—2 раза) имеет место в основном при очень малых и очень больших значениях величины vd. Удо­

влетворительное совпадение соответствует значениям vb =

0,1-*-

-*-0,2 см2/с.

 

 

 

При малых значениях параметра vb выражение (191) можно

приближенно представить

в

виде

 

Н

 

■ " ( - а - ) - 0’877

,19„

0,24*0,5Д/

~

п к (Т пя — Т со) *

{

Результаты расчетов по уравнению (191) и экспериментальные данные для случая проплавления среднеуглеродистой стали не совпадают приблизительно на 3 — 25% для алюминиевых сплавов.

КОМБИНИРОВАННЫЕ МОДЕЛИ ТЕПЛОВОГО ИСТОЧНИКА ПРИ ГЛУБОКОМ ПРОПЛАВЛЕНИИ

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

126

Получение глубокого проплавления можно трактовать как результат дей­ ствия двух источников тепла: точеч­ ного на поверхности и линейного в пла­ стине, не имеющей теплоотдачи, поэтому параметры электронного луча могут быть связаны с глубиной проплавления следующим критериальным уравне­ нием [44]:

„ ( 1 - 0 .5 А )

Я = Cl "

(196)

где величина Сх ^ 1, а величина k оп­

ределена экспериментальной формулой

k — 0,68 (ЯГпл)0-15.

(197)

Рис. 72. Схема теплового воз­ действия электронного лта [162)

Глубину действия линейного источника ограничивают величи­ ной 6. При тепловой мощности луча Q мощность точечного источ­ ника принимают равной qx = kQ, линейного q2 (1 — k) Q, где k — коэффициент распределения энергии [162].

Распределение температуры в полубесконечном теле от дей­ ствия точечного источника мощностью q1 при t —>оо [153]:

 

<198)

где R =

]/V + у г + z2,

здесь х,

у, z — координаты точки в декартовой системе.

Линейный источник рассматривают как совокупность мгновен­

ных точечных источников мощностью dq = q jb d zx при условии

равномерного распределения тепловыделения в канале проплав­ ления на глубину б (рис. 72).

Изменение температуры от воздействия такого источника в про­ извольной точке М (х, у , г) полубесконечного тела с адиабатиче­ ской границей описывают уравнением

dT2(x, у, z, t) .

dq>2

e x p

х Ч -t/2 + (z— г')2

] +

ср (4яа/)3/2

4at

 

 

I

 

+ e x p

X2 + У2 + ( z +

г ')2 Y\ __

q2 dz'

 

 

 

J)

бrp (4naff!2

 

 

 

 

 

x ь

[ -

 

[ - — ’ ' V e + *,)' ])•

(199,

Упрощая выражение (199) для стационарного случая, когда dT не зависит от t, получим

dT2(xt у , z) = q2dzrl(8AnXRf).

(200)

127

Интегрируя выражение (200) по z', легко получить

т*(х ‘

У' г>= -j6

 

=

д-2

,л Г

(г + 6) +

Г г 2 +

(г + б)2 1

4яЛб

I

(г — 6)

V г2 +

(2 0 1 )

(г — б)2 J ’

где

R' == Y х2+ уг + (2 z ' f ; г — V х г + у2-

Это же выражение можно получить из уравнения (199), ин­ тегрируя его по V и г'.

Сучетом формулы (198) и уравнения (201) температура в точке

М(х, у, г) в стационарном режиме может быть представлена уравнением

Т {х, у, г) = 7, + Т 2 =

Qi

Яг

In ~ (г + б) +

Г г 2+ ( г + й)2 '

 

2nXR

4яЯб

. (г — б) +

V г2+ (z — б)2 J

 

 

 

 

(202)

Как видно на изотермах, рассчитанных по уравнению (202), глубина действия линейного источника и соотношение удельных мощностей линейного и точечного источников (qx/q 2) влияют на форму изотерм. Если доля линейного источника незначительна, то форма изотерм близка к получаемым от действия точечного источника. С увеличением доли линейного источника форма изотерм приближается к «кинжальной».

Температуру нагрева металла нормально-распределенным по­ током высокой концентрации можно рассматривать как суммар­

ный результат действия

двух источников теплоты: Т = Тф + Г 2,

где Т х — температура,

определяемая

действием поверхностного

нормально-распределенного источника;

Т2— температура, опре­

деляемая действием движущегося в

глубь металла нормально-

распределенного источника теплоты

[32].

Температуру поверхностного нормально-распределенного ис­ точника определяют по известным уравнениям [153]. Для опре­ деления температурного поля Т 2 от непрерывно-движущегося в глубь тела источника с произвольной ограниченной в простран­ стве и времени функцией теплового потока / (г, /) принимают схему полуограниченного тела с движущейся границей [32].

При перемещении со скоростью v по оси OZ плоской границы с нормально-распределенным источником теплоты уравнение теп­ лопроводности записывают в виде системы двух уравнений,

128

Jr-L #

г

1—\

V

/ 0 / Г

 

Щ

ФГ /

 

V

 

- - -

 

тт^1гооо°

 

Ч Н

 

 

 

1

- 2 мм

 

 

Z Г1500°

 

й)

 

 

 

Ъ)

 

Рис. 73. Температурное поле при воздействии электронного луча на молибден [ 162]:

а Q =

900, qx = 700,

Я* =

200

кал/с,

6 =

0,5

см,

К

=

0.77,

б —

Q = 900, ^

=

300,

Яг ~

 

600

кал/с,

6 =

0,5

см,

&™

0,33, b

Q ~

1500,

qx — 500,

=

1000

кал/с,

6 ~

= 0,5

см,

k =

0,33,

/ —5 — распределение

температур

в

плоскостях

x — q

при

г = 0,

0,125,

0,25, 0,375,

0,5,

0,6 см, а — Q =

2000; ^

=

500;

q2 =

1500

кал/см,

8 = 1

см,

& =

0,25,

/ — 3 —

распределение температур

в плоскостях

х —у

при г =

0, 0,1, 0,2,

0,3,

0,5,

0,7, 1,0, 1,0,

1,1

см

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5 II II Ры^алпн и др

129

соответственно для нижнего и верхнего полупространств в коор­ динатах, движущихся с границей

дт'2

_

 

а2П

,

1

дт:2

,

а2П \

,

а п

dt

~

( д г 2

+ ~Г ~ я Г +

 

 

 

 

 

дг

 

 

 

 

ап

 

/ а2П

,

j_ ап ,

а»ПЛ ,

.. а п

dt

 

 

I

*

г

ал

+

аг! / ^

у

аг2

с краевыми условиями

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при 2Х= 0;

 

 

 

П =

0

 

при 2i = оо;

 

 

П =

Гг;

а п

=

а п

 

при 2j — 22 = 0;

 

 

 

аг.

 

аг2

 

 

 

 

 

 

П = п

=

о

 

при t — 0.

 

 

(203)

(204)

(205)

После применения интегральных преобразований Ханкеля и Лапласа к системе уравнений (203) и (204) с краевыми усло­ виями (205) решения этой системы и последующего обратного преобразования получены уравнения для определения темпера­ туры впереди и позади движущейся границы [32]:

со

Т

f (г, 1 -1)

 

Т.т

 

 

 

2Х/а1/2

 

 

 

 

 

 

X

 

 

 

 

 

 

■erfс f

г1

I

 

 

\2 (a t)^

+

со

t

f(r, x —i)

 

 

 

 

■ -о ио

2Я1/2

 

 

X

 

г ехр (

г2 \

( n t ) 1

\

4ai

/

 

 

erfc (

г2

^

 

 

2 (Ш)1/2

гг,

2а

v

 

 

 

 

 

Р (~5Г

+ Т Н -)

х

2а 1/2

I0(sr)sdsdt‘,

(206)

 

 

 

 

 

 

 

 

 

 

X

и

/

UZo

.

v4

 

2а 1'2 еХР \

2а

+

~4а~

 

га1''2 /. / 0 (ы) s d s d t .

 

(207)

130