applied sciences
Article Assessment of Tire Features for Modeling Vehicle Stability in Case of Vertical Road Excitation
Vaidas Lukoševiˇcius*, Rolandas Makaras and Andrius Dargužis
Faculty of Mechanical Engineering and Design, Kaunas University of Technology, Studentu˛Str. 56, 44249 Kaunas, Lithuania; [email protected] (R.M.); [email protected] (A.D.) * Correspondence: [email protected]
Abstract: Two trends could be observed in the evolution of road transport. First, with the traffic becoming increasingly intensive, the motor road infrastructure is developed; more advanced, greater quality, and more durable materials are used; and pavement laying and repair techniques are im- proved continuously. The continued growth in the number of vehicles on the road is accompanied by the ongoing improvement of the vehicle design with the view towards greater vehicle controllability as the key traffic safety factor. The change has covered a series of vehicle systems. The tire structure and materials used are subject to continuous improvements in order to provide the maximum possi- ble grip with the road pavement. New solutions in the improvement of the suspension and driving systems are explored. Nonetheless, inevitable controversies have been encountered, primarily, in the efforts to combine riding comfort and vehicle controllability. Practice shows that these systems perform to a satisfactory degree only on good quality roads, as they have been designed specifically
for the latter. This could be the cause of the more complicated car control and accidents on the lower-quality roads. Road ruts and local unevenness that impair car stability and traffic safety are not
Citation: Lukoševiˇcius,V.; Makaras, avoided even on the trunk roads. In this work, we investigated the conditions for directional stability, R.; Dargužis, A. Assessment of Tire the influence of road and vehicle parameters on the directional stability of the vehicle, and developed Features for Modeling Vehicle recommendations for the road and vehicle control systems to combine to ensure traffic safety. We Stability in Case of Vertical Road have developed a refined dynamic model of vehicle stability that evaluates the influence of tire tread Excitation. Appl. Sci. 2021, 11, 6608. and suspensions. The obtained results allow a more accurate assessment of the impact of the road https://doi.org/10.3390/app11146608 roughness and vehicle suspension and body movements on vehicle stability and the development of recommendations for the safe movement down the road of known characteristics. Academic Editors: Flavio Farroni, Andrea Genovese and Keywords: directional stability; road profile; road unevenness; vehicle-road interaction; vertical Aleksandr Sakhnevych vehicle excitation; tire models; tire tread
Received: 29 June 2021 Accepted: 16 July 2021 Published: 18 July 2021 1. Introduction
Publisher’s Note: MDPI stays neutral The vehicle-road subsystem includes various aspects of the vehicle-road interac- with regard to jurisdictional claims in tion. This subsystem is emphasized where the technical part of transport operations is published maps and institutional affil- investigated. Considerable emphasis has already been placed on the description of the iations. vehicle-road interaction, vehicle movement under various conditions, and substantiation of the requirements for vehicles both in the theoretical studies and in the experimental investigations where the theory has not provided any viable solutions [1–3]. Basic solutions could mistakenly be considered to have been found, with further inquiry allegedly being limited to observations that are necessary to correct the solutions. Nonetheless, certain Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. limitations have been introduced increasingly by a lot of countries in response to the This article is an open access article aggressive development of this type of transport. These limitations have an effect both distributed under the terms and on individual components of the subsystem and on the subsystem itself and promote the conditions of the Creative Commons improvement of both technical and organizational aspects thereof. The road and vehicle Attribution (CC BY) license (https:// as the elements of the vehicle-road subsystem may be described individually as they are creativecommons.org/licenses/by/ viewed and may be explored as individual elements of the system in an investigation of 4.0/). transport systems [4].
Appl. Sci. 2021, 11, 6608. https://doi.org/10.3390/app11146608 https://www.mdpi.com/journal/applsci Appl. Sci. 2021, 11, 6608 2 of 34
Pavement unevenness values are not completely random. An attempt was made to systematize these values. A method enabling assessment of the unevenness height dispersion and correlation function that has been used for a fairly long time is already more appropriate for the modeling [5–8]. The specifics of unevenness lie in that certain patterns between the unevenness height and unevenness wavelength become established during the pavement laying process and use. Hence, the random process model cannot be applied to the road description, as waviness and ruts are registered in the road pavement, in particular, after the use [9–12]. In the stability tasks, two types of unevenness have the greatest effect: road micro- profile that causes excitation of the suspension operation [13,14] and road roughness due to its relation to the grip ratio [15–17]. The vehicle-road interaction is determined by two intermediate links: the tire and the suspension [18,19]. Vehicle stability modeling tasks involve multiple parameters and need to be split into several individual tasks during the analysis. Nonetheless, a more detailed inquiry has shown that the tasks are interrelated, and integrated analysis is needed when dealing with more complex tasks. Analysis of the evolution of trends in traffic accidents and vehicle design has shown that controllability is one of the key operational properties of a vehicle that influences safety. Vehicle controllability is the ability to move following the trajectory designed by the driver [20,21]. The driver sets the desired direction using the steering wheel mechanism, and the controllability tasks are focused on the steering wheel mechanism [22,23]. Deeper inquiry into the potential excitation caused by road unevenness has shown that it is not enough to know the slip characteristics, as the wheel direction vector is influenced by the wheel spatial position, i.e., camber, toe alignment, etc. The dynamic equations cannot be used directly due to the specifics of the tasks and to the fact that the characteristics of individual elements are clearly nonlinear. Semi-empirical formulae are used to describe the effect of the change of the forces acting on the wheel and the spatial position thereof on the direction vector. Dependencies of the link between the slip angle and lateral force, lateral force and relative wheel slippage, slip angle, and stabilizing moment are required in order to be able to perform the modeling. During the modeling, the dependencies are expressed by either linear or higher-degree dependencies [24–30]. The majority of the currently available methodologies employ simplified models. These models individually assess the kinematics of the steering wheel mechanism, sus- pension kinematics, stiffness of the suspension elements, and forces acting on the vehicle while making a turn [31–35]. This kind of analysis does not enable the identification of the key factors. Hence, the models consisting of multiple elements, the characteristics of which can be known only after specific vehicle design details are available, need to be used for the integrated analysis. This complicates the integrated analysis and generalization in the investigation of directional stability viewed as the criterion for assessment of the technical condition of the vehicle and the condition of the pavement. Hence, the focus was first turned to another driver-vehicle-road system element, i.e., road, prior to designing the integrated directional stability model for the vehicle and performing the analysis. In this context, the road is associated with the effects which may cause vehicle instability, while the purpose of improvement of stability is related to the solution of the vehicle-road interaction. A moving vehicle is subject to external effects and reacts to the driver’s actions (Figure1 ). In the investigation of vehicle directional stability, all these factors need to be considered. However, the vehicle model becomes very complex, making the analysis of the effect of individual factors more difficult. The present paper does not account for the driver’s actions, environmental impact, aerodynamic effects, and vehicles design features. Appl.Appl. Sci. Sci. 20212021, ,1111, ,6608 6608 3 3of of 35 34
FigureFigure 1. 1. InfluenceInfluence of of external external factors factors on on vehicle vehicle directional directional stability. stability. This paper provided a comprehensive review of the key parameters of road-vehicle This paper provided a comprehensive review of the key parameters of road-vehicle interaction that have an effect on vehicle stability and design the movement model, con- interaction that have an effect on vehicle stability and design the movement model, con- sisting of a finite number of parameters, for investigation of stability. The required tire sisting of a finite number of parameters, for investigation of stability. The required tire parameters were determined and revised using the theoretical and experimental methods. parameters were determined and revised using the theoretical and experimental methods. Following the numerical investigation using the dynamic models, experimental studies Following the numerical investigation using the dynamic models, experimental studies were performed to verify the model. The effect of the road parameters on the driving were performed to verify the model. The effect of the road parameters on the driving characteristics of the vehicle was investigated. A revised dynamic model of the vehicle characteristics of the vehicle was investigated. A revised dynamic model of the vehicle directional stability was developed assessing the effect of vehicle suspensions on vehicle directional stability was developed assessing the effect of vehicle suspensions on vehicle stability and traffic safety. The generated results have enabled a more accurate assessment stability and traffic safety. The generated results have enabled a more accurate assessment of the influence of road unevenness and vehicle tire and car body vibrations on vehicle of the influence of road unevenness and vehicle tire and car body vibrations on vehicle stability, and the development of the recommendations on safe movement on the roads stability,using the and available the development characteristics. of the recommendations on safe movement on the roads usingBased the available on the characteristics. topics discussed above, the main contributions of this paper are as follows:Based (1) on wethe examinetopics discussed the vehicle above, stability the main dynamic contributions models of and this develop paper are a refinedas fol- lows:methodology (1) we examine for estimating the vehicle the stability key parameters dynamic of models these models; and develop (2) we a produce refined method- a revised ologytire smoothing for estimating function the key model parameters that takes of into these account models; tread (2) we deflections; produce (3)a revised we investi- tire smoothinggate vertical function dynamics model models that takes and into develop account recommendations tread deflections; for (3)their we investigate use in-vehicle ver- ticalstability dynamics studies; models (4) weand verifydevelop the recommendations validity of the model for their of use the in-vehicle vertical dynamicsstability studies; using (4)experimental we verify the studies; validity and of (5)the we model establish of the the vertical vehicle dynamics stability using model, experimental allowing evaluation studies; andof vehicles’ (5) we establish technical the performancevehicle stability and model, safe drivingallowing conditions evaluation of of vehiclesvehicles’ ontechnical roads withper- formancedifferent and characteristics. safe driving conditions of vehicles on roads with different characteristics.
2.2. Investigation Investigation the the Tire-Road Tire-Road Interaction Interaction 2.1.2.1. Investigation Investigation of of the the Tire Tire Properties Properties WhenWhen dealing dealing with with the the vehicle vehicle stability stability tasks, tasks, the the issue issue of of tire-road tire-road interaction interaction is is one one ofof the the most most complex complex tasks. tasks. Analysis Analysis of of the the vehicle vehicle wheel-to-road wheel-to-road contact contact requires requires taking taking intointo account account the the deformation deformation properties properties of of the the tire. tire. The The complexity complexity of of the the task task stems stems from from thethe uncertainty uncertainty of of the the operating operating conditions conditions of of the the tire tire and and the the complex complex structure structure of of the the tire.tire. The The complexity complexity of of the the problem problem is is defined defined by by the the mere mere fact fact that that no no consensus consensus has has been been reachedreached on on the the shape shape of of the the tire. tire. TheThe performed performed studies studies of of the the tire-road tire-road interaction interaction [36] [36] confirmed confirmed that that it it would would be be necessarynecessary to to consider consider the the deformation deformation properti propertieses of of the the band band in in the the investigation investigation of of the the compensatingcompensating function function of of the the band. band. For For this this purpose, purpose, an an additional additional element element is is added added to to thethe quarter-car quarter-car model model (Figure (Figure 2).2). Additional Additional st studiesudies are are needed in order to identify the deformation properties of the element. deformation properties of the element. The methodology may be used where deformation of an additional element, i.e., the tread, is needed in the quarter-car model (Figure 2). Whereas both tread layers (with and without the pattern) are not made of composite material, the equations applicable to rub- ber parts are used here as they consider that ν = 0.5 [37]. The tread model is comprised of Appl. Sci. 2021, 11, 6608 4 of 35
Appl. Sci. 2021, 11, 6608 4 of 34 two consistently joined layers: the one with a pattern and the one without a pattern. The deformation of the part with a pattern was assessed using three models.
Figure 2. Quarter-car model: Lk—tire-road contact length; m1, c1 and k1—mass, stiffness, and damp- Figure 2. Quarter-car model: Lk—tire-road contact length; m1, c1 and k1—mass, stiffness, and ing ratio of the tire tread part; c2, k2—tire stiffness and damping ratio; c3, k3—suspension stiffness damping ratio of the tire tread part; c2, k2—tire stiffness and damping ratio; c3, k3—suspension and damping ratio; m2—unsprung mass; m3—sprung mass. stiffness and damping ratio; m2—unsprung mass; m3—sprung mass.
TheThis methodology paper analyzes may radial be used tires where of three deformation types that of differ an additional considerably element, in terms i.e., the of tread,structure: is needed passenger in the car quarter-car (175/70R13; model 195/50R15), (Figure2). Whereaslight-duty both truck tread (185/75R14C), layers (with andand withoutheavy-duty the pattern)truck (12.00R20) are not made tires. of In composite the assessment material, of the the vehicle equations tire-road applicable interaction, to rubber it partsis first are necessary used here to describe as they considerthe deformation that ν = pr 0.5operties [37]. Theof the tread tire modelunder the is comprised vertical load. of twoThe consistentlystatic load consists joined layers:of the vehicle the one and with cargo a pattern mass. and The the loaded one withouttire is subject a pattern. to defor- The deformationmation when of forming the part the with contact a pattern area was with assessed the road using pavement. three models. The tire-road contact is depictedThis in paper Figure analyzes 3, and tire radial element tires ofcharac threeteristics types thatare presented differ considerably in Table 1. inThe terms follow- of structure:ing assumptions passenger were car used (175/70R13; when determining 195/50R15), the area light-duty of the contact truck (185/75R14C),with the road: and heavy-duty• the contact truck with (12.00R20) the road tires. has an In elliptical the assessment shape and of the is vehiclelimited tire-roadby the tire interaction, width and it is firstlength necessary of the contact to describe with the the road; deformation properties of the tire under the vertical load.• the The layers static under load the consists tread of are the sufficiently vehicle and stiff, cargo and their mass. length The loaded does not tire change is subject dur- to deformationing deformation when of forming the tire the(ld ≈ contact l0) (Figure area 3a); with the road pavement. The tire-road contact• there is depicted is the linear in Figure dependency3, and tire between element the characteristics radial load are of presentedthe tire and in Tableradial1 .defor- The followingmation assumptions Δh; were used when determining the area of the contact with the road: •• thethe contactpressure with is distributed the road has evenly an elliptical across the shape contact and surface. is limited by the tire width and lengthThe second of the and contact third withassumptions the road; enable determining the length of the tire-road contact: • the layers under the tread are sufficiently stiff, and their length does not change during deformation of the tire (l ≈ l ) (Figure=⋅3a); − d LrkRz0/2 arccos((1 CRr )/ ) (1) • there is the linear dependency between the radial load of the tire and radial deforma- wheretion CR∆—hradial; stiffness of the tire and r—free radius of the tire. • theMore pressure accurate is distributedinformation evenly on the acrosstire element the contact loads surface.is provided by the model ‘elastic ring withThe second the outer and elastic third layer’ assumptions (Figure enable4a,b) [38] determining. The ring is the considered length of to the be tire-road the elastic contact band;: the tensile force N0 caused within it by the inner pressure does not change during deformation; and the ring points are displacedLk/2 only= r in· arccos the radial((1 − directionCRRz)/ byr) the measure w. (1)
where CR—radial stiffness of the tire and r—free radius of the tire. More accurate information on the tire element loads is provided by the model ‘elastic ring with the outer elastic layer’ (Figure4a,b) [ 38]. The ring is considered to be the elastic band; the tensile force N0 caused within it by the inner pressure does not change during deformation; and the ring points are displaced only in the radial direction by the measure w. Appl. Sci. 2021, 11, 6608 5 of 34 Appl. Sci. 2021, 11, 6608 5 of 35
(a) (b)
FigureFigure 3. 3.TireTire load load scheme: scheme: (a ()a tire) tire deformation deformation not not assessed; assessed; (b (b) tire) tire deformation deformation assessed. assessed.
TableTable 1. 1. DataData of of the the tire tire components. components.
TireTire Element Element TireTire Thickness,Thickness, mm mm Fiber Fiber Angle Angle δ1a = 8.0 – 175/70R13 δ1a = 8.0 – 175/70R13 δ1δb1 =b =7.0 7.0 – – ◦ A Breaker, 3 layers 2δ2 = 4.1 2θ2 = 20 A δ = 4.1 θ = 20°◦ Breaker,Cord, 1 layer3 layers δ3 = 1.2 θ3 = 90 δ3 = 1.2 θ3 = 90° Cord, 1 layer δ4 = 1.0 – δ4 = 1.0 – δ1a = 8.1 – δ1a = 8.1 – 195/50R15195/50R15 δ1b = 7.8 – 1b ◦ B Breaker, 2 layers δ δ 2== 7.8 3.9 θ2–= 20 ◦ B Cord, 2 layers δ2δ = =3.9 1.7 θ2θ = =20° 90 Breaker, 2 layers 3 3 δ3δ =4 =1.7 1.3 θ3 = 90°– Cord, 2 layers δδ41 a= =1.3 12.0 – – 185/75R14C δ1b = 11.0 – δ1a = 12.0 – ◦ C Breaker,185/75R14C 2 layers δ2 = 3.2 θ2 = 20 δ1b = 11.0 – ◦ Cord, 2 layers δ3 = 1.8 θ3 = 90 C δ2 = 3.2 θ2 = 20° Breaker, 2 layers δ4 = 1.0 – δ3 = 1.8 θ3 = 90° Cord, 2 layers δ1a = 17.0 4 12.00R20 δδ1 b= =1.0 15.0 – D Breaker, 4 layers δ = 6.4 ◦ δ1a =2 17.0 θ2 = 20 Cord,12.00R20 6 layers δ = 7.2 ◦ δ1b =3 15.0 θ3 = 90
δ4 = 1.0 – D δ2 = 6.4 θ2 = 20° Note: θ—cord angle in relationBreaker, to the rolling 4 layers direction. The indexes of the measures have the following meaning: δ3 = 7.2 θ3 = 90° 1a—tread part with a pattern, 1bCord,—continuous 6 layers tread part, 2—breaker, 3—cord, 4—inner sealing liner. δ4 = 1.0 – Note: θ—cord angle in relation to the rolling direction. The indexes of the measures have the fol- lowing meaning: 1a—tread part with a pattern, 1b—continuous tread part, 2—breaker, 3—cord, 4—inner sealing liner. Appl. Sci. 2021, 11, 6608 6 of 34 Appl. Sci. 2021, 11, 6608 6 of 35
(a) (b)
(c)
FigureFigure 4.4. Model ofof thethe deformabledeformable ringring ofof thethe tire:tire: ((aa)) deformabledeformable ring;ring; ((bb)) deformabledeformable ringring element;element; ((cc)) cross-sectioncross-section schemescheme for assessing the impact of the cord; 1—elastic layer (tread); 2—breaker layers; 3—elastic base (side walls); 4—disk. for assessing the impact of the cord; 1—elastic layer (tread); 2—breaker layers; 3—elastic base (side walls); 4—disk.
WhereWhere thethe impactimpact ofof thethe cordcord threadsthreads isis assessed,assessed, thethe radialradial deformationdeformation modelmodel isis used,used, and the tire reaction is assessed using using single single pressure pressure q z ==C z w .. The The sidewallssidewalls reactreact onlyonly to the displacements displacements in in the the radial radial direct direction,ion, and and thethe reaction reaction per perlength length unit unitequals equals = q s = .C Anzsw elastic. An elastic outer outerlayer layerwith the with thickness the thickness equal equalto H before to H before deformation deformation deforms deforms by ΔH byin the∆H contactin the contactarea. Inelastic area. Inelastic forces are forces not subject are not to subject assessment. to assessment. Tread reaction: Tread reaction: qCH=Δ0 (2) qz =zzsCzs∆H0 (2) Arc length x (Figure 4b) is measured from the center of contact. The main measures Arc length x (Figure4b) is measured from the center of contact. The main measures sought are [38]: sought are [38]: • model deflection hz assuming that the band does not deform in the radial direction: • model deflection hz assuming that the band does not deform in the radial direction: =+Δ() hwHz = (3) hz = (w + ∆H)\0\xx=0 (3)
• contact zone length—2x0; • contact zone length—2x0; • contact pressure distribution—qz(x); • contact pressure distribution—qz(x); • load Fz: • load Fz: x0 xZ0 Fz = 22 qxdxqz()(x)dx (4) zz (4) 00 As a result of integration and rearrangement of Equation (4), the following is obtained: Appl. Sci. 2021, 11, 6608 7 of 34
As a result of integration and rearrangement of Equation (4), the following is obtained: thαx0 2 1 thαx0 Fz = 2kzsx0 a0 1 − − a1x0 − (5) αx0 3 αx0
The unknown measures hz and x0 are determined using the following dependen- cies [38]:
x2 2 β2 1 − (thαx /αx ) = 0 + = p = − 0 0 hz 1 m , β kz/N0, m 1 2 (6) 2R βx0 α 1 + (β/a)thαx0
The described methodology validates the specifics of the tire ring load, namely, that the tire pressure in the contact area is not equal to the internal pressure of the tire due to the breaker and the cord layers. In order to refine the model, the deformation characteristics of the ring consisting of the cord and breaker layers and the tread as a separate deformable element are determined.
2.2. Deformation Properties of the Ring In order to determine the deformation characteristics of the tire, the tire is divided into elements in the band consisting of the tread, breaker, cord, and protective layers attached by a flexible joint to the sides by the cord-reinforced rubber. The deformation characteristics of the individual tire elements were determined by breaking them down into monolayers made of rubber or rubber reinforced with unidirectional fiber. The elastic characteristics of the monolayer were determined according to the composite mechanic’s equations [39–42]:
!−1 Em EC1 = Ef φf + Em(1 − φf ), EC2 = Em 1 − φf + φf (7) Ef 2
!−1 Gm GC12 = Gm 1 − φf + φf , νC12 = ϕ f νf 12 + ϕmνm (8) Gf 12
where Ec and Gc—modulus of elasticity and shear of the composite; Em and Gm—modulus of elasticity and shear of the matrix; Ef and Gf—the modulus of elasticity and shear of the filler; ϕf and ϕm—volumetric parts of the filler and matrix; νC, νf, and νm—Poisson’s coefficients of the composite, filler and matrix; indices 1, 2—the direction along the fiber and the direction transverse to the fiber, respectively. The properties of the cord and breaker layers in the longitudinal and transverse directions were determined using expressions of the plane stress state for the anisotropic monolayer [39–42]: 1 − νyx ηxy εx Ex Ey Gxy σx µ = − v 1 xy εy Ex Ey Gxy σy (9) η µy γxy x 1 τxy Ex Ey Gxy
The components of the elasticity matrix were determined according to [39–42] using the axes provided in Figure5: Appl. Sci. 2021, 11, 6608 8 of 35
−1 cs44 1 ν Ecs()θ =++22 −2 21 , x EE12 G 122 E −1 sc44 1 ν Ecs()θ =++22 −2 21 , y EE12 G 122 E
−1 222 11 ν ()cs− Gcs()θ =+++4222 21 , xy EE12 E 2 G 12 Appl. Sci. 2021, 11, 6608 8 of 34 ν ν yx 44 2211 1 ()θ =−−+−21 ()cs cs, (10) EEy 21212 EEG
η 4 224 c scs 2 2 1 ν21ν −1 1 Exyx(θ) = ( + + c s ( 22− 2 ))21 , ()θ =+−−−2csE1 E2 () cG12 s E2 , 4 4 GEEEGs c 2 2 1 ν21 −1 2 Exyy(θ) = ( + 21+ c s ( − 2 )) 2, 12 E1 E2 G12 E2 2 2 2 2 2 1 1 ν21 (c −s ) −1 Gμxy(θ) = (4c s ( 22+ + 2 ) + ν ) , xy csE1 E2 E2 G12 1 ν θ =++−−22 21 yx () ν221cs4 4 2() c2 1 s 1 1 , E (θ) = E (c − s ) − c s ( E + E − G ), (10) GEEEGyxy 2 211 2 2122 12 η 2 2 xy (θ) = 2cs{ s + c − (c2 − s2)( ν21 − 1 )}, Gxy E2 E1 E2 2G12 µ 2 η 2 E μ E xy (θ) = 2cs{ηc =+ sxyx+ (c2μ− s=2)(xνyy21 − 1 )}, Gxy Ex2 E1 , y E2 2G12 η GE µ EG η = xy xy, µ = xy y xy x Gxy y Gxy wherewherec =c = cos cosθ θand ands = s sin= sinθ; θx-direction; x-direction forming forming an anglean angleθ with θ with the reinforcementthe reinforcement direction direc- (Figuretion (Figure5). 5).
FigureFigure 5. 5.Layered Layered composite composite element: element: 1-the 1-the direction direction along along the the fiber; fiber; 2-the 2-the direction direction transverse transverse to to thethe fiber. fiber.
TheThe modulus modulus of of elasticity elasticity of theof componentsthe componen andts theand Poisson’s the Poisson’s coefficients coefficients were taken were fromtaken the from works the [works36,43]. [36,43]. The methodology The methodology used forused this for purpose this purpose is presented is presented in [39 in– 42[39–]. Direct42]. Direct modulus modulus of elasticity of elasticity calculations calculations does does not not provide provide an an answer answer in in relation relation to to the the characteristicscharacteristics of of the the tread tread and and unreinforced unreinforced ring ring section section that that are are required required for for the the model model providedprovided in in Figure Figure4. 4. To To determine determine the the stiffness stiffness of of the the tread treadC zC,z tread, tread deformations deformations are are analyzedanalyzed additionally. additionally. The The obtained obtained data data are are presented presented in in Table Table2. 2.
2.3. InvestigationTable 2. Modulus of Tread of Deformationelasticity of the tire elements, MPa. Tread deformations were investigated using several models:Tire porous composite and rubberCharacteristics prism models. The research studies analyzing the pneumatic tires specify fairly A B C D wide limits of the modulus of elasticity of rubber Er = 2 ... 20 MPa. In the study [36], a El 1667·103 1667·103 ––– 1367·103 simplified model of a passenger car tire was used, Er = 18 MPa, in order to obtain adequate Breakerdeformation (with nylon) characteristics of the wholeE tire.t In26,736 order to specify26.736 the initial ––– characteristics, 24,568 an experiment was performed to determineEr the26,736 313 mm wide26.736 section of––– the tread 24,568 band of the truck tire corresponding to 12.00R20 tire treads. The height of the tread part with a pattern was δ1a = 17.8 mm, without pattern δ1b = 3.9 mm. At compression, the tread was under the action of the load corresponding to the nominal radial load Rznom of the 12.00R20 type tires in terms of the pressure. An experimental setup for monotonous tension compression tests consisted of a 50 kN testing machine Instron (Norwood, MA, USA) 8801 series and an electronic device that was designed to record the stress-strain curves. The mechanical characteristics were measured with an error not exceeding ±1% of the deformation scale. The data collection was performed by using the Bluehill Universal Software. Compression tests were performed in accordance with ASTM D575-91: Standard Test Methods for Rubber Properties in Compression. The obtained load-deformation dependence (Figure6) is very close to the linear (correlation coefficient rxy = 0.947), proving the applicability of the description of the linear elastic material to the study of tread deformations. The obtained modulus of tread elasticity Appl. Sci. 2021, 11, 6608 9 of 34
determined on the basis of the experimental results was 7.84 MPa, taking into consideration the tread part with a pattern according to the filling factor, Er = 6.16 MPa, which was much lower than in the study [36].
Table 2. Modulus of elasticity of the tire elements, MPa.
Tire Characteristics ABCD 3 3 3 El 1667·10 1667·10 —— 1367·10
Breaker (with nylon) Et 26,736 26.736 —— 24,568
Er 26,736 26.736 —— 24,568
El 73,119 73,119 73,119 73,119
Breaker (with steel) Et 21,737 21,737 21,737 21,737
Er 23,684 23,684 23,684 23,684
El 24,625 24,625 26,819 —— 3 3 3 Cord (with viscose) Et 3553·10 3553·10 4339·10 ——
Er 24,625 24,625 26,819 —— Inner sealing layer E 18 18 18 18 Part of the ring up to the cord E 16,601 16,630 17,118 16,511 Whole ring E 16,890 16,906 17,311 16,541 Sidewall elasticity modulus during radial tension E, 507.5 531.15 788.1 551.7 Equivalent sidewall elasticity modulus during radial bending E flat band, 23.32 24.37 29.77 67.89 ring. 42.69 43.04 61.51 113.4
Equivalent modulus of elasticity of the band bending Ec flat band in the circular direction, 15.89 17.74 19.63 15.92 perpendicular to the circular direction ring, 15.35 17.37 16.34 15.90 in the circular direction. 23.73 24.21 22.50 18.50 Relative position of the neutral layer e/ ∑ δ at bending flat band in the circular direction, 0.288 0.281 0.383 0.377 perpendicular to the circular direction ring, 0.108 0.11 0.088 0.124 in the circular direction. 0.286 0.277 0.376 0.372 Appl. Sci. 2021, 11, 6608 10 of 35 Note: E—in all directions, El—the rolling direction of the tire, Et—perpendicular to the rolling direction of the tire, and Er—in the radial direction (in relation to the tire).
FigureFigure 6.6. SummarySummary experimentalexperimental data:data: 1—experimental1—experimental data;data; 2—generalized2—generalized line.line.
Whereas the tread consists of the parts with and without a pattern and the defor- mation properties of the tread are evaluated in the model as a single generalized layer, the possibilities of assessment of the tread part with a pattern were examined. The simplest model assesses the filling of the part with a pattern (Figure 7), which can be executed using the expression of porous composites [44]: = EkEug (11)
where ku—tread filling factor, where ku = Ai/Ap, Ai—the area of the continuous pattern in the part of the tread part selected according to the repetition of the pattern, Ap—area of the selected part of the tread; and Eg—tread rubber modulus of elasticity.
(a) (b) (c)
Figure 7. Tire tread patterns: (a) for a passenger car; (b) for a light-duty truck; and (c) for a truck.
Since the tread pattern is different across differed parts of the tread (Figure 7), the filling factor for each tread band was calculated separately. In this case, the equivalent modulus of elasticity of the tread was equal to: Σkk EE= ui bi 1 g −−Σ() (12) 11kkkauibi
where Eg—tread rubber modulus of elasticity; ka = δ1a/δ1; δ1 = δ1a + δ1b; kbi—share of the width of the i-th band in the width of the tread, δ1—tire tread thickness, δ1a—tread part with a pattern, and δ1b—continuous part of the tread. The results obtained for different tires are presented in Table 3. The modulus of elas- ticity of the whole tread was calculated using the equation in [44]: Appl. Sci. 2021, 11, 6608 10 of 35
Appl. Sci. 2021, 11, 6608 10 of 34
Figure 6. Summary experimental data: 1—experimental data; 2—generalized line. Whereas the tread consists of the parts with and without a pattern and the deforma- tionWhereas properties the of tread the tread consists are evaluatedof the parts in thewith model and without as a single a pattern generalized and the layer, defor- the mationpossibilities properties of assessment of the tread of are the evaluated tread part in with the amodel pattern as werea single examined. generalized The layer, simplest the possibilitiesmodel assesses of assessment the filling of of the the part tread with part a pattern with a (Figurepattern7 ),were which examined. can be executed The simplest using modelthe expression assesses the of porous filling of composites the part with [44]: a pattern (Figure 7), which can be executed using the expression of porous composites [44]: E ==kuEg (11) EkEug (11)
where ku—tread filling factor, where ku = Ai/Ap, Ai—the area of the continuous pattern in where ku—tread filling factor, where ku = Ai/Ap, Ai—the area of the continuous pattern in the part of the tread part selected according to the repetition of the pattern, Ap—area of the the part of the tread part selected according to the repetition of the pattern, Ap—area of selected part of the tread; and Eg—tread rubber modulus of elasticity. the selected part of the tread; and Eg—tread rubber modulus of elasticity.
(a) (b) (c)
FigureFigure 7. 7. TireTire tread tread patterns: patterns: ( (aa)) for for a a passenger passenger car; car; ( (bb)) for for a a light-duty light-duty truck; truck; and and ( (cc)) for for a a truck. truck.
SinceSince the the tread tread pattern pattern is is different different across across differed differed parts parts of of the the tread tread (Figure (Figure 77),), thethe fillingfilling factor factor for for each each tread tread band band was calculated separately. separately. In In this this case, case, the the equivalent modulusmodulus of of elasticity elasticity of of the the tread tread was was equal equal to: to: Σ kkkuik bi EE= Σ ui bi E1 =1 Eg g −−Σ() (12)(12) 1 −11kakkk(auibi1 − Σkuikbi)
wherewhere EEgg——treadtread rubber rubber modulus modulus of of elasticity; elasticity; ka k=a δ=1a/δ11a; /δδ11 =; δ11a =+ δ1ba; k+bi—δ1bshare; kbi—share of the width of the ofwidth the i-th ofthe bandi-th in band the width in the of width the tread, of the δ tread,1—tireδ 1tread—tire thickness, tread thickness, δ1a—treadδ1a—tread part with part a pattern,with a pattern, and δ1b— andcontinuousδ1b—continuous part of the part tread. of the tread. TheThe results results obtained obtained for for different different tires tires are are presented presented in Table in Table 3. The3. Themodulus modulus of elas- of ticityelasticity of the of whole the whole tread tread was wascalculated calculated using using the equation the equation in [44]: in [44]:
1 δ δ 1 = 1a + 1b · (13) E1 E1a E1b δ1
Table 3. Tread modulus of elasticity obtained by assessment using the first model, MPa.
Tire Type Part with a Pattern E1a, Continuous Part E1b, Tread E1, A 5.72 7.84 6.16 B 5.88 7.84 6.22 C 7.08 7.84 6.51 D 5.33 7.84 5.37
This model does not accurately reflect the operating conditions of the tread. Defor- mation of the tread part with a pattern may be restricted by the contact with the road or a continuous breaker. Moreover, the deformation of the tread part with a pattern does not fully correspond to the ideal compression conditions, as the deformation of the ends at Appl. Sci. 2021, 11, 6608 11 of 34
compression leads to a change in the element shape. The course of deformation would be more accurately estimated using the model of a prism subject for compression by applying a rubber-shock-absorber calculation methodology. The deformations of the tread part with a pattern can be assessed more precisely using the calculation methods applicable to the rubber shock absorbers. Under this approach, the tread element is depicted as a rubber prism, deformed with more or less constrained end deformations depending on the conditions. If the deformations of the both ends of the prism-shaped rubber shock absorber are constrained, the relationship between the load and the deformation is expressed by the equation:
−1 1 1 ∆δ 4 2 2 αδ Rz = βEA − λ , λ = 1 − , β = 1 + η − η 1 − th (14) 3 λ2 δ 3 αδ 2
where A—prism base area; δ—prism height; ∆δ—prism deformation; and α and η—shape change parameters. Whereas α and η are interrelated, a system of equations derived from the energy minimum condition is proposed in the study [45]:
2 2 = 1 48(1+η −η) α b2 η2( a )+(1−η)2 , b (15) αδ−2th αδ + 2− η( a )2+η−1 2 = 1 + 1 η η b 2 αδ αδ 2 1−2η η2( a )2+(1−η)2 αδ(th 2 −1)+2th 2 b where a, b—dimensions of the prism base. Equations (14) and (15) were solved by approximation. This model of the tread part with a pattern was applied to two cases: with the end deformations not constrained and with the end deformations completely constrained. If the deformations of the shock absorber ends were not constrained, β = 1. In case the deformations were constrained, the stiffness of the individual bands of the tread part with a pattern would increase by up to 4.5 times (Table4).
Table 4. The values of the stiffening coefficient of the tread elements for individual tread bands.
Tire Type Band Number Length a, mm Width b, mm β Band 1 39 29 3.83 A 2 32 29.6 3.67 1 37 29 3.75 B 2 31 29 3.67 3 239 24 4.58 1 33 20 2.62 C 2 36 18 2.08 1 116 42 3.33 D 2 90 41 3.17
Under actual conditions, the constraints of the end deformations of the tread part with a pattern were not absolute. In relation to the contact with the road, they are determined by the tire-road friction, while the other end is in contact with the tread part without a pattern. Deformations of the tread part without a pattern are restricted by the breaker layers that have the anisotropic properties (Table2), and the modulus of elasticity of the layer is only 1.2–1.5 times higher than the modulus of elasticity of rubber in the direction transverse to the reinforcement. Under these conditions, the deformation conditions of the tread part with a pattern are closer to those of a shock absorber with freely deforming ends. This option was used for further calculations. Appl. Sci. 2021, 11, 6608 12 of 34
Both models offer different characteristics: the modulus of elasticity E1 or deformation ∆δ corresponding to the radial load Rz value. For comparison of the results, the generalized characteristic describing the relationship between load and deformation—stiffness C1 was used: ∆δ = C1 · Rz (16)
where C1—tread stiffness and Rz—radial load of the tire. When determining the measure C1 according to the first or second models, it should be taken into the account that the road contact area and the volume of the deformable tread part depend on the load. The following assumptions were used when determining the area of the contact with the road: • the contact with the road has an elliptical shape and is limited by the tire width and length of the contact with the road; • the layers under the tread are sufficiently stiff, and their length does not change during deformation of the tire; • there is the linear dependency between the radial deformation and radial load; • the pressure is distributed evenly across the contact surface. The provided assumptions enable determination of the pressure in the contact area. Nonlinear q(Rz) dependence was obtained for all the tested tires (Figure8). The tire load is calculated using the nominal load of the tire Rznom, which is specified in the catalogs, by setting the limits Rzmin = 0.2Rznom and Rzmax = 1.4Rznom. The both models were compared to the first model by calculating the following:
Appl. Sci. 2021, 11, 6608 = q · E 13 of 35 ∆δ δ/ 1 (17)
FigureFigure 8. 8.Radial Radial loadsloads ofof thethe tirestires (marking (marking according according to to Table Table2). 2).
Nonlinear dependence was observed C1a(Rz) and was related to nonlinear dependence between Rz and the length of the contact with the road. This generated the nonlinear q(Rz) dependence that reduced tread deformation in the area of higher loads. The specifics of the tire-road contact are explained by the fact that the type of tire C1a(Rz) has no significant effect on the nature of the dependence (Figure9). The influence of tread wear on its stiffness was determined using the second model that provides more accurate reflection of the deformation conditions of short rubber prisms [36]. The effect of prism end constraint was not assessed when calculating tread stiffness. For the calculation of the tread stiffness of type A tire, the tread height was 8.1 mm for a new tire and 1.6 mm for a worn tire, and an intermediate wear level (4.9 mm) was provided. The tread wear changed the C-Rz dependence slightly (Figure 10) but had a significant effect on the stiffness value. The effect was not monotonic as the stiffness would be subject
to particular change at the beginning of the wear. For the tested tires, tread wear could increaseFigure 9. theDependence stiffness of relative the tread stiffness in relation of the totread the pa nominalrt with a load pattern by on up the to 1.3relative times. radial load for different tires (compressed prism model; marking according to Table 2).
The influence of tread wear on its stiffness was determined using the second model that provides more accurate reflection of the deformation conditions of short rubber prisms [36]. The effect of prism end constraint was not assessed when calculating tread stiffness. For the calculation of the tread stiffness of type A tire, the tread height was 8.1 mm for a new tire and 1.6 mm for a worn tire, and an intermediate wear level (4.9 mm) was provided. The tread wear changed the C-Rz dependence slightly (Figure 10) but had a significant effect on the stiffness value. The effect was not monotonic as the stiffness would be subject to particular change at the beginning of the wear. For the tested tires, tread wear could increase the stiffness of the tread in relation to the nominal load by up to 1.3 times. Depending on the constraint conditions, the second model provided very different results. For more accurate assessment, the effect of the latter constraint, given that the properties of the ring consisting of the cord and breaker are anisotropic, and the influence of deformation characteristics of the tire ring consisting of the tread, breaker, cord and protective layer, on the tread properties were determined, and the road-to-tread contact model was adjusted. Appl. Sci. 2021, 11, 6608 13 of 35
Appl. Sci. 2021, 11, 6608 13 of 34
Figure 8. Radial loads of the tires (marking according to Table 2).
Appl. Sci. 2021, 11, 6608 FigureFigure 9.9. DependenceDependence ofof relativerelative stiffnessstiffness ofof thethe treadtread partpartwith witha apattern patternon onthe therelative relative radial radial14 load ofload 35 for different tires (compressed prism model; marking according to Table 2). for different tires (compressed prism model; marking according to Table2). The influence of tread wear on its stiffness was determined using the second model that provides more accurate reflection of the deformation conditions of short rubber prisms [36]. The effect of prism end constraint was not assessed when calculating tread stiffness. For the calculation of the tread stiffness of type A tire, the tread height was 8.1 mm for a new tire and 1.6 mm for a worn tire, and an intermediate wear level (4.9 mm) was provided. The tread wear changed the C-Rz dependence slightly (Figure 10) but had a significant effect on the stiffness value. The effect was not monotonic as the stiffness would be subject to particular change at the beginning of the wear. For the tested tires, tread wear could increase the stiffness of the tread in relation to the nominal load by up to 1.3 times. Depending on the constraint conditions, the second model provided very different results. For more accurate assessment, the effect of the latter constraint, given that the properties of the ring consisting of the cord and breaker are anisotropic, and the influence FigureFigureof deformation 10.10.Dependence Dependence characteristics of of stiffness stiffness ofof ofB the tireB tire tire tread tread ring part part consisting with with a pattern a pattern of onthe theon tread, treadthe tread breaker, height: height: 1—new cord 1—new tire;and tire; 2—worn by 50%; and 3—worn out. 2—wornprotective by 50%;layer, and on 3—worn the tread out. properties were determined, and the road-to-tread contact model was adjusted. DependingThe tape for on Nokian the constraint passenger conditions, car tires the175/70R13 second 82T model was provided modeled. very Its differentmore de- results.tailed structure For more is accurateprovided assessment,in Table 5. the effect of the latter constraint, given that the properties of the ring consisting of the cord and breaker are anisotropic, and the influence ofTable deformation 5. Layer distribution characteristics of 175/70R13 of the 82T tire tire ring bands. consisting of the tread, breaker, cord and protectiveLayer layer, No. on the tread properties Layer were Designat determined,ion and theThickness, road-to-tread mm contact model was adjusted.1 Tread with a pattern 2.9 The tape for Nokian passenger car tires 175/70R13 82T was modeled. Its more detailed 2 Tread without a pattern 2.1 structure is provided in Table5. 3 Nylon-rubber 0.8 Table 5. Layer distribution4 of 175/70R13 82TSteel-rubber tire bands. 0.9 5 Rubber 0.9 Layer6 No. LayerSteel-rubber Designation Thickness,0.9 mm 17 TreadViscose-rubber with a pattern 2.90.9 28 TreadSealing without layer a pattern 2.10.9 3 Nylon-rubber 0.8 The properties of the cord and breaker layers in the longitudinal and transverse di- 4 Steel-rubber 0.9 rections were determined using expressions of the plane stress state for the anisotropic monolayer [395–42] by applying them to the Rubber axes presented in Figure 5: 0.9 6 Steel-rubber 0.9 ν 7 Viscose-rubber1 21 0.9 − 0 EE 8εσ Sealing12 layer 0.9 xx v 1 εσ=−, 21 yyTT 0 EE γ 22 τ xy1 xy 00 G12
cs22 cs cs22 −2 sc ,22=−=−22 Tsccs, Tsc2 sc (18) −−22 −−22 22(cs cs c s ) sc sc() c s
where c = cosθ and s = sinθ. Deformation properties of the layers were determined according to the methodology set out in Section 2.2. An assumption was made that nylon and viscose fibers were aniso- tropic, with a different modulus of longitudinal elasticity and transverse to the fiber. The issue of identification of the layers emerged in the investigation of the structure of the modeled tire band. There is a rubber layer between the individual bands of the breaker reinforced with steel wire. As a result, the band was investigated additionally by Appl. Sci. 2021, 11, 6608 14 of 34
The properties of the cord and breaker layers in the longitudinal and transverse directions were determined using expressions of the plane stress state for the anisotropic monolayer [39–42] by applying them to the axes presented in Figure5:
1 − ν21 0 εx E1 E2 σx { ε } = [T0]{ − v21 1 0 }[T]{ σ } y E2 E2 y 1 γxy 0 0 τxy G12 (18) c2 s2 cs c2 s2 −2sc [T0] = [ s2 c2 −cs ], [T] = [ s2 c2 −2sc ] −2cs 2cs (c2 − s2) sc −sc (c2 − s2)
where c = cosθ and s = sinθ. Deformation properties of the layers were determined according to the methodology set out in Section 2.2. An assumption was made that nylon and viscose fibers were Appl. Sci. 2021, 11, 6608 anisotropic, with a different modulus of longitudinal elasticity and transverse to the15 fiber.of 35
The issue of identification of the layers emerged in the investigation of the structure of the modeled tire band. There is a rubber layer between the individual bands of the identifyingbreaker reinforced the rubber with layer steel between wire. As the a result,reinforc theed band layers was (Figure investigated 11a; Table additionally 5) and using by simplifiedidentifying separation the rubber into layer two between layers (Figure the reinforced 11b). layers (Figure 11a; Table5) and using simplified separation into two layers (Figure 11b).
(a) (b)
Figure 11. Breaker layer calculation schemes: ( a)) by by identifying identifying the the interm intermediateediate rubber layer and ( b) without identifying the intermediate rubber layer. the intermediate rubber layer.
InIn the the investigation investigation of of the the deformation deformation prop propertieserties of the whole band, the model used inin layered layered composites was was applied initially. The The calculated modulus modulus of of elasticity of the bandband as as a a composite composite was was calculated calculated by by taking taking into into consideration consideration the modulus of elasticity ofof the the layers joined in parallel ( X and Y directions) or consistently ( Z direction). SinceSince the the reinforcement reinforcement angles angles of ofthe the adjacent adjacent breaker breaker layers layers differed differed by the by sign the only sign (+20°only (+20and ◦−20°),and −it 20was◦), assumed it was assumed that the thatresulting the resulting shear forces shear in forcesthe layers in the were layers counter- were balancedcounterbalanced and were and not were transmitted not transmitted to other tolayers. other layers. ConsideringConsidering the the large large difference difference between the modulus of elasticity of the reinforcing fiberfiber and the rubber, it was determined that the modulus of elasticity of the layer in the directiondirection of reinforcementreinforcement waswas not not significantly significantly affected affected by by the the modulus modulus of elasticityof elasticity of the of therubber. rubber. The The method method of identification of identification of theof the layers layers (a and(a and b) alsob) also did did not not have have any any effect effect in inthis this case. case. At At the the same same time, time, it wasit was found found that that the the model model of elasticityof elasticity of theof the layer layer would would be bemore more dependent dependent on on the the modulus modulus of of elasticity elasticity of of the the rubber rubber in in the the breaker breaker fiber fiber direction directionX Xand andZ .Z The. The modulus modulus of of elasticity elasticity in inX Xdirection direction was was less less dependent dependent on on the the total total modulus modulus of ofelasticity elasticityE. ByE. By changing changing the modulusthe modulus of elasticity of elasticity of the of rubber the rubber in the in range the ofrange 2 ... of20 2…20 MPa, MPa,Ex differed Ex differed by 3.7% by 3.7% for model for model (a) and (a) 3.5% and for3.5% model for model (b). E (b).z dependence Ez dependence was particularlywas partic- ularlyevident: evident: rubber rubber stiffness stiffnessEr would Er would change change 10 times, 10 times, while Ewhilez would Ez would change change 9.97 times 9.97 times(Figure (Figure 12). Considering 12). Considering that the that deformation the deformation properties properties of the band of the in bandZ direction in Z direction differed differedconsiderably considerably from the from data the provided data provided by other by authors, other authors, the assumption the assumption was that was stiffer that stiffercord and cord breaker and breaker layers layers may restrict may restrict tread deformations, tread deformations, and in and the statein the of state constrained of con- strained deformation, tread deformation could be significantly influenced by the Pois- son’s ratio close to 0.5. Therefore, tread deformation under compression was investigated additionally under the assumptions that reflected the interaction between the individual layers in the identified tire band more accurately. Two cases were examined for this purpose. In the first case, it was assumed that by de- forming the band in Z direction, deformation of all layers was the same in X and Y directions: εε=== εεε=== ε xx12.... xn, yy12.... yn (19)
The sum of the resulting single axial forces in the individual layers in the X and Y axis directions equaled zero: δσ+++= δ σ δσ 11xxx 22... nxn 0, δσ+++= δ σ δσ 11yyy 22... nyn 0, (20) σσ===== σσ zz12... znzq . Appl. Sci. 2021, 11, 6608 15 of 34 Appl. Sci. 2021, 11, 6608 16 of 35
deformation, tread deformation could be significantly influenced by the Poisson’s ratio close toThe 0.5. numerical Therefore, experiment tread deformation showed underthat th compressione modulus of was elasticity investigated of rubber additionally in Z di- underrection the would assumptions increase thatgreatly reflected and result the interaction in an increase between in the the modulus individual of elasticity layers in of the the identifiedentire band tire in band Z direction more accurately. (Figure 13).
Appl. Sci. 2021, 11, 6608 16 of 35
The numerical experiment showed that the modulus of elasticity of rubber in Z di- rection would increase greatly and result in an increase in the modulus of elasticity of the Figure 12. Tire band elastic modulus in X and Z directions: 1—model a (8 layers) in X direction, 2— Figureentire 12.bandTire in band Z direction elastic modulus (Figure 13). in X and Z directions: 1—model a (8 layers) in X direction, 2—modelmodel a in a inZ direction,Z direction, 3—model 3—model b (7 b (7layers) layers) in inX direction,X direction, and and 4—model 4—model b in b inZ direction.Z direction.
Two cases were examined for this purpose. In the first case, it was assumed that by deforming the band in Z direction, deformation of all layers was the same in X and Y directions: εx1 = εx2 = ... = εxn, εy1 = εy2 = ... = εyn (19) The sum of the resulting single axial forces in the individual layers in the X and Y axis directions equaled zero:
δ1σx1 + δx2σx2 + ... + δnσxn = 0, δ1σy1 + δy2σy2 + ... + δnσyn = 0, (20) σz1 = σz2 = ... = σzn = σz = q.
The numerical experiment showed that the modulus of elasticity of rubber in Z
directionFigure 12. wouldTire band increase elastic modulus greatly andin X resultand Z directions: in an increase 1—model in the a (8 modulus layers) in of X elasticitydirection, 2— of Figure 13. Dependence of the tire band modulus of elasticity on the modulus of elasticity of rubber themodel entire a in bandZ direction, in Z direction 3—model (Figure b (7 layers) 13). in X direction, and 4—model b in Z direction. (σz = 0.2 MPa): 1—7 layers and 2—8 layers.
The numerical experiment confirmed that the models predicting uniform layer de- formation were particularly sensitive to components with the Poisson’s ratio close to 0.5. The conditional modulus of elasticity of the rubber with the layer structure corresponding to the 175/70R1382T tire, depending on the layer, was 86.3 MPa. For this reason, the model with the identified intermediate layer of rubber predicted an unrealistically high stiffness. In the second case, the tread-to-road friction was assessed, and an assumption was made that the tread deformation would be constrained as long as the shear force devel- oping in directions X and Y did not exceed the frictional force. Whereas the control exper- iment was anticipated, the numerical model assumed that the frictional forces would also act on the other side of the band, i.e., in contact with the sealing liner. The constraint pro- vided that the stresses in X and Y directions could not exceed the limit values: σσ==μσ xr yr z (21) FigureFigure 13.13. DependenceDependence ofof thethe tiretire bandband modulusmodulus ofof elasticityelasticity onon thethe modulusmodulus ofof elasticityelasticity ofofrubber rubber (σz = 0.2 MPa): 1—7 layers and 2—8 layers. (σwherez = 0.2 μ MPa):—coefficient 1–7 layers of and friction 2–8 layers. between rubber and steel. The numerical experiment confirmed that friction had an impact on the Ez measure (FigureThe 14). numerical As the coefficientexperiment of confirmed friction incr thateased, the modelsthe stiffness predicting of the uniformband in thelayer trans- de- formationverse direction were particularlywould increase sensitive rapidly. to compon The identifiedents with additional the Poisson’s rubber ratio layers close tohad 0.5. a Thesmaller conditional influence modulus in this case of elasticity (Figure of14). the rubber with the layer structure corresponding to the 175/70R1382T tire, depending on the layer, was 86.3 MPa. For this reason, the model with the identified intermediate layer of rubber predicted an unrealistically high stiffness. In the second case, the tread-to-road friction was assessed, and an assumption was made that the tread deformation would be constrained as long as the shear force devel- oping in directions X and Y did not exceed the frictional force. Whereas the control exper- iment was anticipated, the numerical model assumed that the frictional forces would also act on the other side of the band, i.e., in contact with the sealing liner. The constraint pro- vided that the stresses in X and Y directions could not exceed the limit values: σσ==μσ xr yr z (21)
where μ—coefficient of friction between rubber and steel. The numerical experiment confirmed that friction had an impact on the Ez measure (Figure 14). As the coefficient of friction increased, the stiffness of the band in the trans- verse direction would increase rapidly. The identified additional rubber layers had a smaller influence in this case (Figure 14). Appl. Sci. 2021, 11, 6608 16 of 34
The numerical experiment confirmed that the models predicting uniform layer defor- mation were particularly sensitive to components with the Poisson’s ratio close to 0.5. The conditional modulus of elasticity of the rubber with the layer structure corresponding to the 175/70R1382T tire, depending on the layer, was 86.3 MPa. For this reason, the model with the identified intermediate layer of rubber predicted an unrealistically high stiffness. In the second case, the tread-to-road friction was assessed, and an assumption was made that the tread deformation would be constrained as long as the shear force developing in directions X and Y did not exceed the frictional force. Whereas the control experiment was anticipated, the numerical model assumed that the frictional forces would also act on the other side of the band, i.e., in contact with the sealing liner. The constraint provided that the stresses in X and Y directions could not exceed the limit values:
σxr = σyr = µσz (21)
where µ—coefficient of friction between rubber and steel. The numerical experiment confirmed that friction had an impact on the Ez measure (Figure 14). As the coefficient of friction increased, the stiffness of the band in the transverse Appl.Appl. Sci.Sci. 20212021,, 1111,, 66086608 direction would increase rapidly. The identified additional rubber layers had a smaller1717 ofof 3535
influence in this case (Figure 14).
FigureFigureFigure 14.14.14. DependenceDependence ofof of tiretire tire bandband band deformationdeformation deformation onon on thethe the coefficientcoefficient coefficient ofof friction offriction friction ((σσzz == ( σ 0.20.2z = MPa,MPa, 0.2 MPa, EErr == 66 MPa) using different breaker models: 1—7 layers and 2—8 layers. EMPa)r = 6 MPa)using using different different breaker breaker models: models: 1—7 1—7layers layers and 2—8 and 2—8layers. layers.
ToToTo revise reviserevise the thethe results resultsresults obtained, obtained,obtained, it it wasit waswas attempted attempattemptedted to to deformto deformdeform the thethe band bandband cut cutcut from fromfrom the thethe tire tiretire in theinin theZthedirection. ZZ direction.direction. The The experimentThe experimentexperiment was performedwaswas performeperforme usingdd usingusing the 175/70R1382T thethe 175/70R1382T175/70R1382T Nokian NokianNokian tire band tiretire usedbandband for usedused modeling. forfor modeling.modeling. The summarized TheThe summarizedsummarized experimental expeexperimentalrimental data are datadata presented areare presentedpresented in Figure inin FigureFigure15. 15.15.
Figure 15. Summary experimental data: 1—experimental data and 2—generalized line. FigureFigure 15. 15.Summary Summary experimental experimental data: data: 1—experimental 1—experimental data data and and 2—generalized 2—generalized line. line.
TheThe experimentexperiment confirmedconfirmed thatthat thethe linearlinear relationshiprelationship betweenbetween pressurepressure andand defor-defor- mationmation remainedremained upup toto thethe pressurespressures aboveabove thethe operatingoperating pressurepressure rangerange (tested(tested upup toto qq == 1.51.5 MPa).MPa). TheThe calculatedcalculated modulusmodulus ofof elasticityelasticity ofof thethe wholewhole bandband waswas 2727 MPa,MPa, whichwhich cor-cor- respondedresponded toto thethe numericalnumerical modelmodel ofof thethe simpsimplifiedlified modelmodel structurestructure (7(7 layers)layers) assessingassessing thethe frictionfriction betweenbetween thethe platesplates andand thethe band.band. InIn thethe simplifiedsimplified tiretire models,models, thethe cordcord waswas assessedassessed duringduring thethe calculationcalculation ofof radialradial loads.loads. Therefore,Therefore, thethe breakerbreaker layerslayers couldcould alsoalso bebe consideredconsidered asas conditionallyconditionally involvedinvolved inin thethe radialradial deformationdeformation ofof thethe tread.tread. Thereby,Thereby, thethe calculationcalculation expressionsexpressions ofof thethe layeredlayered compositescomposites withwith sequentiallysequentially arrangedarranged layerslayers wouldwould resultresult inin thethe following:following:
1 δδδδδδ 1 =++=++121233 **** (22) CCCC**** (22) CCCCprzprz123123 z z z z z z
where ∗∗—member of the stiffness matrix accounting for the deformation constraint. where —member of the stiffness matrix accounting for the deformation constraint. IfIf thethe stiffnessstiffness ofof thethe bandband usedused inin thethe experimentexperiment waswas modeledmodeled withwith aa layeredlayered com-com- positeposite withwith sequentiallysequentially arrangedarranged layers,layers, itsits stiffnessstiffness couldcould bebe estimaestimatedted byby EquationEquation (22).(22). ∗∗ InIn thisthis case,case, thethe membermember ofof thethe stiffnessstiffness matrixmatrix ofof thethe testedtested tiretire treadtread would would bebe equalequal toto 2727 MPa.MPa. ThisThis valuevalue waswas usedused inin thethe researchresearch forfor determinationdetermination ofof thethe CCzz measuremeasure inin EquationEquation (22).(22). Appl. Sci. 2021, 11, 6608 17 of 34
The experiment confirmed that the linear relationship between pressure and defor- mation remained up to the pressures above the operating pressure range (tested up to q = 1.5 MPa). The calculated modulus of elasticity of the whole band was 27 MPa, which corresponded to the numerical model of the simplified model structure (7 layers) assessing the friction between the plates and the band. In the simplified tire models, the cord was assessed during the calculation of radial loads. Therefore, the breaker layers could also be considered as conditionally involved in the radial deformation of the tread. Thereby, the calculation expressions of the layered composites with sequentially arranged layers would result in the following:
1 δ1 δ2 δ3 ∗ = ∗ + ∗ + ∗ (22) Cprz C1z C2z C3z Appl. Sci. 2021, 11, 6608 18 of 35 ∗ where Ci —member of the stiffness matrix accounting for the deformation constraint. If the stiffness of the band used in the experiment was modeled with a layered com- posite with sequentially arranged layers, its stiffness could be estimated by Equation (22). ∗ In thisFor case, real the tread member modeling, of the stiffness tread matrix loads of and the testedoperating tire tread conditionsEz would need be equal to be adjusted by bringingto 27 MPa. the This numerical value was model used in cl theoser research to real for conditions. determination The of third the Cz treadmeasure deformation in model Equation (22). was Forused real to tread evaluate modeling, the treadspecifics loads of and the operating tire-road conditions texture need interaction. to be adjusted Road by texture means unevennessbringing the numerical of the road model pavement. closer to real conditions.For asphalt The pavement, third tread deformation the height model of the unevenness was used1–5 tomm. evaluate The theirregularities specifics of the were tire-road assumed texture to interaction. be hemispherical Road texture protrusions means of regular spacingunevenness (Figure of the 16). road This pavement. assumption For asphalt enabled pavement, modeling the height of the of the interaction unevenness of an individual was 1–5 mm. The irregularities were assumed to be hemispherical protrusions of regular unevenness with the tread, assuming that the unevenness acted on the continuous part of spacing (Figure 16). This assumption enabled modeling of the interaction of an individual theunevenness tread. The with coefficient the tread, assuming of filling that with the unevennessspheres necessary acted on the for continuous the conversion part of the refer- enceof the pressure tread. The was coefficient equal ofto fillingπ/4 = with0.785. spheres Based necessary on the geometrical for the conversion conditions, of the it was deter- minedreference that pressure the gaps was equal between to π/4 the = 0.785. spheres Based would on the geometricalbe filled at conditions, the sphe itre was penetration into thedetermined tread equal that the to: gaps Δ = between0.455r, thewhere spheres r—radius would be of filled the sphere. at the sphere The penetrationknown spheres used for into the tread equal to: ∆ = 0.455r, where r—radius of the sphere. The known spheres used thefor theinvestigation investigation of of interaction interaction were were the the solution solution for the for elastic the elastic deformations deformations of the of the plane contact.plane contact. The Theinteraction interaction scheme scheme is shownshown in in Figure Figure 17. 17.
(a) (b) (c)
Figure 16. 16.Road Road unevenness unevenness schemes: schemes: (a)—convexities; (a)—convexities; (b)—concavities; (b)—concavities; (c)—combination (c)—combination of con- of con- vexities and and concavities. concavities.
Figure 17. Unevenness-to-tread interaction scheme.
The hemisphere-shaped unevenness with radius r pressed against the tread located on a rigid base with force F. As a result, the unevenness penetrated into the tread to depth Δ. The radius of the depression was equal to c. The elastic solution for the task was known and taken from [45]. The depression radius c was equal to:
=⋅⋅⋅⋅+3 cFrkk0.721 2 (12 ) (23)