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

А.Г. Галкин - Надежность и диагностика систем

.pdf
Скачиваний:
485
Добавлен:
18.03.2016
Размер:
2.93 Mб
Скачать

1.2. НЕВОССТАНАВЛИВАЕМЫЕ ОБЪЕКТЫ

 

 

x0

 

P xп x0 fп x dx.

 

 

0

 

Нагрузка и прочность - случайные величины

 

fH(x)

fH(x)

 

fП(x)

 

 

 

fП(x)

 

Н

П

x

а) общий вид модели

 

 

fH(x)

fH(x)

fП(x)

fП(x)

fH(x0)

 

fП(x0)

 

x0

x

dx

б) увеличенный масштаб области перекрытия

Рис. 1.1

Но нагрузка xн может меняться в интервале (0, ). Поэтому для определения вероятности отказа необходимо проинтегрировать полученное

выражение на этом интервале, заменив x0 на x

21

1.2.2. Модель отказа: нагрузка и прочность - случайные величины

Q

 

x

 

(1.9)

fн

(x)

(y)dy dx.

 

fп

 

 

 

 

0

 

 

Справедливо и другое уравнение

Q

 

 

 

(1.10)

fп

(x)

(y) dy dx.

 

fн

 

 

 

0

x

 

 

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

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

Прочность и нагрузка имеют распределение Гаусса (нормальное)

 

 

 

 

 

 

 

1

 

 

 

z2

 

 

 

 

 

 

Q 1

 

 

 

 

exp(

 

)dz,

(1.11)

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

a

 

 

 

 

п н

 

 

 

 

 

 

 

 

 

где a =

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п2 н2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п и

 

н -

 

соответственно

параметры центра

распределения

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

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

Вероятность отказа можно выразить через нормальное распределение нормированной случайной величины.

 

 

 

 

 

 

 

 

 

 

 

 

 

п н

 

 

 

 

 

 

 

 

 

 

 

Q Фн

 

 

 

 

,

(1.12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п2 н2

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

2

 

 

 

 

 

 

 

 

 

где

Фн(x)

 

 

 

 

exp(

 

z

 

)dz .

 

 

 

 

 

(1.13)

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

22

1.2. НЕВОССТАНАВЛИВАЕМЫЕ ОБЪЕКТЫ

Прочность и нагрузка имеют экспоненциальное распределение

Q 1

н

,

(1.14)

 

п н

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

Прочность имеет распределение Гаусса, а нагрузка - экспоненциальное

 

 

 

 

п

 

 

 

 

1

 

 

 

 

 

 

 

 

 

2

 

2)

 

Q

 

 

 

 

exp

 

(2

п

 

 

н

 

 

 

 

 

 

 

 

 

Фн

 

п

 

 

 

 

2

 

 

 

 

 

н

 

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.15)

 

 

 

 

 

 

 

 

 

п

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 Фн

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где п и п - соответственно параметр центра (математическое ожидание) и параметр масштаба (среднее квадратическое отклонение) прочности;

п - параметр распределения нагрузки.

Фн(x) - здесь и далее определяется по формуле (1.13).

Прочность и нагрузка имеют гамма - распределение

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

Г(

п

 

н

)

 

(1 r)

 

 

1

 

 

 

1

 

 

Q 1

 

 

 

 

 

1-U

п

 

U

 

н

 

dU,

(1.16)

 

 

 

 

 

 

 

 

 

Г( п) Г( н)

0

 

 

 

 

 

 

 

 

 

где Г(x) - гамма функция, см. прил. 2;

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

r н п ;

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

Нагрузка имеет распределение Гаусса, прочность - Вейбулла

23

1.2.2. Модель отказа: нагрузка и прочность - случайные величины

 

 

 

 

 

 

 

 

x

0п н

 

 

 

 

 

1

 

 

п x0п

 

 

 

 

 

Q 1 Ф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н

 

н

 

 

 

 

 

 

2

 

 

 

 

н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.17)

 

 

 

 

 

 

 

 

1

 

x

 

 

 

 

 

 

x

 

 

 

 

 

2

 

 

 

 

 

 

 

0п

 

 

 

0п

н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp

 

y

 

п

 

 

 

 

п

 

 

y

 

 

 

 

dy,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

н

 

 

 

 

 

 

 

 

 

н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где x0п, п, п - соответственно параметры сдвига, формы и масштаба распределения прочности;

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

Прочность и нагрузка имеют распределение Вейбулла

 

 

 

 

п

1

 

x0п x

 

н

 

 

 

 

 

 

 

 

 

Q

 

exp( y) exp

 

 

y

п (

 

)

 

dy,

(1.18)

 

 

 

 

 

 

 

 

 

 

 

 

 

0

н

 

 

н

 

 

 

 

 

 

 

 

 

 

 

 

 

где x0п, п,

п - соответственно параметры сдвига, формы масштаба

распределения прочности;

 

 

 

 

 

 

 

x, н,

 

н - соответственно параметры центра и масштаба

распределения нагрузки.

 

 

 

 

 

 

 

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

KЗ

 

п

.

(1.19)

 

 

 

 

н

 

 

 

 

 

Покажем, как влияют на вероятность отказа изменения математических ожиданий и коэффициентов вариации (см. прил. 2) распределений. Для простоты будем считать параметры распределения нагрузки неизменными.

Пусть имеются некоторые плотности распределений нагрузки fн(x) и

прочности fп(x), показанные на рис. 1.2, а). При уменьшении математического ожидания прочности (Рис. 1.2, б) площадь перекрытия возрастает, а следовательно, возрастает и вероятность отказа. Фактически это означает, что снижение коэффициента запаса приводит к росту вероятности отказа. Теперь допустим, что прочность имеет такое же математическое ожидание, как на рис. 1.2, а), но коэффициент вариации выше. Эта ситуация показана на рис. 1.2, в). Коэффициент запаса прочности не изменился, но

24

1.2. НЕВОССТАНАВЛИВАЕМЫЕ ОБЪЕКТЫ

вероятность отказа возросла. Это произошло из-за того, что больший коэффициент вариации привел к большему перекрытию распределений. Указанные эффекты наглядно иллюстрирует рис 1.2, г), где показана зависимость вероятности отказа от коэффициента запаса при различных коэффициентах вариации, для случая, когда нагрузка и прочность имеют распределение Гаусса. Из рис. 1.2 можно сделать важный для практики вывод о значении высоких технологий. Высокие технологии позволяют изготовлять объекты с большой степенью точности, и, следовательно, с низким значением коэффициента вариации прочности. Это позволяет значительно снизить коэффициент запаса не увеличивая вероятность отказов и не снижая надежности объектов. Снижение коэффициентов запаса, в свою очередь, позволяет снижать массу, габариты и стоимость объектов.

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

fн(x)

 

Кза

fн(x)

f (x)

 

на

f

(x)

п

 

па

п

 

 

 

 

 

 

 

Qа

 

 

н

п

x

 

 

 

а)

 

 

 

fн(x)

 

Кзв = Кза

 

 

fп(x)

 

нв = на

 

 

 

 

пв > па

 

 

 

 

Qв > Qа

 

 

н п x

в)

Кзб < Кзанб = напб = па

Qб > Qа

 

н

п

 

x

 

 

б)

 

 

Q

 

 

 

 

0,4

 

 

νПН

 

0,3

 

 

 

0,2

 

 

 

νП=0,5

0,1

ν =0,01

 

νП=0,3

 

П

 

 

 

 

1,2 1,4

1,6

1,8 2,0

КЗ

 

 

г)

 

 

Рис. 1.2

Пример:

Электрическая прочность изоляторов электроустановки аппроксимирована распределением Гаусса с параметром центра п = 30 кВ и

25

1.2.2. Модель отказа: нагрузка и прочность - случайные величины

параметром рассеяния п = 2,6 кВ. Коммутационные перенапряжения также распределены по закону Гаусса с параметром центра н = 20 кВ и параметром рассеяния н = 5 кВ.

Найти:

1.Коэффициент запаса прочности.

2.Коэффициенты вариаций прочности и нагрузки.

3.Вероятность отказа (пробоя изоляции электроустановки коммутационными перенапряжениями).

Решение:

1) Коэффициент запаса прочности найдем по формуле (1.19)

30 KЗ20 1,5

2) Коэффициент вариации прочности

п п

2,6п 30 0,086

Коэффициент вариации нагрузки

5н 20 0,25

3) Вероятность отказа по формуле (1.12)

Q Фн(

 

30 20

 

) 3,6 10 2

 

 

 

2,62 52

 

 

 

 

Программа на Бейсике для расчета вероятности отказа, если нагрузка и прочность имеют распределение Гаусса, приведена в прил. 6.

26

1.2.НЕВОССТАНАВЛИВАЕМЫЕ ОБЪЕКТЫ

1.2.3.Модель отказа: нагрузка и прочность - случайные

процессы

В реальных условиях нагрузка и прочность часто могут изменяться со временем. В этом случае можно говорить о случайных процессах (прил. 4) изменения нагрузки и прочности. На рис. 1.3 показаны три характерных

ситуации. Прочность и нагрузку будем считать случайными процессами xп(t)

и xн(t). Под действием нагрузки xн(t) прочность непрерывно уменьшается до тех пор, пока объект не выйдет из строя. Время безотказной работы объекта определяется функцией распределения Q(t). Функция Q(t) зависит от значений прочности и нагрузки, значений параметров их законов распределений, а также от взаимосвязи функций xп(t) и xн(t).

Для получения аналитических зависимостей вероятности отказа объекта Q(t) и других показателей надежности выразим нестационарные

случайные функции xп(t) и xн(t) через более простые типы случайных функций. Для модели на рис. 1.3, а) нагрузка - стационарный случайный процесс с параметрами

t ,

t .

Прочность - нестационарный случайный процесс с монотонно убывающим математическим ожиданием

xп t xп t ,

(1.20)

где xп(t) - случайная величина;

п(t) - неслучайная функция усталости. Тогда п(t) = п(0) (t). Модель (вид) функции (t) зависит от характера физико-химических

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

t exp kt .

(1.21)

Свойства функции усталости:

t 0 1;

lim t 0. t

27

1.2.3. Модель отказа: нагрузка и прочность - случайные процессы

Модель отказа: нагрузка и прочность - случайные процессы

 

Xн(t)

п(0)

Xп(t)

 

 

Xп(t)

 

 

 

 

 

 

fп(x, t1)

 

 

 

 

п(t)

 

 

 

 

 

н(t)

 

 

 

 

Xн(t)

 

 

 

fн(x, t1)

 

 

 

t1

а)

t

Xн(t)

п(0)

 

Xн(t)

п(0)

Xп(t)

Xп(t)

Xп(t)

Xн(t)

 

 

 

 

 

 

п(t)

 

п(t)

 

 

 

 

 

 

Xн(t)

 

Xп(t)

 

н(t)

 

 

(t)

 

 

 

 

н

 

б)

 

t

t

 

 

 

в)

а) нагрузка - стационарный процесс: прочность - нестационарный процесс;

б) нестационарная возрастающая нагрузка;

 

в) нестационарная убывающая нагрузка.

 

Рис. 1.3

В работе Острейковского В.А. приводится ряд моделей физикохимических процессов определяющих деградационные явления.

28

1.2. НЕВОССТАНАВЛИВАЕМЫЕ ОБЪЕКТЫ

Химические реакции

Рост трещин в твердых телах

Изменение относительного удлинения

Старение полупроводниковых ферромагнитных материалов (снижение магнитной проницаемости)

Старение конденсаторов

Старение терморезисторов

Скорость коррозии

Износ

Изменение переходного сопротивления разборных электрических контактов

C C0exp kP t ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ea

 

 

 

 

 

 

 

 

 

 

 

 

 

kP

kP0 exp

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R1T

 

 

 

 

 

 

 

 

 

 

 

 

dl

A kn1;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dN

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

CP

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

b t exp ln bKP

 

0,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

M

 

 

 

 

 

 

t

 

 

 

M

 

l

 

 

 

 

1 exp

 

 

 

 

 

 

 

 

 

t.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

E2

 

 

 

 

 

 

 

 

1 exp utt max 1 exp ut .

JУТ t JУТmax 1 exp( k1 tn1 .

R t Rmax 1 exp( k1 tn1 .

dy

dt c0KP exp( 1T);

dy c0 exp T t 1. dt

x t a2 exp(k2t) 1 , при k2 2; x t a2 exp(k2t) , при k2 0.

Rп t

 

 

 

 

 

Rп0

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

Ea

 

1

 

2 D

0

t exp

 

 

 

R T

 

 

a

0

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

29

1.2.3. Модель отказа: нагрузка и прочность - случайные процессы

Напряжение срабатывания

UCP t UCP0 exp( a tb).

электромагнитного реле

 

Здесь обозначены:

 

t

- время;

 

T

- абсолютная температура;

D0

- коэффициент диффузии при Т = 0 С;

Ea

- энергия активации;

 

R1

- универсальная газовая постоянная;

С- концентрация вещества;

kP

- константа скорости химической реакции;

Э0

- удельная электропроводность при Т = 0 С;

T, Eп, E - температурные коэффициенты;

E

- напряженность электрического поля;

eп0

- начальная диэлектрическая проницаемость при Т = 0 С;

EKP0

- напряженность электрического поля при пробое для Т=0 С;

M

- механическое напряжение;

EпP0

- напряженность электрического поля при пробое для M = 0;

2l

- длина трещины;

N

- число циклов нагружения;

k

- размах коэффициента интенсивности напряжения;

bKP

- критическая глубина трещины;

CP

- среднее механическое напряжение;

0,2

- предел текучести;

e

- деформация;

E1

- модуль упругости;

E2

- модуль упругости при запаздывающей упругой деформации;

 

- постоянная времени;

 

- магнитная проницаемость;

J

- ток утечки;

R

- изменение сопротивления резистора;

m

- масса;

 

- угловая скорость;

e3

- начальный эксцентриситет;

a0

- радиус пятна касания;

U

- напряжение;

30