Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

1658

.pdf
Скачиваний:
7
Добавлен:
07.01.2021
Размер:
1.51 Mб
Скачать

Современные методы прогнозирования по характеру исходной информа-

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

ских данных, полученных в результате предварительных исследований, экс-

пертных оценок, моделировании, включающем физические, физико-

математические и информационные модели.

В зависимости от процедуры прогнозирования износа различают:

– статистическое прогнозирование (по данным о степени рассеивания из-

носа в определенные, т.е. контрольные промежутки времени);

– прогнозирование по реализации (по результатам испытаний одного или нескольких элементов на изнашивание для стационарных процессов, обладаю-

щих свойством эргодичности);

– расчетные методы (по аналитическим выражениям, описывающим про-

цесс изнашивания);

–комбинированные методы (прогнозирование по реализации с учетом ста-

тистических характеристик; расчет износа на основания экспериментальных за-

висимостей и др.).

Наибольшее распространение в настоящее время получили методы про-

гнозирования, основанные на экстраполяции значений прогнозируемого пара-

метра по кривой изнашивания. Процедуру экстраполяции в математике пони-

мают следующим образом. Если известно значение функции в точках х01<…<хn лежащих внутри интервала [x0;xn], то процедуру установления зна-

чения функция f(х) в точках х, лежащих вне интервала [x0;xn], называют экстра-

поляцией. Методы прогнозирования, основанные на экстраполяции, в зависимо-

сти от процедуры построения аппроксимирующей кривой можно разделить на:

прямую экстраполяцию; адаптивные методы экстраполяции; корреляционный метод экстраполяции; метод огибающих кривых; параметрические методы.

Обычно используют метод прямой экстраполяции.

Процедура прогнозирования состоит из анализа исходных данных и по-

строения графика, иллюстрирующего изменения прогнозируемого параметра

(износа, зазора или размеров детали) во времени; определения аналитического

выражения (математической модели), описывающего закономерность измене-

ния прогнозируемого параметра во времени; экстраполяции полученного урав-

нения и прогнозирования изменения параметра на заданный период.

После построения графиков по результатам испытаний или исследований в эксплуатации, которые отражают связь между переменными, подбирают ана-

литическую функцию. Подбор функции составляет важную часть прогнозиро-

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

лом переменных (таблица 6).

Параметры эмпирических формул определяют или методом выбранных точек, или методом средних, или методом наименьших квадратов.

Процедуру прогнозирования износа методом экстраполяции эксперимен-

тальной зависимости можно представить следующим образом.

Износ на прогнозируемый момент времени J( пр), определяется функци-

ей

J( пр) ( 1 ),

(13)

где 1 – период наблюдений или продолжительность испытаний; – прогнози-

руемый период, (0,3–0,5) 1.

Среднее квадратичное отклонение рассчитать по формуле

 

n

 

 

 

 

 

 

yi

yi 2

 

 

i 1

 

,

(14)

 

 

 

 

n

 

где уi – значение износа при выборке (см. таблицу 2); yi – значение износа по теоретическому динамическому ряду; n – количество интервалов наблюдения.

Пример: Рассчитаем износ втулок двигателя дорожной машины на про-

гнозируемый период, равный трем интервалам времени.

В результате экспериментальных исследований получено, что интервалу

времени = 1, 2, 3, 4, 5 соответствует износ J = (7, 10, 12, 15, 20) 10-5 м.

1.По экспериментальным данным строим график функции y = f(x).

2.Выбираем (линейную) модель прогнозирования y a0 a1 x.

3.По формулам (см. таблицу 6).

 

n

n

 

n

 

xi yi n xi yi

a1

i 1

i 1

 

i 1

 

n

2

n

 

 

 

xi

n xi2

 

i 1

 

i 1

 

 

1

n

n

 

a0

 

 

yi a1 xi

n

 

 

i 1

i 1

 

определяем параметры

а0 = 4,54 10-5 м, а1 = 2,6 10-5 м.

4.Линейная модель прогноза имеет следующий вид: y 4,54 10 5 2,6 10 5 x.

5.Вычисляем теоретический динамический ряд:

y1 4,54 10 5 2,6 10 5 1 7,14 10 5 м,

y2 4,54 10 5 2,6 10 5 2 9,74 10 5 м,

y3 4,54 10 5 2,6 10 5 3 12,34 10 5 м,

y4 4,54 10 5 2,6 10 5 4 14,94 10 5 м,

y5 4,54 10 5 2,6 10 5 5 20,14 10 5 м.

6. Прогноз вперед:

на один интервал времени

y6 4,54 10 5 2,6 10 5 6 22,74 10 5 м,

на два интервала

y7 4,54 10 5 2,6 10 5 7 25,34 10 5 м,

на три интервала

y8 4,54 10 5 2,6 10 5 8 27,94 10 5 м.

Таблица 6

Основные зависимости при математическом описании изнашивания

Аналитическая

Формулы для определения

 

График функции

 

запись

 

 

 

 

параметров

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

n

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

xi yi n xi yi

y

 

 

 

a1

i 1

 

i 1

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

n

 

 

 

 

 

 

y a0 a1 x

 

 

 

 

 

xi

 

 

n xi2

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

n

 

 

 

 

n

 

 

 

 

 

 

a

 

 

 

 

 

 

y

 

a

 

x

 

 

 

 

 

 

0

 

 

i

1

i

 

x

 

 

 

 

 

n

i 1

 

 

i 1

 

 

 

 

 

 

 

n

 

n

 

 

y

 

 

 

 

yi xia1 a0 xi2a1 0

 

y a0

xa

1

i 1

 

i 1

 

 

 

n

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yi xia1 lnxi a0 xi2a1 lnxi 0

 

 

 

 

i 1

 

i 1

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

n

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

yi ea1

xi

a0 e2a1

xi

0

 

y a0

ea

1 x

i 1

 

i 1

 

 

 

n

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yi xi ea1

xi

a0 xi e2a1 xi

0

 

 

 

 

i 1

 

i 1

 

 

x

 

 

 

 

 

 

 

 

 

Таблица 7

Значения износа (у), полученные в результате эксперимента

Вариант

 

 

 

 

Интервал

 

 

 

 

 

 

 

 

 

 

 

 

 

х=1

х=2

х=3

х=4

х=5

х=6

х=7

х=8

х=9

 

 

 

 

 

 

 

 

 

 

 

1

2,4083

4,7912

7,1741

9,5569

11,9398

14,3227

16,7056

19,0884

21,4713

2

1,4068

2,7928

4,1789

5,5649

6,9509

8,3370

9,7231

11,1091

12,4951

3

2,8089

5,5905

8,3721

11,1538

13,9354

16,7169

19,4986

22,2802

25,0618

4

3,4099

6,7896

10,1693

13,5439

16,9200

20,3084

23,6881

27,0779

30,4475

5

2,0077

3,9981

5,9759

7,9601

9,9442

11,9284

13,9125

15,8997

17,8808

6

2,6086

5,1909

7,7731

10,3553

12,9376

15,5198

18,1021

20,6843

23,2665

7

2,0062

4,4762

6,5644

8,7757

10,8281

13,1610

14,9940

17,6511

20,0213

8

4,4762

6,5644

8,7751

10,8281

13,1610

14,9911

17,6511

20,0213

21,6706

9

1,0043

2,0645

2,9784

3,9871

4,8736

5,9758

6,6913

8,0497

9,1760

10

1,8051

3,6773

5,3741

7,1845

8,8483

10,7710

12,2296

13,6881

15,9167

11

2,0059

4,0804

5,9729

7,9839

9,8419

11,8697

13,6141

16,0622

18,2201

12

1,4052

2,8708

4,1762

5,5858

6,8609

8,3734

9,4604

11,2547

12,7936

13

3,2074

6,4995

9,5665

12,7068

15,7305

19,0890

22,4479

25,2285

28,7302

14

4,6090

8,7138

12,4365

16,6160

20,3840

24,4843

28,7943

32,6090

36,3843

15

1,0047

2,0691

2,9691

4,0011

4,8053

5,9061

7,0153

8,1130

9,5270

 

 

 

 

 

 

 

 

 

 

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

n 2

yi yi

 

i 1

 

0,21 10 5

 

 

 

 

n

Таким образом, износ втулки на прогнозируемый период с вероятностью

= 0,95 составит:

Ji Ji

,

(15)

t

 

 

 

,

(16)

 

 

 

 

 

 

 

 

 

 

 

 

n

 

где – точность значений; t – квантиль распределения Стьюдента (таблица 8).

Значение t находим по известным таблицам. Так точность значений

2,78 0,21 10 5 0,261 10 5 5

на один интервал времени

J1 = 22,74 10-5 0,261 10-5 м,

на два интервала

J2 = 25,34 10-5 0,261 10-5 м,

на три интервала

J3 = 27,94 10-5 0,261 10-5 м.

 

 

 

 

 

 

 

Таблица 8

 

 

 

Таблица значений t

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,95

0,99

0,999

 

0,95

0,99

 

0,999

 

 

 

 

 

 

 

 

 

5

2,78

4,60

8,61

20

2,093

2,861

 

3,883

6

2,57

4,03

6,86

25

2,064

2,797

 

3,745

7

2,45

3,71

5,96

30

2,045

2,756

 

3,659

8

2,37

3,50

5,41

35

2,032

2,720

 

3,600

9

2,31

3,36

5,04

40

2,023

2,708

 

3,558

10

2,26

3,25

4,78

45

2,016

2,692

 

3,527

11

2,23

3,17

4,59

50

2,009

2,679

 

3,502

12

2,20

3,11

4,44

60

2,001

2,662

 

3,464

13

2,18

3,06

4,32

70

1,996

2,649

 

3,439

14

2,16

3,01

4,22

80

1,001

2,640

 

3,418

15

2,15

2,98

4,14

90

1,987

2,633

 

3,403

16

2,13

2,95

4,07

100

1,984

2,627

 

3,392

17

2,12

2,92

4,02

120

1,980

2,617

 

3,374

18

2,11

2,90

3,97

 

1,960

2,576

 

3,291

19

2,10

2,88

3,92

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Лабораторная работа №4

Тема: «Прогнозирование числа отказов строительных и дорожных

машин по результатам их эксплуатации».

Цель работы: Практическое применение элементов теории надежности применительно к дорожным и строительным машинам.

Для эффективного использования строительных и дорожных машин ин- женеру-механику, осуществляющему их эксплуатацию, необходимы навыки по определению таких показателей надежности, как коэффициент готовности, наработка до отказа, вероятность безотказной работы, среднее время восстановления, вероятность восстановления в заданный срок и т.д.

Особенностью эксплуатации строительных и дорожных машин является огромное разнообразие условий их применения, что приводит к различиям в показателях надежности и, как следствие, к существенным различиям в рекомендациях по использованию конкретного вида машин. Так, при разработке заводами и отраслевыми научно-исследовательскими институтами нормативов по комплектованию машин запасными частями проводится усреднение условий их эксплуатации, а это приводит к расхождению норм и реальной потребности. Различие в условиях эксплуатации приводит к существенным вариациям коэффициента готовности для машин одного и того же типа.

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

Практикум предусматривает выполнение задания: по заданным выборкам наработок структурных единиц машины до отказа осуществить прогноз числа отказов за конкретный промежуток времени.

Данные методические указания позволяют освоить и закрепить навыки по определению и практическому применению элементов теории надежности при-

менительно к дорожным и строительным машинам.

Исходные данные

Для выполнения поставленной задачи проводят контрольные испытания конкретного вида техники в соответствии с планом (N, R, Т). На испытания ста-

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

При определении показателей надежности проводят расчленение машины на составные части, отказ одной из которых не влияет на надежность других

(двигатель, коробка перемены передач, гидросистема, передний мост, задний мост и т.д.).

Считаем, что информация, накопленная за период наблюдения Т по отка-

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

полнения поставленной задачи. Результаты испытаний по данным структурным единицам машины представлены в виде выборок наработок до отказа, (таблица

9). Конкретные варианты выборок при выполнении задания выбираются по но-

меру зачетной книжки.

При восстановлении работоспособности гидросистемы и коробки передач дополнительно фиксировалась частота отказов отдельных деталей, составляю-

щих структурные единицы. Результаты наблюдений представлены в процент-

ном отношении к общему числу отказов структурных единиц и приведены в таблице 10.

Порядок выполнения работы:

1) В соответствии с выборкой наработок до отказа построить вариацион-

ный ряд.

2) Рассчитать, а затем построить гистограмму относительных частот рас-

пределения.

Для определения числа равных интервалов k, на которое следует разбить весь диапазон ряда выборки, можно воспользоваться формулой

k 3,3 lgn 1,

(17)

где n – объем выборки.

 

Шаг выборки h, определяется по формуле

 

h

xmax xmin

,

(18)

 

 

k

 

где xmax, xmin – максимальное и минимальное значение выборки.

Относительную частоту для каждого интервала i, рассчитывают по фор-

муле

 

 

 

 

ni

,

(19)

 

i

n

 

где ni – накопленная частота в каждом интервале.

Плотность относительной частоты pi, определяют по формуле

pi

i

.

(20)

 

 

h

 

Построить гистограмму по плотности относительной частоты (pi) и значе-

ниям выборки (hi)

3) Вычислить статистические оценки параметров распределения: выбо-

рочную среднюю, выборочную дисперсию, среднее квадратичное отклонение,

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

Выборочная средняя xв, определяется по формуле

 

 

 

n

 

 

 

 

xi

 

 

x

в M(x)

i 1

,

(21)

 

 

 

 

 

n

 

где М(х) – математическое ожидание.

Выборочная дисперсия Dв, определяется по формуле

n ni xi xв 2

Dв

i 1

 

,

 

 

 

 

n

где ni – число повторений одинаковых чисел в выборке.

Среднее квадратичное отклонение , определяется по формуле

Dв ,

Коэффициент вариации в, определяется по формуле

в M(x) .

(22)

(23)

(24)

4) По виду гистограммы и численному значению коэффициента вариации выдвинуть гипотезу о законе распределения. Следует помнить, что для нор-

мального закона распределения коэффициент вариации не превышает 0,33; при распределении Вейбулла меняется от 0,33 до 1,0 и для экспоненциального рас-

пределения имеет значение, близкое к 1,0.

5) Проверить гипотезу о законе распределения с использованием критерия Пирсона при уровне значимости =0,05 (таблица 13).

При проверке гипотезы в качестве частичных интервалов можно исполь-

зовать интервалы, принятые для построения гистограммы. Если в интервалах частота много меньше объема выборки (<5), то их следует объединить.

Вероятность попадания случайной величины в частичные интервалы Pi,

определяется зависимостью

Pi F(xi 1) F(xi),

(25)

где F(xi 1) – функция распределения при конечном значении интервала (хi+1);

F(xi) – функция распределения при начальном значении интервала (х).

 

Критерий Пирсона 2 , определяется по формуле

 

2

(ni n Pi)

.

(26)

 

 

n Pi

 

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]