Математическая модель напряженной магнитоанизотропной структуры устройства телеметрии силового воздействия
Автор: Дубинин Александр Ефимович, Дубинин Александр Александрович, Цаплин Николай Николаевич
Журнал: Известия Самарского научного центра Российской академии наук @izvestiya-ssc
Рубрика: Информатика, вычислительная техника и управление
Статья в выпуске: 4-1 т.18, 2016 года.
Бесплатный доступ
Рассматривается математическая модель объемной напряженной магнитоанизотропной структуры (НМАС) устройства телеметрии силового воздействия. При этом обобщенная физическая модель НМАС представлена в виде четверти ферромагнитного кольцевого чувствительного элемента (ЧЭ) в электромагнитном поле при силовом воздействии, которая описывается системой дифференциальных уравнений теорий электромагнитного поля, ферромагнетизма и упругости в сферических полярных координатах. Магнитная цепь НМАС описывается системой алгебраических уравнений. Совместное решение системы дифференциальных и алгебраических уравнений позволяет определить электрические и магнитные параметры ЧЭ и всего устройства, а также построить его выходную статическую характеристику.
Объемная магнитная анизотропная структура, физическая модель, математическая модель, чувствительный элемент, структурная схема, выходная статическая характеристика, устройство телеметрии
Короткий адрес: https://sciup.org/148204645
IDR: 148204645
Текст научной статьи Математическая модель напряженной магнитоанизотропной структуры устройства телеметрии силового воздействия
шары 4 растягивается образуя напряжение растяжения σр [1, 2].
Под действием механических напряжений σс и σр магнитная проницаемость ЧЭ2 уменьшается, а магнитная проницаемость ЧЭ3 увеличивается, согласно теории ферромагнетизма [1, 2], что приводит к изменению частоты выходных сигналов f 1 и f 2 генераторов Г1 и Г2. Эти сигналы с выхода обмоток 7, 8 ЧЭ 2, 3 подаются на входы смесителя См. С выхода смесителя См сигнал разностной частоты поступает на антенну А2 передатчика и принимается антенной А1 приемника, где обрабатывается усилителем разностной частоты УРЧ, амплитудным ограничителем АО, частотным детектором ЧД с гетеродином ГТ и подается на регистрирующий прибор РП, фиксируя силовое воздействие.
Обобщенная физическая модель НМАС ЧЭ сжатия 2 и растяжения 3 представлена на рис. 2 в виде четверти ферромагнитного кольца, сечение которого ориентировано в сферической полярной системе координат таким образом, чтобы вектор нормальных напряжений σ совпадал с направлением магнитного поля Н С, а вектор касательных напряжений σ был с ним перпендикулярен, где Р С – сила воздействия на кольцо направлена по радиусу ЧЭ под прямым углом к полю Н С; а = r н – r в – толщина кольца; r н, r в – наружный и внутренний радиусы кольца; b – ширина кольца; r 0 – средний радиус элемента; φ – угол между осью симметрии сечения и радиусом r в; Er и Eb напряженности электрического поля на гранях сечения; B С – магнитная индукция внешнего электромагнитного поля; Н С – напряженность магнитного поля в сечении кольца; i – ток, протекающий через обмотку возбуждения; i в – вихревой ток; θ – угол

Рис. 3. Схема замещения магнитной цепи упругих элементов

Рис. 1. Структурная схема устройства телеметрии силового воздействия
между осью симметрии сечения и радиусом r 0 по длине кольца с в радианах.

Рис. 2. Обобщенная физическая модель НМАС ЧЭ
-
1. Картина электромагнитного поля в любом сечении кольцевого ЧЭ равномерна и не зависит от угла 9.
-
2. Магнитная индукция B С электромагнитного поля, в котором находится ЧЭ, изменяется по
-
3. Картина силового поля с учетом напряжений сжатия, растяжения и сдвига зависит от угла 9 ЧЭ и может быть рассмотрена с использованием
синусоидальному закону и имеет только нор-
мальную составляющую.
3-й гипотезы прочности [1].
Общее выражение для определения нор-
мальных механических напряжений сжатия и
растяжения в кольцевом ЧЭ под действием силы
нажатия Р записывается в виде [1]:
^ Н, b =
0,5 P ba
±3,82 r - +1 6 r^ +11 Cos 9 a I a )
, (1)
Математическая модель обобщенной физической модели НМАС описывается дифференциальными уравнениями теорий электромагнитного поля, ферромагнетизма и упругости в полярных координатах. При этом принимаются следующие допущения и предположения.
где верхние знаки относятся к наружной поверхности кольцевого ЧЭ, а нижние знаки - к внутренней поверхности ЧЭ.
При 9 = 0 наружная поверхность ЧЭ сжимается под действием стс а внутренняя - растягивается под действием ст р (см. рис. 2). При 9 = п/2 картина меняется на обратную: наружная поверхность растягивается под действием стр, а внутренняя сжимается под действием стс. Картина распределения напряжений повторяется для каждой четвери кольца.
Механические напряжения растяжения стр ЧЭ2 при 9 = п/2 являются причиной возникновения силы растяжения Р Р кольцевого ЧЭ3, который расположен во взаимно перпендикулярной плоскости - рис. 1. При этом внутренняя поверхность ЧЭ3 сжимается, а наружная - растягивается, т. е. картина распределения механических напряжений стс и стр аналогична ЧЭ2.
Используя обобщенную физическую модель ЧЭ на рис. 2, определяются максимальные напряжения сжатия с с при 9 = 0 и растяжения ст р при 9 = π/2 согласно выражению (1):
1,1- P - r 1,91- P - r
° =----г-0 ° р =----2—0
ba 2 ba 2
;
.
Учитывая третью гипотезу прочности для ЧЭ °э = V ° 2 + 4т 2 и максимальные касательные напряжения т = (ст р - сс)/2, определяются эквивалентные механические напряжения при работе ЧЭ2 на сжатие и ЧЭ3 на растяжение:
эс с р с р сэр = V2сР + с2 - 2срсс • (3)
4. Под действием эквивалентных механических напряжений сжатия сэс и растяжения сэр магнитная проницаемость ЧЭ2 уменьшается, а ЧЭ3 – увеличивается [1].
Нэс , р = Рн ±Н (4)
где (–) – при сжатии элемента 2; (+) – при растяжении элемента 3; μн – начальная магнитная проницаемость элементов 2 и 3; ∆μ – изменение магнитной проницаемости.
При принятых допущениях и предположениях при одновременном воздействии силового поля сжатия или растяжения и поля возбуждения состояние соответствующего ферромагнитного ЧЭ описывается дифференциальными уравнениями в сферических полярных координатах [1]:
д Нс д r
= УФ E Ф ;
д Нс дф
= Y rEr ;
д Er д Е ф д Вс ^^^^^^^^^^^^^^^^^^^^^^^^^^. ^^^^^^^^ ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^_ ^^^^^^^^ ^^^^^^^^^^^^^^^^^^^^^^^^^^_ д r дф д t
- jюp ° Hc ;
дВг _ 1 д Вг I д Нг д Нг д Нг
-^- + =C-I ^=C + ^=C + вHc В2 дНс I д° д. Вт
C Cp c
E (w + ^ w L ° .
(д r дф де ) э
где сэ - соответствующие эквивалентные напряжения кольцевых ЧЭ; у r , уф удельные электрические проводимости по направлениям; E r , E ф; H C - комплексы действующих значений напряженностей электрического и магнитного полей по поверхности сечения кольца по переменным r , ф и 9; с - длина кольца; r 0 - средний радиус элемента; ю - круговая частота; цсэ - магнитная проницаемость от механического напряжения сэ; w, , w ф, w 9 - перемещение точек структуры по радиусу r , углам ф и 9; Е - модуль упругости; 9 - угол в радианах.
Граничные условия в этом случае имеют вид (рис. 2).
Нс ( r, ф ) = Нс ( r ) при r = r 0 ± а/ 2, ф = 0; r в< r ^ r н; (8)
Н С (ф, r ) = Н С (ф) при bb - arcsin — < ф < arcsin — (9)
2 r í 2 r â
Ф 2 при малых углах ф b « L = J r 0 d ф .
Ф 1
b
Нс ( r ,ф) = 0 при ф = ±arcsin— ;
2 r â
Нс (ф, r ) = 0 ; при r = r 0 ± а/ 2. (10)
c = J r 0 d 9 = r 09 ; 0 < 9 < 2n. (11)
Уравнение теории упругости (7) рассмотрено с учетом 3-й гипотезы прочности, и его решение представлено в виде формул (3). Уравнение ферромагнетизма (6) для НМАС решено исходя из плотности свободной энергии ферромагнетика и термодинамических соотношений для обратимых магнитомеханических процессов при растяжении-сжатии, что позволило определить магнитные проницаемости ЧЭ цс их изменения при статических и динамических воздействиях P и P(t) [1]:
K 'с р ^as 2 цн’ r0P
Р. = Нн ±ДР; Н° = НН ± • ^2ba2 ^12> где K’с, р – коэффициент при сжатии и растяжении; Xs - изотропная магнитострикция материалаЧЭ; as коэффициент намагниченности материала.
Из решения уравнений электромагнитного поля (5) при выбранных граничных условиях (8) – (11) определены плотности токов, вихревые I В и намагничивающие I Ф токи и сопротивления этим токам Z В и Z Ф [1]. Последнее зависит от силового воздействия и имеет вид:
Z ф= \р. +
У
K 'с,р ^S«S н Н r0рЛ
пВ 2 ba 2
K П ej , (13) )
где S э = 2( a + b ) aВ – сечение эквивалентного контура; aВ – глубина проникновения магнитного поля в материал; a - угол сдвига между намагничивающим током I Ф и магнитным потоком Ф; K П = K Т – коэффициент трансформации.
Выражение (13) позволяет определить электрические потери и угол магнитных потерь в ЧЭ при силовых воздействиях. На рис. 3 показана схема замещения ЧЭ, НМАС которых работает при сжатии-растяжении.
Для схемы замещения по втором у закону Кирхгофа и после перехода к магнитным индукциям и далее к напряженностям поля система алгебраических уравнений преобразуется к виду:
HХДХ - H 2 Д 2 к
^^^^^^^^
T r ^ ) ^2 S' Ck w 2 J
U 1 toSw ^
^^^^^^^в
H х Д х + H 2 Д 2 к
^^^^^^^в
2S^_ mSCw2 J
= 0
.
В результате решения системы уравнений определяется напряженность магнитного поля H 2 и выходные действующие напряжения U Н на Z Н с выхода обмоток 7, 8 – рис. 1.
H = U 2 ^Ckw^
2 ro 9
г
Д н ± к
K 'с,р ^2 Д 2 Г 0 P^
пВ 2 ba 2
J
H 2 r 0 0
U н 7,8 2^н
, w 2
г
= U 2 ^ р.с C К K t Дн ±-к
Zн
K 'с,р <а Д 2 Г о P )
nB 2 ba 2 7
Изменение магнитных проницаемостей (12) магнитных элементов 2 и 3 приводит к изменению индуктивности обмоток возбуждения 5 и 6 и изменению частот генераторов Г1 и Г2 [1].
f Г 1, Г 2
2 n ba 2 w 2 C k
r 0
Д ± к
K с,р ^а .Д . Г о P )
п В 2 ba 2
в 1
, (16)
J_
где w 1 – число витков обмоток возбуждения 5, 6 элементов 2, 3; C к – емкости контурных конденсаторов генераторов Г1 и Г2.
Частоты f Г1 и f Г2 генераторов Г1 и Г2 через измерительные обмотки 7, 8 подаются на входы смесителя См, с выхода которого сигнал ∆ U н = U н7 – U н8 частотой ∆ f = f Г1 – f Г2 поступает в антенну А 2 передатчика и принимается антенной А 1
приемника, где обрабатывается и подается на регистрирующий прибор РП.
Все узлы передатчика и приемника выполнены на интегральных микросхемах.
Предлагаемая методика расчета использовалась для УТСВ со следующими параметрами:
-
1. элемент 2 сжатия 5×3×1,5 мм из феррита 2000 НМ1 w 1 20; w 2 = 10; K Т = 2; Z н = 1×104 Ом, Ck = 90 пФ; μн = 25·104 Гн/м; U 1 = 5B – напряжение генератора Г1.
-
2. элемент растяжения 3 10×6×3 мм из феррита 2000 НМ1; U 1 = 5B – напряжение генератора Г2 w 1 30; w 2 = 15; K Т = 2; Z н = 1×104 Ом, C к = 42 пФ; μн = 28·10–4 Гн/м.
В результате расчета получено:
-
1) при Р = 0 f Г1 = 1,5×106 Гц; U н7 = 210 мВ; f Г2 = 1×106 Гц; U н8 = 73 мВ; на выходе смесителя См в А 2 ∆ f = 0,5 МГц; ∆ U н = 137 мВ; f РП = ∆ f = f ГТ = 0;
-
2) при Р = 30 Н f Г1 = 1,8×106 Гц; U н7 = 150 мВ; f Г2 = 0,95×106 Гц; U н8 = 80 мВ; на выходе смесителя См в А 2 ∆ f = 0,85 МГц; ∆ U н = 70 мВ; f РП = ∆ f – f ГТ = 0,35 МГц; U РП = 2,0 В.
Экспериментальные исследования подтвердили правильность теоретических положений и формул – погрешность не превышает 10%.
Список литературы Математическая модель напряженной магнитоанизотропной структуры устройства телеметрии силового воздействия
- Дубинин А.Е., Дубинин А.А. Математическая модель напряженных магнитоанизотропных структур устройства телеметрического контроля тягового электродвигателя//Вестник транспорта Поволжья. 2013. № 1 (34). С. 34-43.
- Патент на изобретение №2467441, Российская Федерация, МКИ Н01R39/58. Устройство для измерения нажатия щетки на коллектор электрических машин/А. Е. Дубинин, А. А.Дубинин, Р. Е. Капитуров. Заявл. 03.08.11; Опубл. 20.11.12, Бюллетень изобретений № 232.