CC BY
62
НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 - 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
О выборе базиса для моделирования движения провода ЛЭП методом Галеркина # 09, сентябрь 2013 Б01:10.7463/0913.0602290 Иванова О. А.
УДК 533.6.011
Россия, МГТУ им. Н.Э. Баумана [email protected]
Введение
Известно, что провода воздушных линий электропередачи (ЛЭП), находящиеся под воздействием поперечного ветра, подвержены пляске (галопированию) — высокоамплитудным низкочастотным колебаниям [1]. Высокие динамические нагрузки, действующие на провода и арматуру ЛЭП при пляске, могут привести к их повреждению, поэтому в последние десятилетия активно разрабатываются алгоритмы и программы для численного моделирования аэроупругих колебаний проводов. Наиболее распространенным методом решения данной задачи является метод конечных элементов [2, 3], однако применяются также метод конечных разностей [4] и метод Галеркина [5].
Характер движения провода ЛЭП во многом определяется собственными формами и частотами его малых свободных колебаний. Поэтому при моделировании движения проводов решение часто ищется в виде разложения по собственным формам, при этом нахождение числа собственных форм, которые необходимо учитывать для получения качественно и количественно верного результата, является актуальной задачей. Известно [6], что в зависимости от тяжения первой собственной частоте колебаний провода в плоскости его начального провисания отвечает либо симметричная собственная форма, либо антисимметричная. По этой причине можно ожидать, что поведение провода в воздушном потоке в зависимости от тяжения будет иметь качественно различный характер. Однако в некоторых работах (например, [2]) в расчете используется лишь одна собственная форма, при этом неясно, каков будет результат расчета при использовании большего числа форм. В данной работе исследована зависимость параметров пляски провода ЛЭП в воздушном потоке от числа используемых базисных функций при решении задачи методом Галеркина.
1. Постановка задачи
В качестве модели провода принята модель абсолютно гибкого стержня (нити) [7], линейно упругого по отношению к растяжению. Его движение описывается уравнениями
д
( д&т) джг(е,т)\ + _ + .
дЛ 1 + д(е,т+ С + ^) - -
де
+
дХ (е,Т )
де
+
дХз(е,т)
де
1+
д 2хг(е,т) дт2
з(е,т)
0,
1, 2, 3,
а
(1)
Здесь введены следующие безразмерные параметры: е £ [-0,5; 0,5] — дуговая координата на нерастянутом проводе; т — время; ^(е, т) — тяжение; а — жесткость на растяжение; жг(е,т), г = 1, 2, 3, — декартовы координаты точки провода (рис. 1); С — функция, зависящая от скорости движения провода и задающая внутреннее демпфирование; 5^(е,т), г = 1, 2, 3, — распределенная аэродинамическая сила. Ускорение свободного падения д направлено противоположно оси Ох3 и его влияние в уравнениях (1) представлено слагаемым ^¿3 (8¿j — дельта Кронекера). При приведении к безразмерному виду величины, имеющие размерности длины, массы и силы, отнесены соответственно к длине, массе и весу провода; время отнесено к у/Ь/д, где Ь — длина нерастянутого провода.
Рис. 1. Расчетная схема
В качестве начального условия для системы (1) принимается равновесное положение провода хг0(е), г =1, 2, 3, и равновесное тяжение ф0(е) [8]. В дальнейшем будем искать решение в виде
2*(е,т) = хго(е)+«де,т), г = 1,2,3, д(е,т) = а(е) + д$(е,т),
где «¿(е, т) и дд(е, т) —отклонения от равновесной формы провода и равновесного тяжения соответственно. Концы провода считаются закрепленными неподвижно:
«¿(±0,5, т) = 0, г = 1, 2, 3. 2. Решение уравнений движения провода методом Галеркина
В соответствии с методом Галеркина отклонения провода от положения равновесия представляются в виде линейной комбинации некоторого набора базисных функций, в качестве которых целесообразно выбрать собственные формы малых свободных колебаний провода:
(е,т) = £ «(¿)(т(е), г = 1, 2, 3.
(2)
к= 1
2
2
2
Здесь Si, S2, S3 — число базисных функций. Такой подход позволяет, во-первых, корректно моделировать динамику провода, а во-вторых, задать внутреннее демпфирование в соответствии с общепринятым подходом пропорциональным частоте колебаний в виде
ак,т ) = 2в £ ^ Л),
k= 1 Т
где в — коэффициент демпфирования; — k-я собственная частота колебаний в направлении оси ÜXi.
Изменение тяжения AQ(£, т) также представляется в виде линейной комбинации некоторого набора базисных функций, например, тригонометрических:
SQ
AQ(C,T ) = £ akQ)(T )^Q)(£), k=1
^1Q) = 1, viQ) = л/2 sin 2n = /2 cos 2п £ и т. д.
Здесь Sq — число базисных функций. Отметим, что при Sq = 1 изменение тяжения будет постоянным вдоль пролета, как это предполагается, например, в работах [4, 5]: AQ(£, т) = = aiQ)(T). Заметим также, что в работе [4] равновесное тяжение Qo тоже считалось постоянным вдоль пролета.
3. Определение аэродинамической нагрузки
Для определения величины распределенной аэродинамической силы, действующей на провод, применим общепринятый подход, в соответствии с которым обтекание всех сечений провода считается плоскопараллельным, а удельная нагрузка в каждом конкретном сечении равна
(0, ql <¿)T = 1 РвоздЬК2тн(тVCxa(§) + VVСуа(§)),
где рвозд — плотность воздуха; b — хорда (диаметр) провода; Ктн — относительная скорость ветра для данного сечения; тV, vV — орты, направленные вдоль Vотн и перпендикулярно ему; Сха(§), Суа(§) — стационарные аэродинамические коэффициенты профиля сечения, известные из эксперимента или расчета, зависящие от угла атаки § — угла между вектором Vотн и хордой профиля.
Далее в расчетах приняты следующие зависимости аэродинамических коэффициентов профиля сечения провода от угла атаки (рис. 2): Сха(§) = 2,0 + 0,2 cos§, Суа(§) = —1,5 sin 2§. Данные зависимости в окрестности § = 0 количественно и качественно соответствуют аэродинамическим коэффициентам профиля с U-образным обледенением [3]. Кроме того, для выбранных зависимостей Сха(§) и Суа(§) в области § е [—0,37; 0,37] выполнено условие неустойчивости Глауэрта — Ден-Гартога Сха(§) + Суа(§) < 0.
4. Собственные частоты и формы малых колебаний провода
На рис. 3 приведены зависимости квадратов безразмерных собственных частот малых колебаний провода в плоскости начального провисания от стрелы провеса w = |x30(0) обратно пропорциональной тяжению, при значении а = 6906. Из графика видно, что при значении стрелы провеса w ~ 0,022 первая собственная форма меняется с симметричной (четной функции на антисимметричную (нечетная функция £). На рис. 4 для примера приведены две первые собственные формы колебаний в направлении оси Ox3 при значениях стрелы провеса w = 0,018 и w = 0,030. Собственные формы найдены методом сагиттарной функции [9].
0.01 0.02 0.03 0.04
Рис. 3. Зависимость квадратов собственных частот от стрелы провеса
5. Результаты расчетов и их анализ
В расчетах с различным числом базисных функций для отклонений провода и его тяжения принималось, что жесткость провода на растяжение а = 6906; стрела провеса принимала значения w = 0,018 и w = 0,030; число базисных функций для всех перемещений было одинаковым: Si = S2 = S3 = S.
В случае сравнительно малой стрелы провеса w = 0,018, когда тяжение велико и первая форма колебаний является симметричной, во всех расчетах была получена периодическая траектория движения провода; колебания представляли собой стоячую волну. В табл. 1 приведены значения удвоенной амплитуды (размаха) колебаний провода. Частота колебаний во всех расчетах была приближенно равна 13,2, за исключением расчета при S = Sq = 1, в котором она была близка к 13,1.
Таблица 1. Размах колебаний провода при малой стреле провеса
S = 1 S = 2 S = 4 S = 6 S = 8
Sq = 1 0,0086 0,0086 0,0084 0,0083 0,0083
Sq = 3 - 0,0086 0,0084 0,0084 0,0083
Приведенные в таблице данные показывают, что увеличение числа базисных функций приводит к незначительному, на 2-3 процента, уменьшению размаха колебаний. Частота колебаний практически не меняется. Проведенные расчеты демонстрируют, что при моделировании колебаний сильно натянутого провода достаточно использовать по одной базисной функции в каждом направлении.
В случае стрелы провеса т = 0,030 движение провода имеет качественно другой характер (табл. 2). Проведя расчет с одной базисной функцией в каждом направлении, можно, как и в предыдущем примере, получить движение в форме стоячей волны. Увеличение числа базисных функций до двух приводит к решению в форме бегущей волны. Дальнейшее увеличение числа базисных функций вид решения качественно не меняет, добавление третьей и четвертой форм колебаний приводит к уточнению решения на 38%, пятой и шестой — еще на 7 %. Использование большего числа базисных функций представляется нецелесообразным, поскольку оно практически не оказывает влияния на решение. Количество SQ базисных функций для тяжения провода должно быть согласовано с числом S базисных функций для перемещений.
Таблица 2. Размах колебаний провода при увеличенной стреле провеса
S = 2 S = 4 S = 6 S = 8
Sq = 1 0,0115 0,0086 0,0083 0,0054
SQ = 3 0,0112 0,0086 0,0086 0,0073
SQ = 5 0,0110 0,0081 0,0076 0,0076
Заключение
Проведенные расчеты показывают, что для корректного моделирования движения провода ЛЭП методом Галеркина одной базисной функции, как это делается во многих известных работах, в общем случае недостаточно. Это особенно актуально для пролетов ЛЭП с большими стрелами провеса и сравнительно малым тяжением. В этом случае колебания представляют собой суперпозицию бегущих волн и установившееся движение отдельных сечений провода не является периодическим. В то же время для сильно натянутых проводов колебания имеют вид стоячей волны, установившееся движение сечений провода — периодическое, и для описания колебаний провода достаточно использовать по одной базисной функции (собственной форме малых колебаний) в каждом направлении.
Работа выполнена при финансовой поддержке гранта Президента РФ для государственной поддержки молодых российских ученых — кандидатов наук (проект МК-6482.2012.8) и гранта Президента РФ для государственной поддержки ведущих научных школ (проект НШ-255.2012.8).
Список литературы
1. EPRI. Transmission Line Reference Book, Wind Induced Conductor Motion. Palo Alto (California): Electrical Power Research Institute, 1979. 255 p.
2. Desai Y.M., Yu P., Popplewell N., Shah A.H. Finite Element Modelling Of Transmission Line Galloping // Comp. and Struct. 1995. Vol. 57. P. 407-420.
3. Wang X., Lou W.-J. Numerical Approach to Galloping of Conductor // Proc. of the 7th Asia-Pacific Conference on Wind Engineering. Taipei, Taiwan, 2009. 8 p.
4. Сергей И.И., Виноградов A.A. Численное моделирование эксплуатационных статических и динамических режимов проводов ВЛ и кабелей // Электрические станции. 1998. № 1. С. 41-49.
5. Luongo A., Zulli D., Piccardo G. Analytical and Numerical Approaches to Nonlinear Galloping of Internally Resonant Suspended Cables // Journal of Sound and Vibration. 2008. Vol. 315, no. 3. P. 375-393.
6. Irvine H.M., Caughey T.K. The Linear Theory of Free Vibrations of a Suspended Cable // Proceedings of the Royal Society of London. Ser. A. 1974. Vol. 341. P. 299-315.
7. Светлицкий B.A. Механика абсолютно гибких стержней /Под ред. А.Ю. Ишлинского. М.: Изд-воМАИ, 2001. 432 с.
8. Иванова О.А. Приближенные методы определения собственных частот колебаний проводов многопролетных линий электропередач // Вестник МГТУ им. Н.Э. Баумана. Естественные науки. 2011. Спец. выпуск ««Прикладная математика»». С. 34-44.
9. Akulenko L., Nesterov S. High-precision methods in eigenvalue problems and their applications. Boca Raton, CRC Press, 2004. 235 p.
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Choice of basis for simulation of ETL wire by the Galerkin method
# 09, September 2013 DOI: 10.7463/0913.0602290 Ivanova O. A.
Bauman Moscow State Technical University 105005, Moscow, Russian Federation [email protected]
A problem of building a system of base functions suitable for solving motion equations of a wire by the Galerkin method was considered in this paper. Dependence of the peak-to-peak amplitude of oscillations on the number of base functions was studied by a test example. Two sets of wire parameters were considered; wire motion is qualitatively different for these sets. Use of one base function in each direction is reasonable for a taut wire and it leads to an error less than 4%. For a slack wire usage of not less than six base functions is required for correct simulation of wire motion.
References
1. EPRI Transmission Line Reference Book: Wind Induced Conductor Motion. Palo Alto, California, Electrical Power Research Institute (EPRI), 1979. 255 p.
2. Desai Y.M., Yu P., Popplewell N., Shah A.H. Finite Element Modelling of Transmission Line Galloping. Computers and Structures, 1995, vol.57, no. 3, pp. 407-420. DOI: 10.1016/0045-7949(94)00630-L.
3. Wang X., Lou W.-J. Numerical Approach to Galloping of Conductor. Proc. of the 7th Asia-Pacific Conference on Wind Engineering. Taipei, Taiwan, 2009. 8 p.
4. Sergey I.I., Vinogradov A.A. Chislennoe modelirovanie ekspluatatsionnykh staticheskikh i dinamicheskikh rezhimov provodov VL i kabeley [Numerical simulation of static and dynamic operation condition of wires of overhead lines and cables]. Elektricheskie stantsii, 1998, no. 1, pp. 41-49.
5. Luongo A., Zulli D., Piccardo G. Analytical and Numerical Approaches to Nonlinear Galloping of Internally Resonant Suspended Cables. Journal of Sound and Vibration, 2008, vol. 315, no. 3, pp. 375-393.
6. Irvine H.M., Caughey T.K. The Linear Theory of Free Vibrations of a Suspended Cable. Proceedings of the Royal Society of London. Ser. A, 1974, vol. 341, pp. 299-315.
7. Svetlitskiy V.A. Mekhanika absolyutno gibkikh sterzhney [Mechanics of absolutely flexible rods]. Moscow, MAI Publ., 2001. 432 p.
8. Ivanova O.A. Priblizhennye metody opredeleniya sobstvennykh chastot kolebaniy provodov mnogoproletnykh liniy elektroperedach [Approximate methods to determine eigenfrequencies of wire oscillations in multi-span transmission lines]. Vestnik MGTU im. N.E. Baumana. Ser. Estestvennye nauki [Herald of the Bauman MSTU. Ser. Natural science], 2011, spets. vyp. "Prikladnaya matematika" [spec. iss. "Applied mathematics"], pp. 34-44.
9. AkulenkoL., Nesterov S. High-precision methods in eigenvalue problems and their applications. Boca Raton, CRC Press, 2004. 235 p.
