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

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

Gespeichert in:
Bibliographische Detailangaben
Datum:2015
Hauptverfasser: Пилипенко, О.В., Горев, Н.Б., Доронин, А.В., Заболотный, П.И., Коджеспирова, И.Ф.
Format: Artikel
Sprache:Russian
Veröffentlicht: Інститут технічної механіки НАН України і НКА України 2015
Schriftenreihe:Техническая механика
Online Zugang:http://dspace.nbuv.gov.ua/handle/123456789/100791
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Определение комплексного коэффициента отражения электромагнитной волны с помощью зондовых измерений / О.В. Пилипенко, Н.Б. Горев, А.В. Доронин, П.И. Заболотный, И.Ф. Коджеспирова // Техническая механика. — 2015. — № 4. — С. 139-147. — Бібліогр.: 10 назв. — рос.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-100791
record_format dspace
spelling irk-123456789-1007912016-05-27T03:03:47Z Определение комплексного коэффициента отражения электромагнитной волны с помощью зондовых измерений Пилипенко, О.В. Горев, Н.Б. Доронин, А.В. Заболотный, П.И. Коджеспирова, И.Ф. Целью данной статьи является разработка метода измерения комплексного коэффициента отражения с использованием электрических зондов. Метою цієї статті є розробка методу вимірювання комплексного коефіцієнта відбиття з використанням електричних зондів. The aim of this paper is to develop a method for complex reflection coefficient measurement using electrical probes. 2015 Article Определение комплексного коэффициента отражения электромагнитной волны с помощью зондовых измерений / О.В. Пилипенко, Н.Б. Горев, А.В. Доронин, П.И. Заболотный, И.Ф. Коджеспирова // Техническая механика. — 2015. — № 4. — С. 139-147. — Бібліогр.: 10 назв. — рос. 1561-9184 http://dspace.nbuv.gov.ua/handle/123456789/100791 621.002.56 ru Техническая механика Інститут технічної механіки НАН України і НКА України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language Russian
description Целью данной статьи является разработка метода измерения комплексного коэффициента отражения с использованием электрических зондов.
format Article
author Пилипенко, О.В.
Горев, Н.Б.
Доронин, А.В.
Заболотный, П.И.
Коджеспирова, И.Ф.
spellingShingle Пилипенко, О.В.
Горев, Н.Б.
Доронин, А.В.
Заболотный, П.И.
Коджеспирова, И.Ф.
Определение комплексного коэффициента отражения электромагнитной волны с помощью зондовых измерений
Техническая механика
author_facet Пилипенко, О.В.
Горев, Н.Б.
Доронин, А.В.
Заболотный, П.И.
Коджеспирова, И.Ф.
author_sort Пилипенко, О.В.
title Определение комплексного коэффициента отражения электромагнитной волны с помощью зондовых измерений
title_short Определение комплексного коэффициента отражения электромагнитной волны с помощью зондовых измерений
title_full Определение комплексного коэффициента отражения электромагнитной волны с помощью зондовых измерений
title_fullStr Определение комплексного коэффициента отражения электромагнитной волны с помощью зондовых измерений
title_full_unstemmed Определение комплексного коэффициента отражения электромагнитной волны с помощью зондовых измерений
title_sort определение комплексного коэффициента отражения электромагнитной волны с помощью зондовых измерений
publisher Інститут технічної механіки НАН України і НКА України
publishDate 2015
url http://dspace.nbuv.gov.ua/handle/123456789/100791
citation_txt Определение комплексного коэффициента отражения электромагнитной волны с помощью зондовых измерений / О.В. Пилипенко, Н.Б. Горев, А.В. Доронин, П.И. Заболотный, И.Ф. Коджеспирова // Техническая механика. — 2015. — № 4. — С. 139-147. — Бібліогр.: 10 назв. — рос.
series Техническая механика
work_keys_str_mv AT pilipenkoov opredeleniekompleksnogokoéfficientaotraženiâélektromagnitnojvolnyspomoŝʹûzondovyhizmerenij
AT gorevnb opredeleniekompleksnogokoéfficientaotraženiâélektromagnitnojvolnyspomoŝʹûzondovyhizmerenij
AT doroninav opredeleniekompleksnogokoéfficientaotraženiâélektromagnitnojvolnyspomoŝʹûzondovyhizmerenij
AT zabolotnyjpi opredeleniekompleksnogokoéfficientaotraženiâélektromagnitnojvolnyspomoŝʹûzondovyhizmerenij
AT kodžespirovaif opredeleniekompleksnogokoéfficientaotraženiâélektromagnitnojvolnyspomoŝʹûzondovyhizmerenij
first_indexed 2025-07-07T09:22:21Z
last_indexed 2025-07-07T09:22:21Z
_version_ 1836979472724656128
fulltext 139 УДК 621.002.56 О. В. ПИЛИПЕНКО, Н. Б. ГОРЕВ, А. В. ДОРОНИН, П. И. ЗАБОЛОТНЫЙ, И. Ф. КОДЖЕСПИРОВА ОПРЕДЕЛЕНИЕ КОМПЛЕКСНОГО КОЭФФИЦИЕНТА ОТРАЖЕНИЯ ЭЛЕКТРОМАГНИТНОЙ ВОЛНЫ С ПОМОЩЬЮ ЗОНДОВЫХ ИЗМЕРЕНИЙ Измерения комплексного коэффициента отражения широко применяются для контроля различных параметров материалов. Целью данной статьи является разработка метода измерения комплексного коэф- фициента отражения с использованием электрических зондов. Для случая трех зондов, установленных на произвольном расстоянии друг от друга, получено биквадратное уравнение, связывающее модуль ком- плексного коэффициента отражения с токами полупроводниковых детекторов, соединенных с зондами. Показано, что в силу того, что модуль коэффициента отражения не превышает единицу, он однозначно определяется из этого уравнения как его меньший положительный корень. Фаза комплексного коэффици- ента отражения определяется из двух квадратурных сигналов, которые легко находятся из токов полупро- водниковых детекторов при известном модуле коэффициента отражения. Для проведения измерений ком- плексного коэффициента отражения в диапазоне частот межзондовое расстояние удобно выбрать равным одной восьмой длины волны электромагнитного излучения в волноводе при максимальной частоте. Предложенный метод позволяет исключить амплитуду напряженности электрического поля падающей электромагнитной волны, что дает возможность значительно смягчить требования к стабильности выход- ной мощности генератора. Отличительной особенностью предложенного метода является простота его аппаратной реализации, что открывает путь к разработке нового класса измерителей комплексного коэф- фициента отражения. Вимірювання комплексного коефіцієнта відбиття широко використовуються для контролю різнома- нітних параметрів матеріалів. Метою цієї статті є розробка методу вимірювання комплексного коефіцієн- та відбиття з використанням електричних зондів. Для випадку трьох зондів, розташованих на довільній відстані один від одного, одержано біквадратне рівняння, що пов’язує модуль комплексного коефіцієнта відбиття зі струмами напівпровідникових детекторів, з’єднаних із зондами. Показано, що в силу того, що модуль комплексного коефіцієнта відбиття не перевищує одиниці, він однозначно визначається з цього рівняння як його менший додатний корінь. Фаза комплексного коефіцієнта відбиття визначається з двох квадратурних сигналів, які легко знаходяться зі струмів напівпровідникових детекторів при відомому модулі коефіцієнта відбиття. Для проведення вимірювань комплексного коефіцієнта відбиття в діапазоні частот міжзондову відстань зручно обрати рівною одній восьмій довжини хвилі електромагнітного ви- промінювання в хвилеводі при максимальній частоті. Запропонований метод дозволяє виключити амплі- туду напруженості електричного поля падаючої електромагнітної хвилі, що дає змогу значно пом’якшити вимоги щодо стабільності вихідної потужності генератора. Відмітною особливістю запропонованого ме- тода є простота його апаратної реалізації, що відкриває шлях до розробки нового класу вимірювачів ком- плексного коефіцієнта відбиття. Complex reflection coefficient measurements are widely used in materials characterization. The aim of this paper is to develop a method for complex reflection coefficient measurement using electrical probes. For the case of three probes placed an arbitrary distance apart, a biquadratic equation that relates the complex reflection coef- ficient magnitude to the currents of the semiconductor detectors connected to the probes is derived. Due to the fact that the complex reflection coefficient magnitude is no greater than unity, it is unambiguously determined from that equation as its smaller positive root. The complex reflection coefficient phase is determined from two quadrature signals, which are easy to find from the semiconductor detector currents once the complex reflection coefficient magnitude is known. To measure the complex reflection coefficient over a frequency range, it is con- venient that the interprobe distance be equal to one eighth of the guided operating wavelength at the maximum frequency. The proposed method allows one to eliminate the incident wave electric field amplitude, this making it possible to greatly alleviate the requirements for oscillator output power stability. A distinctive feature of the proposed method is the simplicity of its hardware implementation, which opens up the way to the development of a new class of vector reflectometers. Ключевые слова: комплексный коэффициент отражения, электриче- ский зонд, полупроводниковый детектор, выходная мощность генератора, межзондовое расстояние. Для контроля параметров материалов, в частности процентного содер- жания смесей, широко применяются измерения диэлектрической проницае- мости. Одним из методов определения диэлектрической проницаемости яв- ляется измерение комплексного коэффициента отражения. Традиционный  О. В. Пилипенко, Н. Б. Горев, А. В. Доронин, П. И. Заболотный, И. Ф. Коджеспирова, 2015 Техн. механика. – 2015. – № 4. 140 метод измерения комплексного коэффициента отражения основан на исполь- зовании согласованной конструкции двенадцатиполюсного измерительного преобразователя [1]. Основные проблемы метода – оптимизация структуры и характеристик преобразователя, а также разработка способов и средств его калибровки [2]. Константы, необходимые для определения комплексного ко- эффициента отражения, находятся по результатам предварительной калиб- ровки измерителя с использованием не менее четырех эталонов комплексно- го коэффициента отражения. Кроме того, технология изготовления согласо- ванных измерителей, например в волноводном исполнении [3], сложна. Со- вершенствование компьютерной техники позволило сформулировать подход, согласно которому находят применение двухполюсные и четырехполюсные эталоны с параметрами, полученными расчетным путем. В работе [4] пред- ложен векторный рефлектометр для измерения комплексного коэффициента отражения с измерительным преобразователем на основе крестообразного Е-плоскостного волноводного делителя, схематический вид которого пока- зан на рис. 1, где 1 – крестообразный Е-плоскостной делитель с плечами 2, 3, 4, 5; 6 – генератор; 7 – вентиль; 8, 15 – направленные ответвители; 9 – иссле- дуемая нагрузка; 10, 11, 14 – детекторы; 12, 13 – согласованные измерители мощности; 16 – измерительный преобразователь. Рис. 1 Особенностью Е-плоскостного волноводного креста является возмож- ность расчета матрицы рассеяния строгими электродинамическими метода- ми. Элементы матрицы рассеяния имеют весьма плавную зависимость от ча- стоты во всей рабочей полосе. Как показано в работе [5], обработка результа- тов измерений на основе принципа голографии с тремя опорными сигналами [6] позволяет заменить решение системы квадратных уравнений на решение системы линейных уравнений с использованием аппарата регуляризации Ти- хонова [7]. В работе [8] предложен измеритель комплексного коэффициента отра- жения, построенный на основе двухканального двухдетекторного волновод- но-щелевого преобразователя, схематический вид которого приведен на рис. 2, где 1 – сверхвысокочастотный генератор; 2 – двунаправлен- ный ответвитель; 3, 4 – согласованные нагрузки; 5, 6 – детекторные головки; 7, 8 – отрезки волновода прямоугольного сечения; 9 – исследуемый образец. Преимуществом предлагаемого измерителя, по сравнению с рас- смотренным выше измерителем на основе крестообразного Е-плоскостного 141 волноводного делителя, является меньшее количество детекторов – всего два, однако, как видно из рис. 2, его конструкция остается довольно сложной и для определения фазы комплексного коэффициента отражения требуется проводить измерения на двух близких частотах. Рис. 2 С другой стороны, задача определения коэффициента отражения возни- кает также при измерении относительного перемещения объекта интерфе- ренционным методом с использованием электрических зондов. Так, в работе [9] предложена методика измерения относительного перемещения объекта с помощью зондов, установленных в волноводной секции на расстоянии одной восьмой длины волны зондирующего электромагнитного излучения друг от друга. В этой методике перемещение определяется по двум квадратурным сигналам, которые находятся из токов полупроводниковых детекторов, со- единенных с зондами. В выражения, связывающие квадратурные сигналы с токами полупроводниковых детекторов, входит модуль неизвестного ком- плексного коэффициента отражения объекта, который также определяется по токам полупроводниковых детекторов. При известном расстоянии между объектом и зондами из квадратурных сигналов может легко быть найдена и фаза комплексного коэффициента отражения объекта. Аппаратная реализа- ция подхода, лежащего в основе зондовой методики, значительно проще ап- паратной реализации методов измерения комплексного коэффициента отражения на основе крестообразного Е-плоскостного волноводного де- лителя и на основе двухканального двухдетекторного волноводно-щелевого преобразователя. Однако эта методика позволяет проводить измерения толь- ко на одной частоте, а именно на той, при которой межзондовое расстояние равно одной восьмой длины волны зондирующего электромагнитного излу- чения в волноводе. Поэтому для обеспечения возможности измерения ча- стотной зависимости комплексного коэффициента отражения в диапазоне частот лежащий в основе этой методики подход должен быть обобщен на случай произвольного межзондового расстояния. Целью данной статьи является разработка метода измерения комплекс- ного коэффициента отражения с использованием электрических зондов. Рассмотрим три зонда 1, 2, 3, соединенных с полупроводниковыми де- текторами с квадратичной вольт-амперной характеристикой. Зонды располо- жены в волноводной секции между генератором электромагнитного излуче- ния и образцом, коэффициент отражения которого измеряется, на расстоянии 142 l друг от друга, причем ближе к образцу расположен зонд 3. Схема измерения приведена на рис. 3. Как видно из сравнения с рис. 1 и рис. 2, с точки зрения аппаратной реализации зондовые измерения комплексного коэффициента отражения значительно проще методов, предложенных в работах [4] и [8]. Рис. 3 Токи детекторов 321 ,, JJJ следующим образом выражаются через мо- дуль r и фазу  комплексного коэффициента отражения образца   cosпад rrEkJ 21 22 11 , (1)   ]sinsin2coscos2 1[]cos21[ 22 ïàä2 22 ïàä22   rr rEkrrEkJ , (2)   ]sinsin2coscos2 1[]cos21[ 22 ïàä3 22 ïàä33   rr rEkrrEkJ , (3)     gγ 4 L , g   l4 , g   l8 , (4) где падE – модуль амплитуды напряженности электрического поля падаю- щей электромагнитной волны; 1k , 2k , 3k – коэффициент пропорциональности между квадратом модуля амплитуды напряженности электрического поля и током детектора для детекторов 1, 2, 3 соответственно; L – расстояние между образцом и зондом 1; gγ – длина волны зондирующего электромагнитного излучения в волноводе. Задача заключается в том, чтобы найти модуль r и фазу  комплексного коэффициента отражения по измеренным токам детекторов. Соотношения (1) – (3) можно рассматривать как систему трех уравнений с четырьмя неизвестными: r . sin , cos , 2 падE . Для того чтобы число уравнений соответствовало числу неизвестных, эту систему следует допол- нить основным тригонометрическим тождеством 1cossin 22  . (5) Вначале получим уравнение для определения r . Разделив (2) и (3) на (1), получим 143 a Jk Jk rr rrr    12 21 2 2 cos21 sinsin2coscos21 , b Jk Jk rr rrr    13 31 2 2 cos21 sinsin2coscos21 , откуда      r ra a 2 11 sinsincoscos 2  , (6)      r rb b 2 11 sinsincoscos 2  . (7) Соотношения (6) и (7) можно рассматривать как систему двух уравнений с двумя неизвестными sin и cos . Решение этой системы имеет вид               r Ar ba baab r r 2 1 sincossincos cos1cos1 2 1 sin 22      , (8)             B r r ba ba r r 2 1 sincossincos sin1sin1 2 1 cos 22      . (9) Подставив (8) и (9) в (5), получим   22422 214 BArrr  . (10) Из (10) следует, что квадрат модуля комплексного коэффициента отра- жения 2rt  удовлетворяет следующему квадратному уравнению 01 4 2 22 2         t BA t . (11) Так как свободный член этого уравнения положителен, а коэффициент при t отрицателен (последнее следует из положительности t и 12 t ), то корни этого уравнения положительны. Так как коэффициент при t2 и свобод- ный член равны единице, то один корень этого уравнения больше единицы, а другой – меньше. Это позволяет однозначно определить из этого уравнения квадрат модуля комплексного коэффициента отражения 2r как меньший ко- рень этого уравнения, так как r , а следовательно, и 2r не может быть боль- ше единицы. Таким образом, из (11) имеем r 11 2 1 2 2 2222           BABA . При известном модуле коэффициента отражения r содержащие инфор- мацию о фазе коэффициента отражения sin и cos находятся из (8) и (9). 144 При известных sin и cos аргумент этих функций     g L4 можно найти следующим образом n 2 ,                    ,0,0sin,2 sin arctg ,0, sin arctg ,0,0sin, sin arctg cos cos cos cos cos cos где  – неразвернутая фаза; n – целое число. Отсюда получается следующее выражение для искомой фазы  ком- плексного коэффициента отражения n L g     2 4 . (12) Входящее в (12) целое число n выбирается из условия, чтобы фаза  лежала в интервале от нуля до 2 . В настоящее время общепринятым значением межзондового расстояния является 8gl  . При увеличении отношения gl  до 1/4 уравнение (3) пе- реходит в уравнение (1), т. е. количество уравнений в системе (1) – (3), (5) становится меньше количества неизвестных. Поэтому при проведении изме- рений комплексного коэффициента отражения в диапазоне частот межзондо- вое расстояние l удобно выбрать равным одной восьмой минимальной длины волны в волноводе ming , т. е. равным одной восьмой длины волны, соответ- ствующей максимальной частоте, поскольку при этом во всем диапазоне ча- стот будет выполняться условие 81 gl . При проведении измерений в миллиметровом диапазоне длин волн обес- печение межзондового расстояния, равного 8ming , является достаточно сложной задачей. Однако при этом можно поступить по-другому, а именно, с помощью предварительных измерений с использованием короткозамыкаю- щего поршня с микрометрическим винтом определить фактическое межзон- довое расстояние [10]. Пусть расстояние 12l между зондами 1 и 2 установле- но с некоторой погрешностью, т. е.     1 8 min 12 g l . В этом случае при проведении измерений с короткозамыкающим поршнем ( 1r ,  ) урав- нения (1) и (2) принимают вид   cos121relJ , (13)   2sin122 relJ , (14) где введены следующие обозначения 145 2 1 1 1 падEk J J rel  , 2 2 2 2 падEk J J rel  (величины в знаменателях этих формул представляют собой токи соответ- ствующих детекторов в режиме согласованной нагрузки). Формулу для 2relJ можно переписать как    ]2sincos2cossin1[22 relJ , откуда        cos 22cossin1 2sin 2relJ . (15) Пусть minx и maxx – положения короткозамыкающего поршня, в кото- рых ток 1relJ достигает минимума и максимума соответственно. Очевидно, что в этих положениях поршня 1)(cos min  x , 1)(cos max  x , 0)(sin maxmin,  x . (16) Из (15) с учетом (16) имеем   2 2)( 2sin min2   xJ rel , (17)   2 )(2 2sin max2 xJ rel  . (18) При разумном предположении, что 1 , т. е. погрешность установки не превышает заданного межзондового расстояния, величина  определяется из (17) и (18) однозначно. В качестве примера на рис. 4 приведены рассчитанные зависимости нор- мированных токов 1relJ , 2relJ от координаты x короткозамыкающего поршня для g 3 см и  0,2. Расчет проводился по формулам (4), (13), (14) с добавлением моделирующей шум случайной компоненты, равномерно распределенной на промежутке от –0,25 до 0,25. Вертикальные линии 1minx , 2minx соответствуют минимумам тока 1relJ , а вертикальная линия 1maxx – максимуму этого тока (координаты этих экстремумов определялись как по- лусумма координат слева и справа от экстремума, соответствующих одина- ковому значению тока 1relJ ). В этих точках ток 2relJ имеет следующие зна- чения:   1min2 xJ rel 2,60;   2min2 xJ rel 2,76;   1max2 xJ rel 1,56. Соответ- ствующие им значения  2sin  , найденные по формулам (17) и (18), равны 0,30; 0,38; 0,22. Их среднее значение равно 0,30, в то время как истинное зна- чение  2sin  для  0,2 равно 0,31. 146 0,5 1,0 1,5 2,0 2,5 3,0 3,5 0 1 2 3 4 x max1 x min2 x min1 J rel1 , J rel2 x, см J rel1 J rel2 Рис. 4 Точно так же определяется и фактическое расстояние 23l между зондами 2 и 3. При известных расстояниях 12l и 23l входящие в (2) и (3) углы  и  даются выражениями g   124 l ,   g   23124 ll . Следует отметить, что в полученные выше выражения для модуля и фа- зы комплексного коэффициента отражения не входит квадрат модуля ампли- туды напряженности электрического поля падающей электромагнитной вол- ны 2 падE . Следовательно, при проведении измерений амплитуда напряженно- сти электрического поля падающей электромагнитной волны не обязательно должна поддерживаться строго постоянной, что смягчает требования к гене- ратору зондирующего электромагнитного излучения в части стабильности мощности. Таким образом, предложен метод измерения комплексного коэффициен- та отражения в диапазоне частот с использованием трех зондов, аппаратная реализация которого значительно проще аппаратной реализации известных методов, что открывает путь к разработке нового класса измерителей ком- плексного коэффициента отражения. 1. Энген Г. Ф. Успехи в области СВЧ-измерений / Г. Ф. Энген // Труды Института инженеров по электро- технике и электронике. – 1978. – T. 66, № 4. – С. 8 – 20. 2. Опыт разработки ААЦ с двенадцатиполюсными рефлектометрами / А. А. Кабанов, С. М. Никулин, В. В. Петров, А. Н. Салов // Измерительная техника. – 1985. – № 10. – С. 38 – 40. 3. Cronson H. M. A 94-GHz diode-based single six-port reflectometer / H. M. Cronson, R. A. Fong-Tom // IEEE Transactions on Microwave Theory and Technology. – 1982. – V. MTT-30, No 8. – P. 1260 – 1264. 4. Барташевский Е. Л. Векторный СВЧ-рефлектометр на основе четырехплечего делителя мощности / Е. Л. Барташевский, В. А. Карлов // Электронная техника. Серия 1. Электроника СВЧ. – 1989. – № 1. – С. 38 – 44. 147 5. Векторный измеритель коэффициента отражения в 6-мм диапазоне длин волн / О. О. Дробахин, Н. Б. Горев, В. А. Карлов, И. Ф. Коджеспирова, Е. Н. Привалов // СВЧ-техника и телекоммуникацион- ные технологии : 13-я международная Крымская конференция, сентябрь, 2003 г., Севастополь : сб. ма- териалов. – Севастополь : Вебер, 2003. – С. 664 – 665. 6. Bates R. H. T. Image Restoration and Reconstruction / R. H. T. Bates, M. J. McDonnel. – Oxford : Clarendon Press, 1986. – 288 p. 7. Tikhonov A. N. Solution of Ill-Posed Problems / A. N. Tikhonov, V. Y. Arsenin. – New York : Winston-Wiley, 1977. – 258 p. 8. Афонин Ф. И. Измеритель комплексного коэффициента отражения на основе двухканального двухде- текторного волноводно-щелевого преобразователя / Ф. И. Афонин, Г. В. Боков, И. В. Лащенко // Збірник наукових праць Академії військово-морських сил імені П. С. Нахімова. – 2011. – № 3. – C. 52 – 58. 9. Двухзондовая реализация интерференционного метода измерения параметров движения механических объектов / О. В. Пилипенко, Н. Б. Горев, А. В. Доронин, И. Ф. Коджеспирова, Е. Н. Привалов // Техниче- ская механика. – 2013. – № 4. – С. 112 – 122. 10. Displacement measurement using a two-probe implementation of microwave interferometry / A. V. Doronin, N. B. Gorev, I. F. Kodzhespirova, E. N. Privalov // Progress in Electromagnetics Research C. – 2012. – V. 32. – P. 245 – 258. Институт технической механики Национальной Получено 06.10.2015, академии наук Украины и Государственного в окончательном варианте 08.10.2015 космического агентства Украины, Днепропетровск