Глава 1. Показатели технического риска
1.1. Человеческий фактор
Во многих динамических системах человек-оператор был и остается основным звеном в системе контроля и управления. В наибольшей степени это касается таких систем, как «самолет – летчик». Задача оператора-летчика состоит в построении цели и выборе средств ее достижения путем преодоления внешних возмущений. При этом роль человека (летчика) состоит в адекватном отображении полетной ситуации с помощью систем отображения информации, формирования концептуальной модели режима полета, которая создается в его сознании путем анализа информационных потоков и знания моделей реакции ЛА на его управляющие воздействия [8].
На процесс управления, формируемый оператором, сказывается так называемый человеческий фактор, включающий в себя психические, умственные, физиологические, биомеханические, антропометрические и другие свойства человека. Под влиянием различных, в том числе человеческих, факторов могут возникать полетные режимы: штатный, особый, аварийный, критический. Штатный режим возникает, когда параметры xi управляемой динамической системы, подлежащие контролю и управлению, находятся в заданной или расчетной области, которую называют допустимой. Множество всех значений x1 из этой области обозначают Ωдоп. Все остальные ситуации характеризуются выходом параметров xl из области Ωдоп и называются особыми, в частности, аварийными или критическими. Возникновение последних приводит к летным происшествиям, в результате которых происходит полное или частичное разрушение ЛА с возможной гибелью пассажиров или членов экипажа.
Отметим, что в словаре Вебстера слово «риск» определяется как «опасность, возможность убытка или ущерба». В этом смысле особые ситуации, возникающие при полете ЛА, сопряжены с риском.
По оценкам специалистов особые ситуации в процессе эксплуатации ЛА на 20 % обусловлены техническими, а на 80 % – человеческими факторами [8]. Эти обстоятельства приводят к мысли о необходимости устранения человеческого фактора из контроля и управления ЛА. Однако устранить этот фактор невозможно, так как все технические системы любого уровня созданы человеком, и человеческий фактор всегда будет присутствовать в любой технической системе либо на уровне программ управления, либо на уровне технико-технологического совершенства.
Создавая автоматы, позволяющие исключить человека из контура контроля и управления ЛА, мы заменяем одни потери на другие, одни ошибки – другими, и величина потерь должна быть оценена по частоте их возникновения путем сравнения с выигрышем, выгодой. Таким образом, возникает оптимальная задача распределения функций контроля и управления между летчиком и автоматом, решение которой обеспечит максимальные выгоды.
Полет ЛА осуществляется, как правило, с помощью двух информационно-аналитических центров: ИАЦ-1 и ИАЦ-2. Информационно-аналитический центр человека (ИАЦ-1) рассмотрен в работе [8], где указаны следующие ему присущие ограничения: по объему оперативной и долговременной памяти; по скорости обработки информации; по точности обработки информации; по наличию зоны нечувствительности.
Под ИАЦ-2 будем понимать совокупность технических систем, подсистем, блоков, элементов, осуществляющих сбор и обработку информации с целью формирования управлений для достижения целей функционирования ЛА. До недавнего времени бортовое оборудование в авиации, как правило, не включало информационно-аналитический центр (ИАЦ-2).
ИАЦ-1 и ИАЦ-2 могут дополнять друг друга, в полете, тем самым, снижая стоимость бортового оборудования, расширяя область применения возможностей ЛА, что повышает выгоду от его использования. Что касается формирования управления, то ИАЦ-1 – система с заданными свойствами, которая обладает указанными выше ограничениями.
Параметры технической системы ИАЦ-2 можно выбирать и, причем оптимальным образом с учетом ограничений на ее возможности. Эту систему, в отличии от ИАЦ-1, мы можем создавать с нуля, в то время как ИАЦ-1 мы вынуждены воспринимать в основном уже созданной. Мы получаем два антипода, дополняющих друг друга. При этом параметры ИАЦ-1 можно изменять в достаточно широком диапазоне путем тренировок и обучения, однако основные физиологические параметры (быстродействие умственной деятельности, объем памяти и т. п.) мы не в состоянии существенно изменить.
Объединяя ИАЦ-1 и ИАЦ-2 в единый комплекс, мы получаем новый ИАЦ, который лишен недостатков ИАЦ-1 и ИАЦ-2, т. е. каждого в отдельности. Проектирование совместного комплекса связано с определенными трудностями, так как требует
– подбора (в процессе обучения по специальным программам) экипажа с наилучшими возможностями с позиции достижения минимального технического риска и со свойственными ему физиологическими свойствами и ограничениями, параметры которых задают некоторую область Ω1;
– проектирования самолета и двигателя как системы, свойства и параметры которых заполняют некоторую область Ω2, заданную инвестором;
– проектирования бортового оборудования под этот комплекс, способного обеспечить надежное (устойчивое, безопасное, оптимальное) состояние и достижение цели полетного задания в области Ω2, т. е. при расширении области Ω1 до Ω2.
При этом возможны различные ситуации, связанные с совместным проектированием комплекса с позиции минимизации потерь при эксплуатации, начиная от измерения параметров самолета и двигателя до разработки специальных для данного класса самолетов комплексов тренажеров, позволяющих осуществлять необходимую подготовку летчиков (доводку его ИАЦ-1 до необходимого состояния). По-существу, здесь рассматривается идея проектирования с максимальным использованием самого дешевого продукта для эксплуатации ЛА – человека.
В работе [8] рассмотрен частный подход к решению задачи выбора параметров комплекса (рис. 1.1) с использованием одного критерия, связанного с энтропией системы. В общем случае необходимо использовать тот критерий, который связан с выполнением цели и назначения данного класса ЛА. Так, например, при проектировании комплекса следует учесть необходимость исключения летчика из контура управления, когда на вход ИАЦ поступают процессы высокой скорости изменения по амплитуде и частоте. С этой целью необходим анализатор входных сигналов xi
– по амплитуде и удаленности его от Ωдоп;
– по скорости изменения xi;
– по количеству сигналов, которые могут достигнуть Ωдоп.
При этом необходимо осуществлять прогноз достижимости критической области [6].
Рис. 1.1
Несмотря на эти ограничения, рассмотренный подход может быть использован в качестве метода или методики в случае, если критерий уточнен, расширен в зависимости от целей и задач проектируемого ЛА, когда необходимо решать многокритериальные задачи.
1.2. Прибыль и убытки в технике
Прибыль и убытки в авиации, как и везде в технике, экономике, взаимосвязаны, взаимозависимы. Это антиподы, не существующие друг без друга. Таков основной закон среды жизнедеятельности. Как указано в работе [6]: риск – это ситуативная характеристика, представляющая собой величину потерь и частоту их появления, которая связана функционально с прибылью, полученной при этом.
В необходимости учета риска при разработке проекта (создания ЛА и его систем или организации эксплуатационного предприятия) заинтересованы следующие его участники: заказчик, инвестор, исполнитель, страховая компания. При анализе риска любого из участников проекта используются положения, предложенные американским экспертом Б.Берлимером:
– потери от рисков независимы друг от друга;
– потеря по одному направлению «портфеля рисков» не обязательно увеличивает вероятность потери по другому (за исключением форс-мажорных обстоятельств);
– возможный максимальный ущерб не должен превышать финансовые и другие возможности участника проекта.
При анализе характеристик риска выделим два взаимно дополняющих друг друга вида: количественный и качественный. Качественный анализ может быть сравнительно простым, его главная задача определить факторы, влияющие на риск, этапы и работы, при выполнении которых риск возникает. Количественный анализ сводится к численному расчету размеров отдельных компонент риска и риска проекта в целом. Этой проблеме посвящена данная работа.
Все факторы, так или иначе влияющие на рост величины риска в проекте, можно условно разделить на две группы: объективные и субъективные.
К объективным относятся факторы, независящие непосредственно от самой фирмы или авиационного комплекса: это инфляция (старение техники), анархия, политические и экономические кризисы, экология, таможенные пошлины, наличие режима наибольшего благоприятствования.
К субъективным относятся факторы, характеризующие непосредственно данную фирму, данный проект, данный авиационный комплекс; производственный потенциал, техническое оснащение, уровень предметной и технологической специализации, организация труда, уровень производительности.
При разработке нового ЛА или доработке старой модификации (путем установки нового бортового оборудования) возникают как взаимный интерес инвестора и конструкторского бюро, так и противоречия. Задача заказчика состоит в том, чтобы при минимальных затратах создать такой ЛА, который по основным показателям превысил бы известные ЛА. Задача конструкторского бюро в том, чтобы найти возможность удовлетворить требования заказчика. Как правило, не удается полностью достичь того, что хочет заказчик на те средства, которые он выделил. При этом эксплуатационники, а также пользователи услуг, страховые компании, организации типа IKAO требуют ЛА с заданной надежностью.
В качестве основного показателя, предъявляемого к ЛА, являются экономические показатели. Все остальные порождены этим показателем, за исключением показателя, который связан с человеческими жертвами. Так, например, такой показатель, как «регулярность» обеспечивает заданную величину отложенных полетов, учитывая их высокую стоимость. Показатель «безопасность» связан с расходами на поломку или восстановление техники, а также со страховыми выплатами. При этом по-существу из одного показателя экономичности был введен векторный показатель: экономичность, безопасность, регулярность. В последнее время к этим показателям добавился расход топлива.
Экономичность включает в себя чрезвычайно важную проблему авиастроения – оптимальное распределение стоимости производства конструкции ЛА и бортового оборудования, что обеспечивает получение максимальной прибыли при минимальных потерях.
При этом прибыль связана только с этапом эксплуатации, а потери с этапами разработки, производства, эксплуатации. Прибыль также существенно зависит от качества исполнения несущих аэродинамических поверхностей, в том числе от их свойств и возможностей, и от погрешностей δХ функционирования бортового оборудования. Чем выше возможность и меньше δХ, тем дороже самолет, тем большие издержки несет инвестор, которые могут не окупиться.
Введем общее расчетное (максимальное) количество полетов N, которые может совершить самолет за время Т. Пусть из-за погодных условий он не сможет совершить (при его низком показателе регулярности) n1 полетов. Из-за аварийных ситуаций он не завершит n2 полетов (в том числе поломок), а из-за недостоверной информации как бортового оборудования, так и средств управления воздушным движением – n3 полетов. Таким образом, полеты n0 = N – (n1 + n2 + n3) = N – n4 завершены благополучно и могут принести прибыль, а полеты n4 = n1 + n2 + n3 принесут убытки. Каждая из составляющих вектора = {n0, n1, n2, n4} несет в себе определенную информацию с позиции функционирования бортового оборудования:
n0 – выполнение поставленной цели;
n1 – невыполнение поставленной цели, при правильном функционировании систем контроля и управления бортовым оборудованием;
n2 – возникновение аварийных ситуаций, включая катастрофы, обусловленные превышением критических значений параметров состояния ЛА из-за погрешностей δx функционирования систем контроля бортового оборудования;
n3 – недостижение поставленной цели, в том числе отказ ее достижения из-за ложной информации систем контроля и средств управления воздушным движением.
Таким образом, убытки, а, следовательно, технический риск, обуславливают те же события, которые связаны с {n2, n3, n4}, из них {n3, n4} обусловлены погрешностями δx получения и обработки информации.
Качество системы управления и изготовления ЛА будем характеризовать ошибками выдерживания траектории Δx и, в частности дисперсией D(Δx) = σ21(Δx).
Если в качестве цели ставиться полет на дальность L, то в этом случае потери будем характеризовать частотой n2 – невыполнение полетов (в том числе по погодным условиям); безопасность характеризовать частотой n3; экономичность будем оценивать совокупностью {n2, n3, n4}, характеризующих потери в процессе эксплуатации; а с помощью n1 – прибыль, связанную с благополучным выполнением поставленной цели.
Для современной авиации характерны наперед заданные ограничения на компоненты потерь, связанные с безопасностью полетов, т. е. по-существу задана плата за риск эксплуатации и связанная с ним прибыль. Так, на посадке суммарный риск не должен превышать Р = 10–9. При этом предполагается, что современная авиация с современным оборудованием гарантированно имеет потери (убытки). Если абстрагироваться от реальности, то можно добиться от бортового оборудования такого функционирования, при котором нет катастроф, но стоимость такого самолета будет так высока, что доходы за счет n1 могут не покрыть эти расходы.
В общем случае количественные характеристики риска представляют векторные величины, а задача построения и прогноза их чрезвычайно сложна [6]. Таким образом, приступая к проектированию самолета, мы должны учитывать:
– затраты на создание и эксплуатацию;
– прибыль при эксплуатации;
– потери в процессе создания и эксплуатации.
В качестве примера рассмотрим техническую постановку задачи создания новых образцов авиационной техники.
Одной из основных задач, стоящих перед проектировщиками и разработчиками таких сложных и дорогостоящих технических систем, как авиационный комплекс, включающий: самолет и его бортовое оборудование; системы управления воздушным движением; аэродромные средства, является задача выбора и обоснования технических требований к комплексу, в которых отражалось бы целевое назначение его и которые соответствовали бы научно-техническому потенциалу разработчиков. При проектировании авиационного комплекса выбор технических требований к нему должен производиться исходя из целей и задач, стоящих перед проектировщиком самолета, в том числе и его бортовым оборудованием. Такие цели формулируются, как правило, на качественном уровне и позволяют судить лишь об общем направлении работ по созданию авиационного комплекса и его совершенствованию. Для обеспечения необходимой ясности и однозначности формулировок целей последние лучше задать в терминах характеристик авиационного комплекса. Для этого генеральную цель – выполнение самолетом полетного задания – приходится разбивать на совокупность более частных, зато более простых и конкретных подцелей, то есть проводить квантификацию целей. Такими подцелями являются обеспечение регулярности полета, его безопасности и экономичности.
Осуществив квантификацию, получают многоуровневое иерархическое дерево целей, на нижнем уровне которого оказывается полный набор измеримых целей. Для обеспечения полноты в набор целей нижнего уровня приходится включать цели, характеризующие различные стороны процесса функционирования системы. Дерево целей позволяет иметь полный перечень задач для подсистем любого уровня иерархии. При этом проектировщик формулирует технические задания для подсистем, в которых учитываются цели и задачи более высокого уровня.
Процесс квантификации целей завершен, когда получен набор количественно измеримых подцелей, связанных с показателями эффективности функционирования подсистем и системы, т. е. авиационного комплекса в целом. На практике обычно используется следующий, чисто эмпирический подход к построению показателей эффективности и оценке качества систем.
Из множества технических показателей систем авиационного комплекса лицо, принимающее решение, выделяет тот или те, которые, по его мнению, в наибольшей степени характеризуют соответствие системы заданному целевому назначению. Поскольку авиационный комплекс служит для обеспечения регулярности (R), безопасности (Б) и экономичности (Э) полета самолета, последние являются показателем эффективности авиационного комплекса. Отсюда следует, что задача проектирования авиационного комплекса заключается в том, чтобы создать такой авиационный комплекс, который обеспечивал бы самолету значения показателей регулярности, безопасности и экономичности его полета не хуже существующих, и при этом обеспечивал бы прибыль.
Таким образом, целью нового авиационного комплекса или совершенствования старого является, как следует из вышеизложенного, повышение регулярности, безопасности и экономичности полетов самолета. Как правило, реализация этой цели поддается экономической оценке, в результате чего могут быть получены зависимости
J1 = J1(ΔR, ΔБ, ΔЭ, Т), J2 = J2(ΔR, ΔБ, ΔЭ, Т),
где J1 – прибыль за время эксплуатации самолета, оснащенного таким авиационным комплексом; ΔR, ΔБ, ΔЭ – соответственно приращения показателей регулярности, безопасности и экономичности полета нового самолета по отношению к аналогичным показателям старого варианта самолета; J2 – затраты на создание авиационного комплекса; Т – время эксплуатации. Очевидно, что эффект от внедрения
ΔJ = J1 – J2. (11)
Рассмотрим вектор А параметров, полностью характеризующих авиационный комплекс. Тогда R = R(A); Б = Б(А); Э = Э(А), задача заключается в отыскании такого А = А*, при котором показатель (1.1) достигает максимальной величины на множестве значений ΔJ, на границах которого значения ΔJ достигают порога, характеризующего целесообразность создания авиационного комплекса. Таким образом, задача состоит в отыскании А = А*, удовлетворяющего условию
А* = min{J1(R(А), Б(А), Э(А)) – J2(R(А), Б(А), Э(А))}. (1.2)
В результате процесс проектирования авиационного комплекса сводится к построению алгоритма, с помощью которого устанавливается связь между свойствами вектора А* параметров авиационного комплекса и значениями R, Б, Э, а также метода нахождения А*, удовлетворяющего условию (1.2).
Предположим, что показатели регулярности, безопасности и экономичности полета представляют собой вероятности возникновения некоторых событий (например, особых ситуаций, опасных ситуаций, ложных срабатываний [6]). Предположим также, что алгоритм (метод) расчета эффекта J1 в зависимости от значений, указанных показателей известен. В качестве примера такого алгоритма рассмотрим алгоритм, устанавливающий зависимость между эффектом J1 и значениями показателя безопасности полета, под которым будем понимать вероятность или частоту особых ситуаций.
Пусть для всего парка самолетов заданного класса известно общее количество особых ситуаций, имевших место за заданный период времени. Это позволит определить экономические потери П*, обусловленные такими ситуациями. С другой стороны, предположим, что в результате проектирования будет создан такой авиационный комплекс, который обеспечит уменьшение особых ситуаций за тот же период времени, в результате чего потери от них составят величину П**. Тогда экономический эффект P1 от эксплуатации самолета-носителя, имеющего такой авиационный комплекс выразится следующим образом:
Р1 = П* – П**.
С учетом введенных предположений определение вектора А*, характеризующего авиационный комплекса и удовлетворяющего выражению (1.2), сведется к задаче определения затрат J2 на создание комплекса, обеспечивающего самолету значения показателей R, Б, Э полета не хуже заданных (требуемых).
Решение данной задачи может быть сведено к последовательному решению следующих двух задач: задачи синтеза структуры авиационного комплекса, обеспечивающего значения указанным показателям не хуже требуемых, и задачи определения затрат на создание авиационного комплекса, имеющего такую структуру.
Рис. 1.2
В простейшем случае процесс создания нового ЛА или совершенствования старого связан с инвестором (рис. 1.2). Как правило, инвестор, стремясь получить максимальный доход, заказывает КБ проектирование нового самолета с характеристиками R, Б, Э. Назовем их условно характеристиками идеального самолета. В силу ограниченных возможностей КБ создает вариант самолета с характеристиками (R1, Б1, Э1); назовем его проектный вариант ЛА. На последней стадии создания, на стадии производства все изменяется, и мы получаем ЛА с характеристиками (R2, Б2, Э2). Назовем такой самолет реальным или фактическим. В результате идеальный доход (Dн), на который рассчитывал инвестор не получился и стал равен некоторому фактическому значению Dф. В случае, если расхождение ΔD = Dн – Dф велико, в КБ проводятся исследования, направленные на поиск наилучшего соответствия между ΔD и стоимостью оборудования, необходимого для его уменьшения (рис. 1.3).
Представленная на рис. 1.3 структурная модель поиска наилучшего решения включает финансовые вложения (финансовые потоки) и порожденные ими технические совершенства, которые обеспечивают регулирование (изменение) дохода D(t) и прибыль Dп(t) в нужную для заказчика сторону. При этом новая техника внедряется, если происходит увеличение n0/n, т. е. частоты выполненных полетов, при уменьшении n1/n, n2/n, n3/n.
Анализируя все сказанное относительно потерь и прибыли, сопутствующим процессу эксплуатации самолета, в том числе анализируя процессы доработки, эксплуатации и создания самолета, можно сделать вывод, что в самом простейшем случае самолет есть экономическая подсистема с переменными параметрами, подверженная старению, износу, поломке, изменению точностных и надежностных характеристик.
Рис. 1.3
При этом на такую экономическую систему воздействуют внешние и внутренние возмущающие факторы, изменяющие, как правило, в сторону уменьшения, величину дохода от эксплуатации самолета (рис. 1.4). В такой системе, созданной по воле человека и непосредственно его трудом, циркулируют деньги и техника, взаимно преобразовываясь в нужном для человека направлении. Участие человека сводится к выбору параметров системы, например, нелинейной обратной связи, выбора типа и параметров новой техники, а также места и времени применения ЛА.
Рис. 1.4
Так, нелинейная обратная связь (рис. 1.5) имеет зону нечувствительности ε, в которой ΔD = 0, т. е. отчисления ΔD на создание новой техники с целью компенсации потерь, равны нулю. Рассматриваемая экономическая подсистема имеет внешнюю обратную, направленную на расширение бизнеса. В общем случае каждый ЛА при анализе эксплуатационных потерь можно рассматривать как подсистему (объект, элемент) некоторой экономической системы, представляющей собой эксплуатационное подразделение. Это означает, что все основные положения экономических систем могут быть использованы при анализе процессов создания и эксплуатации ЛА.
Рис. 1.5
В зависимости от типа самолета в качестве таких экономических систем выступают эксплуатационные подразделения следующих уровней:
– частное лицо (1–2-х местные самолеты);
– спортивное общество (спортивная авиация);
– фирма (самолеты деловых связей);
– региональные авиаслужбы (самолеты местных авиалиний);
– авиаслужбы страны (самолеты внутренних авиалиний);
– международные авиакомпании (самолеты международных авиалиний).
Каждое из этих эксплуатационных подразделений, в той или иной мере оказывает влияние на создание соответствующей (эксплуатируемой им) техники, а также средств обеспечения ее функционирования (рис. 1.6). Заказчики или инвесторы техники имеют разные итоговые цели ее создания, разные условия функционирования и потому формируют различные требования и назначения. При этом для каждого ЛА характерны свои потери как по величине, так и по виду своего проявления, что обуславливает необходимость анализа потоков технического риска (потерь).
Рис. 1.6
1.3. Истоки технического риска
Технико-технологическая среда, являющаяся одним из основных элементов среды жизнедеятельности человека [6], включает в себя технические и технологические объекты, созданные с использованием научных достижений. Огромная польза науки заключается в ее способности создавать необходимые предпосылки для обеспечения безошибочных действий человека, связанных с созданием новой техники или совершенствованием старой.
На рис. 1.7 показана взаимосвязь: среды жизнедеятельности; науки; технических и технологических объектов. В основе такого взаимодействия лежит обмен информацией Jx от среды к науке Jy, от науки Jy к технике Jz о процессах x, y и z, формируемых средой, наукой и технологией соответственно. В процессе такого обмена возникают погрешности и ошибки, свойственные контролю, обработке и передаче информации. Так, в процессе изучения среды жизнедеятельности, которая является источником информации и базовой основой для построения математических и иных моделей, возникают погрешности δ1, которые зависят от средств измерения и свойств среды жизнедеятельности как носителя информации Jx.
Численные характеристики погрешности δ1 в процессе научных исследований видоизменяются, как правило, увеличиваются, и на входе в подсистему (3) приобретают новое значение δ2. В процессе создания новой техники или технологий δ2 видоизменяется, превращаясь, с одной стороны, в случайные факторы W31(δ1, δ2, δ3, y, z), которые воздействуют на среду жизнедеятельности, а с другой – в δ3, которые наблюдает наука с целью корректировки математических или иных моделей. Среда жизнедеятельности реагирует на W31 и создает на входе (3) случайные факторы W13(δ1, x, t).
Рис. 1.7
В процессе взаимодействия между подсистемами (2) и (3) возникают противоречия, обусловленные трудностями, а подчас и невозможностью выполнения всех требований, сформулированных наукой, подлежащих реализации на этапе создания технико-технологических объектов. Это обусловлено, чаще всего, отсутствием соответствующей производственной базы, например, для обеспечения заданной точности и надежности функционирования самолета, или заданной точности конструктивного исполнения несущих аэродинамических поверхностей.
При этом наука задает расчетную или потребную точность, например, в виде погрешности δ2, а на практике получают возможную погрешность δ*31, которая, как правило, превышает δ31. В результате ухудшаются все показатели, например, тактико-технические данные самолета, такие, как дальность полета, расход топлива, грузоподъемность, ресурс, категория по взлету и посадке.
Ограничимся основными требованиями к технике, которые сформулируем, следующим образом:
1) техника должна оптимальным или наилучшим образом реализовать свое целевое назначение, так например, перемещаться в пространстве с минимальными энергетическими затратами при прочих равных условиях;
2) функционирование техники должно быть безопасным в целом, для биосферы и, в частности, для тех, кто ее использует, так, например, вероятность катастрофы для самолета не должна превышать 10–8.
Невыполнение указанных требований приводит к соответствующим неучтенным или сверхнормативным потерям, которые включают, в частности, потери финансовых средств, направленных на создание новой техники.
Процессу создания и эксплуатации самолета сопутствуют следующие потери:
R1 – от поломки или разрушения самолета;
R2 – топлива, из-за несовершенства бортовой аппаратуры;
R3 – топлива, из-за несовершенства аэродинамическо-прочностных свойств самолета;
R4 – эксплуатационные, из-за несовершенства организации работ по подготовке к эксплуатации;
R5 – на этапе научно-исследовательских работ;
R6 – на этапе опытно-конструкторских работ.
Таким образом, потери R при создании и эксплуатации представляют собой векторную величину, включающую в себя n-компонен-тов, каждая из которых также является, векторной величиной.
Потери R связаны с жизненным циклом, который включает в себя следующие этапы (рис. 1.8): научно-исследовательский, опытно-конструкторский, серийное производство, эксплуатационный. Согласно схеме рис. 1.8, инициатором создания того или иного проекта является инвестор-1, на средства которого в объеме (Z1 + Z2) организуются НИР (научно-исследовательские работы), ОКР (опытно-конструкторские работы), производство. Реализовав объект на рынке-2 по цене Z3, он должен получить прибыль П1 = [Z3 – (Z1 + Z2)] > 0.
Инвестор-2, купив объект по цене Z4, организует его эксплуатацию, представив его на рынок-3 (услуг) для потребителя услуг (10). В процессе эксплуатации ивестор-2 от потребителя получает Z5 – стоимость услуг. При известных эксплуатационных расходах, ресурсе планера инвестор-2 имеет возможность определить прибыль.
Рис. 1.8
Перечень расходов, сопутствующих жизненному циклу, включает в себя: 1) создание самолета, в том числе научно-исследовательские работы, проектирование, производство и испытание опытных образцов; 2) серийное производство; 3) эксплуатационные расходы, в том числе организация управления воздушным движением (управление воздушным движением), аэродромное и техническое обслуживание. Все эти расходы должны окупаться за счет коммерческих перевозок за период эксплуатации, в качестве которого может выступать ресурс работы планера самолета. При этом Z5 есть функция таких факторов, как стоимость тонны километра перевозимого груза, которая изменяется случайным образом во времени; количество километров, налетанных самолетом и др.
Потери инвестора представляют собой превышение расходов над доходами за месяц, год или за весь период эксплуатации данного самолета или всех самолетов данной серии. Каждый из этапов жизненного цикла характеризуется определенными потерями, которые представляют собой этапный риск (рис. 1.9). Потери на каждом из этапов работ зависят от величин погрешностей δiR , допущенных при проведении работ, а также от совокупности средств Аi , с использованием которых проводились работы. В итоге получаем суммарные потери δ4R, которые можно представить в виде
δ4R = ψ(δ1R, δ2R, δ3R, А1 A2, A3),
где ψ – оператор преобразования.
Рис. 1.9
1. На этапе научно-исследовательских работ, включающих: аэродинамические, прочностные расчеты; расчеты по динамике полета; продувки модели самолета и прочие расчеты, потери инвестора обусловлены невозможностью достичь заданную цель. Это приводит к потерям тех финансовых средств, которые были затрачены инвестором на проведение таких работ (рис. 1.9). Обозначим их R(1)1H.
Другой крайностью является ситуация, в которой результаты научно-исследовательской работы показали возможность достижения поставленной цели, а этап ОКР их не подтвердил (R(1)2H). Между этими крайними случаями находится проект, позволяющий достичь заданную цель, но который был отклонен (R(1)3H).
Таким образом, научно-исследовательский риск характеризуется тремя ситуациями, возникающими в процессе проведения работ и которые могут характеризоваться потерями R(1)1H, R(1)2H, R(1)3H.
2. На этапе опытно-конструкторских работ (ОКР), например, для самолета, включающего: проектирование и изготовление опытного образца; проведение аэродинамических, прочностных и летных испытаний, возможны те же ситуации, что и на этапе НИР. Однако потери, а, следовательно, и численная величина риска возрастают за счет более высокой стоимости ОКР. Причиной таких ситуаций являются погрешности δ1R, полученные и не обнаруженные на этапе НИР; а также погрешности ОКР – δ2R.
3. На этапе серийного производства показатели риска увеличиваются по следующим причинам:
– ухудшение показателей объекта, за счет влияния несовершенств технологических процессов производства, обусловленных свойствами металла, станков и инструментов, квалификацией специалистов и т. п.; в результате получаем погрешность δ3R.
– повышение стоимости производства объекта по отношению к заявленной стоимости, что увеличивает численную величину риска, связанную с финансовыми расходами.
Отметим, что изменение характеристик объекта за счет технологических процессов сказывается и учитывается на этапе эксплуатации.
4. Последний этап – эксплуатационный – характеризуется соответствующим риском, связанным, прежде всего, с полной или частичной потерей техники при авиационных происшествиях, а также с фактическими (финансовыми) расходами для обеспечения полета, которые превышают расчетные или оптимальные, например, за счет неоптимального расхода топлива.
Анализ материалов, полученных в процессе эксплуатации самолетов нового поколения, позволяет выявить тенденции изменения причин тяжелых авиационных происшествий в мировой гражданской авиации (рис. 1.10). Из этого рисунка видно, что основная доля авиационных происшествий приходится на этап захода на посадку и посадку (76 %). Внедрение систем автоматического управления и изменение характера летной работы привели к появлению значительной доли авиационных происшествий по причине столкновений с землей на исправном самолете (41 %) и потери способности управлять самолетом в полете (30 %).
Количество происшествий на самолетах СГТ Д-217
1 – заход на посадку и посадка;
2 – столкновение с землей на исправном самолете;
3 – потеря способности управлять самолетом в полете;
4 – прекращение взлета;
5 – полная выработка топлива;
6 – наземные происшествия.
Рис. 1.10. Причины тяжелых авиационных происшествий в мировой гражданской авиации за период с 1989 по 1998 гг. (по материалам Boeing Operations Symposium, 1999)
В следующей таблице приведены основные типы ошибок экипажей приведших к авиационным происшествиям [83].
Типы ошибок экипажей
По данным таблицы видно, что неправильные решения экипажа – основная причина авиационных происшествий. На основе эксплуатации ЛА относительно причин возникновения и способа нейтрализации риска: можно сформулировать следующее:
– уровень риска в полете на современном ЛА определяется уровнем свойств и состоянием всего авиационного комплекса, включая: ЛА; экипаж; бортовые системы управления ЛА и обеспечение жизнедеятельности; наземные средства руководства полета;
– необходимый уровень риска обеспечивается высоким уровнем знаний характеристик ЛА и условий полета, в том числе состояния среды, в которой протекает полет;
– при рассмотрении проблем моделирования риска необходимо учитывать характер многих закономерностей, имеющих место в процессе функционирования авиационного комплекса;
– ведущая роль в благоприятном завершении полета принадлежит системе «ЛА – экипаж».
Эксплуатационный риск – векторная величина, т. е. многомерная характеристика, каждая компонента которой, оценивается разными службами, людьми (организующими полет или участвующими в полете) различным образом. При этом эксплуатационный риск включает в себя следующие виды рисков:
– производственные, связанные с исполнением конструкции ЛА, двигателя и бортового оборудования;
– технические, связанные с бортовым оборудованием, которое обеспечивает предотвращение выхода из области допустимых состояний;
– финансовые, обусловленные превышением фактических экс-плутационных расходов над расчетными или оптимальными, в том числе из-за расхода топлива;
– характеристик ЛА и надежности рекомендаций по управлению им;
– возмущающие факторы, в том числе от среды, в которой протекает полет.
Каждая из указанных компонент риска может быть разбита еще на ряд составляющих. Разбиение зависит от поставленной задачи, степени детализации конечных результатов по оценке влияния отдельных систем ЛА на величину риска.
Важной составляющей риска является уровень знаний характеристик ЛА, его систем бортового оборудования, систем наземного комплекса, с помощью которых осуществляется взлет, посадка, полет по эшелонам. Отказы авиационной техники, ошибки летного состава, непредусмотренные расчетами воздействия среды на ЛА, чаще всего, обусловлены влиянием множества случайных факторов. Неполный учет случайных факторов приводит к случайным непредвиденным исходам в полете.
При решении различных практических задач, принятии различных управленческих (проектных) решений можно учитывать все новые и новые группы факторов от самых существенных до самых ничтожных, тем самым, пытаясь уменьшить роль (предотвратить появление) случайного (ничтожного, неосознанного) фактора. Однако полностью исключить влияние таких факторов невозможно, так как человеческие знания относительны, и каждому уровню познаний соответствует свои погрешности. Кроме того, часто глубина наших знаний ограничивается финансовыми возможностями.
Подводя итог сказанному, сформулируем следующее определение: технический риск есть интегральная характеристика потерь, реализованных за весь жизненный цикл самолета или самолетов данного класса.
1.4. Показатели технического риска
Потери топлива и техники в процессе эксплуатации ЛА являются основными потерями для современных ЛА, о чем свидетельствует анализ материалов их эксплуатации. В качестве основных факторов и параметров, обуславливающих указанные потери при эксплуатации (эксплутационный риск) будем рассматривать параметры:
– траектории полета, в том числе высоту, число Маха, скорость;
– влияющие на расход топлива и безопасность: масса, положение центра тяжести, сопротивление и подъемная сила ЛА.
В процессе разработки показателей технического риска для самолета необходимо учитывать все четыре этапа его жизненного цикла. На этапе научно-исследовательских работ учитывается роль науки; на втором этапе – роль конструкторско-проектных работ; на третьем – роль технологических процессов производства самолета; на четвертом – проблему эксплуатации и, прежде всего, обеспечение безопасности и оптимальности полета. При этом величина технического риска есть интегральная характеристика, обусловленная потерями при создании и эксплуатации технических систем и объектов, т. е. на протяжении всего жизненного цикла, включая НИР, ОКР, производство, эксплуатацию.
Система оптимизации режимов пилотирования позволяет вычислить при заданных свойствах планера и двигателя, конфигурации самолета оптимальные параметры траектории полета, при которых достигается максимальная (оптимальная) дальность полета L или минимальный расход топлива q.
В случае, когда параметры траектории полета отличаются от оптимальных, например, при отсутствии системы оптимизации режимов пилотирования (СОРП), и топливо расходуется не оптимально, происходят потери, что обуславливает соответствующий риск. Величина потерь для современного лайнера может достичь больших величин, если его не оборудовать СОРП.
В реальных условиях на ЛА могут быть установлены:
– только система оптимизации режимов пилотирования;
– только система предупреждения критических режимов;
– одновременно обе системы.
Случай, когда имеется только СПКР, рассмотрен в работе [6]. Назначение СОРП состоит в обеспечении в каждом полете максимальной дальности полета. Рассмотрим, как в процессе проектирования самолета, имеющего СОРП, формируется показатель риска, связанный с достижением или не достижением максимальной дальности полета. Заданное по техническому заданию значение дальности полета (L) это то ее значение, которое имеет место быть, когда математические модели, методы расчетов и проектирования не имеют погрешностей. Обозначим эту величину Lкр – критическая, равная максимальному значению, дальность полета ЛА.
При проектировании мы получаем расчетную величину дальности полета Lр. Если бы методы расчета и математические модели не обладали погрешностями, то Lр = Lкр. В реальности за счет по-грешностей расчета δLp имеем Lp = Lкр + δLр, т. е. Lр ≠ Lкр. Равенство имеет место, когда M{δLр} = 0, где M{δLр} – математическое ожидание случайной величины δLр. При этом имеем M{Lр} = Lкр.
При проектировании с целью обеспечения заданной (максимальной) дальности полета выбираются необходимые параметры самолета (обозначим их через А). За счет ошибок δLр мы можем выбрать параметр ЛА с ошибками δА такими, что дальность полета L0 < Lкр. Прогнозирование и оценка такой ситуации важна для инвестора, так как расходы на работы в случае L0 < Lкр могут не окупиться и возрасти после проведения ОКР, когда обнаружится этот факт, например, за счет того, что они должны быть повторены.
Погрешность расчетов δLр представим в виде: δLр = δL1 + δL2, где δL1 – погрешность, обусловленная расчетами аэродинамических сил Rx, Ry, Rz и другими факторами, влияющими на Rx, Ry, Rz; δL2 – погрешности, обусловленные несовершенством алгоритмов оптимизации, построенных на этапе НИР, включающих погрешности, стабилизации и управления. Отметим, что погрешности существующих методов расчета аэродинамических характеристик достигают 10 %.
На этапе опытно-конструкторских работ уточняется Lкр = Lmax и после экспериментального полета получают величину равную . При этом = Lр + δL0, где δL0 – погрешность дальности полета, обусловленная влиянием погрешностей опытных образцов бортового оборудования, двигателя, конструкции самолета.
В процессе производства каждая из подсистем; планер; двигатель; бортовое оборудование изготавливаются с погрешностью δn и в результате для идеальных условий эксплуатации получаем Lп, которая включает в себя погрешность δLп и тогда Lп = Lр + δLп.
При эксплуатации ЛА за счет влияния внешних δв1 и внутренних δв2 возмущающихся факторов получаем эксплуатационную или фактическую дальность полета вида Lф = Lп + δLэ(δв1, δв2). Погрешности δLр, δLп, δLэ являются случайными, каждая из них характеризуются своей плотностью вероятности. При этом фактическое значение Lф совпадает со значением, полученным при эксплуатации. При таких условиях технический риск, имеющий место при проведении ОКР, будет связан с теми ситуациями, которые возникают в процессе выполнения его, а именно
А1 = (Lф < Lкр, Lр ≥ Lдоп), А2 = (Lф ≥ Lкр, Lр ≤ Lдоп),
А3 = (Lф < Lкр, Lр ≤ Lдоп), А4 = (Lф ≥ Lкр, Lр ≥ Lдоп).
В результате искомые вероятности будут включать: Р1 = Р(А1); Р2 = Р(А2); Р3 = Р(А3). Вероятность Р2 – мы отвергли экономически выгодную конструкцию (объект); Р1 – приняли к исполнению вариант ЛА, который не подтверждает по факту возможность достижения поставленной цели; Р3 – получили отрицательный результат, что обусловило потери, связанные, в данном случае, с расходами проведения ОКР.
Пусть ставятся задачи: обеспечить одновременно оптимизацию режимов пилотирования с целью снижения потерь расхода топлива, а также предупреждение критических режимов пилотирования с целью минимизации потерь техники и расходов на ремонт ее в после аварийный период.
Пусть контролю и ограничению подлежат n параметров х1, х2, …, хn траектории полета самолета и состояния внутренних систем. В этом случае с помощью системы предупреждения критических режимов мы планируем осуществить предупреждение режимов полета, при которых параметры хi принимают критические значения.
При этом финансовые (материальные) потери, обусловленные возникновением событий Bi = (xi > xi кр), где xi кр – критическое значение параметра xi, обозначим через Gip. В процессе проектирования мы можем оценить стоимость Цi и массу mi такой системы. Управление расходом топлива q будем осуществлять с помощью системы оптимизации режимов пилотирования.
Суммарная стоимость СПКР, СОРП (ЦΣ) и масса (mΣ) зависят от объема обрабатываемой информации, в том числе от количества контролируемых, ограничиваемых и управляемых параметров xj.
Итак, задача состоит в оптимизации ψ = (Р*, q, ЦΣ, mΣ). Трудность решения такой задачи заключается в том, что для минимизации (Р*, q) необходимо контролировать все параметры хi траектории полета, в том числе возмущающие факторы, например, ветер. С другой стороны, для минимизации (ЦΣ, mΣ) требуется ограничивать количество параметров хi . Это противоречие сегодня решается конструкторами ЛА; СПКР; СОРП на интуитивном уровне, путем анализа потребностей рынка и возможностей производственно-технологической базы авиастроения. Возможно, такой подход оправдан. Однако необходимо иметь аналитические методы, позволяющие проводить анализ технического риска.
Согласно приведенной схеме процесса создания и эксплуатации пассажирского самолета, необходимо обеспечить заданную величину прибыли Z5 = Z (рис. 1.8). В качестве критической величины прибыли выступает Zкр = 0, т. е. когда прибыль равна нулю. В рамках рассматриваемой задачи прибыль Z зависит от расхода топлива (q) и потерь техники (α) в процессе эксплуатации самолета:
Z = f(q, α). (1.3)
Часто последнюю зависимость можно представить в более простом виде:
Z = k1φ1(q) + k2φ2(α),
где f, φ1, φ2 – непрерывные функции своих аргументов.
Если отклонение фактического расхода q от оптимального qopt мало, то возможно дальнейшее упрощение зависимости (1.3):
Z = k3Δq + k4Δα,
где Δq = q – qopt; Δα = α – αopt.
Отметим, что частота Pi выхода параметров траектории xi в критическую область позволяет вычислить αi = PiЦi, где Цi – стоимость работ по восстановлению техники, обусловленных выходом параметра xi в критическую область.
В процессе эксплуатации самолета, оборудованного СОПР и СПКР, возможны следующие события [6]:
на выходе СПKР: AX,BX,CX,DX;
(1.4)
на выходе СОПР: Aq,Bq,Cq,Dq;
где
AX = (x < xкр, xизм < xпрдоп), BX = (x < xкр, xизм > xпрдоп),
CX = (x > xкр, xизм > xпрдоп), DX = (x > xкр, xизм < xпрдоп),
Aq = (qф < qкр, qр < qпрдоп), Bq = (qф > qкр, qр > qпрдоп),
Cq = (qф < qкр, qр > qпрдоп), Dq = (qф > qкр, qр < qпрдоп).
При этом qкр есть критическое значение расхода топлива, начиная с которого полет совершается без прибыли, бизнес несет финансовые потери, прибыль Z или равна нулю, или отрицательная; qпрдоп – допустимая величина расхода, отличающаяся от qкр на некоторый запас, обусловленный влиянием погрешностей бортовых измерительных систем и систем оптимизации.
Введем событие A1(Z), которое происходит тогда, когда выполняются одновременно два события А1 = (Ах,Аq). Обозначим через Р1 вероятность такого события и тогда
P1 = P(Ax, Aq) = P(Ax) + P(Aq)
в силу независимости Ах и Аq из соотношения (1.4).
На практике возможны различные сочетания событий из соотношений (1.4) вида
А1 = (Ах, Aq); А2 = (Ах, Вq); А3 = (Ах, Cq); А4 = (Ах, Dq);
А5 = (Bх, Aq); А6 = (Bх, Вq); А7 = (Bх, Cq); А8 = (Bх, Dq);
А9 = (Cх, Aq); А10 = (Cх, Вq); А11 = (Cх, Cq); А12 = (Cх, Dq);
А13 = (Dх, Aq); А14 = (Dх, Вq); А15 = (Dх, Cq); А16 = (Dх, Dq).
Для всех Аi характерно событие Z < Zкр, где прибыль Z в соответствующем режиме полета будет меньше критической. Кроме того, эти соотношения включает события, при которых мы отказались от выгодного (оптимального) режима. Каждому событию Аi соответствует своя вероятность Pi = P(Ai). В силу независимости Ai между собой вероятность суммы всех событий Ai имеет вид
Задача проектировщика состоит в том, чтобы подобрать такой комплекс: ЛА, двигатель, СПКР, СОПР, для которого вероятность Р1 достигала бы максимума, а вероятность – минимума.
При расчете потерь прибыли следует включить стоимость по величине расхода q; стоимость восстановительных работ, обусловленных разрушениями при аварии самолета; стоимость работ по созданию самолета, СПКР, СОПР, двигателя; прибыль, полученную в процессе эксплуатации самолета. При этих условиях вероятность Р(А) представляет собой вероятность невыполнения поставленной цели, подлежащую выполнению, т. е. технический риск. С целью анализа, прогнозирования и управления риском необходимы математические модели, позволяющие установить явные зависимости между Pi(Ai) и параметрами контроля и управления на различных этапах жизненного цикла ЛА.
1.5. Этап научно-исследовательских работ
1.5.1. Показатели риска
На этапе научно-исследовательских работ, cвязанных с созданием самолета, мы работаем, как правило, с математическими моделями. При этом мы имеем дело с двумя видами моделей, М1 и М2 (рис. 1.11), проек тируемого самолета. Модель М1 описывает функционирование реального объекта без погрешностей, М2 – модель, принятая при расчетах.
Рис. 1.11
Для рассматриваемой ситуации
x(t) = Ψ(z, A, t), y(t) = Ψ1(z, A, δ, t),
где Ψ, Ψ1 – операторы описывающие модели М1 и М2 соответственно; δ(t) – погрешность модели или метода разработанного в теории; А – вектор заданных параметров, в том числе случайных возмущающих факторов.
О векторе x(t) на этапе НИР, как правило, мы имеем мало информации, но можем предположить, что в общем случае хi(t) ≠ yi(t) . При этом x(t) есть фактическое значение параметров состояния объекта хф в силу того, что в модели М1 учтены все внутренние и внешние возмущающиеся факторы и особенности систем его насыщающих. При этом процессе y(t) получен на выходе модели М2, которую нам дала теория (наука). Именно этой моделью мы владеем и пользуемся при проведении НИР, а на выходе ее имеем y(t) – расчетный или оценочный процесс. В этой модели учтена только часть внешних и внутренних возмущений (факторов), следовательно, модель М2 приближенно описывает изменение параметров ее состояния во времени.
Имея ввиду сказанное, получим показатель научно-исследовательского риска.
Рис. 1.12
Рис. 1.13
Пусть параметр х проектируемого объекта ограничен величиной по максимуму xвдоп (рис. 1.12). При этом справедливо x ≤ xвдоп, а G1 есть область допустимых значений х. В качестве такого параметра могут выступать, например, километровой расход топлива q, перегрузка, угол атаки. Ограничение на х может быть и снизу, т. е. ограничения для минимального значения, например, дальности полета L, прибыли инвестиционного проекта. При этом должно выполнятся условие x ≥ xндоп, и область допустимых значений х есть G2 (рис. 1.13).
При проведении научно-исследовательских работ параметры динамической системы выбираем таким образом, чтобы выполнялось условие x ≥ xвдоп или x ≥ xндоп. Такой подход возможен, если х известно и детерминировано. Однако это не так в силу того, что объект (самолет) подвержен возмущающим факторам как внутреннего, так и внешнего происхождения. В этих условиях параметр состояния х представляет собой случайный процесс. При этом у имеет вид y(t) = x(t) + δ1x и также является случайным процессом, δ1x – погрешность оценки.
За счет δ1x мы обязаны ввести запас Δв = xвдоп – (x0доп)в или Δн = xндоп – (x0доп)н (рис. 1.14).
Рис. 1.14
В силу случайности δ1х, у, х возможен ряд событий, представляющих собой потери инвестора (случай ограничения сверху, например, по расходу топлива):
A1 = (x ≥ xвдоп; y ≤ (x0доп)в), A2 = (x ≤ xвдоп; y ≥ (x0доп)в),
A3 = (x ≥ xвдоп; y ≥ (x0доп)в).
Вероятности этих событий
P1 = P(A1) = P(x ≥ xвдоп; y ≤ (x0доп)в),
P2 = P(A2) = P(x ≤ xвдоп; y ≥ (x0доп)в),
P3 = P(A3) = P(x ≥ xвдоп; y ≥ (x0доп)в).
Вероятность Р1 представляет собой численную величину того, что расчетный параметр фактически выполненного объекта будет превышать его допустимое значение, что неприемлемо либо по экономическим признакам, либо из условий безопасности эксплуатации объекта.
Вероятность Р2 представляет собой число, характеризующее следующее событие: объект забракован из-за невозможности удовлетворить ограничениям на параметр х, хотя эти ограничения для реального объекта выполняются.
Вероятность Р3 представляет собой событие, связанное с невозможностью удовлетворения ограничений, накладываемых на параметр х, т. е. разрабатываемый проект не может быть реализован, а средства, затраченные на проведение НИР, есть убытки, которые представляют риск инвестора.
Вероятность Р1 назовем вероятностью принятия убыточного проекта, вероятность Р2 – отклонением прибыльного проекта, Р3 – отклонением убыточного проекта. В совокупности Р1, Р2 и Р3 являются показателями научно-исследовательского риска, т. е. Рнр = (Р1, Р2, Р3). Запись в явной форме Р1, Р2, Р3, с помощью плотностей вероятностей W(x,y,t) рассматривается в работе [6].
В частном случае вероятности Р1, Р2, Р3 могут быть вычислены с помощью плотностей вероятностей W(Δx,δx), где Δx = х – mх; mх – математическое ожидание х; y = x + δх; δх – погрешность оценки х. Если Δx и δх независимы, то получаем еще более простые соотношения: W(Δх, δх) = W1(Δх)W2(δx).
1.5.2. Управление техническим риском на этапе проектирования
Развитие методов проектирования с учетом требований экономичности и безопасности полетов получило развитие только в последнее время и в настоящее время не достигло должного уровня. В результате иногда наблюдается либо перерезервирование в некоторых системах, т. е. утяжеление и удорожание самолета, либо, наоборот, вводятся конструкционные несовершенства, которые снижают эксплутационные качества, самолета. И те, и другие недостатки приходится устранять в процессе испытаний и доводки самолета, что существенно увеличивает цикл его создания и стоимость разработки, т. е. ведет к большим потерям и соответственно увеличивает технический риск.
Управление техническим риском, как правило, осуществляется с целью минимизации его численной величины. Процесс минимизации осуществляется различными способами и в различной степени на всех стадиях жизненного цикла летательного аппарата: научно-исследовательские работы; проектирование; опытно-конструкторские работы; серийное производство; эксплуатация. Существенно здесь то, что величина потерь характеризуется сочетанием этапа полета, типа или вида внутренних и внешних возмущающих факторов, характерных для ожидаемых условий эксплуатации. Под характерными ожидаемыми условиями эксплуатации подразумевается состояние внешней среды (барометрическое давление, температура наружного воздуха, направление и скорость ветра и т. п.) и эксплуатационные факторы (возможные конфигурации самолета, масса и центровка самолета, режим работы двигателя и т. п.). При этом суммарный риск обусловлен также внутренними возмущающими факторами, источниками которых являются пилот, функциональные системы обеспечения полета, их отказы и погрешности.
На рис. 1.15 представлено разделение суммарного технического риска Р(z) на две категории: Р(А) – потери, обусловленные появлением опасных режимов пилотирования и связанным с ними выходом из строя (полного или частичного) двигателя, самолета, его бортовых систем; Р(В) – экономические потери, обусловленные, например, перерасходом топлива. Каждая из составляющих Р(z) = P(A) + P(B) разбивается на P(A) = Р1i и P(B) = P2i согласно количеству этапов полета.
Рис. 1.15
Как указано выше, наиболее опасным этапом является посадка. Этот этап мы рассмотрим в плане управления техническим риском. Одним из главных моментов планирования и координации работ по летной годности самолета является выбор тех показателей, которые будут критериями обеспечения требуемого качества, свойств и условий летной годности на рассматриваемом этапе проектирования. Для пассажирского и транспортного самолетов такими показателями могут быть допустимые и расчетные для данного этапа полета вероятности возникновения различных ситуаций (безопасности и экономичности) i-й функциональной системы: (Pj)доп, (Pj)расч.
Порядок расчета, контроля и корректировки перечисленных показателей можно охарактеризовать как процесс управления риском при проектировании на трех уровнях:
1. ОКБ в целом Рдоп и Ррасч;
2. подразделения-разработчики функциональных систем Pjдоп, Pjрасч;
3. конкретные исполнители: конструкторы узлов, блоков, агрегатов проектируемой системы Pjдоп, Pjkрасч.
При эскизном проектировании в этом процессе можно выделить четыре основных этапа:
1. распределение заданных (допустимых, расчетных) показателей между функциональными системами и комплексами проектируемого самолета;
2. распределение показателей между подразделениями-разработчиками функциональных систем с целью более оперативного управления показателями внутри подразделения, ответственного за проектирование нескольких систем;
3. выполнение экспресс-анализа показателей риска;
4. определение допустимых вероятностей, в том числе функциональных отказов.
Рассмотрим этап посадки. Отметим что, как правило, на этапе посадки потерями топлива пренебрегают. При анализе риска в процессе посадки выявляются возможности реализации опасных условий полета. Этот фактор риска при посадке учитывает возможность приземления за пределы взлетно-посадочной полосы, наличие избыточной вертикальной, боковой или продольной составляющей скорости полета и ситуацию, требующую уход на второй круг после срыва захода на посадку.
В нижеследующей таблице приведены значения вероятностей выхода параметров движения за критические значения для самолета ТУ-154М с АБСУ-154-3.
В колонке «в среднем» даны вероятности событий в предположении, что все случайные параметры замкнутого контура управления изменяются в соответствии со своими законами распределения. В колонке «в предельном случае» даны вероятности событий, когда один из случайных параметров находится на своем предельно допустимом значении, а все остальные изменяются в соответствии со своими законами распределения.
На рис. 1.16 приводится детализация составляющих суммарной вероятности катастрофы, обусловленной наличием систем (посадка в автоматическом режиме) при посадке. При этом роль отказов бортового оборудования (внутренние возмущающие факторы) чрезвычайно велика и сравнима со всеми остальными возмущающими факторами, вызывающими реализацию опасных траекторий движения с выходом в критическую область.
Рис. 1.16. Составляющие средней вероятности катастрофы при посадке
На рис. 1.17 представлены составляющие суммарного риска при заходе на посадку самолета по категории II. На этом этапе полета нет особой необходимости учитывать потери в расходе топлива. Однако отказы систем играют здесь важную роль, и потому их выделяют отдельной величиной, которая при проектировании должна быть заданной для всех подсистем обеспечения захода и выполнения посадки.
Рис. 1.17. Составляющие суммарного риска при заходе на посадку по категории II
1.6. Модели для анализа технического риска
Развитие среды жизнедеятельности, осуществляемое человеком, включает следующие этапы:
– построение моделей, адекватных исследуемым процессам, и разработку соответствующих методов исследования свойств этих процессов, новых физических явлений и принципов;
– создание новой более совершенной техники, технологических процессов на основе новых физических явлений и принципов, описываемых соответствующими моделями, так, например, в космической связи, в передаче информации по каналам сети Internet.
Как правило, при создании технических, технологических объектов и систем разработчики стремятся использовать известные математические модели. В зависимости от уровня знаний в рассматриваемой области возможны следующие ситуации:
1) математическая модель изучаемого процесса отсутствует, но есть возможность использовать известные модели для ее построения;
2) нет возможности использовать известные модели, но есть возможность для ее построения применить известные законы и закономерности;
3) необходимо построить технический или иной объект, в основу работы которого закладываются новые законы и закономерности (ракеты, и т. п.); эта ситуация характеризуется максимальной неопределенностью и связана с наибольшими потерями и соответственно наибольшей величиной технического риска.
Известные математические модели были созданы и создаются для исследования
– технико-технологических процессов с тем, чтобы создать управляемые процессы с более низким уровнем затрат или более высокой эффективностью в сравнении с существующими;
– динамики и эффективности современных воздушных и космических летательных аппаратов и грузопассажирского транспорта;
– функционирования человека в человеко-машинных системах, в том числе при управлении;
– биологических систем, в том числе подсистем жизнеобеспечения: центральной нервной; сердечно-сосудистой.
История создания моделей связана с процессом развития научной мысли, познания среды жизнедеятельности и включает в себя ряд этапов. Первичными в этом процессе всегда были результаты наблюдений процессов и явлений, лежащих в основе построения модели через абстрактное мышление. Модель создается для практической деятельности человека и, в частности, для создания развития технико-технологической среды.
Построение модели требует большого количества исходной информации, которая создается в результате отдельных наблюдений, которые являются случайными в силу: 1) невозможности учесть все многообразие фактов, в том числе возмущающих, действующих на изучаемый объект, даже если он простой; 2) сложных взаимосвязей между этими факторами и объектом; 3) несовершенства естественных и искусственных средств наблюдения и контроля. В процессе таких наблюдений и контроля производится формализация для построения основных признаков, связей, закономерностей, присущих изучаемому объекту.
Как правило, для одного и того же объекта в зависимости от конкретных требований практики и типа решаемой задачи строят несколько моделей, с помощью которых осуществляют формализацию различных функций этого объекта или внешних воздействий на него. Это позволяет получить достаточно простые модели, что упрощает процедуру анализа и синтеза технико-технологических систем и объектов, уменьшая экономические затраты.
Процесс построения моделей для осуществления формализации при их описании основан на большом разнообразии подходов и языков. Ранее на первом этапе формализации осуществлялось, в основном, содержательное описание явления и, чаще всего, использовался разговорный язык (Греция), удобный как для изложения результатов исследований, так и их передачи для использования и развития. Однако разговорный язык не мог описать, зафиксировать и надежно передать большое число качественных и количественных сторон явлений. Развитие, расширение и углубление процессов исследования привело к необходимости создания различных языковых форм, более экономных по объему передаваемой информации, понятных и заполнявших содержательное описание результатов наблюдений. Так появились многие языки схем, уравнений, карт, чертежей, которые используются и по сей день для построения моделей и их формального описания. Это связано с тем, что при их построении применяются различные способы наблюдения, обобщения, формализации и представления результатов наблюдения.
Во второй половине XX в. в связи с развитием технико-технологической среды возникла необходимость в новых принципах построения моделей формализации результатов наблюдений, связанных с управлением целенаправленной деятельностью человеческого общества, живыми организмами и искусственными, созданными человеком, машинами и механизмами, которые представляют математические модели.
Существенно, что при решении указанных проблем, модель, принятая для целей проектирования, на этапе научно-исследовательских работ часто сильно отличалась от реального объекта, включающего внешние и внутренние возмущающие факторы, отказы техники, старение и деградацию, что значительно уменьшало или сводило на нет эффективность разработанной системы. В настоящее время математические модели включают в себя следующие модели:
– детерминированные;
– стохастические, в том числе с неопределенностью и с идентификацией параметров;
– оптимальных процессов и т. п.
Согласно требованиям практики, в теории вероятностей и математической статистике получили развитие два способа построения математических моделей исследуемых физических объектов.
Первый – априорный способ введения модели. Он основан на том, что физический объект описывается некоторой вероятностной или статистической моделью с заданной структурой, включающей в себя неизвестные параметры. Эта модель получена из дополнительных исходных соображений, например, на основе опыта описания аналогичных объектов той или иной физической природы с учетом физико-химических и других особенностей объекта, а также физически обоснованных и экспериментально проверенных моделей элементарных составляющих объекта.
На втором этапе исследований с объектом проводят эксперимент, в результате которого получают экспериментальные данные, содержащие информацию о реальных значениях исследуемых параметрах объекта.
На третьем этапе производится соответствующая обработка результатов измерений неизвестных параметров модели. На четвертом – проверка адекватности полученной М2 (рис. 1.11) и фактической М1 моделей объекта, например, путем сравнения выходных сигналов y модели М2 и х модели М1, если задан входной сигнал z(t). При этом предполагается, что физический объект существует, и на нем можно проводить необходимые исследования.
На пятом этапе, в случае необходимости, корректируется модель М2, а точнее ее параметры.
При таком подходе, основанном на использовании аксиоматики А.Н.Колмогорова и являющимся априорным, есть ясная физическая интерпретация параметров модели, возможность широкого использования этой модели для анализа, синтеза, численного моделирования, прогнозирования и управления поведением объекта.
Вторым способом построения моделей является априорный способ, при котором структура модели не задана. При этом построение модели основано на подборе и подгонке по экспериментальным данным как вида модели из заданного класса моделей, так и ее параметров, которые наилучшим или приемлемым образом соответствуют результатам статистической обработки экспериментальных данных. При этом на модель не накладывается ограничений. Так, например, если известно, что случайный процесс x(t) – стационарный и эргодический с нормальной плотностью вероятностей и корреляционной функцией kx(τ) = σе–|τ|, то в этом случае мы можем поставить процессу x(t) в соответствие процесс y(t), являющийся выходом модели М, описываемой дифференциальным уравнением вида
где η(t) – случайный процесс типа белого шума, моделируемый с помощью генератора случайных чисел [21].
Апостериорный подход используется, например, для кибернетики или информационной модели физического объекта, когда объект рассматривается как система, осуществляющая преобразование входного воздействия {η(t)} в выходной y(t). При этом она считается полностью описанной, если имеется математическая модель, т. е. модель в абстрактной форме, не зависящей от физической природы входных и выходных переменных и процессов, происходящих в системе.
В отличие от вероятностных моделей в статистических (Фон Мизеса) модель строится на основе характеристик, получаемых предельным переходом от эмпирических характеристик при неограниченном увеличении объема экспериментальных данных. Адекватность использования таких моделей для описания физических объектов также проблематична, как и вероятностных, а в нестационарных условиях – не очевидна. Ответ на вопрос, какие модели и какие способы их получения использовать, в конкретном случае зависит от поставленной цели исследования, вида решаемой задачи, априорной неопределенности. Общим для этих методов является применение единых вероятностно-статистических методов, используемых в статистических измерениях, для характеристик и параметров реальных физических объектов. Основной вклад в технический риск происходит на этапе обработки материалов испытаний, т. е. при систематическом анализе, когда решаются следующие задачи:
1) выбор вида модели из заданного класса;
2) статистические описания характеристик модели или измерения ее параметров;
3) проверка статистических гипотез;
4) построение итоговой модели М, включая подгонку модели;
5) анализ свойств модели М;
6) статистическое моделирование на полученной модели М.
1.7. Упрощение математических моделей
Одна из причин отличия реальной математической модели М2, принятой в расчетах, от идеальной М1, которая описывает функционирование реального объекта без погрешностей, состоит в том, что идеальная модель для расчетов упрощается. Это упрощение происходит по различным причинам и, в частности, из-за невозможности применить М1 при расчетах на этапе НИР из-за ее сложности. Рассмотрим пример такого упрощения и допускаемые при этом погрешности.
Пусть создана математическая модель М1, которая представляет собой систему из n уравнений, заданная в пространстве RN с координатами x = (x1, x2, …, xN), вида
где 1 ≤ i ≤ n; 1 ≤ l, m ≤ N, 1 ≤ k ≤ α; x, t – соответственно пространственные и временные переменные.
Таким образом, (1.5) представляет собой систему интегро-дифференциальных уравнений для объекта с распределенными параметрами. Решения ψj = ψj(t,x;Ck) описывают состояние объекта или системы в некоторый момент времени t, поэтому их называют переменными состояниями. При этом предполагается, что уравнения Fi = 0 зависят от параметров Ck, например, числа Маха (М), числа Рейнольдса в авиации, структурных констант, которые могут влиять на свойства решений ψj, поэтому их называют управляющими параметрами.
Проблема исследования решений системы уравнений (1.5), даже если речь идет лишь о том, как зависят эти решения от управляющих параметров Сk, является исключительно сложной. По этой причине вводится ряд последовательных предположений, каждое из которых вносит соответствующие погрешности. К числу таких предположений относятся изменения в структуре модели М:
1. Погрешность δ11, обусловленная пренебрежением элементов объекта, описываемых интегралами. В результате модель (1.5) заменяется новой моделью вида
которая представляет собой систему нелинейных уравнений в частных производных.
2. Погрешность δ12, обусловленная новым упрощающим предположением, согласно которому модель не содержит пространственных производных любого порядка, т. е.
которая описывает объект с сосредоточенными параметрами.
3. Погрешность δ13, обусловленная переходом от модели (1.7) к модели вида
когда модель содержит производные не выше первого порядка и, кроме того, эти производные входят в упрощенную функцию Fi специальным каноническим образом. Систему уравнений Fi = 0 называют динамической системой.
4. Погрешность δ14, обусловленная предположением, что функция fi в (1.8) полностью не зависит от времени, тогда получаем модель
которая называется автономной динамической системой. Для таких систем, когда k ≤ 4, имеются некоторые интересные результаты в области свойств решений.
5. Погрешность δ15, обусловленная переходом от (1.9) к системе вида
когда все функции fi могут быть заданы антиградиентом по отношению к ψi некоторой потенциальной функцией V(ψj; Ck). Свойства таких систем изучены довольно глубоко.
6. Погрешность δ16, обусловленная переходом к модели, описывающей состояния равновесия dψi/dt = 0 градиентных динамических систем. В этом случае модель изучаемого процесса принимает вид
Эти уравнения могут не иметь решений, если V(ψ) = ψ; иметь одно решение при V(ψ) = ψ2; более, чем одно решение, если V(ψ,c) = ψ4 + cψ2: при c < 0 – имеется три решения, а при c > 0 – одно решение. При этом задача исследования заключается в том, каким образом состояние равновесия ψj(Ck) потенциальной функции V(ψj;Ck) изменяется при изменении управляющих параметров Ck.
Среди других методов исследования системы, основанных на переходе от (1.5) к более простым моделям, что обуславливает погрешности типа δ1j, рассмотренные выше, укажем:
1) метод линеаризации;
2) вариационные методы;
3) метод прогонки или метод начальных параметров.
Важным моментом в этих методах является приближенная замена дифференциальных уравнений (1.6) в частных производных системой обыкновенных дифференциальных уравнений, характеризующих динамические свойства системы с конечным числом степеней свободы.
1.7.1. Стохастические модели и их погрешности
Пусть рассматриваются абстрактные вероятностные или статистические модели физических объектов. Модель исследуемого объекта и модель полученных данных отличны. Эти отличия обусловлены тем, что экспериментальные данные отражают не только состояния изучаемого объекта, значения его характеристик и параметров, но и условия проведения, связанные с внешними возмущающими факторами δ1, а также упрощения и допущения, использованные при обработке эксперимента δ2. Если первые связаны с внешней средой, то вторые – со структурой объекта, с ее упрощением, которые проявляют себя в моделях, алгоритмах, построенных на основе экспериментальных данных.
Рассмотрим подробнее этот процесс, в основу которого положен морфологический анализ, широко используемый при исследовании и проектировании различных технических объектов [11]. Разберем задачу измерения характеристик случайных процессов, т. е. первый этап построения модели, например, в форме дифференциальных уравнений, позволяющих прогнозировать состояние изучаемого объекта, исследовать его устойчивость, прогнозировать численную величину показателей риска. Этапы таких работ включают следующие погрешности:
1. Получение исходных данных: принцип получения исходных данных (ансамблевый или траекторный) и оператор усреднения, а также количество обращения к исходным данным – погрешности δ21.
2. Схема сравнения с образцовой величиной: количество сравнений (одно или несколько) представляет собой этап, в котором производится сравнение измерительных преобразований; вычислительные преобразования – погрешности δ22.
3. Схема дискретизации траектории по времени: непрерывная; дискретная – погрешность δ23.
4. Схема квантования значений реализаций по уровню: без квантования (непрерывнозначный; аналоговый); с квантованием всех значений (цифровой); с квантованием части значений – погрешность δ24.
5. Добавление вспомогательных сигналов – погрешность δ25.
6. Преобразование по формулам оценивания – погрешность δ26:
– количество значений, используемых в качестве результата измерения одной величины (одноинтервальный – две точки) – δ261;
– способ оценивания (параметрический, непараметрический) – δ262;
– объем исходных данных (фиксированный или переменный) – δ263;
– принцип реализации вычислительных операций (детерминированный, стохастический) – δ264;
– принцип организации вычислительного процесса (формул связи текущего значения с предшествующим): одношаговый; многошаговый – δ265;
– метод вычислений: с фиксированным числом операций; с переменным числом операций – δ266;
– вид измеряемой величины (скалярный, векторный): детерминированная, постоянная; реализация случайного элемента (байесовский) – δ267;
– суммарная погрешность на данном этапе имеет вид: .
7. Схема изменения аргумента модели в процессе усреднения: аргумент остается постоянным; аргумент изменяется по заданному закону – погрешность δ27.
8. Получение конечного результата: непосредственно по исходным данным (прямой); пересчетом результатов измерения других характеристик, т. е. косвенным методом – погрешность δ28.
9. Аналитическое описание результатов измерения: без аналитического описания; с аналитическим описанием – погрешности δ29.
10. Учет априорных и апостериорных данных: наличие адаптации (без и при наличии) к априорным и апостериорным (первичным, вторичным) данным – погрешности δ210.
Суммарная погрешность модели физической системы, которую мы заложили при построении ее математической модели в процессе эксперимента, записывается в виде
Предположим, что построена математическая модель случайного процесса, т. е. построен, например, алгоритм, связывающий выходной процесс x(t) с входным процессом y(t), характеристики которого заданы в виде x(t) = ψ(y(t),δ2i), где ψ – оператор преобразования. Наличие погрешностей δ2i заставляет нас искать показатели качества алгоритма, которые являются характеристикой соответствия алгоритма его назначению, т. е. пригодность алгоритма для получения решения поставленной задачи и близость достижения цели.
Рассмотрим критерий применимости упрощенных математических моделей изучаемых динамических систем. Зададим Р1, Р2 и Р3 Предположим, х = mх + Δx, где mх – математическое ожидание х; Δx – отклонение х от его среднего значения (mx). В этом случае вероятность Р1 можно записать так [6]:
где а = хвдоп – mx; b = ходоп – mx – δx; W1, W2 – плотности вероятностей Δх, δх.
Выбирая ту или иную модель М2, мы изменяем W2(δх), оставляя W1(Δх) неизменной. Можно показать, что при возрастании дисперсии σ2(δх) условие Р1 = const выполняется, если запас Δ = хдоп – ходоп увеличивается. Очевидно, что существует то значение σ2(δх), при котором Δ очень велико и эффективность применения такого объекта недопустимо мала (рис. 1.18).
Рис. 1.18
Таким образом, переход к упрощенной модели при проведении НИР приводит к уменьшению области допустимых состояний ограничиваемого параметра х.
Задача разработчика проекта: наложить ограничения на Δ и принять ее равной Δ*. Все модели, для которых Δ ≤ Δ*, будут допустимы для целей анализа и синтеза. Таким образом, проблема анализа риска и управления им включает в себя: построение плотностей вероятностей фактических хф = mx(t) + Δx(t) и измеренных (оценочных) хизм = хф(t) + δх(t) значений процессов, подлежащих контролю и ограничению. Для построения плотностей вероятностей W1(xф) или W1(Δx), а также W2(xизм) или W2(δx) необходимы соответствующие математические модели. Построению таких моделей посвящена данная работа.
1.7.2. Модели плотностей вероятностей случайных процессов
Для расчета технического риска необходимо иметь плотность вероятностей ограничиваемых параметров, представляющих собой, как правило, случайные процессы. Для аналитического построения искомых плотностей вероятностей необходима математическая модель случайного процесса, а также внешних и внутренних возмущающих факторов. Для решаемой задачи математическая модель движения самолета представляет, в достаточно общем виде, систему нелинейных дифференциальных уравнений, на вход которой поступает белый шум [6]. При этом на выходе системы (модели) получают многомерный марковский процесс, что существенно упрощает нахождение плотностей вероятностей W(·) аналитическим методом. С этой целью используется математический аппарат уравнений Фоккера-Планка-Колмогорова (ФПК). Однако при размерности вектора x(t) ограничиваемых параметров больше трех, решение ФПК-уравнения представляет собой трудноразрешимую задачу. В этом случае используют различные процедуры, упрощающие вычисления W(xi,t), где xi – процесс, подлежащий контролю, ограничению и прогнозированию, т. е. вычислению xi(t + Δt), где Δt > 0, а также моделированию его с помощью известных технических средств.
Пусть задана техническая система (объект), на выходе которой возможна регистрация процесса xi = z(t). Требуется построить математическую модель процесса z(t). Пусть поставлен эксперимент, в котором измерены значения z(t) с погрешностями δ2Σ, в результате получено z(t) = zф + δ2Σ, где zф – фактическое значение z(t); для z(t) построены одномерная плотность распределения ω(z) и корреляционная функция Bz(τ). Построение математической модели процесса z(t) осуществляется в три этапа.
На первом этапе производится аппроксимация одномерной плотности распределения ω(z) изучаемого процесса z(t) плотностью S-распределения Джонсона [19, 20]. При этом исключаются из рассмотрения те процессы, которые не могут быть аппроксимированы указанным образом. Поскольку класс плотностей S-распределения Джонсона образован тремя семействами (SB, SU и SL-распределения), то аппроксимация распределений ω(z) заключается в выборе соответствующего семейства плотностей S-распределения и в определении его параметров. Процедура такой аппроксимации описана в работе [6]. Следует только добавить, что исходной информацией для аппроксимации являются первые четыре момента распределения процесса z(t), которые, в общем случае, следует предварительно вычислить. Для некоторых видов плотности распределения ω(z) эти моменты указаны в справочной литературе. После того, как для ω(z) определено соответствующее S-распределение, для построения математической модели искомого процесса z(t) необходимо описать нормированный гауссовский процесс y(t) и подвергнуть его нелинейному преобразованию Джонсона. Указанное преобразование будет иметь вид
если плотность ω(z) представима в виде плотности SB-распределения Джонсона;
если плотность ω(z) пред ставима в виде плотности SU-распределения Джонсона;
если плотность ω(z) представима в виде плотности SL-распределения Джонсона.
Входящие в выражения (1.10)÷(1.12) параметры λ, γ, η и ε являются параметрами, соответствующими S-распределению Джонсона. Таким образом,
z(t) = ψj(y(t)), (113)
где 1 ≤ j ≤ 3 в зависимости от вида плотности S-распределения Джонсона, аппроксимирующей заданную плотность ω(z).
Для того, чтобы случайный процесс z(t) имел заданную корреляционную функцию Bz(τ), необходимо подобрать соответствующую функцию (обозначим ее ρy(τ)) нормированному гауссовскому процессу y(t), подвергаемому нелинейному преобразованию ψj(·) (1 ≤ j ≤ 3). Поэтому на втором этапе проводится расчет корреляционной функции ρy(r), который осуществляется по формуле
где
Нn – полином Эрмита n-й степени. Расчет заключается в получении набора значений ρy(τi) корреляционной функции случайного процесса y(t), соответствующих выбранным значениям τi ее аргумента и удовлетворяющих (1.14). При проведении конкретных расчетов ряд, стоящий в правой части этого равенства, необходимо ограничить несколькими первыми членами; количество оставленных членов должно быть таким, чтобы обеспечивалась требуемая точность выдерживания указанного равенства.
На третьем этапе осуществляется представление случайного процесса y(t) в виде компонент двумерного линейного марковского процесса. При этом полученные на предыдущем этапе значения корреляционной функции ρy(τ) процесса y(t) сглаживаются выражением
ρy(τ) = e–h|τ|(cosωτ + βsinω|τ|). (115)
Обеспечение хорошей точности сглаживания позволяет, используя приведенную в [6, 18, 44, 47, 71, 74, 75] методику, представить случайный процесс y(t) в виде линейного двумерного марковского процесса (y1(t), y2(t)), порождаемого следующей системой стохастических дифференциальных уравнений:
где ξ(t) – гауссовский белый шум с нулевым математическим ожиданием и дельтаобразной корреляционный функцией.
Отметим, что часто встречающиеся на практике корреляционные функции вида
представляют собой частные случаи корреляционной функции (1.15), следовательно, процесс y(t), корреляционная функция которого принадлежит семейству (1.17), может быть представлен в виде компоненты y1(t) марковского процесса (1.16) после его соответствующего упрощения.
Таким образом, генерация стационарного процесса z(t) с заданными одномерной плоскостью распределения ω(z) и корреляционной функцией Bz(τ) может быть осуществлена путем генерации двумерного линейного марковского процесса (y1(t), y2(t)), описываемого системой уравнений (1.15), и последующего преобразования с помощью соотношения (1.13) компоненты y1(t) данного процесса, которое необходимо проводить после установления ее стационарности. При этом двумерный марковский процесс (y1(t), y2(t)) может быть получен, например, как выходной процесс с формирующего фильтра [69], возбуждаемого белым шумом с указанными статистическими характеристиками и построенного в соответствии с системой (1.16).
1.8. Цели и задачи работы
Оправдывает ли себя в эксплуатации тот или иной ЛА, окупятся ли расходы разработки, изготовления и эксплуатации тех или иных новых систем бортового оборудования или новых конструкций планера и двигателя? Технологии их изготовления (рис. 1.19) зависят от соотношения прибыли и потерь при эксплуатации, характеризуемые вероятностями P1 и P = (P2, P3, P4). При анализе P1 и P необходимо учитывать целевое назначение ЛА в силу различия, например, спортивных самолетов от самолетов гражданской авиации, для которых, в отличие от первых решающее значение имеет место выход параметров траектории полета в критическую область.
Рис. 1.19. (Здесь 1 – система встроенного контроля бортового оборудования; y = (x1 опт, …, хl опт, х1, …, хι, δ11, …, δ1α), W2(δ2) – внешние возмущающие факторы, в том числе от аэродромных средств посадки, управления воздушным движением, состояния атмосферы, СУ – система управления, ИИС – информационно-измерительная система)
Согласно сказанному выше, имеет место несколько групп параметров траектории полета, которые при эксплуатации ЛА подлежат контролю и управлению. Одна из групп А принадлежит области G1 оптимальных значений. Как правило, G1 узкая, а параметры, принадлежащие G1, должны измеряться и выдерживаться с высокой точностью. Вторая группа B должна принадлежать области G2 (достаточно большой) допустимых значений. Выход параметров B из области G2 связан с потерями техники, в том числе с ее восстановлением.
Параметры, принадлежащие A и B при эксплуатации ЛА, подлежат контролю и управлению. Кроме управляющих воздействий они изменяются под действием возмущающих факторов внутреннего и внешнего происхождения, которые обуславливают их выход из G1 и G2 и соответствующие потери. Как правило, указанные факторы не поддаются контролю или неконтролируемы по каким-либо причинам. Как только тот или иной возмущающий фактор начинают контролировать с помощью бортового оборудования, он переходит в разряд параметров состояния ЛА. Среди таких параметров можно рассматривать массу m и центровку xT ЛА.
Некоторые возмущающие факторы, например, параметры восходящих потоков, сдвиг ветра, мы можем только контролировать. Среди внутренних возмущающих факторов выделим следующие (рис. 1.19):
1) технологические факторы производства планера и двигателя, которые обуславливают некоторые колебания аэродинамических характеристик планера на величину ΔY, ΔX, ΔM, где Y, X, M – подъемная сила, сопротивление и момент ЛА соответственно, а также колебания газодинамических характеристик двигателя на величину ΔP, где P – тяга двигателя, которые обозначим C1, C2, C3, C4 соответственно;
2) погрешности систем контроля δ1 и управления, которые обусловливают отклонения параметров состояния ЛА от заданных значений на величину Δxi ;
3) погрешности систем оптимизации, обусловленные влиянием Δxi, отклонениями фактических параметров траектории (xф)i от их оптимальных значений (xopt)i, т. е. Δxi = (xф)i – (xopt)i, а также собственно их методическими и инструментальными погрешностями.
Задача состоит в разработке таких математических моделей, которые описывают по возможности всю совокупность параметров x и возмущающих факторов, влияющих на плотность вероятностей W(x,δ), служащей для расчета вероятностей Pi . Отметим, что как только подсчитаны вероятности Pi, им с достаточной степенью приближения можно поставить соответствующие частоты рассматриваемых событий и для данного класса ЛА в данных условиях его эксплуатации (стоимость одного рейса, прибыль от этого рейса) можно подсчитать потери и прибыль (в стоимостном выражении) от совокупности рейсов (например, 109) и оценить целесообразность (окупаемость) установки тех или иных новых систем повышения регулярности, безопасности, экономичности, в том числе расхода топлива.
Согласно сказанного выше, цели и задачи количественного анализа технического риска включают в себя
– построение математических моделей погрешностей функционирования оптимизатора расхода топлива ЛА, производства планера ЛА, вычислителя массы и центровки ЛА, измерителя скоростей полета и тяги несущего винта вертолета;
– построение с помощью математических моделей погрешностей, возникающих на всех этапах создания и эксплуатации ЛА, искомых плотностей вероятностей W(x,δ);
– вычисление количественных показателей технического риска с использованием полученных плотностей вероятностей W(x,δ);
– расчет потерь и прибыли по известным количественным показателям технического риска.
Указанным целям и задачам и посвящена данная работа. Методический подход достижения цели и решения задач иллюстрируется на ряде новых систем контроля, алгоритмов оптимизации на материалах расчетов W(x,δ) и летных испытаний систем.