Методы исследования надежности в эксплуатации
Так как отказы имеют случайный характер, т. е. момент их наступления не может быть указан заранее, для изучения надежности используются методы теории вероятностей и математической статистики (теории случайных величин). За меру количественной оценки надежности принимается вероятность Р(Т) или P(L) безотказной работы конструкции в течение заданного времени Т или пробега L, т. е. вероятность того, что время t или пробег I исправной работы конструкции будет больше заданного времени Т или пробега L:
Р (Т) = Вер (t>T) или P(L)= Bep (l>L). (3.1)
Практически определение вероятности P(Т) или P(L) основывается на результатах проведения статистических испытаний определенного числа N0 деталей или систем в заданных условиях эксплуатации. Эти испытания заключаются в подсчете числа отказов исследуемых деталей или агрегатов по интервалам времени наблюдений ∆t или пробега . При этом период наблюдения t или пробег I разбивают на интервалы ∆t или ∆l и в каждом из этих интервалов подсчитывается число отказов ni. Надежность определяют как математическую вероятность безотказной работы по формуле:
(3.2)
где N0 — число объектов наблюдения; m = l/∆l-число интервалов пробега, в которых подсчитываются отказы; ni-число отказов в i-ом интервале;
- вероятность отказов конструкции
Расчет надежности по формуле (3.2) предполагает, что независимо от природы и способа оценки (классификации) отказов вагон или троллейбус в целом, и каждый из его элементов может находиться только в двух состояниях: исправном или неисправном. Поэтому сумма вероятностей Р и Q равна единице. Вероятность, равная единице, соответствует достоверному событию, равная нулю - невозможному событию.
В реальных условиях число N0 объектов наблюдения всегда конечно, что исключает возможность непосредственного использования формулы (3.2). Практически вероятности Q и Р определяют по приближенным формулам:
(3.3)
Применение формул (3.3) вместо (3.1) связано с определенной ошибкой, величина которой-тем больше, чем меньше партия N0.
Если при этом выводы исследования распространяются на партию Nr>N0, то вероятность ошибки возрастает еще больше. Поэтому практическим исследованиям надежности должна предшествовать оценка допустимой ошибки и определение на основе ее размера партии N0, которую необходимо поставить под наблюдение. Эту партию называют выборочной совокупностью. Величина ее определяется по формуле теории математической статистики:
(3.4)
где Nr - партия, на которую распространяются выводы исследования (так называемая генеральная совокупность),∆ - допустимая (принятая) величина отклонения (ошибки) выборочной средней от генеральной средней; Р — заданная (ожидаемая) вероятность безотказной работы объектов наблюдения. Если вероятность Р неизвестна даже приблизительно, то принимают Р = 0,5, т. е. исходят из максимально возможной дисперсии (разброса) отказов; k — величина, выбираемая в зависимости от строгости требований к отклонению ∆ вероятности отказа выборки от вероятности отказа генеральной совокупности. Принимают k = 2—3, причем, чем больше k, тем меньше вероятность отклонения.
Принимая, например, Nr= 1000, ∆ = 0,05 (5%), k = 2 и Р = 0,5, получим N0 = 286. При тех же данных, но P = 0,95 соответственно найдем N0 = 70 и т. д.
Распределение отказов во времени или пробеге характеризуется частотой отказов fi = ni/N0 и интенсивностью отказов ,λ(t) или λ(L), которая определяется как число отказов в рассматриваемом интервале времени или пробега, отнесенное к величине этого интервала и числу элементов, оставшихся исправными к его началу:
(3.5)
Характерная гистограмма (ступенчатая кривая) и ее огибающая (кривая распределения) интенсивности отказов λ для элементов, поступающих в эксплуатацию после заводского изготовления или ремонта, показана на рис. 3.1.
На этой кривой видны три характерных периода работы, которым соответствуют три типа отказов: 1 — период приработки (приработочные отказы), 2 — период нормальной эксплуатации (случайные отказы) и 3— период аварийного износа (аварийные или износовые отказы).
Рис.3.1. Гистограмма (ступенчатая кривая) и ее огибающая (кривая распределения) интенсивности отказов λ для элементов, поступающих в эксплуатацию после заводского изготовления или ремонта
К концу периода 2, при пробеге L0 параметры системы достигают предельных значений (исчерпывается эксплуатационный допуск) и в период 3 ее износ будет уже запредельным (аварийным). Элементы с приработочными отказами нормально должны отбраковываться на заводе в процессе обкатки или контрольных испытаний и не выпускаться в эксплуатацию. Наличие их в эксплуатации указывает на недостатки технологии производства и некачественный контроль изделий на заводе-изготовителе или ремонтном заводе. Отказы в период нормального эксплуатационного износа характеризуются примерно постоянной интенсивностью, а длительность этого периода в обычных условиях значительно превышает период приработки. Запредельные (аварийные) отказы периода 3 недопустимы. Точка перехода периода 2 в 3 определяет максимальный пробег подвижного состава до очередного ремонта рассматриваемого элемента.
По характеру кривой распределения интенсивности отказов (λ-характеристике) элемента или системы можно сделать ряд важных выводов — оценок качества заводского изготовления или ремонта и правильности построения системы технического обслуживания ПС. Примеры таких характеристик приведены на рис. 3.2.
Рис.3.2. Примеры λ-характеристик (по результатам исследований отказов различных элементов тяговых двигателей)
Рассматривая, в частности, характеристику, показанную на рис. 3.2, а, видим, что область 1 выражена здесь слабо. Это говорит о том, что детали с приработочными отказами отбракованы на заводе почти полностью. Качество заводского изготовления характеризует и длина интервала L0. Присутствие в кривой области 3 говорит о том, что межремонтные сроки слишком завышены и их нужно снизить, по крайней мере, до величины L0. На рис. 3.2, б область 1 растянута на весь период испытаний. Она говорит о некачественности заводского изготовления и контроля продукции или недостаточной длительности выбранного периода испытаний. Бугристость λ-характеристики, показанной на рис. 3.2, в, указывает на то, что на нее оказывали влияние несколько видов отказов, и на необходимость отдельного анализа каждой из этих причин. Таким образом, но виду λ -характеристики можно судить о качестве изготовления и эксплуатации элементов и разрабатывать мероприятия, направленные на повышение надежности.
Вероятность Р безотказной работы и вероятность отказов Q, как видно из формул (3.1) и (3.2), характеризуются накопленным числом отказов по всем т интервалам периода l или t наблюдения вплоть до рассматриваемого момента L или Т. Кривые зависимости P = F1(L) и Q=F2(L), приведенные на рис. 3.3, называют интегральным законом распределения вероятностей отказов или безотказной работы.
Рис.3.3. Интегральные (а) и дифференциальная (б) кривые распределения вероятности отказов (Q, q) и безотказной работы (P)
Эти кривые показывают вероятность того, что сумма отказов в заданный момент L не превысит Q, а вероятность безотказной работы будет не ниже Р. Дифференцируя интегральную кривую вероятности Q, можно получить так называемый дифференциальный закон или кривую плотности распределения вероятности:
q(L) = dQ(L)/dL = -dP(L)/dL. (3.6)
По кривой плотности распределения вероятностей можно определить пробег Lмакс, при котором можно ожидать максимальное число отказов, а также пробеги, отвечающие разной плотности распределения. Кривая плотности распределения -очень важная характеристика закона распределения. Ее основными параметрами являются:
математическое ожидание (среднее значение)
(3.7)
дисперсия (рассеяние, разброс относительно среднего)
(3.8)
среднеквадратичное отклонение
, (3.9)
где А — В — интервал интегрирования.
Когда зависимости Q=F2(L) и q(L) представляют собой ступенчатые кривые (гистограммы), параметры т и D кривой плотности распределения определяются формулами:
(3.10)
Задачу анализа надежности различных конструкций и схем при наличии всех данных о надежности их элементов решают вполне однозначно. Отказы классифицируют и анализируют по каждой группе, исходя из принципа, что отказ любого элемента приводит к отказу всей системы. Вероятность Рс надежной работы системы, состоящей из n элементов, при таком потоке отказов определяется теоремой умножения вероятностей:
, (3.11)
где Рi — вероятность надежной работы i-го элемента.
Так как каждая из частичных вероятностей Рi меньше единицы, можно сделать вывод: чем больше в системе элементов, т. е. чем она сложнее, тем ниже и ее надежность. Этим объясняется, что подвижной состав со сложной схемой менее надежен, чем подвижной состав с простой схемой. Из формулы (3.11) следует также, что надежность системы всегда меньше надежности самого ненадежного ее элемента.