
книги из ГПНТБ / Хитрик, В. Э. Методы динамической оптимизации механизмов машин-автоматов
.pdfПреобразуем выражение (11.95) интегрированием по частям. Обозначим
|
и |
|
Г |
\Р п- т 2 - ^ |
|
d v = |
dv. |
|
||||
|
|
|
J |
\ |
О’ + с, |
|
|
|
|
|
||
|
|
|
с, |
|
|
|
|
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
Тогда в соответствии с |
изложенным |
находим |
|
|
||||||||
|
|
|
|
|
|
|
|
ct |
|
|
|
|
- 0 , 5 |
= у |
(p* + PdP\1:2 |
|
|
|
|
|
|
||||
\ |
|
Cj J |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|||
откуда |
|
|
с, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
|
|
И |
Pi-\- Р\Уг |
1/2 |
|
|
|
|
||||
|
|
|
|
|
dy = 0,5. |
|
|
|||||
|
|
|
|
|
Ху |
С\ |
|
|
|
|
|
|
Кроме того, |
граничное |
условие _у(0) = О приводит к coot- |
||||||||||
ношению |
|
|
с . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
|
P^-t Pi.У2 |
1/2 |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
Ху + С1 |
fify = 0,5. |
|
|
||||
|
|
|
о |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
Обозначив 7 = |
— С4 X и исключив множитель |
получим соот |
||||||||||
ношение |
|
|
|
|
|
|
|
|
|
|
|
|
|
I у ( - * |
|
|
|
= j ( |
i " - ) 1 V |
(Н .9 6 ) |
|||||
|
О |
|
|
|
|
|
|
О |
|
|
|
|
Сделаем замену переменных y = |
j z |
и обозначим |
v-— {P2lP\)1/2 |
|||||||||
у = рь2Д 2. |
Тогда |
уравнение (Н.96) |
можно |
привести |
к виду |
|||||||
|
|
1 |
|
|
|
19 |
|
|
1 |
22 |
12 |
|
|
[Л Г_ / vТ Z2 |
\ 1,2 |
dz |
р j ч+ |
(11.97) |
|||||||
|
/ |
'>SО' |
1 - |
2 |
|
|
1 |
|
|
|||
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
Параметр р, отражает соотношение допустимых величин сил инерции и динамической мощности и является известной вели чиной. Неизвестная величина v может быть определена из (П.97). Далее определяем
>= |
c, = - x i . |
Таким образом, найдены все неизвестные и задача решена.
Вкачестве примера рассмотрим типичный случай: р1= р 2= 1, р,= 1, соответствующий равноопасности и равнозначности мини
мизируемых динамических факторов. Решение (11.97) на ЭЦВМ «Минск-22» дает значение для v: v = 0,5127, причем
62
Я^-)1!* =1,9973.
О
После ряда последовательных вычислений получим
Т = j max= |
1,3966, \ = — 43,4657, С, = |
60,7024. |
||
Найденный закон движения имеет вид |
|
|||
|
|
1,3966 |
|
|
.V = |
0,5 + |
1+32 |
d% |
|
V - 60,7024 - 43,46.573 |
||||
|
|
Таблица 6
Инварианты подобия для закона комплексной динамической оптимизации
5 |
X |
£ |
; |
|
0 |
0 |
7,791 |
0 |
|
1 |
—7,791 |
1 |
||
|
||||
0,05 |
0,0033 |
7,641 |
0,000082 |
|
0,9967 |
-7,641 |
0,999918 |
||
|
||||
0,10 |
0,0099 |
7,470 |
0,00057 |
|
0,9901 |
—7,470 |
0,99943 |
||
|
||||
0,15 |
0,0169 |
7,279 |
0,00144 |
|
0,9831 |
-7,279 |
0,99856 |
||
|
||||
0,20 |
0,0231 |
7,072 |
0,00252 |
|
0,9769 |
-7,072 |
0,99748 |
||
|
||||
0.25 |
0,0306 |
6,849 |
0,00417 |
|
0,9694 |
-6,849 |
0,99583 |
||
|
||||
0,30 |
0,0383 |
6,613 |
0,00629 |
|
0,9617 |
-6,613 |
0,99371 |
||
|
||||
0,35 |
0,0461 |
6,367 |
0,00889 |
|
0,9539 |
—6,367 |
0,99111 |
||
|
||||
0,40 |
0,0542 |
6,111 |
0,01189 |
|
0,9458 |
-6,111 |
0,98811 |
||
|
||||
0,45 |
0,0622 |
5,849 |
0,01529 |
|
0,9378 |
-5,849 |
0,97871 |
||
|
||||
0,50 |
0,0712 |
5,584 |
0,01961 |
|
0,9288 |
—5,584 |
0,90440 |
||
|
||||
0,55 |
0,0805 |
5,315 |
0,02488 |
|
0,9195 |
-5,315 |
0,97572 |
||
|
||||
0,60 |
0,0901 |
5,046 |
0,03003 |
|
0,9099 |
-5,046 |
0,96997 |
||
|
||||
0,65 |
0,1002 |
4,776 |
0,03628 |
|
0,8998 |
-4,776 |
0,86372 |
||
|
||||
0,70 |
0,1109 |
4,508 |
0,04303 |
|
0,8891 |
-4,508 |
0,95697 |
||
|
5
0,70
0,75
0,80
0,85
0,90
0,95
1,00
1,05
1,10
1,15
1,20
1,25
1,30
1,35
1,3965
| |
-V |
|
|
|
0,1109 |
4,5077 |
0,04303 |
|
0,8891 |
-4,5077 |
0,95697 |
|
0,1226 |
4,241 |
0,05128 |
|
0,8774 |
—4,241 |
0,94872 |
|
0,1342 |
3,973 |
0,06028 |
|
0,8658 |
-3,976 |
0,93972 |
|
0,1477 |
3,714 |
0,07138 |
|
0,8523 |
—3,714 |
0,92862 |
|
0,1611 |
3,453 |
0,08308 |
|
0,8389 |
—3,453 |
0,91692 |
|
0,1767 |
3,194 |
0,09748 |
|
0,8233 |
-3,194 |
0,90252 |
|
0,1927 |
2,936 |
0,11311 |
|
0,8073 |
-2,936 |
0,98692 |
|
0,2108 |
2,6767 |
0,13168 |
|
0,7892 |
—2,677 |
0,86832 |
|
0,2302 |
2,415 |
0,15228: |
|
0,7698 |
—2,415 |
0,84772 |
|
0,2521 |
2,148 |
0,17708 |
|
0,7479 |
-2,148 |
0,82292 |
|
0,2772 |
1,871 |
0,20668 |
|
0,7228 |
—1,871 |
0,79332 |
|
0,3063 |
1,577 |
0,24248 |
|
0,6937 |
—1,577 |
0,75752 |
|
0.3418 |
1,249 |
0,26768 |
|
0,6582 |
-1,249 |
0,71232 |
|
0,3896 |
0,8468 |
0,39968 |
|
0,6104 |
—0,847 |
0,60032 |
|
0,5 |
0 |
0,5 |
|
0,5 |
0 |
0,5 |
63
5
Рис. 8. Кривые закона движения, синтези рованного по комплексному критерию оптимальности, учитывающему уровень динамической мощности и ускорений ме ханизма.
64
где Ь=у — инвариант скорости. В табл. 6 приведены значения инвариантов ускорения I, скорости б и пути %полученного зако на движения. Соответствующие графики приведены на рис. 8.
Если по условиям работы «мягкие» удары недопустимы, то по лученный закон движения следует откорректировать в начале
иконце интервала.
§7. Оптимизация механизмов по комплексному энергетическому критерию
Вряде случаев представляет интерес улучшение энергетиче ских характеристик проектируемого механизма как в динамиче ском отношении, так и в отношении затраченной работы дви
жущих сил. В этом случае становится целесообразным требова ние минимизации некоторого комплексного критерия, характе ризующего сумму затраченной работы и нормы работ сил инерции системы за период. При этом обобщенный момент тех нологических сопротивлений Мс полагается известной функцией
положения механизма. Улучшение энергетических характери стик механизмов может быть достигнуто рациональным выбо ром передаточной функции системы, которая, как обычно, долж на обеспечивать заданный ход ведомого звена на рассматривае мом интервале и удовлетворять условиям непрерывности и без ударного движения.
При решении этой задачи предполагаем, что ведущее звено имеет постоянную угловую скорость. Такая постановка пред ставляет интерес для задачи оптимизации силовых механизмов, у которых с ведущим звеном связаны значительные маховые массы, так что его скорость близка к постоянной.
Критерий оптимальности можно записать в виде
|
R = |
j М с (?) П' (9) d f + ||л г I- |
|
(11.98) |
|||
|
|
6 |
|
|
|
|
|
Здесь первое |
слагаемое — работа |
за |
период |
сил |
полезного'со |
||
противления |
(или, |
что то же, работа |
движущих |
сил), |
второе |
||
слагаемое — норма |
работы сил инерции за |
период. Оба |
этих |
||||
слагаемых (и величина Л2) зависят |
от выбора |
передаточной |
|||||
функции П'(ф). Если каким-либо |
образом |
найдена функция |
|||||
1Г(ф), минимизирующая величину |
|
то для механизмов с ма |
лыми инерционными нагрузками эта функция будет минимизи ровать потребную работу (и мощность движущих сил), а для быстроходных механизмов с малыми технологическими нагруз* ками — норму работы сил инерции за период.
К числу дополнительных соотношений, которым должна удо влетворять искомая передаточная функция П', относится требо вание, чтобы механизм имел заданные углы рабочего и холо стого хода (фр и фх), а также заданное значение хода ведомого звена (По).
5 Зак. 117 |
65 |
Обозначив передаточные функции на рабочем и холостом ходу Пр и Пх, запишем указанные выше условия:
П р ( 0 ) = П р (<рр) = 0 , |
| П рс Ь = П 0, |
|
о |
|
2- |
n x(tpp) = Пх(2^) = |
0, |'Пхй?Ф=: — П0. |
Обычно момент сопротивления |
на холостом ходу М с. х значи |
тельно ниже, чем на рабочем ходу Мс. р. Учитывая это, соот
ношение (11.98) можно записать |
в виде |
|
|
|
||||
|
?р |
|
|
|
|
2j |
|
|
л , - |
г [Мс. Рп ; + |
/со2(Прп;)2]d* + |
j /о)2(П;п ;)2^ , |
|||||
|
о |
|
|
|
|
?р |
|
|
где / — момент инерции |
(для вращательного движения ведомо |
|||||||
го звена) |
или масса |
(для поступательного движения ведомого |
||||||
звена); со — угловая скорость ведущего |
звена. |
|
|
|||||
Рассмотрим два случая. |
|
|
|
М с. Р = const. |
||||
1) П о с т о я н н ы й |
мо ме нт |
с о п р о т и в л е н и я |
||||||
Обозначив |
|
м |
, |
Пр=_у, ПР= У , |
придем |
к |
изопериме- |
|
М 0 -- |
||||||||
трической |
вариационной задаче |
о минимизации |
функционала |
|||||
|
|
|
|
р |
|
|
|
(11.99) |
|
|
R = |
\ (МоУ + У2У’2)й? |
|
||||
при наличии |
условия |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
f V f = |
n , |
|
|
(11.100) |
и граничных |
условий |
|
и |
v |
|
|
|
|
У (0)=^(<рр) = 0. |
|
|
(11.101) |
|||||
|
|
|
|
|
Эта задача допускает решение точными методами. Уравне ние Эйлера имеет для функционала (11.99) вид
y y - + y'2y - - ^ f ^ = 0, |
(П.102) |
где /- — подлежащий определению множитель Лагранжа. По лученное уравнение совпадает с уравнением задачи о мини муме динамической мощности (§ 6, гл. II), что является след
ствием предположения о постоянстве момента сопротивления. Решая его, получим
<?— С2 ± 2(^ :2Cl) V a y + С„ |
(11.103) |
где а = Ж0-(-Х. Постоянные интегрирования Q и С2 опреде
ляются из граничных условий (11.101) Сх = ^ - а у 9а?‘ , С . =
66
= 0,5 срр. Теперь легко найти наибольшее значение передаточ
ной фуНКЦИИ На рабочем |
ХОДУ Пргаах=.Утах = — У Oa'fp. |
Для |
|||||
окончательного |
решения задачи нужно определить параметр а. |
||||||
Он найдется из |
условия |
|
|
по |
|
||
(11.100) а = 3,350—=-. |
|
||||||
Для |
случая |
|
®р = 210°, |
П0= |
1 |
|
|
|
ф= |
1,8326 У; 5°>8у ~~3д1,7468 V 50,8у + 10,8. |
|
||||
2) |
П е р е м е н н ы й |
м о м е н т с о п р о т и в л е н и я . |
По-преж |
||||
нему считаем, |
что момент сопротивления |
на холостом |
ходу |
||||
пренебрежимо |
|
мал. Тогда для |
построения оптимальной пере |
||||
даточной функции на холостом |
ходу Пх |
можно руководство |
|||||
ваться изложенными выше результатами. |
Задача об установле |
нии оптимальной передаточной функции на рабочем ходу Пр
заключается |
в следующем. Учитывая, что при 0 < ?<С?Р функ |
||||
ция Пр должна |
быть знаконоложительной, и |
обозначая |
|||
|
т = т (?) = |
Л4 (ф! |
|
||
|
— |
— = km0(®), |
|
||
поставим задачу |
о минимизации |
искомой функцией функцио |
|||
нала R: |
|
*р |
|
|
|
|
|
|
|
|
|
|
|
R = j ( m n p+ n pn p ) d ? |
(11.104) |
||
|
|
6 |
|
|
|
при условии, |
что |
|
|
|
|
|
|
Пр(0) = |
Пр (ч»р) = 0 , |
(11.105) |
|
|
|
*р |
, |
|
(11.106) |
|
|
|Прй?? = П0. |
|||
|
|
о |
|
|
|
Здесь ш(<р) |
и т0(ср) — безразмерные функции |
сопротивления; |
параметр k характеризует соотношение технологических и инер
ционных нагрузок в проектируемом механизме.
Для функционала (11.104) соответствующее уравнение Эй лера не сводится к квадратурам, поэтому .тля решения постав ленной задачи мы применим прямой мелид Ритца. Будем искать
передаточную функцию |
|
Пр = ®(?р — <p)(«?2- f by + c), |
(11.107) |
причем параметры а, Ь, с находятся из условий относительного
минимума функционала (11.104). В соответствии с (11.106) ко эффициенты а, Ь, с должны удовлетворять соотношению
?р |
?р |
|
¥р ■ |
f ==a~ i + b '\1 |
+ |
с " г — По== °- |
|
5* |
|
|
67 |
Подставив выражение (11.107) |
в (11.104), |
получим зависимость |
R — R(a, b, с), после чего в соответствии |
с правилом неопреде |
|
ленных множителей Лагранжа |
будем искать минимум функции |
Ф = Я + X/.
Нетрудно видеть, что определение коэффициентов а, Ь, с сво
дится к решению нелинейной системы уравнений
(II. 108)
Так как одно из соотношений системы (11.108) |
(/= 0) |
ли |
|
нейно относительно искомых коэффициентов, |
то |
оказывается |
|
возможным методом исключения неизвестных |
получить |
одно |
|
уравнение 4-й степени относительно какой-либо |
неизвестной |
или систему двух уравнений 2-го порядка, что удобнее с расчет
ной точки зрения. После ряда преобразований получим общую
систему уравнений |
для |
определения |
параметров оптимизации |
|||||
а, Ь |
|
|
|
|
|
|
|
|
P i а 2 Л - Р ф г + Р г Р Ь - \- P ifi - \ - р ъЬ - f Рч = |
0, |
(II. 109) |
||||||
Яid2+ ЯФ2+ |
ЯгаЬ + |
?*а + |
Я-Ф+ ?6 = |
0, |
||||
|
||||||||
где |
|
|
|
|
|
|
|
|
Pi = |
0,274701?$, |
p t = |
0,009523?*, |
|
|
|||
|
|
|
|
|
P> |
|
|
|
РP%i = 0,455238??,,238?p p tk = |
-4,144762- |
|
|
|||||
|
|
/?5= 0,800000 — , |
|
|
||||
96 |
По |
. 1 |
20 |
f |
|
|
|
|
|
|
|
|
j* m y3 (cpp — cp)d y — |
|
0
q, = - 0,725298?4p, ?2=0,009096?*,
Яг = - 0 , 64761 V p, = ~ 17,276192-^,
<75 = 2 ,0 0 0 0 0 0 -^ ,
Yp
68
После того как система (11.109) решена, величина с может быть
определена из соотношения
Система (11.109) имеет обычно несколько решений. Поэтому результатом является, как правило, семейство передаточных функций, сообщающих экстремум функционалу (11.104) и удо влетворяющих условиям (II.105) и (11.106). Среди этих функций возможны и такие, которые меняют знак на отрезке [0, <рр]. Эти функции должны быть исключены из рассмотрения, так как их невозможно использовать при синтезе реальных механизмов. Из числа оставшихся (допустимых) функций непосредственным вы числением функционала (11.104) выбирается функция, сооб щающая этому функционалу наименьший минимум.
Пусть заданы угол рабочего хода механизма <рр = 210° и наи большее значение функции положения ведомого звена (ход или угол качаний) П0=1. Момент сопротивления на рабочем ходу
является квадратичной параболой, так что
Отношение параметра технологических сопротивлений к
энергетическому |
уровню |
ведомого звена |
Л»2 обозначим |
k = |
|||||
= - ^ - , |
тогда |
rti — k — |
^? — |
. Система |
уравнений (11.109) в |
||||
данном |
случае имеет вид |
|
|
|
|
|
|||
49,573125а2 + |
(22,317441 6 — 1,130845) а + (0,12892S62 + |
|
|||||||
|
+ 0,0595526 — С,002836/г — 0,001131) = |
0, |
|
||||||
130,888815а2 + |
(31,8866886 - f 4,638990) а - (0,12219262 + |
|
|||||||
|
+ 0,1488806 - |
0,002836* — 0,001131) = |
0. |
|
|||||
Решим задачу |
для |
случаев kt — 0,1, |
k2 = \. |
1) *! = |
0,1. |
||||
После |
вычислений получим |
|
|
|
|
||||
|
П! = |
<р (<Рр — ?) (0,0107<р2 ■+- 0,0259? + 0,0302), |
|
П2= ? (?Р — ?)(0,0259? + 0,0745),
Пз = ? (?Р — ?) ( - 0,0107?г + 0,0259? + 0,1165),
П4 = ? (?р — ?) (0,0235?2 — 0,0304? - f 0,0854),
Пб = ? ( ? р — ?) (0,0131?2 — 0,0304? + 0,1269).
Этим функциям соответствуют следующие значения функцио налов (11.104): /?! = 0,643, R 2 = 0,354, R3 = 6,46, Я4= 2,4, # 5 =
— 18,8. Очевидно, что наилучшей в данной постановке функ цией П' (?) является Пг:
П0пт = Пг = ? (3,665 — ?)(0,0259? + 0,0745).
. Совершенно аналогично решая задачу для k —\, получим
ПоПХ= ® (3,665 — <р) (0,28<р2 + 0,010).
Графики этих функций приведены на рис.-9.
Рис. 9. График закона движения, син тезированного из условий минимума комплексного энергетического критерия.
§ 8. Приближенный метод построения оптимальных безударных законов движения
I
Как уже указывалось, применение точных методов, связан ных с интегрированием дифференциального уравнения Эйлера для данной вариационной задачи, ограничивается тем, что по условиям работы механизма искомый закон движения должен удовлетворять дополнительным граничным условиям. Поэтому полное число граничных условий превышает число постоянных интегрирования уравнения Эйлера.
Например, функция у(х), сообщающая минимум величине
среднеинтегральных ускорений ведомого звена |^ |
1
||^||еее/? = j y '2d x
о
и не имеющая „жестких" ударов, должна |
удовлетворять |
гра |
|
ничным условиям следующего типа: |
|
|
|
У(0) = я, У (\) = $. |
|
(11.110) |
|
Однако если по условиям работы механизм не должен |
|||
иметь „мягких" ударов и должны быть |
удовлетворены |
соот |
|
ношения |
|
|
|
У ( 0 ) = т. У ( 1) = |
8, |
(11.111) |
70
то использование точных методов невозможно, так как число граничных условий (4) превышает порядок дифференциального уравнения Эйлера п = 2.
Наряду с основным функционалом |
|
||
|
1 |
|
|
R = |
^ F (х, |
у, y')d.x |
(11.112) |
|
О |
|
|
рассмотрим вспомогательный функционал |
|
||
1 |
|
|
|
JF(x, |
У, y') + |
j r f ( y " ) dx, |
(11.113) |
где / — функция, удовлетворяющая всем необходимым требова ниям непрерывности и гладкости; р — вещественное число.
Функция, сообщающая минимум функционалу (11.113), мо жет быть найдена из дифференциального уравнения Эйлера— Пуассона, которое для данной задачи имеет вид
6 F ____d |
I dF |
_1_ _сР_ \дНГУ\ |
= |
0, |
||
ду dx |
( ду’ |
dx2 |
ду" |
|||
|
|
т. е. является уравнением 4-го порядка, общий интеграл кото рого содержит 4 постоянные интегрирования. Эти постоянные интегрирования могут быть определены из граничных условий (11.110), (11.111). Таким образом определяется функция у(х,р),
удовлетворяющая 4 граничным условиям и сообщающая мини мум функционалу (11.113). Можно показать, что при достаточно больших р по модулю функция у(р,х), удовлетворяющая гра
ничным условиям |
(11.110)— (11.111), сообщает функционалу |
R Pзначение, сколь |
угодно близкое к минимальному, сообщаемо |
му ему функцией у(х), которая является решением вариацион ной задачи (11.110) — (11.112).
На практике проще всего выбирать функцию / ( / ' ) в виде /( У ') = У '2. В этом случае уравнение Эйлера—Пуассона имеет вид
Принимая в качестве критерия оптимальности величину среднеинтегральных ускорений ведомого звена, находим, что искомая функция у(х) должна сообщать минимум функционалу
1
R= J y '2dx
О
при выполнении однородных граничных условий
у (0) = 0-, |
(11Л14) |
71