книги / Методы расчета ресурса работы элементов машин
..pdf
|
|
Находим среднее значение: |
|
|
|
|
|
|
||||||||||
|
|
ki |
|
185,941 |
37,188. |
|
|
|
|
|
|
|||||||
kср |
|
n |
5 |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
Находим итоговое среднее значение: |
|
|
|
|||||||||||||
kср |
|
ki |
|
37,068 37,139 37,171 37,188 37,142. |
|
|||||||||||||
|
|
|
n |
|
|
|
|
|
|
4 |
|
|
|
|
|
|
||
|
|
Находим расчетные значения выходного параметра при вновь |
||||||||||||||||
рассчитанном значении kср , используя формулу |
|
|||||||||||||||||
Tрасч |
|
|
|
kср |
|
|
|
|
37,142 |
|
|
. |
|
|||||
X 0,436 |
0,584 |
0,429a0,438 |
X 0,436 |
0,584 |
0,429a0,438 |
|
||||||||||||
|
|
|
p |
|
ср |
|
|
0 |
|
|
p |
ср |
|
0 |
|
|
||
|
|
Расчетные значения Т, ч, в зависимости от переменных Xp , cp , |
||||||||||||||||
, a0 следующие: |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Tрасч |
|
|
|
|
|
|
|
8236 |
|
|
|
|
5867 |
|
4316 |
|
|
3361 |
2719 |
||||||
|
5133 |
|
|
|
|
3866 |
|
4227 |
|
|
4264 |
3298 |
||||||
|
5104 |
|
|
|
|
5021 |
|
3701 |
|
|
2884 |
3330 |
||||||
|
4572 |
|
|
|
|
4905 |
|
3576 |
|
|
2767 |
2942 |
||||||
|
4419 |
|
|
|
|
3319 |
|
3352 |
|
|
3712 |
2852 |
Вычислим ошибку по формуле
Tэксп Tрасч.
Сведем в таблицу значения ошибок. Значения ошибки следующие:
488 |
96 |
63 |
–11 |
–82 |
438 |
68 |
35 |
–185 |
–112 |
172 |
202 |
101 |
0 |
–50 |
214 |
–4 |
114 |
148 |
–73 |
70 |
4 |
29 |
–8 |
81 |
101
Для того чтобы выбрать из двух рассчитанных разными способами значений k необходимо оценить значения ошибки . Сравнивая значения ошибки, приведенные на этой странице и на с. 97, невозможно однозначно сказать, в каком случае ошибка является меньшей, поэтому для более рациональной оценки необходимо рассчитать отклонение расчетного значения от экспериментального в процентах.
Воспользуемся формулой
100 .
Tэксп
Составим таблицы вычисленных отклонений. |
|
|||||
Значения отклонения , |
%, при k 37,652 следующие: |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4,3 |
0,25 |
|
0,09 |
1,7 |
|
4,51 |
6,59 |
0,38 |
|
0,54 |
5,98 |
|
4,93 |
1,93 |
2,55 |
|
1,31 |
1,39 |
|
2,93 |
3,15 |
1,45 |
|
2,55 |
3,77 |
|
3,94 |
0,2 |
1,26 |
|
0,5 |
1,59 |
|
1,43 |
Значения отклонения , |
%, при k 37,142 следующие: |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5,59 |
1,6 |
|
1,44 |
0,33 |
|
3,11 |
7,86 |
1,73 |
|
0,82 |
4,53 |
|
2,96 |
3,26 |
3,87 |
|
2,66 |
0 |
|
1,52 |
4,47 |
0,08 |
|
3,06 |
5,08 |
|
2,54 |
1,56 |
0,12 |
|
0,86 |
0,22 |
|
2,76 |
Сопоставляя эти значения, делаем вывод, что при расчете коэффициента k первым способом общая ошибка получается меньшей, соответственно, этот способ расчета является более точным и оптимальным.
Можно сказать, что в результате представления расчетной зависимости (35) в виде степенной (48) ошибка хоть и появляется, но является не слишком большой.
102
Глава 6. ОПРОБОВАНИЕ РАСЧЕТНЫХ ЗАВИСИМОСТЕЙ ДЛЯ ОЦЕНКИ РЕСУРСА ЭЛЕМЕНТОВ МАШИН
ПРИ НЕПОСТОЯННОЙ СКОРОСТИ ПРОЦЕССА ИЗНОСА
В проведенных исследованиях мы рассматривали скорость процесса старения как постоянную величину ( ср const ). Это ха-
рактерно для процессов равномерного изнашивания. Однако существуют процессы, скорость которых меняется во времени. Примером могут служить процессы изнашивания поверхностей, когда износ протекает нелинейно.
На рис. 7 приведены графики, характеризующие износ во времени.
По ним мы видим, что процесс со временем постепенно затухает. Если по рис. 7 дать аналитическое описание зависимости износа от времени:
U aT b , |
(49) |
то скорость будет являться производной от этой функции: |
|
dU abT b 1. |
(50) |
dT |
|
Найдем коэффициент a и показатель степени b . Для этого составим систему из двух уравнений износа для наработок в 2000 и 9000 ч. Рассматриваем второй условный номер станка (см. рис. 7):
U1 aT1b ;
(51)
U 4 aT4b .
Прологарифмируем эти уравнения:
lgU1 lg a blgT1 ;
(52)
lgU4 lg a blgT4 .
103
Подставим значения U и T . Для этого воспользуемся табл. 9:
lg11,9 lg a blg 2000; lg34,1 lg a blg9000;
1,075 lg a 3,301b;
1,533 lg a 3,954b.
(53)
(54)
Для того чтобы решить систему, выразим из обоих уравнений lg a и приравняем полученные выражения:
1,075 3,301b 1,533 3,954b
0,653b 0,457 b 0,7.
Подставим полученное значение b в первое уравнение системы (54), выразим lg a и найдем значение a :
lg a 1,075 3,301 0,7 lg a 1,235
a 0,058.
При найденных значениях a и b найдем значения U |
при раз- |
|||||||
ном времени наработки: |
|
|
|
|
|
|
||
|
|
U 0,058T 0,7 . |
|
|
(55) |
|||
Составим таблицу значений: |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
T , ч |
|
2000 |
|
4000 |
|
6000 |
|
9000 |
U , мкм |
|
12 |
|
20 |
|
26 |
|
34 |
Следует |
отметить, что |
значения износа |
для наработок 2000 |
и 9000 ч практически совпадают с подобными значениями табл. 9, полученными экспериментальным путем.
Построим новый график зависимости износа от времени наработки, он будет отличаться от графика на рис. 7, так как мы приняли
104
для износа функцию (49) как наиболее подходящую для описания процесса, а на самом деле ее вид может отличаться и быть более сложным (рис. 24).
Рис. 24. График зависимости износа от наработки при непостоянной скорости
Полученный график практически совпадает с графиком 2 (см. рис. 7). Это означает, что мы подобрали верную функцию для описания процесса износа.
Найдем при разном времени наработки:
0, 058 0, 7T 0,7 1 0, 04T 0,3 .
Составим таблицу значений:
Т, ч |
2000 |
4000 |
6000 |
9000 |
γ, мкм/ч |
0,00416 |
0,00338 |
0,00299 |
0,002 65 |
По графику из рис. 25 видим, что зависимость является нелинейной и с увеличением наработки скорость процесса старения снижается.
Сравним полученные с помощью дифференцирования значениясо значениями, полученными прямым способом (см. табл. 9).
105
Рис. 25. График зависимости скорости процесса старения от наработки
Рассматривая непостоянной, мы получили меньшие значения
скорости процесса в период наработки до 9000 ч (2,65–4,16 мкм на 1000 ч), тогда как в исследованиях ресурса работы при постоянной скорости мы принимали в соответствии с прямым расчетом скорости ср 4,45 мкм на 1000 ч.
ср U2000 / 2000 U9000 |
/ 9000 |
12 / 2000 |
34 / 9000 |
0,004 45 мкм/ч. |
2 |
|
2 |
|
|
Если принять Umax 40 мкм, то при нелинейном изменении износа ресурс, ч, в соответствии с формулой (55)
T 0,7 0,058Umax 0,7 0,05840 11 355.
Если вернуться к расчету при постоянной средней скоростиср 0,004 45 мкм/ч (линейный износ), то ресурс, ч, составляет
T Umax |
|
40 |
8989. |
|
0,004 45 |
||||
|
|
|
106
Зависимость U f (T ) носит линейный характер. Представим ее на графике для более наглядного изображения (рис. 26).
Рис. 26. График зависимости износа от наработки при постоянной скорости
Это говорит о том, что при рассмотрении затухающего процесса старения, при котором увеличение износа происходит нелинейно и после некоторой наработки становится менее значительным, мы получаем большие значения срока службы.
Преобразуем формулу (30) с учетом выражения для скоро-
сти (49):
|
X |
p |
2 2T 2 |
U |
max |
a |
abT b 1 T; |
(56) |
|||||||
|
|
|
a |
|
|
|
|
0 |
|
|
|||||
|
|
|
X |
p |
2 |
2T 2 |
U |
max |
a abT b . |
(57) |
|||||
|
|
|
|
a |
|
|
|
|
|
|
0 |
|
|||
Подставим численные значения a и b : |
|
||||||||||||||
X |
p |
|
2 |
2T 2 U |
max |
a 0,058 0,7T 0,7 |
; |
||||||||
|
|
|
a |
|
|
|
|
|
0 |
|
|
|
|||
|
X |
p |
|
2 2T 2 |
U |
max |
a 0,04T 0,7 . |
(58) |
|||||||
|
|
|
|
a |
|
|
|
|
|
|
0 |
|
|||
Исследуем зависимость ресурса от |
Xp , a0 , a , |
по форму- |
ле (58).
107
Принимаем исходные данные этих параметров из аналогичного расчета как при ср const:
1,5 10 3 мкм/ч;
Umax 40 мкм;
a0 12 мкм;
a 1мкм;
P(t) 0,96.
Для расчета значений T воспользуемся программой Excel, командами «подбор параметра» и «поиск решения».
Преобразуем выражение (56). Возведем обе части уравнения в квадрат и перенесем все слагаемые в одну часть:
Xp2 a2 Xp2 2T 2 Umax a0 2 |
(59) |
||
0,08 Umax a0 T 0,7 0,0016T1,4 |
0. |
||
|
Решение уравнения с помощью подбора параметра и поиска решения основано на подборе (поиске) такого значения T , при котором выполняется равенство (59). Поиск решения используем для проверки численных значений T , полученных подбором параметра.
Например, найдем значение T , используя все исходные данные с помощью команды «подбор параметра»:
1,82 12 1,82 0,00152 T 2 40 12 2
0,08 40 12 T 0,7 0,0016T1,4 0.
Записываем левую часть уравнения (59) в ячейку B1. Вместо T подставляем адрес ячейки А1, в которой впоследствии и получим ответ (значение T , при котором функция, записанная в ячейке В1, обращается в нуль). Запускаем команду «подбор параметра»:
108
Получаем подобранное значение Т.
T 4762,544.
При этом значении T текущее значение функции близко к нулю.
Команда «поиск решения» имеет аналогичный подход. Найдем зависимость T f (P(t)) при abT b 1.
При Xp 1,282 P(t) 0,9 , а T 5683,968 ч 5684 ч. При Xp 1,6 P(t) 0,94 , а T 5078,653 ч 5079 ч. При Xp 1,8 P(t) 0,96 , а T 4762,544 ч 4762 ч. При Xp 2,1 P(t) 0,98 , а T 4358, 279 ч 4358 ч.
При Xp 2,326 P(t) 0,99 , а T 4097,607 ч 4098 ч.
Составим таблицу значений ( abT b 1 ):
109
Р(t) |
0,9 |
0,94 |
0,96 |
0,98 |
0,99 |
Т, ч |
5684 |
5079 |
4762 |
4358 |
4098 |
Исследуем аналогичную зависимость при ср 0,004 45 мкм/ч. Расчет производим по формуле (30), также используя Excel.
При Xp 1,282 |
P(t) 0,9 , а T 4378,307 ч 4378 ч. |
|
||||||
При Xp 1,6 |
P(t) 0,94 , а T 4068,582 ч 4069 ч. |
|
||||||
При Xp 1,8 |
P(t) 0,96 , а T 3894,693 ч 3895 ч. |
|
||||||
При Xp 2,1 P(t) 0,98 , а T 3659,245 ч 3659 ч. |
|
|||||||
При Xp 2,326 |
P(t) 0,99 , а T 3499,236 ч 3499 ч. |
|
||||||
Составим таблицу значений ( ср const ): |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
Р(t) |
|
0,9 |
|
0,94 |
0,96 |
0,98 |
|
0,99 |
Т, ч |
|
4378 |
4069 |
3895 |
3659 |
|
3499 |
Сравнив значения ресурса таблиц на с. 105 и на этой странице (верхняя), можно сказать, что при рассмотрении скорости как переменной величины значения T получаются, как и следовало ожидать, большими, так как процесс износа со временем затухает. По сравнению с прямым расчетом ресурса (54), (55), значения в два раза меньше, что объясняется учетом сложной зависимости (30), (58) от большого числа негативных факторов, снижающих ресурс работы.
Для наглядного подтверждения этих выводов построим график
(рис. 27).
По графику видим, что зависимости носят аналогичный характер, но при изменении вероятности от 0,9 до 0,99 ресурс приabT b 1 изменяется на 27 %, тогда как при ср const на 20 %.
Сопоставляя график зависимости при ср const с зависимо-
стью, полученной при многофакторном эксперименте (см. рис. 20), можем сказать, что они очень похожи по характеру и полученным значениям.
110