CC BY
32
ФУНДАМЕНТАЛЬНЫЕ ОСНОВЫ ПРОБЛЕМ НАДЕЖНОСТИ И КАЧЕСТВА
УДК 621. 396
/
К ПРОБЛЕМЕ СИНТЕЗА ЦЕЛЕВОЙ ФУНКЦИИ ПАРАМЕТРИЧЕСКОЙ ОПТИМИЗАЦИИ СЛОЖНЫХ ТЕХНИЧЕСКИХ СИСТЕМ
А. В. Саушев
Введение
Оптимизация сложных технических систем (СТС) на этапе параметрического синтеза предполагает решение двух основных задач - определение номинальных значений внутренних параметров системы и допустимых пределов их изменения. Внутренние параметры - это параметры элементов ТС, которые характеризуют состояние и свойства самой системы. При проектировании они определяют вектор X управляемых (варьируемых) параметров. Математическая функциональная модель СТС представляет собой алгоритм вычисления вектора выходных параметров У при заданных векторах внутренних параметров X и внешних параметров V. Внешние параметры характеризуют свойства внешней по отношению к СТС среды и оказывают влияние на ее функционирование. Выходные параметры характеризуют свойства ТС, интересующие потребителя. Они представляют собой параметры-функционалы, т.е. функциональные зависимости фазовых переменных СТС, и параметры, являющиеся граничными значениями диапазонов внешних переменных, в которых сохраняется работоспособность системы. К выходным параметрам на стадии параметрического синтеза относятся показатели назначения, параметрической надежности и экономичности [1, 2].
Показателем параметрической надежности при ограниченных статистических данных о законах распределения внутренних параметров СТС во времени является запас работоспособности [1, 3, 4]. Область работоспособности О = Р П М определяет множество допустимых значений внутренних параметров, при которых выполняются все требования, предъявляемые к выходным параметрам СТС [5]. Эта область определяется условиями работоспособности, которые в случае двухсторонних ограничений на параметры имеют следующий вид:
У ■ < У = Е(X) < У / = Гт'
А ] Ш1П —А у ' — ] шах'7 (1)
X ■ < X < X I = Гп
шах''
где У/ шах (XI шах), У/ш1п (XI ш1п) - соответственно максимально и минимально допустимые значения /-го выходного У'(/-го внутреннего X) параметра; ^ - оператор, устанавливающий связь между внутренними и выходными параметрами; Ви Р - допусковые области, определяемые соответ-
ственно первым и вторым неравенствами (1). При этом области D в пространстве внутренних параметров будет соответствовать допусковая область M.
Трудность решения задач параметрического синтеза обусловлена тем, что эти задачи приходится решать в условиях многокритериальности [1, 6]. Кроме того, область работоспособности, определяющая допустимые пределы изменения внутренних параметров, имеет, как правило, достаточно сложную конфигурацию.
В статье рассматривается поисковый метод оптимизации, который предполагает, что каждая из функций-ограничений неравенства (1) Yjmax - Fj (X) > 0 и Fj (X) - Yjmin > 0 аппроксимирована конечным множеством линейных гиперповерхностей fj, и допусковая область M задана в виде следующей системы неравенств:
2m n
X fj (X) > 0, fj (X) = bj0 + X b]1X1 > 0 .
j=1 1=1
Методы аппроксимации области работоспособности системой линейных неравенств рассмотрены, например, в работах [5, 6]. В отличие от известных методов оптимизации [6] предлагаемый метод позволяет сформировать такую целевую функцию, для которой возможно применение любого известного алгоритма поисковой оптимизации. При этом достигается максимально возможный или заданный запас работоспособности СТС.
Выбор целевой функции. Рассмотрим методологические аспекты выбора целевой функции при оптимизации СТС. В основу методологии положен аксиоматический принцип.
Постулат 1. Основополагающим постулатом методологии является констатация того факта, что любая СТС характеризуется двумя обобщенными параметрами - эффективностью (полезностью) и затратами (платой за полезность). Постулат логически вытекает из фундаментального философского закона диалектического единства противоположностей. Следствием постулата является формулировка двух критериев оценки СТС - эффективности (Q1) и цены (Q2). Критерий цены характеризует СТС значением параметра Q2 при фиксированном значении параметра Q1, а критерий эффективности - значением параметра Qi при фиксированном значении параметра Q2.
Постулат 2. Обобщенные параметры СТС Q1 и Q2 могут быть представлены в виде функциональных зависимостей от ее выходных и ресурсных параметров, под которыми понимаются параметры системы, определяющие величину затрат и характеризующие потребление различных ресурсов при ее создании и эксплуатации [7-11].
Постулат 3. Два варианта S1 и S2 решения задачи оптимизации считаются равно эффективными Q1(S1) = Q2(S2) лишь в том случае, когда равны попарно их соответствующие выходные параметры Yj(S1) = Yj(S2). На основе постулата 3 функция Q1(Y) конкретизируется в форме вектора.
Постулат 4. Ресурсы СТС тождественны товарам. При этом любой товар как экономическая категория принципиально имеет цену, и она может быть установлена. Следствием постулата является конкретизация функции цены в форме Q2 = YbT, где b - стоимость (цена) единицы ресурса (товара) T. Эта формула по своей сути принципиально отличается от похожей по форме аддитивной функции записи обобщенного показателя качества, которую вводят с целью сведения к скалярной векторную оценку СТС [6, 12].
Задача оптимизации критерия цены налагает на обобщенный параметр Q1 условие его постоянства. Единственно возможным способом удовлетворить это требование является фиксирование всех выходных параметров СТС, т.е. задание определенных потребительских свойств системы. При этом само выражение и, следовательно, методология оптимизации инвариантны относительно величины фиксированного значения цели системы и определяют лишь форму ее задания.
Задача оптимизации критерия эффективности Q1 предполагает нахождение его максимума. Вместе с тем операция оптимизации вектора не имеет смысла. Единственным логически возможным способом разрешения этого противоречия является фиксация всех компонент вектора,
кроме одного, который и подлежит максимизации. При этом Q1 = i^Yq, Yj = constj, j = 1, m, j Ф q ,
где q - нефиксированный выходной параметр СТС.
Выражение для критерия Q1 накладывает на параметр Yj, j Ф q, только требование фиксиро-ванности, но оставляет полную свободу в выборе его конкретных значений.
Анализ литературных источников показывает, что для большинства СТС на первое место выдвигается требование высокой надежности. Применительно к решаемой задаче это означает, что в качестве параметра Yq целесообразно выбрать вероятность безотказной работы или запас
работоспособности системы.
Непосредственное использование важнейшего показателя параметрической надежности -вероятности безотказной работы в качестве целевой функции при оптимизации внутренних параметров системы не всегда эффективно вследствие его малой чувствительности [6] вдали от границ области работоспособности и большой трудоемкости вычислений. Кроме того, при отсутствии статистических данных о распределении параметров элементов СТС вероятность безотказной работы принципиально не может быть использована в качестве такой функции.
Целевой функцией при решении задачи параметрической оптимизации СТС предлагается выбрать запас работоспособности X(X) или минимальный запас работоспособности Xj (X), т.е. minXi(X). Эти критерии позволяют получать в отличие от других известных критериев оптимизации любое Парето-оптимальное решение [6, 13, 14].
Доказано, что для разных законов распределения выходных параметров максимизация выбранной целевой функции, т.е. максимум минимального запаса работоспособности
X(X0) = max X(X) = max, (2)
XeG XeG
дает близкую к максимальной вероятность безотказной работы системы [6].
Решение, получаемое по алгоритму (2), является единственным и в максимальной степени обеспечивает выполнение всех условий работоспособности (1). При этом автоматически учитываются показатели назначения СТС, стоимость ее изготовления и, кроме того, чувствительность и возможные уходы выходных параметров, связанные с вариациями параметров комплектующих элементов в процессе изготовления, хранения и эксплуатации системы.
Для некоторых ЭТС в силу специфики их работы в качестве целевой функции целесообразно выбрать не запас работоспособности X(X) или его минимальное значение Xt (X), а какой-либо показатель назначения. Однако при этом одним из критериев, на которые накладываются ограничения, как правило, должен быть требуемый уровень запаса работоспособности.
Формирование целевой функции. В пространстве Rn внутренних параметров введем метрику l, которая является функцией координат двух любых точек этого пространства, например
rn 2
точек A и B. При этом l = (A) - Xt(B)) , где Xt(A), Xi(B) - координаты векторов то-
чек A и B соответственно; ^ - нормирующий множитель по i-й координате параметров X. Если одна из точек, например точка A, является граничной точкой области работоспособности, а точка B находится внутри этой области и ее координаты характеризуют состояние СТС в рассматриваемый момент времени, то данная метрика будет определять запас работоспособности системы и служить критерием определения координат оптимальной точки.
Для формирования целевой функции представим область работоспособности в виде единого аналитического выражения. С этой целью введем в рассмотрение и воспользуемся свойствами логических R-функций [5].
Пусть Y = F(X) есть функция, определенная всюду в пространстве Rn. Согласно определению R-отображения [5] данная функция является R-функцией, если в каждой из областей Hj (j = 1, 2, ..., 2n) она сохраняет постоянный знак. При этом область Hj представляет собой совокупность всех точек пространства Rn, для которых хотя бы одна координата Xi равняется нулю. В результате использования R-функций область работоспособности может быть задана следующим неравенством:
d
G = (...((Ф1 a£i Ф2) л^2 Фз) л^з ...) g-1) 9d = Akag фg > 0, (3)
g=1
где ag, g = 1, d - величины, принадлежащие интервалу age [-1; 1].
В том случае, когда все ограничения (1) являются двухсторонними d = 2(m + n). При этом Фg = Фg (X) = fj (X), g = j = 1,2m . Функции-ограничения Ximax - Xi > 0 и Xi — Ximin > 0 описываются гиперплоскостямиf, ф (X) = f i(X), g = 2m, d, i = 1,2n .
Для построения Л-коньюнкции удобно воспользоваться формулой [14, 15]:
Ф ЛаФ2 = 0,5( +Ф2 — ^/ф? +ф2 -2«Ф1Ф2 )л(ф1,Ф2), (4)
где Л(фь ф2) - функция, обеспечивающая наличие k производных Л-коньюнкции.
В том случае, если не требуется, чтобы Л-коньюнкция была дифференцируема, формула (3) может быть упрощена. Принимая а = 1 , получим
Ф1 ЛФ2 = 0,5 ( + Ф2 — |Ф1 -Ф21). (5)
В формуле (3) могут быть опущены скобки, и конечный результат не будет зависеть от последовательности свертки.
На основании использования свойств Л-функций можно записать: G=M л^ P > 0 . Принимая значение а = 1, а также следуя формулам (3) и (5), получим аналитическое описание области работоспособности в виде следующего рекуррентного соотношения:
G = 0,5( M+P—| M+P|>0; (6)
M = M 2m = 0,5 ((2m—1 + Ф2т) — |M2m —1 — Ф2т)|);
M2m —1 = 0,5 ((2m—2 + Ф2m —1 — |M2m—2 — Ф2m —1I);
< Mg = 0,5(—1 + Фg — \Mg—1 —ФgI);
M 2 = 0,5 ((M1 + ф 2 — \M1 —Ф21);
M1 =ф1.
Допусковая область P описывается аналогичным образом, только функции ограничения будут иметь вид фg = фg (X) = f(X).
В случае если m = n = 2, d = 2(m + n), а внутренние параметры заданы в относительных единицах, причем X 1min = X2min = —1, X1max = X2max = 1, области M и P запишутся в виде следующих неравенств:
M = 0,25(7jmax + Y2max — Y1min — Y2min — |2F1(X1,X2) — Y1max — Y1min | — — I2F2(X1,X2) — Y2max — Y2min | — |Y1max + Y2min — Y1min — Y2max + + I2FL(^^ 1,X2) — Y2max — Y2min| — |2F2(X1,X2) — Y1max — Y1min||) > 0;
P = 0,5(2 — IX1I — IX 2\ —I X 2\ — I X1|) > 0.
Используя основное свойство Л-функций, заключающееся в том, что логические операции и простейшие арифметические операции над Л-функциями образуют новую функцию, которая также принадлежит к классу Л-функций, можно заключить, что аналитические описания областей M, P и G также будут являться Л-функциями.
Получим уравнение границы области G^, расположенной эквидистантно области G и внутри нее. При этом множество граничных точек областей G и Gß будут располагаться относительно друг друга по направлению градиента к функции G на одну и ту же величину l.
Рассмотрим две граничные точки N е /. (X) е М и N"е /" (X) е М". Координаты точки Ым определяются выражением
х--х+1.-х+/1 • -^(х) . Откуда
/" (X) - / (X) - вгаа/. (X)/. (7)
Уравнение границы области М" на основании использования свойств ^-функций примет следующий вид:
м" - (X) - Т^з: (X)+/1^1(X) - /" (X) - /'^+1(x)|)}-о, /".(X) = /" (X). (8)
В том случае, если граничная точка принадлежит области Р, координаты точки N определяются выражением X" - Х1 ± / . Аналитическое описание областей Р" и О" производится аналогично описанию областей Р и О по формуле (6).
Доказано, что при любом числе ^-функций ф^ (X), g - 1, й значение ^-функции О, аналитически описывающей границу области работоспособности в виде конъюнкции этих функций, для любой внутренней точки области определяется значением функции фр^), которое является наименьшим среди всех других ^-функций фg (X). Кроме того, для всех точек, определяющих
линейную гиперповерхность /" области О", вычисленное значение ^-функции б" есть величина
постоянная [6].Отсюда следует, что функция О принципиально может быть использована в качестве целевой функции при параметрическом синтезе СТС по критерию запаса работоспособности. Недостатком такой целевой функции является невозможность использования поисковых методов оптимизации, поскольку для произвольной точки, принадлежащей границе области работоспособности, функция О" не является постоянной, а принимает значение из множества возможных
значений {у/1, /2" /2"т,/}.
Построим ^-функцию, которая будет принимать единственное значение для любой точки, находящейся на одинаковом расстоянии от границы области работоспособности, т.е. от любой гиперповерхности/. С этой целью на основании (6) сформируем ^-функции вида
<pf (X) = fj (X)/|grad f j (X )|.
Легко увидеть, что полученные ^-функции обладают требуемым свойством, причем для любой внутренней точки области работоспособности, находящейся на одинаковом расстоянии от ее ближайшей граничной точки, вычисленное значение функции будет равно /. При этом искомая функция будет иметь следующий вид:
Gf= 0,5 (Mf + Pf-
Mf + if
где pf = P, Mf = 2УЧ1'{0,5(ф£у(X) + ф£у+1(Х)-|ф^у(X)-qjf,y+1(X)|) = о}, ф^) = фГ(X).
Таким образом, функция Gf может являться целевой функцией при оптимизации СТС на максимум запаса работоспособности, причем для вычисления координат оптимальной точки по критерию max Gf может быть использован любой поисковый метод оптимизации. Кроме того, при поиске исключается зацикливание алгоритма в независимости от формы границы области ра-
ботоспособности. Важным свойством полученной функции является возможность распознавания состояния СТС. Если вычисленное значение функции положительное, то система находится в работоспособном состоянии. Если результат окажется отрицательным, то система находится в неработоспособном состоянии. В том случае, если значения внутренних параметров выражены в относительных единицах, вычисленное в любой внутренней точке области работоспособности значение функции будет характеризовать относительное значение запаса работоспособности СТС, принадлежащее интервалу [-1; 1].
При использовании для оптимизации градиентных методов, характеризующихся наибольшим быстродействием, при построении функции Gf следует использовать формулу (3), которая позволяет осуществлять операции дифференцирования.
Вывод
Выбор оптимальных значений внутренних параметров СТС по критерию максимального или заданного запаса работоспособности обеспечивает работоспособное состояние системы на предстоящий период времени. Это особенно актуально для СТС водного транспорта, большинство из которых характеризуются параметрической нестабильностью. В том случае, если известна информация о границе области работоспособности, заданная в виде системы линейных ограничений на значения ее внутренних параметров, возможно построение целевой функции, обеспечивающей поиск оптимума на основе известных алгоритмов. При этом исключается зацикливание в процессе поиска, а полученный результат в относительных единицах характеризует запас работоспособности системы. Аналитическое описание области работоспособности на основе использования логических ^-функций позволяет достаточно просто идентифицировать текущее состояние СТС и решать задачи прогнозирования. Рассмотренный метод был апробирован при решении задач параметрического синтеза электротехнических систем и устройств объектов водного транспорта.
Список литературы
1. Саушев, А. В. Основы управления состоянием электротехнических систем объектов водного транспорта / А. В. Саушев. - СПб. : ГУМРФ им. адмирала С. О. Макарова, 2015. - 222 с.
2. Саушев, А. В. Структура процесса управления состоянием сложных электротехнических систем / А. В. Саушев // Надежность и качество сложных систем. - 2013. - № 3. - С. 23-30.
3. Абрамов, О. В. Параметрический синтез стохастических систем с учетом требований надежности / О. В. Абрамов. - М. : Наука, 1992. - 176 с.
4. Норенков, И. П. Основы автоматизированного проектирования / И. П. Норенков. - М. : МГТУ им. Н. Э. Баумана, 2002. - 336 с.
5. Саушев, А. В. Области работоспособности электротехнических систем / А. В. Саушев. - СПб. : Политехника, 2013. - 412 с.
6. Саушев, А. В. Параметрический синтез электротехнических устройств и систем / А. В. Саушев. - СПб. : ГУМРФ им. адмирала С. О. Макарова, 2013. - 315 с.
7. Yurkov, N. K. Synthesis of a Conceptual Model of a Subject Domain. Characteristic Features of Modeling Complex Systems / N. K. Yurkov // Measurement Techniques. - 2004. -Vol. 46, № 2. - February. - P. 128-133.
8. Баранов, Н. А. Управление состоянием готовности системы безопасности к отражению угрозы / Н. А. Баранов, Н. А. Северцев // Труды Междунар. симп. Надежность и качество. - 2012. - Т. 1. - С. 8-10.
9. Лысенко, А. В. Способ снижения величины вибрационных нагрузок в несущих конструкциях ЭС и методика его реализующая / А. В. Лысенко // Надежность и качество сложных систем. - 2013. - № 4. -С. 41-44.
10. Дедков, В. К. Компьютерное моделирование характеристик надежности нестареющих восстанавливаемых объектов / В. К. Дедков, Н. А. Северцев // Труды Междунар. симп. Надежность и качество. - 2010. -Т. 1. - С. 368-370.
11. Yurkov, N. K. Characteristic Features of the Control of Complex Systems Utilizing Conceptual Models / N. K. Yurkov // Measurement Techniques. - 2004. - Vol. 47, № 4. - April. - P. 339-342.
12. Саушев, А. В. Параметрический синтез технических систем на основе линейной аппроксимации области работоспособности / А. В. Саушев // Автометрия. - 2013. - Т. 49, № 1. - С. 61-67.
13. Саушев, А. В. Планирование эксперимента в электротехнике / А. В. Саушев. - СПб. : СПГУВК, 2012. -272 с.
14. Саушев, А. В. Структура, метод и алгоритмы оптимального параметрического синтеза динамических систем / А. В. Саушев // Труды Междунар. симп. Надежность и качество. - 2013. - Т. 1. - С. 214-217.
15. Саушев, А. В. Метод синтеза многопараметрических динамических систем на основе информации о границе области работоспособности / А. В. Саушев // Труды Междунар. симп. Надежность и качество. -2014. - Т. 1. - С. 120 - 123.
Саушев Александр Васильевич кандидат технических наук, профессор, кафедра электропривода и электрооборудования береговых установок, Государственный университет морского и речного флота им. адмирала С. О. Макарова (198035, Россия, г. Санкт-Петербург, ул. Двинская, 5/7) E-mail: [email protected]
Аннотация. Рассматривается метод определения оптимальных значений внутренних параметров технических систем, основанный на аналитическом описании области работоспособности, заданной в виде совокупности линейных ограничений. Формируется выражение для целевой функции, которая позволяет использовать любой известный поисковый метод оптимизации системы по критерию запаса работоспособности. Основой для построения целевой функции является аппарат логических Л-функций.
Ключевые слова: параметрический синтез, запас работоспособности, Л-функции.
Saushev Aleksandr Vasil'evich
candidate of technical sciences, professor,
sub-department of electric drive
and electric equipment of coastal installations,
State University of sea and river fleet
named after admiral S. O. Makarov
(198035, 5/7 Dvinskaya street, Saint-Petersburg, Russia)
Abstract. The method of determination of optimum values of internal parameters of technical systems based on the analytical description of area of the working capacity set in the form of set of linear restrictions is considered. Expression for criterion function which allows to use any known search method of optimization of system by criterion of a stock of working capacity is formed. A basis for creation of criterion function is the device of logical R-functionsand real electro-technical systems at the solution of a problem of their parametrical optimization.
Key words: parametrical synthesis, stock of working capacity, R-function.
УДК 621. 396
Саушев, А. В.
К проблеме синтеза целевой функции параметрической оптимизации сложных технических систем / А. В. Саушев // Надежность и качество сложных систем. - 2015. - № 3 (11). - С. 3-9.
