Calculation of normal stiffness and contact patch of ultra-low pressure tires for agricultural machinery
- Authors: Balakina E.V.1, Godzhaev Z.A.2, Konshin A.A.1, Kochetov M.S.1
-
Affiliations:
- Volgograd State Technical University
- Federal Scientific Institution Agroengineering Center VIM
- Issue: Vol 90, No 1 (2023)
- Pages: 49-58
- Section: Theory, designing, testing
- Submitted: 24.01.2023
- Accepted: 11.04.2023
- Published: 23.05.2023
- URL: https://journals.eco-vector.com/0321-4443/article/view/133657
- DOI: https://doi.org/10.17816/0321-4443-133657
- ID: 133657
Cite item
Abstract
BACKGROUND: The impact of tires on a deformable ground surface determines not only the ability of a vehicle to move, but also the pressure on ground through the contact patch geometric characteristics. The pressure on soil is normalized and limited by the conditions for the preservation of its agrophysical properties. Therefore, when creating vehicles for various purposes for their movement or movement of their trailers on deformable ground surfaces, it is required to calculate the pressure on ground in the design process. However, currently there are no universal dependencies for calculation of the mentioned parameters of ultra-low pressure tires. This applies to both vertical and inclined tire arrangement.
AIMS: Development and implementation of universal methods for calculation of normal stiffness and con-tact patch parameters of ultra-low-pressure tires for agricultural machinery.
MATERIALS AND METHODS: Software products such as MS Excel, MATLAB and Curve Expert Professional were used in modeling.
RESULTS: The universal analytical-experimental dependence has been obtained to determine the normal stiff-ness of ultra-low pressure tires, including a possible inclination of the wheel rotation axis. The dependence takes into account tire pressure, normal wheel load, an inclination angle of the rotation axis and tire geometric characteristics. The formula is obtained for determining the contact patch contour area of an ultra-low pressure tire. It helps to calculate the current pressure on ground, which should not exceed the specified value.
CONCLUSIONS: The method for calculation of normal stiffness and contact patch parameters such as contact length, contour area and contact pressure for ultra-low pressure tires has been developed and implemented. The method is universal for all ultra-low pressure tires.
Full Text
ВВЕДЕНИЕ
Шины транспортных средств по различному воздействуют на твердую и деформируемую опорные поверхности. Их воздействие на твердую опорную поверхность [1–3] определяет их контактные сцепные и геометрические характеристики. Эти характеристики влияют на безопасность движения транспортного средства (устойчивость, управляемость, тормозные свойства).
Воздействие шин на деформируемую опорную поверхность [4–6] определяет не только возможности машины по передвижению, но и давление на грунт через геометрические характеристики пятен контакта. Давление на грунт нормируется и ограничивается условиями сохранности его агрофизических свойств.
Поэтому, при создании машин различного назначения для передвижения по деформируемым опорным поверхностям или их прицепов требуется в процессе проектирования рассчитывать давление на грунт. Как известно, этот параметр связан не только с нормальной нагрузкой на колесо, но и с контурной площадью пятна контакта шины [7–10]. Чем эта площадь больше, тем ниже давление на грунт. Контурная площадь определяется упругими свойствами шины, а именно, нормальной жесткостью. Для увеличения контурной площади при сохранении прочностных свойств и долговечности шин применяют шины сверхнизкого давления [11, 12]. Для определения длин и контурных площадей пятен контакта при проектном моделировании машины требуются расчетные зависимости. Такие зависимости существуют для легковых и грузовых шин обычного назначения [13–15]. Упругие свойства шин обычных и сверхнизкого давления изучались в работах [16–22]. Однако, на сегодняшний день отсутствуют универсальные зависимости для расчета указанных параметров шин сверхнизкого давления. Это относится и к вертикальному, и к наклонному расположению шин.
ЦЕЛЬ
Целью данного исследования является разработка и реализация универсальных методик расчета нормальной жесткости и параметров пятен контакта шин сверхнизкого давления для сельскохозяйственной техники.
МЕТОДИКА ИССЛЕДОВАНИЯ
С участием авторов проведено экспериментальное определение упругих свойств шин сверхнизкого давления. Теоретическое исследование основано на результатах этих экспериментов. Были измерены 7 шин сверхнизкого давления, фотографии которых показаны на рис. 1.
Рис. 1. Фотографии исследуемых шин.
Fig. 1. Photographs of the studied tires.
Под общим руководством чл.-корр. РАН З.А. Годжаева на экспериментальном стенде определялись нагрузочные характеристики этих шин при разных значениях внутреннего давления. Эти зависимости между нагрузкой и деформацией имели вид пар экспериментальных точек ( ; Z ) при каждом значении давления.
С помощью программного продукта Curve Expert Professional измеренные нагрузочные характеристики были аппроксимированы, и построены расчетно-экспериментальные зависимости вида для каждого значения давления в шине. Для примера, полученные зависимости представлены на рис. 2 в графическом виде.
Рис. 2. Расчетно-экспериментальные нагрузочные характеристики шины 1300х600-533: 1 – р =20 кПа; 2 – р =30 кПа; 3 – р =60 кПа; 4 – р =70 кПа.
Fig. 2. Analytical-experimental load characteristics for the 1300х600-533 tire: 1 – р =20 kPа; 2 – р =30 kPа; 3 – р =60 kPа; 4 – р =70 kPа.
Авторами разработана методика построения универсальной расчетно-экспериментальной зависимости для определения нормальной жесткости произвольной шины сверхнизкого давления. В основе методики лежат исходные экспериментальные нагрузочные характеристикам различных шин сверхнизкого давления. Суть предлагаемой методики можно охарактеризовать следующими этапами:
- Аппроксимация для каждой опытной шины сверхнизкого давления экспериментальных нагрузочных характеристик функциями одинакового вида для каждого давления в данной шине. Здесь , – постоянные коэффициенты.
- Аппроксимация зависимости полученных коэффициентов и от давления в шине для каждой исследуемой шины сверхнизкого давления. В результате получаем дополнительные функции:
и .
Они имеют вид:
; .
В итоге для каждой шины сверхнизкого давления получены выражения:
или .
- По определению, коэффициент нормальной жесткость шины (нормальная жесткость) – есть первая производная от нормальной нагрузки колеса по нормальному прогибу шины. Поэтому, необходимо дифференцирование полученной формулы нормальной нагрузки по нормальному прогибу. Оно позволяет получить выражение для нормальной жесткости каждой шины сверхнизкого давления:
.
Вводя новые обозначения ; получим выражение:
. (1)
Пример расчета нормальной жесткости показан на рис. 3 в графическом виде.
Рис. 3. Нормальная жесткость шины 1300х600-533: 1 – р =20 кПа; 2 – р =30 кПа; 3 – р =60 кПа; 4 – р =70 кПа.
Fig. 3. Normal stiffness for the 1300х600-533 tire: 1 – р =20 kPа; 2 – р =30 kPа; 3 – р =60 kPа; 4 – р =70 kPа.
- Избавление от неизвестной величины в правой части уравнения (1) для каждой шины сверхнизкого давления:
; ;
; ; .
Обозначим
; ; .
В результате выражение (1) примет общий вид:
. (2)
где – давления в шине; – нормальная нагрузка на шину; – постоянные коэффициенты, характеризующие конструкцию данной шины.
Таким образом, удалось избавиться от неизвестной величины нормального прогиба в формуле (1) для расчета коэффициента нормальной жесткости. Теперь в этом выражении присутствуют условия: давление в шине , нормальная нагрузка на шину и постоянные коэффициенты , характеризующие конструкцию данной шины.
В табл. 1 приведены расчетные постоянные коэффициенты для 7 моделей шин сверхнизкого давления.
Таблица 1. Постоянные коэффициенты уравнения (2)
Table 1. Constant coefficients of the equation (2)
Полученные постоянные коэффициенты | № шины | ||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | |
0,829 | 0,738 | 0,650 | 0,675 | 0,642 | 0,695 | 1,152 | |
1,358 | 1,263 | 1,420 | 1,380 | 1,302 | 1,336 | 1,282 | |
1,126 | 0,934 | 0,923 | 0,931 | 0,836 | 0,929 | 1,476 | |
0,358 | 0,263 | 0,420 | 0,380 | 0,302 | 0,336 | 0,282 | |
1,091 | 0,947 | 0,945 | 0,949 | 0,871 | 0,946 | 1,255 | |
0,736 | 0,792 | 0,704 | 0,725 | 0,768 | 0,749 | 0,780 | |
0,264 | 0,208 | 0,296 | 0,275 | 0,232 | 0,251 | 0,220 |
- Получение универсальной расчетно-экспериментальной зависимости для нормальной жесткости шины сверхнизкого давления.
В табл. 2 приведены геометрические характеристики, рассматриваемых шин сверхнизкого давления, и полученные постоянные коэффициенты уравнения (2) для расчета нормальной жесткости.
Таблица 2. Характеристики шины и полученные коэффициенты
Table 2. Tire properties and obtained coefficients
№ Шины | Характеристики шины | Полученные коэффициенты | ||||||
, H | , мм | , мм | , мм | |||||
1 | 4000 | 1300 | 380 | 530 | 0,717 | 1,091 | 0,736 | 0,264 |
2 | 6000 | 1300 | 380 | 600 | 0,633 | 0,947 | 0,792 | 0,208 |
3 | 8000 | 1300 | 345 | 700 | 0,493 | 0,945 | 0,704 | 0,296 |
4 | 8000 | 1300 | 345 | 700 | 0,493 | 0,949 | 0,725 | 0,275 |
5 | 6000 | 1200 | 330 | 600 | 0,550 | 0,871 | 0,768 | 0,232 |
6 | 8000 | 1300 | 345 | 700 | 0,493 | 0,946 | 0,749 | 0,251 |
7 | 8000 | 1300 | 345 | 700 | 0,493 | 1,255 | 0,780 | 0,220 |
Анализ соответствия полученных постоянных коэффициентов разным характеристикам шин выявил взаимосвязи вида:
; ; ,
где – высота профиля шины в мм; – ширина профиля шины в мм.
Аппроксимация указанных зависимостей позволила получить следующую систему универсальных коэффициентов:
(3)
На рис. 4 представлены результаты расчета значений безразмерных коэффициентов по выражениям (3) в графическом виде. Точки показывают значения соответствующих коэффициентов для разных шин из табл. 2.
Рис. 4. Коэффициенты для расчета нормальной жесткости шин сверхнизкого давления: линии – расчет по зависимостям (3); точки – эксперимент.
Fig. 4. Coefficients for calculation of normal stiffness of ultra-low pressure tires: lines – values calculated with the formulae (3); points – experimental values.
По результатам расчета коэффициентов для шины с заданным соотношением появляется возможность вычисления нормальной жесткости шины сверхнизкого давления по выражению (2):
.
Реализация разработанной математической модели выполнялась в средах MatLab и Excel.
Формулы (2) и (3) относятся к случаю вертикального расположения главной плоскости колеса. При его наклонном расположении требуется корректировка уравнения (2):
, (4)
где – коэффициент нормальной жесткости при наклонном расположении; – угол наклона оси вращения колеса; – коэффициента коррекции нормальной жесткости шины.
Для расчета авторами получена приближенная зависимость
. (5)
ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ ИССЛЕДОВАНИЯ
Таким образом, получены расчетно-экспериментальные выражения (2) и (3) для определения нормальной жесткости шин сверхнизкого давления. Формулы (2) и (3) являются универсальными для всех шин сверхнизкого давления. Они учитывают давление в шине, нормальную нагрузку на нее в конкретных условиях и ее геометрические характеристики: высоту и ширину профиля. Полученные коэффициенты являются функциями отношения геометрических характеристик профиля . Зависимостей указанных коэффициентов от других геометрических параметров шины не выявлено. Коэффициенты можно вычислить по выражению (3).
По значению нормальной жесткости можно рассчитать нормальную деформацию и геометрические характеристики пятна контакта шины сверхнизкого давления с грунтом, а также давление на грунт.
, (6)
где – нормальная деформация (прогиб); – нормальная нагрузка.
Длину пятна контакта можно вычислить согласно работе [13] по следующей формуле
, (7)
где – экспериментальный коэффициент уменьшения длины пятна контакта, рассчитанный по формуле Хедэкеля ( ).
Формула (7) справедлива при паспортном давлении в шине или отличающемся от него не более чем на 20%.
Контурную площадь пятна контакта шин больших размеров с грунтовым основанием можно приближенно вычислить по полученной авторами зависимости от паспортной нагрузки по выражению
, (8)
где – коэффициент учета грунтового основания, зависящий от наружного диаметра шины колеса (табл. 3); = 0,335; = -0,0000645; – нормальная нагрузка в кГс (0,1Н); – контурная площадь пятна контакта в см2 (0,0001 м2).
Давление на грунт
. (9)
Воспользовавшись формулами (6), (8) и (9) можно рассчитать давление на грунт. Оно не должно превышать нормируемого значения, ограничиваемого условиями сохранности агрофизических свойств грунта. В случае превышения указанного значения, требуется выбирать шину с меньшей нормальной жесткостью.
ВЫВОДЫ
- Разработана и реализована методика расчета нормальной жесткости и параметров пятна контакта для шин сверхнизкого давления: длины контакта, контурной площади и давления в контакте. Методика является универсальной для всех шин сверхнизкого давления.
- Получена расчетно-экспериментальная зависимость для определения нормальной жесткости шин сверхнизкого давления, в том числе, при возможном наклоне оси вращения колеса. Зависимость учитывает давление в шине, нормальную нагрузку на колесо, угол наклона оси вращения и геометрические характеристики шины: высоту и ширину профиля. В полученном выражении три безразмерных коэффициента являются функциями отношения высоты и ширины профиля шины.
- Получена универсальная расчетно-экспериментальная зависимость для определения контурной площади пятна контакта шины сверхнизкого давления. Зависимость позволяет рассчитывать текущее давление на грунт, которое не должно превышать нормируемой величины.
ДОПОЛНИТЕЛЬНО
Вклад авторов. Е.В. Балакина — разработка концепции, разработка теоретической методики, общее руководство теоретическими исследованиями, подготовка статьи; З.А. Годжаев — разработка концепции, разработка экспериментальной методики, общее руководство экспериментальными исследованиями, подготовка статьи; А.А. Коньшин — обработка результатов экспериментов, теоретические расчеты, подготовка статьи; М.С. Кочетов — участие в получении формулы (5), библиографический поиск, подготовка статьи. Авторы подтверждают соответствие своего авторства международным критериям ICMJE (все авторы внесли существенный вклад в разработку концепции, проведение исследования и подготовку статьи, прочли и одобрили финальную версию перед публикацией).
Конфликт интересов. Авторы декларируют отсутствие явных и потенциальных конфликтов интересов.
Источники финансирования. Исследование выполнено за счет гранта Российского научного фонда № 23-21-00004, https://rscf.ru/project/23-21-00004/.
ADDITIONAL INFORMATION
Authors contribution. E.V. Balakina — concept development, development of a theoretical methodology, general management of theoretical research, writing the text of the manuscript; Z.A. Godzhaev — concept development, development of an experimental methodology, general management of experimental studies, writing the text of the manuscript; A.A. Konshin — processing of experimental results, theoretical calculations, writing the text of the manuscript; M.S. Kochetov — participation in obtaining the formula (5), search for publications, writing the text of the manuscript. All authors made a substantial contribution to the conception of the work, acquisition, analysis, interpretation of data for the work, drafting and revising the work, final approval of the version to be published and agree to be accountable for all aspects of the work.
Competing interests. The authors declare the absence of obvious and potential conflicts of interest.
Funding source. The study was supported by the grant from the Russian Science Foundation № 23-21-00004, https://rscf.ru/en/project/23-21-00004/.
About the authors
Ekaterina V. Balakina
Volgograd State Technical University
Author for correspondence.
Email: fahrgestell@yandex.ru
ORCID iD: 0000-0001-5152-9340
SPIN-code: 2595-2802
Professor of the Technical Operation and Maintenance of Vehicles Department
Russian Federation, VolgogradZahid A. Godzhaev
Federal Scientific Institution Agroengineering Center VIM
Email: fic51@mail.ru
ORCID iD: 0000-0002-1665-3730
SPIN-code: 1892-8405
Corresponding Member of the Russian Academy of Sciences, Dr. Sci. (Tech.), Deputy Director for Innovational and Implemental Activities
Russian Federation, MoscowAlexander A. Konshin
Volgograd State Technical University
Email: alex.tiger.vd.1999@gmail.com
ORCID iD: 0000-0002-7786-5913
SPIN-code: 3503-9285
postgraduate student of the Technical Operation and Maintenance of Vehicles Department
Russian Federation, VolgogradMikhail S. Kochetov
Volgograd State Technical University
Email: kochetov_m.s.1995@mail.ru
ORCID iD: 0000-0003-2527-6890
SPIN-code: 9405-4356
Researcher of the Technical Operation and Maintenance of Vehicles Department
Russian Federation, VolgogradReferences
- Pacejka HB. Tire and Vehicle Dynamics. Waltham: Butterworth-Heinemann; 2012. doi: 10.1016/C2010-0-68548-8
- Tarasik VP. The theory of car movement. Saint Petersburg: BKhV-Peterburg; 2022. (in Russ).
- Balakina EV, Kochetkov AV. Coefficient of adhesion of the tire to the road surface. Moscow: Innnovatsionnoe mashinostroenie; 2017. (in Russ).
- Brunotte J, Lorenz M. Anpassung der Lasteinträge landwirtschaftlicher Maschinen an die Verdichtungsempfindlichkeit von Böden – Wunschtraum oder bereits Realität? In: Tagungsband „Jahr des Bodens“. Würzburg: Bundesministerium für Ernährung und Landwirtschaft, 2015. P. 11–17. Accessed: 24.01.2022. Available from: https://www.lfl.bayern.de/mam/cms07/iab/dateien/jahr-des-bodens-2015_tagung_tagungsband.pdf
- Pryadkin VI, Shapiro VYa, Godzhaev ZA, et al. Transport and technological means on tires of ultra-low pressure. Voronezh: Izd-vo VGLTU im. G.F. Morozova; 2019. (in Russ).
- Godzhaev ZA, Goncharenko SV, Artyomov AV. Calculation and experimental assessment of the impact on the soil of tires of ultra-low pressure of mobile power equipment. Tractors and agricultural machinery. 2020;87(3):35–47. (in Russ). doi: 10.31992/0321-4443-2020-3-35-47
- Revenko VYu, Godzhaev ZA, Rusanov AV. Methods for estimating the contact area of wheel propellers with a support base. Tractors and agricultural machinery. 2019;86(5):48–54. (in Russ). doi: 10.31992/0321-4443-2019-5-48-54
- Leiva-Villacorta F, Vargas-Nordcbeck A, Aguiar-Moya JP, et. al. Influence of Tire Footprint Area and Pressure Distribution on Pavement Responses. In: The Roles of Accelerated Pavement Testing in Pavement Sustainability. Cham: Springer; 2016:685–700. doi: 10.1007/978-3-319-42797-3_45
- Minca C. The determination and analysis of tire contact surface geometric parameters. Review of the Air Force Academy. 2015;13(1(28)):149–154. Accessed: 24.01.2022. Available from: https://www.afahc.ro/ro/revista/2015_1/149.pdf
- Ivanov V. Analysis of tire contact parameters using visual processing. Advances in Tribology. 2010. Vol. 2010. P. 491723. doi: 10.1155/2010/491723
- Kotlyarenko VI, Goncharenko SV, Godzhaev ZA, et al. Ultra-low pressure tire – the optimal propulsion vehicle for vehicles on weakly bearing soils. Tractors and agricultural machinery. 2014;81(2):17–21. (in Russ). doi: 10.17816/0321-4443-65619
- Shapiro VYa, Dmitriev VN, Godzhaev ZA, et al. Study of ultra-low pressure tires for snow and swamp vehicles. Trudy NAMI. 2009;241:123–131. (in Russ).
- Balakina ЕV, Zadvornov VN, Sarbaev DS, et. al. The calculation method of the length of contact of car tires with the road surface. IOP Conf. Ser. Mat. Sci. Eng. 2019;632. doi: 10.1088/1757-899X/632/1/012022
- Khaleghian S, Ghasemalizadeh O, Taheri S. Estimation of the tire contact patch length and normal load using intelligent tires and its application in small ground robot to estimate the tire-road friction. Tire Science and Technology. 2016;44(4):248–261. Accessed: 24.01.2022. Available from: http://hdl.handle.net/10919/75187
- Tomaraee P, Mardani A, Mohebbi A, et al. Relationships among the contact patch length and width, the tire deflection and the rolling resistance of a free-running wheel in a soil bin facility. Spanish Journal of Agricultural Research. 2015;13(2). doi: 10.5424/sjar/2015132-5245
- Goncharenko SV, Godzhaev ZA, Artyomov AV, et al. Elastic characteristics of ultra-low pressure tires. Vertical loads. Avtomobilnaya promyshlennost. 2020;8:18–21. (in Russ).
- Goncharenko SV, Godzhaev ZA, Pryadkin VI, et al. Elastic characteristics of ultra-low pressure tires. lateral and tangential loads. Avtomobilnaya promyshlennost. 2020;10:32–36. (in Russ).
- Kristalny SR, Balakina EV, Popov NV. Friction in contact of a studded wheel with a hard icy bearing surface. Journal of Friction and wear. 2022;43(1):92–104. (in Russ). doi: 10.32864/0202-4977-2022-43-1-92-104
- Balakina EV, Zadvornov VN, Berezovsky MS, et al. Investigation of tire stiffness coefficients. Radial stiffness coefficient. Avtomobilnaya promyshlennost. 2020;8:7–10. (in Russ).
- Balakina EV, Zadvornov VN, Berezovsky MS, et al. Investigation of tire stiffness coefficients. Lateral stiffness coefficient. Avtomobilnaya promyshlennost. 2020;9:18–20. (in Russ).
- Salnikov VI, Barashkov AA, Zadvornov VN, et al. Calculation and experimental universal dependences for determining the radial stiffness of tires. Avtomobilnaya promyshlennost. 2014;7:13–14. (in Russ).
- Zadvornov VN, Balakina EV, Mishchenkov NA. Predicting tread wear from tire stiffness. Journal of Friction and wear. 2020;41(4):485–490. (in Russ).