яді:
.
В
Рис. 2 Визначення сумарної реактивної сили
Тоді
В В
Вираз для визначення сумарної реактивної сили має вигляд:
.
Вираз для визначення сумарної сили опору обертанню має вигляд:
В
Приведений коефіцієнт тертя ковзання дорівнює:
.
Представлення Pi у вигляді вираженні Pi = P2 (li) В· P1 (?) вимагає визначення вхідних у нього компонентів P2 (li) і P1 (?).
P1 (?) найчастіше може бути описано залежностями P1 (?) = const (рис. 1 а) або P1 (?) = cos? (Рис. 1 б). p> Величина P2 (li) В· визначається видом епюри розподілу тиску в поздовжньому напрямку.
Оскільки Qni = Qn (li)-зосереджена сила, діюча в i-му поперечному перерізі системи вал-підшипник, те:
В
Звідки
В
2.2 Математична модель чисельного інтегрування методом трапеції
В
Рис.3 Метод трапеції
Інтеграл оцінюється обчисленням суми площ елементарних трапецій зі сторонами, дорівнювати значенням на початку і наприкінці елементарного відрізка. Це наближення рівносильно заміні функції відрізком прямої, що з'єднує значення в початковій і кінцевій точках відрізку (рис.3). p> Площа кожного елементарного сегмента розбиття вважається за формулою:
В
де
В
Просуммируем елементарні площі:
В
2.3 Математична модель чисельного інтегрування методом Сімпсона (парабол)
Більш високу точність розрахунків забезпечує використання параболічної (квадратичної) інтерполяції за трьома сусіднім точкам відрізка. Рівняння полінома другого ступеня, що проходить через точки можна записати у вигляді:
В
Проінтегрував формулу з урахуванням того, що отримаємо - площа під параболою на відрізку
В В
Підсумувавши всі елементарні площі, отримаємо:
В
де
В
3. Алгоритм
3.1 Алгоритм головний програми
алгоритм програма знос підшипник тертя
1. Введення вихідних даних
.1 Введення зосередженої сили Q
.2 Введення коефіцієнта тертя f
.3 Введення граничного кута контакту ?
.4 Введення числа обертів вала m
.5 Введення радіусу вала r
.6 Введення кутової швидкості обертання вала w
.7 Введення довжини підшипника Ln
.8 Введення мінімального значення сили Qmin
.9 Введення числа розбиттів n
.10 Запис функції P 1 (