книги из ГПНТБ / Хитрик, В. Э. Методы динамической оптимизации механизмов машин-автоматов
.pdfгде Х; — неопределенные |
множители. |
|
Необходимое |
условие |
||||
минимума функции |
Ф имеет вид |
|
|
|
||||
|
# - |
= |
0, |
т = 1,2, . . . |
, |
N. |
|
|
|
vam |
|
|
|
|
|
|
|
Обозначив для |
симметрии |
Xx = a v+1, |
\ 2 = aN+2, \ 3 = a N+3> |
|||||
после преобразований |
получим |
|
|
|
||||
N |
(да |
|
1) (fe ~t~ |
1) |
|
|
|
|
дФ |
|
|
|
|
||||
~да^ |
т + k -J- 1 |
a k + a N + l + (m + 1) a N + 2 " Ь |
||||||
|
“ ,V+3 = |
0, |
m — 1, 2, |
. . . , N . |
(11.26) |
|||
|
да -у 2 |
|
|
|
|
|
||
Полученные N уравнений (11.26) совместно с соотношениями |
||||||||
(11.22) — (11.24) |
образуют |
линейную |
алгебраическую |
неодно |
||||
родную систему уравнений N + 3 порядка. |
|
|||||||
Для решения этой системы целесообразно использование стандартных вычислительных программ для ЭЦВМ. В частно сти, полученная система при N = 9 решалась на цифровой вы
числительной машине «Минск-22» для случаев: а) выстой—
перемещение—выстой |
(или полный ход шарнирного |
механиз |
||
м а ) — а = р = 0; б) разбег системы — а = 0; р = 0,05; |
0,10; |
0,15; |
||
0,20; в) торможение |
системы — р = 0; |
а = 0,05; 0,10; |
0,15; |
0,20. |
Выбор численных |
значений а и р |
для разбега и торможе |
||
ния охватывает характерные случаи в работе механизмов ме таллорежущих станков-автоматов (автоматов продольного то чения и токарно-револьверных автоматов). В табл. 2 приведе ны значения коэффициентов оптимизации аи для рассмотрен
ных вариантов. |
распространенного |
случая |
а = р = 0 искомый |
||||
Для |
широко |
||||||
закон движения имеет вид |
|
|
|
|
|||
у' г 8 ( л = |
89,7957л 2 — 563,2818л3 + |
1542,021л:4 - |
|||||
- |
1929,021л3 + |
718,2745л6 + |
77,1585л:- + |
1066,516л8- |
|||
|
|
- |
1620,031л:9 + |
618,5676л10. |
(11.27) |
||
Соответственно зависимости |(х ) |
и |
£(*) |
определяются выра |
||||
жениями |
|
|
|
|
|
|
|
5 (л) = 179,5914л — 1689,8454л;2 + 6168,0840л:3 -
-9645,1050л4 - f 4309,6470л3 + 540,1095л6
+8532,1280л7 - 14580,2790л8 + 6185,6760л9. . (11.28)
С(л) = 29,93190л3 — 140,82045л4 -ф 308,40420л3 —
- |
321,50350л6 + 102,61064л7 + |
9,64481л8 + |
|
+ |
118,50170л9 — 162,0031л1О+ |
56,23342л11. |
(11.29) |
Сравним этот закон с законом движения, сообщающим по ставленному критерию минимум в более широком классе функ ций, а именно в классе функций, допускающих в граничных
32
117 .Зак 3
|
|
|
|
|
|
|
|
|
Т а б л и ц а 2 |
|
|
Коэффициенты полиномиального закона движения без „мягких" ударов |
|
|
|||||
|
|
|
при различных граничных условиях |
|
|
|
|||
О/Э |
«1 |
а2 |
я9 |
ах |
«5 |
Й,1 |
«7 |
|
о.» |
0/0 |
89,796 |
2563,282 |
1542,021 |
-1929,021 |
718,275 |
77,152 |
1066,52 |
-1620,03 |
618,57 |
0/0,05 |
88,081 |
—550,955 |
1504,197 |
-1874,965 |
695,089 |
58,529 |
1075,75 |
-1605,27 |
609,59 |
0,0,10 |
86,366 |
-538,627 |
1466,373 |
—1820,910 |
671,901 |
39,897 |
1084,99 |
—1590,51 |
600,62 |
0/0,15 |
84,651 |
-526,300 |
1428,549 |
-1766,853 |
648,718 |
21,262 |
1094,22 |
21575,74 |
591,65 |
0/0,20 |
82,937 |
-513,973 |
1390,724 |
—1772,797 |
625,529 |
2,6304 |
1103,46 |
—1560,98 |
582,67 |
0,05/0 |
87,021 |
—547,446 |
1502,744 |
—1886,616 |
705,525 |
91,966 |
1003,93 |
21553,79 |
596,61 |
0,10/0 |
84,247 |
-531,620 |
1463,527 |
-1844,377 |
692,980 |
106,70 |
941,240 |
—1487,42 |
574,62 |
0,15/0 |
81,471 |
—515,777 |
1424,216 |
—1801,895 |
680,147 |
121,52 |
878,712 |
—1421,22 |
552,68 |
0,20/0 |
78,696 |
-499,936 |
1384,918 |
-1759,444 |
667,352 |
136,31 |
816,197 |
—1355,03 |
530,74 |
|
|
|
Т а б л и ц а 3 |
|
Инварианты полиномиального закона |
||
|
при однородных |
граничных условиях |
|
X |
|
6 |
с |
0 |
0 |
0 |
0 |
0,05 |
5,4670 |
0,1631 |
0,00934 |
ОЛО |
6,3086 |
0,4703 |
0,01867 |
0,15 |
5,1961 |
0,7620 |
0,05615 |
0,20 |
3,7256 |
0,9844 |
0,09363 |
0,25 |
2,6450 |
1,1414 |
- 0,14937 |
о'.зо |
2,1015 |
1,2581 |
0,20710 |
0,35 |
1,8923 |
1,3572 |
0,29189 |
0’40 |
. 1,6923 |
1,4474 |
0,37667 |
0 45 |
1,2351 |
1,5220 |
0,43553 |
050 |
0,4259 |
1,5650 |
0,49438 |
■0*55 |
—0,6300 |
175604 |
0,57194 |
0^60 |
-1.6818 |
1,5020 |
0,64950 |
0,65 |
—2,4595 |
1,3969 |
0,71911 |
0,70 |
—2,8337 |
1,2630 |
0,78871 |
0,75 |
—2,9513 |
1,1180 |
0,84455 |
0,80. |
—3,2529 |
0,9648 |
0,90038 |
0,85 |
—4,2515 |
0,7808 |
0,93883 |
0',90 |
-5.9068 |
0,5278 |
0,97727 |
0,95 |
- 6 3912 |
0,2076 |
0,98864 |
1,00 |
0 |
0 |
1 |
Рис. 2. Корректировка закона движений с разрывным гра фиком ускорений достаточно гладкой функцией.
34
точках отрезка [0, 1] разрывы непрерывности инварианта уско
рений («мягкие» удары). Такой закон был синтезирован в § 1, гл. II настоящей работы, где показано, что в случае постоян ной скорости ведущего звена он вырождается в известный за кон равноубывающего ускорения, для которого максимальное значение инварианта скорости 6тах=1,5, а максимальное зна
чение инварианта ускорений |шах= —£min = E(0) = 6. |
|
В настоящем же |
законе (11.27) — (11.29) коэффициенты ма |
ксимальной скорости |
и ускорения равны (см. табл. 2): 6тах = |
= 1,5650, £тах= 6,3086, s,„in=—6,3912. |
|
Таким образом, в настоящем законе ценой увеличения ко эффициента максимальной скорости и средних ускорений на 4,3% и коэффициента максимального ускорения на 6,5% уда лось ликвидировать скачок ускорений Д£, который в законе
равноубывающего ускорения определялся соотношением |
Д£= |
|
~ ^шах= 6. Численные значения инвариантов |
ускорения |
£(*)» |
скорости б(х) и пути £(х) приведены в табл. |
3. На рис. 2 |
при |
ведены графики безразмерного позиционного коэффициента ускорения %(х), скорости 6 (х) и пути £(х) для базового слу чая а = р = 0.
Полученный закон движения целесообразно применять в сравнительно быстроходных механизмах'при равномерном вра щении ведущего звена в целях уменьшения инерционных на грузок, получения благоприятных углов давления в кулачко вых механизмах, уменьшения крутящих моментов на главном валу.
§ 3. Оптимизация движения по критерию средневзвешенных ускорений ведомого звена
В предыдущих задачах динамически оптимальный закон движения находился из условия равномерной минимизации ускорений ведомого звена на заданном интервале при извест ной скорости ведущего звена. Иногда возникает задача о бо лее выгодном распределении сил инерции по ходу ведомого, звена при одновременном уменьшении сил инерции на всем ходу. Например, при синтезе тяжело нагруженных кулачковых механизмов в зоне удаления (подъема) более выгодным явля ется уменьшение сил инерции в начале подъема, когда усилие замыкающей пружины, усилие трения и силы инерции нагру жают пару кулачок—толкатель. Напротив, в конце участка удаления, когда силы инерции разгружают контактную пару, можно допустить более высокий уровень сил инерции. В этом и в других подобных случаях возникает задача о минимиза ции средневзвешенных ускорений ведомого звена. Полагая, что ведущее звено вращается с постоянной угловой скоростью, для решения поставленной задачи используем форму безраз мерных позиционных коэффициентов пути £, скорости б и уско рения С использованием этих коэффициентов кинематиче
3' |
35 |
ские функции q, q, q ведомого звена на отрезке [0, ф0] можно
записать в виде |
Я£(х), |
- д = |
|
9 = |
-5Lco8(x), |
<Ро
Искомый закон движения 3(л) должен сообщать минимум функционалу оптимальности, который в данном случае имеет вид (
R = ^p( x)?(x)dx, (Н.ЗО)
о
где р (х) > 0 — весовая (управляющая) функция.
Искомый закон движения должен удовлетворять изопериметрическому условию (11.31), которое выражает требование,
чтобы ход ведомого звена был |
равен заданной |
величине qo |
на заданном интервале: |
|
|
J Srfjc = |
1. |
(11.31) |
о |
|
|
Граничные условия в общем случае имеют вид |
|
|
8 (0) = 80, 8(1) = 8,. |
(11.32) |
|
В дальнейшем мы рассмотрим случаи, когда р(х) непрерывна
на рассматриваемом интервале или является кусочно-постоян ной функцией.
Обозначим
&(■*)= У (*). £(■*) = У '(■*)•
Тогда соотношения (Н.ЗО)—(11.32) будут |
иметь вид |
1 |
(11.33) |
R = ^p(x)y'2dx, |
|
О |
|
1 |
|
Jyrfx = l, |
(И.34) |
У(0) = 80, y(l)=8i. |
(11.35) |
1. С л у ч а й н е п р е р ы в н о й в е с о в о й ф у н к ц и и . В
этом случае оптимальный закон движения может быть найден в результате интегрирования уравнения Эйлера для поставлен ной вариационной задачи (11.33) — (11.35). Искомая функция
у(х) должна сообщать безусловный минумум функционалу R*:
1
R * = |
j [р (х)у'2+ Ху] dx, |
(11.36) |
где X — неопределенный |
о |
|
множитель Лагранжа. |
|
36
Уравнение Эйлера для функционала (11.36) имеет вид
i - z - l h - l p w y w 1 = ° . |
(н.з7) |
Находим общий интеграл уравнения (11.37)
Постоянные интегрирования Сь Сг и неопределенный множи
тель А определяются из граничных условий (11.35) и изопериметрического условия (11.34).
Вкачестве примера рассмотрим случай, когда весовая
функция линейно зависит от х: р(х) = 1 + рх. |
Граничные усло |
||
вия полагаем однородными: 6o = 6i= 0, |
что |
соответствует ти |
|
пичной схеме работы цикловых механизмов. |
|
||
После ряда несложных преобразований из соотношений |
|||
(11.34) , (11.35) |
получим |
|
|
г _ 2:П1п (1 +ц) — и] |
р __п |
||
1 |
2,а - (;*+■ 2) In (1 + 1*) ’ |
°2 — и, |
|
у = ___V M L + ! 4 _
*2р- ( ! * + 2) In(!+ (!)•
Находим закон движения в рассматриваемом случае:
8 (■*) = 2^ —0Г+~2) In ('1+ ,х) Iх ln (1 + !А) ln (1 + Р*)] •
Соответственно выражения для инварианта ускорения 1(х) и
инварианта пути £(х) имеют вид
£ = |
2,а — (|х + '2)1п(1 +|х) [ 1п ( 1 + |
I1) |
1 +JJUC ’ |
’ |
= 2i » - ( i i + 2,)iii{ i+ ^ { 'T |
' " ( Н - : * ) - |
|
[in <! + , « ) - - ! - ] } .
Ha рис. 3, а приведены графики инвариантов подобия для ■—0,5 и р = 0, что соответствует условию равномерной оп
тимизации.
2. С л у ч а й к у с о ч н о - п о с т о я н н о й в е с о в о й ф у н к
ции. Рассматриваемый |
отрезок {0, 1] разбит на |
п участков |
|
точками х0= 0, х\, |
Xi, |
хи= 1. , На отрезке |
x ^ ^ x ^ X i |
весовая функция имеет постоянное значение pi. Искомый опти
мальный закон движения по-прежнему удовлетворяет изопериметрическому соотношению (П.34) и граничным условиям (11.35) , причем в данном случае принято 60= 6i = 0. Величина
средневзвешенного квадратического ускорения ведомого звена, которая в соответствии с постановкой задачи является крите
37
рием оптимальности движения, принимает вид функционала
(11.38)
я = 2 |
л / ‘ * d x ' |
(1К38) |
i=1 |
xi-l |
|
где у{(х) = ^ (х )— безразмерный позиционный коэффициент скорости на i-м участке при x t.x^ х < x t. Критерий оптималь
ности на г'-м участке имеет вид
Рис. 3. Оптимальные
а- в случае непрерывной весовой функции: / — для равномерной оптимизации (/> (х) = 1); 2—
дни: 1 — для кривой p (jr) =■ / |
5- 0 < х < ° '5’ г - для кривой |
I |
1. 0,5 <х< 1,0; |
38
R i = |
f |
y ' f d x . |
(П.39) |
|
%! |
|
|
|
xl- 1 |
|
|
На [отрезке [x,_b x,] функция |
_y,(.*:) должна |
удовлетворять |
|
изопериметрическому условию |
|
|
|
*1 |
|
|
(11.40) |
j |
у id x |
= ai , |
|
xi-1
5
39
где af— неопределенный положительный параметр. С учетом соотношения (11.40) уравнение Эйлера на /-м участке для функ ционала (11.39) можно записать в виде
Хг — 2у; = 0,
где X, — неопределенный множитель Лагранжа для /-го участка. Функция y t {x), сообщающая условный минимум функцио
налу (11.39), имеет вид
У,(х)=±1-х* + |
Сих + С21. |
(И.41) |
Постоянные интегрирования |
на Z-м участке определяются |
|
по граничным условиям на этом участке. Полагаем, что пози ционный коэффициент скорости у (л:) на / — 1 участках построен.
Тогда на левом |
конце /-го участка при л = |
функция y t (х) |
|||
в силу |
условии |
непрерывности должна иметь |
значение _уг-ь |
||
которое |
принимает функция у Н1(х) при |
х — х ^ . |
Значение |
||
функции у{(х) |
при x = x t не определено |
для |
всех |
участков, |
|
за исключением |
последнего, на котором |
из граничных усло |
|||
вий следует, что у„(х„) = 0. Для определения постоянных ин тегрирования на /-м участке привлечем кроме условий непре
рывности условие трансверсальности на правом конце |
для всех |
||
Z<1 п — 1. |
|
|
|
В данном случае условие |
трансверсальности имеет |
вид |
|
|
dy'i X- |
= 0 , |
(11.42) |
|
|
|
|
где / T = j ';2+ |
X,jv |
|
|
С учетом |
соотношения (11.41) условие (11.42) может оыть |
||
приведено к виду |
|
|
|
|
■у x i + |
Сц — 0. |
(11.43) |
Условие трансверсальности (11.43) совместно с условием пери одичности (11.44) и изопериметрическим условием (11.40) по
зволяет определить постоянные интегрирования Clh C2i и па
раметр X,-:
|
|
У/ (x i-i) — У;-1, |
|
|
|
(11.44) |
гдеУы — значение функции на левом конце /-го |
участка |
(при |
||||
x = x t_i). Представим соотношение (11.44) |
в виде |
|
|
|||
- у x 2i-\ -)- Сих Н1 -f- C2i = |
v,-_i- |
(11.45) |
||||
Изопериметрическое |
условие (11.40) |
с |
учетом выражения |
|||
(11.41) приводится |
к уравнению |
|
|
|
|
|
“ПГ (*?— x i-i) + |
~2 Си |
— х 1~0 + С-и(x i — *i-i) = |
?•;. |
(11.46) |
||
40
|
Из системы уравнений, образованных соотношениями (11.43), |
|||||||||||
(11.45), (11.46), получим, |
обозначив |
предварительно |
k i—x l—x i _ly |
|||||||||
|
|
Си — — 2>х1Уг-А - < |
|
|
|
|
|
|||||
|
|
C 2i= y i - 1 - |
|
v*2 |
У/-А — Ч qv - у,-А — ч |
|
||||||
|
|
|
Л 1 - \ |
|
|
|
|
|
|
|||
|
|
xf = |
6 |
|
а< |
|
|
|
|
|
|
|
|
|
|
|
д? |
|
|
|
|
|
|
|
|
Тогда |
закон |
движения |
на г-м участке при |
^ .л :^ |
л:г |
при |
||||||
мет вид для |
всех г-<га— 1: |
|
|
|
|
|
|
|||||
|
У< = У/-1 + |
4 - (y" |
|
T af) (Л - |
* м ) {х - |
2х‘+ * '« )• |
(п-47) |
|||||
Из |
граничных |
условий |
(11.35) |
и условия ^ = 0 |
следует, |
что |
||||||
у0= |
0- |
Тогда |
при i = 1 |
|
имеем |
|
|
|
|
|
|
|
|
|
|
|
J'i = |
- |
4 - 7 Г -*(-* — 2JCi). |
|
|
|
|
||
Таким |
образом, |
по формуле |
(11.47) составляется |
оптимальная |
||||||||
функция у{х) |
для всех |
|
— 1. |
На последнем, |
п -м |
участке |
||||||
постоянные интегрирования определяются из следующих усло вий:
Уп (■*„-1) = |
J^-i. У„ (А) = 0. |
(11.48) |
Условия (11.48) в развернутой |
форме с учетом того, |
что -*„=1, |
имеют вид |
|
|
-^ -^ п -\-\-Сых п_^-{-С2п — уп-{, |
(11.49) |
|
4 f- + Cln+ C2„= 0.
Из системы уравнений, образованных соотношениями (11.49) и (П.46) при i = n , определим
Ап — 2 |
Уп- |
л3„ |
■ 2а„ |
■ (1 + ДСд-l) ' |
У л- 1 |
|
|
|
|
|
|
Г __о У п - i A n |
2ал |
. уП- 1 |
> |
||
'-'2п — |
'У |
.я |
|
"*71-1 -Г д |
|
|
>-„=12 |
Уп~ \^ п |
|
||
|
|
|
|
Дп |
|
где А „ = 1 |
— х п. х. |
|
Учитывая приведенные выше соотношения, получим выра |
||
жение для искомого закона движения на последнем |
участке |
|
^ _ x< . v < |
1 |
|
Уя |
= (1 — х) У П -1 ■3 |
(П.50) |
41
