Файл: Лебедев И.В. Элементы струйной автоматики.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.04.2024

Просмотров: 238

Скачиваний: 10

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Исходная конфигурация элемента выбирается на основании гидродинамических исследований, расчетов, интуитивных пред­ ставлений и т. п. Единственным требованием, которому должна удовлетворять эта конфигурация, является возможность выпол­ нения элементом основных функций и возможность определения величин параметров, входящих в условия работоспособности. Так, для дискретного элемента достаточно, чтобы он переклю­ чался при подаче сигнала управления и возвращался в исходное состояние после снятия сигнала.

Понятно, что построить конфигурацию, удовлетворяющую таким общим требованиям, можно без больших затрат времени, методом проб и ошибок. Желательно, конечно, чтобы эта исход­ ная конфигурация была возможно ближе к оптимуму, так как это позволит существенно уменьшить объем экспериментов, не­ обходимых для достижения оптимума. Поэтому, чем полнее изу­ чена картина течений в элементе, тем меньше опытов необходи­

мо для достижения оптимума.

 

Варьируемые факторы выбираются

на основании данных

о механизме явлений. Роль различных

факторов уточняется

в ходе факторного эксперимента на основе анализа значимости факторов с помощью критерия Стыодента. Незначимые факторы отбрасываются. Составление матрицы планирования описано выше. Отметим, что для получения оценок только линейных чле­ нов целесообразно ставить две дробные реплики, отличающиеся одна от другой знаками.

Интервалы варьирования 7; должны быть достаточно малы­ ми, для того чтобы получить более точную линейную аппроксима­ цию. С другой стороны, они должны существенно превышать возможную ошибку установки положения пластин.

После того, как проделана эта предварительная работа, мо­ жно приступать к реализации факторного эксперимента.

Особенности статистического анализа. Особенностью фак­ торного эксперимента при оптимизации элементов струйной ав­ томатики является то, что с ошибкой измеряются не только зна­ чения функций отклика, как это обычно принимается, но с ошиб­ кой устанавливаются и значения факторов (положения пластин). Однако, как показал Бокс [75], для линейной модели можно поль­ зоваться обычными методами статистического анализа. При этом ошибка опыта, определенная из повторных измерений зна­ чения функции отклика, будет учитывать не только ошибку из­ мерения функции отклика, но и ошибки установки уровней фак­ торов.

7. П рим ер получения функций отклика

Рассмотрим пример получения функций отклика, связавающих параметры струйного элемента с внутренней обратной связью с его геометрическими размерами. Ниже эти функции

22 за к. 935

337


будут использованы для решения задач оптимизации струнного элемента по нагрузочной способности.

Пример. И сходную конфигурацию выбираем на основании предваритель ­

ных исследований

(рис. 157, б): Ь„= 0,4 мм; и\=

11°;

CI2

=

10°;

«і

=

0,23 мм;

а 2 0,4 мм; I= 4

мм; R =

0,3 мм.

 

 

 

 

 

 

 

И сследование

рабочего

процесса элементов

(см.

гл.

I ll

и

V)

и

предва ­

рительные эксперименты позволили установить, что на парам етры элемента

наиболее

существенно влияю т

следую щ ие факторы (рис.

157,

а ):

положение

разделителя,

характеризуемое координатами

хі, уі, ширина

приемных каналов,

определяемые

смешением і/ 2

и

ysпластин 2 и 8 ; ширина

атмосферного

к ан а ­

ла, за д а в а е м а я координатой

х 2, смещением

у3,у7стенок. Эти

семь

факторов

выбираем в качестве варьируемых.

 

 

 

 

 

Д л я

получения

линейных

приближений

функций отклика

д о л ж н а

быть

реализована дробная

реплика

2 '~ 4, позволяю щ ая определить

коэффициенты

при семи линейных членах. Учитывая то, что точность установки пластин не превыш ает 0,005—0,01 мм, принимаем интервалы варьирования всех ф акто ­ ров л = 0,04 мм. Таким образом, верхнему уровню значений факторов будет

отвечать смешение топ или иной пластины

на

0,04 мм

в

положительном

на ­

правлении

соответствующей

оси

координат,

а

нижнему — на

0,04 мм

в отри ­

цательном

направлении топ ж е оси.

 

 

 

 

 

 

 

 

 

 

М атрицу планирования

(табл. 12)

получаем

следующим образом .

С о ­

ставляем

матрицу

планирования

полного

факторного

эксперимента

23.

Д л я

этого

в первые

три колонки

записываем

все

возм ож ны е

сочетания

значений

-г 1

и

— 1

для

.чу,

х2,

х3 (для

простоты

цифры

1

опущены).

Следующие

ко ­

лонки

получаем перемножением

первых

трех

колонок;

Х]Х2,

Х|.ѵ3, х 2 х а, .ѵ 1 2 Х3 .

Эти

обозначения

взаимодействий

заменяем

новыми переменными; х\= ХіХ2,

х 5 =

Х1 Х3 ,

х 6 =

х2 х3, X? = Х|Х2 х3. К аж д о й

строке матрицы планирования соот­

ветствует

один

опыт,

т. е.

определенное

сочетание

положений

пластин.

В каж до м

опыте

определяем

значения

нагрузочных способностей К„л, К дл.

Д ив, давления отпускания рот, остаточного давления р0ст и границы коле­

баний

Q]jCT

. Д л я

определения

коэффициента нагрузочной способности изме­

ряем

расход

на

выходе

Q B

при

минимальном давлении

единичного уровня

давлений р>У) =

40

мм

вод.

ст.

и расход срабаты вания

Q cp. К оэффициент

нагрузочной способности вычисляем как отношение расхода на выходе к рас ­ ходу срабатывания;

 

^ИА = ^И тах/^А ср’

*ДА =

тах/^Аср!

*НВ = QjT шах/^Вср•

 

где <3И т ш ,

Яд ш а х — расход на вы ходах

И и

Д

соответственно;

Q A C P ,

Q BCP — расход срабаты вания

по в х о д а м /1 и В соответственно.

 

 

 

Д авлени е

срабаты вания

во всех

опы тах

поддерж ивалось

постоянным

Pop =

30 мм

вод. ст. за счет

регулирования сопротивления

дросселя в к ан ал е

управления,

а давление питания поддерж ивалось

равным

150

мм

вод.

ст.

Результаты

измерений

указанны х

парам етров

в

каж до м из

опытов

за п и ­

саны в соответствующ их граф ах табл.

1 2 . В нижней

части

таблицы приведе­

ны значения

коэффициентов уравнений связи д л я

каж дого

из

параметров,

вычисленные

по

формулам

(435). Д л я

оценки

уровня

значимости

коэфф ици ­

ентов

Ь0,Ьі,...,

6 7 после каж до го опыта устанавливаем

исходную

конф игура ­

цию и

измеряем

значения

параметров. П роделав таким

образом

10 повторны х

опытов, определяем в центре эксперимента среднее арифметическое и среднее

квадратичное отклонение полученных значений

для каж до го парам етра

(см.

табл.

12). З атем

по формуле (438)

определяем

уровень

значимости коэф ф и ­

циентов. Значение

критерия Стыодеита tвыбирается по

статистическим

т а б ­

лицам

д л я 5 % -ной доверительной вероятности

при девяти степенях свободы .

После отбрасы вания незначимых членов в соответствии с формулой

(438)

получаем следующие вы раж ения для

функции отклика:

 

 

338


&ид = 3,43 + 0, lX] — 0, 12х2 + 0 ,15х.,— 0,075х5 — 0 ,15.ѵ6;

Ад А = 4,17 — 0,23х2 + 0,0Эл-4 — 0,11х5 — 0.09х6;

ротп = 5,2 + 0,75л-! + 0,25х2 — 0,25х3 — 0,25х., + 0,25xä —

— 0,25х6 + 0,25х7;

(440)

QycT =

13,75 + 6,75х[ + 2 ,5 х 2— 1,2 5 х 3— 1 , 5хц — 2 х ; ;

Р ост =

1 + 0,5Х[ + 0 , 5х3 + 0 ,5 х 4 + x ä +

0 , 5х6;

АііВ = 2 . 5 — 0,09хі — 0 ,З х2+ 0 , 2 5 х4— 0 ,2 х5.

Чтобы оценить линейные члены отдельно

от парных взаимодействии,

была реализована м атрица планирования, в которой знаки изменены на об­ ратные по сравнению с матрицей табл. 12. Полученные при этом коэфф ици ­

енты отличались от коэффициентов, полученных

в

табл.

12

на

величины,

меньшие уровня значимости. Это служ ит подтверждением

того,

что

коэф ф и ­

циенты системы являю тся оценками только линейных членов.

 

 

 

 

Функции

отклика

(440) позволяю т оценить

влияние различных

факторов

на изменение

того или

иного параметра. Н и ж е

эти

функции

отклика будут

использованы для решения задачи оптимизации

элемента

 

по

нагрузочной

способности.

 

 

 

 

 

 

 

 

8. Метод решения задач оптимизации по рабочим параметрам

Задача оптимизации по рабочим параметрам формулируется следующим образом: отыскать в факторном пространстве точку, принадлежащую области работоспособности и характеризую­ щуюся максимальным значением критерия качества. Как говори­ лось выше, основной особенностью этой задачи является нали­ чие ограничений (условий работоспособности), которые необхо­ димо учитывать при выборе направления крутого восхождения.

Для отыскания направления крутого восхождения наиболее удобен алгоритм выпуклого программирования [19], позволяю­ щий рассчитать направление крутого восхождения, ведущее строго в глубь области работоспособности. Допустим, выбрана исходная конфигурация струйного элемента. Выбрав варьируе­ мые факторы и интервалы варьирования, реализуем факторный эксперимент.

В результате проведения факторного эксперимента получаем линейные приближения функций отклика в виде:

П

П

А = й і + 2

в і = ь й' і + 2 Ь ц Х і '

£-1

£=І

Подставив функции отклика в условия работоспособности, можно записать

Ф, = ± ( в , - в ; )= Ь'оі +

> 0.

!» 1

22*

339



ю

н

е?

с;

\£>

,. X о л

>>а

О*

с

ь

ft.о

К(

ю

2,9

,62

СО

,62

,33

2,2

2,0

 

 

 

 

 

 

 

 

 

со" см

 

 

 

 

 

 

о

о

о

о

о

см

о

Tf

СО

о

О0

о

оо

 

 

--,

см

 

СМ

 

 

со

■чт

о

о

о

-D

со

■cf

со

 

Г--

_

г-

СО

со

со

00

ю

со

-■

ю

—1

со

■Cf

со

LO

 

■D

-?f

со

Tf

 

со

со

<

с»

CD

см

ІО

о

ю

•cf .

тг

СО

•—1

со

00

со

о

со

 

со

СО

со

со

СО

со

со

со

 

 

 

 

 

+

I

I

_|_

 

 

 

 

 

 

 

-j_

О

 

:

+

I

+

О

 

 

 

 

 

о - к -f

+

+

 

 

+

I

I I

+

+

+

+

°- * I I + +

 

 

_|[_

о

 

 

 

 

 

°- * I +

і +

I +

I +

 

+

+

+

+

+

0,075

0,05

1

 

U0

■4f

о

о

со

со

о

см

о

о

со

со

о

о

о

о

г-

ю

о

о

о

о

Ю

СО

Ю

О

О

О

-М О

О О О <М О

О

!I

U0

CD

ю

НО ю

о

о

см

о

о

о

о

о

о

I"-

I

смI

!

 

см

_

ІО

ю

о

о

о

 

 

о

о

о

о

о

ю

CD

ю

 

ю

. со

о

<М Ю

см

о

о

о

о

о

 

оI

 

 

 

о

м

со

ю

 

 

I

 

 

со

СМ (М

 

 

 

 

 

 

 

о

о

ГМ ю

ю

ю

CD

о

h-

f-

о

о

о

о

СО о

о

со

со

см

ю

о

ю

■^f

О)

 

со

со

ю

со

см

•е*-

о

 

 

 

 

 

 

 

 

Си

 

 

Си

 

 

 

 

 

 

 

 

с

 

 

 

о

 

 

 

 

 

 

 

 

 

н S

 

 

 

 

 

 

 

 

 

 

 

*5

ю

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

оа

со

Tf

ю

со

is:

00

 

 

 

а

 

 

 

си

?-

fr-

fr-

fr-

{-

fr-

н

н

 

 

о

<ь>

3

cf)

 

ю s

н

3

3

3

3

3

3

3

's

 

с

с

с

с

с

е

о

СО CQ Оз £Q «з з

 

о

о

о

о

о

о

о

о

 

 

о