книги из ГПНТБ / Оптимизация процессов грузовой работы
..pdfПлощадь поперечного сечения и |
массу |
груза в вагоне легко найти из |
|||
построений, приведенных на рис. 2 1 , а: |
|
|
|
||
т (Ф) = |
V c t§ (Ф — Фо) — . |
(ГѴ .65> |
|||
|
z |
|
|
s |
|
Уравнение движения на третьей фазе после преобразований |
|
||||
[/<, + % ctg (ер— фо) Я 2 ( ф ) ] ^ |
= М В— M g- |
|
|||
— P /у sin ф + É»! ctg (ф— фо) у (ф)— сг ctg (ф— фо), |
(ІѴ .6 6 ) |
||||
где |
|
|
|
|
|
а ,= |
H 2ly d |
bi |
Н Ч у |
|
|
Сі==^ |
г |
2 |
|
||
2 g |
|
|
|||
В четвертой фазе угол поворота ротора изменяется в границах — ^ ф
<+ фо (рис. 21, б). В начале этой фазы меняется знак статического
момента т ( ф ) г / ( ф ) £ , который начинает действовать в направлении, противо положном вращающему моменту. В конце фазы заканчивается разгрузка вагона. Выражение для определения массы груза т(ф) имеет вид (IV .65). Д ля этой фазы уравнение движения системы ротор— вагон следующее:
[/o + Ö! ctg ('Ф--- Фо) Я 2 (Ф )1 -2 - = Ма----М о: — Ргг sin Фі — dv
— bLctg (ф— Фо)у (ф)— Сх ctg (ф— Фо). |
(ГѴ .67) |
Значения постоянных аг, Ьх и сг прежние.
Чтобы установить законы изменения радиуса инерции R(ф) и координа ты г/(ф) от угла поворота ротора, необходимо получить уравнение траекто рии центра тяж ести поперечного сечения груза в процессе выгрузки полува гона. Анализ качественной стороны процесса свидетельствует о том, что при выгрузке центр тяж ести массы груза О' занимает последовательно положения О' — 0[ — 0'2— Оз — 0 4', перемещаясь постепенно вправо внизотносительно центра вращения ротора О (см. рис. 2 0 , 2 1 ).
Точно определить траекторию Д(ф) и зависимость у(ф) довольно слож но. Сделаем это приближенным способом. Задаваясь различными характерными положениями ротора фг (например, относя их к границам фаз), можно для каждого из них на основе несложных геометрических построений рассчитать
130
величины Л?(cp) и у(ср), а затем в границах каждой фазы построить аппрокси мирующие их зависимости. На рис. 22 построены указанные зависимости для случая, когда угол естественного откоса груза р = 35-4-40° и ср0 = 45°. Анализ показывает, что с достаточной степенью точности можно принять ку сочно-линейную зависимость R(ср) и у(ср):
В первой фазе выгрузки с ..................... |
#(ср) = г0 л |
Во второй и третьей фазах ...................................... |
R(q>) = аср-{-ß МІѴ.68) |
В четвертой фазе .......................................................... |
Л (ф )~ «іф + Рі |
Постоянные а и р можно определить, совместно решая два уравнения с двумя неизвестными, подставив для второй фазы в первое уравнение си
стемы (IV .6 8 ) начальные условия ср =<р0, ІДф) = |
і?(ф0) и конечные условия |
||||||||||||||
|
|
|
Ф = Фо + |
arctg — , |
#(< р)= Я (фо + |
a rc tg — |
J ; |
|
|
||||||
для |
третьей |
фазы |
начальные |
условия |
ф = |
|
|
н |
R (ф) = |
||||||
ф0 - j- a r c t g — , |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
В |
|
|
= R ^Фо + arctg |
|
и конечные |
условия ф = Д - , |
R (ф) = R |
. |
Ана |
|||||||||
логично |
определяем |
величины |
ссх |
и |
ßx |
во |
втором |
уравнении |
систе |
||||||
мы |
(IV .6 8 ), |
для |
которого |
начальные |
условия |
ф = |
, R (ф) = |
R |
|
||||||
конечные |
ф = - у + |
фо, R (ф) = |
R f ф0 + |
~ |
j . |
Д л я |
у ( ф) |
получим |
следую- |
||||||
а) |
|
|
|
|
|
|
|
б) |
|
|
|
|
|
|
|
Рис. 22. Аппроксимирующие зависимости:
а —ЯСф); б — і/(ф)
5* |
131 |
щие зависимости: в первой фазе у — r0sin (я + ф0); во второй, третьей и четвертой фазах принимаем одну аппроксимирующую зависимость
|
у(ф) = Д с р + 5 і . |
(ІѴ |
.6 9 ) |
Величины А г и Вх определяем по начальным ср = |
ср0, у(ср) = у(ср0) |
и ко |
|
нечным условиям ф = |
+ фоі г/(ф) = У ^фо + |
. |
|
П осле подстановки |
соотношений (IV .6 8 ) и (ІѴ |
.69) в (ІѴ .64), (ІѴ |
.6 5 ) |
и (ІѴ .67) получим развернутые уравнения движения системы ротор— вагон во второй, третьей и четвертой фазах работы вагоноопрокидывателя. Перей дем к методу и алгоритму поиска закона оптимального управления роторным вагонорпрокидывателем. Интегральный функционал, выражающий эксплуа тационные расходы за один полуцикл работы вагоноопрокидывателя, имеет следующий вид
т |
|
/ = I (a0+ b 0y \u \)d t, |
(IV .70) |
о |
|
где а0, Ь0— постоянные коэффициенты — удельные расходы на заработную плату и электроэнергию.
П ервая составляющ ая этого выражения — расходы на заработную плату, вторая — на электроэнергию и смазочные материалы. Если состояние управляемой системы описывают уравнения (IV .59), (ІѴ .64), (ІѴ .65) и
(ІѴ .67), то задача заклю чается в том, |
чтобы найти в каждой фазе такой за |
||
кон оптимального |
управления u(t), при котором интегральный функционал |
||
(IV .70) приводится |
к минимуму. При |
этом на управление u(t) налагаю тся |
|
ограничения: |и |^ |
етах. Заметим, |
что могут быть такж е наложены ограни |
|
чения на угловую скорость ротора |
ф ^ |
со0. Величина и на отдельных фазах |
|
разгрузки выражается следующим образом:
первая фаза и1 = Мв— Мс— dm; вторая фаза и2 = М в— М0;
третья и четвертая фазы и3= щ = М в— М
Поэтому в общем случае иг Ф и2 Ф (и 3, и4) и в каждой фазе система ограничений, налагаемых на и, имеет вид
I I ^ ех ш ах, I « 2 1 ssC е2 ш ах и |и3|^ е3 ш ах.
Д ля решения задачи целесообразно воспользоваться численными метода ми— алгоритмом пошаговой оптимизации, например динамическим програм мированием. Д л я этого представим уравнения (V I .59), (V I.64), (V I.65) и
132
(IV .6 7 ) в конечных разностях, перепишем их в фазовых координатах ф1; ф2,
а функционал (IV .70) |
запишем |
в |
виде интегральной суммы: |
|
|
|||||
|
|
|
|
/ = |
|
2 |
( О о + Ь о Ф я |“ ь | ) Д ( !ѵ |
|
(IV .71) |
|
|
|
|
|
|
к= |
1 |
|
|
|
|
где k = |
1, 2, ..., N — |
число |
ш агов итерационного процесса. Д ля |
первой |
||||||
фазы из системы (IV .59) следует: |
|
|
|
|||||||
|
|
|
|
|
|
|
Афг --= фзДг1; |
I |
(IV .72) |
|
|
|
|
|
|
|
|
|
|
||
|
|
Дф2 ф2 = Дфху [% — Prx sin фг — mgr0sin (ф + я)]. | |
|
|||||||
Д л я |
второй |
фазы, |
учитывая |
уравнение (IV .64), |
|
|
||||
|
|
|
|
|
Д ф і = |
ф 2 А ^ ; |
|
|
||
Дф2 ф2= Дф |
ца+ с tg (Фі— Фо) — Р г г sin y x+ g ( a — b tg (Фі— Ф))(Лх Ф + Д Д |
(IV . 73) |
||||||||
|
|
/o+ |
Ia'Pi + ß)2 (а — b tg (q>— Фо)) |
|
|
|||||
|
|
|
|
|
|
|||||
для |
третьей |
фазы из уравнения (IV .65) |
|
|
||||||
Дф ф=Дф ( из~ |
|
|
Дфх = |
ф2Д^; |
|
|
||||
Pri sin фі~^bl ctg f o - V o W 1 фі+g0~gictg (Фі—Фоі |
(ІѴ .74) |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
\ |
|
^о+ |
аі сі£ (Ф і— Фо) (афі + ß)3 |
|
|
|||
для четвертой фазы из уравнения (IV .67)
|
|
|
|
Дфі = ф2At |
|
Дф |
ф = Д ф |
( “ i — P ri sin Фі— 1ь і ctg (Фі— Фо) (Лі Фі+Ді)— Ді c tg (Фі—Фо)' |
(IV .75) |
||
|
|||||
2 |
2 |
\ |
+ |
9 o )(a i cp1 -{-ßi)2 |
|
|
Уравнения (ІѴ |
.72)— (ІѴ |
.75) выражают состояние системы ротор— вагон. |
||
Вид их свидетельствует о том, что оптимальная фазовая траектория движе ния изображающей точки, моделирующей работу вагоноопрокидывателя в рабочем полуцикле, представляет собой кусочногладкую кривую, состоя щую из четырех участков. Каждый участок соответствует определенной фазе работы вагоноопрокидывателя.
Оптимальное управление u(t) и уравнение фазовой траектории целесооб разно искать одновременно от начального и конечного положений изобра жающей точки, навстречу друг другу до тех пор, пока на последнем этапе
133
|
скорости |
не сравняю тся, а |
сумма |
||||||
|
отрезков |
пройденных путей |
не бу |
||||||
|
дет равна заданному углу поворо |
||||||||
|
та ротора ерт- Процедура эта от |
||||||||
|
начальной точки а с координатами |
||||||||
|
t — О, ф10 |
= |
Фт, ф2 г = |
0 и конечной |
|||||
|
точки |
О |
с |
координатами |
t = |
Т, |
|||
|
срІГ |
= |
0 , |
ф2 г = 0 |
изображена |
на |
|||
|
рис. |
23. |
|
|
|
|
|
|
|
|
|
Оптимальное управление ищем |
|||||||
|
в направлении от точки а к началу |
||||||||
|
координат |
|
при |
положительном |
|||||
Рнс. 23. Схема многоэтапного поиска опти |
значении |
Uj, а в направлении от |
|||||||
мального закона управления вагоноопро- |
точки |
О — |
при отрицательном |
|
|||||
кипывателем |
В данном |
|
случае |
/ — количество |
|||||
|
|
||||||||
отсчитываемых от точки а (J = 1 , |
шагов |
многоэтапного |
процесса, |
||||||
2 , ..., |
|
/У), |
і — количество |
ш агов, |
|||||
отсчитываемых от точки О (і — 1 , 2, .... М). И сходя |
из описываемой проце |
||||||||
дуры поиска, вытекает, что М — N. Поэтому индексы / и і обозначают лишь |
|||||||||
направление поиска. Чтобы повысить |
точность расчетов и улучшить сходи |
||||||||
мость итерационного процесса на завершающих этапах, может оказаться целесообразным уменьшить величину шага Дфд Рассмотрим процедуру поис ка на первом этапе.
|
1. |
Задаемся начальным |
приращением |
угла поворота |
ротора |
Дфід |
|||||
и величиной « Л в границах, определяемых условием (|«д| = |«і|) ^ е ^ іа х . |
|||||||||||
Подставляем ДфіЯ |
и Дф2Л в уравнения (ІѴ .72) |
и решаем их относительно |
|||||||||
Дtj\ |
и Дфзд. В данном случае Д/д— приращение времени на первом ш аге |
||||||||||
итерационного процесса, |
Дф2д — средняя скорость поворота ротора в ин |
||||||||||
тервале (0, Дфід) на участке разгона. П ервая цифра в индексе величин |
М |
||||||||||
и и — номер шага при переборе значений и на каждом этапе; вторая цифра — |
|||||||||||
номер этапа. Аналогичные значения |
имеют вторая и третья цифры в индек |
||||||||||
сах величин Дфхд, Дф2д и ф2д . П осле подстановки Дфід |
и Ыд |
уравнения |
|||||||||
(ІѴ .72) |
переписываем следующим |
образом: |
|
|
|
|
|||||
|
|
|
|
|
ДфіЛ— Фгд А^д; |
|
|
|
|
||
|
|
ДфІд |
|
Кд — PrxSin Дфід — m qr0 sin (Дфід + я) |
(IV . 76) |
||||||
|
|
Дфхд -!■----------------------------------=----------------------• |
|
|
|||||||
|
Функционал записывается так |
|
|
|
|
|
|
||||
|
|
|
І 1}-= |
(а0-j- Ь0А(р Д I Ujy I) Д(д, |
|
(IV .77) |
|||||
|
|
|
|
|
|
|
2 |
|
|
||
причем |
ид > 0. |
Варьируя |
« д |
в |
границах, |
определяемых |
условием |
||||
( j |
|= I цд I ) ^ 6 j m ax |
и воспользовавш ись системой |
(IV .7 6 ), |
находим |
|||||||
134
такие |
и*,, а |
такж е |
соответствующ ие им |
значения |
А**,, |
Асрі/і = фг/ь |
ко |
||
торые минимизируют функционал (IV .77). |
и ДАъ |
|
|
|
|||||
2. |
Д л я определения параметров |
Acp1 (1 |
и ф2 і1 воспользуемся |
||||||
итогами расчетов предыдущего полуэтапа, принимая во |
внимание, |
что |
|||||||
ФаЛ ~ |
Фгп. и уравнениями |
(ІѴ .75), которые переписываем следующим об |
|||||||
разом: |
|
|
|
|
|
|
|
|
|
|
|
|
|
" / 1 |
— Р 'г sin Д ф іа — 6i clg (ДфШ — фо) X |
|
|
||
|
Дф2 |
= Дф |
|
X (ѵ 4і А ф щ + Д і) — CiCtg(A fi a — Фо) |
|
|
|||
|
m |
|
і а |
Л> + Яі ctg (Д ф ш — фо) (й! Дфк'і + ß])2 |
(ІѴ |
.78) |
|||
^Фш = Фіі& п -
Ввы раж ении функционала (IV .77) индекс j заменяем на индекс і;
ограничение, |
налагаемое |
на |
и, |
записываем так (| ип |= |
|и4 \) ^ |
е3 m ax, |
при этом ип ^.О. Варьируя |
в |
заданных границах величиной ип и вос |
||||
пользовавш ись |
системой |
(ІѴ .78), |
находим такие и*ь Дф*,ь |
А/*ь |
которые |
|
приводят функционал (IV .77) к минимуму. Определяем величину этого функционала.
Рассчитываем сумму Д фн і +Д ф і /і и проверяем условия Аф?п +
+ АфТл < ф Г И ф2П <со0.
3. Затем последовательно переходим ко второму, третьему.........k эта пам расчетов. Ф азовые координаты для любого этапа k расчетов {к = 1, 2, ..,
.... N) |
определяем при |
помощи следующих |
рекуррентных соотношений |
|
Ф 2 ih = Ф г л й - і ) + А Ф 2 ^ ; |
Ф ш = Ф и ( * — и + |
А Ф ш ; |
к = Е ( А ^ д + A ^ / J ; |
|
Філ= 2 ( а Фіій+ а Фш )- |
|
|
|
|
k |
|
|
|
|
Если величина Дф = |
фГ — ф1д достаточно мала |
и условия ф2;-й < 1 со0 |
||
и \и \ |
emax соблю даю тся, то следует перейти к последнему этапу расчетов, |
|||
который, по-видимому, соответствует третьей или четвертой фазе разгрузки вагона. Е сли ж е при проверке неравенство ф1Л < фГ не соблюдается, то надо уменьшить величину шага АфіЛуѵ _ и на предпоследнем этапе и повторить процедуру расчета.
В итоге итерационного процесса будет определен закон u{t) и по извест ным значениям ф1;-й, ф2;-к построена оптимальная фазовая траектория дви жения изображающей точки, найдены фазовые координаты и моменты точек переключения ускорений при переходе из одной фазы разгрузки в другую , а такж е точка переключения с режима разгона на режим торможения ротора. Если в процессе поиска текущие значения ф22-й (или ф2;й) достигнут величи ны со0 ранее, чем пересекутся между собой интегральные кривые разгона и торможения (по-видимому, этот случай наиболее типичен для работы ваго
135
ноопрокидывателя), то следует определить координаты точек переключения фазовых траекторий с прямой cp2fe = с£>0 (рис. 23). Таким образом заверш ает ся синтез оптимального управления роторным вагоноопрокидывателем при реализации рабочего полуцикла. Процесс поиска оптимальных параметров управления холостой части полуцикла при возвращении ротора с порожним вагоном в исходное положение значительно упрощен, так как уравнение (ІѴ .72) при т = 0 описывает все фазы работы вагоноопрокидывателя, ме няется только знак момента Е М .
Алгоритм поиска оптимальных параметров управления можно реализо вать на Э ВМ или в схеме автоматического оптимизатора. Основная функция такого оптимизатора состоит в решении уравнений (ІѴ .7 2 )— (ІѴ .75) на от дельных этапах процесса разгрузки вагона, определении моментов переклю
чения управлений и выработке управляющ их |
воздействий для остановки |
и переключения электропривода ротора. Д ля |
этого оптимизатор должен |
получать постоянную информацию от специальных датчиков о текущ их зн а чениях угла поворота ротора фг и его скорости ср.
Существенное отличие роторного вагоноопрокидывателя от подъемно поворотного то, что у второго ось вращения расположена за габаритами по лувагона выше его кузова. На рис. 24 изображены расчетная схема роторно
подъемного вагоноопрокидывателя |
и отдельные фазы его |
работы. |
Введем |
|||||||
обозначения: |
R(ф) — по-прежнему |
радиус инерции |
массы |
груза в |
вагоне; |
|||||
у(ср) — плечо силы тяжести груза относительно оси вращения ротора: |
|
|||||||||
|
|
|
У(ф) = Я (ф) sin (ф + |
а„). |
|
|
(IV .79) |
|||
|
Заметим, |
чго |
радиус инерции R(ср) изменяется |
от максимального |
зн а |
|||||
чения Ro при ср = |
0 и ср = сро до минимального R k при ф = -у - + |
ф0. |
Д о |
|||||||
пускаем, что при изменении угла в пределах от ф = |
ф0 (конец первой и на |
|||||||||
чало |
второй |
фазы разгрузки) до ф = ф0 -j- у |
(окончание поворота ротора) |
|||||||
R(ср) |
изменяется |
линейно в зависимости |
от ф. Учитывая, |
что — ^ |
|
|
||||
^ 10 ч - 12%, а по абсолютной величине |
|
|
|
АО |
|
|||||
R 0— R h не превышает 0 ,6 м, при |
||||||||||
веденное допущение можно считать корректным, а выражение радиуса инер ции записывается так
я
( У - Ф о). |
( I V .80) |
если
Фо < Ф < Фо +
136
/
1.37
Следовательно, |
|
|
|
У(Ф )= i^ o — - (Ron Rh) ) (У — 9 o)sin((p + |
a B). |
(IV .8 1 ) |
|
В соотношениях (IV .80) и (IV .81) R h представляет собой |
радиус |
инерции |
|
массы груза в конечный момент разгрузки, когда ср — <р0 |
-)- |
. |
|
Примем еще одно допущение, вытекающее из того, что у подъемно-по воротного вагоноопрокидывателя частично статически уравновеш ивается вся система ротор— вагон, а именно в течение всех фаз работы вагоноопрокиды вателя суммарный момент тары вагона и веса вагоноопрокидывателя направ
лен в |
сторону |
поворота ротора. Полагаем такж е, что статический момент |
|
массы |
груза во |
всех фазах выгрузки |
направлен в сторону, обратную вра |
щению ротора, за исключением, быть |
может, конечного момента третьей за |
||
вершающей фазы. Однако в этот момент масса оставш егося груза в вагоне незначительна, ее статический момент, направленный в сторону вращ ения, мал, поэтому им можно пренебречь.
Перейдем к математическому описанию отдельных фаз работы подъем но-поворотного вагоноопрокидывателя. Учитывая специфику конструк ции, процесс разгрузки полувагона можно разделить на три фазы. В первой фазе при изменении угла поворота в границах 0 ^ <р ф0 груз из вагона не высыпается, положение центра тяж ести массы груза не меняется; радиус инерции такж е не изменяется. Д л я этой фазы уравнение движения системы вагоноопрокидыватель— вагон следующее:
[(mp + mB)r? -f |
= М В—М с + |
|
+ Pr1sin (ф + я 4- ßo) — mgR0sin (ф + а 0), |
(IV .82) |
|
где гх — расстояние от оси вращения до центра тяж ести суммарной массы вагоноопрокидывателя и тары вагона (радиус инерции этой массы груза);
а 0 — угол между радиусом инерции массы груза в вагоне и вертикалью , проходящей через ось вращения, в исходном положении вагоно опрокидывателя при t = 0 и ф = 0 ;
ß0 — угол между радиусом инерции суммарной массы вагоноопрокиды вателя и тары вагона и вертикалью , проходящей через ось вра щения, при t = 0 и ф = 0
В о второй фазе работы вагоноопрокидывателя при Фо ^ ф ^ arctg -щ4
4 фо происходит высыпание груза с нарастающей скоростью ; непрерывно изменяется радиус инерции и плечо веса груза относительно оси вра
138
щения |
по закону, |
определяемому выражениями (IV .80) и (IV .81). И з |
меняется |
такж е и |
масса груза в зависимости от угла поворота ср. Заметим, |
что для подъемно-поворотного вагоноопрокндывателя зависимость площади поперечного сечения и массы груза от ср такая ж е, как и для роторного. Д ля второй фазы работы вагоноопрокидывателя уравнение движения системы вагоноопрокидыватель— вагой в общем виде записывается так:
І(/Лр + тв) г\-Ь т (ср) R 2 (ср)] |
■= М в — Мс0 — |
|
|
||
— dm (ср) — т (ф) gR (ф) sin (ф + а0) + Pi\ sin (ф + |
л - f ß0). |
(IV .83) |
|||
Если ось вращения вагоноопрокидывателя опирается на подшипники |
|||||
скольж ения, то d = ріг0, где г0— радиус |
цапфы. Величина і?(ф) |
и масса |
|||
груза т(ф ) определяются соответственно |
выражениям |
(ІѴ |
|
.80) и |
(IV .63). |
н |
|
я |
\ |
вагоноопроки |
|
(при arctg — + фо ^ ф ^ фо + |
|
|
|||
дывателя характеризуется интенсивным процессом высыпания груза из ва гона (в начале и в середине фазы), который замедляется в конце фазы. Конец третьей фазы протекает в режиме торможения. Уравнение движения для третьей фазы выгрузки имеет вид (ІѴ .83), масса груза т(ф) определяется формулой (ІѴ .65). Д л я поиска закона оптимального управления фі(ф2) урав нения движения (ІѴ .82) и (ІѴ .83) переписываются в фазовых координатах, а затем интегрируются численными методами. Алгоритм поиска оптимально
го |
управления ф^фа) такой лее, как для роторного вагоноопрокидывате |
ля |
(см. рис. 23). |
