RENEWABLE ENERGY
СОЛНЕЧНАЯ ЭНЕРГЕТИКА
SOLAR ENERGY
СОЛНЕЧНЫЕ ГОРОДА
SOLAR CITIES
Статья поступила в редакцию 31.08.15. Ред. рег. № 2331
The article has entered in publishing office 31.08.15. Ed. reg. No. 2331
УДК 62-523.2
ЭНЕРГОЭФФЕКТИВНОЕ ТЕРМИНАЛЬНОЕ УПРАВЛЕНИЕ СЕРВОПРИВОДАМИ СОЛНЕЧНЫХ МОДУЛЕЙ
Е.В. Капля
Филиал НИУ «МЭИ» в г. Волжском 404110 г. Волжский, Волгоградская обл., пр. Ленина, д. 69 Тел.: (8443) 210160, e-mail: ev-kaple@yandex.ru
Заключение совета рецензентов: 03.09.15 Заключение совета экспертов: 06.09.15 Принято к публикации: 09.09.15
Получена новая формула, позволяющая вычислить оптимальную продолжительность поворота солнечного модуля. Установлена зависимость оптимальной по расходу энергии продолжительности поворота солнечного модуля от величины углового перемещения и параметров сервопривода. Увеличение углового перемещения приводит к увеличению оптимальной продолжительности поворота. Представлена новая структурная схема автоматизированного сервопривода модуля с блоком терминального управления. Для поворота солнечного модуля использован сервопривод с двигателем постоянного тока, червячной передачей и датчиком углового положения. Построены графики зависимости расхода энергии при терминальном управлении серводвигателем модуля от продолжительности поворота модуля. Оптимальная продолжительность поворота модуля соответствует минимуму энергетических затрат на процесс поворота.
Ключевые слова: солнечный модуль, солнечная электростанция, оптимальное потребление энергии, терминальное управление, оптимальная продолжительность, сервопривод, серводвигатель, двигатель постоянного тока.
ENERGY EFFICIENT TERMINAL CONTROL BY THE SOLAR SERVO MODULES
E.V. Kaplya
Volzhsky Branch of the National Research University «Moscow Power Engineering Institute» 69 Lenin str., Volzhsky, Volgograd reg., 404110, Russia Tel.: (8443) 210160, e-mail: ev-kaple@yandex.ru
Referred: 03.09.15 Expertise: 06.09.15 Accepted: 09.09.15
The new formula that allows to compute the optimal duration of rotation of the solar module. The dependence of optimal consumption duration of rotation of the solar module on the magnitude of angular displacement and servo parameters. The increase in the angular displacement increases, the optimal duration of the rotation. Presents a new structural scheme of an automated servo module with the terminal control unit. To rotate the solar module used a servo with a DC motor, worm gear and a sensor of angular position. Plots of power consumption at terminal servo motor control module on the duration of rotation of the module. The optimal duration of rotation of the module corresponds to the minimum energy consumption for the turning process.
Keywords: solar module, solar power plant, optimal consumption energy, terminal control, fixed-time control, optimal termination time, optimal duration, servo-driver, servo motor, DC motor.
Егор Викторович Капля Egor V. Kaplya
Сведения об авторе: канд. физ.-мат. наук, доцент кафедры «Автоматизация технологических процессов и производств» ВФ МЭИ, ведущий научный сотрудник лаборатории диагностики и прогнозирования.
Образование: Волгоградский государственный технический университет (2000).
Область научных интересов: энергетика, автоматизация.
Публикации: 62.
Information about the author: candidate of physical and mathematical sciences, docent faculty «Automation of technological processes and productions» of the Volgzhky branch of Moscow power institute, leading researcher in laboratory of diagnostics and prediction of the Volgzhky branch of Moscow power institute.
Education: Volgograd State Technical University (2000).
Research area: power engineering, automation.
Publications: 62.
Введение
Современные солнечные электростанции (СЭС) -сложные автоматизированные объекты энергетики. Повышение энергоэффективности СЭС достигается совершенствованием фотоприемных устройств, использованием автоматических систем позиционирования и оптимизацией процессов управления системами позиционирования модулей [1-4]. Создание устойчивых, надежных и энергетически экономичных автоматизированных систем позиционирования солнечных модулей является перспективной задачей.
Автоматическое позиционирование фотоэлектрических модулей, гелиостатов и параболоцилиндриче-ских концентраторов обычно осуществляют с помощью электродвигателей. В современных СЭС каждый модуль имеет индивидуальный сервопривод с контроллером (управляющим устройством).
Редуктор с червячной передачей позволяет исключить возможность прокручивания солнечного модуля под действием воздушного потока и собственного веса модуля. Для поворота солнечных модулей можно использовать сервопривод с двигателем постоянного тока (ДПТ), червячной передачей и датчиком углового положения. При угловом позиционировании солнечных модулей с помощью ДПТ целесообразно использовать систему терминального управления. Терминальное управление сервоприводом реализуется на конечном отрезке времени [5-7]. При заданных энергетических ресурсах исполнительных элементов системы минимально допустимое значение Т продолжительности поворота солнечного модуля зависит от начальных и конечных условий. От продолжительности поворота зависят требуемые значения сил, осуществляющих движение, и количество затраченной энергии. Поиск аналитических выражений, позволяющих определить оптимальное значение Т, - актуальная задача, которую целесообразно решить для оптимизации систем позиционирования.
Цели исследования:
1. Анализ зависимости энергии, затраченной сервоприводом на заданное угловое перемещение солнечного модуля, от продолжительности поворота модуля.
2. Формулировка выражения, определяющего оптимальную (по расходу электроэнергии) продолжительность поворота солнечного модуля.
Геометрия задачи поворота солнечного модуля
Поворот солнечного модуля вокруг горизонтальной и вертикальной осей (рис. 1) предполагает использование двух серводвигателей. Алгоритмы управления двигателями аналогичны, поэтому рассмотрим систему позиционирования фотоэлектрического модуля при повороте вокруг одной горизонтальной оси ОХ.
Рис. 1. Система позиционирования фотоэлектрического солнечного модуля Fig. 1. The positioning system photovoltaic solar module
Начальное положение модуля характеризуется углом 60 между вектором нормали к поверхности модуля и горизонтальной плоскостью. Необходимо повернуть модуль на угол Д6 с минимальными затратами энергии и остановить его в положении с вектором нормали N, направленным под углом 6 относительно горизонтальной плоскости (рис. 1).
Структурная схема сервопривода солнечного модуля
Структурная схема автоматизированного сервопривода с ДПТ и блоком терминального управления (БТУ) представлена на рис. 2.
На схеме использованы следующие обозначения: 6(0 - угол поворота солнечного модуля; а(0 -угол поворота якоря ДПТ; ю(0 - угловая скорость
вала ДПТ; ф(0 - угловая скорость солнечного модуля; £(0 - угловое ускорение солнечного модуля; е(0 - угловое ускорение вала ДПТ; - управляющее напряжение; i(t) - ток в цепи якоря ДПТ; L - индуктивность цепи якоря; R - активное сопротивление цепи якоря; E(t) = - противодействующая ЭДС якоря; ^ - коэффициент противоЭДС, связывающий E(t) с угловой скоростью; Mвр(t) = = ^Щ) - вращающий момент, создаваемый ДПТ; ^ - коэффициент пропорциональности, связывающий ток в цепи якоря и развиваемый двигателем вращающий момент; Mc(t) - момент сопротивления нагрузки, приведенный к валу двигателя; J - суммарный осевой момент инерции якоря ДПТ, колес редуктора и солнечного модуля; п - передаточное число редуктора сервопривода.
Рис. 2. Структурная схема автоматизированного сервопривода с блоком терминального управления Fig. 2. The structural diagram of automated servomotor with the terminal control unit
Жесткая механическая связь вала двигателя и ведомого вала нагрузки посредством зубчатых колес редуктора позволяет связать фазовые координаты вала двигателя и вала модуля:
6(0 = а(0Л?; ф(0 = ю(0/п; £(0 = Е(0/П.
Измерение угла 6(0 ориентации солнечного модуля осуществляется с помощью энкодера (датчика углового положения). Мгновенная угловая скорость ф(0 и угловое ускорение £(0 солнечного модуля вычисляются контроллером в результате дифференцирования функции 6(0 по времени. Численное дифференцирование по двум узлам может привести к некорректным результатам. Поэтому для вычисления производных целесообразно использовать многоузловую разностную схему, которая выполняет функцию сглаживающего фильтра низких частот для измеренных значений 6изм(0.
На вход БТУ подаются заданные конечные значения фазовых координат вала солнечного модуля: 6зад & + ^ = 6, ¥зад ^ + ^ = 4, £зад ^ + ^ = £, где Т0 - момент времени начала поворота; Т - продолжительность поворота. Модуль переводится из одного неподвижного состояния в другое неподвижное состояние, поэтому начальную и конечную угловую скорость солнечного модуля примем равными нулю: у ^) = у = 0 . Текущие значения фазовых координат сравниваются с заданными. БТУ вычисляет и формирует управляющий сигнал ^0 на основе выбранного закона терминального управления.
Терминальное управление сервоприводом
Известный [5, 6] закон терминального управления при трех конечных условиях применительно к сервоприводу представим в следующем виде:
'(', T) =
Л
k.,n
12(9-9(0) 6 ( + V(0)
T -1
+ I
T -t) при t e[T0; T0 + T]; 0 - в остальных случаях,
(1)
0
Q(T)= J и (t,T)i (t,T)dt .
(5)
где Т = Т0 + Т + АТ - момент времени окончания
переходного процесса; АТ - жесткость управления (временной интервал между ведущей и ведомой фазовыми точками); АТ = Т; ки - коэффициент пропорциональности между напряжением на обмотке ДПТ и угловым ускорением ротора; ки = Ж / кт .
Ток в цепи якоря ДПТ представим решением (3) относительно /(/, Т) для / > Т0:
/ (, Т ) = ехр Г-| (-Т. ))х
х 7 J ехр f Т To)
<(', T)- £0
d а(', T)
dt'
dt' +
Т
+ exp|-7(t - To )|i (To,T),
(6)
Математическая модель сервопривода постоянного тока
В основе известной [8] математической модели сервопривода постоянного тока - система дифференциальных уравнений:
J
d 2a(t)
dt
2 = kmi(t) - M c(t);
d a(t)
di(t)
L-^- + Ri(t) = и (t) - k dt
dt
(2)
(3)
Момент сил сопротивления будем считать линейно зависящим от угловой скорости вращения ротора ДПТ:
Х1ЦО + Х0, ю()> 0;
%1ю(()-%о, ю()< 0;
М () = Мвр (), если ю () = 0 Л |Мвр ()| < Хо; (4)
Хо, ю() = 0 л МвР ()>Х0;
^ ю() = 0 Л мвр () <
где х0 - момент трогания ротора двигателя и механической нагрузки сервопривода, приведенный к валу двигателя; х1 - коэффициент сопротивления движению.
Уравнения (2)-(4) соответствуют структурной схеме на рис. 2.
Расход энергии при терминальном управлении
Основной энергетической характеристикой процесса управления сервоприводом является величина электрической энергии Q, затраченной на питание двигателя в процессе управления (от момента Т0 до момента Т0 + Т конца процесса управления).
Энергия, потребляемая якорем ДПТ от источника электрической энергии в процессе терминального управления сервоприводом, выражается интегралом
где /(Т0, Т) - ток в цепи якоря в момент времени Т0. При пуске двигателя можно принять 1(Т0, Т) = 0, т.к. ток через индуктивный элемент не может измениться мгновенно.
1000
800
600
400
200
30
20е
10
Q(V. \ Дж 1 5° ^0^____
Т. с
0 5 Д6. рад -0.5 10 1 5 20 25 до/
-0,4
-0.3
0.2
-0,1
Тт.С
10
15
20
25
Рис. 3. Графики зависимостей Q(T) и Тт(Д0) Fig. 3. Dependency graphs Q(T) and Тт(Д0)
Функция Q(T) имеет минимум (рис. 3). Оптимальная продолжительность Tm поворота модуля соответствует минимуму целевой функции Q(T): Q (Tm )= min Q (T), где Q(Tm) - минимальная энер-
V ' 0<T
гия, которую необходимо затратить на питание двигателя в процессе поворота модуля.
Убывание функции Q(T) в области 0 < T < Tm связано с уменьшением максимума напряжения
max \u (t,T )| при увеличении T. Возрастание
T0 <t <T0 +T1 4 'I
функции Q(T) при T > Tm обусловлено ростом продолжительности простоя ротора и увеличением энергии, затраченной за время простоя [7]. Величина Tm зависит от заданных начальных и конечных условий, а также от характеристик сервопривода.
Оптимальная продолжительность поворота
В начале процесса терминального управления при T = Tm управляющее напряжение примерно равно |u (To, Tm )| = 2Up, где Up = Rx J km - пусковое
напряжение ДПТ. Полярность напряжения в момент пуска ДПТ определяется знаком углового перемещения: sign[u(T0,T)] = sign(де), где A9 = 9-9(T0), следовательно
sign(A9). u {T0,Tm ) = 2Up .
(7)
Уравнение (7) для (^(T0) = \¡/ = o)n(| = о) с учетом (1) принимает вид
sign (A0)- kun-
6A9
,2 = UP .
(Tm +AT)
Решив (8) относительно Tm, получаем формулу
(8)
Тт = ^ 6Д9| "АТ ^ 6Д9| -ДТ . (9)
Формула (9) позволяет вычислить приближенное значение Тт оптимальной продолжительности поворота модуля для заданного угла поворота Д6, представленного в радианах. Увеличение углового перемещения приводит к увеличению оптимальной продолжительности процесса терминального управления.
Результаты моделирования
Численные эксперименты проведены с моделью сервопривода, обладающего характеристиками: Ь = = 0,01 Гн; R = 2 Ом; ^ = 0,08 Н-м/А; = 0,2 В-с/рад; J = 1,5 кг-м2; п = 10; х = 0,1 Н-м-с/рад; х0 = 0,2 Н-м при следующих начальных и конечных условиях: у(Т>) = ю(Т0 = 0; 6(Т0) = 20°; 6 = 6(Т0) + Д6; ф = ю = 0; £ = е = 0.
При моделировании были приняты: ^ = 0; Т0 = 1 с; ДТ = 0,01 с; Дt = 10-3 с.
Угловое перемещение солнечного модуля варьировалось в диапазоне 0°< Д6 < 34°.
Графики зависимости энергии Q(T), затраченной в процессе поворота, от продолжительности Т поворота построены с помощью (1)-(6) при Д6 = 5, 10, 15, 20, 25° и показаны на рис. 3.
В результате минимизации функции Q(T) для 0°< Д6 < 34° построен график зависимости Тт(Д6), представленный в нижней части рис. 3 сплошной линией. Результаты расчета Тт(Д6) по (9) представлены на рис. 3 наклонной пунктирной кривой. Близость сплошной и пунктирной расчетных кривых Тт(Д6) подтверждает корректность расчета оптимальной продолжительности поворота модуля с помощью (9).
Заключение
Получена новая формула для расчета оптимальной (по расходу энергии) продолжительности поворота солнечного модуля: Тт = ^ 6(Уп/х0)|Д0| -ДТ,
где Д6 - заданное угловое перемещение солнечного модуля (в радианах); п - передаточное число редуктора сервопривода; У - суммарный осевой момент инерции вращающихся деталей сервопривода и солнечного модуля; %о - момент трогания сервопривода; ДТ - жесткость управления. Предлагаемая формула позволяет вычислить оптимальную продолжительность Тт процесса поворота солнечного модуля из одного неподвижного состояния в другое.
- с-
Ж
%
Список литературы
1. Виссарионов В.И., Дерюгина Г.В., Кузнецова В. А., Малинин Н.К. Солнечная энергетика. М.: Издательский дом МЭИ, 2008.
2. Капля Е.В., Кузеванов В.С., Шевчук В.П. Моделирование процессов управления в интеллектуальных измерительных системах. М.: ФИЗМАТЛИТ, 2009.
3. Camacho E.F., Berenguel M., Rubio F.R., Martinez D. Advances in industrial control: Control of solar energy systems. Springer-Verlag London Limited, 2012.
References
1. Vissarionov V.I., Derugina G.V., Kuznecova V.A., Malinin N.K. Solnecnaa energetika. M.: Izdatel'-skij dom MEI, 2008.
2. Kapla E.V., Kuzevanov V.S., Sevcuk V.P. Modelirovanie processov upravlenia v intellektual'nyh izmeritel'nyh sistemah. M.: FIZMATLIT, 2009.
3. Camacho E.F., Berenguel M., Rubio F.R., Martinez D. Advances in industrial control: Control of solar energy systems. Springer-Verlag London Limited, 2012.
4. Ochieng R.M. Solar collectors and panels, theory and applications. Croatia, Rijeka: Sciyo, 2010.
5. Батенко А.П. Системы терминального управления. М.: Радио и связь, 1984.
6. Батенко А.П. Управление конечным состоянием движущихся объектов. М.: Советское радио, 1977.
7. Капля Е.В. Энергоэффективное терминальное управление серводвигателем постоянного тока // Автоматизация процессов управления. 2015. № 1 (39). С. 27-33.
8. Фираго Б.И., Павлячик Л.Б. Теория электропривода. Минск: ЗАО «Техноперспектива», 2004.
4. Ochieng R.M. Solar collectors and panels, theory and applications. Croatia, Rijeka: Sciyo, 2010.
5. Batenko A.P. Sistemy terminal'nogo upravlenia. M.: Radio i svaz', 1984.
6. Batenko A.P. Upravlenie konecnym sostoaniem dvizusihsa ob"ektov. M.: Sovetskoe radio, 1977.
7. Kapla E.V. Energoeffektivnoe terminal'noe upravlenie servodvigatelem postoannogo toka // Avtomatizacia processov upravlenia. 2015. № 1 (39). S. 27-33.
8. Firago B.I., Pavlacik L.B. Teoria elektroprivoda. Minsk: ZAO «Tehnoperspektiva», 2004.
Транслитерация по ISO 9:1995