Математическая модель и анализ импульсного режима работы линейного электродвигателя при независимом питании его обмоток Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

УДК 621.313

В.И. Милых, С.В.Ткаченко

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

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

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

Введение. Линейный импульсный электродвигатель (ЛИЭД) рассматриваемого типа [1, 2] предназначен для создания в земной коре слабых сейсмических колебаний для поиска полезных ископаемых. Его конструкция и принцип действия представлены в [1-3], а здесь его упрощенная схема показана на рис. 1.

При подаче импульсов токов іг и іа в обмотки взаимодействие реактора и якоря порождает электромагнитную силу Егт. Якорь движется со скоростью уа и через опорную плиту производит ударное воздействие на грунт, а инерционный реактор со скоростью уг отскакивает вверх. За рабочий интервал взаимное перемещение якоря и реактора составляет примерно половину полюсного деления (полюсное деление тр - шаг между осями соседних зубцов якоря или реактора).

Рис. 1. Схема активной силовой части ЛИЭД: 1, 2, 3 - обмотка, сердечник и инерционная масса реактора; 4, 5, 6 -обмотка, сердечник и опорная плита якоря; 7 - грунт

Математическая модель и расчетный анализ ударного режима работы ЛИЭД были представлены в [2]. Там рассматривался вариант его питания от одного источника при последовательном соединении обмоток реактора и якоря. Однако с практической точки зрения может быть целесообразным и вариант с их независимым питанием от индивидуальных источников. Это позволяет понизить уровень напряжения ЛИЭД и всего электрооборудования в два раза и, тем самым, повысить электробезопасность в условиях эксплуатации и упростить систему электроизоляции.

Электрическая схема замещения ЛИЭД при независимом питании обмоток представлена на рис. 2,

где индексы г и а в обозначениях параметров соответствуют реактору и якорю.

Рис. 2. Электрическая схема замещения ЛИЭД

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

В схеме на рис. 2 представлены также сопротивления Ксг, Яса с температурной нелинейностью и индуктивности Ьсг, Ьса подводящих электрокабелей.

Конденсаторные накопители энергии имеют емкости СЕг, СЕа и изменяющиеся напряжения

І , І

иСг - иСго ~ С 1іг^І ; иСа

— }Іа^І , (1)

. г“* ’ “Са иСао | ‘а

СЕг о СЕа о

где иСго, иСао - их начальные значения; t - время.

Инициализация работы ЛИЭД производится одновременной коммутацией тиристоров УБ1 и ¥Б2 в цепях обмоток реактора и якоря, они же исключают изменение направления токов в процессе работы.

Электромагнитные процессы в ЛИЭД неразрывно связаны с силовыми взаимодействиями.

Источником силовых процессов является упоминавшаяся сила Еет, которая в равной мере действует и на реактор, и на якорь, только в разные стороны (рис. 1). Эта сила также является трехпараметрической функцией, завися от токов 1Г, 1а и взаимного аксиального смещения реактора и якоря С,.

На реактор и якорь действуют также силы инерции Етг = таг, Ета = тааа и их весовые силы Р, = т^^, Ра = т^,., где тг, та, аг, аа - массы и ускорение движения реактора и якоря; ggr - гравитационная постоянная.

К важнейшим силам относится сила реакции упругого грунта Е^, которая в условиях эксплуатации зависит от целого ряда слабо прогнозируемых факторов, поэтому ее адекватное определение в процессе расчета является практически невыполнимой задачей. Но для обеспечения расчетов ЛИЭД в [2] была принята гипотетическая зависимость силы Е^(С,а) от осевого смещения якоря. В исходном положении = 0, а при

ударном режиме работы имеет место максимальная деформация грунта на С,ат. При этом Е^ достигает максимального значения Е^т, которое задается на основе тестовых расчетов ЛИЭД в зависимости от максимально возможной силы Еет.

Для прямого хода якоря с углублением опорной плиты в грунт была принята линейная зависимость Е^ = £,-ыСа, гДе коэффициент упругости грунта кге!!\ = Е^/^ат. После ослабления силы Еет упругая сила Е^ отбрасывает якорь в исходное положение и на обратном ходе принята параболическая функция Е^ = £ге.д(Са _ Сао)2 с частично необратимой деформацией грунта С^ао, где коэффициент упругости грунта кгеЛ =

Еgrm / (Сат Сао) .

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

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

Математическая модель ЛИЭД.

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

№

ж

Са ж

(2)

(3)

где Тг, - результирующие магнитные потокосцеп-

ления (МПС) обмоток реактора и якоря, складывающееся из собственных и взаимных МПС; иЯ, = 1Г(ЯГ + Ясг), иЯа = 1а(Яа + Яса) - падения напряжения на соответствующих сопротивлениях.

С учетом трехпараметрических зависимостей МПС их производные преобразуется к виду

(№г СТ. + Жа • (4)

Ж Л д1г Ж д1а Ж

Жа ^а + ^ . (5)

Ж Л д1а Ж д1г Ж

В этих выражениях выделим дифференциальные собственные и взаимные индуктивности

г Я* д¥а „ д% „ д^а

Ь, =—- ; Ьа =——; Мга =—- ; Маг =-

(6)

' и ' ги ' иг

01, 01а 01а 01,

а также коэффициенты, которые для удобства опери рования назовем условно индуктивностями движения

(7)

с1^ ъ Сс

Обозначим временные производные токов через

I,' = &,/&, 1а' = &а/А, а также скорость взаимного перемещения реактора и якоря уга = сСуА, которая составляется из их индивидуальных скоростей уга = V, + уа.

Теперь можно оперировать с ЭДС самоиндукции и взаимоиндукции

еь, = Ь,1, ; еЬа = Ьа1а ; еМт = Мт‘а ; еМаг = Ма,‘, , (8) а также с ЭДС движения

evr ~ ЬrCVra ; eva ~ ЬаСУ,а . (9)

В итоге уравнения (2), (3) представляются в виде

(Ь, + Ьс, ) • I, + М,а 'га = иС^ ~ иЯ^ ~ evr ; (10)

Ма, '+ (Ьа + Ьса ) ' 1а = иСа ~ иЯа ~ ^а . (11)

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

Т = Т +-

£е £ео 1

(12)

сретСие о

Яе = Яео[1 + аре (Те - Тео)] , (13)

где индекс е подразумевает соответствующие параметры проводника каждого из названных элементов; тсие, сре, Тео, Яео, Орс - их массы, удельная теплоемкость, исходные температуры и соответствующие сопротивления, температурный коэффициент сопротивления.

Баланс сил в механической системе ЛИЭД (рис. 1) описывается на основе уравнений

Е = Р + Е ; Е = Е + Е - Р (14)

1 ет 1 , ~1 т, ? ет х та 1 х gr х а ?

из которых получаются ускорения и скорости движения реактора и якоря и их аксиальные смещения:

а, ~ Еет / т, ~ ggr ; аа ~ (Еет ~ Еgr ) / та ggr , (15)

г г

V, =| а,Сг; va =| ааСг , (16)

о о

г г

(17) о о

что и дает их взаимное смещение

С=С, +Са +Со , (18)

С , = | vrdt ; С а = | vadt,

где - начальное смещение реактора в сторону, противоположную его движению из положения "нейтрали", когда зубцы и пазы реактора и якоря расположены посредине друг друга, т.е. < о [1, 2].

Принципы подготовки и работа с трехпараметрическими функциями

Для решении системы всех представленных уравнений и соотношений необходимы трехпараметрические функции МПС ^(С, /,, /а), %(С, /,, /а) и электромагнитной силы Еет(С1,/8, /а). Их подготовка, а также оперирование с частными производными МПС, оказались наиболее сложной частью расчета переходного процесса, характеризующего импульсный рабочий режим ЛИЭД. Принципы преодоления этой проблемы заложены в [2, 4], но теперь задача еще многократно усложнилась ввиду перехода от двух- к трехпараметрическим функциям.

Как И В [2, 4], функции ¥Г(С, /,, /а), %(С, /,, /а) и Еет(С>, /& /а) были подготовлены численно-полевыми расчетами методом конечных элементов по программе ГБММ [5] с учетом реальной геометрии, "коллекторной" шихтовки и насыщения сердечников. Было определено, что для получения достаточно "гладких" функций должно хватить структуры сочетаний из Кг = 1о значений С, и по К= К/а = 12 значений /, и /а, т.е. общее число сочетаний составило (Юх 12x12) или 144о, на что сравнительно эффективному компьютеру (2,8 ГГц) потребовалось около 22 часов расчетного времени. Естественно, что "ручная" организация перебора вариантов и "извлечения" результатов расчета были неприемлемы, а с задачей справилась управляющая программа, написанная на алгоритмическом языке Ьиа [5].

Для получения в ходе расчета переходного процесса конкретных значений Еет и дифференциальных индуктивных параметров Ь,, Ьа, Мт, М^, Ь^ и Ь^ при текущих значениях С, /, и /а использовалась локальная аппроксимация табличных функции кубическими полиномами отдельно по каждому из аргументов при фиксировании двух других аргументов, подобно тому, как это делалось для двухпараметрической функции [4]. Рассмотрим это на примере более трудного варианта -МПС, так как для него надо определять не сами значения, а частные производные по трем параметрам (4)-(7).

Как уже представлено, имеется трехмерный массив МПС (индексы , или а при Т уберем, т.к. для МПС реактора и якоря действия, в принципе, аналогичны):

%,и , п = 1.К, / = 1..К/,; / = 1 ..К/а , (19)

которому соответствовали одномерные массивы перемещений и токов

Сп , п = ; /, , ] = 1..-К/, ; /а , / = 1..К/а . (2о)

Для иллюстраций структуру полученного массива (19) представим фрагментарно на рис. 3 в виде трехкоординатной сеточной структуры, где координаты линий соответствуют (2о), а узлы - (19).

Предположим, что в процессе расчета определены значения С,, /, и /а, а для продолжения расчета необходимы частные производные (6), (7). Для этого предлагается следующая процедура, основанная на аппроксимации функции Т(^, /,, /а) однокоординат-

ными кубическими полиномами. На рис. 3 точка с указанными координатами помечена кружочком и попала внутрь ячейки в виде параллелепипеда, левый нижний ближний угол которого идентифицируется тройным номером П,І, І.

І I _ І- I „

\-¥------Т\Г*-------Ї 7*------Г^--

_ І “I--Г7+-1-------Т ^ - 1------Г-*-Г

/ + 2 і “і—1~г “і—і— _г“гг ~г і

I ” ~Г -Г"1"-?Т*ГТ-р^-ÓÓ1----------

Г'У+2 і і і і і і і і і і і і і і і і

і і -Аґі. _ і_ _і_ -Ъ?-± -1 -1 _±-і -1 -1-

/ +11 -+Ґ-І- -і_|_ .

II

І I 1 1 1 1 1 1- -# • _і_ і_

і I _

І І I 1^1 I _ |_ Іґ_ І І I _ І І _1Ґ П+ - *

1 1 І-ТГ Г і„Т 1~1 І^тт Т ЬТн + ГгЬя+2

1/Т і 1-ГГ і 1" Т і ь-ТП ГЬл+1

_ _1 __I _ .1_

III

-1,01 і і ^

__і I 1~Т7~Т

I 1 1^Т _ і_ 1_ I

_ 1_ _1_ 1 _ 1111 _ 1 _ 1_________

*Т і I гт і і і

*а,1+1 1-1 1 ' + -

Рис. 3. Условная структура массивов для определения производных МПС и электромагнитных сил

Сначала из структуры на рис. 3 возьмем одну плоскость с фиксированным номером п (рис. 4,а). На горизонтальной линии с номером / в пределах отрезка между номерами / и /+1 предполагаем кубический полином

^ = а + Ь • 1а + с • + С • . (21)

Значения коэффициентов а, Ь, с, С получим, используя четыре значения МПС из массива (19) в прилегающих узлах на пересечении линии / с линиями /-1, /, /+1 , /+2 и на основе (21) составим систему алгебраических уравнений

а + Ь ■ І

а,і-1

+ С • І

а,і-1

+ ^ * І(Я,і-1 %,І,і-1;

а + Ь ' іа,і + С' ^ + ^ ' іа,і = і і ;

а + Ь 'іа,і+1 + С' і'а,і+1 + й ' іа,і+1 = ^п,І,і+1; а + Ь ' іа,і+2 + С' іа,і+2 + й ' іа,і+2 = ^п, /,і+2.

(22)

. 9і

J + =*.----■*---¥—*—

І І

І

.4!..

/і І І

-------4

і

/-11 1 і

V

+Уг

г,п 2

І

*г,у+1

/ \ \ / п

/1

. ^1! .

І -І- ^ -■О—!— ! /Л, ,

" !/ 0 ■ *' і+и / . ■* ---♦-----'

V-1

І

; + 2

, л 1 . , ^ . І / І

'і ' ж

І Ки+1

и + 2

і

/ — 1 г 1 + 1

7-1,*и-1; а б

Рис. 4. Частн^1е одноплоскостн^1е структуры элементов массивов

Решив систему алгебраических уравнений (22) по (21) находим значение МПС в точке при значении /а (крестик) и здесь же производную по /а

(23)

%, і,іа а ^ Ь * іа С * іа й * іа ;

(24)

(^п,],/а )/а Ь 2 'с ' ^а ^ 3 'С '/а .

Проделав все то же на линиях _/-1, у+1, _/'+2, получим значения МПС и его производной в наборе точек, отмеченных на рис. 4,а крестиками, а именно,

%,у-1,10 ; %,у,1а ; %,у +1,1а ; %,у+2,10 ; (25)

(^п,у—1,1а ) 1а ; (^п,],1а )/а ; (^п,у+1,1а )га ; (^п,у+2,Ш )га (26)

Теперь на вертикальной линии (рис. 4,а) с точками, отмеченными крестиками, в пределах отрезка между номерами у и у+1 со значением /V предполагается свой полином со своими коэффициентами а, Ь, с, й

Т = а + Ь • 1, + с • 1, + й • 1, . (27)

Для четырех указанных точек составляется система алгебраических уравнений

а + Ь ■ і а + Ь • і

г, У-1 + ,г 2 + ё • 13 _ іг, У-1 п,

; . + г, У 2 с • іг у ^ ё • і3 ■ г, У = Т п, У іа ;

г, У+1 + с • іг,у+1 + ё • 13 - V, У+1 - = ^п,

г, у+2 + с'ІІу+2 + ё 13 іг, У+2 =%

(28)

После ее решения получаются значения а, Ь, с, ё и подстановкой их и конкретного значения тока іг в (27) получаем МПС Тпігіа в точке, помеченной на рис.4,а квадратиком. Здесь же из (27) получается производная МПС по іг

(^п,іг,іа )іг = Ь + 2 • с • іг + 3 • ё ■ іГ . (29)

Полином и система уравнений, аналогичных (27) и (28), составляются для производной МПС по іа, т.е.

(30)

а + Ь ■ Іг ,■ _і + С ■ Іг, ,■ _і + ё ■ Іг ,■ _і у_!,іа )іа;

Т У -

а + Ь ■ іг, у + с • іг, У + ё •і3, У = (^п, і.іа )іа;

г, У

г, У

п, У,іа іа

(31)

а + Ь ■ іг, У+1 + с ■ іг, У+1 + ё ■ іг, У+1 _ (^п, У+1,іа )іа ;

а + Ь ■ Іг,у+2 + с • Іг,у+2 + ё • Іг,у+2 — (^п,У+2,іа )іа.

В итоге по (30) в точке (квадратик) на рис. 4,а при фиксированных значениях іг и іа получается значение ПРОИЗВОДНОЙ МПС ПО Іа, Т.е. (¥'п,іг,іа)іа.

Действия от (21) до (31) выполняются на рис. 3 во всех плоскостях того же типа, как на рис. 4.а при замене номера п на п-1, п+1 и п+2. Таким образом, в точках, отмеченных на рис. 4,6 квадратиками, получатся четверки значений МПС и его производных

^п-1,іг,іа ; ^п,іг,іа ; ^п+1,іг,іа ; ^п+2,іг,іа ; (32)

(^п—1,іг,іа )іа ; (^п,іг,іа )іа ; (^п+1,іг,іа)іа ; (^п+2,іг,іа )іа ;

(33)

(% - 1,іг,іа)іг ; (^п,іг,іа)іг ; (^п+1,іг,іа)іг ; (%, + 2 ,іг,іа )іг .

(34)

На линии с этими значениями (рис. 4,6) в точке с координатой С, (точка обозначена увеличенным кружочком) с использованием уже изложенных действий в итоге получаются значения частных производных по трем параметрам іг, іа и С, а именно

,іа )? ; (^,іг,іа )іг ; (^іма )іа . (35)

Если в качестве функции МПС Т использовать МПС реактора х¥г, то в соответствии с (6) и (7) получим Ьг (ХТ г,С,іг,іа)іг; -Мга ОТ г,£,іг,іа)іа; Ьг^ 0^ г,С,іг,іа)С;

если Та, ТО Маг = (У'а&іг,іа)іг'; Ьа = (^'а,^іг,іа)іа; Ьа^ =

(^а,Сіг,іак. Очевидно, что эти значения для конкретного сочетания параметров іг, іа и С, получены по 64 зна-

чениям каждой из функций ВМП в окружающих узлах из массивов (19).

Оперирование с функцией силы ^ет(^, 1а) про-

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

Примеры зависимости МПС и силы Еет от тока обмоток и величины С, есть в [2, 3]. Иллюстрация частных зависимостей индуктивных параметров обмоток ЛИЭД от взаимного аксиального смещения якоря и реактора дана на рис. 5-7. Здесь каждая из кривых соответствует конкретному сочетанию токов обмоток, хотя вместо токов использованы более универсальные величины - МДС, приходящиеся на паз реактора Епг и якоря Еп а, которые связаны с токами простыми соотношениями

Еп,т ~ п / а р ; Еп,а ~ /aNn|ap , (36)

где Ып, ар - одинаковые для обеих обмоток числа витков в пазу и параллельных ветвей.

Рис. 5. Зависимость индуктивных параметров обмоток ЛИЭД от взаимного смещения реактора и якоря: 1 - ¥пг = 0;

мГн

0.2

0.1

0

1-г —І.............*;........!

1*’ / І .і.

М,

■І-/-*---Маг :

Ш___________—-у..........ъ...

-12 0 12 24 36 мм

а - г = 100 кА; К, а = 0;

...& |.л?+....Л(гаі-Маг ..

мГн

0.1 0

-\2 0 12 24 36 мм

б - Рпг = 100 кА; ¥П} а = 100 кА Рис. 6. Варианты зависимости индуктивных параметров обмоток от взаимного смещения С,

Рис. 7. Варианты зависимости индуктивностей движения обмоток от взаимного смещения С, при сохранении К г = 100 кА и 1 - К а =0; 2 - К а = 50 кА; 3 - К а = 100 кА

(37)

Указанные значения МДС 100 кА соответствуют достигаемым при работе ЛИЭД максимумам токов. Значение Q = -12 мм соответствует Qo [6], Q = 0 - положение на геометрической "нейтрали", Q = 42 мм -совмещение зубцов реактора и якоря.

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

Представленная математическая модель, с учетом подготовленных данных и принципов оперирования с нелинейными дифференциальными индуктивными параметрами, позволила перейти к созданию алгоритма и программы расчета и проведению расчетного анализа рабочего режима ЛИЭД. В основу этого положено численное решение системы дифференциальных уравнений (10), (11) с использованием метода Рунге-Кутта 4 порядка в совокупности с комплексом остальных уравнений и соотношений. В целом это подобно представленному в [2], но теперь надо оперировать с парой производных токов из (10), (11)

Да Ar

ia ~ ; ir ~ ,

Aar Aar

где общий и частные определители имеют вид

Aar — (Lr + LCr ) •(La + Lca ) - M ra *M ar ;

Дг = (uCr ~ uRr ~ evr )(La + Lca )-(uCa _ uRa ~ eva )Mra ;

Aa = (uCa - uRa - eva )(Lr + Lcr )-(uCr _ uRr ~ evr ar .

Тестовый расчет и его результаты.

Конструктивные параметры ЛИЭД в целом соответствуют представленному в [1-6], то есть зазор между сердечниками реактора и якоря - 2 мм; диаметр по зазору - 416 мм; длина сердечников - 533 мм; Nn = 14; ap = 3; ma = 750 кг; mr = 1500 кг, общая высота

- около 2,6 м. Но в связи с изменением схемы электропитания обмоток для обеспечение того же режима используются новые параметры конденсаторов: UCro = UCao = 600 В; CEr = CEa = 0,17 Ф, т.е. соединенные ранее последовательно конденсаторы и обмотки теперь работают независимо.

Функции, характеризующие протекание переходного процесса в рабочем интервале, даны на рис.8-9. В этом интервале графики Qa и va показывают, что якорь осуществил свое ударное действие и вернулся в исходное состояние: = 0 и va = 0.

Функции даны в относительных единицах (о.е.), но для числовой оценки назовем их максимальные значения: токи /V и ia - Irmax =21,3 кА, Iamax = 19,2 кА; сила Fem - Femm = 1013 кН, сила Fgr - Fgrm =1695 кН; скорости реактора vr - vrm = 4,26 м/с и якоря va - vam = 2,65 м/с; ускорения реактора ar - arm = 665 м/с2 и якоря aa - aam = 822 м/с2; сдвиг якоря Qa - C,am = 11,3 мм; сдвиг реактора еще продолжается. Для прочих величин дадим масштабы - значения на 1 о.е.: Lr и Mra -0,1 мГн; для всех напряжений и ЭДС - 150 В. Заметим, что ряд графиков дан только для величин, относящихся к реактору. Для якоря они близки к ним, отличаясь примерно в той степени, как uCa и uCr.

Главными функциональными параметрами ЛИЭД считаются амплитуда механической силы Femm, дли-

тельность действия tiem и импульс этой силы [i, 2, б]

~iem

Sem ~ ^ Fem (t) *dt .

X 7 \ .

«іД Р

і

X-.......

Ik ~“Са ......

*10 ■--------

14

Рис. 8. Временные функции электрических величин

F ' Рц гет v 5 ту vr Л

У Js\ «г

1

2 4 \ \ 10 1 /. 2 14 мс 16 |і

v\ V / |і

\ У -Ji

\ /

\

Рис. 9. Временные функции механических величин

Но более объективную оценку дает введенный нами [2, 6] импульс силового воздействия непосредственно на грунт

Sgr = І Fgr dt

(39)

где timg+ - время действия силы Fgr до того момента, пока якорь еще движется в грунт (т.е. на рис. 1 и рис. 9 пока еще уа > 0).

Условный КПД ЛИЭД ^ = Wgr/Wel, где механиче-

timg+

ская энергия деформации грунта - Wgr = | Fgrvadt;

о

отданная конденсаторами электрическая энергия

Wel = 0,5[СЕг (иСго ~ иСгге1) + СЕа (иСао ~ иСаге1)] ,(40)

исГГе1, иСаге1 - остаточные их напряжения.

Получено Бет = 6557 Н-с; = 5469 Н-с; ^ =

0,198, длительности импульсов (по [2]) для токов 1Г и 1а - 12,5 и 15 мс; для сил Fem и Fgr 13,1 и 9,4 мс.

За один рабочий режим нагрев обмоток оказывается менее 1°С, несмотря на значительный максимум плотности тока - порядка 134 А/мм2. Но при многократном повторении импульсов с периодом Ттр =6 с нагревы становятся значительно больше.

Отметим, что полученные электросиловые параметры рабочего ударного режима ЛИЭД при последовательной схеме включения обмоток и двух конден-

t

саторов [2, б] по значениям и характеру изменения достаточно близки к полученным здесь. В частности общий для обмоток ток достигал 2G, 1 кА; сила Fem -Femm = 1GG8 кН, сила Fgr - Fgrm = 1б90 кН.

Выводы.

Представленная математическая модель является достаточно уникальной, так как основана на использовании сложных нелинейных трехпараметрических функций МПС, дифференциальных индуктивностей и взаимоиндуктивностей и электромагнитной силы. Эта модель позволяют эффективно исследовать динамический импульсный режим работы ЛИЭД с независимым включением обмоток реактора и якоря, предназначенного для возбуждения в земной коре слабых сейсмических колебаний. Ее усложняющим фактором также является совокупность взаимосвязанных электромагнитных, механических и тепловых величин.

СПИСОК ЛИТЕРАТУРЫ

1. Бару Ю.А., Мельник А.К., Ткаченко С.В., Милых В.И. Выбор типа и основных параметров линейного импульсного электродвигателя для невзрывных сейсмоисточников // Електротехніка і електромеханіка. - 2GG8. - № 3. - С. 41-44.

2. Милых В.И., Ткаченко С.В. Алгоритмизированная математическая модель линейного импульсного электродвигателя как ударного источника слабых сейсмических колебаний // Електротехніка і електромеханіка. - 2G11. - № 1. - С. 33-38.

3. Милых В.И., Ткаченко С.В. Силовые взаимодействия в линейном электродвигателе для источников сейсмических колебаний // Електротехніка і електромеханіка. - 2GG8. - № б.

- С. Зб-З9.

4. Милых В.И., Ткаченко С.В. Математическое обеспечение численного анализа импульсного возбуждения линейного электродвигателя с двухпараметрической функцией магнитного потокосцепления // Вісник Кременчуцького державного університету ім. М. Остроградського. - Кременчук: КДУ. - 2G1G. - Вип.3/2010(б2). Ч. 2. - С. 7G-73.

5. Meeker D. Finite Element Method Magnetics. Version 4.2. User’s Manual, September 2б, 2GG6 // http://femm.berlios.de, 2GG5.

6. Милых В.И., Ткаченко С.В. Анализ влияния механических параметров линейного электродвигателя на эффективность его работы как источника сейсмических колебаний // Гірн. елекгромеханіката автоматика: Наук.-техн. зб. - 2G1G.

- Вин. 85 - С. 124-13G.

Bibliography (transliterated): 1. Baru JU.A., Melnik A.K., Tkachenko C.V., Mil'ih V.I. V'ibor tipa i ocnovn'ih parametrov line'nogo im-pulcnogo elektrodvigatelja dlja nevzr'ivn'ih ce'cmoictochnikov // Elektrotehnika і elektromehanika. - 2GG8. - № 3. - C. 41-44. 2. Mil'ih V.I., Tkachenko C.V. Algoritmizirovannaja matematicheckaja model line'nogo impulcnogo elektrodvigatelja kak udarnogo ictochnika clab'ih ce'cmicheckih kolebani' // Elektrotehnika і elektromehanika. - 2G11. -№1. - C. 33-38. 3. Mil'ih V.I., Tkachenko C.V. Cilov'ie vzaimode'ctvija v line'nom elektrodvigatele dlja ictochnikov ce'cmicheckih kolebani' // Elektrotehnika і elektromehanika. - 2GG8. - № б. - C. Зб-З9. 4. Mil'ih V.I., Tkachenko C.V. Matematicheckoe obecpechenie chiclennogo analiza impulcnogo vozbujdenija line'nogo elektrodvigatelja c dvuhparametrichecko' funktsie' magnitnogo potokoctseplenija // Vicnik Kremenchutskogo derjavnogo univercitetu im. M. Octrogradckogo. -Kremenchuk: KDU. - 2G1G. - Vip^GiG^). CH. 2. - C. 7G-73. 5. Meeker D. Finite Element Method Magnetics. Version 4.2. User's Manual, September 2б, 2GG6 // <http://femm.berlios.de/>, 2005. б. Mil'ih V.I., Tkachenko C.V. Analiz vlijanija mehanicheckih parametrov line'nogo elektrodvigatelja na effektivnoct ego rabot'i kak ictochnika ce'cmicheckih kolebani' // Girn. elektromehanika ta avtomatika: Nauk.-tehn. zb. - 2G1G. - Vip. 85 - C. 124-13G.

Поступила 10.03.2011

МтыхВладимирИванович, д.т.н., проф.,

Ткаченко Сергей Валериевич Национальный технический университет "Харьковский политехнический институт" кафедра "Электрические машины"

6iGg2, Харьков, ул. Фрунзе, 21

тел. (G57) 707-б5-14; e-mail: mvikpi@kpi.kharkov.ua

Milykh V.I., Tkachenko S.V.

Mathematical model and analysis of the impulsive mode of operations of linear electric motor at the independent feed of his windings.

A linear electric motor is intended for excitation of vibrations of the earth's crust at the search of minerals. His mathematical model is formulated. She corresponds the mode of shock power impulse. Its feature is a presence of nonlinear functions of differential own and mutual inductances windings. These functions depend on three parameters. The same is function of electromagnetic force which operates on a reactor and anchor.

Key words — linear impulsive electric motor, mathematical model, own and mutual inductance, electromagnetic force, functions of three parameters.