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

Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. Т. 61, № 2 (2018), с. 97-107 Energetika. Proc. CIS Higher Educ. Inst. and Power Eng. Assoc. V. 61, No 2 (2018), pp. 97-107 97

DOI: 10.21122/1029-7448-2018-61-2-97-107 УДК 621.315.176

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

Ю. В. Бладыко1*

''Белорусский национальный технический университет

© Белорусский национальный технический университет, 2018 Belarasian National Technical University, 2018

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

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

Для цитирования: Бладыко, Ю. В. Механический расчет гибких токопроводов при замене сосредоточенных сил распределенной нагрузкой / Ю. В. Бладыко // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2018. Т. 61, № 2. С. 97-107. DOI: 10.21122/10297448-2018-61-2-97-107

Адрес для переписки

Бладыко Юрий Витальевич Белорусский национальный технический просп. Независимости, 65/2 220013, г. Минск, Республика Беларусь Тел.: +375 17 292-71-93 eie@bntu.by

Address for correspondence

Bladyko Yurii V.

Belarusian National Technical University 65/2 Nezavisimosty Ave., 220013, Minsk, Republic of Belarus Tel.: +375 17 292-71-93 eie@bntu.by

Mechanical Calculation of Flexible Wires

when the Concentrated Forces

are Replaced by a Distributed Load

Y. V. Bladyko1*

'-Belarasian National Technical University

Abstract. The objective of mechanical calculation of flexible wires of substations and overhead lines is to determine the sag and tension in different climatic conditions. A wire with a uniformly distributed load is considered as a homogeneous flexible thread having the form of a parabola. Concentrated loads from spacers, barrier balls, stubs, taps to electrical apparatus and other elements are replaced with the one distributed over the span. On behalf of a span without tension springs of insulators the action of concentrated loads on the wire is considered, an error is determined when replacing the concentrated forces with the one uniformly distributed along the span by the load. The sag for the equivalent wire is determined with the aid of the ratio of the increase of the sag, due to the presence of concentrated forces. An equation of state has been drawn up which makes it possible to determine the tension after changing the number of concentrated loads, e.g., after installing spacers, hanging the barrier balls, fixing the taps. The dependence of the maximum sag on the number of concentrated forces is given. The definition of the coefficient of concentrated forces as the ratio of the sum of the concentrated loads to the weight of the wire in the span is presented. A relationship between the load factors, the increase of the sag and the coefficient of concentrated forces is established. The formula has been deduced for determining the error in the replacement of concentrated forces by a uniformly distributed load along the span as a function of the number of concentrated forces and the coefficient of concentrated forces. A decrease in the error with an increase in the number of concentrated forces has been demonstrated. A more accurate calculation of mechanical tensions and sag is possible with the use of a vector-parametric method for calculating the flexible busbar of substations and air-line wires, where the design model of wires in the form of a flexible elastic thread is placed, taking into account the spatial disposition of all structural elements.

Keywords: sag, tension, flexible thread, chain line, parabola, equation of state, span, equivalent wire, load factor

For citation: Bladyko Y. V.

Задачей механического расчета гибких проводов распределительных устройств (РУ) и воздушных линий (ВЛ) является определение стрел провеса и тяжений в различных режимах климатических воздействий. Актуальная проблема современной электроэнергетики - проектирование РУ с уменьшенными габаритами и компактных ВЛ, что повышает требования к их расчету, обеспечению в любой точке пролета необходимых по ПУЭ минимальных расстояний до различных объектов при всех возможных видах нагрузок, принятых при проектировании (гололед, ветер и др.). В механическом расчете гибких проводов сосредоточенные нагрузки от распорок, заградительных шаров, гирлянд изоляторов, шлейфов, отпаек к электрическим аппаратам и других элементов заменяются распределенной по пролету [1]. Цель статьи - определение погрешности такой замены.

В приближенных методах провод с равномерно распределенной нагрузкой обычно рассматривается как однородная гибкая нить, имеющая очертание параболы. Расчет такой нити производится по уравнению состояния, вывод которого основан на сравнении длин нити при различных климатических условиях и разном нагружении [1].

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

В методике [1] уравнение состояния провода не отличается по структуре от обычного уравнения состояния, по которому рассчитываются провода ВЛ с равномерно распределенными нагрузками. Но для того чтобы уравнение было применимо для неоднородной нити, ее заменяют приведенной (эквивалентной) однородной нитью, что достигается путем введения приведенной (эквивалентной) равномерно распределенной по длине пролета нагрузки. Приведенная нагрузка определяется из условия равенства длины однородной нити с равномерно распределенной приведенной нагрузкой длине комбинированной нити с распределенными нагрузками разной интенсивности и в общем случае также длине нити с сосредоточенными нагрузками разной величины, приложенными в различных местах пролета.

Длина неоднородной (комбинированной) нити определяется путем интегрирования дифференциального уравнения нити по участкам. Для перехода от неоднородной нити к эквивалентной однородной при любой схеме загружения пролета определяется переходный коэффициент ^ называемый коэффициентом нагрузки. Приведенная нагрузка вычисляется как произведение величин равномерно распределенной нагрузки для данного провода и коэффициента нагрузки.

В большинстве случаев трудоемкого непосредственного интегрирования по участкам можно воспользоваться приемом перемножения эпюр «балочных» поперечных сил - приемом Верещагина [2]. Могут быть использованы и другие приближенные способы [3, 4].

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

Гибкая однородная нерастяжимая тяжелая нить с закрепленными концами в однородном гравитационном поле принимает форму цепной линии [5].

По уравнению цепной линии стрела провеса в середине пролета

/о=* [ ^ \)-1

где а = H/q; I - длина пролета; Н - тяжение в проводе (горизонтальная составляющая); q - погонный вес провода. Длина провода

Ь _ 2а  sh \ — |.

12а )

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

I2 /4 ql2 qъlA

/0 = —2-+ "7-Г + ... _ +-Г +...

22  2! а 24  4!а3 8Н 384Н3

при длине провода

_ . 13 15 . q2l3 Ь _ I +—;-+ ^-г + ... = I + - ' '

22  3! а1 24  5! а4 24Н2 1920Н4

Учитывая, что для встречающихся на практике случаев 2а >> I, пользуются только первым слагаемым, соответствующим представлению провода параболой [1]:

/ _ ^ /о = 8Н'

что дает погрешность порядка 0,2 %.

Длина провода определяется при этом как

ь = I +

24Н 2 '

Стрела провеса для эквивалентного провода

/ _ ql2К

/о _ "

8Н

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

Длина эквивалентного провода

-213 К 2

Ь _ I + -

где К - коэффициент нагрузки.

24Н 2

Расчетная модель пролета показана на рис. 1а.

а

RA _ _ _ ' 1 \ 1 \ 1 \ '< с l > /f\RB / 1 / 1 / 1

A B

HA A Pi у x - 1 1 УВ Hb q(x) \ 1 \<p„ YP

b

Рис. 1. Расчетная схема пролета: а - провод с вертикальными нагрузками, распределенными q(x) и сосредоточенными P; b - провод в виде простой разрезной балки с шарнирными опорами, загруженной так же, как и провод

Fig. 1. Estimated span scheme: a - wire with vertical distributed loads q(x) and concentrated loads P; b - wire in the form of a simple split beam with hinged supports, loaded in the same way as the wire

Вертикальные реакции А и В (рис. 1b) можно назвать балочными реакциями, как это принято, например, при определении вертикальных реакций трехшарнирной арки или вантовых систем [5].

Составим уравнение моментов относительно какого-нибудь произвольно выбранного сечения провода с координатами x и y от всех внешних сил, расположенных слева от этого сечения. Учитывая, что изгибающий момент в любом сечении провода (как гибкой нити) равен нулю, получим [1]

XM = Ax - XMq -XMp -Hy = 0,

где A - опорная балочная реакция; x, y - координаты провода; 7LMq - сумма моментов распределенных нагрузок; IMP - сумма моментов сосредоточенных сил.

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

как провод. Обозначая этот балочный изгибающий момент через M(x), получим

M

y = —(1) H

Так как M(x) есть переменная величина балочного изгибающего момента в зависимости от абсциссы x, то по формуле (1) может быть определена стрела провеса в любом месте пролета.

В середине пролета x = //2, и максимальная стрела провеса f0 = = y = M(//2)/H.

Приведенная (эквивалентная) нагрузка q = qK .

Напряжение в проводе с = H/F, где F - площадь сечения провода.

Определив для всех режимов эквивалентные погонные q и удельные у = q/F нагрузки, составляется уравнение состояния. При этом следует учитывать упругое и температурное удлинения провода на длине /1. Тогда получим разность длин провода в двух режимах, из которых один исходный (с индексом 0):

AL = Ш-Ш = ^(с-со) + аА (f-0),

24с2 24с0

откуда

L^L = с0-î^-aE

24с2 24с2 V '

где у = уК; у0 = у0К0; E - модуль упругости провода; а - коэффициент температурного удлинения провода; t - температура провода.

При неизменной температуре t = t0 = const уравнение состояния (2) можно записать в виде

н _ q2l2EFK2 = н _ q20l2EFKl

24Н2 0 24Н2

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

Коэффициент нагрузки при отсутствии натяжных гирлянд изоляторов и одинаковых высотах подвеса проводов на опорах [1]

K2 =1 + 12YsPab + 1

2 ;2

Q/2 Г Q2/

ZP2ab

RaY Pbl + Pa У Pbl + ... + Pn 2 а 2 У Pbl + Pn a Pb

1 1 ^^ il 2 2^^ ï ï n-2 n-2 ^^ ï ï n-1 n-1 n n

n-1

где Q = ql - вес провода в пролете без учета провиса; ai = li/(n + 1) - расстояние до опоры А i-й сосредоточенной силы при их равномерной расстановке вдоль пролета; bi = l - ai = l(n + 1 - i)/(n + 1) - расстояние до опоры В i-й сосредоточенной силы; Pi = P/n - вес i-й сосредоточенной нагрузки; P - суммарный вес всех сосредоточенных нагрузок; n - число сосредоточенных сил в пролете.

После подстановки получим

K\n) = 1 + 2 KP + Kp, n +1 n

где KP = P/Q - коэффициент сосредоточенных сил.

При большом числе сосредоточенных сил коэффициент нагрузки

К = lim К

После расчета опорных балочных реакций

A = B = (Q + P)/2 = Q(1 + Kp)/2 по (1) можно определить максимальную стрелу провеса

/о (n ) = H

ql2 Pl _ PI ff Г1__т_

8 + 4 n [ 2 n +1

После упрощения получим

. . ql2Kf

/о

где

Кг (п ) = 1 + КР [1 +1 j.

При п ^ да провод можно считать загруженным распределенной нагрузкой q + Р/1 = д(1 + КР), т. е. коэффициент нагрузки К = К(да) = 1 + КР, что подтверждает ранее полученное выражение (3).

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

При совместном действии вертикальных и горизонтальных нагрузок стрелы провеса и горизонтальные отклонения находят независимо в двух взаимно перпендикулярных плоскостях по формуле (1).

Расчеты выполнялись для следующих исходных данных: длина пролета ВЛ I = 400 м; провод 264/34 ACSR с погонным весом q = 0,98 даН/м,

площадью поперечного сечения ¥ = 297,8 мм , модулем упругости Е = = 7400 даН/мм2; тяжение Н = 1608 даН; суммарный вес сосредоточенных нагрузок в пролете Р = 35 даН. Рассматривался пролет без натяжных гирлянд изоляторов, с одинаковыми высотами подвеса, без ветра и гололеда, с неизменным тяжением. Для РУ длина пролета I = 40 м. Зависимость максимальной стрелы провеса пролета ВЛ от числа сосредоточенных сил приведена на рис. 2, то же для пролета РУ - на рис. 3.

14,5

fo(n), м

14,0

13,5

13,0

1

2 е—е—в—в—

0

5

10

15

Рис. 2. Зависимость стрелы провеса пролета воздушной линии от числа сосредоточенных сил: 1 - расчет по формуле (4); 2 - стрела провеса эквивалентного провода

Fig. 2. Dependence of the boom of the sag of overhead line span sag on the number of concentrated forces: 1 - calculation by the formula (4); 2 - boom of the sag of the equivalent cable

п

Рис. 3. Зависимость стрелы провеса пролета распределительного устройства от числа сосредоточенных сил: 1 - расчет по формуле (4); 2 - стрела провеса эквивалентного провода

Fig. 3. Dependence of the boom of the sag span of the switchgear on the number of concentrated forces: 1 - calculation by the formula (4); 2 - boom of the sag of the equivalent cable

Как видно из рис. 2, наибольшая абсолютная погрешность при п = 1 составляет около одного метра, что не позволительно для проектирования.

Для частного случая, когда в середине пролета действует одна сосредоточенная сила, стрела провеса

f - J_

f = H

ql + P l 2 2

ql2

ql2 Kf 8H

где

Kf - 1 + 2KP.

Коэффициент нагрузки при одной сосредоточенной силе в середине пролета определяется по формуле

2 P

K2 - 1 + 3— + 3

Г P V

ql V ql

-1 + 3KP (1 + KP).

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

5 f -

fo

ql 2

8 H

8H

ql2 Kf (n ) 8H

1 + n

v KP J

Как следует из рис. 2, погрешность от замены п = 1 сосредоточенной силы равномерно распределенной нагрузкой по длине пролета ВЛ составляет 7,5 %, при п = 3 - до 2,7 %, при п = 5 - до 1,6 %, при п = 7 - до 1,2 %. Для пролета РУ (рис. 3) при п = 1 погрешность возрастает до 32 %. Таким образом, самые большие погрешности замены сосредоточенных сил распределенной нагрузкой будут при их малом числе. Выражение (5) позволяет оценить погрешность при любом числе сосредоточенных сил п и коэффициенте КР.

При наличии натяжных гирлянд изоляторов, а также при учете ветровой нагрузки погрешность возрастает из-за искажения формы системы «гирлянды изоляторов - провод», отклонения составляющих проходят в разных плоскостях. В [6] показано увеличение погрешности с ростом длины гирлянды изоляторов. Наличие отпаек требует учета их реального расположения в РУ. Методика [1] может быть распространена на пролеты с гирляндами изоляторов, заградительными шарами, ветром, гололедом. Однако погрешность замены сосредоточенных сил равномерно распределенной нагрузкой возрастает.

Более точный расчет механических напряжений возможен при применении векторно-параметрического метода расчета гибкой ошиновки РУ и проводов ВЛ, где положена расчетная модель проводов в виде гибкой упругой нити [6-8].

1

Дифференциальные уравнения второго порядка, описывающие статику гибкой упругой нити [6], являются нелинейными. Их численное решение производится разностным методом с помощью ЭВМ [8-10]. Система конечно-разностных уравнений решается на основе вложенных друг в друга итераций: относительно координат, тяжения или длины провода. При наличии сосредоточенных нагрузок в виде ответвлений к электрическим аппаратам добавляется еще одна итерация по координатам точек их крепления к гибким шинам РУ. Эти координаты являются функцией тяже-ния в отпайке, а оно в свою очередь зависит от положения провода. Произвольное расположение отпаек дает составляющие сосредоточенных сил, действующие во всех направлениях. Как показано в [6], расхождение результатов по [1] можно снизить при более точном задании этих составляющих. Разработанная программа позволяет учесть до трех отпаек, до двух шлейфов и гирлянд изоляторов, в то время как методика [1] дает при малом числе сосредоточенных нагрузок наибольшую погрешность.

ВЫВОДЫ

1. Дано определение коэффициенту сосредоточенных сил ^ как отношению суммы сосредоточенных нагрузок к весу провода в пролете.

2. Установлена связь между коэффициентом нагрузки K и коэффициентом сосредоточенных сил

3. Установлена связь между коэффициентом увеличения стрелы провеса .K, обусловленным наличием сосредоточенных нагрузок, и коэффициентом сосредоточенных сил

4. Выведены формулы определения погрешности замены сосредоточенных сил равномерно распределенной нагрузкой по длине пролета как функции от числа сосредоточенных сил п и коэффициента сосредоточенных сил ^. Показано уменьшение погрешности при увеличении числа сосредоточенных сил п.

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

ЛИТЕРАТУРА

1. Бошнякович, А. Д. Расчет проводов подстанций и больших переходов ЛЭП / А. Д. Бош-

някович. Л.: Энергия, 1975. 248 с.

2. Качурин, В. К. Гибкие нити с малыми стрелками / В. К. Качурин. М.: Гостехиздат, 1956.

224 с.

3. Розанов, Г. Н. Некоторые вопросы расчета механической части воздушных линий / Г. Н. Розанов. М.; Л.: Госэнергоиздат, 1954. 224 с.

4. Тер-Мкртчан, Л. Н. Механический расчет проводов для открытых подстанций / Л. Н. Тер-

Мкртчан // Электричество. 1935. № 4. С. 35-45.

5. Глазунов, А. А. Основы механической части воздушных линий электропередачи. Т. 1: Работа и расчет проводов и тросов / А. А. Глазунов. М; Л.: Госэнергоиздат, 1956. 192 с.

6. Стрелюк, М. И. Численный метод расчета статики гибкой ошиновки ОРУ в различных режимах климатических воздействий / М. И. Стрелюк, И. И. Сергей, Ю. В. Бладыко // Энергетика. Изв. высш. учеб. заведений. 1983. № 8. С. 8-14.

7. Анализ действия гололедно-ветровых и электродинамических нагрузок в пролетах с произвольным расположением проводов / И. И. Сергей [и др.] // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2012. № 1. С. 38-44.

8. Стрелюк, М. И. Расчет на ЭВМ статики гибкой ошиновки РУ / М. И. Стрелюк, И. И. Сергей, Ю. В. Бладыко // Научные и прикладные проблемы энергетики: межвед. сб. Вып. 12. Минск: Вышэйш. шк., 1985. С. 75-79.

9. Стрелюк, М. И. Программа механического расчета гибкой ошиновки РУ в различных режимах климатических воздействий / М. И. Стрелюк, И. И. Сергей, Ю. В. Бладыко // Инв. № ГосФАП - П007594, инв. № РФАП БССР - 00143. Минск, 1984. 12 с.

10. Стрелюк, М. И. Численный метод механического расчета гибких токопроводов электроустановок энергосистем / М. И. Стрелюк, И. И. Сергей, Ю. В. Бладыко // Энергетика. Изв. высш. учеб. заведений. 1995. № 3-4. С. 21-29.

Поступила 28.11.2017 Подписана в печать 05.02.2018 Опубликована онлайн 30.03.2018

REFERENCES

1. Boshnyakovich A. D. (1975) Calculation of the Wires of Substations and of Large Transitions of Transmission Lines. Leningrad, Energia Publ. 248 (in Russian).

2. Kachurin V. K. (1956) Flexible Threads with Small Sags. Moscow, Gostekhizdat Publ. 224 (in Russian).

3. Rozanov G. N. (1954) Some Issues of Calculating the Mechanical Part ofAir Lines. Moscow; Leningrad, Gosenergoizdat Publ. 224 (in Russian).

4. Ter-Mkrtchan L. N. (1935) Mechanical Calculation of Wires for Open Substations. Elektrich-estvo [Electricity], (4), 35-45 (in Russian).

5. Glazunov A. A. (1956) Fundamentals of the Mechanical Part of Overhead Transmission Lines. Vol.1. Work and Calculation of Wires and Cables. Moscow; Leningrad, Gosenergo-izdat Publ. 192 (in Russian).

6. Strelyuk M. I., Sergey I. I., Bladyko Y. V. (1983) Numerical Method for Calculating the Statics of a Flexible Busbars in Various Modes of Climatic Influences. Energetika. Izvestiya Vys-shikh Uchebnykh Zavedenii [Energetika. Proceedings of the CIS Higher Educational Institutions], (8), 8-14 (in Russian).

7. Sergey I. I., Bladyko Y. V., Panamarenka Y. G., Tsemekhman B. D., Tarasov V. E. (2012) Analysis of the Effect of Ice-Wind and Electrodynamic Loads in Spans with Arbitrary Arrangement of Wires. Energetika. Izvestiya Vysshikh Uchebnykh Zavedenii i Energeti-cheskikh Ob'edenenii SNG = Energetika. Proceedings of the CIS Higher Educational Institutions and Power Engineering Associations, (1), 38-44 (in Russian).

8. Strelyuk M. I., Sergey I. I., Bladyko Yu. V. (1985) Calculation on the Computer of Statics of a Flexible Busbars of the Switchgear. Scientific and Applied Problems of Power Engineering. Minsk, Vysheishaya Shkola [Higher School] Publ., (12), 75-79 (in Russian).

9. Strelyuk M. I., Sergey I. I., Bladyko Yu. V. (1984) Program of Mechanical Calculation of the Flexible Busbars of the Switchgear in Various Types of Climatic Influences. Accession Number GosFAP - P007594, Accession Number of RFAP of the BSSR - 00143. Minsk. 12 (in Russian).

10. Strelyuk M. I., Sergei I. I., Bladyko Y. V. (1995) Numerical Method of Mechanical Calculation of Flexible Current Conductors of Power Plants of Power Systems. Energetika. Izvestiya Vysshikh Uchebnykh Zavedenii [Energetika. Proceedings of the CIS Higher Educational Institutions], (3-4), 21-29 (in Russian).

Recеived: 28 November 2017 Accepted: 5 February 2018 Published online: 30 March 2018