Периодическое изменение скорости ветра при постоянном моменте нагрузки
В действительности скорость ветрового потока в результате перемещений воздуха, связанных как с общими атмосферными условиями, так и с вихреобразова - ниями, вызванными местными условиями у поверхности
земли, подвержена непрерывным пульсациям— порывам. Структурная характеристика порывов ветра чрезвычайно разнообразна, но, несмотря на это, в аэрологии установлены некоторые приближенные закономерности, характеризующие ход порывов во времени.
Для уяснения характера движения системы при пульсациях скорости ветра рассмотрим случай, когда среднее значение максимальной • скорости ветра и частота ее колебаний (частота порывов) могут быть выражены приближенной зависимостью [JI. 10]:
V = 1,2и 4-7,2:
,2,28 ср ' |
Где |
На |
У 1/кек |
Макс ' CD I ' '
= 0,0038а' " і
Среднее значение порывов ветра; v — число порывов в секунду. На основании этих соотношений на рис. 2-23 построе - диаграмма, показывающая, в какой степени могут
Изменяться частота и кратность порывов для заданного диапазона рабочих скоростей ветра. Представим далее внешнее возмущение со стороны ветра в виде периодического изменения скорости воздушного потока с частотой, равной частоте порывов ветра, и предположим, что при постоянном моменте Мт колебания скорости - ветра происходят по синусоидальному закону Да = v sin v^,
Где v — амплитуда колебания скорости ветра; v — частота колебания. При этом условии уравнения (2-34) запишутся в виде:
(.р3 - f а^2 - f а2р - j - а3) х2 = b3 cos v* - f 64 sin чі |
(р3 - f - ахрг + а2р - f а3) х1 = Ьг cos v* -f b2 sin vf,
Из максимумов скорости
?дё |
Найдем частное решение для координаты В силу линейнос. ти уравнений (2-43) можно аналогично (2-35) положить, что
(2-44) |
Ю> |
Хі—х.
Представив, здесь частное решение л:10 в форме
-£10 = <3 COS Vt 4 х sin
Где о и х — постоянные коэффициенты, подбираемые так, чтобы при их подстановке уравнения (2-43) обращались в нуль.
Заменяя в (2-44) л:1о его выражением (2-45), дифференцируя последовательно 3 раза и подставляя в (2-43), получим:
Р3х^ -}- v3a sin vt — tv3cosv/-|-
-f - at (p2Xtl — v2e cos v* — tv2 sin it) - f
4 a2 (px x — 5v sin 4t 4 cos vf) 4
4 «з j 4 з cos v^ 4 x sin vf) = bx cos чі 4 b2 sin vf.
Путем несложных преобразований это равенство приводится к виду:
(.Р3 + а, р2 4 агр + а3) хщ1 = [(v3 — а2у) % 4
+ (a, v2 — a3) з 4- bt] cos vf 4 [(ay — a3) % — (v3 —
(2-45) |
— aav)<7-f-62]smv/.
Для того чтобы правая часть этого выражения обра - тилась в нуль, необходимо одновременно:
(v3 — a2v) x 4 (аУ — a,) a 4 bx = 0; (a. v2 - a,) (v3 - a2v) a 4 = 0.
Отсюда находим значения постоянных коэффициентов
_ (ч3 — a2i) Ь2 — (а, у2 — а,) 6, _ ° (V3 ~ a24f + (a, v« - а,)2 '
(у* —а, у)»1 + (а1у«—
Х (>3 — е2у)2 + (е, у2 — а3)2 '
Чтобы получить выражение для частного решения в раскрытом виде, надо пол ученные значения коэффициентов а и т подставить в уравнение (2-45).
Если до начала действия возмущающей силы система находилась в установившемся движении, то можно принять прежние начальные условия [*і(0)=0, *г(0) = — О, рхг(0)=0]. Исходная система уравнений (2-20) § 2-4 дает:
/tt,(0) = 0, р'хЛ 0) = ^-.
Подставляя эти значения производных в уравнения (2-37) и замечая, что
*1С(0) = о, Pxl0 (0) = vt, /Ле10(0) = - оЛ -
Получим следующую систему уравнений для определения произвольных постоянных:
С + Л = -3;
— СЯ — Aw - f - Bq = — vt;
СЯ2 -f- А К - ?2) — 2Bqw=^- - f - v2a.
Определив произвольные постоянные, можно вычислить общий интеграл движения, который в данном случае имеет вид:
Хг = Се~и-f e~wt (A cos qt--B sin qt) - f
+ (a cos v/ - f x sin v^). (2-47)
Первые два члена характеризуют свободные колебания системы с затуханием, а третий член ■— вынужденные колебания от периодически меняющейся скорости ветра. Амплитуда свободных колеба«ий зависит только от начального состояния движения. В первый момент
Времени после нарушения установившегося движения амплитуда колебания системы будет определяться разностью частот обоих колебаний, но после того как свободные колебания исчезнут, она будет определяться величиной амплитуды колебания вынуждающей силы.
Напишем выражение для частного решения уравнений (2-43) по координате х2:
Х20 = Sj COS v/ - f - tj sin v/.
Применяя тот же метод, найдем значения коэффициентов и
__ (у3 —а2у)64 —(д, у2 —д3)63 .
• 1 (v3_a24)2+(aiv2_a3)2 .
_ (v3 —а2у)63 + (д1у2 —д3)64 1 (У3 — Д2У)2 + (д, У2 — Д3)2 '
При подстановке в исходную систему (2-20) § 2-4 начальных условий, принятых для координаты xt, получим /72х,(0) = 0, и для определения произвольных постоянных будем иметь следующие уравнения:
С,+ 4 =-о,;
— С^, — AjW - f Btg ~ — vx,;
С^2 + А (ш2 — д2) — 2B, gw = v23l.
Общий интеграл движения по координате х2 может быть запиЬан в следующем виде:
Х2 = C, e~xt - f e~wt (At cos qt + Bt sin qt) - f
-f (^cosv^-f-^sinv*). (2-48)
Для наглядного представления о характере возмущенного движения системы от периодически меняющейся скорости ветра рассмотрим несколько характерных примеров.
Пусть скорость ветра изменяется по синусоидальному закону с амплитудой, равной v= ~ v0 и периодом
Г^Ю сек (соответствует частоте v = 0,628). Характер движения системы после начала действия возмущающей силы показан на рис. 2-24. Как следует из этой диа-
9—2412 129
граммы, несмотря на некоторое запаздывание действия регулятора по отношению к изменению скорости ветра, существенного увеличения движущего момента не наблюдается. Максимальное значение движущего момента не превышает 1,35 (Мв)0.
Учитывая, что подобного рода системы содержат ряд существенно нелинейных функций, при выборе того или
Иного значения внешнего возмущения надо строго следить за тем, как при этом могут измениться коэффициенты уравнений упрощенной линейной системы (2-20) § 2-4.
Нельзя задавать такие внешние возмущения, которые вызывают изменение этих коэффициентов более чем на 25—30%. В тех же случаях, когда требуется получить характеристику точности регулирования при больших внешних возмущениях, вызывающих значительные изменения коэффициентов, следует применять метод численного интегрирования исходной нелинейной системы (2-14) § 2-4.
Увеличение скорости ветра после начала регулирования связано с переводом лопастей на другие углы
Установки, а следовательно, и с изменением параметров Ул, Угр, k0 и др. В результате изменяется частота собственных колебаний системы. Чтобы выяснить, возможны ли резонансные колебания, надо проследить изменение частоты собственных колебаний системы и частоты порывов ветра в функции скорости ветра.
Случай, показанный на рис. 2-25, предполагает увеличение частоты колебаний скорости ветра, равной теперь v=2,l.
0.3 0.2 Ф О -01 -0.2 -Q3 - Щ |
1.2 |
0.5 |
■ О |
-12 |
-05 |
Характерным здесь является то, что наблюдается увеличение запаздывания действия регулятора, вслед-
М ij
-007 |
Рис. 2-£5. Движение системы при изменении скорости ветра по закону [а=-"2_ sin 2,1/. |
Ствие чего движущий момент ветродвигателя растет быстрее, чем это имело место в предыдущем примере. В этом случае движущий момент достигает своего максимального значения 2,5 (-Мв)о в первую четверть секунды, после чего, несмотря на дальнейшее нарастание возмущающей силы, движущий момент быстро уменьшается. Установившиеся периодические колебания движущего момента в данном случае происходят с амплитудой, равной 0,8
9* |
Как видно из приведенных по ветродвигателю ID-18 примеров, для системы прямого аэродинамического ре-
131
-гулирования динамическая неравномерность скорости вращения будет составлять не более ±'5%., а максимально возможное увеличение движущего момента не будет превосходить 2,5 {Мв)й.