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

-►

ЭЛЕКТРОТЕХНИКА

DOI: 10.18721/ JEST.230109 УДК 621.315.176

Ю.Н. Бочаров, В.В. Титков, Р.Ш. Абитаева, А.Б. Бекбаев, М.А. Джаманбаев

СВОБОДНОЕ КРУТИЛЬНОЕ КОЛЕБАНИЕ РАСЩЕПЛЕННОГО ПРОВОДА

ЛИНИЙ ЭЛЕКТРОПЕРЕДАЧ

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

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

Ссылка при цитировании:

Ю.Н. Бочаров, В.В. Титков, Р.Ш. Абитаева, А.Б. Бекбаев, М.А. Джаманбаев. Свободное крутильное колебание расщепленного провода линий электропередач // Научно-технические ведомости СПбГПУ. 2017. Т. 23. № 1. С. 90-97. DOI: 10.18721/ JEST.230109

Y.N. Bocharov, V.V. Titkov, R.S. Abitaeva, A.B. Bekbaev, M.A. Jamambaev

TORSION VIBRATIONS OF SPLITTED WIRE OF OVERHEAD TRANSMISSIONS

Free torsional oscillation of a split wire of an overhead transmission line is considered. Torsional oscillations of a split wire were calculated based on a non-linear differential equation obtained by the authors using the Bubnov-Galerkin method. Calculated formulas for torsional stiffness of a split wire have been derived. Ratios between the characteristics of torsional and lateral oscillations of a split wire were obtained, which allow optimal selection of parameters of the devices for vibration and oscillation dampening. The results of calculations agree with the experimental data as well as known the theoretical results. Obtained calculation formulas are recommended for use in design of new overhead transmission lines with split wires and reconstruction of lines already in operation.

OVERHEAD TRANSMISSIONS LINES; SPLITTING OF WIRE; TORSION MOTION EQUATION; TORSION HARDNESS; EXPERIMENTAL DATA AND CALCULATIONS COMPARISON.

Citation:

Y.N. Bocharov, V.V. Titkov, R.S. Abitaeva, A.B. Bekbaev, M.A. Jamambaev, Torsion vibrations of splitted wire of overhead transmissions, St. Petersburg polytechnic university journal of engineering sciences and technology, 23 (1) (2017) 90-97, DOI: 10.18721/ JEST.230109

На воздушных линиях электропередачи сверхвысокого напряжения, как правило, применяются расщепленные провода (РП), позволяющие существенно снизить потери энергии, вызываемые коронным разрядом [1]. Конструкция расщепленного провода представляет собой совокупность одиночных проводников, центры сечения которых размещаются на окружности радиуса Я (рис. 1). При таком расположении универсальными геометрическими характеристиками расщепленного провода становятся радиус кругового сечения одиночного проводника (г0), упомянутый выше радиус расщепления Я и количество одиночных проводов в фазе — п [1, 3].

Одной из проблем эксплуатации воздушных линий электропередачи высокого напряжения является обеспечение минимально допустимых расстояний между фазными проводами, а также между фазными проводами и элементами конструкции опор (стойки, траверсы) [2]. В условиях ветровых нагрузок наиболее опасный режим — это так называемая пляска проводов, сопровождающаяся их масштабными взаимными перемещениями [3].

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

Синхронизированные вертикальные и крутильные колебания расщепленного провода вызывают периодические изменения угла атаки воздушного потока, воздействующего на провод [3, 4, 6]. Для линий с расщепленными проводами крутильные колебания проявляются в периодических закручиваниях фазы в целом.

Как известно, одно из предложений по борьбе с пляской проводов заключается в активном вмешательстве в этот процесс. Один из подходов основывается на достижении максимальных различий в частотах поперечных и крутильных колебаний, что повышает устойчивость РП к пляске. Широко используются гасители, принцип действия которых направлен на принудительную расстройку частот поперечных и крутильных движений РП (например, маятниковые гасители) [3, 5]. Для более эффективной реализации указанных мер интересно проанализировать крутильные движе-

Рис. 1. Конструкция расщепленного провода

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

Кроме того, при расчете параметров некоторых типов гасителей, воздействующих на крутильные движения РП, необходимо иметь сведения о его крутильной жесткости. Сведения о крутильной жесткости полезно знать и при моделировании пляски РП с помощью эквивалентного провода, эквивалентных пружинных устройств и т. д. Для оценки крутильной жесткости применяются аналитические [7] и экспериментальные методы. Последние были, в частности, реализованы на опытном полигоне Казахстанского научно-исследовательского института энергетики.

В задачу данной работы входит определение частотных характеристик свободных колебаний и крутильной жесткости РП.

Уравнение крутильного движения РП

Уравнение крутильного движения РП определяем исходя из уравнения Лагранжа [8]:

йг

д1А дф (г X

дф (г)

= о,

(1)

где ЬА = Ек - Ей — функция Лагранжа; Ек—кинетическая энергия РП; — энергия деформаций проводов расщепленной фазы.

При аппроксимации форм колебаний расщепленного провода функциями вида

П7

Ф( г,г) = ф(г )8ш—

(2)

кинетическая энергия от вращательного движения РП равна

i Т

Ек =í

дФ( z,t) dt

dz =

. nPplR2 .2

4 g

Ф 2(t),

(3)

где — обобщенная координата; г — расстояние от опоры до произвольного сечения РП; I — длина пролета; Р0 — масса одного метра провода; п — число расщеплений (число проводов в фазе); Я — радиус расщепления; /ф — момент инерции РП,

Jф = P R2.

(4)

Edi = T0 L0 ) + ^L0 ) ,

2i

(5)

Lo =J

1+1 fdy( z) dz

dz;

Lp=J

1 +1 2

дУф( z,t) dz

dz,

(6)

y( z) = A- z(i - z).

21a

(7)

Уф (z,t) =

/X 1 Л • 2/^-2 nz = У (z) + 2Rsin Ц; Ф (t)sin -y -

-R cos ц; ф(t)sin у.

(8)

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

Ц; = Ц +

2я(/ -1)

(9)

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

где I = '1 - п.

Например, для РП, состоящего из трех проводов (п = 3), согласно (9) имеем: ^ = 30° ;

=150°; ц3 = 270° .

Длины проводов РП до (Ь0) и после (Ьф) закручивания согласно формулам (6), (7), (8) равны следующим:

где Е — модуль Юнга; Г — площадь поперечного сечения провода; Т0 — начальное натяжение провода.

Длину провода Ь0 в положении статического равновесия и длину Ьф, соответствующую закрученному состоянию г-го провода РП, определяем по приближенным формулам:

P 2i3

Lo = i + р0±~1 24Т2

(10)

P i2

L= L0 + Ф 0 2nT0

4

—j

n

k — k2 | +

i

+—

4

(( + k2)-ü k,k2;

(11)

где у( г) и уф (г, ^) — функции, описывающие конфигурацию (провес) провода в пролете до и после закручивания РП.

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

где k = у Rsin Ц; ф2^); k2 = ПRcos ц ф(t).

Подставляя разности длин проводов (Lф - L0)

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

Ed=X Edi =

n2R 2T0

4i

1 +

8EFP02 i2 ^ n t0 j

X cos2 Ц; Ф2(t) + 1

После геометрических расчетов кручения РП вдоль пролета определено выражение для функций у^г,г):

n2R 2T0

16i

EFP02i2 ^ 2n2t03 j

о

и

о

о

<Ё 81П2 Иг +

п4Я 4 ЕЕ Л 4

-5—> ео8 и,

32£3 1 г

ф4(г)+

к4Я 4 ЕЕ 64£3

X 81и2и,- ео82 и, ф6(г) +

я4Д4ЕЕ п

512£3

х sin4 и, ф8(г).

(12)

О2 s2 + П- Лл3

4

\ 3_4

ф3(г)+ф5(г)+ 16

+^ DS4 ф7(г) = о,

где Л = ^; 02 = п2 .ТО

РоГ

2Ро г

2Й2\

1 +

ЕЕРо2£

2п ТО у

(13)

; юл =

2

Ро £

1 +

8ЕЕРо2£2 ^

_4т" 3

П То У

В табл. 1 приведены отличные от нуля значения тригонометрических коэффициентов.

Решение нелинейного дифференциального уравнение движения

Решение уравнения (13) осуществляется приближенным методом Бубнова — Галеркина[8]. Вариационная форма уравнения имеет вид

2п

По выражениям (3) и (12) определяем функцию Лагранжа Ьл = Ек - Ей . Подставляя Ьл в уравнение (1), получим дифференциальное уравнение крутильного колебания:

ф (г)+®л2у1 ф(г)+

4

ф(г)+ ®л 51ф(г)+

О2л2 + П- Лл3

\

У

+^Пт °*5ф5(г)ш4ф7(г)

16 32

ф3(г)+ Ьфйг = о.

Положим

ф(г) = фоео8ш^; 5ф = 5фоео8ш^,

где ф 0 — амплитуда крутильного колебания РП, шк — искомая частота крутильного колебания.

Интегрируя вариационное уравнение с учетом 5ф, получим частотное уравнение

так = ®к

15п4

1+

3

4®|

О2 л2 + П- Лл3

\

фо +

у

128®;

"Л^5фо +

35п4

2о48®к

Л^фо

(14)

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

поперечного колебания линеаризованной системы [10].

®к

= ®л>Й".

(15)

Таблица 1

Числовые значения тригонометрических коэффициентов

ш

1

о

Число расщепления п

Xи,

Xsin2 и,

Xеos4 и,

Xsin4 и,

Е. 2 2

sin и, еos и,

1

0,5 0,5 0,5 0,5 0,5

1

0,5 0,5 0,5 0,5 0,5

1

0,37 0,25 0,37 0,37 0,37

0 0,37 0,25 0,37 0,37 0,37

0 0,12 0,25 0,12 0,12 0,12

л =

л =

2

3

4

5

п

п

п

п

п

Как следует из формулы (15), отношение собственных частот поперечного и крутильного движения РП — величина постоянная и составляет

* =72.

Ок л/51

(16)

ОЬ

= Юг,

1 +

3

2

4«2

О2 + Д?3

V

Фо;

^^ 1 +

3О2 52 4«2

(17)

Согласно [8] крутильная частота вала определяется по формуле

Р =

(19)

Формула (16) не учитывает влияния амплитуды поперечного и крутильного колебания на частотные характеристики. Частота поперечного колебания практически не изменяется с увеличением амплитуды [10], а частота кручения чувствительна к изменению угла закручивания (рис. 1, а). Следовательно, с увеличением амплитуды кручения отношение (16) несколько уменьшается. Для определения отношения частот с учетом влияния амплитуды крутильных колебаний ограничимся двухчленным приближением формулы (14) (для реальных линий электропередачи коэффициентом Б можно пренебречь):

где р — крутильная частота вала; К— крутильная жесткость вала; / — момент инерции вала. Заменим РП эквивалентным валом:

пЕ\

р=Ок; К = крф; 1 = 1рф = Я2. Решив (19) относительно искомой величины

Крф, находим

=про 2 у

К РФ" 2 к

п g

пР0£2 Я2 Г 2 3 2п g V 4 .

1 +

3О2 ^ л

4О

к У

(20)

Для РП из двух горизонтально расположенных проводов (п = 2) выражение (20) несколько упрощается:

К_=п^ о2 = п^ . (21)

рф

П2 g

2п2 g

В формуле (17) при упрощении значение ф0 принято равным 1. Используя формулу (17), находим (с учетом = S2 = 0,5)

О

со к

8юд

[4ю2 + 3О2

(18)

Крутильная жесткость расщепленной фазы

Крутильная жесткость РП может быть определена исходя из зависимости между жесткостью и крутильной частотой эквивалентного вала.

Результаты расчетов

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

Экспериментальные данные по частотным характеристикам, а также по крутильной жесткости РП получены на опытном полигоне Каз-НИИ энергетики имени академика Ш.Ч. Чоки-на. Все расчетные и экспериментальные данные относятся к проводу марки АС0-330/39 сечением 330 мм2.

Vk, Гц

0,39 0,38 0,37 0,36 0,35 0,34 0,33

п = 3

1 = 200 м

0

20

40

60

Ф0

'а / Vк 1,2 1

0,8 0,6 0,4 0,2 0

0

^—♦—а-

С - * - *-1

▼

п = 5

1 = 288 м

10 а0, даН/мм2

Рис. 2. Зависимость частоты крутильного колебания РП от амплитуды кручения (а); отношение частоты поперечного и крутильного колебания РП (б)

V*, Гц 0,4

0,3 0,2 0,1 0

0

♦ ♦♦

п = 5

/ = 288 м

ю а0, даН/мм2

"^ь ГЦ 0,4

0,3

0,2

0,1

0

♦

п = 8

/ = 354 м

О

а0, даН/мм2

Рис. 3. Зависимости частоты крутильных колебаний РП от механического напряжения провода: а — РП из пяти проводов, / = 288 м; б — РП из восьми проводов, / = 354 м

На рис. 2,а приведены зависимости частоты крутильного колебания от амплитуды кручения. Расчет выполнен при растягивающем напряжении провода а0, равным 8 даН/мм2. Число расщепления — п = 3. Длина пролета — / = 200 м. Как следует из рис. 2,а, при амплитуде кручения Ф0, равной 60° (характеризует максимальный угол поворота РП при пляске), частота крутильного колебания может увеличиваться до 17 %. С увеличением натяжения провода и длины пролета это отличие несколько сглаживается.

На рис. 2,б приведены результаты сопоставления экспериментальных и расчетных данных об отношении частот поперечных и крутильных колебаний РП при различных а0. Длина опытного пролета — / = 288 м, число расщепления п = 5. Как видно из рисунка, при малых напряжениях между

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

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

В табл. 2 приведены результаты сопоставления экспериментальных и теоретических данных по

Таблица 2

Сопоставление экспериментальных и теоретических данных крутильной жесткости РП

а0, даН/мм2 Крф, даН-м2/рад

Эксперимент Теоретич. по ф-ле (20) Теоретич. по ф-ле (22)

РП с 4,5 5.4 8.5 13,2 п = 5;/ = 288 м 2659 2426 2493 3177 4087 (53 %) 3145 (29 %) 2123 (17 %) 2186 (45 %) 1323 (100 %) 1008 (240 %) 655 (380 %) 651 (488 %)

РП с 5,7 6,6 7,1 8,4 п = 8;/= 354м 10365 9371 9467 9862 15144 (46 %) 12482 (33 %) 11491 (21 %) 9890 (0,2 %) 3049 (340 %) 2490 (376 %) 2280 (415 %) 1934 (510 %)

крутильной жесткости РП. Теоретические расчеты выполнены по различным формулам. В скобке приведены отличия между экспериментальными и теоретическими данными в процентах.

Согласно [7] для определения эквивалентной крутильной жесткости РП принимается формула

Крф = m+T0R

1 + EFP02l

2,2 Л

12T3

(22)

V ~ 0 у

где т — крутильная жесткость одиночного провода (при расчетах можно пренебречь, посколькУ т << Крф ).

Как показывает сравнительный анализ, расчетные формулы, полученные в [7], дают для крутильной жесткости существенно заниженные значения.

Выводы

1. Частота крутильного колебания расщепленного провода всегда меньше частоты поперечного колебания. При минимальном расщеплении (число проводов в фазе п = 2) частоты

крутильных и поперечных колебаний совпадают. При увеличении числа расщепления (п>2), применяемом в России для линий классом выше 330 кВ, частота крутильных колебаний всегда ниже, чем поперечных.

2. С увеличением натяжения и с уменьшением длины пролета разности между частотой поперечного и крутильного колебания уменьшается.

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

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

5. Сопоставление расчетных и экспериментальных данных подтверждает надежность полученных формул. Они могут быть использованы при решении прикладных задач, связанных с ограничением пляски проводов.

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

1. Александров Г.Н. Передача электрической энергии. СПб.: Изд-во Политехнического университета, 2007.

2. Кучинский Г.С., Кизеветтер В.Е. Пинталь Ю.С.

Изоляция установок высокого напряжения. М.: Энер-гоатомиздат, 1987.

3. Электрические сети сверх- и ультравысокого напряжения. Теоретические и практические основы. Т. 1 / Ред. А.Ф. Дьяков. М.: НТФ «Энергопрогресс» Корпорации «ЕЭЭК», 2012.

4. Пустыльников Л.Д., Шкапцов В.А. Аэродинамически неустойчивые колебания проводов воздушных линий электропередачи с гололедными отложениями // Известия АН СССР. Энергетика и транспорт. 1991. № 2. С. 103-106

5. Крылов С.В., Шкапцов В.А. Улучшенная система подвески проводов для больших воздушных переходов с промежуточными опорами // Электрические станции. 1999. №3. С. 36-42.

6. Бекметьев Р.М., Жакаев А.Ш., Ширинских Н.В.

Пляска проводов воздушных линий электропередачи. Алма-Ата: Наука КазССР, 1979.

7. Wang I., Lilien J. L. Overhead electrical transmission line galloping. AfUllmulti-Span 3—DOF—Model, some Application and design recommendations //IEEE Transactions on Power Delivery. 1998. Vol. 13. №3. P. 909-916.

8. Светлицкий В.А., Стасенко И.В. Сборник задач по теории колебаний. М.: Высшая школа, 1973.

9. Бекметьев Р.М., Джаманбаев М.А. Методика расчета динамических нагрузок при пляске проводов // Сборник докладов советских специалистов на международном совещании по проблемам пляски проводов ЛЭП. Сочи. 1985. С. 56-68.

10. Бекбаев А.Б., Джаманбаев МА., Акпанбетов Д.Б., Токенов Н.П. Исследование влияние амплитуды колебаний проводов линий электропередач на собственную частоту // Вестник Казахстанской Национальной академии естественных наук. 2012. №1. С. 64-66.

СВЕДЕНИЯ ОБ АВТОРАХ

БОЧАРОВ Юрий Николаевич — доктор технических наук профессор Санкт-Петербургского политехнического университета Петра Великого. 195251, Россия, г. Санкт-Петербург, Политехническая ул., 29. E-mail: [email protected]

ТИТКОВ Василий Васильевич — доктор технических наук профессор, заведующий кафедрой Санкт-Петербургского политехнического университета Петра Великого. 195251, Россия, г. Санкт-Петербург, Политехническая ул., 29. E-mail: [email protected]

АБИТАЕВА Рахимаш Шанракбаевна — ассистент Казахского национального технического университета (КазНТУ). 050013, Казахстан, г. Алматы, ул. Сатпаева, 22а. E-mail: [email protected] БЕКБАЕВ Аманкельды Бекбаевич — доктор технических наук профессор Казахского национального технического университета имени К.И. Сатпаева. 050013, Казахстан, г. Алматы, ул. Сатпаева, 22а. E-mail: [email protected]

ДЖАМАНБАЕВ Мураткали Абенович — кандидат технических наук (PhD) доцент Казахского национального технического университета имени К.И. Сатпаева. 050013, Казахстан, г. Алматы, ул. Сатпаева, 22а. E-mail: r [email protected]

REFERENCES

1. Aleksandrov G.N. Peredacha elektricheskoy ener-gii. SPb.: Izd-vo Politekhnicheskogo universiteta, 2007. (rus.)

2. Kuchinskiy G.S., Kizevetter V.Ye. Pintal Yu.S. Izo-lyatsiya ustanovok vysokogo napryazheniya. M.: Ener-goatomizdat, 1987. (rus.)

3. Elektricheskiye seti sverkh - i ultravysokogo napryazheniya. Teoreticheskiye i prakticheskiye osnovy. T. 1 / Red. A.F. Dyakov. M.: NTF «Energoprogress» Korporat-sii «YeEEK», 2012.

4. Pustylnikov L.D., Shkaptsov V.A. Aerodinami-cheski neustoychivyye kolebaniya provodov vozdushnykh liniy elektroperedachi s gololednymi otlozheniyami. Izves-tia ANSSSR Energetika i transport. 1991. № 2. S. 103-106. (rus.)

5. Krylov S.V., Shkaptsov V.A. Uluchshennaya sistema podveski provodov dlya bolshikh vozdushnykh perek-hodov s promezhutochnymi oporami. Elektricheskiye stantsii. 1999. №3. S. 36-42. (rus.)

6. Bekmetyev R.M., Zhakayev A.Sh., Shirinskikh N.V.

Plyaska provodov vozdushnykh liniy elektroperedachi. Alma-Ata: Nauka KazSSR, 1979.

7. Wang I., Lilien J. L. Overhead electrical transmission line galloping. AfUllmulti-Span 3—DOF—Model, some Application and design recommendations . IEEE Transactions on Power Delivery. 1998. Vol. 13. №3. P. 909-916.

8. Svetlitskiy V.A., Stasenko I.V. Sbornik zadach po teorii kolebaniy. M.: Vysshaya shkola, 1973.

9. Bekmetyev R.M., Dzhamanbayev M.A. Metodika rascheta dinamicheskikh nagruzok pri plyaske provodov. Sbornik dokladov sovetskikh spetsialistov na mezhdunarod-nom soveshchanii po problemam plyaski provodov LEP. Sochi. 1985. S. 56-68. (rus.)

10. Bekbayev A.B., Dzhamanbayev M.A., Akpanbetov D.B., Tokenov N.P. Issledovaniye vliyaniye amplitudy kolebaniy provodov liniy elektroperedach na sobstvennuyu chastotu. Vestnik Kazakhstanskoy Natsionalnoy akademii yestestvennykh nauk. 2012. №1. S. 64-66. (rus.)

AUTHORS

BOCHAROV Yurii N. — Peter the Great St. Petersburg Polytechnic University. 29 Politechnicheskaya St., St. Petersburg, 195251, Russia. E-mail: [email protected]

TITKOV Vasilii V. — Peter the Great St. Petersburg Polytechnic University. 29 Politechnicheskaya St., St. Petersburg, 195251, Russia. E-mail: [email protected]

ABITAEVA Rakhimash S. — Kazakh National Technical University (KazNTU). Satpaev str. 22 Almaty, Kazakh Republic. E-mail: [email protected]

BEKBAEV Amankeldy B. — Kazakh National Technical University named after K.I. Satpayev. Satpaev str. 22 Almaty, Kazakh Republic. E-mail: [email protected]

JAMAMBAEV Muratkali A. — Kazakh National Technical University named after K.I. Satpayev. Satpaev str. 22Almaty, Kazakh Republic. E-mail: [email protected]

Дата поступления статьи в редакцию: 22.11.2016.

© Санкт-Петербургский политехнический университет Петра Великого, 2017