The mathematical model of the system of orientation of the rotor of the wind turbine by a weathervane plane

In certain criteria for the design of low-power windturbines should consider them on the one side as an energyobject that needs most effectively use the energy of the windflow, and the other side as a high risk, requiring to ensure reliableoperation of the installation during its service life.For or...

Full description

Saved in:
Bibliographic Details
Date:2017
Main Authors: Holovko, V., Kokhanievych, V., Shykhaylov, M.
Format: Article
Language:Ukrainian
Published: Institute of Renewable Energy National Academy of Sciences of Ukraine 2017
Online Access:https://ve.org.ua/index.php/journal/article/view/38
Tags: Add Tag
No Tags, Be the first to tag this record!
Journal Title:Vidnovluvana energetika
Download file: Pdf

Institution

Vidnovluvana energetika
_version_ 1871103087609905152
author Holovko, V.
Kokhanievych, V.
Shykhaylov, M.
author_facet Holovko, V.
Kokhanievych, V.
Shykhaylov, M.
author_institution_txt_mv [ { "author": "V. Holovko", "institution": "Institute of Renewable Energy, NAS of Ukraine" }, { "author": "V. Kokhanievych", "institution": "Institute of Renewable Energy, NAS of Ukraine" }, { "author": "M. Shykhaylov", "institution": "Institute of Renewable Energy, NAS of Ukraine" } ]
author_sort Holovko, V.
baseUrl_str https://ve.org.ua/index.php/journal/oai
collection OJS
datestamp_date 2026-07-18T06:32:07Z
description In certain criteria for the design of low-power windturbines should consider them on the one side as an energyobject that needs most effectively use the energy of the windflow, and the other side as a high risk, requiring to ensure reliableoperation of the installation during its service life.For orientation of the rotors with plane orientation energylosses during orientation of rotor depends on the angular speedof orientation rotor, rotor deflection others. On the other side ofthe rotor, angular velocity orientation affects the gyroscopic,inertial and other loadings in the blade, in accordance, determinesthe lifetime of the blade and wind turbine as a whole.Today the mathematical models of the orientation systemsof the rotor in the direction of air flow, taking into account anumber of parameters and system performance orientation of therotor. It is necessary to take into account the continuousimprovement of designs of wind turbines rotor orientation,which requires amending the already proven mathematicalmodels of wind turbines rotor orientation, and some assumptionsin making mathematical models can dramatically change thepicture of the physical process.In this paper an improved mathematical model basedorientation rotor gyroscopic torque that occurs during orientationof the rotor in the direction of air flow. This mathematical modelallows obtaining an equation for calculating the velocity of theorientation of the rotor depending on the wind speeds, the angleof deflection of the rotor from the direction of air flow and anumber of system design parameters orientation of the rotor,which in turn allows to identify energy losses and gyroscopicloads on structural elements of wind turbines in the processorientation rotor.
first_indexed 2025-07-17T11:37:04Z
format Article
fulltext ВІТРОЕНЕРГЕТИКА ISSN 1819-8058 Відновлювана енергетика. 2017. № 3 48 УДК 621.548 В.М.Головко1, докт.техн.наук, В.П.Коханевич2, канд.техн.наук, М.О.Шихайлов3 (Інститут відновлюваної енергетики НАН України, Київ) Математична модель системи орієнтації ротора вітроустановки флюгерною площиною Удосконалено математичну модель системи орієнтації ротора з урахуванням гіроскопічного моменту.Отримано рівняння для розрахунку швидкостей орієнтації ротора в залежності від швидкості вітру, кута відхилення ротора від напрямку повітряного потоку та конструктивних параметрів системи орієнтації ротора. Бібл. 9, рис. 1. Ключові слова: вітроенергетика, вітроустановка, орієнтація ротора, флюгерна площина. Orcid: 10000-0003-0195-9654; 20000-0003-0033-1355; 30000-0003-1845-9904. При визначенні критеріїв для проектування вітроустановок (ВУ) малої потужності необхідно розглядати їх, з однієї сторони, як енергетичний об’єкт, що потребує максимально ефективно ви- користовувати енергію вітрового потоку, а з ін- шої – як об’єкт підвищеної небезпеки, що вима- гає забезпечення надійної роботи установки впродовж її терміну експлуатації [1, 2]. В роботі [3] для систем орієнтації роторів флюгерною площиною було показано, що енер- гетичні втрати у процесі орієнтації ротора зале- жать від кутової швидкості орієнтації, кута від- хилення ротора тощо. З іншого боку, кутова швидкість орієнтації ротора впливає на гіроско- пічні, інерційні та інші навантаження в лопаті [4, 5], що відповідно визначає термін експлуатації лопаті та вітроустановки в цілому. На сьогодні запропоновано математичні моделі систем орієнтації ротора за напрямком повітряного потоку, які враховують ряд параметрів та характе- ристик системи орієнтації ротора. При цьому необ- хідно враховувати постійне вдосконалення конс- трукцій систем орієнтації роторів ВУ, що потребує внесення змін до вже апробованих математичних моделей систем їх орієнтації. А деякі допущення, внесенні в математичні моделі, можуть кардиналь- но змінити картину фізичного процесу, що заклада- ється у вищезгадані математичні моделі системи орієнтації роторів вітроустановок. В даній роботі пропонується в математичній моделі системи орієнтації ротора вітроустановки врахувати ряд сил та моментів, що вони ство- рюють, які раніше не враховувались розробника- ми подібних моделей. Розвертання ротора за вітром при зміні його напрямку в загальному випадку описується ди- ференціальним рівнянням: 2 2Г Z dJ M dt    , (1) де ГJ – момент інерції гондоли відносно осі її повороту (осі Z); γ – кут повороту ротора; ZM – алгебраїчна сума моментів, що діють на ротор та елементи гондоли в процесі орієнта- ції ротора за напрямком повітряного потоку від- носно осі її повороту. Розглянемо схему сил, що діють на ротор та елементи гондоли в процесі орієнтації ротора за напрямком повітряного потоку (рис. 1), та про- аналізуємо сили, що входять у вираз (1). Z B Хв Р Г Тр ДM M M M M M M      , (2) де BM – поворотний момент від сили лобового ти- ску; ХвM – поворотний момент від аеродинамічних сил, що діють на флюгерну площину хвоста; РM – реактивний момент; ГM – гіроскопічний момент; ТрM – момент від сил тертя в опорно-поворотному пристрої; ДM – момент демпфування. © В.М.Головко, В.П.Коханевич, М.О.Шихайлов, 2017 ВІТРОЕНЕРГЕТИКА ISSN 1819-8058 Відновлювана енергетика. 2017. № 3 49 Рис. 1. Схема сил, що діють на ротор та елементи гондоли в процесі орієнтації ротора в напрямку повітряного потоку. Розглянемо вищенаведені моменти та сили, що привели до їх виникнення. Поворотний момент від сили лобового тиску спричиняється зміщенням центру лобового тиску відносно центру ротора при набіганні повітряно- го потоку під кутом γ до площини ротора. Дане питання було детально розглянуто в [6] та отри- мано вираз BM . Тоді, з урахування того, що в даному випадку відсутнє зміщення осі ротора відносно осі повороту, будемо мати:    2 3 3 0 04B TM k C V R cos        , де 0 – кут набігання повітряного потоку на ро- тор при зміні напрямку вітру (або кут відхилення ротора від напряму повітряного потоку при зміні напрямку вітру); ρ – питома густина повітря; TC – коефіцієнт лобового тиску; k  0,0015, якщо кут γ задається в градусах; V – швидкість повітряного потоку; R – радіус ротора. З урахуванням значень ρ, k, TC (для γ у раді- анах) вираз для поворотного моменту від сили лобового тиску набуде вигляду:    2 3 3 0 00 075BM , V R cos      . (3) Поворотний момент від аеродинамічних сил, що діють на флюгерну площину хвоста, створю- ється складовою RY від сумарної аеродинамічної сили R, яка виникає при набіганні повітряного потоку на площину хвоста під кутом  0  . Для виразу ХвM використаємо наступну залеж- ність [7]:     2 0 2 0 1 2 1 Хв C Хв Хв Хв M k F l V k sin ,                      де Ck – тангенс кута нахилу прямої СY для виб- раного типу флюгерної площини хвоста; ХвF – площа флюгерної площини хвоста; Хвl – відстань від осі обертання гондоли до центру тиску аеро- динамічних сил на флюгерну площину хвоста; Хвk – коефіцієнт зменшення швидкості повітря- ного потоку в зоні флюгерної площини внаслідок затінення її площиною ротора; χ – граничний кут, при якому ще діє ефект затінення флюгерної площини. Кут χ розраховується за таким виразом [7]: Р Хв Rarc tg l l    , де Рl – відстань від осі обертання гондоли до центру ротора. З урахуванням величини ρ запишемо вираз для поворотного моменту від аеродинамічних сил, що діють на флюгерну площину хвоста:     2 0 2 0 0 63 1 Хв C Хв Хв Хв M , k F l V k sin .                     (4) Реактивний момент наявний у конструкціях вітроустановок з вертикальним валом. Такий тип конструкцій у даний час практично відсутній, тому можна прийняти 0РM  . Гіроскопічний момент становить: ВІТРОЕНЕРГЕТИКА ISSN 1819-8058 Відновлювана енергетика. 2017. № 3 50 Г P PМ J     , (5) де JР – момент інерції ротора, вала ротора та ін- ших з’єднаних з валом ротора елементів; ωР – кутова швидкість обертання ротора; ω – кутова швидкість обертання гондоли. Введення гіроскопічного моменту в рівняння (1) є одним із принципових моментів, тоді як ряд авторів вважають його несуттєвим та ігнорують [4, 7, 8]. Момент від сил тертя в опорно- поворотному пристрої складається з наступних складових моментів:         Г X Тр Тр Тр Тр ТрB M G R M M M M M ,    (6) де  Тр B M – момент тертя від сили лобового ти- ску на ротор;   Г Тр M M – момент тертя від дії гіроскопічного моменту;  Тр G M – момент тертя від маси гондоли;   X Тр R M – момент тертя від складової RX аеродинамічної сили, що діє на флюгерну площину хвоста. При розрахунку моменту тертя від сили ло- бового тиску на ротор у загальному випадку по- трібно враховувати зміщення центру тиску у вер- тикальній площині в напрямку обертання ротора на величину εВ. Оскільки величина даного змі- щення незначна, то нею можна знехтувати. Тоді вираз для  Тр B M прийме такий вигляд:   В П В Тр В В Н НB П П l l lM B f r f r l l        , де  2 2 2 02 TB C V R cos     – сила лобового тиску на ротор; Вf та Нf – коефіцієнти тертя у верхньому та нижньому підшипниках; Вr та Нr – середні значення радіусів верхнього та нижнього підшипників; Пl – відстань між підшипниками; Вl – відстань від верхнього підшипника до осі обертання ротора. Зазвичай верхній та нижній підшипник встановлюються в опорно-поворотному при- строї одного типорозміру, тобто можна записа- ти, що В Нf f f  та В Нr r r  . З урахуван- ням даних припущень вираз для  Тр B M набу- де такого вигляду:    2 2 2 02 2 Тр TB В П П M C V R cos l lf r . l               Відповідно з урахуванням значень ρ, TC мо- мент тертя від сили лобового тиску на ротор складе:    2 2 2 0 21 74 В П Тр B П l lM , V R cos f r l           . Момент тертя від маси гондоли буде становити:  Тр G M G f r   , де G – маса гондоли, що включає масу ротора, хво- ста та інших конструктивних елементів гондоли. В загальному випадку центр маси гондоли не співпадає з центром її обертання, але у конструк- тивних схемах орієнтації з допомогою хвоста ве- личина зміщення незначна і нею нехтують. Момент тертя від дії гіроскопічного моменту складає:     2 Г Г Г Тр В В Н НM П П M MM f r f r f r l l    . З урахуванням виразу для гіроскопічного моменту вираз для моменту тертя від його дії на- буде такого вигляду:   2 Г P P Тр M П JM f r l    . Момент тертя від складової RX аеродинаміч- ної сили, що діє на флюгерну площину хвоста, незначний у порівнянні з іншими моментами тер- тя, і ним також можна знехтувати, тобто:   0 X Тр R M  . Підставимо отримані вирази для  Тр B M ,   Г Тр M M ,  Тр G M ,   X Тр R M у вираз (6) та за- пишемо в кінцевому вигляді формулу для розра- ВІТРОЕНЕРГЕТИКА ISSN 1819-8058 Відновлювана енергетика. 2017. № 3 51 хунку моменту від сил тертя в опорно- поворотному пристрої:  2 2 2 01 74 2 2 Тр В П P P П П M , V R cos f r l l J f r G f r. l l                (7) Момент демпфування у загальному випа- дку має вираз:      Д Д Д ДХв Р П M M M M   , де  Д Хв M – момент демпфування хвоста;  Д Р M – момент демпфування ротора;  Д П M – момент демпфування додаткових пристроїв. Момент демпфування хвоста знаходиться з виразу [7]:   2 30 5Д O Хв ХвХв M , k F l  , де Ok – коефіцієнт лобового опору вибраного типу флюгерної площини хвоста. Момент демпфування ротора в конструкціях вітроустановок із системами орієнтації ротора за допомогою хвоста значно менший від моменту демпфування хвоста, а тому допустимо прийняти   0Д Р M  . Також необхідно зауважити, що в даних системах орієнтації ротора використання додаткових пристроїв демпфування практично не зустрічається, і для спрощення приймемо   0Д П M  . З урахуванням прийнятих допущень та чис- лового значення ρ вираз для моменту демпфу- вання можна записати так:   2 30 63Д Д O Хв ХвХв M M , k F l  . (8) Як зазначалося вище, рівняння (1) описує рух ротора за напрямом повітряного потоку при зміні напряму вітру. Розкриємо праву частину даного диференціального рівняння, підставимо замість моментів BM , ХвM , ГM , ТрM , ДM вирази (3), (4), (6), (7), (8) та отримаємо рівняння (1) у на- ступному вигляді:           2 2 3 3 0 02 2 3 2 2 0 0 2 2 2 0 0 075 0 63 0 63 1 1 74 2 2 Г O Хв Хв C Хв Хв Хв P P В П P P П П dJ , V R cos dt , k F l , k F l V k sin J G f r , V R cos l l Jf r f r. l l                                                   (9) Дане рівняння є вихідною математичною мо- деллю, що дозволяє отримати рівняння кутової швидкості гондоли вітроустановки в залежності від швидкості повітряного потоку V, що набігає на ротор, кута відхилення ротора від напрямку повітряного потоку та ряду конструктивних па- раметрів. При цьому необхідно зауважити, що дана математична модель базується на відповід- них допущеннях. Частина цих допущень була проаналізована вище. Крім цього, необхідно обумовити інші допущення, а саме: – кутова швидкість обертання ротора P в процесі установлення ротора в напрямку повіт- ряного потоку залишається незмінною, оскільки час орієнтації ротора за напрямом повітряного потоку значно менший від часу, що обумовлює суттєве зниження обертів ротора; – коефіцієнт лобового тиску TC є незмінним і у відповідності з теорією Бетца для номінально- го режиму складає 8 9TC  ; – відстань від осі обертання гондоли до центру тиску аеродинамічних сил на флюгерну площину хвоста залишається незмінною. У зага- льному випадку центр тиску зміщується, але ве- личина даного зміщення складає декілька відсот- ків відстані від осі обертання гондоли до центру тиску аеродинамічних сил, тобто даним зміщен- ням можна знехтувати. Для отримання аналітичного виразу для ку- тової швидкості орієнтації ротора за напрямком повітряного потоку при зміні напрямку вітру ви- користаємо таке перетворення [9]: 2 2 2 1 2 d d ddt      . (10) ВІТРОЕНЕРГЕТИКА ISSN 1819-8058 Відновлювана енергетика. 2017. № 3 52 З урахуванням виразу (10) запишемо дифе- ренціальне рівняння (9) у наступному вигляді:           2 3 0 0 2 0 0 2 2 0 2 1 B G Г Г B T Г T Хв Хв Д Т J d C cos C d cos C C C k sin C C ,                                            (11) де 2 30 075BC , V R ; 2 2 21 74 B В П T П l lC , V R f r; l        Г P PC J  ; GTC G f r ; 2 Г P P Т П JC f r; l   20 63Хв C Хв ХвC , k F l V ; 30 63Д O Хв ХвC , k F l . Для отримання рівняння кутової швидкості орієнтації ротора в напрямку повітряного потоку при зміні напрямку вітру проінтегруємо вираз (11) та отримаємо: 2 2 0 1 2 Г B B G Г B B Хв Хв Г Д Т T T T J C D C D C C C C D C C ,                (12) де           0 0 0 0 1 33 36 4 1 31 3 12 4 BD cos cos sin sin ;                         0 1 1 2 2 4BTD sin     ;         2 2 2 0 0 0 2 2 02 2 0 2 2 0 1 1 4 2 2 1 2 8 2 1 2 4 Хв Хв Хв Хв Хв Хв Хв Хв D k k cos k k cos sin k k k sin .                                                                               Для визначення вільного члена можливі два варіанти: – початок відліку кута γ співпадає з віссю Х, тоді ω=0 при γ=0; – початок відліку кута γ співпадає з напрям- ком швидкості вітру, тоді відповідно ω=0 при γ=γ0. Оскільки гондола вітроустановки починає свій рух від осі Х, то перший варіант ближче від- повідає фізичному процесу, що проходить у сис- темі, тоді вираз для С0 набуде такого вигляду:      0 0 00B BT T B B Хв ХвC C D C D C D   , де   00 1 2 4BTD sin   ;   0 00 0 0 1 33 36 4 1 33 12 4 BD cos cos sin sin ;                2 0 0 00 2 2 2 0 02 2 1 2 2 1 2 2 8 Хв Хв Хв Хв Хв D k k cos k cos k sin .                                        При цьому необхідно зауважити, що коефіці- єнти BC , BTC та ХвC розраховуються для конк- ретної величини швидкості вітру, а коефіцієнти BD , BTD , ХвD та  0BD ,  0BTD ,  0ХвD для конкретного відхилення ротора від напрямку по- вітряного потоку. З рівняння (12) можна отримати аналітичний вираз для ω для розрахунку втрат при орієнтації ротора за напрямком повітряного потоку. Оскі- льки даний вираз буде надто громіздким, тому більш доцільним є розрахунок числових значень ω з певним інтервалом кута γ і подальший розра- хунок втрат за виразом, наведеним в [3], та гіро- скопічних навантажень у лопатях та інших еле- ментах конструкції вітроустановки [4, 5]. Висновки. 1. Удосконалено математичну модель орієнтації ротора з урахуванням гіроско- пічного моменту, що виникає у процесі орієнтації ротора за напрямком повітряного потоку. 2. Отримано рівняння для розрахунку швид- костей орієнтації ротора в залежності від швид- костей вітру, кута відхилення ротора від напрям- ку повітряного потоку і ряду конструктивних па- раметрів системи орієнтації ротора, що дозволяє визначити енергетичні втрати та гіроскопічні на- ВІТРОЕНЕРГЕТИКА ISSN 1819-8058 Відновлювана енергетика. 2017. № 3 53 вантаження на елементи конструкції вітро- установки в процесі орієнтації ротора. 1. Кудря С.О., Пепелов О.В., Ціцікан Р.М. Сучасний стан розвитку вітроенергетики в світі // Матеріали ХI Між- народної науково-практичної конференції «Відновлювана енергетика ХХІ соліття». – АР Крим, смт. Миколаївка, 13 – 17 вересня 2010 року. – С. 118 – 120. 2. Кудря С.А., Кузнецов Н.П., Сахно Б.Г. Опыт использо- вания ветровых електрических установок в Украине // Альтер- нативная энергетика и экология. – 2011. – № 8. – С. 50 – 53. 3. Головко В.М., Коханєвич В.П., Шихайлов М.О. Ана- ліз систем орієнтації ротора вітроустановки малої потужно- сті // Відновлювана енергетика. – 2015. – № 2. – С. 55 – 60. 4. Фатеев Е. М. Ветродвигатели. – Л. : ВИМЭ, 1946. – 244 с. 5. GERMANISCHER LLOYD: Vorschriften und Richtli- nien, Teil 1 Windenergie, Richtlinien fur die Zertifizierung von Windenergieanlagen. – K. 1-10. – 1993 mit Ergдnzung 1994 (German). 6. Головко В.М., Коханєвич В.П., Шихайлов М.О., Мар- ченко Н.В. Вплив параметрів системи орієнтації ротора з використанням конструктивної схеми підпружиненого хво- ста на статичні характеристики вітроустановки // Віднов- лювана енергетика. – 2015. – № 3. – С. 30 – 39. 7. Мосалев В.Ф. Об изгибе упругой лопасти гироско- пическим моментом при ориентации поворотной головки ветродвигателя на ветер хвостом // Труды ВНИИЭМ. – 1970. – Т. № 34. – С. 78 – 93. 8. Кузьо І.В, Корендій В.М. Динаміка систем орієнтації вітроколіс // Вісник Нац. ун-ту «Львівська політехніка». – 2012. – № 730: Динаміка, міцність та проектування машин і приладів. – С. 51 – 57. 9. Федотов В.Е. Исследование работы ветродвигателя с ветроколесом за башней, кандидатская диссертация, ЭНИН АН СССР, 1953. REFERENCES 1. Kudrya S.О., Pepelov O.V., Tsitsikan R.M. The current state of development of wind power in the world // Proceedings of the XI International Scientific Conference «Renewable energy ХХІ century». _ Crimea, smt. Mykolayivka, 13 – 17 September 2010. – Р. 118 – 120. 2. Kudrya S.О., Kuznetsov М.P., Sakhno B.G. Experience of using wind power plants in Ukraine // Alternative energy and ecology. – 2011. – № 8. – Р. 50 – 53. 3. Golovko V.M., Kokhanievych V.P., Shykhailov M.O. Rotor orientation system analysis of a low capacity wind turbine // Renewable Energy. – 2015. – № 2. – P. 55 – 60. 4. Fateev E.M. Wind turbines. – L.: VIME, 1946. – 244 p. 5. GERMAN LLOYD: Regulations and directives, Part 1 Wind energy, Guidelines for the certification of wind energy installations. – K. 1 – 10. – 1993 with addition 1994 (German). 6. Golovko V.M., Kokhanievych V.P., Shykhailov M.O., Marchenko N.V. The effect of rotor orientation system parameters with a spring-loaded tail constructive scheme over static characteristics of the wind turbines // Renewable Energy. – 2015. – № 3. – P. 30 – 39. 7. Mosaliev V.F. On the bending of an elastic blade with a gyroscopic moment when the rotating head of a windmill is oriented to the wind with a tail // Proceedings of VNIIEM. – 1970. – T. No. 34. – P. 78 – 93. 8. Cousin I.V., Korendiy V.M. Dynamics of orientation wind wheels // Journal of Nat. Univ «Lviv Polytechnic». – 2012. – № 730: Dynamics, strength and design of machines and devices. – P. 51 – 57. 9. Fedotov V.E. Research of wind turbine operation with a wind wheel behind the tower, Ph.D. thesis, ENIN of the USSR Academy of Sciences, 1953. В.М.Головко, докт.техн.наук, В.П.Коханевич, канд.техн.наук, Н.А.Шихайлов (Институт возобновляемой энергетики НАН Украины, Киев) Математическая модель системы ориентации ротора ветроустановки флюгерной плоскостью Усовершенствована математическая модель системы ориентации ротора с учетом гироскопического момента. Получены уравнения для расчета скоростей ориентирования ротора в зависимости от скорости ветра, угла отклонения ротора от направления воздушного потока и конструкционных параметров системы ориентации ротора. Библ. 9, рис. 1. Ключевые слова: ветроэнергетика, ветроустановка, ориентация ротора, флюгерная плоскость. Holovko V., prof., dokt.tekhn.nauk, Kokhanievych V., kand.tekhn.nauk, Shykhaylov M. (Institute of Renewable Energy National Academy of Sciences of Ukraine, Kyiv) The mathematical model of the system of orientation of the rotor of the wind turbine by a weathervane plane The mathematical model of the rotor orientation system is improved with allowance for the gyroscopic moment and equations are obtained for calculating the rotor orientation speeds as a function of wind speed, the angle of the rotor deviation from the direction of the air flow, and the structural parameters of the rotor orientation system. References 9, figures 1. Keywords: wind power, wind turbine, rotor orientation, weathervane plane. SYNOPSES In certain criteria for the design of low-power wind turbines should consider them on the one side as an energy object that needs most effectively use the energy of the wind flow, and the other side as a high risk, requiring ensure reliable operation of the installation during its service life . For orientation of the rotors with plane orientation energy losses during orientation of rotor depends on the angular speed of orientation rotor, rotor deflection others. On the other side of the rotor angular velocity orientation affects the gyroscopic, ВІТРОЕНЕРГЕТИКА ISSN 1819-8058 Відновлювана енергетика. 2017. № 3 54 inertial and other loading in the blade, in accordance determines the lifetime of the blade and wind turbine as a whole. Today the mathematical models of the orientation systems of the rotor in the direction of air flow, taking into account a number of parameters and system performance orientation of the rotor. It is necessary to take into account the continuous improvement of designs of wind turbines rotor orientation, which requires amending the already proven mathematical models of wind turbines rotor orientation, and some assumptions in making mathematical models can dramatically change the picture of the physical process. In this paper an improved mathematical model based orientation rotor gyroscopic torque that occurs during orientation of the rotor in the direction of air flow. This mathematical model allows to obtain an equation for calculating the velocity of the orientation of the rotor depending on the wind speeds, the angle of deflection of the rotor from the direction of air flow and a number of system design parameters orientation of the rotor, which in turn allows to identify energy losses and gyroscopic loads on structural elements of wind turbines in process orientation rotor. Стаття надійшла до редакції 23.05.17 Остаточна версія 23.08.17 УДК 621.519.2 П.Ф.Васько1, докт.техн.наук, А.П.Вербовий2, канд.техн.наук, С.Т.Пазич3 (Інститут відновлюваної енергетики НАН України, Київ) Реалізація стохастичної двопараметричної моделі поздовжньої складової швидкості вітру для задач вітроенергетики Розроблено реалізацію стохастичної моделі поздовжньої складової швидкості вітру для використання в сучасних математичних пакетах. Моделювання поривів швидкості вітру здійснено шляхом генерування двох послідовностей випадкових величин. Розподіл амплітуд поривів підпорядкований нормальному закону, а тривалості поривів – закону Пуассона. Бібл. 24, табл. 2, рис. 5. Ключові слова: вітротурбіна, вітер, випадкова величина, модель, порив вітру, швидкість вітру. Orcid: 10000-0001-8807-7173; 20000-0003-2838-6032; 30000-0002-6310-9842 Вступ. При розв’язанні науково-технічних задач, пов’язаних із проектуванням, дослідженням робочих режимів, оптимізацією параметрів, техніко-економічним обґрунтуванням використання вітроенергетичних установок (ВЕУ) різного призначення широкого вжитку набуло математичне моделювання процесів перетворення енергії вітрового потоку в механічну енергію обертового руху, теплову та електричну енергію. Дослідження режимів роботи складових частин ВЕУ та установки в цілому потребує застосування адекватних моделей зміни швидкості вітру в часі відповідно до постановки задачі. Динаміка механічних процесів ВЕУ моделюється нелінійним диференціальним рівнянням руху обертових частин [1, 2]: ( ( ), , ) ( ),c dJ M v t M dt      (1) де  – кутова швидкість обертання валу вітротурбіни; J – сумарний момент інерції всіх обертових мас, приведений до валу вітротурбіни; t – поточний час; ( )v t – залежність швидкості вітру від часу;  – кут повороту лопаті вітротурбіни; ( ( ), , )M v t   – момент, що розвиває вітротурбіна; ( )cM  – сумарний момент опору від сил тертя і навантаження. Момент, що розвиває вітротурбіна, істотно залежить від швидкості вітру, тому точність, з якою задається ( ),v t впливає на результати мо- © П.Ф.Васько, А.П.Вербовий, С.Т.Пазич, 2017
id veorgua-article-38
institution Vidnovluvana energetika
keywords_txt_mv keywords
language Ukrainian
last_indexed 2026-07-19T01:02:37Z
publishDate 2017
publisher Institute of Renewable Energy National Academy of Sciences of Ukraine
record_format ojs
resource_txt_mv veorgua/3b/f4b4bcfd37c67ceb75cdc74caf3ae93b.pdf
spelling veorgua-article-382026-07-18T06:32:07Z The mathematical model of the system of orientation of the rotor of the wind turbine by a weathervane plane Математическая модель системы ориентации ротора ветроустановки флюгерной плоскостью Математична модель системи орієнтації ротора вітроустановки флюгерною площиною Holovko, V. Kokhanievych, V. Shykhaylov, M. In certain criteria for the design of low-power windturbines should consider them on the one side as an energyobject that needs most effectively use the energy of the windflow, and the other side as a high risk, requiring to ensure reliableoperation of the installation during its service life.For orientation of the rotors with plane orientation energylosses during orientation of rotor depends on the angular speedof orientation rotor, rotor deflection others. On the other side ofthe rotor, angular velocity orientation affects the gyroscopic,inertial and other loadings in the blade, in accordance, determinesthe lifetime of the blade and wind turbine as a whole.Today the mathematical models of the orientation systemsof the rotor in the direction of air flow, taking into account anumber of parameters and system performance orientation of therotor. It is necessary to take into account the continuousimprovement of designs of wind turbines rotor orientation,which requires amending the already proven mathematicalmodels of wind turbines rotor orientation, and some assumptionsin making mathematical models can dramatically change thepicture of the physical process.In this paper an improved mathematical model basedorientation rotor gyroscopic torque that occurs during orientationof the rotor in the direction of air flow. This mathematical modelallows obtaining an equation for calculating the velocity of theorientation of the rotor depending on the wind speeds, the angleof deflection of the rotor from the direction of air flow and anumber of system design parameters orientation of the rotor,which in turn allows to identify energy losses and gyroscopicloads on structural elements of wind turbines in the processorientation rotor. Усовершенствована математическая модель системыориентации ротора с учетом гироскопическогомомента. Получены уравнения для расчета скоростейориентирования ротора в зависимости от скоростиветра, угла отклонения ротора от направлениявоздушного потока и конструкционных параметровсистемы ориентации ротора. Удосконалено математичну модель системи орієнтації ротора з урахуванням гіроскопічного моменту.Отримано рівняннядля розрахунку швидкостей орієнтації ротора в залежності від швидкості вітру, кута відхилення ротора від напрямкуповітряного потоку та конструктивних параметрів системи орієнтації ротора. Institute of Renewable Energy National Academy of Sciences of Ukraine 2017-09-11 Article Article application/pdf https://ve.org.ua/index.php/journal/article/view/38 Vidnovluvana energetika ; No. 3 (50) (2017): Scientific and Applied Journal Vidnovluvana energetika; 48-54 Возобновляемая энергетика; ##issue.no## 3 (50) (2017): Научно-прикладной журнал Возобновляемая энергетика; 48-54 Відновлювана енергетика; № 3 (50) (2017): Науково-прикладний журнал Відновлювана енергетика; 48-54 2664-8172 1819-8058 uk https://ve.org.ua/index.php/journal/article/view/38/28 Copyright (c) 2017 V. Holovko, V. Kokhanievych, M. Shykhaylov https://creativecommons.org/licenses/by-nc-nd/4.0
spellingShingle Holovko, V.
Kokhanievych, V.
Shykhaylov, M.
The mathematical model of the system of orientation of the rotor of the wind turbine by a weathervane plane
title The mathematical model of the system of orientation of the rotor of the wind turbine by a weathervane plane
title_alt Математическая модель системы ориентации ротора ветроустановки флюгерной плоскостью
Математична модель системи орієнтації ротора вітроустановки флюгерною площиною
title_full The mathematical model of the system of orientation of the rotor of the wind turbine by a weathervane plane
title_fullStr The mathematical model of the system of orientation of the rotor of the wind turbine by a weathervane plane
title_full_unstemmed The mathematical model of the system of orientation of the rotor of the wind turbine by a weathervane plane
title_short The mathematical model of the system of orientation of the rotor of the wind turbine by a weathervane plane
title_sort mathematical model of the system of orientation of the rotor of the wind turbine by a weathervane plane
url https://ve.org.ua/index.php/journal/article/view/38
work_keys_str_mv AT holovkov themathematicalmodelofthesystemoforientationoftherotorofthewindturbinebyaweathervaneplane
AT kokhanievychv themathematicalmodelofthesystemoforientationoftherotorofthewindturbinebyaweathervaneplane
AT shykhaylovm themathematicalmodelofthesystemoforientationoftherotorofthewindturbinebyaweathervaneplane
AT holovkov matematičeskaâmodelʹsistemyorientaciirotoravetroustanovkiflûgernojploskostʹû
AT kokhanievychv matematičeskaâmodelʹsistemyorientaciirotoravetroustanovkiflûgernojploskostʹû
AT shykhaylovm matematičeskaâmodelʹsistemyorientaciirotoravetroustanovkiflûgernojploskostʹû
AT holovkov matematičnamodelʹsistemioríêntacíírotoravítroustanovkiflûgernoûploŝinoû
AT kokhanievychv matematičnamodelʹsistemioríêntacíírotoravítroustanovkiflûgernoûploŝinoû
AT shykhaylovm matematičnamodelʹsistemioríêntacíírotoravítroustanovkiflûgernoûploŝinoû
AT holovkov mathematicalmodelofthesystemoforientationoftherotorofthewindturbinebyaweathervaneplane
AT kokhanievychv mathematicalmodelofthesystemoforientationoftherotorofthewindturbinebyaweathervaneplane
AT shykhaylovm mathematicalmodelofthesystemoforientationoftherotorofthewindturbinebyaweathervaneplane