Порядок ввода и состав входных данных представлен в табл.1.
Таблица 1
Идентификатор | Тип | Размерность | Формат | Подпрограмма |
L | Целочисленный массив | 7I10 | FORMMA | |
CX | Вещественный массив | L(2) < 50 | 7F10.0 | FORMMA |
BA | Вещественный массив | L(3) < 2500 | 7F10.0 | FORMMA |
LX | Целочисленный массив | 5I10 | PODXKT | |
H | Целочисленный массив | 7I10 | ENTRD | |
KP | Целочисленный массив | 35I2 | ENTRD | |
GR | Вещественный массив | H(2) < 250 | 7F10.0 | ENTRD |
AR | Целочисленный массив | H(3) < 150 | 7I10 | ENTRD |
MAC | Вещественный массив | H(3) < 150 | 7F10.0 | ENTRD |
BMP | Вещественный массив | H(6) < 250 | 7F10.0 | BMPBMH |
BMH | Вещественный массив | H(7) < 600 | 7F10.0 | BMPBMH |
KPR | Вещественный массив | - | 7F10.0 | REG1 |
KPA | Вещественный массив | - | 7F10.0 | REG1 |
KF | Целочисленный массив | < 7 | 7I10 | ENTRD |
FORM | Вещественный массив | < 28 | 7F10.0 | ENTRD |
AY | Вещественный массив | < 21 | 7F10.0 | ENTRD |
ADR, NEW | Целочисленная и вещественная переменные | I10 | ENTRD | |
KDR | Целочисленный массив | F10.0 | ENTRD |
Состав входных данных математической модели приведен в [1]. При решении данной задачи он имеет следующие особенности:
1) Первый элемент целочисленного массива "Н" всегда равен единице, т.к. проведение завязки предусмотрено на одном режиме.
2) Программа управления содержит один закон, т.е. массив "ВМР" состоит из трех элементов - номера программы, количества элементов в массиве "ВМР" (равного 3) и номера закона управления.
3) Вещественные массивы "КPR" и "КРА" из состава входных данных исключаются, т.к. используется простая программа управления.
Входные данные задачи
Порядок ввода входных данных представлен в табл.2.
Таблица 2
Идентификация | Тип | Размерность | Формат | Подпрограмма |
KZ, KY | Целочисленные переменные | 2I 10 | ENTRZ | |
Дополнительная информация | - | - | - | ENTRZ |
NY1, ix1NY2, ix2NYky, ixky | Целочисленные переменные | 2I10 | ENTRZ |
2.2.1. Целочисленная переменная "KZ" определяет тип выполняемого расчета. Она может принимать следующие значения:
KZ = 0 - выполняется завязка методом вариантного расчета;
KZ = 1 - завязка с поиском оптимального значения одного из параметров;
KZ = 2 - завязка с многомерной оптимизацией.
Целочисленная переменная "KY" определяет количество узлов из массива "СХ", требующих подключения типовых характеристик после выполнения завязки.
2.2.2. Дополнительная информация в зависимости от значения целочисленной переменной “KZ” имеет следующее значение и состав:
При KZ = 0 в состав входных данных не включается.
При KZ = 1 состоит из четырех переменных, вводимых на одной перфокарте по формату 2I10, 2F10.0.
iA - программный адрес параметра, минимальное или максимальное значение,
которого определяется при проведении расчета. При отыскании минимума
имеет знак "минус";
iV - программный адрес варьируемого параметра;
VMAX - максимальная граница изменения варьируемого параметра;
VMIN - минимальная граница изменения варьируемого параметра.
При KZ = 2 состав и порядок ввода дополнительной информации представлен в табл.3.
Таблица 3
Идентификатор | Тип | Размерность | Формат | Подпрограмма |
KO | Вещественный массив | 4I10 | ENTRZ | |
LF | Целочисленный массив | KO(2)+1< 50 | 7I10 | ENTRZ |
AYT | Целочисленный массив | KO(1) < 200 | 7F10 | ENTRZ |
XMA | Вещественный массив | KO(1) | 7F10.0 | ENTRZ |
XMI | Вещественный массив | KO(1) | 7F10.0 | ENTRZ |
FO | Вещественный массив | KO(1) | 7F10.0 | ENTRZ |
Целочисленный массив "КО" содержит информацию о размерности решаемой задачи.
КО(1) - количество переменных оптимизации;
КО(2) - количество функциональных ограничений;
КО(3) - максимальное количество шагов поиска;
КО(4) - период печати результатов поиска.
Целочисленный массив "LF" имеет следующий состав:
LF(1) - программный адрес параметра, являющегося функцией цели. При поиске минимума задается со знаком "минус".
Последующие КО(2) элементов массива содержат программные адреса параметров, используемых в качестве функциональных ограничений. Если задано ограничение типа "не более заданного значения", то соответствующий программный адрес задается со знаком "минус".
Целочисленный массив “AYT” содержит программные адреса переменных оптимизации (варьируемых параметров).
Вещественные массивы "ХМА" и "XMI" содержат, соответственно, максимальные и минимальные границы изменения переменных оптимизации.
Вещественный массив "FO" содержит значения параметров, принятых в качестве функциональных ограничений. Порядок их следования должен совпадать с порядком следования программных адресов в массиве "LF", начиная со второго.
2.2.3. Целочисленные переменные "NY" и "iX" вводятся отдельной строкой для каждого из узлов, требующих подключения типовой характеристики.
NY - условный номер узла в массиве "СХ";
iX - признак подключаемой характеристики.
Признак подключаемой характеристики может принимать следующие значения:
Для узла типа "Компрессор"
iX = 0 - подключается типовая характеристика осевого нерегулируемого компрессора;
iX = 1 - подключается типовая характеристика одноступенчатого центробежного компрессора.
Для узла типа "Разделитель потоков"
iX = 1 - подключается характеристика sв= f(lв) 2);
iX = 2 - подключается характеристика sс= f(lс);
iX = 0 - подключаются обе характеристики.
Для узла типа "Основная камера сгорания"
iX = 1 - подключается характеристика s = f(l);
iX = 2 - подключается характеристика h= f(a);
iX = 0 - подключаются обе характеристики.
Для узла типа "Турбина"
iX = 0 - подключается типовая характеристика осевой нерегулируемой турбины;
Для узла типа "Форсажная камера"
iX = 1 - подключается характеристика sг= f(l);
iX = 2 - подключается характеристика hф= f(aS);
iX = 0 - подключаются обе характеристики.
Для узла типа "Выходное устройство"
iX = 1 - подключается характеристика j = f(pс*);
iX = 2 - подключается характеристика m = f(pс*);
iX = 0 - подключаются обе характеристики.
Для остальных узлов полагается iX = 0.
__________________________
2)sbи sc- коэффициенты восстановления полного давления, соответственно, во внутреннем и наружном контурах.
ВЫХОДНАЯ ИНФОРМАЦИЯ
3.1. Результатом выполненной задачи формирования облика двигателя являются его основные параметры, хранимые в массиве "ВА". Он содержит уточненные в процессе завязки значения искомых параметров, заданных ранее в нулевом приближении, а так же подгруппы, содержащие типовые характеристики отдельных узлов расчетной схемы двигателя. Узлы, требующие подключения типовых характеристик, определяются пользователем.
Данный массив записывается в базу данных с индексом L(1) и идентификатором iD = 1.
3.2. Печать результатов расчета и условная печать осуществляются в соответствии с заданными массивами "КР", "АR", "МАС" и "КОR".
3.3. При выполнении завязки с оптимизацией (КZ = 2) предусмотрен вывод дополнительной печати, информирующей о ходе поиска.
В процессе поиска с периодом КО(4) шагов осуществляется печать следующих величин.
1) Порядковый номер - NS;
2) Условный номер сделанного шага, т.е. номер проработавшего алгоритма - NO;
3) Условный номер очередного шага - NZ;
4) Общий масштаб поиска - FMO, определяющий величину шагов поиска. Он служит критерием окончания поиска. В начале поиска он задается равным единице и постепенно уменьшается до предельно допустимого - FMO = 0,001;
5) Константа поиска - FKP, которая определяет величину случайного вектора при вычислении очередного шага поиска. Значение “FKP” изменяется при поиске в пределах 0,4...0,85.
6) Одиннадцать признаков, свидетельствующих о возникшей ситуации - М1...М11.
М1 - признак способа проверки функциональных ограничений:
М1 = 0 - учет ограничений в явном виде;
М1 = 1 - поиск методом штрафных функций;
М2 - признак качества шага:
М2 = 1 - удачный.
М2 = -1 - неудачный
М2 = 0 - выход за границу;
М3 - счетчик числа неудачных шагов подряд при поиске в области;
М4 - признак величины “овражных” шагов;
М4 = 0 - шаг "единичной" длины;
М4 = 1 - увеличивающиеся шаги;
М4 = -1 - уменьшающиеся шаги;
М5 - счетчик шагов по алгоритму № 3 ("случайный шаг вне области");
М6 = 0
М7 - счетчик общего числа шагов по алгоритмам № 3 и № 4 (алгоритм № 4 - "экстраполяция в область");
М8 = 0
М9 - признак пустой области:
М9 = 0 - точка в области поиска не найдена;
М9 = 1 - точка в области уже найдена;
М10 - признак прекращения поиска:
М10 = 1 - при поиске выполнено заданное предельное число шагов;
М10 = 2 - при ошибке во входных данных, т.е. когда начальная точка находится вне границы изменения переменных оптимизации, или верхняя граница меньше нижней;
М11 = 0.
7) Коэффициент интерполяции – FLMIN при выходе за границы изменения переменных оптимизации.
8) Значение функции цели в текущей, лучшей, старой, запасной, оптимальной точках.
9) Печать массивов координат (т.е. значений переменных оптимизации) и значений функциональных ограничений в указанных выше точках (по столбцам). В правом столбце печатается массив индивидуальных масштабов для каждой переменной.