Математическая модель гибридных систем магнитной левитации энерговырабатывающего оборудования автономных систем электроснабжения

УДК: 621.313

В. Е. Вавилов

Уфимский государственный авиационный технический университет

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ГИБРИДНЫХ СИСТЕМ МАГНИТНОЙ ЛЕВИТАЦИИ ЭНЕРГОВЫРАБАТЫВАЮЩЕГО ОБОРУДОВАНИЯ АВТОНОМНЫХ СИСТЕМ ЭЛЕКТРОСНАБЖЕНИЯ

Дата поступления: 30.08.2016 Решение о публикации: 01.09.2016 Дата публикации: 30.09.2016

Аннотация. Введение. Обычно при разработке математических моделей магнитолевитационных систем (магнитных подшипников, электростатических подшипниковых опор, гибридных магнитных подшипников (ГМП и т.д.), рассматривается непосредственно сама подшипниковая опора, как отдельная, изолированная электромеханическая система. Данный подход позволяет довольно точно исследовать процессы в магнитолевитационных системах, но практически не позволяет оценивать процессы, происходящие в системе магнитная опора-объект установки. Для решения данной проблемы автором предлагается иной подход к анализу процессов в магнитолевитационных системах в целом, и в ГМП в частности- рассмотрение системы магнитная опора-объект установки, как единого комплекса.

Цель. Работа направленна создание обобщенной аналитической модели высокооборотных ЭМПЭ с высококоэрцитивными постоянными магнитами (ВПМ) на упругих подшипниковых опорах, учитывающей взаимовлияние процессов в ЭМПЭ и подшипниковых опорах. Данная задача является новой и актуальной и имеет важное значение для современного энергомашиностроения. Для решения данной задачи в статье разрабатывается обобщенная математическая модель роторной системы на гибридном магнитном подвесе. Оценивается влияние гибридных магнитных подшипников на общее поведение роторной системы. Производиться анализ процессов в электромеханическом преобразователей энергии и механических процессов, протекающих в роторной системе.

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

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

Ключевые слова: гибридные магнитные системы левитации, электромеханические преобразователи энергии

V. E. Vavilov

Ufa State Aviation Technical University

MATHEMATICAL MODEL OF THE HYBRID SYSTEM OF MAGNETIC LEVITATION ENERGY PRODUCTION EQUIPMENT AUTONOMOUS POWER SUPPLY SYSTEMS

Abstract: Introduction. Typically, when developing mathematical models magnetocavitation systems (magnetic bearings, electrostatic bearings, hybrid magnetic bearings (GMP, etc.) is considered just the very bearing as a separate, isolated Electromechanical system. This approach allows us to accurately explore the process magnetocavitation systems, but practically does not allow to evaluate the processes occurring in the system of magnetic bearing-object position. To solve this problem, the author proposes a different approach to the analysis of the processes in magnetocavitation systems in General and GPC in particular considering the magnetic bearing-object position, as a single complex.

Goal. The work aimed the creation of a generalized analytical model of high-speed, AMPE with coercivity permanent magnet (VPM) on an elastic bearing supports, taking into account the mutual influence of processes in AMPA and bearing supports. This task is new and relevant and is essential to modern engineering. To solve this problem this paper developed a generalized mathematical model of the rotor system on a hybrid magnetic suspension. Evaluate the impact of hybrid magnetic bearings on the overall behavior of the rotor system. Performed analysis of processes in Electromechanical energy converters and mechanical processes occurring in the rotary system.

Method and methodology. The research methodology is based on the joint solution of Maxwell equations and equations describing the mechanical processes of a rotor system with five degrees of freedom.

Conclusion. The generalized mathematical model of high-speed, AMPE with VPM on a non-contact bearing supports and conducted her research. Based on research of the developed mathematical model, the authors developed an original control algorithm for the rotor position in a hybrid magnetic bearings, which allows for the design of high-speed, AMPE with VPM to abandon the position sensors of the rotor. In addition, on the basis of the results of calculations, a method was developed for diagnostics of eccentricity of rotor are high-speed, AMPE with VPM, as well as new methods of calculation of high-speed, AMPE with VPM, past experimental verification.

Key words: hybrid magnetic levitation system, Electromechanical converters of energy

Введение

Развитие авиакосмической и машиностроительной отрасли, робототехники и автономной энергетики ставит перед мировым электромашиностроением задачи по созданию высокооборотных электромеханических преобразователей энергии с высококоэрцитивными постоянными магнитами (ЭМПЭ с ВПМ) мощностью от 10 до 700 кВт, с частотами вращения ротора до 200 000 об/мин и сверхвысокооборотных

ЭМПЭ мощностью до 1 кВт и частотами вращения до 1 000 000 об/мин [1, 2].

Применение данных ЭМПЭ с ВПМ является одним из перспективных направлений, которое позволяет значительно расширить функциональные возможности объектов, в которых они установлены, а также минимизировать массогабаритные показатели этих объектов и повысить их вырабатываемую мощность. Поэтому различными производственными компаниями и корпорациями ведутся исследования, разработки и производство подобных ЭМПЭ [3-9].

Так компанией Calnetix разрабатываются высокооборотные ЭМПЭ с ВПМ мощностью от 3 до 700 кВт для авиационной промышленности, кораблестроения, нефтяной и газовой промышленности. В ЭМПЭ с ВПМ компанией Calnetix для минимизации потерь на трение, обусловленных высокой частотой вращения ротора ЭМПЭ используются магнитные подшипник.

Компанией Capstone Turbine Corporation серийно выпускается несколько типоразмеров микротурбинных установок в которых используются высокооборотные ЭМПЭ с ВПМ (частоты вращения от 60000 до 96000 об/мин, мощности от 35 до 1200 кВт). Для минимизации потерь на трение в ЭМПЭ с ВПМ компании Capstone используются аэродинамические подшипниковые опоры.

NASA (Glenn Research Center, Cleveland, Ohio) работает над созданием высокооборотного ЭМПЭ с ВПМ для малых космических аппаратов. Мощность данного ЭМПЭ будет составлять 7 кВт, частота вращения ротора 50000 об/мин. Для снижения потерь на трение применяются электромагнитные подшипник и магнитные подшипник на постоянных магнитах [3].

Аналогичные работы ведутся ABB, ETH Zurich, Cryostar Group, E+A Elektromaschinen und Antriebe AG и другими компаниями [3-7].

Из анализа промышленно производимых высокобортных ЭМПЭ с ВПМ [3-7], большинство данных ЭМПЭ, с целью повышения их эффективности и минимизации потерь выполняются на бесконтактных подшипниковых опорах (БПО, аэродинамических или магнитных подшипниковых опорах) или упругих подшипниковых опорах (например гидродинамических подшипниках).

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

инструментов для их разработки, которые позволили бы учесть взаимовлияния БПО и ЭМПЭ с ВПМ в процессе эксплуатации.

Постановка задач исследования

В работе [9] представлена математическая модель ЭМПЭ с ВПМ, а также разработана ее программная реализация в Matlab Simulink. В данной модели рассмотрены только электромагнитные процессы в высокооборотных ЭМПЭ. В [10] произведен многодисциплинарный анализ тепловых и механических процессов в высокооборотных ЭМПЭ с помощью программного комплекса Ansys. В [11] также с помощью программного комплекса Ansys решена задача многодисциплинарных исследований сверхвысокооборотных ЭМПЭ. В [12] разработана математическая модель для беспазовых ЭМПЭ в которой учитываются только электромагнитные процессы в высокооборотных ЭМПЭ. В [13] решается многодисциплинарная задача оценки теплового состояния высокооборотных ЭМПЭ. В [14] приведены совместные исследования тепловых и электромагнитных процессов в ЭМПЭ. В [15] разработана математическая модель магнитных подшипников, но при этом не рассматривается их влияние на высокооборотный ЭМПЭ. Таким образом, математических моделей и междисциплинарных исследований высокооборотных отдельно ЭМПЭ с ВПМ и отдельно БПО в литературе представлено значительное множество, а многодисциплинарных исследований системы высокооборотный ЭМПЭ - бесконтактные подшипниковые опоры, с учетом их взаимовлияний в литературе не представлено. Хотя данная задача является весьма актуальной, как для фундаментальной науки (на ее основе может сформулирована обобщенная целевая функция для оптимизации высокооборотных ЭМПЭ с ВПМ), так и для прикладных задач (решение данной задачи позволит значительно повысить точность проектирования высокооборотных ЭМПЭ с ВПМ, по результатам решения данной задачи могут быть сформулированы новые алгоритмы управления положением ротора высокооборотных ЭМПЭ с ВПМ).

Поэтому данная работа направленна создание обобщенной аналитической модели высокооборотных ЭМПЭ с ВПМ на упругих подшипниковых опорах, учитывающей взаимовлияние процессов в ЭМПЭ и БПО. Данная задача, как видно из представленного обзора, является новой и актуальной и имеет важное значение для современного энергомашиностроения. При аналитическое решение данной задачи будет иметь общий, масштабный характер, в отличи от решения данной задачи методами компьютерного моделирования (Femm, Ansys).

Обобщенная математическая модель высокооборотного ЭМПЭ на упругих подшипниковых опорах

При разработке обобщенной математической модели, описывающей высокооборотные ЭМПЭ на БПО, рассматривается жесткий, недеформируемый ротор на упругих подшипниковых опорах (в качестве упругих подшипниковых опор могут выступать активные магнитные подшипники, газовые подшипники, гибридные магнитные подшипник), при этом ротор массой тр установлен на двух опорах (рисунок 1). Статическая неуравновешенность массы ротора характеризуется эксцентриситетом е, динамическая - угловым параметром р.

В качестве ЭМПЭ при разработке математической модели рассматривается неявнополюсная ЭМПЭ с ВПМ, рисунок 1 Для упрощения математических преобразований используются следующие допущения:

- магнитная проницаемость окружающей среды и воздушного зазора равна магнитной проницаемости вакуума;

- величина немагнитного зазора намного меньше радиуса его кривизны, кривизна;

- обмотка статора представляется в виде тонкого электропроводящего слоя, распределенного по диаметру расточки магнитопровода статора; ;

- плотность индуцированных токов по толщине тонкого медного слоя постоянна;

- теплопроводность активных элементов ЭМПЭ постоянна по всем координатным осям (активные элементы ЭМПЭ изотропны);

- аксиальная составляющая напряженности магнитного поля в торцевых поверхностях ротора равна 0, то есть рассматривается ЭМПЭ бесконечной длины;

- температура в воздушном зазоре считается известной.

При решении поставленной задачи ищется совместное решение нескольких систем уравнений:

- системы уравнений, описывающих электромагнитные процессы в высокооборотном ЭМПЭ (системой уравнений Максвелла для медленно движущихся сред):

Ы '

7 = сз\Е + [УхВ)\

д\\] = 0, Н = \10В,

<НуЯ = ре

где В- вектор магнитной индукции результирующего магнитного поля; Е, Н - вектор напряженности электрического и магнитного полей; V - вектор скорости движения ротора; а - электрическая проводимость обмотки статора; у - вектор плотности индуцированных токов; /ст - вектор плотности сторонних токов;

Рис. 1. Обобщенная конструктивная схема ЭМПЭ на ГМП

- математической модели, описывающей механическое движение ротора:

m1x1 - m12 x2 + h0

m1y1 - m12y2 - hü

- m12Х=1 + m2x2 - h0

- m12 y1 - m2y2 + h0

где Q1m -Q4m - силы реакции опор; Q1 - Q4 - внешние силы, определяющие характер движения ротора;

2 9

mpeю (J1 - J3)ero2/. Qlv (t) = -(COS ЮГ) - ( 1 j3)-(Sin Юt) ;

2 9

трею (J - J3)ero2, Q2 v (t) = —p-(z2 sin rot) - -^-(- cos rot);

1в 1в

mnea

2 ,тт\„2

Q3v (t) = (-zj cos at) - J—J3 )gL° (-sin rot);

2 2 q4v (t) = (-z sin at) -(J- J3)ea (cos at); e - эксцентриситет ротора; тъ-

lB lB

масса вала с ротором; 4 - длина вала ротора; a - частота вращения ротора;

J1 - экваториальный момент инерции вала; J3 - осевой момент инерции

2 2 J + mz2 J + mz z9 J + mzi

вала; m1 = ——т-^-2, m12 = —-r 1 2 , m2 = ——r-^-1— изменение массы

l l l

1в 4в 4в

ротора; к0 = -3— гироскопический коэффициент; х1? х2, у1, у2 - радиальные

I 2

смещения вала в упругой подшипниковой опоре; ^, - осевая координата первого и второго БПО; Х1? Х2, у1? у2 - вторые производные от перемещений по времени.

Для общего анализа физических процессов в ЭМПЭ с ВПМ разрабатываемая математическая модель должна учитывать зависимость параметров магнитного поля в рабочем зазоре ЭМПЭ от температуры.

Учитывая, что магнитная индукция в воздушном зазоре ЭМПЭ пропорциональна остаточной индукции ВПМ, то зависимость магнитного поля ЭМПЭ с ВПМ от температуры может быть описана в виде:

B(0) = B011 -

kBr

100 )

где 5(0) - действующие значения остаточной индукции и коэрцитивной силы ВПМ соответственно; 0ВПМ - температура ВПМ; квг-температурный коэффициент остаточной индукции ВПМ.

Данное выражение может быть представлено в следующем виде:

Впкя т

дБ - Б -

'0'* Br ВПМ

100

Разделив обе части выражения на dt:

дБ _ Bo Bo0kBr д0впм

& & 100 дt Из системы уравнений Максвелла следует, что

- = го1£.

дt

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

Н = Й1+Й2,

В = В1 + В2.

Тогда система уравнений Максвелла может быть переписана в виде:

^ . . дВ

гоШ = ] + 7ст + —:

дг

гоЕ = —

дВ дВ

2

дг дг

Г Во

ч дг

Во квг д© 100

У = °

ВПМ

дг .

сНУ5 = 0, сПУ/ = 0,

ё^О = р е,

Из выражения (15) видно, что если

д© ВПМ

дВ1 = 0, дВ2 = о

(ротор

неподвижен), но при этом

дг

дг дг

^ 0 (то есть температура ЭМПЭ

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

Данный вывод является одним из свойств ЭМПЭ с ВПМ, которое нехарактерно для других типов электрических машин.

Решая систему уравнений (1)—(11) относительно напряженности вторичного поля, имеем следующее:

ЛЯ.

Ш2 ёг

той

При решении уравнения (21) обычно принимается, что ротор ЭМПЭ обладает одной степенью свободы, и скорость ротора не зависит от координат. При рассмотрении ЭМПЭ на упругих подшипниковых опорах ротор обладает 5 степенями свободы, в связи с чем, данное допущение при решении для ЭМПЭ на упругих подшипниковых опорах будет не верным.

Представим гов виде:

к*(Гх Н2) = (Й2  -(V У)//2 + ШуЙ2 - Н2(ИУ¥ . п*(Кх Й,) = • У)К - {V  + ШуНх - Й^хУ¥ ,

где V- оператор Гамильтона. Так как &\уЙ2 = 0, divЙ1 = 0:

тог(? хя2) = (я2-у)г-(к-у)я2- Й2&\?,

гс*(г X нх) = (я2  у)Г - (V  у)я2 - Й2$ыУ.

АЯ2-ц0о

ая, а/

-

^-{{й^у-^-^щ-й^у)

Из векторного анализа известно, что:

/ - \- д¥ д¥ д¥

дг

дх у ду 2 дг При проектировании ЭМПЭ с упругими опорами более целесообразно оперировать переменными в декартовой системе координат. Это обусловлено тем, что практически все системы управления БПО построены именно на использовании декартовой системы координат.

Тогда с учетом выражений (27), (28) можно записать следующее:

АН2-[10а

а/

/

Н4\уУ +

дУ

дУ

дУ

V

дх ду дг

У

/г/ сН2 г/ сН2 г/ с//2ЛЛЛ Ух —2 + V —2 + У2 —2

дх у ду дг

-ц 0°

а/

V

ая

'и + н 89 + н

дх ду дг

.

V

/

дх у ду

2 г / дН2

+ V,—-

дг

\\

)))

Допускается, что в направлении оси г ротор жестко закреплен:

А#2-ц0ст

'6Н2 '

а/

-Н2дл\У -

ГнЛ + н.ЛЛ Г

дх

у 2

= -Ц 0°

^ г

а/

-НЛ\У-

дУ_ дх

Н

у1

дУ_

ду

ду

Н

т. дН2 дН2 К—2 + V —-

х дх у ду

дУ_

\ (

IV) дЙ,

дН,

\\\

V*— + ку— + К—

дх ду дг

Совокупность представленных выражений является обобщенным математическим описанием высокооборотного ЭМПЭ с ВПМ на БПО.

Заключение

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

Работа выполнена при поддержки Российского Фонда Фундаментальных Исследований (проект 16-38-60001).

Библиографический список

1. Hendershot J. & Miller T. Design of brushless permanent-magnet motors // Oxford University Press, 1996.

2. Borisavljeviс A. Limits, Modeling and Design of High-Speed Permanent Magnet Machines // Printed by Wormann Print Service. - Zutphen, the Netherlands, 2011. - P. 209.

3. Nagorny A. S., Jansen R. H. & Kankam D. M., Experimental performance evaluation of a highspeed permanent magnet synchronous motor and drive for a flywheel application at different frequencies, in Proceedings // Electrical Machines-ICEM, 2006. - pp. 1-8.

4. Binder A. & Schneider T. Permanent magnet synchronous generators for regenerative energy conversion - a survey // Power Electronics and Applications, 2005. - pp. 110-115.

5. Malmquist A. Mikrogasturbinen als Wegbereiter der dezentralen Wärme- und Stromversorgung // ABB Technik.- 3/2000. - p.22-30.

6. Nagorny A., Dravid N., Jansen R. & Kenny B. Design Aspects of a High Speed Permanent Magnet Synchronous Motor/Generator for Flywheel Applications // NASA/TM-2005-213651, 2005. - pp.1-7.

7. Bailey C., Saban D. & Guedes-Pinto P. Design of High-Speed Direct-Connected Permanent-Magnet Motors and Generators for the Petrochemical Industry // IEEE Transactions on Industry Applications, 2009. - vol. 45. № 3. -pp. 1159-1165.

8. Isomura K., Murayama M., Teramoto S., Hikichi K., Endo Y., Togo S. & Tanaka S. Experimental Verification of the Feasibility of a 100W Class Micro-scale Gas Turbine at an Impeller Diameter of 10 mm, J. // Micromech. Microeng, 2006. - no. 16. - pp. 254-261.

9. Kulkarni S. S. & Thosar A. G. Mathematical Modeling and Simulation of Permanent Magnet Synchronous Machine // International Journal of Electronics and Electrical Engineering. - vol. 1. - №. 2. - pp. 66-71.

10. Kolondzovski Z., Arkkio A., Larjola J. & Sallinen P. Power limits of high-speed permanent-magnet electrical machines for compressor applications. Report Series on Electromechanics // Report 76, Aalto University School of Science and Technology, Espoo 2010. - 21 p.

11. Li W., Zhang X. & Cheng S. Thermal Optimization for a HSPMG Used for Distributed Generation Systems // IEEE Transactions on Industrial Electronics, 2013. - vol. 60. - no. 2. - pp. 474-482.

12. Pierre-Daniel Pfister Very High Speed Slotless Permanent Magnet Motors: Analytical Modeling, Optimization, Design and Torque Measurement Methods // IEEE Transactions on industrial electronics. - vol. 57. - no. 1.

- pp. 45-49.

13. Tuysuz A., Steichen M., Zwyssig C. & Kolar J. W. Advanced cooling concepts for ultra-high-speed machines // Power Electronics and ECCE Asia (ICPE-ECCE Asia), pp. 2194-2202.

14. Vong P. & Rodger D. Coupled electromagnetic-thermal modeling of electrical machines // IEEE Trans. on Magn. - vol. 39. - no 3. - pp. 1614-1617.

15. Shen J. & Fabien B. Optimal control of a flywheel energy storage system with a radial flux Hybrid Magnetic Bearing // Journal of the Franklin Institute-Engineering and Applied Mathematics, 2002. - vol. 339. - no 2.

- pp. 189-210.

References

1. Hendershot J. & Miller T. Design of brushless permanent-magnet motors, Oxford University Press, 1996.

2. Borisavljevic A. Limits, Modeling and Design of High-Speed Permanent Magnet Machines. Zutphen, 2011. 209 p.

3. Nagorny A. S., Jansen R. H. & Kankam D. M. Electrical Machines-ICEM, 2006, pp. 1-8.

4. Binder A. & Schneider T. Power Electronics and Applications, 2005, pp. 110-115.

5. Malmquist A. ABB Technik, 2000, no. 3, pp.22-30.

6. Nagorny A., Dravid N., Jansen R. & Kenny B. NASA/TM,, 2005, 213651, pp.1-7.

7. Bailey C., Saban D. & Guedes-Pinto P. IEEE Transactions on Industry Applications, 2009, vol. 45, no. 3, pp. 1159-1165.

8. Isomura K., Murayama M., Teramoto S., Hikichi K., Endo Y., Togo S. & Tanaka S. Micromech. Microeng, 2006, no.16, pp. 254-261.

9. Kulkarni S. S. & Thosar A. G. International Journal of Electronics and Electrical Engineering, vol. 1, №. 2, pp. 66-71.

10. Kolondzovski Z., Arkkio A., Larjola J. & Sallinen P. Power limits of high-speed permanent-magnet electrical machines for compressor applications. Report Series on Electromechanics, Report 1976, Aalto University School of Science and Technology, Espoo 2010, 21 p.

11. Li W., Zhang Х. & Cheng S. IEEE Transactions on Industrial Electronics, 2013, vol. 60, no. 2, pp. 474-482.

12. Pierre-Daniel Pfister IEEE Transactions on industrial electronics, vol. 57, №. 1, pp.45-49

13. Tuysuz A., Steichen M., Zwyssig C. & Kolar J. W. Power Electronics andECCE Asia (ICPE-ECCE Asia), pp. 2194-2202.

14. Vong P. & Rodger D. IEEE Trans. on Magn., vol. 39, no. 3, pp. 1614-1617.

15. Shen J. & Fabien B. Journal of the Franklin Institute-Engineering and Applied Mathematics, 2002, vol. 339, no. 2, pp. 189-210.

Сведения об авторе:

ВАВИЛОВ Вячеслав Евгеньевич E-mail: s2_88@mail.ru Information about author:

Viacheslav VAVILOV E-mail: s2_88@mail.ru

© ВАВИЛОВ В. Е., 2016