книги / Методы обеспечения надежности изделий машиностроения
..pdfДля нахождения параметров распределения Ро, ^тр, Щ используют метод максимального правдоподобия, для чего необ ходимо составить функцию правдоподобия:
(5.16)
где /г, — число испытаний между г-й и (г + 1 )-й доработками; т / — число отказов в я,- испытаниях, после которых проводилась доработка.
Подставляя в формулу (5.16) вместо Р, его значение из выражения (5.15), получим искомую функцию правдоподобия, зависящую от трех параметров модели:
Уравнение (5.17) включает в себя известные опытные дан ные mif я/ и неизвестные параметры Ртр, Ро, а*. Оценками мак
симального правдоподобия Ртр, Ро, а этих параметров будут такие значения Ртр, Ро, а, при которых функция (5.17) обра щается в максимум. Для решения этой задачи используют стан дартные программы минимума нелинейной функции. Для нахожде ния параметров перейдем к логарифмической отрицательной функции правдоподобия, минимум которой совпадает с макси мумом выражения (5.17). Сокращаем эту функций на постоян ный множитель:
— InL =
Поскольку функция (5.18) является нелинейной, то необ ходимо убедиться в том, что она имеет только один минимум.
Функция |
(5.18) |
задана |
на выпуклом множестве, поскольку по |
|||
логике область |
изменения |
параметров |
находится |
в пределах |
||
O ^ a ^ l , |
O ^ P o ^ l, 0 |
^ Р |
тр^ 1 . Так |
как функция |
дифферен |
цируема, то необходимым и достаточным условием выпуклости является положительная полуопределенность матрицы вторых производных, т. е. в нашем случае матрицы
А = V 2 [—1п(Л, Л)1 |
(5.19) |
121
где
а
П = — вектор параметров модели;
Х = Щ — вектор опытных значений mh щ.
Функция (5.18) является дважды дифференцируемой, в чем можно непосредственно убедиться, взяв производные по парамет рам. Для положительной полуопределенности матрицы А следует показать, что ее квадратичная форма неотрицательна. По опре делению квадратичная форма с учетом введенных обозначений может быть представлена зависимостью [5]
В = IV АП, |
(5.20) |
где Пт— транспонированная матрица параметров Я. Развернутый вид выражения (5.20):
|
|
cPlnL |
дЧпL |
<?ЧпL |
a |
|
|
|
|
да2 |
дадРа |
дадРтр |
|
||
|
|
|
|
||||
В = |
[а/50Ртр] |
дР\пL |
(?MnL |
&\nL |
Po |
(5.21) |
|
дР0да |
дГ* |
дР0дРтр |
|||||
|
|
|
|
||||
|
|
d?\nL |
cP\nL |
(T'ln L |
|
|
|
|
|
dP-tfda |
дРтрдР0 |
дР% |
|
|
После преобразований получим
В =
(5.22)
, 2--^— аР ■
дадРтр ТР
Проведенные расчеты показывают, что частные производные, входящие в выражение (5.22), являются отрицательными в обла сти заданных параметров 0 ^ (а , Ро, Ртр)^ 1 , что подтверждает положительность квадратичной формы, а следовательно, и вы пуклость функции (5.18). Таким образом, функция (5.18) имеет единственный минимум в точке
у [ - 1 п ( Я , Х)] = 0, |
(5.23) |
являющейся решением системы уравнений правдоподобия:
dinL |
__ « |
d\nL |
л |
d in t |
__ |
(5.24) |
|
да |
- U ; |
дР0 |
- U ; |
дРтр |
~ |
||
|
122
Для случая, когда доработки проводят только лишь после появления отказов и если причина явно установлена, имеем уравнение правдоподобия вида
- I n L |
= - |
min (1 - |
Р0) - |
mH |
| + !)] ln( |
1 - |
f |
-) - |
|
m |
|
|
|
|
|
|
|
- |
Y, fa |
- 0 4 |
л» - |
(p„ - |
f>.)( i - |
- Я |
] |
<5-25) |
|
i= 1 |
|
|
|
|
|
|
|
Эта функция правдоподобия отражает процесс отработки изделий, при этом в качестве исходных данных могут быть ис пользованы значения параметров, полученные по результатам ис пытаний аналогичных изделий.
5.5. МНОГОПАРАМЕТРИЧЕСКАЯ СТАТИСТИЧЕСКАЯ МОДЕЛЬ РОСТА НАДЕЖНОСТИ
Так как в процессе проведения испытаний причины обнаруженных отказов не всегда достоверны, то проводимые до работки также не одинаково влияют на изменение надежности изделия. Вместе с тем, принимая допущение о том, что в среднем наблюдается определенная тенденция к росту уровня надежности от одной доработки к другой, процесс изменения надежности изделия в результате проведения доработок носит случайный характер. Для восстанавливаемых изделий весь процесс испыта ний можно представить в виде некоторого числа циклов.
Под циклом испытаний понимают функционирование изделия в течение некоторого времени в условиях, соответствующих эксплуатационным при выполнении определенной работы. До работку изделия проводят по всем отказам, причина которых выявлена и по ним требуется изменение конструкции. В этом случае изменение функции надежности происходит в зависимости от номера цикла испытания (/=1, 2 , ..., п). Следовательно, если имеются данные о результатах испытаний, то можно найти ста
тистические оценки Я, хотя ничего неизвестно о параметрах, определяющих поведение этой функции.
Для описания изменения функции надежности Р, может быть использован метод, основанный на представлении Pi в виде функ
ции случайных аргументов, значения которых могут быть опреде лены по результатам испытаний данного изделия или изделия-ана
лога. Такой подход позволяет установить аналитическую форму зависимости, которая является функцией надежности. Таким об разом, статистический метод заключается в том, что по опыт ным данным строят кривую роста надежности произвольной формы, используя экспериментальные точки графика, а затем методами математической статистики находят оценки параметров, входящих в эту функцию. В общем виде функцию можно записать следующим образом:
123
Pi = P (a l,a 2, . . . , a s;i), |
(5.26) |
где as — неизвестные параметры.
Выбрав несколько различных форм функции А, можно оценить, какая из них наилучшим образом описывает экспери ментальную кривую. Так как испытания изделия проводят в циклическом режиме, то можно, как наиболее простую, принять биномиальную модель распределения вероятности доработок.
Для нахождения параметров распределения воспользуемся
методом максимального правдоподобия. С этой целью построим
функцию правдоподобия вида |
|
||
L = п |
|
(1 - а Г ^ - Ч |
(5.27) |
1= |
1 |
(я, —/П,)! v |
|
где mi — число отказов, зафиксированных в щ циклах испытаний; Р, _ вероятность успеха в i-й доработке; k — число проводимых доработок.
Подставляя в (5.27) модель изменения надежности, получим уравнение (5.26), описывающее правдоподобие в зависимости от искомых параметров:
L (OL, nh т.) = |
П:\ |
Pi (al,a 2, . . . , a s-,i)]m‘ X |
||
П —-т------- гг- [1 - |
||||
V *’ |
t=i mf! (я/ - |
т.)! |
|
|
|
X [A (a,, |
oj,.. .,o s; 0 f ~ m' |
(5.28) |
Выражение (5.28) представляет собой функцию правдоподо бия, в которой функция надежности связана с номером доработ ки. Оценками максимального правдоподобия параметров as яв ляются такие значения as, при которых функция правдоподобия принимает максимальное значение:
L (a,, Oj,..., as, nit m) = max. |
(5.29) |
При расчете удобнее использовать не саму функцию, а ее логарифмы, так как максимум и минимум функции L совпадают соответственно с InL и —InL. Для нахождения оценок параметров распределения решают систему уравнений правдоподобия вида
din L |
= 0 . |
(5.30) |
да5 |
Метод максимального правдоподобия позволяет получить несмещенные, состоятельные и эффективные оценки [1]. При нахождении оценок параметров распределения может быть использован также метод наименьших квадратов, который при больших объемах испытаний практически совпадает с методом максимального правдоподобия. В этом случае функция правдо подобия выражается зависимостью вида
124
— inL = |
£ [Д — Pt (fli. 02. • • •. <v. 0]2> |
(5.3i) |
|
|
i=l |
|
|
где Pi — оценка надежности, |
полученная по результатам |
||
испытания в /-й серии. |
квадратов |
хорошо описывается, |
когда |
Метод наименьших |
в /-й серии пС^ 20-=-50 циклов испытаний, т. е. биномиальный закон распределения приближается к нормальному:
Р , = 1 - - ^ - |
(5-32) |
Оценки параметров распределения, обращающие функцию правдоподобия в минимум, находятся из выражения вида
_ |
= о. |
(5.33) |
|
да5 |
|
Таким образом, статистический поход к определению мо дели роста надежности состоит в интуитивном выборе функцио
нальной зависимости Р{а\уаг, |
as\ i) и соответственно нахожде |
|
нии тем или иным способом |
параметров распределения |
а*. |
В качестве функциональной зависимости часто используют |
ал |
гебраические выражения, включающие несколько неизвестных параметров.
5.6. ГИПЕРБОЛИЧЕСКАЯ МОДЕЛЬ РОСТА ФУНКЦИИ НАДЕЖНОСТИ И ДРУГИЕ ВИДЫ ЗАВИСИМОСТЕЙ
Основанием гиперболической модели может служить предположение того, что вся программа испытаний представляет собой k этапов, на каждом из которых производят определен ное число испытаний и в каждом фиксируется отказ или ус пешная работа [10]. По результатам испытаний строят кривую роста надежности — гиперболу, описываемую зависимостью
Pi = P,Р- - р |
(5-34) |
где P i----- истинная надежность изделия на /-м этапе испытаний; Ртр — требуемое значение надежности; а — параметр, характери
зующий скорость роста надежности; / = 1 , 2 , 3, |
k — этапы |
испытаний. |
|
Модель зависит от двух параметров Ятр и а. Для нахождения оценок, параметров распределения воспользуемся методом мак симального правдоподобия. Тогда для /-го этапа отработки бу дем иметь
Ц = |
(1 — P tf‘, |
(5.35) |
где rii — число испытаний |
на /-м этапе; т , — число |
отказов на |
/-м этапе. |
|
|
125
Полагая, что результаты всех этапов статистически неза висимы, получим
к |
к |
т‘(1 - |
т |
‘ |
(5.36) |
L = П L; = |
П Р-‘ |
|
|||
i=i |
i=i |
|
|
|
|
или после логарифмирования |
|
|
|
|
|
к |
к |
|
|
|
|
InL = £ (nt - m ) In (PTp - f ) + X |
m‘ln 0 |
- |
+ f ) |
(537) |
|
/=1 |
i = 1 |
|
|
|
|
Уравнения правдоподобия можно записать в виде |
|
|
|
|
i - i |
|
яТрТР - т i- |
+ l |
|
|
+ |
=0; |
<538> |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
п{- |
т. |
|
|
|
|
|
|
d\nL |
_ _ V* |
/___ |_V* |
|
0. |
(5.39) |
||||||
|
da |
|
|
|
|
п |
а |
Z-i |
|
|
||
|
|
|
|
|
|
|
|
|
||||
|
|
|
/ = 1 |
^ т р |
I |
i = 1 1 |
^ т р + |
^ |
|
|||
В результате решения уравнений (5.38) и (5.39) получаем |
||||||||||||
оценки параметров: |
|
|
|
|
|
|
|
|
|
|||
|
|
k_ |
|
|
|
|
|
|
|
|
|
|
|
|
х п'Тх *^~т») ~ * х1х |
~т)] |
|
||||||||
|
А |
| = 1 |
L t = j |
|
|
|
/ = 1 |
J |
(5.40) |
|||
|
ос ------------------------------------------------------------------------ |
|||||||||||
|
|
|
|
|
|
|
A + L c _ * |
|
|
|
||
|
|
|
k_ |
|
Г |
k |
|
|
к |
|
|
|
|
|
|
к |
|
|
|
"1 |
|
||||
|
|
X |
п' |
L |
т Х |
‘(я- - т- ) - Х |
|
-J |
|
|||
|
тр |
<•=1 |
|
1 = 1 |
л. + 1 |
|
/•= 1 |
|
(5.41) |
|||
|
|
|
|
|
|
с — k |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|||
где г = |
£ - L « ln |
( * + |
4 - ) |
+Е, Е = 0,57721 |
— iпостоянная |
|||||||
Эйлера. 1 - 1 1 |
|
|
|
|
|
|
|
|
|
|
|
|
В качестве аппроксимирующих зависимостей можно исполь |
||||||||||||
зовать также модели вида |
|
|
|
|
|
|
||||||
или Р = |
1; |
|
Pi = |
Р~ ~ ( Р „ |
— |
Р0) е |
а' |
|
(5.42) |
|||
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
Р, = 1 - |
(1 - Р |
0)е - а', |
|
(5.43) |
126
где Pi — надежность на /-м этапе испытаний; Р — предельно достижимое значение надежности; Р0 — начальное значение на дежности; е — основание натурального логарифма; щ—параметр, характеризующий изменение скорости роста надежности.
Применяют также модели роста надежности вида [10]
_ |
Я0 Р„ |
|
(5.44) |
|
Ро + |
(/>» - |
Ро)е~°‘ ’ |
||
|
||||
Р = Р |
6 |
° , |
(5.45) |
|
|
+ С’ |
|
где а, Ьу с — неизвестные параметры.
Наиболее распространенной является модель, описываемая
показательным законом [1 0 ]: |
|
р. = 1 — |
(5.46) |
где А = рР; с = In — |
|
Функция (5.46) зависит от трех параметров р, Р и а, кото рые могут быть найдены методом максимального правдоподо бия или методом наименьших квадратов.
Более общая модель роста надежности [10] имеет вид
Pt = P „ - a f ( i ) , |
(5.47) |
|||
где f(i) — положительная, |
монотонно убывающая |
функция /. |
||
Для нахождения оценок Р^ |
и а |
методом наименьших квад |
||
ратов вычисляют суммы: |
°° |
|
|
|
k |
|
|
k |
|
с, = £ |
ДО; |
<* = |
£/*(*■)• |
|
/= 1 |
|
|
/= 1 |
|
Для удобства вычислений С\ и С2 [10] можно использовать приближенные формулы для нескольких видов f(i\ приведенные в табл. 5.2.
Оценки параметров распределения находят, решая уравнения
|
|
|
1 |
|
к |
« |
- I |
», |
8 |
+ |
|
• |
|
|
|||
|
1 = 1 |
|
|
|
|
dQ |
о у |
Г Щ~ Щ |
|
+ «Д0] = 0 ; |
|
dP„ |
~ L |
L ", |
- - P - |
/ = 1
(5.48)
(5.49)
dQ _ |
о V Г n< ~ m‘ |
Л» + «ДО] ДО = о. |
(5.50) |
да |
L, L Щ |
/=1
127
h
/-1/2
г'
г2
г3
5.2. Приближенные формулы для вычислений с\ и с2
|
С1 |
|
С2 |
|
/ |
1\ |
1/2 |
|
|п(й+ т ) + 0,577 |
2(к+т) |
- |
1-46 |
||
|п(* + у ) |
+ |
0.577 |
т -(‘+ т)" |
|
|
т-(*н-4)"
' ■ - 1 ( * + ! ) "
_ 2 |
4 |
6 |
П р и м е ч а н и е . — = |
1,645; -г— = |
1,082; -гг— = 1,017. |
о |
90 |
945 |
В результате |
решения уравнений (5.49) и (5.50), проведен |
||
ных в работе |
[1 0 ], получено |
|
|
р |
= |
. != 1________________ ini_________________• |
(5.51) |
00 |
|
kc2 - c \ |
|
|
|
k |
|
|
|
f(0 |
|
* |
_ |
i = 1 |
(5.52) |
ОС = |
------------- |
||
|
|
kc2 — c\ |
|
где /лi — число |
отказов на i-м этапе; nt — число испытаний на |
||
/-м этапе. |
|
|
|
Изменение надежности изделия при его поэтапном испыта нии можно также оценить статистически [1 0 ].
Оценку надежности после /-й серии испытаний определяют
зависимостью вида |
|
/>.= аЯ/ + (1 -а)Я ,_ ,, |
(5.53) |
где Р! — оценка надежности по результатам только i-й серии испытаний; а — коэффициент, характеризующий вес каждой оценки; Я, - 1 — оценка надежности по результатам всех испыта ний, проведенных до /-й серии.
128
Примем следующее допущение — надежность испытываемых образцов в каждой серии одинакова. Особенности этой модели таковы, что при а = 1 , 0 не будут учитываться все испытания до /-й седин, а при а = 0 — результаты последней серии. Величину оценки Р\ определяют по результатам испытаний по формуле
3 = |
1 “ |
. |
(5.54) |
Значением параметра а |
можно |
лишь задаваться, |
исходя |
из практических соображений, что снижает эффективность использования этой модели.
5.7. ВЕРОЯТНОСТНЫЕ МОДЕЛИ, ОСНОВАННЫЕ НА ЛОГИЧЕСКИХ ПРЕДПОСЫЛКАХ
Рассмотрим построение логико-вероятностной модели, основанное на предположении того, что изделие состоит из опреде
ленного числа k блоков, при |
этом |
вероятность |
отказа каждо |
го блока qi= 1 —Pi\ отказы |
блоков |
являются |
независимыми,и |
отказ любого из них приводит к отказу изделия |
в целом; отказ |
каждого блока можно устранить с некоторой вероятностью г*. С учетом принятых допущений вероятность отказа 1-то блока после п испытаний определяется соотношением
|
|
Q, = (1 - |
Я,г,Т> |
|
|
(5-55) |
||
а вероятность безотказной работы изделия |
после п |
испыта |
||||||
ний — соотношением |
|
|
|
|
|
|
||
|
К = |
П (1 - Q,<7,) = |
П |
[1 - q, (1 |
- |
qtr$]. |
(5.56) |
|
|
|
i'= 1 |
|
i= 1 |
|
|
|
|
|
Принимая допущение о постоянстве величин qi = q, rt = r, имеем |
|||||||
|
|
Pn = |
[ l - q ( l - q r ) n]k |
|
|
(5.57) |
||
Так как произведение ^г<С 1, то можно записать: (1—qr)n~ |
||||||||
« е |
чгп Подставим полученное |
выражение |
в |
формулу |
(5.57): |
|||
|
|
/>„ |
« (1 - |
|
|
|
|
(5.58) |
|
Значение qe |
^ " < 1 , поэтому |
|
|
|
|
|
|
|
|
Рп « |
- ( 1 - kq)e~qrn, |
|
|
(5.59) |
||
где |
1 —kq — вероятность |
успеха |
до первой |
доработки |
(началь |
|||
ная надежность). |
qr = R, |
1 — kq = Po |
и подставим их |
|||||
|
Введем обозначения: |
в формулу (5.59). В результате получим двухпараметрическую экспоненциальную модель роста надежности:
Ря = 1 - Р 0 е -Лп, |
(5.60) |
129
Полученная модель, как указывалось выше, основана на предположениях и должна быть подтверждена статистическими
данными.
При построении моделей роста надежности сложных техни ческих систем была использована теория стохастических моде лей обучаемости. Изменение надежности изделия при отработке можно представить в виде случайного процесса. Рассмотрим частные свойства случайного процесса, которые могут быть использованы при построении кривой роста надежности. Так, на пример, может быть принято допущение, что: 1 ) выявление при чины отказа и устранение его доработкой проводят только после появления отказа; 2 ) надежность изделия может изменяться только после доработок, причем каждая доработка либо не из меняет надежность, либо увеличивает ее.
Практика отработки сложных технических систем показыва ет, что доработку изделий многоразового использования прово дят только после появления отказа, когда причина явно уста новлена. В этом случае процесс отработки представляет собой следующую логическую основу. Пусть Я, — вероятность успеха в /-м испытании; Л1 — вероятность проведения доработки после успеха; Л2 — вероятность проведения доработки после отказа. После каждого испытания могут быть два несовместных события: успех или отказ. Так как принято допущение, что доработку про водят только в случае появления отказа, то вероятность отсут ствия доработки после успеха равна ( 1 — jii)= l, а после отказа соответственно ( 1 — лг). Отсюда приращение надежности в резуль тате доработок после успеха равно нулю, т. е. ДРн = 0, а после отказа равно ДЯг/. Схема рассматриваемого процесса приведена в табл. 5.3.
Соотношение для определения приращения надежности за пишется в виде
5.3. Схема процесса изменения надежности
Исходы |
Вероят |
|
|
/-го |
|
||
ность |
Последствие |
||
испы |
|||
исхода |
|
||
тания |
|
||
|
|
||
Успех |
Pi |
Доработка |
|
А и |
проводится |
||
|
|
Доработка |
|
|
|
не проводится |
|
Отказ |
i - Pi |
Доработка |
|
А21 |
проводится |
||
|
|
Доработка |
|
|
|
не проводится |
Вероятность
последствия
—
Pi
л2(1 - P i )
0-**)(!-Р<)
Вероятность
успеха в (/+ 1)-м
испытании
—
Pi
Р<+ДР2,-
Pi
130