Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Хитрик, В. Э. Методы динамической оптимизации механизмов машин-автоматов

.pdf
Скачиваний:
25
Добавлен:
19.10.2023
Размер:
4.88 Mб
Скачать

где Х; — неопределенные

множители.

 

Необходимое

условие

минимума функции

Ф имеет вид

 

 

 

 

# -

=

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)

получим

 

 

г _ 21п (1 +ц) — и]

р __п

1

2,а - (;*+■ 2) In (1 + 1*) ’

°2 — и,

у = ___V M L + ! 4 _

*2р- ( ! * + 2) In(!+ (!)•

Находим закон движения в рассматриваемом случае:

8 (■*) = 2^ —0Г+~2) In ('1+ ,х) 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

Соседние файлы в папке книги из ГПНТБ