MODELING OF THERMAL AND WEAROUT PROCESSES OF BRAND - NEW FRICTION CLUTCH FOR MECHANICAL ENGINEERING

Abstract


The brand-new friction clutch for mechanical engineering, with the complex relative movement of the effective area points on the friction pairs has the following advantages: higher smoothness of turning on and off, lower noise and vibration levels, as well as lower uniform wearout of the effective areas on the friction pairs. This requires reducing the contact temperatures in the zone of the frictional interaction. The in-depth research of the thermal and wearout processes is required. Objective: the development of the complex mathematical model for the thermal and wearout processes of the brand-new friction clutch for mechanical engineering (with the complex relative movement of the effective area points on the friction pairs) by improving the existing mathematical models, considering the kinematic and constructional features of the brand-new friction clutch. Methods Applied: the mathematical modeling of the thermal and wearout processes, the experimental design technique, the laboratory experimentation, the comparative and statistical data analysis of the experimental results. The complex mathematical model, which describes the regularities of the energy release, the heat transfer, and the wearout of the friction clutch with the complex relative movement of the effective areas on the friction pairs, is proposed. The new analytic dependences were obtained: the dependence between the coefficient of sliding friction (for the effective areas on the friction pairs of the friction clutch), the wearout of friction materials, the temperature and the relative speed of the effective area points; the dependence between the thermophysical parameters of friction materials and the temperature. The universal correlation expression for the dependence between the relative speed of the effective area points on the friction pairs of the friction clutch and the heat transfer coefficients for the internal and external convection is proposed. It was found that the results of the mathematical modeling cognate with the experimental results. The developed complex mathematical model stand by the numerical experimentation for the brand-new friction clutches of various designs. The experimental results can be used to obtain the optimal parameters for the friction clutches of various designs.

Full Text

Введение В машиностроении, на транспорте, в легкой промышленности достаточно широко применяются фрикционные муфты различного назначения и конструктивного исполнения. Несмотря на такие различия, фрикционные муфты традиционно используют простое относительное движение рабочих поверхностей фрикционных пар, расположенных соосно [1-3]. Это обстоятельство в ряде случаев ограничивает дальнейшее улучшение эксплуатационных характеристик таких узлов трения. Важную роль в совершенствовании конструкции и повышении эксплуатационных характеристик фрикционных муфт играют современные математические среды и современные математические модели, позволяющие сократить время на разработку и затраты на испытания таких устройств за счет проведения численных экспериментов [4-8]. Повышение качества существующих и получение новых конструкционных и фрикционных материалов для фрикционных муфт также способствует уменьшению массы и габаритов таких устройств, при этом не снижая либо улучшая их эксплуатационные характеристики [9, 10]. Наиболее важными эксплуатационными характеристиками фрикционных муфт, предназначенных для машиностроения, которые имеют резервы для улучшения, являются: уровень вибрации и шума, износ рабочих поверхностей фрикционных пар, плавность включения и выключения, условия охлаждения и удобство технического обслуживания [11, 12]. Дальнейшее развитие техники и технологии, повышение требований надежности, малошумности, ресурсосбережения, внедрение новых, перспективных, конструкционных и фрикционных материалов способствует усовершенствованию фрикционных муфт для машиностроения. Фрикционные муфты нового типа, предназначенные для машиностроения, которые основаны на использовании сложного относительного движения рабочих поверхностей, при некотором усложнении их конструкции, по сравнению с фрикционными муфтами традиционной конструкции, обладают такими важными преимуществами, как малошумность и низкий уровень вибрации при включении и выключении, высокая плавность включения и выключения, сравнительно меньший и более равномерный износ рабочих поверхностей фрикционных пар таких муфт [13]. Уменьшение линейного износа рабочих поверхностей фрикционных элементов муфты со сложным относительным движением рабочих поверхностей в значительной степени определяется сравнительно более низкими контактными температурами в зоне фрикционного взаимодействия. Исходя из этого, при разработке конкретных конструкций фрикционных муфт нового типа для решения задач машиностроения, в частности станкостроения, важным этапом является проведение численного эксперимента, экономящего ресурсы и время, необходимое для разработки фрикционной муфты конкретного назначения. Для проведения численного эксперимента необходима адекватная математическая модель, описывающая тепловые процессы и процессы изнашивания фрикционной муфты нового типа со сложным относительным движением точек рабочих поверхностей фрикционных пар. Существуют математические модели, предназначенные для описания тепловых процессов, возникающих при простом относительном движении точек рабочих поверхностей фрикционных пар. В качестве основы для разработки комплексной математической модели может быть взята известная математическая модель [14]. Разработка комплексной математической модели тепловых процессов и процессов изнашивания применительно к сложному относительному движению рабочих поверхностей фрикционных пар фрикционной муфты нового типа является актуальной научно-технической задачей. Целью исследования является разработка комплексной математической модели тепловых процессов и процессов изнашивания фрикционной муфты нового типа со сложным относительным движением точек рабочих поверхностей фрикционных пар путем усовершенствования существующих математических моделей, с учетом кинематических и конструкционных особенностей фрикционной муфты нового типа. Полученные результаты и их обсуждение В предложенной автором фрикционной муфте нового типа для машиностроения передача вращающего момента от ведущего диска к ведомому диску фрикционной муфты осуществляется посредством фрикционных дисков, оси вращения которых связаны с ведомым диском [13]. При включении фрикционной муфты фрикционные диски начинают вращаться, участвуя в двух вращательных движениях одновременно. В результате этого с момента включения муфты и до полного сцепления, а также при выключении фрикционной муфты, возникает сложное относительное движение точек рабочих поверхностей фрикционных дисков и ведущего диска фрикционной муфты нового типа, которое обеспечивает плавность включения и выключения, малошумность, низкий уровень вибрации, равномерный и сравнительно меньший износ рабочих поверхностей фрикционных пар по сравнению с фрикционной муфтой традиционной конструкции, использующей простое относительное движение точек рабочих поверхностей фрикционной пары. Внешний вид одного из конструктивных исполнений фрикционной муфты нового типа показан на рисунке. Комплексная математическая модель тепловых процессов и процессов изнашивания разработана на основе известной математической модели для описания тепловых процессов при простом относительном движении рабочих поверхностей фрикционной пары [14]. Исходные данные для моделирования приведены в таблице. Дополнительно к ним использованы аппроксимационные выражения для зависимости теплофизических величин фрикционных и конструкционных материалов от температуры T, которые получены автором по справочным данным [15-18]. Рис. Опытный образец фрикционной муфты: 1 - ведущий вал; 2 - ведущий диск; 3 - ведомый (прижимной) диск; 4 - ведомый вал; 5 - ось одного из шести фрикционных дисков; 6 - один из пяти фрикционных дисков Для материала ведущего диска и ведущего вала фрикционной муфты: - коэффициент теплопроводности, Вт/(м·°С), - плотность, кг/м3, - теплоемкость, Дж/(кг·°С), - твердость, Па, - модуль Юнга, Па, - коэффициент Пуассона, Для материала фрикционных дисков муфты: - коэффициент теплопроводности, Вт/(м·°С), Исходные данные к моделированию тепловых процессов и процессов изнашивания фрикционной муфты Параметр Обозначение Значение Размерность Число срабатываний фрикционной муфты в единицу времени n0 2 с-1 Шаг по времени τ 1 с Толщина стенки ведущего диска d 0,003 м Толщина фрикционного диска 0,005 м Безразмерная координата по толщине стенки диска X 0,5 - Первый корень характеристического уравнения η 3,381 71 - Температура ведущего диска в начальный момент моделирования 20 °С Температура окружающей воздушной среды 20 °С Коэффициент черноты материала ведущего диска муфты 0,3 - Коэффициент черноты материала фрикционного диска 0,6 - Характеристический размер вентиляционных каналов между элементами фрикционной муфты 0,006 м Характеристический размер каналов между фрикционным и прижимным дисками 0,001 м Радиус ведущего диска 0,085 м Средний радиус по прижимному диску 0,06 м Характеристический размер для ведущего диска 0,17 м Характеристический размер для фрикционных дисков 0,048 м Коэффициент трения скольжения фрикционной пары µ 0,3 - Коэффициент влияния скорости на коэффициент трения 0,184 - Коэффициент влияния скорости на коэффициент трения 0,1 с/м Коэффициент влияния температуры на коэффициент трения 0,105 - Коэффициент влияния температуры на коэффициент трения 0,014 1/°С Объем ведущего диска 68∙10-6 м3 Площадь ведущего диска 26,7∙10-3 м2 Площадь фрикционной поверхности фрикционного диска 1,8∙10-3 м2 Открытая поверхность фрикционного диска 2,1∙10-3 м2 Объем фрикционного диска 9∙10-6 м3 Коэффициент пропорциональности в формуле износа 1∙10-6 - Средняя мощность электродвигателя P 11∙103 Вт Момент инерции приводимого во вращение механизма J 0,125 кг∙м2 - плотность, кг/м3, - теплоемкость, Дж/(кг·°С), - твердость, Па, Внешняя среда: - теплоемкость воздуха, Дж/(кг·°С), - коэффициент теплопроводности воздуха, Вт/(м·°С), - динамическая вязкость воздуха, Па·с, - кинематическая вязкость воздуха, м2/с, Разработанная комплексная математическая модель тепловых процессов имеет следующий вид. Удельная изобарная теплоемкость материала ведущего диска, Дж/(м3·°С), Температуропроводность материала ведущего диска, м2/с, Для расчета средней скорости относительного движения фрикционных поверхностей ведущего диска и фрикционного диска, совершающего вращение относительно поверхности ведущего диска, в момент включения муфты автором предложено выражение где - шаг разбиений по осям x и y, м; r - радиус фрикционного диска, м; Rd - радиус ведущего диска, м; - угловые скорости вращения ведущего диска и фрикционных дисков, с-1; и - координаты элемента площади фрикционного взаимодействия в системе координат, связанной с осью вращения фрикционного диска муфты, м. Зависимость коэффициента трения от средней скорости относительного движения фрикционных поверхностей ведущего и фрикционного дисков и от температуры [19] где T(t) - средняя температура в зоне фрикционного контакта, °С. Работа фрикционного взаимодействия при включении муфты в интервале времени где - время включения без учета сил сопротивления движению, Дж, где P - средняя мощность электродвигателя, Вт; J - момент инерции приводимого во вращение механизма, кг∙м2; - угловая скорость вращения механизма, с-1, Доля работы фрикционного взаимодействия, переходящая в нагрев ведущего диска [20], Плотность теплового потока из зоны фрикционного взаимодействия в ведущий диск фрикционной муфты, Вт/ , Безразмерное время Безразмерная относительная температура [9] Число Рейнольдса для внешней поверхности ведущего диска Число Рейнольдса для внутренней поверхности ведущего диска Число Прандта В расчетах используются корреляционные соотношения между характеристической длиной внутренних каналов между элементами муфты и скоростью воздушного потока для внутренней конвекции [21]. Коэффициент теплопередачи, Вт / (м2∙°С), Для внешней конвекции [22] коэффициент теплопередачи, Вт / (м2∙°С), Внутренний конвективный тепловой поток в направлении нормали к поверхности ведущего диска, Вт/ м2, Внешний конвективный тепловой поток в направлении нормали к поверхности ведущего диска, Вт/ м2, Тепловой поток в направлении ведущей оси муфты, Вт/м2, Поток мощности теплового излучения с поверхности ведущего диска, Вт/м2, где σ - постоянная Стефана - Больцмана, кг/ Тепловой поток, идущий на нагрев ведущего диска, Вт/м2, Температура ведущего диска в момент времени , °С, Переход от безразмерного времени к t, с, Энергия, рассеиваемая внутренней конвекцией ведущего диска за время Дж, Энергия, рассеиваемая внешней конвекцией ведущего диска за время Дж, Энергия, рассеиваемая ведущей осью муфты за время Дж, где - площадь поперечного сечения ведущей оси, м2. Энергия, рассеиваемая тепловым излучением с поверхности ведущего диска за время Дж, Суммарная энергия, рассеиваемая внутренней конвекцией ведущего диска за время включения муфты, Дж, где N - число разбиений временного интервала. Суммарная энергия, рассеиваемая внешней конвекцией ведущего диска за время включения муфты, Дж, Суммарная энергия, рассеиваемая ведущей осью муфты за время включения муфты, Дж, Суммарная энергия, рассеиваемая тепловым излучением с поверхности ведущего диска за время включения муфты, Дж, Из закона Арчарда [23] получим выражение для линейного износа фрикционной поверхности ведущего диска за время включения муфты м: где - коэффициент износа, экспериментальное значение которого для разных материалов фрикционных пар находится в пределах - нормальное давление, Па; - средняя скорость относительного движения точек фрикционных поверхностей в интервале времени от t до м/с. Получим закон износа в удобной для проведения сравнительного анализа форме. Для этого выразим нормальное давление, действующее в зоне фрикционного взаимодействия, через силу прижатия фрикционного диска к ведущему диску и через площадь контактной поверхности фрикционного диска: где - сила прижатия фрикционного диска к ведущему диску, Н. Перепишем выражение для линейного износа в виде Скорость линейного износа, м/с, Суммарный линейный износ фрикционной поверхности ведущего диска за время включения муфты, мм, Суммарный объемный износ фрикционной поверхности ведущего диска за время включения муфты, см3, Суммарный объемный износ фрикционной поверхности ведущего диска за время включения муфты, см3/кВт, Работа фрикционного взаимодействия, которая идет на износ фрикционной поверхности ведущего диска за время включения муфты, Дж, где - эмпирический коэффициент пропорциональности между работой фрикционного взаимодействия, которая идет на износ, и величиной линейного износа. Для различных фрикционных материалов Дж/м. Суммарные диссипативные потери энергии ведущего диска за время включения муфты, Дж, Математическая модель содержит аналитические зависимости для описания тепловых процессов и процесса износа фрикционных дисков, аналогичные приведенным зависимостям для ведущего диска. Доля работы фрикционного взаимодействия, переходящая в нагрев каждого фрикционного диска, Температура фрикционного диска в момент времени t с учетом объема фрикционного диска, °С, Полученная математическая модель была дополнительно усовершенствована путем замены корреляционных соотношений между характеристическими размерами ведущего и фрикционного дисков и скоростью воздушного потока, определяющими коэффициенты теплопередачи и следующими полученными автором эмпирическими корреляционными соотношениями: где соответствует ведущему диску, - фрикционному диску; для внутренней конвекции, для внешней конвекции; - число слагаемых в корреляционном соотношении; - эмпирические коэффициенты, зависящие от конструкции фрикционной муфты. Выводы Результаты численных экспериментов, проведенных с применением математического пакета MATLAB, показали, что при одинаковых начальных условиях при многократном циклическом включении и выключении фрикционных муфт, предназначенных для машиностроения, традиционной и новой конструкций наблюдаются различия как тепловых процессов, так и процессов изнашивания. Фрикционная муфта новой конструкции может быть охарактеризована более плавным ростом контактных температур в зоне фрикционного взаимодействия до существенно, в 2-3 раза, меньших максимальных значений, что также способствует уменьшению скорости линейного износа, линейного и объемного износа рабочих поверхностей ведущего диска и фрикционных дисков муфты нового типа, согласно зависимости твердости фрикционного материала от контактной температуры. Понижение температуры приводит к реализации более высоких значений коэффициента трения скольжения, в соответствии с зависимостью коэффициента трения скольжения от температуры, что улучшает работу фрикционной муфты. Наличие воздушных каналов между конструктивными элементами фрикционной муфты нового типа, а также вращение вокруг своих осей фрикционных дисков в интервалы включения и выключения муфты способствует лучшей внутренней и внешней конвекции ведущего диска и фрикционных дисков, что также способствует снижению рабочих значений температуры. Полученные в ходе исследования аппроксимационные зависимости будут полезны при проектировании и постановке численных экспериментов работы узлов и механизмов, использующих в своей работе фрикционные пары, например таких, как фрикционные приводы, преобразующие вращение в линейное перемещение и наоборот, а также дисковые тормоза для машиностроения, включая станкостроение. В указанных узлах трения и фрикционных механизмах также может быть весьма эффективно применено сложное относительное движение рабочих поверхностей фрикционных пар.

About the authors

V. A Voytenko

Vladimir Dahl Lugansk State University

References

  1. Бабаджанов С.Х., Дустова Ф.Х. Фрикционные муфты на швейных машинах, процессы сцепления и расцепления // Молодежь и наука: шаг к успеху: сб. науч. ст. 4-й Всерос. науч. конф. перспективных разработок молодых ученых, г. Курск, 19-20 марта 2020 г. - Курск: Изд-во Юго-Запад. гос. ун-та, 2020. - Т. 5. - С. 28-32.
  2. Friction behavior of pre-damaged wet-running multi-plate clutches in an endurance test / T. Schneider, K. Voelkel, H. Pflaum, K. Stahl // Lubricants. - 2020. - Vol. 8(7). - P. 68. doi: 10.3390/lubricants8070068
  3. Thermal stress analysis of dry friction clutches / O.I. Abdullah, J. Schlattmann, H. Jobair, N.E. Beliardouh, H. Kaleli // Industrial Lubrication and Tribology. - 2020. - Vol. 72, no. 2. - P. 189-194. doi: 10.1108/ILT-05-2018-0176
  4. Static and dynamic analysis of single plate clutch in four wheeler application using ANSYS / P. Viswabharathy, G. Vigneshwar, M. Pragadhishwaran, M. Gopalakrishnan // Int. J. of Emerging Techn. in Eng. Res. - 2017. - Vol. 5, iss. 3. - P. 222-233.
  5. Gouse seema begum S., Santosh Priya K., Noor Mahammed S. Design and transient analysis of friction clutch plate for two wheeler by using ANSYS // Int. J. of Adv. Sci. and Techn. - 2019. - Vol. 28(19). - P. 181-187. - URL: http://sersc.org/journals/index.php/IJAST/article/view/ 2510 (accessed 7 August 2021).
  6. Barve N.A., Kirkire M.S. Analysis of single plate friction clutch using finite element method // Int. J. of Adv. Sci. Res. and Eng. Trends. - 2017. - Vol. 2, iss. 11. - P. 273-276.
  7. Parameters distribution characteristics of material fatigue life based on improved bootstrap method / M. Zhang, X. Liu, Y. Wang, X. Wang // Int. J. Damage Mech. - 2019. - Vol. 28(5). - P. 772-793. doi: 10.1177/1056789518792658
  8. Modepalli G. Clutch assembly modeling and dynamic analysis // Int. J. of Mech. and Prod. Eng. - 2015. - Vol. 3, iss. 11. - P. 102-107.
  9. Harish K., Kumar Y.D. Optimization of friction clutch for various friction materials // Int. J. of Res. - 2017. - Vol. 4, no. 17. - P. 3493-3498.
  10. Dhengre N., Mogra A., Gupta A. Investigating behaviour of multi-clutch plate frictional materials using ANSYS // Int. J. of Appl. Eng. Res. - 2018. - Vol. 13, no. 16. - P. 12657-12662.
  11. Life cycle prediction and evaluation of clutch friction plate considering wear models and thermal stress / M. Zhu, X. Liu, F. Kan, Z. You //j. of Tribol. - 2021. - Vol. 143(4). - P. 041701. doi: 10.1115/1.4048383
  12. Cakmak T., Kilic M. Estimation of friction surface temperature of a dry clutch // Int. J. of Heavy Vehicle Sys. - 2020. - Vol. 27, no. 5. - P. 683-702. doi: 10.1504/IJHVS.2020.111259
  13. Войтенко В.А. Плавный пуск дискового режущего инструмента за счет применения фрикционной муфты со сложным относительным движением поверхностей фрикционной пары // Ресурсосберегающие технологии производства и обработки давлением материалов в машиностроении. - 2021. - № 2 (35). - С. 39-49.
  14. Heat Transfer Module User's Guide. Version: 5.4. 2018.COMSOL. - URL: https://doc.comsol.com/5.4/doc/com.comsol.help.heat/HeatTransferModuleUsersGuide.pdf (accessed 01 July 2021).
  15. Чиркин В.С. Теплофизические свойства материалов ядерной техники: справ. - М.: Атомиздат, 1968. - 485 с.
  16. Теплофизические свойства веществ [Электронный ресурс]: справ. - URL: http://temperatures.ru/pages/ teplofizicheskie_svoistva_veschestv (дата обращения: 01.07.2021).
  17. Варгафтик Н.Б. Справочник по теплофизическим свойствам газов и жидкостей. - М.: Наука, 1972. - 721 с.
  18. Чиркин В.С. Теплофизические свойства материалов: справ. рук-во. - М.: Физматгиз, 1959. - 356 с.
  19. Lee N.J., Kang C.G. The effect of a variable disc pad friction coefficient for the mechanical brake system of a railway vehicle // PLoS ONE. - 2015. - Vol. 10(8). - e0135459. - URL: https://doi.org/10.1371/journal.pone. 0135459 (accessed 5 July 2021).
  20. Сахин В.В. Теплообмен в однородной среде (теплопередача) [Электронный ресурс]: учеб. пособие // БГТУ «Военмех» им. Д.Ф. Устинова. - СПб.: Лань, 2017. - 121 с. - URL: https://e.lanbook.com/book/121813 (дата обращения: 01.07.2021).
  21. Chen J., Gao F. Temperature field and thermal stress analyses of high-speed train brake disc under pad variations // The Open Mech. Eng. J. - 2015. - No. 9. - P. 371-378.
  22. Codina R.Comparison of some finite element methods for solving the diffusion-convection-reaction equation // Comp. Meth. Appl. Mech. Engrg. - 1998. - Vol. 156. - P. 185-210.
  23. Archard J.F. Wear theory and mechanisms. Wear Control Handbook. - N.Y.: ASME, 1980. - 80 p.

Statistics

Views

Abstract - 92

PDF (Russian) - 49

Refbacks

  • There are currently no refbacks.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies