Рис.1. Гистограмма середины интервалов, кривая распределения плотностейа вероятностей э(l_i), теоретическую кривую распределения т(i) и выравнивающая(огибающая)а кривая.

п.3.1.5. Проверка согласия между эмпирическим и теоретическим (нормальным) законом распределения апо критерию 2 Пирсона :

.) Определим меру расхождения 2 между эмпирическим и теоретическима распределениями:

_r

c²= S (i - n_i`)<sup>2</sup> / n<sub>i</sub>` , где

sup>i=1

n_i аи i` <-- соответствие эмпирической и теоретической частоты попадания случайной величины в

Для добства вычислений критерий 2 определим по формуле:

_{r <_ <_ <_}

c² = N _* D* S [ э(i) - f_т(i) ]² / f_т(i) ,

c² =5,12

б.) Вычислим число степеней свободы i меньше 5-ти объединим с соседними интервалами):

m = r₁ - k - 1, агде

r₁ --а число интервалов полученное при объединении;

k - количество параметров закона распределения.

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

m = 4-2-1 = 1

в.) По значениям 2 и 2) теоретического и эмпирического измерения

2) = P(5,12) = 0,0821; Р(2 ) <> 0,05, значит эмпирическое распределение согласуется с нормальным законом распределения.

п.3.1.6. Определение оценок показателей надёжности детали:

) рассчитаема значение среднего ресурса R при нормальном законе распределения, который численно равен математическому ожиданию а, поэтому R<= = 188,73 (тыс. км)

б) рассчитаем вероятность безотказной работы детали по интервалам наработки по формуле:

_{а<_ <_ rа}

P(l_i) = (N - S i / N),

P(l₁) = (66-3)/66 = 0,95;...

7) =(66-66)/66 = 0

в) построим кривую вероятности безотказной работы детали

i) в зависимости от ее наработки

Рис.2 График

i) кривая вероятности безотказной работы детали в зависимости от наработки

п. 3.2. Расчёт параметров распределения ресурсов детали по корреляционным равнениям долговечности.

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

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

Для деталей двигателя данный подход реализован в виде корреляционных равнений долговечности:

К = А+В(R - С*-1 , где

К<- критерий нагруженности;

, В, С -- коэффициенты;

R -- средний ресурс детали;

n = Т-Т₀=1980-1970=10 - прогнозируемый период (Т- год начала выпуска двигателя, Т₀- 1970 год точка отсчета прогнозируемого периода).

Критерий нагруженности рассчитывается по формуле:

К_к = м^к_*k_т*S_к(

R + 0.1D²_*p_i*b^-_1*r^-1),

средний ресурс рассчитывается равнением: К_к = - 25,2 + 81840 / (R_к - 2,75

k_м^к -- дельный критерий физико-механических свойств кольца;

k_та <-- дельный критерий тепло напряженности;

p_R -- дельное давление на стенку цилиндра от сил пругости кольца Па;

Dа <-- диаметр цилиндра, дм;

p_iа <-- среднее значение индикаторного давления, Па;

bа <-- высота верхнего компрессионного кольца, дм;

r <= 0,5(D -

t <-- радиальная толщина кольца, дм;

S_к -- путь трения кольца, м/км;

lа <-- отношение радиуса кривошипа к длине шатуна;

Sа <-- ход поршня, м;

r а<--а плотность материала кольца, Н/м³.

п.3.2.1. Расчет критерия нагруженности детали двигателя включает следующие этапы:

) Находятся значения сопротивлений дороги Р_yij, воздуха P_wij, разгона

g_ij автомобиля при заданных вариантах дорожно-транспортных условиях эксплуатации:

Р_yij = (G_a + н )iа (H), где <{1}

G_a а<-- сила тяжести снаряженного автомобиля, Н;

q_на <-- номинальная грузоподъемность, Н;

g <-- коэффициент использования грузоподъемности,

y_i -- коэффициент сопротивления движению.

G_a = 15125*9.8 = 148225 (Н),

q_н = 8100*9.8 =79380а (Н),

(79380+148225)*0.04=9104,23175,21587,62964,54552,1

P_wij = (2_aij)/13 (H), где <{2}

k -<- фактор обтекаемости автомобиля, Н*с²/м² ;

F - лобовая площадь автомобиля, м²_;

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

P_gij = i [( M_e<× i_o<× i_wj <×k] (H), где <{3}

k_i а<- коэффициент, учитывающий инерционные нагрузки(междугородние перевозки - i<=0, город и подъездные пути i = 0,2, карьеры i = 0,3);

M_e - максимальный крутящий момент M_e = 700, Н_*м;

i_o<× - передаточное число главной передачи o = 7,22;

i_wj - передаточное число коробки передач в

ТАБЛИЦА 5.

Значения рассчитанных сил сопротивлений дороги, воздуха и разгона.

б) Рассчитываются средние значения эффективного давления

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

P_eij =<ååi j[(1.25r_k<×10^-2)/(V_h<×o<×kij<×т)]<×[(Цk_i)(P_yij+P_wij)+P_gij], где

r_k - динамический радиус колеса, м; на дорогах с твёрдым покрытием r_k < r_ст ;

V_h - рабочий объем цилиндров двигателя, л;

i_o - передаточное число главной передачи;

i_kij - средневзвешенное передаточное число коробки передач;

h_т - к.п.д. трансмиссии автомобиля;

a_i, j - коэффициенты, учитывающие распределение пробега автомобиля по типам дорог

åi = 1 и использование пробега <åj <= 1;

P_yij , P_wij P_gij <- соответственно сопротивления дороги, воздуха и разгона в

i_kij<= 0,6<× V_max(<å j <åi <× V_ij)^-1,где

V_max - максимальная скорость автомобиля, км/ч;

V_ij - средняя скорость автомобиля в

V_ij<=1×V_а1г+2×V_а2г+3×V_а3г)+(1-1×V_а1п+2×V_а2п+3×V_а3п), где

b - коэффициент использования пробега.

V_ij<=0.68*(0.5*25+0.46*35<+0.04*5)+(1-0.68)*(0.5*30+0.46*40<+0.03*10)=30,368, (км/ч),

i_kij <= 0.6*90

0.5*0.68*((1.25*0.488*0.01)/(9*7.22*1.778*0.9))*((1-0.02)*(3175.2+31.108)+04474.033)=0,1519127

Значение средневзвешенного эффективного давления Р_е определяется по формуле :

Р_е = 1×Р_е1г+2×Р_е2г+3×Р_е3г)+(1-1×Р_е1п+2×Р_е2п+3×Р_е3п), где

b <- коэффициент использования пробега;

a_i - коэффициент, учитывающий процент движения автомобиля по типам дорог;

Р_егi, Р_епi - среднее эффективное давление при движении автомобиля в груженом и

порожнем состоянии по различным типам дорог.

0.68*(0.5*0.2080473+0.46*0.09642883+0.04*0.293379)=0.1088789

(1-0.68)*(0.5*0.1257443+0.46*0.04076188+0.04*0.15191270)=0.02806372; P_e=0,1369426 , Па

ТАБЛИЦА 6.

а Таблица рассчитанных значений давления.

S 0,1369426 0,2067211 0,3436637

Для определения Р_м используется зависимость:

P_м = А+В*с_m, где

, В --а коэффициенты, станавливаемые экспериментально;

с_m = (2S <× o<× 0,6V_max)/(60<× 0.377<×r_k)

c_mа -- средняя скорость поршня, м/с;

c_mа <= (2*0.12*7.22*0.6*90)/(60*0.377*0.488) =8,476757, м

P_м = 0.105+0.012*8.476757 =0,2067211 (Па).

Определим среднее индикаторное давление.

Р_i = P_e + P_м , (Па)

Р_i= 0,1369426<+0.2067211<= 0,3436637 Па

п.3.2.2. Рассчитаем значение дельного давления, возникающего от сил пругости компрессионного кольца:

P_R = (0.424*E*A)/[(3--1-1)³], (Па), где

Е - модуль упругости, МПА;

x - постоянная, зависящая от эпюры давления (

- зазор в замке кольца в свободном состоянии.

P_R<=(0.424*1.2*100*0.170)/((3-0.196)*1.2*(1.2*(1/0.05)<-1)*23*23)=0,2112775 Мпа,

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

) гильза - компрессионное кольцо :

К_м^к =(0,2^t_*НВ_к^m_*НВ_гⁿ)/(НВ_к+НВ_г), где

e_0,2^t - коэффициент микрорезания;

НВ_к, НВ_г - соответственно, твердость по Бринеллю кольца и гильзы, ед.;

m, n <- показатели степени, при расчете ресурса кольца принимаются

К_м^к =(1,59*700^1,5_*230²)/(700+230) = 1675008

Удельное значение критерия найдем из соотношения :

k_м = 1/ м

k_м = 1/

п.3.2.4. Оценивается критерий теплонапряженности детали:

К_т = D^0.38_* c_m^0..5 [ (632p_i)/(H_H*i)]^0.88, где

H_H - низшая теплотворная способность топлива, для дизельного топлива H_H<=42496кДж/кга К_т = 0,12^0.38*8,476757^0.5*((632*0,3156)/(42496*0,45))^0.88 = 0,023458596

Определим дельное значение критерия теплонапряженности:

k_т = К_т / К_т^max, где

К_т^max - предельное значение критерия теплонапряженности для рассматриваемой конструкции двигателя:

с_m = (2S*n_e)/60

с_m<= (2*0.12*2600)/60 =10.4;

р_е = [(0.314*e)/V_h ] *10^-2 ,

p_e<= ((0.314*4*700)/9)*0.01 = 0,976, Па

К_т^max <= 0.12^0,38*10.4^0,5*((632*0.976)/(42496*0.45))^0.88 = 0,0702317

k_т = 0.023458596

п.3.2.5. Рассчитаем путь трения компрессионного кольца за один километр пути:

S_т = (100*S*o*i_kij)/(

k),

S_т= (100*0.12*7.22*1.778 )/(3.14*0.488) =100,5312, м/с

На основании рассчитанных параметров определим критерий нагруженности :

К_к = м^к_*k_т*S_к(

R+0.1D²_*p_i*b^-1_*r^-1)

К_к<=0.16067*0.3340*1005.312*(0.2112775<+(0.1*1.2*1.2*0.34367*(1/0.03)*(1/(0.5*(1.2<-0.05)))))<=166,1719

Из корреляционного равнения долговечности:

К_к = -25,2+81840/(R_к-2,75

R_к = 81840 / (К_к + 25,2) + 2.75n

R_к= (81840/(166.1719+25.2) <+13*2.75)<=463,399, (тыс.км).

п.3.2.6. Определим среднеквадратичное отклонение распределения ресурсов детали:

Вычислим коэффициент вариации по корреляционной зависимости

V= 16,507 R^-0,807,

V = 16,507*463.399^-0,807 = 0,1165

среднеквадратичное отклонение вычисляется из соотношения :

s_R =V <× R

s_R <= 0.1165*463.399 =53,98598, (тыс.км)

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

<

f(l) = 1/(R×<Ö 2

*exp.(-(l_i - R_k )2 /2R²)

ТАБЛИЦА 7.

Таблица рассчитанных значений для кривой распределения плотности вероятности

По результатам расчетов построим кривую распределения ресурсов детали по КУД на рис.4..

п.3.3. Расчет параметров распределения ресурсов детали автомобильного двигателя по комбинированному прогнозу.

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

п.3.3.1. Комбинированный прогноз составляется с четом параметров плотности распределения ресурсов, полученных в результате их расчета по КУД и обработки статистических данных распределения ресурсов детали автомобильных двигателей в эксплуатации. Для нормальных законов распределений с параметрами и R (определение по КУД) параметры распределения ресурсов по комбинированному прогнозу определяются следующими зависимостями.

f_S (t) = (2

S)^-0.5 exp(-((t-t_S)² аS)), где

Математическое ожидание определяется по формуле:

t_å = 1*R+ а2*a,

t_å <= 0.5772487*463.399<+0.4227513*188.73 =347,2823, (тыс. км)

Среднеквадратичное отклонение вычисляется по формуле:

D₁ = 2 ; D₂ = R²

D_S <= 1² D₁ + 1²D₂ ;

<

s_å = <Ö 1²*R²+ 2²*2, где

x_1, x₂ - весовые коэффициенты, определяемые по формуле:

x₁=2/(R²+2) ;

x₂=R² / (R²+2);

x₂= 53.98598² / (46.2²+53.98598²) =0,5772487;

x₁= 46.2² / (46.2² +53.98598²) =0,4227513.

D₁ = 46.2² =2134,44 ; D₂ = 53.98598² =2914,486

D_S = 0.5772487² _*2914.486 +0.4227513² _*2134.44=1352,618

<

s_å <=а <Ö 0.5772487² *53.98598² +а 0.422751²* 46.2² = 36.772

Рассчитываем значения для теоретической кривой распределения плотности вероятности нормального закона с параметрами полученными по комбинированному прогнозу и по полученным данным построим кривую н рис 4..

ТАБЛИЦА 8.

Таблица рассчитанных значений для теоретической кривой распределения плотности вероятности.

t(т.км))	84	120	156	192	228	264	300	336	372	408		480	516	552
f(t)	0.00	0.00	0.00	0.1	0.56	0.83	0.00475	0.01084	0.00865	0.00276	0.34	0.17	0.00	0.00

Рис. 4. Графики распределения плотности вероятности.

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

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

Однако лимитирует надежность ЦПГ, как правило, износ верхнего компрессионного кольца по радиальной толщине. Глубина диагностирования определяется ровнем, при котором оценивается значение параметра технического состояния предопределяющего ремонт узла. Для деталей ЦПГ, с четом изложенного, в качестве структурного параметра может быть выбрана площадь зазора в замке верхнего компрессионного кольца (F_2-i).

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

F_{2-i/sub> = F2-0+i}t^a , где

F_2-0 - среднее значение начальной площади в замке компрессионного кольца, мкм²;

g_i - средняя скорость изменения F_{2-i/sub> мкм²/ тыс.км;}

t <- наработка, тыс.км;

a - показатель степени функции изменения параметра.

Для определения доверительных границ используется зависимость среднего квадратического отклонения структурного параметра F2-i от наработки:

s_F2-i² = F2-i²+g_i² t^2a, где

s_F2-i, g_i - среднее квадратическое отклонение F_2-0 и i.

Расчет проводится по следующим этапам.

1.   Определяется значение:

g_i = (F_2-п - F_2-0)/R^a , где

F_2-п - предельное значение структурного параметра, мкм²;

g_i <= ((42.6-9.5)*1)/ 463.399^1.4=а 61.304305

2.   На основании метода линеаризации после преобразования равнений оценивается g_i:

s_gi = [2((i^(2+2a)/a/(F_2-п-F_2-0)^2/a))R²-(i²/(F_2-п-F_2-0)²)F2-0²]^1/2.

s_gi <= (1,4²_* ((61,304^3,43 /331^1,43) _*53,98598²- (61,304² /331²)_*5175²)^0,5 =

=(1.96*((1352342.7 0.5 = 9.846

3.   Находятся доверительные границы изменения параметра, используя в качестве F_2-0,

s_F2-0, i, g_i их оценки:

F_2-i^BH = (F_2-0 bs_F2-0)+( i <
bs_gi)a, где <{26}

F_2-i^B, F_2-i^H - текущие значения верхнего и нижнего доверительных пределов структурного параметра, мкм²;

t_b - статистика Стьюдента для

R²(tsub>1,t₂) =0,8 - нормированная корреляционная функция, деталей ЦПГ;

ТАБЛИЦА 9.

Таблица рассчитанных значений границ изменения параметров.

l (i)	84	120	156	192	228	264	300	336	372	408		480	516	552	588	624
FB, 104	13.97	11.4	18.74	21.52	24.51	27.7	31.07	34.61	38.3	42.13	46.11	50.21	54.44	58.79	63.26	67.83
FH ,104	11.09	12.77	14.67	16.76	19.01	21.41	23.94	26.6	29.37	32.25	35.24	38.33	41.51	44.78	48.13	51.57

<

f <=<Ö 1-R²(tsub>1,t₂) - коэффициент перемешивания реализаций;

На основании расчетов, для 5-6 значений структурного параметра в диапазоне от

t_b F_2-0 до F_2-п изображаются на рис. 5, кривые нижней и верхней границ в таблице 9..

4.Определяются минимальное R^в и максимальное R^н значения ресурса деталей. Для этого в равнение {26} подставляются F_2-i^B= F_2-п , тогда:

R^вн = <{[ F_{2-п а}<-(F_2-оb F2-o)] / (i bg_i)}^1/a, <{27}

R^в = ((42.6*10⁴-(9.5*10⁴+1.96*5175))/ (61.3+0.45*1.96*9.846))^1/1.4= 412.31 , мкм²

R^н <= ((42.6*10⁴-(9.5*10⁴-1.96*5175))/ (61.3+0.45*1.96*9.846))^1/1.4= 430.76, мкм²

Рис.5. Графики верхней и нижней границ изменения параметра.

5. Оценивается наработка до первого ресурсного диагностирования:

tg₁ = R_в - L_тo, где <{28<}

Lтo - периодичность TO-2, станавливается с четом марки и словий эксплуатации автомобиля,

1 <= 430 -12 = 418, т. км

5. Прогнозирование остаточного ресурса детали ЦПГ автомобильного двигателя на основе результатов диагностирования.

Прогнозная оценка остаточного ресурса осуществляется на основе математической модели изменения параметра в функции наработки. Значение структурного параметра при tg_i определяется на основе результатов диагностирования ЦПГ.

п.5.1. Определение структурного параметра на основе результатов диагностирования.

В качестве средства ресурсного диагностирования ЦПГ может быть использован пневмотестер модели К-272. Принципиальная схема измерения площади в замке верхнего компрессионного кольца по величине падения сжатого воздуха, подаваемого в цилиндр, представлена на рис. 6.

Значение структурного параметра рассчитывается нa основании следующей зависимости:

F_2-1 = K(2D

2 / p_i²]p_i))^1/2, где <{29}

К=(1/2)_*F₁/3,13,

К <- коэффициент, учитывающий соотношение коэффициентов сопротивления истечения через входное сопло 1 и зазор кольца 2, также площадь входного сопла (К=0,54Х10⁶ мкм²);

6

5

2 3

1

4

Рис. 6. Принципиальная схема диагностирования ЦПГ пневмотестером модели К-272 :

1 - фильтр; 2, 3 - блок питания; 4 - входное сопло; 5 - измерительный блок; 6 - манометр.

р₂ - атмосферное давление;

Dр = р_о - p_i,

D