Надежность трубобуров

Автор работы: Пользователь скрыл имя, 15 Декабря 2011 в 13:52, курсовая работа

Краткое описание

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

Содержимое работы - 1 файл

курсач%20фарахова[1] (восстановлен).doc

— 890.50 Кб (Скачать файл)

Определим математическое ожидание: 

                      

 

Определим среднеквадратичное отклонение: 

                       

 

Определим значение дисперсии: 

                   

 

Определим значение коэффициента вариации: 

               

 

Проверим член 124 и 127 по критерию Романовского по формуле: 

            

Проверяем t=127: 

                Проверяем t=124: 

 
 

Следовательно член 124 и 127 исключаем из дальнейшего рассмотрения. 

3.3 Проверка резко выделяющихся величин по критерию Ирвина по формуле: 

                         

      

     Следовательно,  анализируемые величины оставляем при дальнейшем рассмотрении.

 

3.4 Проверка резко выделяющихся величин по критерию Груббса по формуле: 

                            

     

     Так как для обеих точек при n=189 заведомо (таблица 5 приложения), то оставляем крайние точки в рассматриваемой совокупности.

      Таблица №7 – расчетная таблица

    Интервал,

          ч

    ∆t Середина n*i       p*i
    1 0-10 10 5 23 0,121693
    2 10-20 10 15 20 0,10582
    3 20-30 10 25 19 0,100529
    4 30-40 10 35 20 0,10582
    5 40-50 10 45 20 0,10582
    6 50-60 10 55 17 0,089947
    7 60-70 10 65 15 0,079365
    8 70-80 10 75 23 0,121693
    9 80-90 10 85 8 0,042328
    10 90-100 10 95 10 0,05291
    11 100-110 10 105 7 0,037037
    12 110-120 10 115 7 0,037037
 
 

Определим частность в i-ом интервале: 

                 

                                

 

     4 Построение графиков статистических функций(гистограмм)

4.1Определим значение интегральной функции распределения: 

              

 
 

4.2Определим эмпирическую плотность распределения: 

              

    

 

4.3Определим значение обратной интегральной функции распределения: 

             

 
 

4.4Определим значение функции интенсивности: 

             

 

Сведем полученные данные в таблицу:

Таблица №8

      
N ∆t ticp n*i
1 10 5 23 0,121693 0,012169 0,121693 0,878307 0,013855
2 10 15 20 0,10582 0,010582 0,227513 0,772487 0,013699
3 10 25 19 0,100529 0,010053 0,328042 0,671958 0,014961
4 10 35 20 0,10582 0,010582 0,433862 0,566138 0,018692
5 10 45 20 0,10582 0,010582 0,539683 0,460317 0,022989
6 10 55 17 0,089947 0,008995 0,62963 0,37037 0,024286
7 10 65 15 0,079365 0,007937 0,708995 0,291005 0,027273
8 10 75 23 0,121693 0,012169 0,830688 0,169312 0,071875
9 10 85 8 0,042328 0,004233 0,873016 0,126984 0,033333
10 10 95 10 0,05291 0,005291 0,925926 0,074074 0,071429
11 10 105 7 0,037037 0,003704 0,962963 0,037037 0,1
12 10 115 7 0,037037 0,003704 1 0
 

Рисунок №1 –  гистограмма эмпирической плотности  распределения

 

 

Рисунок №2 – гистограмма интегральной функции распределения

 

Рисунок №3 – гистограмма обратной интегральной функции 

 

Рисунок №4 – гистограмма функции интенсивности

 

5. Выбор теоретического закона распределения. 

При обработке  статистического материала важной задачей является подбор теоретического

закона распределения, наилучшим образом описывающего статистическое распределение. 

Теоретический закон подбирают, принимая во внимание: 

- физическую  природу отказов;

- опыт отработки  деталей и изделий аналогичного назначения;

- форму кривой  плотности распределения;

- совпадение  опытных точек с теоретической  кривой интегральной функции  безотказности;

- коэффициент  вариации. 

Значение коэффициента вариации, характеризующего рассеивание  показателя надежности,

уже позволяет  судить об условиях эксплуатации машин  и их технологии изготовления. 

5.1 Определим  коэффициент вариации: 

                  
 

Так как V=0,64, то выбираем закон распределения Вейбулла. 

 
 
 
 

5.2 Определим параметры распределения Вейбулла. 

Так как коэффициент  вариации V=0,64, то по таблице 2 приложения определим находим: 

b=1.6 

5.3 Для определения  параметра a нужно найти Кb. 

                 

По таблице  находим

0,640=V;

0,897=Kb;

 

Следовательно уравнения закона распределения  Вейбулла примут вид: 

 

 

     6.Рассчитаем  все теоретические функции.

6.1 определим значение интегральной функции распределения:

             

 

 

6.2 Определим дифференциальную функцию распределения: 

        

 

 

  6.3 Расчет обратной интегральной функции распределения: 

                   

 
 
 
 
 

 

6.4  Расчет значение функции интенсивности отказов: 

                    

 

 

 Найдем  разность между функциями F*(t) и F(t): 

        

 
 

 Определим  вероятность попадания случайной  величины в i-ый интервал pi: 

        

 
 

Сведем полученные данные в таблицу:

 

Таблица №6

t f(t) P(t) pi F*(t) F(t) λ(t)
0   0 1     0 0  
10 23 0,001795 0,936 0,0635 0,121693 0,064 0,001917 0,058
20 20 0,004762 0,819 0,1169 0,227513 0,181 0,005812 0,047
30 19 0,007596 0,683 0,1362 0,328042 0,317 0,011119 0,011
40 20 0,009633 0,547 0,1364 0,433862 0,453 0,017618 0,019
50 20 0,010624 0,422 0,1247 0,539683 0,578 0,025178 0,038
60 17 0,010619 0,315 0,1069 0,62963 0,685 0,033707 0,055
70 15 0,009837 0,228 0,086 0,708995 0,772 0,043135 0,063
80 23 0,008565 0,160 0,0676 0,830688 0,840 0,053409 0,009
90 8 0,007075 0,110 0,05 0,873016 0,890 0,064485 0,017
100 10 0,005581 0,073 0,0365 0,925926 0,927 0,076326 0,001
110 7 0,004225 0,048 0,0255 0,962963 0,952 0,0889 0,010
120 7 0,003081 0,030 0,0173 1 0,970 0,102179 0,030
 
 

      f(t) - дифференциальная функция распределения (эмпирическая плотность распределения)

    P(t) - обратная интегральной функции распределения

    pi - вероятность попадания случайной величины в i-ый интервал

    F(t) - интегральная функция распределения

    λ(t) - функция интенсивности отказов

    |D| - разность между функциями F*(t) и F(t) 
    Рисунок №5 – гистограмма и график теоретической эмпирической плотности

      Рисунок №6 – гистограмма и график теоретической интегральной функции

 

Рисунок №7 – гистограмма и график теоретической обратной интегральной функции распределения

Рисунок №8 – гистограмма и график теоретической функции интенсивности отказов

 

7. Проверка гипотезы о соответствии эмпирического и теоретического распределений с помощью критериев согласия. 

Назначение критерия χ2 - критерия Пирсона 

Критерий χ2 применяется  в двух целях:  

1) для сопоставления  эмпирического распределения признака  с теоретическим - равномерным, нормальным или каким-то иным;  

2) для сопоставления  двух, трех или более эмпирических  распределений одного и того  же признака (в скрипте до 10).  

Описание критерия 

Критерий χ2 отвечает на вопрос о том, с одинаковой ли частотой встречаются разные значения признака в эмпирическом и теоретическом распределениях или в двух и более эмпирических распределениях.  

Преимущество  метода состоит в том, что он позволяет  сопоставлять распределения признаков, представленных в любой шкале, начиная от шкалы наименований. В самом простом случае альтернативного распределения "да - нет", "допустил брак - не допустил брака", "решил задачу - не решил задачу" и т. п. мы уже можем применить критерий χ2.  

Информация о работе Надежность трубобуров