Аннотация: В статье рассмотрено применение программной среды Mathcad для реализации оценочного расчёта электрических параметров ЛЭП методом зеркальных отображений.
Ключевые слова: метод зеркальных отображений, расчёт параметров ЛЭП, программная среда Mathcad.
Фізико-математичні науки
УДК 004.94
Мірошник Костянтин Артурович
магістрант кафедри передачі електричної енергії
Національний технічний університет
«Харківський політехнічний інститут»
Мирошник Константин Артурович
магистрант кафедры передачи электрической энергии
Национальный технический университет
«Харьковский политехнический институт»
Miroshnyk Kostiantyn
master of the department of electric power transmission
The National Technical University
"Kharkiv Polytechnic Institute"
РОЗРАХУНОК ПАРАМЕТРІВ ЛЕП МЕТОДОМ ДЗЕРКАЛЬНИХ ВІДОБРАЖЕНЬ У ПРОГРАМНОМУ СЕРЕДОВИЩІ MATHCAD
РАСЧЁТ ПАРАМЕТРОВ ЛЭП МЕТОДОМ ЗЕРКАЛЬНЫХ ОТОБРАЖЕНИЙ В ПРОГРАММНОЙ СРЕДЕ MATHCAD
COMPUTATION OF ELECTRICAL PARAMETERS OF ELECTRICITY TRANSMISSION LINES BY THE METHOD OF MIRROR CHARGES IN THE MATHCAD SOFTWARE ENVIRONMENT
Анотація: У статті розглянуто використання програмного середовища Mathcad для реалізації оціночного розрахунку електричних параметрів ЛЕП методом дзеркальних відображень.
Ключові слова: метод дзеркальних відображень, розрахунок параметрів ЛЕП, програмне середовище Mathcad.
Аннотация: В статье рассмотрено применение программной среды Mathcad для реализации оценочного расчёта электрических параметров ЛЭП методом зеркальных отображений.
Ключевые слова: метод зеркальных отображений, расчёт параметров ЛЭП, программная среда Mathcad.
Summary: The article describes the using of Mathcad software environment for the computation of the electric parameters of electricity transmission lines by the method of mirror charges.
Key words: the method of mirror charges, computation of the parameters of electricity transmission line, Mathcad software environment.
Теоретичні відомості
Метод дзеркальних відображень
Для розрахунку електростатичних полів, особливо обмежених якоюсь провідною поверхнею правильної форми або в яких є геометрично правильної форми межа між двома діелектриками, широко застосовується розрахунковий прийом, який називають методом дзеркальних відображень. Розрахунок поля заряджених провідників, розташованих поблизу плоских поверхонь, які обмежують провідне середовище, зводиться за допомогою методу дзеркальних відображень до розрахунку поля декількох провідників при відсутності провідного середовища.
Розглянемо поле прямолінійного дроту, розташованого на відстані h від плоскої поверхні провідного середовища. Це відповідає, наприклад, дроту, підвішеному на висоті над поверхнею землі. Усі лінії напруженості поля, які починаються на позитивно зарядженому дроті, закінчуються на поверхні провідного середовища, де з'являється індуктований негативний заряд. Поле визначається як зарядом дроту, так і усім зарядом, розподіленим по поверхні провідного середовища. Розподіл індуктованого заряду з умов задачі невідомо і так само підлягає визначенню.
На перший погляд завдання розрахунку поля в такій системі здається досить складним. Проте воно вирішується дуже просто за допомогою методу дзеркальних відображень. Усунемо провідне середовище і замінимо його дротом, що є дзеркальним відображенням реального дроту у поверхні розділу, та яке має заряд тієї ж величини, що і заряд реального дроту, але протилежного знаку.
Дійсний дріт і його дзеркальне відображення становлять двохпровідну лінію. Площина, розташована посередині між дійсним дротом і його дзеркальним відображенням, є поверхнею рівного потенціалу. У дійсних умовах поверхня провідного середовища якраз збігається з цією площиною і також є поверхнею рівного потенціалу. Звідси випливає, що якщо замінити провідне середовище дзеркальним відображенням дроту зі зміною знаку заряду, то в області над провідним середовищем поле залишиться таким самим, як і в дійсних умовах. В цьому і полягає метод дзеркальних відображень.
Цей метод можна застосовувати і при будь-якій кількості дротів, протягнутих паралельно один одному і паралельно плоскій поверхні, що обмежує провідне середовище. Кожен провід повинен бути дзеркально відображений у поверхні провідного середовища зі зміною знаку заряду, після чого провідне середовище може бути подумки вилучено, та розглянуто поле сукупності дійсних дротів і їх дзеркальних відображень. В такому полі площина, розташована на місці поверхні провідного середовища, є поверхнею рівного потенціалу, так як заряди протилежних знаків розміщені симетрично відносно цієї площини (рис. 1). Знайдене таким шляхом поле і буде дійсним полем в області над поверхнею провідного середовища.
Рис. 1. Схема заміщення трифазної ЛЕП з двома складовими у фазі
В якості системи заряджених тіл розглянемо багатопровідну лінію з n досить довгих дротів з зарядом τk на одиницю довжини (індекс заряду відповідає номеру проводу), протягнутих паралельно поверхні землі. Висота підвісу і радіус кожного дроту відомі. Відома також електрична проникність е середовища, поблизу дротів.
Візьмемо у діелектрику деяку довільну точку М (рис.2) і знайдемо її потенціал. Потенціал точки М дорівнює сумі потенціалів, створених кожним дротом і його дзеркальним зображенням.
Складова потенціалу точки М від дроту 1 і його дзеркального зображення запишеться наступним чином (постійну, з точністю, до якої визначається потенціал, опускаємо):
де - є відстань точки М до дзеркального зображення першого дроту; - є відстань точки М до першого дроту.
Рис 2. Модель розрахунку за методом дзеркальних відображень
Будемо вважати, що висоти підвісу дротів над землею на багато більше радіусів цих дротів. При цьому електричні осі практично співпадають з геометричними.
Складова потенціалу точки М від другого дроту і його дзеркального зображення:
Таким чином:
Потенційні коефіцієнти
Перша група формул Максвелла. Точка М може бути поміщена на поверхню першого дроту. При цьому, - відстань першого дроту до дзеркального зображення другого дроту, - відстань першого дроту до другого дроту і т. д.:
Коефіцієнти при зарядах τ1, τ2 та ін. залежать тільки від геометричних розмірів тіл, їх взаємного розташування та від властивостей середовища. Вони не залежать ні від кількості, ні від знаку зарядів і потенціалів.
Для скорочення запису останню строку та інші аналогічні їй запишемо наступним чином:
Цю систему рівнянь прийнято називати першою групою формул Максвеллу.
Тут:
Коефіцієнт . Так як і то .
Коефіцієнти називають потенційними коефіцієнтами. Розмірність їх дорівнює розмірності одиниці довжини, поділеної на фарад.
Так як у всіх коефіцієнтів під знаком логарифма стоїть дріб, чисельник якої завжди більше знаменника, то всі коефіцієнти додатні.
Коефіцієнтам може бути дано таке тлумачення: нехай заряди всіх дротів, крім першого, дорівнюють нулю
Тоді, тобто чисельно дорівнює потенціалу першого дроту, якщо на першому проводі знаходиться одиничний заряд, а заряди на інших дротах відсутні. Аналогічно чисельно дорівнює потенціалу другого дроту в тих самих умовах. Система дозволяє підрахувати потенціали заряджених тіл за відомими загальними зарядами тіл.
Для розрахунку ємності фази необхідно мати на ній максимальну фазну напругу. Тому розрахунок повинен проводитися для кожної фази окремо, при максимальній напрузі на фазі для якої проводиться розрахунок. Визначаємо напругу на інших фазах в момент, коли на розрахунковій U = UФ.
Простіше всього це зробити виходячи з векторної діаграми напруги (рис. 3). Якщо зафіксувати момент часу, то напруги на фазах в цей момент часу будуть рівні проекції на дійсну вісь. Як видно з векторної діаграми напруга на розрахунковій фазі дорівнює UФ, а на двох інших -0,5·UФ. Тому для розрахунку ємності систему рівнянь необхідно вирішити три рази змінюючи U відповідно до таблиці 1.
Рис. 3. Векторна діаграма напруги на фазах
Таблиця 1. Значення напруги на кожній з фаз ЛЕП
|
UA |
UB |
UC |
CA |
Uф |
-0,5Uф |
-0,5Uф |
CB |
-0,5Uф |
Uф |
-0,5Uф |
CC |
-0,5Uф |
-0,5Uф |
Uф |
Де Uф - значення фазної напруги.
Опис алгоритму функції СЕ, що реалізує метод дзеркальних відображень
Кафедрою передачі електричної енергії НТУ «ХПІ» метод дзеркальних відображень було реалізовано у вигляді функції користувача СЕ у програмному середовищі Mathcad, на вхід якої подається крок розщеплення фази (а), еквівалентна висота підвісу дротів (Не), відстань між осями фаз (do), число складових у фазі (n), радіус дроту поділений на 20 (ro), номінальна напруга ЛЕП (U).
Реалізація методу здійснюється у наступній послідовності (алгоритм методу):
Спочатку в програмі задається кількість складових у кожній фазі, також загальна кількість складових в двох фазах, потім загальна кількість складових в лінії. Наступний крок програми – розрахунок кута α і еквівалентного радіусу розщеплення дроту (rr). Залежно від кількості складових у фазі складається геометрична модель фазного дроту, виконується розрахунок власних і взаємних коефіцієнтів, визначається значення довжин відрізів Hi і Di, формується система алгебраїчний рівнянь (СЛАР), яка вирішується три рази при різних значеннях потенціалу на дроті. Це викликано необхідністю адаптації методу дзеркальних відображень до розрахунку трифазного поля. В результаті отримуємо вектор стовпець (К) зі значеннями зарядів, далі визначаємо величини зарядів на кожному дроті і значення максимального заряду, визначаємо значення максимального заряду і номер дроту з максимальним зарядом. Після цього виконується розрахунок ємності кожної фази, розраховується напруженість на дроті, далі розраховується напруга на дроті. Потім визначаємо значення максимальної напруженості (Emx), робочу ємність лінії, хвильовий опір (zv) і натуральну потужність (Pn). Потім задається вектор - стовпець з результатами, і задається команда виводу результатів.
На виході функції СЕ маємо вектор-стовпець з наступними результатами:
Далі наведено текст функції користувача СЕ, реалізованої у середовищі Mathcad та пояснення до неї (рис. 4).
Рис. 4. Текст функції користувача СЕ, реалізованої у середовищі Mathcad
У якості прикладу визначимо електричні параметри ЛЕП з наступними конструктивними параметрами: крок розщеплення дротів фази (а) = 40 см, еквівалентна висота підвісу дротів (Не) = 10 м, відстань між осями фаз (do) = 12 м, число складових у фазі (n) = 3, радіус дроту поділений на 20 (ro) = мм, номінальна напруга ЛЕП (U) = 500 кВ.
У результаті розрахунку було отримано наступні показники:
Отримані значення цілком відповідають параметрам реальних ЛЕП схожого конструктивного виконання.
Висновок: Запропонована функція користувача дозволяє просто, ефективно та з достатньо високою точністю проводити оціночний розрахунок електричних параметрів повітряних ліній електропередачі виходячи з їх конструктивних параметрів.
Література