Надёжность авиационной техники

Тип работы:
Контрольная
Предмет:
Физико-математические науки


Узнать стоимость

Детальная информация о работе

Выдержка из работы

Надёжность авиационной техники

1. Задача 1

Выполнить анализ данных эксплуатационных наблюдений за отказами изделий СКВ парка ЛА, для чего:

— определить вид случайных величин наработки изделий (реализаций);

— построить ранжированную временную диаграмму;

— выбрать размах и число интервалов временной диаграммы.

Решение

Т.к. исходные данные варианта не содержат сведений о наработке изделий «ВВР 4487T» до цензурирования, необходимо использовать только данные наработок изделий до отказа (таблица П. 2. 3). Поэтому в данном случае выполним построение ранжированной временной диаграммы наработок до отказа ?i, i=1, n; (смотри Приложение 1).

Проведённые крайние сечения ранжированной диаграммы левее минимального значения t=1025 ч. и правее максимального значения t=1088 ч. дают значение размаха ?=1900−1000=900 (ч), полученное значение которого разбиваем на L=9 неравных интервалов? ti с сечениями, соответствующими границам интервалов, из которых первые 6 интервалов с шагом 75 ч, а последние 3 — с шагом 150 ч. Правое крайнее сечение диаграммы будем считать также границей периода эксплуатационных наблюдений Т=1900 ч.

2. Задача 2

Выполнить оценку показателей безотказности параметрическим методом для однократно цензурированной выборки, для чего:

а. Выполнить оценку и построение статистической плотности распределения f*(t) и статистической интенсивности отказов ?*(t);

б. Выполнить оценку параметров распределения для однократно цензурированной выборки: Тср. *; а* и б*; mt* и ?t*; (в зависимости от принятого закона распределения наработки до отказа).

в. Выполнить проверку гипотезы о законе распределения для однократно цензурированной выборки;

г. Выполнить оценку показателей безотказности для принятого закона распределения наработки до отказа.

Решение

1. Для оценки показателей безопасности параметрическим методом, выполним построение ранжированной временной диаграммы для однократно цензурированной выборки.

2. Выполним построение диаграмм плотности вероятности наработки до отказа f*(t) и интенсивности отказов ?*(t);

Искомые величины определим по формулам:

f*(t)=?n/(Ni·?t),

?*(t)=?n/([Ni-n (t)]·?t), где:

?ni — число отказавших изделий в интервале? ti;

Ni — число изделий, наблюдаемых в интервале? ti;

ni(t) — число отказавших изделий до начала i-го интервала;

Ni= N-?mi — общее число всех реализаций гистограммы, за исключением неполных реализаций. Результаты расчётов представим в виде таблицы 1 и гистограмм.

Таблица 1.

№ п/п

1

2

3

4

5

6

7

8

9

ti(ч)

1075

1150

1225

1300

1375

1450

1600

1750

1900

?ti(ч)

75

75

75

75

75

75

150

150

150

?ni

1

1

1

3

2

1

1

1

1

Nj

20

19

18

17

14

12

11

10

9

f*(t) (·10-4)

6,667

7,018

7,407

23,529

19,048

11,111

6,061

6,667

7,407

ni(t)

0

1

2

3

6

8

9

10

11

N

20

20

20

20

20

20

20

20

20

mi

8

8

8

8

8

8

8

8

8

Ni

12

12

12

12

12

12

12

12

12

Ni-ni(t)

12

11

10

9

6

4

3

2

1

?*(t) (·10-4)

11,111

12,121

13,333

44,444

44,444

33,333

22,222

33,333

66,667

По результатам расчётов строим гистограммы f*(t) и ?*(t):

f*(t)

(·10-4)

23,529

19,048

11,111

6,667

7,018

7,407

6,061

6,667

7,407

t

Гистограмма плотности вероятности наработки до отказа f*(t)

? *(t) (·10-4)

66,667

44,444

44,444

33,333

33,333

22,222

11,111

12,121

13,333

t

Гистограмма интенсивности отказов ?*(t)

При сравнении полученных гистограмм с теоретическими кривыми f (t) и ?(t) по их виду, предполагаем, что в нашем случае имеется наибольшее сходство с нормальным законом распределения.

Выполним оценку параметров распределения методом максимального правдоподобия.

Т.е. построим функцию правдоподобия L (?, ?), зависящую от результатов наблюдений выборки из N изделий и параметра? неизвестного закона распределения F (?, ?). Найдём оценку ?* при максимальной вероятности наблюдаемого результата Р {?, ?}> max, т. е. построим L (?, ?)=П•Р {?, ?} и определим max из условия? L (?, ?)/??=0 и соответствующее ?*.

Для предполагаемого закона распределения определим математическое ожидание mt и среднее квадратическое отклонение ?t, для чего воспользуемся выражениями:

?t=(T-T1)/[c•?(k) — k] и mt=k•?t+T, где:

T — период наблюдения, T1 — среднее арифметическое выборки наработок до отказа,

T 1=? ti/n, T22=? ti2/n, С= ?(Ni-n)/n•L,

(здесь L-число интервалов группирования L=9).

Используя метод максимума правдоподобия, воспользуемся двумя уравнениями для определения искомых параметров распределения:

Y1=(T2+T12)/(T-T12) и Y2=(1+ck?(k) — [c?(k)]2)/[c?(k) — k]2,

для чего, определив Y1,

приравняем его значение Y2, т. е. Y1= Y2;

В качестве периода наблюдения выберем значение крайней правой границы размаха временной диаграммы ?=1900−1000=900 (ч), т. е. Т=1900 ч. Для определения Y1 определим значения T1, T2 и С, для чего, подставив значения ti, получим:

T1=(1025+1110+1180+1250+1275+1280+1310+1340+1405+1560+1720+1880)/12=1361. 25;

T22=(10252+11102+11802+12502+12752+12802+13102+13402+14052+15602+17202+18802)/12=(1 050 625+1232100+1 392 400+1562500+1 625 625+1638400+1 716 100+1795600+1 974 025+2433600+2 958 400+3534400)/12=1 726 961. 25, откуда T2=1314,139;

Тогда Y1=(1314,139+1361,252)/(1900−1361,252)= 1 854 315. 7015/(-1 851 101. 5625)=-1,0017;

Из условия равенства Y1 и Y2, получим:

Y1= Y2=(1+ck?(k) — [c?(k)] 2)/[c?(k) — k] 2= -1,0017;

Или 1+c•k•?(k) — [c•?(k)]2 =-1,0017•[c•?(k) — k]2=-1,0017•{[c•?(k)]2-2c•k•?(k)+k2}=-1,0017•[c•?(k)]2+2,0034•c•k•?(k) — 1,0017•k2; или, перенося члены полученного уравнения:

1+c•k•?(k) — [c•?(k)]2+1,0017•[c•?(k)]2-2,0034•c•k•?(k)+1,0017•k2=0;

0,0017•[c•?(k)]2-1,0034•c•k•?(k) +1,0017•k2+1=0; Так как k-задаваемое значение, полученное выражение представляет собой квадратное уравнение:

[c•?(k)]2-716,2857•c•k•?(k)+ 715,2857•k2+714,2857=0, подставляя в которое значение с и одно из табличных значений k=1,3 (см. таблицу 2), получим:

[0,2037•?(k)]2-716,2857•0,2037•1,3•?(k)+ 715,2857•1,32+714,2857=0, или, упрощая выражение, получим: [?(k)]2-3790•?(k)+38 470=0, решая которое, получим:

?(k)1; ?(k)2;=(3790±v3790?-4•1•38 470)/2•1=(3790±3769,6445)/2= 3779,8223/10,1778;

С определим, используя данные Таблицы 1:

С=[(20−12)+(19−12)+(18−12)+(17−12)+(14−12)+(12−12)+(11−12)+(10−12)+(9−12)]/12/9=0,2037;

Если округлить значение Y1=-1,0017?-1 и подставить в выражение для Y2, после преобразования получим:

1+ck?(k) — [c?(k)]2 = - [c?(k)]2+2ck?(k) — k2 или k2 — ck?(k)+1=0, откуда: ck?(k)=k2+1 или ?(k)=(k2+1)/ck;

После подстановки в полученное выражение с=0,2037 и k согласно таблицы 2, получим значения ?(k), близкие по значению с вторым корнем, полученным при решении квадратного уравнения. Следовательно мы можем принять допущение для Y1 и воспользоваться упрощённым уравнением.

Тогда ?(k)=(k2+1)/ 0,2037k;

Для определения Кнач вычислим значение:

F0нач)=1−1/?(Ni/n•L)=1−1/(20+19+18+17+14+12+11+10+9)/(12·9)= 0,1692;

Тогда по таблицам 5. 1, 5.2 приложения 5 определим Кнач=1,31;

Задавая значения k, близкие Кнач, построим зависимости ?р(k)=f (К) и ?Т(k)=f0(k)/F0(k);

Для построения воспользуемся Таблицей 2.

Таблица 2.

К

1,30

1,4

1,5

1,6

1,7

1,8

1,9

2,0

f0(k)

0,9032

0,9192

0,9332

0,9452

0,9554

0,9641

0,9713

0,9772

F0(k)

0,1714

0,1497

0,1295

0,1109

0,0940

0,0790

0,0656

0,0540

?р(k)

10,1578

10,3794

10,6366

10,9230

11,2334

11,5638

11,9112

12,2730

?Т(k)

5,2687

6,1402

7,2062

8,5230

10,1638

12,2038

14,8064

18,0963

По полученным данным построим графики зависимости ?р(k)=f (К) и ?Т(k)=f (k) (смотри приложение).

Точка пересечение указанных графиков даёт следующие значения:

?(k)=11,4; k=1,765;

Тогда ?t=(T-T1)/[c•?(k) — k]=(1900−1361. 25)/(0,2037•11,4−1,765)= 538. 75/0. 55 718=966,9227;

Соответственно: mt=k•?t+T=1,765•966,9227+1900=3606,6186;

Выполним проверку гипотезы о законе распределения для однократно цензурированной выборки по критерию Пирсона ? 2, для чего воспользуемся формулой:

? 2= L·?(?ni-N·Рi)/N·Рi, где:

L — число интервалов группирования;

?ni — число наблюдаемых статистических данных, попавших в i-й интервал;

NPi — среднее число данных, попавших в i-й интервал при условии, что гипотеза о законе распределения верна,

Pi = F (ti) — F (ti-1).

Для подтверждения гипотезы о характере закона распределения необходимо соблюдения условия:

?2расч? ?21-o. o1?, где:

?21-o. o1? — левая граница интервала критической области, квантиль ?2 распределения с r=L-1-S степенями свободы, отвечающий вероятности 1−0. 01?, где S — число наложенных связей, зависимых от числа параметров предполагаемого закона распределения наработки до отказа. Значения квантилей определим по таблице приложения 5.

? — принятый уровень значимости в%.

Наиболее употребительные уровни значимости — 1; 5; 10%.

Для нормального закона распределения число независимых условных связей, накладываемых на выбранный закон распределения, S=2. Следовательно r=9−1-2=6;

Расчёт критерия Пирсона сведём в Таблицу 3.

Таблица 3

L

1

2

3

4

5

6

7

8

9

?ni

1

1

1

3

2

1

1

1

1

F (ti) (·10-4)

6,667

7,018

7,407

23,529

19,048

11,111

6,061

6,667

7,407

Pi(·10-4)

6,667

0,351

0,389

16,122

-4,481

-7,937

-5,05

0,606

0,74

NPi

133,34

7,02

7,78

322,44

-89,62

-158,74

-101

12,12

14,8

(?ni-NPi)/NPi

-0,993

-0,858

-0,871

-0,991

1,022

1,006

1,01

-0,917

-0,932

?2расч

-22,716

1−0. 01? (?=1%)

0,99

1−0. 01? (?=5%)

0,95

1−0. 01? (?=10%)

0,9

?21-o. o1?(?=1%)

22,5

?21-o. o1?(?=5%)

16,8

?21-o. o1?(?=10%)

12,6

По данным таблицы можно сделать вывод, что при любом из трёх наиболее употребимых уровней значимости ?, выполняется условие: ?2расч? ?21-o. o1?, подтверждающее гипотезу о характере закона распределения, соответствующем нормальному распределению случайных величин.

3. Задача 3

Выполнить оценку показателей безотказности по полным данным, для чего:

а. Определить число невосстанавливаемых изделий N и число отказавших изделий n. Найти доверительные границы для вероятности безотказной работы Рн и Рв при двухсторонней доверительной вероятности ?=0,95;

б. Обосновать необходимость определения доверительного интервала (нижнюю Рн и верхнюю Рв границы) для доверительной вероятности ?;

Выполнить оценку показателей безотказности по многократно цензурированным выборкам Р (t), Тср. , Т?, для чего:

а. Определить точечные оценки вероятности безотказной работы за непрерывный беспосадочный полёт самолёта t и за период наработки до профилактики t2=300 ч;

б. Вычислить гамма-процентную наработку Т? при ?=85%.

Решение

Оценку вероятности безотказной работы определим по формуле Р=1-q, где q — оценку вероятности отказа получим методом максимального правдоподобия. Функция правдоподобия будет иметь вид: L (q)=CNn•qn•(1-q)N-n; После преобразований функции правдоподобия получены точечные оценки для вероятности отказа при наработке Т:

q*= n/N;

Соответственно для вероятности безотказной работы p*=1 — n/N;

Так как общее число невосстанавливаемых изделий N=20 (в соответствии с объёмом парка и количеством изделий на одном самолёте — согласно таблице П. 2. 9), а число отказавших изделий за период наблюдения n=12, то q*=12/20=0,6 и p*=1−0,6=0,4;

Для более точной оценки вероятностей отказа и безотказной работы изделий определим границы доверительного интервала для генеральной характеристики q, т. е. получим выражение: ?*= Р (qн?q? qв), где ?* - двусторонняя доверительная вероятность, qн и qв — соответственно нижняя и верхняя доверительные границы характеристики.

Для их определения воспользуемся формулами:

qн= n/N•R1, qв= n/N•R2,

т.к. в нашем случае n?0;

Здесь R1 и R2 — коэффициенты, определяемые по таблицам П. 5. 5, П. 5.6 приложения (…).

Определяем R1=1,59 и R2=0,76, тогда qн= n/N•R1=12/20•1,59=0,3774;

qв= n/N•R2=12/20•0,76=0,7895; то есть можно записать ?*= Р (0,3774?q? 0,7895) или окончательно оценка доверительных границ вероятности безотказной работы будет представлена так: Рн=1 — qв=1−0,7895=0,2105; Рв=1 — qн=1−0,3774=0,6226, с доверительной вероятностью ?=0,95.

В связи с отсутствием в задании данных по многократно цензурированным выборкам рассматриваемого изделия, оценку непараметрическим методом по многократно цензурированным выборкам не выполняем.

безотказность самолет вероятность выборка

ПоказатьСвернуть
Заполнить форму текущей работой