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

книги / Надёжность технических систем и техногенный риск. Структурно-энергетическая теория отказов

.pdf
Скачиваний:
3
Добавлен:
12.11.2023
Размер:
1.47 Mб
Скачать

ГЛАВА 7. МЕТОДИКА ОЦЕНКИ НАДЕЖНОСТИ ЭЛЕМЕНТОВ С УЧЕТОМ ИНТЕНСИВНОСТИ ЭНЕРГЕТИЧЕСКОГО ВОЗДЕЙСТВИЯ

7.1. Общие положения

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

n1

It)

i

 

P(t) = exp(−σ It)

(σ

 

.

(7.1)

 

i!

 

i =0

 

 

 

Для определения параметров модели (7.1) необходимо разработать и провести апробирование методики, позволяющей определить параметры n и α . Ранее было отмечено, что указанные параметры являются функциями размеров дефектных структур. Это обстоятельство позволяет прогнозировать вероятность отказов элементов, обеспечивающую заданный уровень надежности, по результатам одной или двух партий элементов с известными размерами дефектных структур. Решение данной задачи вполне возможно с использованием методов неразрушающего контроля.

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

кэнергетическому воздействию.

7.2.Методика определения параметров n и α

1.Отобрать для испытаний партию из N элементов с известными размерами дефектных структур материалов. Для повышения

161

достоверности результатов испытаний целесообразно брать не менее 100 образцов.

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

3.По результатам испытаний определить среднее значение t и

дисперсию σ t2 наработки элементов до отказа:

 

 

 

 

 

1

 

N

 

 

 

 

1

 

 

 

 

 

N

 

 

 

 

 

=

ti

; σ t2=

 

 

 

 

 

(ti

 

)2 ,

(7.2)

 

 

t

t

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

N i=1

 

 

 

 

 

1 i =1

 

где ti – наработка i-го элемента до отказа.

 

4. По формулам

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n =

t 2

;

 

 

 

(7.3)

 

 

 

 

 

 

 

 

 

 

 

σ t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α =

 

 

 

t

 

 

 

 

 

 

(7.4)

 

 

 

 

 

 

 

 

 

 

 

 

I σ

 

t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

определить параметры n и α

модели отказов (6.1).

 

5. Поскольку испытания элементов проводятся при I = const, то

среднее значение

 

, и дисперсию σ 2E энергии возникновения отка-

E

за определить по формулам

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

; σ 2E= I σ2

 

 

t2 .

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

= It

 

 

 

6. По формулам

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n =

E

 

;

 

 

 

(7.5)

 

 

 

 

 

 

 

 

 

 

 

σ 2E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α =

E

σ 2E

 

определить параметры n и α модели (7.1) для I = вающие интенсивность энергетического воздействия.

(7.6)

const, учиты-

162

Если материалы элементов

имеют

порог

чувствительности

к энергетическому воздействию,

то параметры n и α

необходимо

определять по формулам

 

 

 

 

 

 

 

 

 

n =

(

 

t0 )2

; α =

 

 

 

 

t0

 

 

t

t

,

(7.7)

 

 

 

 

 

 

 

t2

 

 

 

σ

t2

 

 

Iσ

 

 

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

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

Отметим, что с увеличением n кривая плотности вероятности отказа f(t) стремится к нормальному закону, поэтому формально можно считать, что имеет место нормальный закон распределения, если отношение

M (t)

=

n α

 

= n > 3,5,

 

 

 

D(t)

n α

2

 

т.е. если n > 12.

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

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

163

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

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

Рассмотрим некоторые возможные случаи зависимости среднего значения и среднего квадратического отклонения энергии возникновения отказа от величины характерного размера l дефектных структур материалов.

1. Зависимости среднего значения E = It и среднего квадратического отклонения σ энергии возникновения отказа задаются выражениями

 

 

 

 

 

 

 

 

l b

 

 

 

 

 

 

 

 

E = It = k1 exp

 

 

 

;

(7.8)

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k2

 

 

l c

 

 

 

σ E=

exp

 

 

 

,

 

(7.9)

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где k1 и k2 – константы для данного материала; L – характерный размер материала элемента.

Подставляя формулы (7.8) и (7.9) в формулы (7.5) и (7.6), получаем следующие выражения для параметров структурноэнергетической модели отказов:

164

 

t

2

 

{k1

exp (l / L)b }2

n =

 

=

 

 

 

σ t2

{k2

exp

(l / L)c }2

 

 

 

 

 

 

 

k1

2

 

 

l b

=

 

 

 

exp 2

 

 

k

 

 

 

2

 

 

 

L

 

 

 

 

 

 

 

 

+l c

; (7.10)

L

 

 

 

 

 

 

 

{k1

exp

(l / L)b }

 

 

k1

 

 

 

t

 

 

 

 

 

α =

 

 

 

=

 

 

 

 

 

 

 

=

exp

I σ

t2

 

{k2

 

 

 

 

c 2

 

Ik22

 

 

 

exp

(l / L)

 

 

 

 

 

 

 

 

 

 

}

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

l +L

l c

2 . (7.11)

L

1. Зависимости среднего значения и среднего квадратического отклонения энергии возникновения отказа от величины характерного размера дефектных структур материалов элементов задаются выражениями

 

 

 

 

 

 

k3

 

 

 

 

 

 

 

 

 

 

 

 

E = It

=

 

;

(7.12)

(l / L)b

 

 

σ t=

 

 

k4

,

 

(7.13)

 

 

 

(l / L)c

 

где k3 и k4 – константы для данного материала.

Тогда параметры n и α структурно-энергетической модели отказов будут определяться по формулам

 

t 2

k3

2

 

l 2c2b

 

n =

 

 

=

 

 

 

 

 

;

σ

2

 

 

 

t

k4

 

 

L

 

 

 

 

 

 

k3

 

l 2cb

α =

 

t

 

=

 

 

 

 

 

 

.

I σ

2

2

 

 

t

Ik4

 

L

(7.14)

(7.15)

2. Зависимость среднего значения энергии возникновения отказа от величины характерного размера дефектных структур материалов элементов задается выражением (7.8), а среднего квадратического отклонения – выражением (7.13).

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

165

 

t 2

k1

2

 

l 2c

n =

 

 

=

 

 

 

 

 

σ

2

 

 

 

t

k4

 

 

L

 

 

l b

 

exp

2

 

 

;

(7.16)

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

k1

 

l 2c

 

 

 

t

 

 

α =

 

 

 

=

 

 

 

 

exp

I σ

2

2

 

 

t

Ik4

 

L

 

 

 

 

 

 

 

 

 

 

 

 

l b

 

 

L

.

(7.17)

 

 

 

 

 

 

 

 

 

 

 

 

 

3. Зависимость среднего значения энергии возникновения отказа от величины характерного размера дефектных структур материалов элементов задается выражением (7.12), а среднего квадратического отклонения – выражением (7.9).

В этом случае выражения для определения параметров n и α будут иметь вид

 

t 2

 

 

 

k3

2

l

2b

l

c

n =

 

 

 

=

 

 

 

 

 

 

 

 

 

exp 2

 

 

 

 

;

σ

2

 

 

 

 

 

 

 

t

 

 

 

k2

 

 

L

 

 

 

L

 

 

 

 

 

 

 

 

 

k3

 

 

l b

 

l c

 

α =

 

 

t

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp 2

 

 

 

.

 

I σ

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

t

 

Ik2

 

 

L

 

 

L

 

(7.18)

(7.19)

Анализ формул (7.10), (7.11) и (7.14)–(7.19) показывает, что в зависимости от значений b и c возможны следующие закономерности изменения параметров структурно-энергетической модели отказов:

1)n и α – монотонно возрастающие функции размеров дефектных (чувствительных) структур материалов элементов при I = const.

2)n и α – монотонно убывающие функции размеров дефект-

ных (чувствительных) структур материалов элементов при I = const. 3) n = const, α – монотонно возрастающая функция размеров дефектных (чувствительных) структур материалов элементов при

I = const.

Эти закономерности получаются, если в формулах (7.10),

(7.11), (7.14), (7.15) взять b = c = d. Тогда получим

 

k1

2

 

k1

 

l d

 

n =

 

= const;

α =

 

exp

 

 

(7.20)

 

2

 

k2

 

 

Ik2

 

L

 

166

 

 

 

 

 

 

 

 

для экспоненциального закона изменения среднего значения и среднего квадратического отклонения энергии возникновения отказа и

k3

2

k3

 

l d

 

n =

 

 

= const; α =

 

exp

 

 

(7.21)

 

2

 

k4

 

 

Ik4

 

L

 

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

4)n = const, α = const. Этот случай будет иметь место при c = 0

вформулах (7.10), (7.11), (7.14)–(7.19).

Однако одновременное постоянство параметров n и α на практике маловероятно. Это следует хотя бы из того, что при b = c = 0, как видно из формул (7.8), (7.9) и (7.12), (7.13), должны быть постоянными величинами среднее значение и среднее квадратическое отклонение энергии возникновения отказа, но это возможно только в случае полного отсутствия дефектов в материалах элементов, т.е. при m = 0, или в том случае, когда число микроскопических дефектов стремится к бесконечности.

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

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

167

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

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

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

7.4. Методика оценки и прогнозирования надежности элементов

Экспериментальные исследования механического и электрического разрушения различных материалов показывают, что при интенсивностях энергетического воздействия І > 0,7Iр (Iр – интенсивность воздействия, при котором происходит практически мгновенное разрушение элементов) зависимость среднего значенияE и среднего квадратического отклонения σE энергии возникновения отказа от размеров дефектных структур материалов может быть аппроксимирована гиперболой, а при І < 0,7Iр – экспонентой. Следо-

168

вательно, по режиму испытаний элементов можно судить о типе закономерности изменения параметров n и α.

Для простоты дальнейшего изложения допустим, что испытания партии элементов с известными размерами дефектных структур проводятся при І > 0,7Iр, т.е. имеет место гиперболическая зависимость среднего значения и среднего квадратического отклонения энергии возникновения отказа от размеров дефектных структур. Тогда оценка и прогнозирование надежности элементов и критических размеров дефектных структур материалов будет проводиться следующим образом:

1. Из соотношения (7.12)

 

 

 

l b

 

It

 

 

= k3 = const,

(7.22)

 

 

 

 

L

 

 

для заданного I и известного размера l определяется показатель степени

 

 

b =

 

lg[k3 / (It)]

 

 

 

 

 

 

 

.

(7.23)

 

 

 

 

 

 

 

 

 

lg(l / L)

 

 

 

2. Поскольку согласно формулам

 

 

 

 

 

 

n

 

 

n

 

It =

и I 2σ t2=

 

 

 

 

 

α

 

α 2

 

и σt изменяются пропорционально 1/ α

, то для соблюдения этого

t

условия показатель степени с в формуле (7.13) должен быть равен b. Тогда константа k4 в формуле (7.13) может быть определена сле-

дующим образом. Из формулы

n =

t

2

находим выражение для σt и

 

 

 

σ

t2

 

 

 

 

 

 

приравниваем его к формуле (7.13):

 

 

 

 

 

 

t 2

 

 

 

k4

σ t=

 

 

=

 

 

.

 

n

(l / L)c

С учетом равенства (7.22) получаем

169

 

k4

 

 

=

 

 

 

k3

 

 

 

2 1

=

 

 

k3

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

(l / L)

c

 

 

 

 

 

 

 

b

 

 

 

 

b

n

 

 

 

 

 

 

(l / L)

 

n (l / L)

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k4 = k3

 

 

1 l

 

cb

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

(7.24)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n L

 

 

 

 

 

 

 

 

 

 

 

 

 

Но поскольку b = c, то окончательно будем иметь

 

 

 

 

 

 

 

 

 

 

 

k4 = k3

 

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.25)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, выражение для σt будет иметь вид

 

 

 

 

 

 

σ t=

 

 

k3

 

 

 

 

 

 

 

1

.

 

 

 

 

 

 

 

 

 

 

(7.26)

 

 

 

 

 

 

(l / L)b

 

 

n

 

 

 

 

 

 

 

 

 

 

3. Задавшись значениями характерных размеров дефектных

структур материалов элементов l1, l2, …, lk по формуле

 

 

 

 

k3

l

b

 

 

 

 

 

k3

 

 

 

l b

 

 

n

 

 

 

l b

 

 

α =

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

,

(7.27)

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

Ik3

 

 

 

 

Ik4

L

 

 

Ik3 / n

L

 

 

 

 

L

 

определить значение параметра α i для каждого li и построить график зависимости α = f(l) (рис. 7.1).

Рис. 7.1. Зависимость параметра α i от размера дефектных структур

170