Надежность последовательных систем при нормальном распределении нагрузки по однотипным подсистемам

Если рассеяние нагрузки по однотипным подсистемам пренебрежимо мало, а несущие способности элементов независимы друг от друга, то отказы элементов статистически независимы и поэтому вероятность p(R> F0) безотказной работы последовательной системы с несущей способностью R при нагрузке F0 равна произведению вероятностей безотказной работы элементов:

где p(Rj> F0) — вероятность безотказной работы у-го элемента при нагрузке F0; п — число элементов в системе; FKj(F0)— функция распределения несущей способности у'-го элемента при значении случайной величины Rj= F0.

В большинстве случаев нагрузка на конкретные системы имеет существенное рассеяние при применении систем, например одни и те же универсальные машины (станки, автомобили и др.) могут эксплуатироваться в различных условиях. При рассеянии нагрузки по системам оценку вероятности безотказной работы системы p(R > F) в общем случае следует находить по формуле полной вероятности, разбив диапазон рассеяния нагрузки на интервалы ДF, определив для каждого интервала нагрузки произведение вероятности безотказной работы p(Rj> F()у'-го элемента при фиксированной нагрузке на вероятность этой нагрузки f(F,)AF, а затем, просуммировав эти произведения по всем интервалам,

или, переходя к интегрированию,

где f(F) — плотность распределения нагрузки; FRj(F) — функция распределения несущей способности у'-го элемента при значении несущей способности Rt = F.

Расчеты по формуле (4.14) в общем случае трудоемки, так как предполагают численное интегрирование, поэтому при большом п возможны только на ЭВМ. Чтобы не вычислять p(R> F) по формуле (4.14), на практике часто оценивают вероятность безотказной работы систем p(R > Fmax) при нагрузке Fwax максимальной из возможных. Принимают, в частности, Fmax = /лД1 + 3iv), где mFматематическое ожидание нагрузки; vF — ее коэффициент вариации. Это значение Fmax соответствует наибольшему значению нормально распределенной случайной величины F на интервале, равном шести средним квадратическим отклонениям нагрузки. Такой метод оценки надежности существенно занижает расчетный показатель надежности системы.

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

Сравнительные расчеты на ЭВМ по формуле (4.14) (точное решение) и по приведенному далее упрощенному методу показали, что точность упрощенного метода достаточна для инженерных расчетов надежности систем, у которых коэффициент вариации несущей способности не превышает 0,1... 0,15, а число элементов системы не более 10... 15. Суть упрощенного метода заключается в следующем.

1. Задаются двумя значениями FA и FB фиксированных нагрузок. По формуле (4.13) проводят расчет вероятностей безотказной работы системы при этих нагрузках. Нагрузки подбирают с таким расчетом, чтобы при оценке надежности системы вероятность безотказной работы системы получилась в пределах p(R > FA) = = 0,45...0,60 и p(R> FB) = 0,95...0,99, т.е. охватывали бы интересующий интервал. Предварительно ориентировочные значения нагрузок принимают близкими к значениям

  • 2. По табл. 4.1 находят квантили нормального распределения ирА и ирВ, соответствующие найденным вероятностям.
  • 3. Аппроксимируют закон распределения несущей способности системы нормальным распределением с параметрами математического ожидания mR и коэффициента вариации vR. Пусть SR — среднее квадратическое отклонение аппроксимирующего распределения, тогда mR - FA + upASR = 0 и mR - FB + upBSR = 0. Из этих выражений находят выражения для mR и vR = SR/mR.

Таблица 4.1

n

M

?

up = -0,525 0 = 0,700)

ир = -1,281 (/> = 0,900)

ир = -1,645 0? = 0,950)

up = -2,326 (/»= 0,990)

ир = -2,326 (p = 0,995)

1

0,000

1,000

1,000

1,000

1,000

1,000

2

0,545

0,830

0,848

0,857

0,873

0,878

5

1,129

0,684

0,707

0,723

0,752

0,761

10

1,499

0,594

0,632

0,650

0,684

0,684

15

1,694

0,558

0,596

0,642

0,652

0,662

20

1,824

0,533

0,573

0,593

0,630

0,643

30

1,998

0,505

0,544

0,565

0,603

0,618

40

2,116

0,484

0,527

0,550

0,586

0,599

50

2,204

0,469

0,513

0,533

0,574

0,588

60

2,274

0,463

0,503

0,526

0,560

0,581

80

2,371

0,450

0,498

0,516

0,558

0,570

100

2,462

0,436

0,474

0,503

0,541

0,558

150

2,604

0,419

0,457

0,484

0,523

0,538

200

2,701

0,408

0,452

0,474

0,511

0,520

300

2,833

0,394

0,434

0,454

0,493

0,515

400

2,923

0,380

0,427

0,448

0,489

0,507

500

2,9923

0,377

0,420

0,436

0,485

0,496

1000

3,20

0,36

0,40

0,42

0,46

0,47

2000

3,39

0,35

0,38

0,40

0,43

0,44

5000

3,64

0,33

0,36

0,38

0,41

0,42

10000

3,81

0,32

0,35

0,36

0,40

0,41

20000

3,99

0,30

0,33

0,35

0,37

0,38

50000

4,21

0,27

0,32

0,33

0,36

0,37

100000

4,35

0,27

0,32

0,32

0,33

0,34

4. Вероятность безотказной работы системы p(R> F) для случая нормального распределения нагрузки F по системам с параметрами математического ожидания mF и коэффициента вариации vF находят обычным способом по квантилю нормального распределения ир. Квантиль ир вычисляют по формуле, отражающей факт, что разность двух распределенных нормально случайных величин (несущей способности системы и нагрузки) распределена нормально с математическим ожиданием, равным разности их математических ожиданий, и средним квадратическим, равным корню из суммы квадратов их средних квадратических отклонений:

где п — условный запас прочности по средним значениям несущей способности и нагрузки, п = mR / mF.

Использование изложенного метода рассмотрим на примерах.

Пример 4.2. Оценить вероятность безотказной работы одноступенчатого редуктора, если: условные запасы прочности по средним значениям несущей способности и нагрузки составляют: для зубчатой передачи л, = 1,5; для подшипников входного вала Я2 = = 1,4; для подшипников

выходного вала л4 = п5 = 1,6, выходного и входного валов пь = п7 = 2,0. Это соответствует математическим ожиданиям несущей способности элементов mRl=l,5mF; mRi = mR} = l,4mf ; = mRi = ,6mF ;

= mR. = 2mF. Часто в редукторах значения пь и п7 и соответственно и w/?7 бывают больше указанных. Задано, что несущие способности передачи, подшипников и валов нормально распределены с одинаковыми коэффициентами вариации vR{ =vRl = ... = vRl =0,1, а нагрузка по редукторам распределена также нормально с коэффициентом вариации

“0,1.

Решение. Задаемся нагрузками FA и FB. Принимаем FA = 1,3/Иуг, FB = = l,lmyr, предполагая, что эти значения близки к требуемым значениям вероятностей безотказной работы систем при фиксированных нагрузках p(R 5 Fa) и p(R > F„).

Вычисляем квантили нормального распределения всех элементов, соответствующие их вероятностям безотказной работы при нагрузках FA и FB:

Аналогично

По квантилям из табл. 4.1 находим вероятности безотказной работы элементов:

• при нагрузке FA

PxiR * FA) = 0,9099; p2(R > FA) = p}(R > FA) = 0,7627;

PAR* Fa) =ps(R> Fa) = 0,9699; p6(R> FA) = p7(Rz FA) = 0,9998;

• при нагрузке FB

P{R * FB) = 0,9967; p2(R > FB) = Pi(R > FB) = 0,9842;

Pa(R г FB) = ps(R > FB) = 0,9991; p6(R * FB) = p7(R> FB) = 1,0000.

Вероятности безотказной работы редуктора p(R> FA) и p(R > > FB) при фиксированных нагрузках FA и FB оцениваем по формуле (4.1) как произведение вероятностей безотказной работы элементов. Получаем p(R ? Fa) = 0,4979; p(R > FB) = 0,9637. Из сопоставления этих значений с диапазонами допустимых значений, приведенными ранее, следует, что вероятности находятся внутри этих диапазонов. Поэтому изменять значения FА и FB не будем.

По формулам (4.15), (4.16) вычисляем математическое ожидание mR и коэффициент вариации vR несущей способности редуктора:

где upA и upB квантили нормального распределения, соответствующие вероятностям p(R> FA) и p(R> FB).

По формуле (4.17) рассчитываем квантиль ир нормального распределения, соответствующую вероятности p(R> F) безотказной работы редуктора:

По табл. 4.1 находим искомую вероятность, соответствующую полученной квантили:

Пример 4.3. По исходным данным примера 4.2 найти вероятность безотказной работы редуктора по максимальной нагрузке в соответствии с методикой, применявшейся ранее для практических расчетов. Принимаем максимальную нагрузку Fmax =mF (1 + ?>vF )=mF (1 +3 0,1 )= l,3mf.

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

По табл. 4.1 находим соответствующие квантилям вероятности:

Вероятность безотказной работы редуктора при нагрузке Fmax вычисляем по формуле (4.13). Получаем p(R > Fmax) = 0,496.

Из сопоставления результатов решения двух примеров следует, что первое решение дает оценку надежности значительно более близкую к действительной и более высокую, чем во втором примере. Действительное значение вероятности, рассчитанное на ЭВМ по формуле (4.14), равно 0,9774.

 
Посмотреть оригинал
< Пред   СОДЕРЖАНИЕ   ОРИГИНАЛ     След >