
книги из ГПНТБ / Паньков, Н. П. Ремонтопригодность автомобильной техники учебное пособие
.pdfвы'м начальным моментам, определяет собой среднее время без отказной работы. Ето можно найти но зависимости
|
00 |
|
|
|
|
|
|
(1.29) |
||
|
М{%}= J t f ( t ) d t . |
|
|
|
|
|||||
|
О |
|
|
|
|
|
|
|
|
|
Подставив в зависимость (1.29) |
выражение |
(1.28), |
получим |
|||||||
|
оо |
|
|
|
-I |
|
|
|
|
|
М {*} = J t \ e ~ Xt dt = l f |
te~lt d t = - f |
|
|
(1.30) |
||||||
О |
0 |
|
|
|
|
|
|
|
|
|
Плотность экспоненциального распределения для различных |
||||||||||
значений .параметра |
X. приведена на |
рис. 1.7. |
|
|
|
|
|
|||
|
|
|
Рассмотренная схема |
|||||||
|
|
|
минованных |
отказов |
по |
|||||
|
|
|
зволяет обратить .внима |
|||||||
|
|
|
ние |
на |
следующие |
два |
||||
|
|
|
момента: |
|
|
|
|
|||
|
|
|
1) отказ детали на |
|||||||
|
|
|
ступает |
не в |
результате |
|||||
|
|
|
постепенного |
|
изменения |
|||||
|
|
|
ее |
свойств |
|
|
(износа, |
|||
|
|
|
структуры и т. д.), а |
|||||||
|
|
|
лишь как следствие внеш |
|||||||
|
|
|
него |
случайного |
воздей |
|||||
|
|
|
ствия, имеющего величи |
|||||||
|
|
|
ну |
больше |
допустимой;' |
|||||
|
|
|
2) уровень |
предельно |
||||||
|
|
|
допустимой |
|
|
нагрузки |
||||
|
|
|
Qn |
остается |
постоянным |
|||||
|
|
|
в течение всего |
времени |
||||||
|
|
|
эксплуатации |
|
машины. |
|||||
|
|
|
Отсюда 'следует важ- |
|||||||
|
|
|
т яый для практики вывод. |
|||||||
|
|
|
При |
|
экспоненциальном |
|||||
Рис. 1.7. Плотность экспоненциального рас |
распределении |
|
времени |
|||||||
безотказной |
работы |
нет |
||||||||
пределения для различных значений пара |
||||||||||
метра X. |
|
смысла |
прибегать |
к |
||||||
|
|
|
предварительной |
заме- |
не деталей с целью повышения технической готовности машины, поскольку отказ наступает в результате внешнего воздействия. Путями повышения безотказности в данном случае будут либо конструктивное улучшение детали, либо улучшение условий ее работы за счет снижения действующих нагрузок.
Рассмотрим теперь случай появления мгновенных отказов вследствие технологических дефектов производства. В этом слу чае влияние пиковых нагрузок будет проявляться, прежде всего, на деталях, имеющих технологические дефекты.
20
Обозначим |
через д долю |
дефектных деталей, а через |
1 |
— д |
долю высококачественных. |
|
|
|
|
Вели из общей совокупности изготовленных деталей выбрать |
||||
одну, то она |
с вероятностью |
1—д будет .принадлежать ко |
II |
и с |
вероятностью д к I группе деталей.
.Вероятность же безотказной работы детали в течение времени
Т у |
наугад выбранного |
экземпляра |
машины |
определится |
по |
||
формуле |
Р(Т) = ( 1 - д ) Р 1(Т) + дР,(Т) , |
|
(1.31) |
||||
|
|
|
|||||
где |
Рг (Г) — вероятность безотказной работы |
детали, |
не |
име |
|||
|
|
ющей технологических дефектов; |
|
с техно |
|||
|
Р3 (Т) — вероятность безотказной работы деталей |
||||||
|
|
логическим дефектом. |
|
|
|
|
|
Плотность распределения времени безотказной работы в этом |
|||||||
случае определяется по |
формулам: |
|
|
|
|
||
|
при |
Т < О |
f ( T ) = 0 |
; |
|
|
|
|
при |
Г > О |
|
|
|
||
|
|
|
|
(1.32) |
|||
|
|
f(T) = (\ — д)11е - 1'т+ дХ2е-'^т , |
где — параметр деталей, не имеющих дефектов; 12— параметр деталей с дефектами.
Ощенить значения указанных величин можно следующим обра зом. Обозначим:
N — общее количество отказавших деталей;
—количество деталей из общего числа, которые отказали по технологическим дефектам;
Ti. т2... — время безотказной работы деталей с техноло
|
|
|
гическими дефектами; |
|
|
||||
Т2 • **т N-Ni — время |
безотказной |
работы деталей |
|
||||||
Тогда |
|
|
тов. |
|
|
|
|
|
|
|
|
|
N - N , |
|
|
|
|||
д |
|
• |
ч ... |
Л2 |
|
(1.33) |
|||
д/ |
jv- ^ |
5 |
JV, |
||||||
|
> |
||||||||
|
|
|
|
2 |
^ |
|
7=1 |
|
|
|
|
|
|
(=1 |
|
|
|
На рис. 1.8 приведены кри вые плотности распределения, построенные на основании фор мул 1.28 и 1.32. Анализ их по казывает, что кривая плотнос ти распределения, построен ная по формуле 1.32, имеет белее крутой спуск.
Рис. 1.8. Плотность распределения времени.безотказной работы деталей, имеющих дефекты.
21
|
|
СХЕМА ВОЗНИКНОВЕНИЯ |
ИЗНОСНЫХ |
ОТКАЗОВ |
|
|
||||||||||
И ЗАКОНЫ РАСПРЕДЕЛЕНИЯ ВРЕМЕНИ БЕЗОТКАЗНОЙ РАБОТЫ |
||||||||||||||||
Износ удовлетворительно 'сконструированного сопряжения |
||||||||||||||||
характеризуется, как правило, тремя периодами: |
|
I —периодом |
||||||||||||||
приработки, |
II —периодом |
установившегося или |
нормального из |
|||||||||||||
|
|
|
|
|
|
|
|
носа |
и |
III—периодом |
||||||
|
|
|
|
|
|
|
|
ф орсир01Ванно/го |
износа. |
|||||||
|
|
|
|
|
|
|
|
На рис. 1.9 приведе |
||||||||
|
|
|
|
|
|
|
|
на типичная кривая из |
||||||||
|
|
|
|
|
|
|
|
носа. |
|
|
приработки I |
|||||
|
|
|
|
|
|
|
|
Период |
||||||||
|
|
|
|
|
|
|
|
характеризуется |
тем, |
что |
||||||
|
|
|
|
|
|
|
|
трущиеся |
|
|
поверхности |
|||||
|
|
|
|
|
|
|
|
приспосабливаются |
друг |
|||||||
|
|
|
|
|
|
|
|
к другу и к условиям на |
||||||||
|
|
|
|
|
|
|
|
гружения. В них проис |
||||||||
|
|
|
|
|
|
|
|
ходят изменения в на |
||||||||
|
|
|
|
|
|
|
|
правлении |
создания |
оп |
||||||
|
|
|
|
|
|
|
|
тимальной |
для |
данных |
||||||
|
|
|
|
|
|
|
|
условий |
работы шерохо |
|||||||
|
|
|
|
|
|
|
|
ватости |
и |
микротвердо |
||||||
|
|
|
|
|
|
|
|
сти. Скорость износа в |
||||||||
|
|
|
|
|
|
|
|
этот |
период |
постепенно |
||||||
|
|
|
|
|
|
|
|
уменьшается |
и |
отличает |
||||||
Рис. 1.9. |
Типичная кривая износа и соответ |
ся |
взаимной |
зависимо |
||||||||||||
стью |
приращений |
из |
||||||||||||||
|
ствующая ей кривая отказов. |
|
||||||||||||||
|
|
|
|
|
|
|
|
носа. |
|
|
|
|
|
|
||
На рис. 1.10 показано изменение интенсивности износа в ходе |
||||||||||||||||
приработки. |
В начале |
приработки |
(на участке |
Tv Тх-f- Д Т) |
при |
|||||||||||
ращение |
износа |
состави |
|
|
|
|
|
|
|
|
|
|
||||
ло |
A v 1(t), |
а |
в конце |
|
|
|
|
|
|
|
|
|
|
|||
приработки |
(на |
участке |
|
|
|
|
|
|
|
|
|
|
||||
Т2, r 2 + |
A T ) - A v 2(t). |
|
|
|
|
|
|
|
|
|
|
|||||
При этом приращение из |
|
|
|
|
|
|
|
|
|
|
||||||
носа |
Д v2{t) |
зависит |
от |
|
|
|
|
|
|
|
|
|
|
|||
того, какова была его ве |
|
|
|
|
|
|
|
|
|
|
||||||
личина |
в предшествую |
|
|
|
|
|
|
|
|
|
|
|||||
щий момент времени. |
|
|
|
|
|
|
|
|
|
|
|
|||||
Период установивше |
|
|
|
|
|
|
|
|
|
|
||||||
гося или нормального из |
|
|
|
|
|
|
|
|
|
|
||||||
носа |
II |
характеризуется |
|
|
|
|
|
|
|
|
|
|
||||
постепенным |
ростом |
за |
|
|
|
|
|
|
|
|
|
|
||||
зора |
в |
сопряжении |
без |
|
|
|
|
|
|
|
|
|
|
|||
изменения |
физической |
|
|
|
|
|
|
|
|
|
|
|||||
картины |
взаимодействия |
|
|
|
|
|
|
|
|
|
|
|||||
трущихся деталей. Вэтой |
Рис. 1.10. |
Изменение! |
интенсивности износа |
|||||||||||||
зоне |
величина |
прираще- |
|
|
в ходе |
приработки. |
|
|
22
ния износа Д©!^) |
не влияет сколько-нибудь ощутимо |
на |
вели |
чину приращения |
Изменения, которые протекают |
в соп |
|
ряжении, носят в основном количественный характер. |
|
|
Период форсированного износа III характеризуется тем, что медленные количественные накопления -износа «а II участке при водят ж скачкообразным качественным изменениям .в состоянии трущихся пар. На интенсивность износа начинают оказывать .воз действие -новые факторы, которые ранее не проявлялись.
. Изменения в состоянии трущихся деталей носят -направленный характер, приращения изн-о-са взаимно зависимы.
На рис. 1.11 приведена реализация процесса износа в период
разрушения |
сопряжения. |
|
Из |
-рис. |
1Л1 следует, чт-о .приращение износа на участ |
ке 7\, |
Тх А- Д Т вызывает еще .большее увеличение интенсивности |
|
износа на участке Т2, Г2 -f- АТ. |
.-■Возьмем -сопряжение «вал—подшипник» двигателя. Форсиро ванный износ возникает при большом зазоре за -счет ударных на грузок и .нарушения условий смазки. С появлением ударных .на грузок (стуков) идет .направленное ухудшение .качества поверх
ностей, рост шероховатостей и, в |
конечном |
счете, |
разрушается |
|||
подшипник. |
|
т |
(рис. 1.12) при равномерном |
|||
Межремонтный -срок -службы |
||||||
износе сопряженной пары определится по выражению |
||||||
|
|
с |
_ |
с |
|
|
|
|
° т а х |
|
° н |
|
(1.34) |
|
|
2tg а |
|
|||
|
|
|
|
|||
Следовательно, |
для определения |
t необходимо |
знать Д’н — на |
|||
чальный |
зазор по-сле приработки; |
|
— предельно допустимый |
|||
V(i)i |
|
|
|
|
|
|
dVzk) |
|
|
|
|
|
|
^ |
|
|
|
|
|
|
iVM |
|
|
|
|
|
|
Рис. 1.11. |
Реализация |
процесса износа |
Рис. 1.12. |
Изменение зазора в со |
||
в период разрушения сопряжения. |
|
пряжении в зависимости от времени |
его работы.
23
зазор и tg а —величину, характеризующую степень нарастания износа в процессе эксплуатации.
Изучая закономерности нарастания износа на различных ста диях работы детали, можно предупредить появление форсирован
ного |
износа своевременным |
восстановлениемсопряжения. |
|
|||||
/Рассмотрим, каким законам подчиняется время безотказной |
||||||||
работы детали в |
период ее нормальной работы. |
|
|
|
||||
|
|
|
|
На рис. 1.13 приве |
||||
h |
|
|
|
дена |
расчетная схема, |
|||
|
|
|
|
из |
которой |
следует, |
||
|
|
|
|
что |
размер |
предельно |
||
|
|
„Ал |
|
допустимого |
и з но с а |
|||
|
|
|
можно |
разделить |
на |
|||
|
|
А |
|
„/■“ одинаковых поясов |
||||
5 |
|
|
|
„у“, х а р а к т е р и з у ю |
||||
А |
_______ |
|
щих |
собою |
прирост |
|||
4 |
|
|||||||
3 |
7л |
|
|
износа. Н а к л о н н а я |
||||
2 — А^ |
|
|
пунктирная |
линия со |
||||
0 ^ |
г |
|
|
ответствует |
средней |
|||
|
|
величине износа, |
на |
|||||
~ Т+АТ ^ |
|
|
копленного за время t. |
|||||
|
б |
\ |
|
|||||
|
|
|
Случайные |
изменения |
||||
|
|
|
|
величины т) (t) около |
||||
|
|
|
|
этой |
прямой |
обуслов |
||
Рис. 1.13. Расчетная схема износа деталей, |
лены случайностью мо |
|||||||
скачкообразного |
изменения |
износа. |
мента |
возникновения |
||||
|
|
|
Г до Г -f- |
|||||
|
Вероятность |
возникновения скачка за время от |
||||||
+ Д Т будет равна |
|
|
|
|
|
|
||
|
|
Т(Г) = |
Т = ХЛ7. + |
0(Д Т) |
|
|
(1.35) |
|
и ,не будет зависеть от того, |
сколько таких скачков было на участ |
ке от 0 до Т. Как следует из предыдущего рассмотрения, предпо ложение о независимости 'вероятности прироста износа от разме ра ранее /накоплепного износа отвечает зоне нормального износа. Так ли это?
Обозначим через W:(T) размер износа к моменту времени Т.
Тогда величина износа в момент времени Т будет |
|
■n{T)=yW(T) . |
(1,36) |
Прирост износа за время Д Т определится как разность |
|
7](Г+ДГ) — M T) = y [ W { T + L Т ) - W(T)] . |
(1.37) |
(Определим математическое ожидание этой разности. Из урав нения (il.36) следует, что износ детали за время Д Т /может полу чить либо приращение, равное у с вероятностью f , либо прира
24
щение, равное 0 с вероятностью 1 — f. Поэтому |
математическое |
|
ожидание приращения износа |
будет равно |
|
М Ы Г + Д П - т ] ( 7 ' ) } = тз/ + |
(1 -_ у )0 = \ кАТ +0 { АТ) } у . (1.38) |
|
Бели принять ,во внимание, |
что -математическое |
ожидание раз |
ности случайных величин равно разности математических ожида ний, то тол учим
M{- q( T+A Г)}— M h ( 7 ) } |
= [ХД7’+ 0 ( ДГ ) ] у . |
(1.39) |
||
Разделив |
обе части равенства |
на |
А Т и перейдя к |
пределу |
при А Т -> 0, |
получим |
|
|
|
|
dM |т] Г} |
= 1у |
|
(1.40) |
|
dT |
|
|
|
Из зависимости (1.40) следует, что |
Ху от времени не зависит, |
что подтверждает идею о постоянстве средней скорости износа на участке нормальной работы.
Эволюция износа детали (см. рис. |
1.13) |
описывается цепочкой |
|
•$0 |
-> S 2 - > . . . . |
Sr-> |
... , |
где каждое единичное приращение износа происходит по схеме
мпновенного |
повреждения. |
|
|
|
|
Поэтому время безотказной работы должно исчисляться до мо |
|||||
мента получения г-го по счету приращения износа. |
|||||
В этом случае вероятность перехода |
Sk — Sk+ 1 за время А Т |
||||
определяется |
формулой |
|
|
|
|
|
Ч(Т) = \ А Т + 0(АТ) . |
(1.41) |
|||
Причем |
состояние |
износа |
детали |
будет |
характеризоваться |
функциями |
{Рь (Т)ХК — 0,1, 2. ..) , где |
Рк (Т) |
есть не что иное, |
||
как вероятность того, |
что к моменту Т деталь |
будет находиться |
|||
в состоянии |
Sk .Чему равно |
Я* (Г)? |
|
|
Рассмотрим стохастический процесс накопления износа, в ко тором система функций (ЯЙ(Г)) задает распределение некото рой целочисленной случайной величины х( Т), равной количеству
накопленных к моменту Т повреждений |
(в |
нашем |
случае число |
||||||
ступеней износа) на детали. |
Т) — х ( Т) |
на |
интервале времени |
||||||
Приращение |
х ( Т А |
||||||||
от Г до Т 4- А Т |
является |
суммой |
приращений |
х (А Т) — х (Т) |
|||||
и х ( Т 4- А Т) — х (А Т), соответствующих интервалам от Г до |
А 7 |
||||||||
и о т Д 7 ' д о 7 + А 7 ' . |
|
|
того, |
что |
х(Т-{-А Т) — х(А Т) |
||||
Пусть hn{T) — вероятность |
|||||||||
примет значение |
«(где |
« = 1 , 2 , 3 . . . ) . |
х ( Т + А Т ) — х ( АТ) |
w |
|||||
Аналитически |
зависимость |
величин |
|||||||
х ( АТ) — х(Т) может быть выражена |
системой равенств |
|
|||||||
М Т + д Т) = 2 |
hj (А Т) |
(Т) . |
|
(1.42)' |
|||||
|
i=o |
|
|
|
|
|
|
|
25
Единственным распределением, которое удовлетворяет урав нению (1.42), является сложное распределение Пуассона
Р*(Т) = Ч { ! |
6 |
при ЛС>0 . |
(1.43) |
Сумма вероятностей Р0 (Т) + (Т) |
Pr l |
(Т) есть не что |
иное, как вероятность того, что к моменту Г число ступеней износа (накопленных повреждений) равно либо 0, либо 4, либо г — 1, т. е. меньше предельно допустимого г.
. Поэтому время безотказной работы т должно быть не мень ше, чем Т, т. е.
P{ *>T} = % P k(T) . |
(1.44) |
fc= 0 |
|
Подставив вместо Рк(Т) его значение из формулы (1.43), оп ределим величину 1 — Р { т < 7 ’) по формуле
■ |
( 1 ' 4 5 ) |
S - 0
Выражение (1.45) есть функция распределения времени т, подчиняющаяся гамма-распределению при целых значениях г.
Таким образом, схема наносных накапливающихся поврежде-, ■ний может быть описана посредством гамма-раопределения .вре мени безотказной работы х.
В общем виде функция распределения F(Т) гамма-распреде ления задается соотношением
т |
1 |
|
|
F ( T ) = \ f ( T ) d t = - ~ ^ T r- 1 ех 7 |
dt . |
(1.46) |
|
о |
о |
|
|
Плотность же этого распределения при |
Г < 0 равна 0, |
а при |
|
Т > 0 имеет вид |
|
|
|
/( П = гТТ) xr |
е~ХТ |
' |
(1-47) |
где г —число ступеней износа, необходимых для возникнове ния отказа (достижения предельно допустимого изно са);
* Доказательство этого приведено в книге В. Феллер «Введение в теорию вероятностей и ее применение», т. 1, «Мир», 1967, с. 291—293.
26
Г (>) — гамма-функция, |
определяемая |
формулой |
|
||
Г (г) = J |
x r ~l |
е~х |
dt |
, |
(1.48) |
о |
|
|
|
|
|
Для целых г справедливо соотношение |
|
||||
|
Г (г) = |
(г — 1) ! |
|
(1.49) |
|
В наиболее общем виде плотность гамма-распределения может |
|||||
быть записана так |
|
|
|
|
|
___ ^ |
|
|
, |
(1.50) |
|
) |
Г (а 4-1) |
|
|
|
|
где а и р — параметры |
распределения, |
которыми |
могут быть |
||
любые положительные числа. |
|
||||
Математическое ожидание |
М {т} |
и дисперсия гамма-распреде |
|||
ления D {т} равны |
|
|
|
|
|
|
f - |
; |
|
. |
(1.51) |
На рис. 1.14 приведены значения плотности гамма-рашреДе- ления при различных значениях г и X.
Из приведенных кривых можно сделать следующие выводы, имеющие большое практическое значение и объясняющие разно
27
бой в литературе о законах распределения времени безотказной работы:
1.При r= 1 плотность гамма-распределения совпадает с плот ностью экспоненциального распределения. Это хорошо согласует ся с рассмотренной ранее схемой возникновения экспоненциаль ного распределения при мгновенных отказах (единичное повреж дение) .
2.При малых значениях г кривые ассиметричны, а с ростом г они становятся все более симметричными. Другими словами, на блюдается переход от гамсма-распределения к нормальному.
При возрастании г кривая плотности f(T)гамма-распределе
ния,определяемаяуравнением |
(1.47), стремится |
<к |
виду |
|
= в |
. |
(1.52) |
|
г |
|
|
|
U |
|
|
Если в формулу (1.52) ввести значения математического ожи дания нормального распределения параметра С = М {т} и дис персию D {т} = о2, то получим
|
|
|
(Г-С)2 |
|
|
f ( T) |
= |
е |
232 |
. |
(1.53) |
т/2 7Г О |
|
|
|
|
|
Сопоставляя уравнения (1,52) |
и |
(1.53), маходам, что |
|||
C = - f |
, а |
= |
у |
. |
(1.54) |
3. При малых значениях а |
|
по сравнению со средним |
|||
временем безотказной работы |
С — М (т) значения |
плотности f(T) |
близки к 0 на большом интервале времени. Из этого следует, что вероятность появления отказа по причине износа в начале рабо ты детали мала.
Физически этот вывод отражает тот факт, что при большом значении уровня предельного износа и малой величине накоплен ного износа (деталь новая) вероятность отказа мала.
Это обстоятельство должно быть положено ib основу принуди тельных замен деталей, не достигших предельного износа, с целью резкого уменьшения вероятности отказов между ремонтами или техническими обелуживаниям«.
Периодичность таких замен может быть установлена при ис пользовании нормального закона распределения времени безот казной работы (г велико).
28
Здесь же отметим, что (время безотказной работы детали в пе риод приработки подчиняется экспоненциальному (распределению, имеющему максимум плотности при Т=Ю.
4. Разнобой в литературе по вопросу законов распределения времени безотказной работы объясняется там, что различные авто ры начинают исследовать отказы деталей на различных стадиях их износа и получают несравнимые результаты.
ПОРЯДОК ОБРАБОТКИ ЭКСПЕРИМЕНТАЛЬНОГО МАТЕРИАЛА
Для оценки параметров распределений времени безотказной работы могут быть применены методы моментов, квантилей и максимального правдоподобия. Одним из необходимых условий использования этих методов считается определенность выборки
и сведения об отказах элементов, которые подвергаются статисти ческому анализу.
Методика обработки статистических данных о моментах отка зов деталей следующая:
1.Время безотказной работы детали т и функция его распре
деления F (t) рассматриваются в фиксированном интервале рабо ты. Из рассмотрения исключаются те значения т , которые нахо дятся вне границ рассматриваемого интервала времени, т. е.
■с< Тх и т > Ту .
2. Рассматриваемые наблюдения % дают усеченное распреде ление .с функциями распределения
при t < T x F{t) — 0; при £ > 7 у F(t) — 1; при Tx<^t < Ту
Плотность распределения имеет вид
(1.55)
Таким |
образом, функции f (t) и |
F(t) |
задают распределение |
|
Случайной .величины т |
при условии, что ее значения не выходят |
|||
за пределы заданного |
промежутка |
Тх и Ту, т. е. они являются |
||
функциями условных распределений. |
|
|
||
3. |
Определяются |
значения начальных и центральных момен |
||
тов. При |
этом,® первую |
очередь, находятся |
значения первых двух |
29