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

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

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

 

q(ε )= Ci qi (ε );

C1= 1; Ci0.

(5.17)

i=1

i =1

 

Определим функцию распределения следующим образом:

 

)=

l

1

l

q

 

 

q

C q)+

C

).

(5.18)

(l )

 

i i

 

i

i

 

 

 

 

i =1

 

i=1

 

 

 

 

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

 

 

l

 

 

 

 

l

 

 

 

 

S(l ) Ci

Si

+ 1

Ci Si .

(5.19)

 

 

i=1

 

 

 

 

i =1

 

 

Заметим, что

 

 

 

 

 

 

 

 

 

 

l

 

l

 

 

 

 

 

lim Ci Si

+ 1

Ci

Si

= Ci Si + 0.

(5.20)

l →∞

i=1

 

 

i =1

 

 

 

i =1

 

Будем обозначать впредь через m(K) и mi(K) K-й момент распределений q(ε) и qi(ε) соответственно. Из условия (5.17) следует, что

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

ε

K dq)=

 

 

i

 

 

 

 

 

 

 

 

 

 

C

 

ε K dq

 

),

(5.21)

 

 

 

 

0

 

 

 

i=1

 

 

0

 

 

 

 

следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m( K )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= Ci mi( K ) .

 

 

(5.22)

 

 

 

 

 

 

 

 

i =1

 

 

 

 

 

 

Если m(K) есть K-й момент распределения q(l)(ε), то согласно

формуле (5.18)

 

 

 

 

 

 

 

 

 

 

 

 

 

m( K ) =

ε K dq (ε) =

l

C

ε

K dq(ε) +

1

 

l

C ε

K d(ε)q ,

(5.23)

(l )

(l )

i

i

 

 

 

 

i

i

 

 

0

 

i =1

 

0

 

 

 

 

 

i =1

 

0

 

 

следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

limm((lK) ) = Ci mi(K ) + 0.

 

 

(5.24)

 

 

 

 

l →∞

 

i =1

 

 

 

 

 

 

 

Сопоставим между собой выражения (5.22) и (5.24) и получим

91

m(K ) = m(K ) .

(5.25)

lim

(l )

 

l →∞

 

 

Аналогично

 

 

 

limS (l) = S.

(5.26)

 

l →∞

 

Из соотношения (5.19) с учетом выражений (5.20) и (5.26) сле-

дует, что

 

 

 

 

S Ci Si .

(5.27)

i =1

Таким образом, утверждение 1 доказано в общей форме для функции q(ε), определяемой выражением (5.17).

Справедливость доказанного утверждения подтверждается и результатами испытаний на разрушение образцов из алюминия, латуни и алюминиевой фольги, представленными в табл. 5.1–5.3.

Таблица 5. 1

Результаты испытаний на механическое разрушение образцов из дюралюминия марки Д16АТ при скорости нагружения v = 20 мм/мин

 

 

 

 

 

σE,

 

 

 

 

2

 

K

 

 

 

,

 

 

 

 

E

 

 

E

 

 

 

Состав партии

 

 

 

2

Дж/мм

2

S =

 

 

 

 

Ci Si

Дж/мм

σ

2

 

 

 

 

E

 

i =1

N1 = 106 образцов с отв. dср = (8 ± 1) мм

1,38

 

0,130

 

112,67

 

112,67

N2 = 105 образцов с отв. dср = (10 ± 2) мм

1,20

 

0,136

 

77,85

 

77,85

N3 = 99 образцов с отв. dср = (12,5 ± 3) мм

0,78

 

0,116

 

45,21

 

45,21

N4 = 104 образцов с отв. dср = (14 ± 4) мм

0,50

 

0,101

 

24,51

 

24,51

N1 + N2

1,29

 

0,164

 

61,87

 

95,34

N1 + N3

1,09

 

0,325

 

11,25

 

80,09

N1 + N4

0,94

 

0,464

 

4,11

 

 

69,01

N2 + N3

0,99

 

0,269

 

13,54

 

62,01

N1 + N2 + N3

1,13

 

0,269

 

17,65

 

79,33

N1 + N2 + N4

1,03

 

0,397

 

6,73

 

 

71,96

N2 + N3 + N4

0,83

 

0,310

 

7,17

 

 

49,35

N1 + N2 + N3 + N4

0,97

 

0,367

 

6,98

 

 

65,56

92

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 5. 2

Результаты испытаний на механическое разрушение образцов из латуни марки Л63М при скорости нагружения v = 20 мм/мин

 

 

 

 

 

 

σE,

 

 

 

 

2

 

K

 

 

 

 

E

,

 

 

E

 

 

 

Состав партии

 

 

 

2

S =

 

 

 

 

Ci Si

 

 

Дж/мм2

σ

2

 

 

 

Дж/мм

 

E

 

i =1

N1 = 98 образцов с отв. d = 2 мм

2,53

0,402

39,61

 

39,61

N2 = 162 образца с отв. d = 3 мм

2,04

0,317

41,41

 

41,41

N3 = 98 образцов с отв. d = 4 мм

1,76

0,285

38,14

 

38,14

N4 = 112 образцов сотв. d = 5 мм

1,60

0,257

38,76

 

38,76

N1 + N2

 

2,22

0,429

26,78

 

40,73

N1 + N3

 

2,13

0,505

17,79

 

38,87

N1

+ N4

 

2,03

0,557

13,28

 

39,16

N1

+ N2

+ N3

2,1

0,443

22,47

 

40,02

N2

+ N3

+ N4

1,83

0,348

27,65

 

39,75

N1

+ N2

+ N3 + N4

1,98

0,458

18,69

 

39,72

Таблица 5. 3

Результаты испытаний на электрическое разрушение образцов из алюминиевой фольги при W = 13,2 Вт

 

 

 

 

 

 

 

σE,

 

 

 

2

 

K

 

 

 

 

 

E

,

 

 

E

 

 

 

 

Состав партии

 

 

 

2

S =

 

 

 

 

Ci Si

 

 

 

 

 

 

σ

2

 

 

 

 

2

Дж/мм

 

E

 

i 1

 

 

 

 

Дж/мм

 

 

 

 

 

=

N1 = 100

образцов с отв. d = 1 мм

884,86

377,43

5,5

 

 

5,5

N2 = 100

образцов с отв. d = 1,5 мм

536,34

238,96

5,04

 

 

5,04

N3 = 100

образцов с отв. d = 2 мм

299,26

134,64

4,94

 

 

4,94

N4 = 100

образцов с отв. d = 3 мм

97,40

42,9

5,15

 

 

5,15

N1 + N2

 

 

710,6

379,15

3,51

 

 

5,27

N1 + N3

 

 

592,06

426,05

1,93

 

 

5,22

N1 + N4

 

 

491,13

479,84

1,05

 

 

5,32

N1

+ N2

+ N3

572,70

382,95

2,24

 

 

5,16

N1

+ N2

+ N4

506,20

413,44

1,5

 

 

5,23

N2

+ N3

+ N4

311

240,56

1,67

 

 

5,04

N1

+ N2

+ N3 + N4

453,68

389,91

1,35

 

 

5,15

93

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

5.4.Минимальное число квазичастиц, вызывающих отказ элементов

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

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

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

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

94

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

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

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

Обозначим через ri число событий поглощения квазичастиц класса i (i = 1, 2, …, K), происходящих в материале элемента испытуемой совокупности. Пусть вектор Rε есть вектор с K компонентами ri. Индекс ε у вектора R указывает на то, что он является случайной величиной, зависящей от величины энергетического воздействия ε.

Пусть B – случайная переменная, которая принимает значение 0 при отсутствии отказа и значение 1 при возникновении отказа. Вероятность отказа отдельного элемента в результате энергетического воздействия является функцией числа происходящих событий поглощения квазичастиц различных видов, и поэтому обозначим ее как P(B = 1/Rε = R).

95

Если через P(Rε = R) обозначим вероятность того, что при энергетическом воздействии ε в определенном элементе произошло точно ri событий поглощения квазичастиц класса i (i = 1, 2, …, K), то вероятность отказа согласно формуле полной вероятности будет определяться выражением

q(ε) = P(Rε = R) P(B= 1|) / Rε = R.

(5.28)

R A

 

При этом суммирование происходит по множеству A всех векторов R с целочисленными неотрицательными компонентами.

Будем считать, что события поглощения квазичастиц являются статистическими независимыми. Пусть α i ( α i ≥ 0) есть среднее значение числа событий поглощения квазичастиц i-го вида на единицу объема и на единицу энергетического воздействия. Тогда α i ε есть

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

Тогда вероятность осуществления точно ri событий поглощения квазичастиц класса i при энергетическом воздействии ε есть

P(r

= r ) = exp(−α ε) (α εi

)ri .

(5.29)

i,ε

i

i

ri

!

 

 

 

 

 

Отсюда следует

P(Rε = R) = K

P(ri,ε = ri ) =

i=1

 

= exp(−α εi ε) n K

i=1

где

K

α = α i и

i=1

K

exp(−α εi

) (−α εi

) i=

 

 

 

 

r

i=1

 

 

ri !

 

 

α iri ,

ri !

K

n = ri = R .

i =1

Очевидно, что

96

exp(−αε ε)

n

 

α i=

exp(−αε

)

(αε )

 

nn!

 

α

i .

 

 

K

ri

 

 

 

 

n

 

K

 

ri

 

 

i=1 ri !

 

 

n! α

i =1 ri !

На основании этой зависимости формально можно говорить о наложении K простых процессов Пуассона. Разумеется, что в этом случае изменяется не время, а величина энергетического воздействия ε. Если ввести обозначение

 

f (R) =

n!

K

α iri

,

 

 

 

n

 

 

 

 

 

α

i =1 ri

 

 

 

то

 

 

 

 

 

 

 

 

P(R

= R) = exp(−αε

) (αε

)n

f (R) .

(5.30)

ε

 

 

 

n!

 

 

 

 

 

 

 

 

 

Теперь вычислим условную вероятность наступления отказа в том случае, когда произошло точно n событий поглощения квазичастиц. Используя уравнение (5.30), можно показать, что эта вероятность не зависит от ε:

 

 

 

 

 

P(B =1, nε = n)

 

P(B =1, Rε = R)

P(B =1 nε = n) =

=

R

=n

 

=

 

 

 

 

 

 

 

 

 

P(nε = n)

 

 

 

P(Rε

 

 

 

 

 

 

 

 

 

 

= R)

 

 

 

 

 

 

 

 

 

 

R

 

=n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P(Rε = R) P(B =1 Rε = R)

 

 

=

 

R

 

=n

 

 

 

 

.

(5.31)

 

 

 

 

 

 

 

 

 

P(Rε

= R)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R =n

Вероятность наступления отказа зависит только от поглощенной энергии и ее пространственно-временного распределения. Следовательно, P(B = 1/Rε = R) не зависит от величины ε. Если для краткости ввести обозначение

h(R) = P(B = 1/Rε = R),

(5.32)

то уравнение (5.31) с учетом соотношения (5.30) можно записать в следующем виде:

97

 

 

exp(−αε

) (αε )n f (R)h(R)

 

 

 

f (R) h(R)

 

P(B =1 n

= n) =

 

 

n!

 

R

 

=n

=

 

 

R

 

=n

. (5.33)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(αε

 

)n

 

 

 

 

 

 

 

ε

 

 

f (R)

 

 

 

 

 

 

f (R)

 

 

 

exp(−αε

)

 

 

 

 

 

 

 

 

 

 

R

 

=n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n!

R

=n

 

 

 

 

 

 

 

 

 

 

Тем самым мы показали, что вероятность отказа элемента при условии, что произошло точно n событий поглощения квазичастиц, не зависит от ε. Обозначим эту вероятность для краткости qn:

qn = P(B = 1 nε = n).

(5.34)

Вероятность qn изменяется дискретно и с ростом n (n = 0, 1, 2…) монотонно возрастает от 0 до 1. Это очевидное утверждение вытекает из предположения о том, что дополнительное поглощение энергии материалом не может понизить вероятность отказа элемен-

та, т.е.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qn+1 qn, n = 0, 1, 2… .

 

 

 

(5.35)

Таким образом, с учетом выражения (5.34) можно записать

следующее уравнение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q(ε )= P(B= 1 nε = n) P(nε = n)=

 

 

 

 

 

 

n=0

 

 

 

 

 

 

 

 

 

 

(5.36)

 

 

 

= qn exp(−αε

)

(αε )

.

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

n=0

 

 

 

 

n!

 

 

 

 

 

 

Производная этого выражения есть плотность вероятности

энергии разрушения:

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ (ε) = δθ(ε) =

 

ξ π(−αε) ε)

 

ν

 

ξ π(−αε)(αε)

 

ν−1

 

 

 

θ

−αε

 

+ αε

 

 

 

=

 

 

 

 

 

 

δε

ν

 

 

ν!

 

 

 

(ν −1)!

 

 

 

ν=1

 

 

 

 

 

 

 

(5.37)

 

 

 

(αε)

ν −1

 

 

 

ε)

 

ν−1

 

 

 

 

 

 

 

 

= ν − θ ν−1 αε)

ξ π(−αε)

 

=

ϕν αε ξ π(−αε)

,

 

ν=1

 

 

 

(ν −1)!

ν=1

 

 

 

(ν −1)!

 

где fn = qn – (qn

– 1) есть разность между вероятностями отказа по-

сле n событий поглощения квазичастиц и после n – 1 событий поглощения квазичастиц; ξ – постоянный коэффициент, ξ = α n.

98

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

 

 

 

 

 

= ndq =nfn .

(5.38)

n

0

n=1

 

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

Плотность вероятности энергии разрушения, согласно уравнению (5.37), является выпуклым наложением плотностей гаммараспределений, для которых, как известно, среднее значение ε i и

дисперсия σ i2 равны

ε i =

i

; σ i2=

i

,

α

α 2

следовательно, относительная крутизна

S

 

=

ε i2

= i.

i

σ i2

 

 

 

Согласно утверждению 1 в его общей форме из уравнения (5.37) следует, что

 

S fi Si =

fii,

(5.39)

i =0

i=0

 

причем S есть относительная крутизна кривой q = f(ε). А на основании уравнения (5.38) получаем

S

 

,

(5.40)

n

что и требовалось доказать.

Если используемая в опытах совокупность элементов неоднородна, то ее можно разделить на L классов однородных элементов. Разделение на конечное число классов является вполне обоснованным, так как L может быть выбрано сколь угодно большим, кроме того, мы всегда имеем дело с конечным числом испытуемых элементов. Зависимость вероятности отказа от величины энергетического воздействия ε является, таким образом, выпуклым наложени-

99

ем распределений qi(ε) (i = 1, 2,…, L), соответствующих отдельным классам:

L

 

q(ε) = Ci qi(ε) ,

(5.41)

i =1

где

L

Ci =1 и Ci ≥ 0,

i=1

Ci – относительная доля элементов данной структуры материалов i (i = 1, 2, …L) в совокупности.

Выше было показано, что соотношение

S

 

=

ε i2

 

 

 

n

 

σ i2

 

i

 

 

i

справедливо для каждой кривой qi(ε) идентичных элементов. Но тогда, согласно утверждению 1 для распределения q(ε), являющегося выпуклым наложением распределений qi(ε), справедливо соотношение

S =

ε 2

 

.

n

σ 2

 

 

 

 

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

Материал элемента, поведение которого зависит от величины энергетического воздействия, можно рассматривать как некий измерительный прибор, который регистрирует превышение величины энергетического воздействия над некоторым значением ε. Таким измерительным прибором, очевидно, можно измерить энергию, поглощенную в определенном объеме материала, иначе говоря, определить вид и число событий поглощения квазичастиц. Здесь мы должны оставить открытым вопрос, каким образом происходит измерение. В частности, измерительный прибор может «ошибаться», т.е. содержать элемент случайности; на точность измерений могут влиять и восстановительные процессы.

100