CC BY
26
ИЗВЕСТИЯ
ТОМСКОГО ОРДЕНА ТРУДОВОГО КРАСНОГО ЗНАМЕНИ ПОЛИТЕХНИЧЕСКОГО Том 87 ИНСТИТУТА имени С. М. КИРОВА 1957 г.
К РАСЧЕТУ РАДИУСА РАВНОВЕСНОЙ ОРБИТЫ
БЕТАТРОНА
М. Ф. ФИЛИППОВ
(Представлено научным семинаром физико-технического факультета)
Общие соображения о выборе основных параметров электромагнита бетатрона
Расчет электромагнита индукционного ускорителя электронов-—бетатрона, как и трансформатора или электрической машины, представляет задачу с возможными многими решениями. Но только некоторые из этих решений являются технически и экономически целесообразными.
В релятивистской механике соотношение между двумя основными параметрами бетатрона —радиусом равновесной орбиты (г0) и плотностью магнитного потока на орбите (В0) с некоторым приближением определяется известным уравнением:
А = 3.10~4 .г0В0 , (1)
где А — 0,51 полная энергия ускоренных электронов, а их кинетическая энергия в Мэв;
В0 плотность магнитного потока на равновесной орбите в гс\
га радиус равновесной орбиты в см.
Из бесконечного числа нар значений г0 и В0, удовлетворяющих уравнению, могут быть приняты только значения, отвечающие выполнению следующих условий:
1. Распределение магнитного потока в зазоре между полюсами должно обеспечить в процессе ускорения соблюдение так называемого бетатронного условия (Вк = 2 В0), где Вк — средняя плотность магнитного потока в круге равновесной орбиты.
2. Электромагнит бетатрона должен иметь возможно меньший вес.
3. Нагревание стали магнитопровода и обмоток электромагнита при длительной работе в номинальном режиме не должно превышать допустимого по нормам.
Чтобы получить возможно меньший радиальный размер полюсных сердечников, следует стремиться к меньшему значению г0, но при этом возрастает необходимая для получения заданной энергии ускоренных электронов плотность магнитного потока на равновесной орбите.
При больших значениях плотности магнитного потока в междуполюсном пространстве электромагнита бетатрона заметно сказывается влияние насыщения стали центральной части магнитопровода. Влияние насыщения приводит к нарушению бетатронного условия.
В данном случае нарушение бататронного отношения вызывает смещение орбиты ускоряемых электронов к внутренней стенке ускорительной камеры.
Кроме того, увеличение плотности магнитного потока в полюсных сердечниках приводит к возрастанию потерь от вихревых токов и гистерезиса. При неравномерном распределении плотности магнитного потока по сечению магнитопровода наличие больших потерь в стали может создать недопустимую разность по фазе между ведущим и ускоряющим магнитными полями бетатрона, что также вызывает нарушение бетатронного условия. Повышенные потери в стали создают чрезмерный нагрев магнитопровода.
Для выявления наиболее удачного решения при выборе радиуса равновесной орбиты г0, если пользоваться уравнением (1), необходимо рассматривать несколько вариантов и по конечным результатам расчета судить о пригодности того или иного из них.
Чтобы внести определенность в выборе радиуса равновесной орбиты, ниже предлагается расчетное уравнение, связывающее основные параметры и радиальные размеры междуполюсного пространства электромагнита бетатрона при условии соблюдения бетатронного отношения.
Вывод расчетного уравнения
На рис. 1 представлены профиль полюсных наконечников бетатрона и соответствующая кривая распределения плотности магнитного потока В —/ (г) в средней плокости междуполюсного пространства.
Ряс. 1
В рабочей зоне междуполюсного пространства, занятого полезным сечс-нием ускорителькой камеры (участок аб), кривая распределения плотности магнитного потока вдоль радиуса должна удовлетворять соотношению:
В„ г
63
В,
или Вх - В,
г —п ' О
Г п ' X
(-2
где В>„ и Дг - мгновенное значение плотности магнитного потока на радиусах г0 и гх, п - - показатель спадания магнитного поля по радиусу, который должен быть в пределах 1 > п > О, чтобы обеспечить радиальную и аксиальную фокусировки пучка ускоряемых электронов.
Показатель спадания магнитного поля п на разных радиусах в рабочей зоне ускорительной камеры имеет неодинаковые значения (рис. 1), но в кольце, ограниченном радиусами г0 и гс, магнитный поток которого учитывается расчетным уравнением, п можно считать постоянным и равным значению его на радиусе г0.
В рабочей зоне магнитный поток через элементарную кольцевую поверхность с радиусом /*г и шириной ¿¿г (рис. 1) будет равен:
ёФ — 2^гхВхс1г — 2пВог0п.гх1-п(1г,
откуда магнитный поток в кольце, ограниченном произвольными радиусами гх и г2 в пределах рабочей зоны
га
Ф12 = 2 1Г в о г0" I' Гл-1-". йг = - ^ г22-п ___г 2-П ^ . (3)
Чтобы выдержать бетатронное отношение, центральная часть полюсных .наконечников снабжается цилиндрическими выступами или между полюсами помещается несколько цилиндрических дисков (галет).
Действительное распределение плотности магнитного потока вдоль радиуса на рис. 1 представлено кривой абвгд. Для упрощения расчетов с достаточной для практики точностью действительную кривую можно заменить кривой абези, тем более, что возможная неточность, обусловленная таким допущением, легко корректируется после сборки электромагнита регулированием зазора между полюсами.
Полагая, что плотность магнитного потока в кольце, ограниченном радиусами га и гс (рис. 1), изменяется в соответствии с уравнением (2), а в круге радиуса гс она постоянна и равна Вс, найдем из уравнения (3) магнитный поток в круге равновесной орбиты
ф ^ 2 *Д0г0"
^(^«-^Лг^Д,. (4)
2 п V /
Используя, кроме полученного выражения, бетатронное условие, найдем
2 во = (го2 п - ^ )+— ве,
г0Ц2 — п) V / Гл-
от к уда
5« =
Вс
1 /, _\1 (5)
Н' Г"^ О]
где
(6)
Величина в0 зависит от радиального размера поперечного сечения ускорительной камеры, т. к.
в0 = в -{- (0, 5 1, 0), см,
тде в размер радиальной полуоси поперечного сечения ускорительной камеры (рис. 1.).
Из уравнения (6)
О "О
(7)
Подставив в уравнение (1) значение г0 из уравнения (7) и В0 из уравнения (5), получим
А
1,5.10Все0
2-я
1 - Т
п —2
или
А
[,5-Ю-"4 Всв()
(В)
где
ТгО ' )1
IV К-
(9)
Зависимость ¡3 от -(■ для различных значений п (по уравнению 9) изображена в виде кривых (рис. 2).
Р
1.0
0,5
Уравнение (8), как модификация уравнения (5), связывает основные размеры и параметры междуполюсного пространства бетатрона при условии соблюдения бетатронного отношения.
Постоянная бетатрона
В уравнение (8), если учесть и уравнение (9), входят три независимых, неизвестных величины: Вс, п и ви.
При питании обмотки электромагнита бетатрона переменным синусоидальным током мгновенное значение плотности магнитного потока на орбите в конце процесса ускорения
70 Вор — Воа.ъта.^
где Воа -амплитуда плотности магнитного потока на орбите а — 65-=- 75е — рекомендуемая фаза сброса электронов на мишень.
Следовательно, максимальная расчетная энергия ускоренных в бетатроне электронов согласно уравнению (8) равна
л_ \,5A0~*.BcMKeeosinap
где ВсМ максимальная допустимая плотность магнитного потока в стали центральных вставок, кс~0,83 -г-0,87 коэффициент заполнения центральных вставок сталью.
При выборе значения Всм следует учитывать кривую намагничивания, удельные потери стали магнитопровода и способ охлаждения электромагнита.
В бетатронах с максимальной энергией от 5 до 25 Мэв конструкции Томского политехнического института при работе на переменном токе промышленной частоты (/"—50 гц) с естественным охлаждением электромагнита
Всм— 13500-:- 14000 гс •
Для изготовления магнитопровода бетатронов использовалась стальЭ- 42 толщиной 0,35 мм.
Расчетное значение показателя спадания поля п обычно принимается равным п — 0,5-ь- 0,75.
Остается неизвестным значение ви. Как указывалось ранее, ва зависит от размеров поперечного сечения ускорительной камеры бетатрона.
Размеры поперечного сечения ускорительной камеры зависят от многих факторов. В принципе, чем больше размеры поперечного сечения камеры, тем вероятнее захват в ускорение большего числа электронов, испускаемых электронной пушкой в момент подачи импульса напряжения инжекцни, а следовательно, получение более высокой интенсивности гамма-излучения. Но с увеличением размеров поперечного сечения камеры увеличиваются размеры электромагнита и мощность батареи конденсаторов для компенсации реактивной мощности.
На размеры поперечного сечения камеры в сильной степени влияет конструкция магнитопровода, от которой зависит фазовая структура магнитного поля в рабочей зоне междуполюсного пространства.
В бетатронах с симметричным магнитопроводом поперечное сечение камеры относительно меньше, чем при Ш-образной конструкции магнитопровода.
Вопрос выбора размеров поперечного сечения ускорительной камеры здесь не рассматривается и при расчете радиуса равновесной орбиты бетатрона, размер в0 должен быть задан.
Из уравнения (10) можем написать следующую зависимость:
_в0 __ 6,7-Ю-__
АЬ ВсМ.кс.$\
где А выражено в Мзв.
Правая часть этого уравнения при выбранных значениях Всмл, и кс представляет постоянную величину. Обозначив эту величину через Б, можем написать
ЛГ1 i
6,7.103
где Ь ---- — -------
BCM.Kc.smap
может быть названа „постоянной бетатрона".
(П) (12)
Определение радиуса г„
11ри шдбраниых значениях Б, п и в0 из уравнения (И)
8 — •
Определив по кривым рис. 2 соответствующее значение найдем из уравнения (7)
Го~ 6 О ■ ,- •
V .... 1
Как видно из рис. 2, кривые ¡3=/(г) для различных показателей степени спадания поля п при значениях у < 1,5 совпадают, стремясь к прямой линии, соответствующей п — 1. Следовательно, при значениях у <1,5 без ущерба для точности расчета можно пренебречь влиянием п на величину расчетного радиуса равновесной орбиты. Принимая п = 1, из уравнения (9) получим:
Р = Т~1. (13)
Подставив полученное значение ¡3 в уравнение (11) и учитывая равенство (6), на идем
= АБ + в0\. (14)
Уравнение (14) дает вполне достаточную для технических расчетов точность. если -(<1,5, что соответствует бетатронам с максимальной энергией электронов порядка 10 Мэв и выше.
Обычные конструкции магнитопровода бетатрона позволяют осуществлять, в небольших пределах юстировку радиуса равновесной орбиты, если при испытании электромагнита обнаружено отклонение радиуса орбиты от расчетного.
Учитывая возможность юстировки, уравнение (14) можно использовать и для определения радиуса равновесной орбиты в бетатронах с максимальной энергией ускоренных электронов меньше 10 Мэв.
Если максимальная заданная кинетическая энергия ускоренных в бетатроне электронов V? много больше собственной энергии (г0 —0,51 Мэв), то вместо уравнения (14) можно пользоваться уравнением
г0^ЪБ + в0. (14й)
Уравнения (14) и (14д) позволяют быстро и для заданных условий в известной мере однозначно решить вопрос о выборе радиуса равновесной орбиты бетатрона.
Расчет радиуса г0 в бетатронах с кольцевыми пазами
на полюсах
Значительная часть энергии магнитного поля бетатрона, а следовательно, и реактивной мощности, потребляемой намагничивающей обмоткой электромагнита, приходится на долю воздушного зазора центральной части междуполюсного пространства в круге с радиусом гс, где плотность магнитного потока имеет наибольшее значение.
Емкость батареи конденсаторов, необходимых для компенсации реактивной мощности электромагнита бетатрона, можно значительно снизить, если уменьшить или совсем устранить всздушный зазор в центральной части междуполюсного пространства.
В этом случае необходима дополнительная обмотка, компенсирующая действие главной намагничивающей обмотки в центре (в круге с радиусом гс).
Когда применяется компенсационная обмотка, полюсные сердечники должны иметь кольцевой паз для размещения обмотки (рис. 3).
Кольцевые пазы на полюсах необходимы также для размещения компенсационной обмотки при подмагничивании электромагнита бетатрона постоянным током.
Кривая зависимости B—j(r) в электромагните бетатрона с кольцевыми пазами на полюсах приведена на рис. 3.
этом случае для определения радиуса равновесной орбиты получаются уравнения, аналогичные приведенным выше, но вводится новая величина вк- - ширина кольцевого паза, которая может быть вычислена предварительно (стр .74).
ЯР-МО
Рис. 3
Для бетатронов с кольцевыми пазами на полюсах уравнение (14) принимает вид
f
го = (*о + *к) + 0,5АБ+у (0ч5АБ)* + вкАБ .
(15)
При точных расчетах этим уравнением и следует пользоваться. В приближенных расчетах, особенно для выявления закономерностей при использовании метода подобия, уравнение (15) неудобно ввиду его сложности.
Разложим в ряд по формуле бинома Ньютона последний член правой части уравнения (15):
^0,5 АБ^
2 i вкАБ
°-5 = 0,5АБ + вк-
АБ
подставив первые два члена ряда в уравнение (1о), получим
\г0^АБ + {е0 + 2вк) .
(15«)
73
Предварительное определение ширины кольцевого паза
полюса
При расчете радиуса равновесной орбиты в электромагните бетатрона с кольцевыми пазами на полюсах ширина кольцевых пазов должна быть известной. Но в то же время окончательные размеры пазов могут быть определены только после расчета обмотки, которая размещается в кольцевых пазах. В этом случае размеры паза необходимо с достаточной точностью предугадать.
Для предварительного выбора ширины кольцевого паза вк (рис. 3) можно руководствоваться следующими соображениями.
При отсутствии зазора в средней части междуполюсного пространства бетатрона амплитуду ампервитков компенсирующей обмотки в приближенных расчетах можно считать равной магнитному напряжению зазора В0 при Воа< т. е.
А^к = 0,8ооВ„а.
Если зазор сохраняется, но имеет размер меньше против расчетного без компенсирующей обмотки, то
А — 0,8с,- В()аЬ0,
о — Ьк
где ск — —---------степень компенсации,
ое — расчетный зазор в центре, без компенсирующей обмотки, зазор при наличии компенсирующей обмотки.
Площадь поперечного сечения кольцевого паза, необходимая для размещения одной катушки компенсирующей обмотки
Л^ = °'283 с*'Во«-?'» (Ш)
|/ •_) ^ Л'
где '<з ~~ плотность тока в пазу,
]к — плотность тока в меди компенсирующей обмотки, к3 коэффициент заполнения паза медью, зависящий от изоляции паза и обмотки и от системы охлаждения.
Примерные значения А1С в зависимости от напряжения на зажимах обмотки при^ непрерывном питании обмотки электромагнита переменным током указаны в таблице 2. При импульсном питании обмотки допустимая плотность тока, в соответствии со скважностью импульсов, увеличивается.
Т аблица 2
у к кв
~л,г
3 6
6 12
50 — 40 . 40 35
12 24
35 30
Для обмотки с проводом прямоугольного сечения в небольших бетатронах при отсутствии воздушных каналов в катушках обмотки плотность тока в пазу может иметь значения порядка 100-г-120 а/см2.
Воздушный зазор на радиусе равновесной орбиты зависит от размера вп. В приближенном расчете можно принять
^-Ч, I 1 л -* п
Если компенсирующая обмотка размещается только в кольцевых пазах, то площадь, занимаемая поперечным сечением одной катушки при к, ~т$к, равна 8к—твк2 и ширина кольцевого паза
- 07)
V т
где АЛ. - высота (глубина) кольцевого паза.
При выборе коэффициента т следует учитывать, что чем больше тем меньше расчетный радиус равновесной орбиты, т. е. меньше радиальные размеры полюсов. В то же время не может быть больше высоты полюса (1гп). В зависимости от конструкции магнитопровода /те—1.5--:-3.
В случае, когда компенсирующая обмотка располагается по всей высоте центрального сердечника, площадь поперечного сечения, занимаемая всег компенсирующей обмоткой, будет равна
откуда ширина кольцевого паза
ви
(ib)
4 т
т
В уравнение (16) входит неизвестная пока величина Вшг амплитуда плотности магнитного потока на равновесной орбите. В зависимости от размеров поперечного сечения ускорительной камеры, от степени компенсации и от постоянной бетатрона Воа может иметь значения в пределах Воа = 2200н-4500 гс. Меньший предел относится к небольшим бетатронам, работающим на переменном токе повышенной частоты, высший соответствует бетатронам с максимальной энергией ускоренных электронов л. 100 Мэв и выше.
& 0« ¿Ojee
Чтобы определить площадь поперечного сечения кольцевого паза полюса., необходимо предварительно задаваться значением Воа и, определив ширин1, паза вк, а затем г0 и соответствующее ему значение Воа, сравнить с предварительно принятым. Если вычисленное значение В,м отличается от пред-
дарительного более, чем на 5°ч>, следует повторить расчет Sh-, вк и г() при ;ыпи /i' ином значении Воп.
Для ориентировки при выборе предварительного значения Воа можно пользоваться кривыми (рис. 4), соответствующими постоянной бетатрона В = '\в. Если постоянная проектируемого бетатрона (Б) не равна 0,6, то
Г» -л 0'6
значение Ьоа, наиденное по кривон рис. 4, следует умножить на----•
Б
Описанный метод дает приближенное, ориентировочное значение вк. После расчета обмотки размеры паза иногда приходится корректировать, что лучше всего осуществлять за счет изменения глубины паза hK. В этом случае никаких изменений в других размерах магнитопровода не потребуйся.
