CC BY
500
Динамика проводов ЛЭП
229
УДК 62-742
DOI: 10.17586/0021-3454-2015-58-3-229-235
ДИНАМИКА ПРОВОДОВ ЛЭП
ПРИ ДВИЖЕНИИ ПО НИМ ДИАГНОСТИЧЕСКОЙ МАШИНЫ
В. В. Елисеев, М. Р. Бахрами
Санкт-Петербургский государственный политехнический университет, 195251, Санкт-Петербург, Россия
E-mail: yeliseyev@inbox.ru
Рассматривается движение диагностической машины-автомата по проводам линии электропередачи. При математическом моделировании движения используется лагранжева механика твердых и упругих тел с применением аппарата компьютерной математики: выявлены пилообразные колебания в вертикальной плоскости и поперечные движения типа параметрического резонанса.
Ключевые слова: диагностика проводов ЛЭП, колебания, лагранжева механика, параметрический резонанс, компьютерная математика.
Для исследования динамики проводов линий электропередачи создаются специальные диагностические машины-автоматы (ДМА) [1, 2]. Так, М. Р. Бахрами разработана оригинальная ДМА [3], включающая в себя моторизованные ведущие колеса, „руку“ с захватом и центральную коробку (рис. 1, а). Машина может проходить через различные препятствия на проводах (амортизаторы, зажимы, предупреждающие шары).
Однако работа подобных устройств может быть нарушена вследствие реакции проводов на различные воздействия [4]; например, даже при медленном движении ДМА по проводам возможны опасные колебания, приводящие к ее отказам или разрушению. Это обусловливает необходимость математического моделирования движения ДМА на стадии проектирования, чтобы выбрать приемлемые варианты конструкции машины, снизить затраты на эксперимент и обеспечить безопасность работы персонала.
В настоящей статье представлены результаты по разработке методики моделирования движения ДМА, включающей вывод уравнений динамики машины и решение их средствами компьютерной математики (MathCad). Уравнения выведены на основе механики Лагранжа для твердых и деформируемых тел с применением вариационного исчисления [5].
Рис. 1
Допустим, что провод является натянутой струной, ДМА — точечной массой, и движение происходит в вертикальной плоскости x, y (рис. 1, б).
В известной постановке задачи выражение для прогиба струны u (x, t) имеет следующий
вид:
Tu" + P5(x - ^(t)) = pu; x = 0, l: u = 0; t = 0: u = 0, it = 0, (1)
где T — сила натяжения струны; P = mg — вес ДМА; p — погонная масса струны; l — длина струны; штрих и точка означают дифференцирование по координате x и времени t; ^(t) — заданный закон движения нагрузки, причем ^(0) = 0.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2015. Т. 58, № 3
230
В. В. Елисеев, М. Р. Бахрами
Выражение (1) представляет собой упрощенную постановку задачи, поскольку в динамике нагрузка может быть направлена не вертикально, и соответственно струна колеблется не в вертикальной плоскости.
Статическому приближению ( р = 0 ) соответствует условие
u = US <
, х < &
%
l - х
, х > &
Us (%) = -
P (1
,-1
& l - &
Динамическое решение задачи (1) можно найти методом собственных функций. При % = vt (движение ДМА с постоянной скоростью)
12 nn
u(хt) = Z un(t)фп (х); фп =Ji •sin (Кх), К =-r;
n=1 ' 1 l
P
К(v - c2) v Tpl
2 T
[vsin(knct)-csin(knvt)], c = —.
(2)
Как показали расчеты (в среде MathCad), прогиб в точке нагружения u (vt, t) отличается от квазистатического Us (vt) только пилообразными колебаниями струны.
Учтем массу струны. При законе движения %(t) и прогибе струны u (х, t) имеем следующие выражения для координат, скоростей и ускорений движущейся массы: х = %(t), y = u(%(t), t) ^ х = %, y = u"% + u, х = %, y = u"%2 + u"% + 2u"% + u.
Используем вариационную методику Лагранжа — Ритца — Канторовича [5]. Обобщен-
k
ные координаты U., i = 1...k, вводятся аппроксимацией прогиба u(х, t) = щ. (t )ф( (х) =
i=1
= U (t)ф(х) (в матричных обозначениях). Координатные функции ф. заданы такими же, как
в уравнениях (2). Удвоенная кинетическая энергия струны при р^0 определяется как
l
2K = m( х2 + у2) + р J u = т%2 + UT [M + p(%)]U + 2%UT n(%)U + % 2UT o(%)U, (3)
0
l
где введены матрицы M = J рфф'' Сх, p = тфф , о = тф"ф" , n = тфф" .
Определив потенциальную энергию
1 1 1
П = - T Ju"2dc - Pu(%, t) = -UTCU - UTF(%), C = T J ф'ф^Сс, F = Рф(%)
О J О J
(4)
и приняв % = v = const, составим уравнения Лагранжа: dK Л dK = дП ;
jU J dU ~ dU ’
dK ■ T
^ = (M + p)U + Щ U; dU
dK T • 2
— = vnTU + v2oU, dU
дП
dU
= CU - F ^ (M + p)U + [C + v2 (n"(%) - o)]U + v[p"(%) + nT - n]U = F.
(5)
х
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2015. Т. 58, № 3
Динамика проводов ЛЭП
231
Здесь производные по столбцу — это столбцы из производных по компонентам. Для квадратичных функций имеем d(UTCU)/dU = 2CU . При этом заметим, что п"©Т - ° = к, к(^) = тфф"Т, ц"(£) + nT - n = 2n .
Правую часть последнего уравнения в системе (5) преобразуем к каноническому виду для последующего численного интегрирования ( Л = (M + ц)-1):
U = V,
V = Л[ F-(с + v2 к )u - 2vnV ].
Важным для последующих расчетов является ускорение точечной массы a = v2u" + 2vu" + ii = v2ty"TU + 2vф,TU + фТ U =
= [v2ф"T - фТЛ С + v2к)] U + 2v ( - фТЛп)v + фТ AF.
Система (6) решается численно (программа Rkadapt в пакете MathCad) при нулевых начальных условиях.
Заданы следующие значения параметров: Т = 10 кН, P = 1 кН, v = 1м/с, р = 5 кг/м, l = 200 м, число степеней свободы к = 15 . На рис. 2, а—в соответственно представлены результаты расчета — временные зависимости прогиба струны в нагруженной точке Um (t) = u(vt, t), скорости Um и ускорения Um . Величина прогиба в данном случае сопоставима
с решением (2), но скачкообразные изменения скорости, характерные для решения (2), отсутствуют, что является следствием аппроксимации при вариационном подходе.
а) ит, м 5 4 3 2 1
0 100 200 t, c
Рис. 2
Резким изменениям скорости соответствуют ударные импульсы сил инерции (переносного движения), возбуждающие параметрические колебания ДМА. Рассмотрим более сложную модель с движущимся по струне физическим маятником (рис. 3).
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2015. Т. 58, № 3
232 В. В. Елисеев, М. Р. Бахрами
Координаты и скорость центра масс системы, при прогибах струны uy (x, t), uz (x, t) и угле отклонения 9(t) маятника, определяются выражениями
x = 4(t); y = uy (4, t) + s cos 9; z = uz (4, t) + s sin 9 ^ x = 4;
y = uy 4 + Uy - s9 sin 9; Z = uz4 + U z + s9 cos 9,
где s — расстояние от центра масс до точки подвеса.
Кинетическая энергия тела
Kт = 2(mvC + Ic92); vС = x2 + y2 + z2 = 42(1 + u’2 + u'2) + ii2y + u2z + s292 +
+2^(uyuy + u'zuz) + 2^s9(-uy sin 9 + u'z cos 9) + 2s9(-uy sin 9 + uz cos 9), где Ic — момент инерции относительно центра масс.
(7)
Д
А
Рис. 3
Уравнения динамики, базирующиеся на выражении вида (7), слишком сложны. Поэтому положим величину uy заданной и линеаризуем задачу, считая uz, 9 малыми величинами одного порядка. Тогда уравнение (7) преобразуется к виду
2 Кт = m
42u'2 + ^ + 2by'zuz + 2^s9(-u!, 9 + u'z ) + 2s9(-uy 9 + uz ) + (ms2 + Ic )92 +...
Используем, аналогично вышеизложенному, вариационный подход:
n l
uz(у t)=Z U(t)9i(x) = uT (tMx);2K=pj (uy2 + uz2 )dx+mv2+Ic9'1 =
i=1 0
= UT (M + ^(4))U + 4 2UT o(4)U + 24UT n(4 )U +
+2ms9(4UT ф '(4) + UT ф(4) - vm9) + (ms2 + Ic )92 +...,
где vm = uy + vu’y, а матрицы инерционных характеристик аналогичны принятым в выражении (3).
Потенциальная энергия
1 1 1
П =— Т j (u y2 + u’2)dx - mg(uy (4, t) + s cos 9) = — (UTCU + mgs92) +...
2 •* y y 2
0
Составим уравнения Лагранжа ( 4 = v = const):
ck)dK = an (ck ] ck an
dU) CU ~ CU, lc9 ) C9 = C9 .
Определив производные
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2015. Т. 58, № 3
Динамика проводов ЛЭП
233
5K T 5K 2 ■
—- = (M + ц)U + vnTU + meQф, ---= v2cU + v(gU + meQф');
dU dU
^ = (me2 + Ic )Q + шефT ф(^} + vuT ф'(£) - vmQ);
SQ
8K д 6И 6И e
---= -mevjd; = CU, = mgeQ,
ae m su se
перепишем уравнения следующим образом:
(M + p,)U + 2vn U + (v2 к + C }U + meQ ф = 0, (me2 + Ic)Q + me\qTU + 2vф,TU + v2ф''TU + (g -am}Q
= 0,
am = vm = u y + 2vu'y + v2u’y ■
Далее разрешим систему относительно вторых производных; обозначив U = V, Q = ю, запишем систему обыкновенных дифференциальных уравнений в каноническом виде:
(U > ( 0 I 0 0 Л
Y = F(t, Y) = A(t)Y; Y - V e , A = A21 0 A22 0 A23 0 0 1
lю у 1A41 A42 A43 0 ,
A21 = М0 C0, A22 = 2avM0 п A23 =eeM0 g0ф,
A41 = Р (ФТМ0-1C0 - v2ф''T ), A42 = -2vP ( - афТМ-1П), ) = -Pg0 l1 + mPeфTмo'1ф),
М0 - М + ац, а -
2 2 e2 + r/
в -
2 2 e2 + С
C0 - C + v2
ак,
r
где rc — радиус-вектор центра масс.
Система решается в среде MathCad. На рис. 4, а, б представлены результаты расчета угла Q(t) и поперечного смещения uz (vt, t) системы при у- e = rc = 0,02 м и принятых выше
значениях остальных параметров. а) 0, о.е.-!
0,8 ■
0,4
0
-0,4
-0,8
50 100 150 t, c
б) Uz, о.е. 0,01
0,005
0
-0,005
-0,01
-0,015
0
50 100 150 t, c
0
в) 0, о.е.-0,8 -
0,4
0
-0,4
-0,8
mimm miiir мим
М1Г!! ' 1 TMIII'I J11‘Т'fi||
|| тн 'I ''М 1 •'
ii ii у mi iiLiiiuj iiijnui
0
50
100
150 t, c
Рис. 4
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2015. Т. 58, № 3
234
В. В. Елисеев, М. Р. Бахрами
Характер приведенных зависимостей позволяет говорить о параметрическом резонансе. Это — не классический тип резонанса, поэтому его можно исследовать, только используя аппарат компьютерной математики. В отличие от классического варианта [6], в данном случае отсутствует периодическое возбуждение. Однако при других значениях параметров системы результаты могут качественно измениться. Например, при у = 0,2 наблюдаются ограниченные параметрические колебания, причем значение угла 9 не превышает начального (рис. 4, в, г). Изменяя значение параметра у, найдем, что при у « 0,06 появляется резонанс: т.е. это значение у — есть верхняя граница зоны неустойчивости. Несмотря на то, что нижнюю границу определить не удалось, можно утверждать — при уменьшении значения у резонанс становится более выраженным.
Итак, при математическом моделировании движения ДМА по проводам линии электропередачи выявлены колебания системы — вертикальные пилообразные и поперечные параметрические. Многовариантные расчеты показывают, что параметрического резонанса можно избежать. Представленная методика позволяет построить области устойчивости в пространстве параметров.
СПИСОК ЛИТЕРАТУРЫ
1. Разумовский И. Т. Пирометр для контроля перегрева скруток проводов линий электропередач // Изв. вузов. Приборостроение. 2011. Т. 54, № 11. С. 43—48.
2. Toussaint K., Pouliot N., Montambault S. Transmission line maintenance robots capable of crossing obstacles: state-of-the-art review and challenges ahead // J. of Field Robotics. 2009. N 26(5). P. 477—499.
3. Bahrami M. R. A novel design of electrical transmission line inspection machine // Современное машиностроение. Наука и образование. 2014. № 4. С. 292—298.
4. Гимадиев Р. Ш., Динмухаметов Ф. Ф., Галимуллин Н. Р. Расчет пространственной динамики ЛЭП при совместном воздействии ветровой и весовой нагрузки // Изв. вузов. Проблемы энергетики. 2010. № 3—4.
С. 28—37.
5. Елисеев В. В. Механика деформируемого твердого тела. СПб: Изд-во Политехи. ун-та, 2006. 231 с.
6. Зевин А. А. К теории параметрических колебаний // Прикладная математика и механика. 2014. Т. 78, № 1.
С. 46—59.
Сведения об авторах
Владимир Васильевич Елисеев — д-р физ.-мат. наук, профессор; СПбГПУ, кафедра машиноведения и
основ конструирования; E-mail: yeliseyev@inbox.ru
Мохаммад Реза Бахрами — аспирант; СПбГПУ, кафедра машиноведения и основ конструирования;
E-mail: mr.bahrami1982@gmail.com
Рекомендована кафедрой машиноведения и основ конструирования
Поступила в редакцию 22.12.14 г.
Ссылка для цитирования: Елисеев В. В., Бахрами М. Р. Динамика проводов ЛЭП при движении по ним диагностической машины // Изв. вузов. Приборостроение. 2015. Т. 58, № 3. С. 229—235.
DYNAMIC OF ELECTRICAL TRANSMISSION WIRES AT A DIAGNOSTIC MACHINE MOVEMENT
ALONG THE LINE
V. V. Eliseev, M. R. Bahrami
Saint Petersburg State Polytechnical University, 195251, Saint Petersburg, Russia
E-mail: yeliseyev@nbox.ru
Traveling of an automatic diagnostic machine along electrical transmission line is considered. Mathematical modeling based on Lagrangian mechanics of rigid and flexible bodies is carried out with
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2015. Т. 58, № 3
Динамика проводов ЛЭП
235
the use of computer mathematics software. Saw tooth oscillations of the wires in the vertical plane and transversal movements of the type of a parametric resonance are revealed.
Keywords: diagnostic of electrical transmission line, oscillation, Lagrangian mechanics, parametric resonance, mathematical software.
Data on authors
Vladimir V. Eliseev — Dr. Sci., Professor; Saint Petersburg State Polytechnical University, De-
partment of Machine Science and Mechanical Design;
E-mail: yeliseyev@inbox.ru
Mohammad R. Bahrami — Post-Graduate Student; Saint Petersburg State Polytechnical University,
Department of Machine Science and Mechanical Design;
E-mail: mr.bahrami1982@gmail.com
Reference for citation: Eliseev V. V., Bahrami M. R. Dynamic of electrical transmission wires at a diagnostic machine movement along the line // Izvestiya Vysshikh Uchebnykh Zavedeniy. Priborostroenie. 2015. Vol. 58, N 3. P. 229—235 (in Russian).
DOI: 10.17586/0021-3454-2015-58-3-229-235
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2015. Т. 58, № 3
