О влиянии геометрических параметров системы электродинамического подвеса на силовые характеристики

Размещено на http://www.allbest.ru/

О влиянии геометрических параметров системы электродинамического подвеса на силовые характеристики

Кочубей Т.В., Майборода А.О.

Введение

В настоящее время происходит активное освоение космоса и актуальной задачей развития космонавтики является радикальное снижение затрат на доставку грузов в космос. На данный момент доставка осуществляется при использовании космических летательных аппаратов (КЛА). Однако, двигательные установки КЛА должны иметь высокий уровень удельной мощности, работая при этом на пределе прочности. Поэтому стоимость единицы массы конструкции получается чрезвычайно высокой.

Одним из направлений удешевления КЛА является сокращение нынешней высокой удельной мощности двигательных установок (ДУ) за счет отказа от технического совершенствования конструкции, в том числе замены турбонасосной системы подачи топлива в двигатели вытеснительной системой. В этом случае стоимость двигательной установки снижается в 50-100 раз. Однако, такой вариант удешевления КЛА и повышения ресурса их работы может применяться только в межорбитальных операциях, но не при выводе КЛА с Земли в космос. Это связано с дополнительными гравитационными потерями ракеты при использовании тяжелых термохимических двигателей и топливных баков с вытеснительной подачей топлива в процессе разгона до орбитальной скорости в условиях старта в земном гравитационном поле. Данную проблему можно решить, если поместить КЛА (на время разгона) на специальную опорную поверхность, которая проходит вне плотных слоев атмосферы и в несколько раз увеличить время разгона с пропорциональным сокращением тяги двигателя и, соответственно, массы. Подобной заатмосферной опорой могут быть различные варианты орбитальных тросовых систем, которые ориентированы вдоль траекторий своего движения и движения разгоняемых на них КЛА. В общем виде они представляют собой орбитальную платформу с электроракетными двигателями (ЭРД), которая тянет за собой монорельсовую опорную поверхность. Тяга ЭРД уравновешивает силу аэродинамического торможения троса-монорельса в верхних слоях атмосферы и при условии превышения сил торможения над локальными силами микрогравитации, действующими на трос, движение монорельса стабилизируется. Стабилизации движения троса вдоль самого себя по низкой круговой орбите тем эффективнее, чем больше отношение силы торможения к силе микрогравитации. Тросовый монорельс может быть создан непосредственно на орбите, путем непрерывной формовки из запасов алюминиевого сырья, недорогая доставка которого на орбиту теперь возможна благодаря системе космических аппаратов-накопителей сырья (патент RU 2398717).

Такие тросовые системы представляют собой гибкие стержни, в форме тонкостенных труб, движущиеся вдоль самих себя по круговой орбите, и носят название орбитальной путевой структуры (ОПС). Необходимо учесть, что относительная скорость КЛА и ОПС варьируется от первоначальной 8000 м/с до финишных 800-500 м/с и при этом должен быть исключен механический контакт между КЛА и ОПС. Данное условие может быть реализовано на основе магнитного подвеса КЛА над или под ОПС. Перспективным видом подвеса является электродинамический подвес на основе постоянных магнитов КЛА, движущихся над направляющим проводником ОПС.

Принцип работы электродинамического подвеса основан на левитации (свободном парении) твердых тел в магнитном поле. Когда тело, несущее ток, перемещается над направляющим проводником, то в проводнике появляются вихревые токи с собственным электромагнитным полем. Это поле, взаимодействуя с первичным источником, создает подъемную силу, направленную вверх, и тормозящую силу, используемую для разгона КЛА. Чтобы конструкция подвеса отвечала поставленным задачам и имела оправданную стоимость, важно знать зависимости упомянутых сил от геометрических параметров подвеса и расположения магнита относительно проводника. Для этого необходимо рассчитать упомянутые выше вихревые токи.

В данной работе предлагается модель расчета для вихревых токов, возбуждаемых в направляющем проводнике электродинамического подвеса, в виде интегро-дифференциального уравнения второго рода. Для решения уравнения используется метод преобразований Фурье. С помощью полученного решения выполняется анализ влияния геометрических параметров и расположения магнитов относительно направляющего проводника на силовые характеристики электродинамического подвеса.

1. Постановка задачи

Рассмотрим электродинамический подвес, состоящий из направляющего проводника и движущегося вдоль него несущего магнита. Длина проводника много больше его диаметра, поэтому можно пренебречь концевым эффектом и рассматривать проводник как бесконечно длинную цилиндрическую оболочку, выполненную из однородного материала с проводимостью . Поперечное сечение проводника представлено на рис. 1. На рисунке отмечены направления компонент вектора силы, действующей на несущий магнит со стороны проводника, а также гравитационной силы , действующей на груз. Магнит движется вдоль проводника с постоянной скоростью .

Рис. 1 Поперечное сечение направляющего проводника (1) и несущего магнита (2)

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

Для расчета магнитного поля и силовых характеристик подвеса необходимо выполнить расчет вихревых токов направляющего проводника. Сделаем еще одно допущение: будем полагать проводник тонким (толщина проводника много меньше его диаметра и электрического клиренса ) и перейдем к идеализированной (бесконечно тонкой) оболочке. Причем под проводимостью и линейной плотностью поверхностных вихревых токов можно понимать произведение соответствующих величин в объеме на толщину рассматриваемого проводника, если поверхностный эффект мал. Способ учета этого эффекта будет указан ниже.

Запишем уравнения Максвелла, материальные уравнения и граничные условия для системы подвеса с учетом принятых идеализаций:

, , ,

,

на ,

,

на ,

, на ,

где знаками «+» и «-» отмечены предельные значения поля со стороны положительного и отрицательного направлений нормали , соответственно; - срединная поверхность проводника; , . Стоит заметить, что данную задачу удобно рассматривать в цилиндрической системе координат .

Индукцию поля можно представить в виде:

,

где - индукция магнитного поля, создаваемого магнитами, расчет которого будет описан ниже, а индукция поля магнитной реакции поверхности по закону Био-Саварра-Лапласа [1] выражена через плотность вихревых токов .

Введем скалярную функцию тока , связанную с плотностью следующим равенством:

.

Такое представление возможно в силу соленоидальности поля . Тогда можно преобразовать к виду

.

Применяя операцию к закону Ома из векторной задачи и используя , получим интегро-дифференциальное уравнение относительно функции потока

на ,

где - поверхностный лапласиан функции .

2. Численное решение уравнения

Решить уравнение можно несколькими методами, одним из которых является метод интегрального преобразования Фурье [2]. Это преобразование раскладывает искомую функцию в интегральную сумму тригонометрических гармоник, круговая частота которых принимает непрерывные значения от до и называется параметром преобразования.

Формально преобразование Фурье задается следующими правилами:

и

,

устанавливающими соответствия между оригиналом и его Фурье образом . Это прямое и обратное преобразования Фурье соответственно. Здесь - параметр преобразования.

С помощью первого правила задача для искомой функции преобразуется в задачу для ее Фурье-образа. Поскольку последняя обычно оказывается проще, ее решают и находят Фурье-образ ответа. Остается применить к нему обратное преобразование и получить оригинал. Однако, в некоторых случаях необходимость в обратном преобразовании отсутствует. Например, если искомая функция необходима только для вычисления интегральных характеристик. В этом случае удобно использовать равенство Планшереля [2]

.

Знак «» в формуле означает комплексное сопряжение.

Выполним преобразование Фурье по координате с параметром преобразования , а также разложим функцию в тригонометрический ряд Фурье

.

После этого уравнение можно записать в виде

на ,

Где ;

- функция Макдональда, получившаяся в результате использования табличного интеграла [3].

Применив теорему сложения для функции Макдональда [4], получим

,

где «?» обозначает производную по всему аргументу.

Раннее отмечалось, что под можно понимать произведение на толщину проводника. Однако, в этом случае мы не учитываем явление поверхностного эффекта. Для его учета [5] под следует понимать комплексную величину

,

где , - толщина проводника.

С учетом всех преобразований, можно записать аналитическое решение уравнения . А именно,

.

Остается нерешенным вопрос о способах задания и расчета поля, создаваемого несущим магнитом. Рассмотрим его далее.

проводник силовой электродинамический подвес

3. Расчет поля несущего магнита

На практике магниты изготавливают из магнитотвердых материалов, например из сплава самария и кобальта, которые в условиях эксплуатации практически сохраняют неизменным вектор намагниченности. Это делает возможным их эквивалентную замену токовым слоем с плотностью , связанной с вектором намагниченности (рис. 2) следующим образом:

.

Здесь - внешняя нормаль к поверхности магнита.

Рис. 2 Геометрия и расположение несущего магнита относительно направляющего проводника

Рассмотрим прямоугольный виток с током (рис. 3). Индукция магнитного поля такого витка вычисляется по формуле

Рис. 3 Геометрия и расположение витка с током относительно направляющего проводника

Для расчета требуется нормальная компонента индукции, которая может быть вычислена следующим образом

.

Применяя к обеим частям равенства (3.1) преобразование Фурье по координате , а затем теорему сложения для функции Макдональда, получим

При рассмотрении пары витков с током, расположенных симметрично относительно оси , возможны две ситуации, представленные на рис. 4.

а) б)

Рис. 4 Расположение и варианты полярности витков с током а) одинаковая полярность, б) чередующаяся полярность

На рис. 4а токи в витках однонаправлены, а на рис. 4б - имеют противоположные полярности. Формулы для расчета Фурье-образа в этих случаях имеют вид

а) ,

б) .

Для учета толщины магнита, чтобы получить индукции , и , остается проинтегрировать по толщине магнита величины , и соответственно.

4. Расчет силовых характеристик электродинамического подвеса

При конструировании электродинамического подвеса важно знать зависимости его силовых характеристик от различных величин (радиуса и толщины проводника, высоты подвеса, расположения и массы магнита). Исследование этих зависимостей позволяет минимизировать затраты на его изготовление, например, за счет снижения расхода материала. При этом подвес будет обладать необходимыми силовыми характеристиками.

Формула для вычисления электромагнитной силы, действующей на несущие магниты подвеса со стороны направляющего проводника, имеет вид

С учетом представления , тождеств векторного анализа [6] и соленоидальности магнитной индукции получим, что

.

Используя равенство Планшереля , окончательно получим

где

Сделаем важное замечание, необходимое для дальнейших расчетов. А именно, обратим внимание, что вещественная часть Фурье-образа всякого вещественного оригинала обладает четной симметрией относительно начала оси параметра и имеет своим оригиналом функцию, обладающую четной симметрией относительно начала оси координат . Аналогично, мнимая часть обладает нечетной симметрией, а ее оригинал сохраняет эту симметрию на оси координат . Данное свойство является очень важным при численных расчетах, так как позволяет упростить и существенно сократить объем вычислений. Из него следует, что интеграл в симметричных пределах от произведения вещественной части одного из сомножителей в на мнимую часть другого сомножителя равен нулю, а интеграл в тех же пределах от одноименных частей сомножителей, благодаря четной симметрии подынтегральной функции, равен удвоенному интегралу на половине интервала интегрирования. Таким образом, компоненты электромагнитной силы, действующей на магниты, могут быть вычислены по формулам:

а) ,

,

;

б) ,

,

Здесь при вычислении в формуле в качестве необходимо подставить значение или для случая а) и б) соответственно.

В качестве силовых характеристик электродинамического подвеса можно использовать зависимости сил левитации и торможения , действующие на несущие магниты, от геометрических и физических параметров. Значения этих сил выражаются следующим образом

,

Рассмотрим электродинамический подвес, удовлетворяющий нашим идеализациям. Направляющий проводник подвеса выполнен из материала с удельной проводимостью . Радиус проводника м. Конструкция имеет одну пару магнитов, соответствующих случаю а), со следующими геометрическими параметрами: толщина магнита - м, м. Магниты движутся вдоль проводника со скоростью , Тл. На рис. 5 и 6 изображены результаты расчетов сил левитации и торможения в зависимости от высоты подвеса и толщины проводника соответственно при фиксированных значениях скорости. Здесь и далее графики 1, 2, 3, 4, 5 и 6 рассчитаны для значений скорости равных 1, 4, 7, 10, 13 и 16 км/с соответственно.

Рис. 5 Зависимость сил левитации и торможения от высоты подвеса (толщина проводника м)

Выполнив анализ этих графиков можно найти оптимальную конфигурацию электродинамического подвеса, отвечающего необходимым условиям. Например, из рис. 5 видно, что увеличение высоты подвеса приводит к уменьшению как силы торможения, так и силы левитации.

Рис. 6 Зависимость сил левитации и торможения от толщины проводника (высота подвеса м)

Еще одной характеристикой системы является отношение , где - подъемная сила, а - сила тяжести, действующая на магнит. Ниже на рис. 7 приведена зависимость этого отношения от скорости при толщине проводника м и высоте подвеса м, а на рис. 8 - зависимость от массы магнита, которая изменяется путем увеличения его толщины, при тех же условиях.

Рис. 7 Зависимость отношения от скорости

Рис. 8 Зависимость отношения от массы магнита

Из рис. 8 видно, что, начиная с определенного момента, увеличение массы магнита приводит к ухудшению качественных характеристик электродинамического подвеса. Это особенно заметно для малых скоростей (1-4 км/с).

Таким образом, учет свойств электродинамического подвеса и исследование его силовых характеристик позволяют создавать конструкции с параметрами, необходимыми для выполнения поставленных задач и имеющие при этом оптимальную стоимость.

Литература

1. Шимони К. Теоретическая электротехника /Пер. с нем. М.: Мир, 1964. 773 с.

2. Титчмарш Е. Введение в теорию интегралов Фурье / Пер. с англ. М.: ОГИЗ

3. Двайт Г.Б. Таблицы интегралов и другие математические формулы. / М.: Наука, 1973. 228 с.

4. B.G. Korenev Bessel Functions and their Applications, Taylor&Francis, London, 2002. 276 p.

5. Астахов В.И. Математическое моделирование инженерных задач в электротехнике / Учебное пособие. Новочеркасск: НГТУ, 1994. 192 с.

6. Кочин Н.Е. Векторный анализ и начала тензорного исчисления. М.: Изд. АН СССР, 1961. 425 с.

Размещено на Allbest.ru