Введение Человекомашинные системы (ЧМС) составляют значительную часть промышленного и иного оборудования, в частности корабельной техники. Обеспечение безотказности функционирования здесь практически всегда связано с обеспечением безопасности человека как важной и неотъемлемой части указанных систем. При этом существует взаимозависимость между техникой и обслуживающим ее персоналом. Нарушение работоспособности человека-оператора может приводить и приводит к негативным последствиям, порой катастрофическим. Соответственно, отказы технической части ЧМС могут подвергать существенному риску для жизни не только обслуживающий персонал, но и не связанных с данной техникой людей. Однако с точки зрения природы нарушений работоспособности и развития таких нарушений как в технической (ТС), так и в человеческой (ЧС) составляющих ЧМС существуют значительные различия. Указанные различия обусловливают возможность выполнения мероприятий по обеспечению безотказности функционирования ЧМС отдельно для человека-оператора и для обслуживаемой им техники. Такой подход является типичным и широко распространенным на практике. Вместе с тем игнорирование взаимосвязи обеих частей ЧМС может приводить к существенным искажениям при оценке состояния всей ЧМС в целом. Более целесообразным здесь может быть применение подхода, базирующегося на признании человека-оператора и техники равноправными партнерами, вклад каждого из которых в безотказность функционирования зависит от другого. При этом, если исходить из принципа минимакса (расчета на «наихудший» случай), надежность и безотказность всей системы в целом может быть определена как надежность и безотказность части, более склонной к отказам на заданном интервале эксплуатации и, следовательно, менее надежной. Постановка задачи Среди различных показателей надежности и безотказности технических систем вообще наиболее общим следует считать показатель вероятности безотказной работы на заданном интервале времени. Указанный показатель является одним из самых распространенных на практике. Однако его вычисление всегда связано с необходимостью сбора значительных объемов данных (априорных и апостериорных). Реалии эксплуатации ЧМС редко позволяют рассчитывать на получение потребной информации как в качественном, так и количественном отношениях. Иначе говоря, применение показателя вероятности безотказной работы и других подобных показателей статистического характера (среднее время безотказной работы, коэффициент готовности и т. п.) на практике затруднено. Более приемлемыми здесь могут быть показатели, пригодные для использования в условиях существенной ограниченности и неопределенности исходной информационной базы, в частности связанные при своем построении с минимаксным подходом. К таким показателям можно отнести показатель остаточного эксплуатационного ресурса . Указанный показатель, вообще говоря, был разработан для оценки потенциального временного ресурса безотказного функционирования технических систем с учетом параметрических отказов [1], поэтому сначала рассмотрим конструкцию при его использовании в задаче обеспечения безотказности функционирования ТС ЧМС. Пусть – вектор параметров, определяющих состояние ТС. Под воздействием дестабилизирующих факторов вектор y изменяется на интервале эксплуатации Т как неопределенный процесс . При нахождении где D – область работоспособности, ТС функционирует в заданном режиме. , где А, В – заданные ограничения; D – ортогональный параллелепипед. В результатах измерений параметров на интервале возможно присутствие аддитивной и статистически неопределенной ошибки : . (1) Относительно ошибки известно только, что , (2) где D – известный вектор, величины компонентов которого связаны со свойствами измерительной техники (и обслуживающего персонала), т. е. определяются паспортными данными средств измерения, а также условиями выполнения соответствующих работ. Остаточный ресурс эксплуатируемой ТС опишем как интервал времени где – момент времени последнего измерения на интервале – первый момент времени достижения границ области D, т. е. . Определение остаточного ресурса технической системы Для формирования критерия, характеризующего остаточный ресурс ТС, используем критерий «запаса работоспособности» , определяемый по каждой компоненте y(t) как расстояние в текущий момент времени между и ближайшей точкой на [1]. При этом оказывается возможным рассматривать изменения S, связанные с изменением y под воздействием дестабилизирующих факторов, как равноускоренное прямолинейное движение S, т. е. (3) где – значение j-го запаса работоспособности на начало эксплуатации ТС (соответствует номинальному значению состояния ТС по j-й компоненте y и задается техническими условиями на эксплуатацию ТС); – скорость изменения j-го запаса работоспособности в направлении к ближайшей точке на – ускорение изменения j-го запаса работоспособности в направлении к ближайшей точке на . В принципе зависимость (3) можно рассматривать как некоторую аппроксимацию найденную из вполне понятных физических соображений. Соответственно, результаты контроля , описываемые соотношением (2) и (3), правомерно считать и результатами контроля . При этом нахождение оценки , будет соответствовать нахождению по результатам контроля оценки т. е. С учетом зависимости (3) поиск практически сводится к поиску по результатам контроля оценок и . По известным несложно найти , причем по каждой координате или Ограниченность и неопределенность исходных данных в задаче нахождения определяет целесообразность использования здесь принципа минимакса [2]. C учетом указанного принципа для ТС в целом остаточный эксплуатационный ресурс можно определить как (4) В дополнение к зависимости (4) на основе принципа минимакса могут быть найдены гарантированные по достоверности относительно имеющейся совокупности исходных данных и без принятия каких-либо допущений или гипотез по статистическим свойствам указанной совокупности как интервальные, определяемые в виде неравенств (5) так и точечные оценки определяемые соотношением (6) Несложно убедиться, что оценка (6) является оптимальной в смысле минимума максимально возможной ошибки оценивания (7) Расширение спектра применений показателя остаточного ресурса на ЧС ЧМС Ключевым вопросом при построении показателя является вопрос о возможности формирования вектора измеряемых параметров , определяющих состояние ТС и чувствительных к изменениям этого состояния. Для ТС, пусть и достаточно сложных, ответ на такой вопрос несложно найти при анализе технической документации (инструкций) по эксплуатации практически любой, в том числе измерительной технической системы. Для биологических объектов, а следовательно, и для ЧС ЧМС, формирование вектора y может составить отдельную и сложную задачу. Указанная задача фактически является задачей функциональной медицинской диагностики, когда по измерениям некоторых определенных параметров необходимо определить наличие или отсутствие той или иной стадии заболевания, связанного с утратой работоспособности человеком. При решении такой задачи можно говорить и о том, что близость или отдаленность значений измерений рассматриваемых параметров от границ допустимых пределов (фактически от ) позволяет судить и о возможности дальнейшего участия данного персонала (людей) в процессе эксплуатации ЧМС. В большей своей части искомые параметры связаны с контролем сердечно-сосудистой системы человека как основы обеспечения жизнедеятельности и работоспособности всего человеческого организма в целом [3]. Одним из таких параметров является, например, показатель оптической плотности миокарда [4]. Для измерения данного показателя существует целая линейка интраскопических приборов, включая дешевые малогабаритные и общедоступные. Современные средства контроля других параметров сердечно-сосудистой системы человека во многих случаях также могут быть построены (или уже построены) как малогабаритные и общедоступные. При этом, при осуществлении контроля за состоянием ТС, целесообразно выполнять и контроль параметров ЧС ЧМС. Иными словами, задача формирования вектора в настоящее время получила достаточное разрешение. Следовательно, существует и очевидная возможность расширения спектра применений показателя на ЧС ЧМС. При этом формирование минимаксных оценок (4)–(7) как для ТС, так и для ЧМС сводится к одинаковым и описанным выше процедурам. Однотипность процедур формирования минимаксных оценок (4)–(7) позволяет рассмотреть и один общий алгоритм решения задачи определения остаточного эксплуатационного ресурса ЧМС с учетом функционирования обеих ее частей (ТС и ЧС). Указанный алгоритм может быть представлен следующим образом. С учетом зависимостей (1) и (2) по результатам контроля состояния ТС (и ЧС) можно сформировать систему неравенств: (8) Исходя из физики процессов изменения состояния ТС (и ЧС), систему (8) можно дополнить ограничениями на непринадлежность экстремума зависимости (3) интервалу эксплуатации ТС (и ЧС), определяемому как где и – соответственно заданные начальный и конечный моменты времени эксплуатации ТС, т. е. (9) Объединение (8) и (9) составляет ограничения в следующих задачах линейного программирования для нахождения минимаксных оценок неизвестных и (10) где – произвольная фиксированная точка в области . На основе результатов решения задач (10) несложно построить два комплекта оценок (4)–(7) для ТС – и ЧС ЧМС (каждый по отдельности). Учитывая партнерский характер взаимосвязей между рассматриваемыми составляющими (ТС и ЧС) ЧМС и исходя из принципа минимакса, можно найти оценки показателя для всей ЧМС в целом из следующего выражения: . Сравнение критериев управления эксплуатацией сложных технических систем Минимаксный подход, вообще говоря, является пессимистическим и поэтому позволяет получить пессимистические оценки состояния и выработать соответствующие мероприятия по управлению состоянием сложной технической (человекомашинной) системы. Затраты на такое управление могут превышать затраты на оптимистическое управление эксплуатацией или управление по статистическим критериям, и в частности на управление по критерию вероятности безотказной работы. При этом оказывается целесообразным сравнить критерии управления эксплуатацией, имеющие различный физический смысл: критерий вероятности безотказной работы на заданном интервале времени и критерий остаточного (гарантированного или минимаксного) эксплуатационного ресурса. Вообще говоря, показатель для проведения сравнения с показателем вероятности безотказной работы на заданном интервале времени следует рассматривать как родственный P(t). При этом можно рассматривать как некоторый случайный процесс, характеристики которого могут быть заданы или определены по результатам наблюдений. Кроме того, показатель следует пронормировать на интервал [0,1]. После нормировки значение показателя дает возможность рассматривать его как интегральную выборочную функцию распределения времени безотказной работы ТС, пусть как H. При этом значения можно интерпретировать как коэффициент присутствия любого отрезка на интервале Т и в области работоспособности объекта D – . Исходя из монотонного убывания присутствия отрезка пути , в области работоспособности объекта D, для любого момента времени , верно неравенство где – значения выборочной интегральной функции распределения времени пребывания некоторой компоненты в области D. Полученное неравенство позволяет оценить вероятность безотказной работы эксплуатируемого объекта на некотором заданном интервале времени T. Показатель , исходя из расчета на наихудший случай, можно считать нижней границей среднего времени наступления отказа ТС. Соответственно, полученное неравенство можно использовать для оценки вероятности безотказной работы эксплуатируемой ТС на интервале T. Таким образом, легче вычисляемый и более надежный по достоверности (как соответствующий реалиям информационного обеспечения эксплуатируемых ТС) показатель может быть использован для решения задач управления эксплуатацией ТС. Его использование целесообразно на этапе эксплуатационной активности ТС, когда отсутствует возможность выяснения значений таких показателей надежности ТС, как вероятность безотказной работы . Для сравнения пронормируем показатель на интервал [0,1], при этом оказывается возможным интерпретировать нормированный показатель , после подстановки найденных значений оценок неизвестных параметров (после решения задач (10)) интерпретировать как интегральную выборочную функцию распределения времени безотказной работы ТС. Далее, можно использовать один из статистических критериев согласия, в частности критерий А. Н. Колмогорова, для подбора некоторого теоретического распределения времени безотказной работы ТС. Подобранное распределение дает возможность находить значения показателя вероятности безотказной работы ТС на заданном интервале времени. Использование алгоритма определения остаточного эксплуатационного ресурса для оценки состояния судовых ЧМС Рассмотрим примеры приложений вышеизложенного подхода к определению показателя элементов судовых человекомашинных систем. Пример 1. Износ дейдвудных подшипников может быть оценен путем измерений величины зазора и проведения расчетов. Погрешность измерений 3 %. Срок службы подшипника гребного вала T = 20 000 часов, допуск износа составляет мм. В табл. 1 приведены результаты измерений зазора, а в табл. 2 – прогноз состояния подшипника. Таблица 1 Результаты измерения зазора в подшипнике Время, t×103, чЗазор, мм 0,50,1 1,00,18 2,00,32 3,00,5 5,00,8 6,00,96 7,01,2 8,01,36 9,01,54 10,01,8 Таблица 2 Результаты прогноза изменения зазора в подшипнике Время, t×103, чmin , ммmax ммСреднее значение мм 12,02,052 4762 2488 13,02,132 80082 4653 14,02,173 18292 6765 15,022 2453 59472 9096 Пример 2. Техническое состояние теплообменных аппаратов, применяемых для охлаждения двигателей, можно оценить по температуре охлаждающей жидкости, измеряемой с погрешностью ±2 %. В табл. 3 приведены данные температурного контроля, а в табл. 4 – прогноза состояния системы охлаждения судового дизеля. Срок службы такой системы 10 000 ч, диапазон рабочих значений температуры от 40 до 80 °C. Таблица 3 Результаты измерения температуры охлаждающей жидкости двигателя Время, t× 103, чТемпература, t, °C 1,058,0 2,061,1 4,064,8 5,067,0 Таблица 4 Прогноз температуры охлаждающей жидкости двигателя Время, t × 103, чmin t, °Cmax t, °CСреднее значение t, °C 6,066,579,6373,07 7,067,1487,077,07 8,067,8395,581,67 9,068,56105,1386,84 10,069,3115,992,6 Пример 3. Показатель оптической плотности миокарда у экипажа судна определялся методом ультразвуковой денситометрии с погрешностью ±10 %. В табл. 5 приведены данные оптической плотности миокарда межжелудочковой перегородки (МЖП), а в табл. 6 – прогноза состояния оптической плотности МЖП. Возраст обследуемого электромеханика на момент начала обследования – 32 года. Частота обследования – 1 раз в два года. Всего проведено 4 обследования. Таблица 5 Результаты измерения оптической плотности МЖП Возраст, летОптическая плотность МЖП, ед. 3011,7 3211,7 3411,9 3611,95 Таблица 6 Результаты прогноза оптической плотности МЖП Возраст, летmin y(t), ед.max y(t), ед.Среднее значение y(t), ед. 3811,8612,412,13 4011,8712,8512,36 4211,8813,412,64 Заключение Анализируя результаты прогнозов, следует отметить, что применение предложенного подхода позволяет гарантированно устанавливать возможность безотказного функционирования ЧМС на основе малого количества результатов контроля состояния (желательно одновременного для ТС и ЧС).