Assessment of Tire Features for Modeling Vehicle Stability in Case of Vertical Road Excitation

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 improved 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 possible 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 avoided even on the trunk roads. In this work, we investigated the conditions for directional stability, the influence of road and vehicle parameters on the directional stability of the vehicle, and developed recommendations for the road and vehicle control systems to combine to ensure traffic safety. We have developed a refined dynamic model of vehicle stability that evaluates the influence of tire tread and suspensions. The obtained results allow a more accurate assessment of the impact of the road 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.


Introduction
The vehicle-road subsystem includes various aspects of the vehicle-road interaction. This subsystem is emphasized where the technical part of transport operations is investigated. Considerable emphasis has already been placed on the description of the 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][2][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 limitations have been introduced increasingly by a lot of countries in response to the aggressive development of this type of transport. These limitations have an effect both on individual components of the subsystem and on the subsystem itself and promote the improvement of both technical and organizational aspects thereof. The road and vehicle as the elements of the vehicle-road subsystem may be described individually as they are viewed and may be explored as individual elements of the system in an investigation of transport systems [4]. 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, consisting of a finite number of parameters, for investigation of stability. The required tire parameters were determined and revised using the theoretical and experimental methods. 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 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 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 stability, and the development of the recommendations on safe movement on the roads using the available characteristics.
Based on the topics discussed above, the main contributions of this paper are as follows: (1) we examine the vehicle stability dynamic models and develop a refined methodology for estimating the key parameters of these models; (2) we produce a revised tire smoothing function model that takes into account tread deflections; (3) we investigate vertical dynamics models and develop recommendations for their use in-vehicle stability studies; (4) we verify the validity of the model of the vertical dynamics using experimental studies; and (5) we establish the vehicle stability model, allowing evaluation of vehicles' technical performance and safe driving conditions of vehicles on roads with different characteristics.

Investigation of the Tire Properties
When dealing with the vehicle stability tasks, the issue of tire-road interaction is one of the most complex tasks. Analysis of the vehicle wheel-to-road contact requires taking into account the deformation properties of the tire. The complexity of the task stems from the uncertainty of the operating conditions of the tire and the complex structure of the tire. The complexity of the problem is defined by the mere fact that no consensus has been reached on the shape of the tire.
The performed studies of the tire-road interaction [36] confirmed that it would be necessary to consider the deformation properties of the band in the investigation of the compensating function of the band. For this purpose, an additional element is added to the quarter-car model ( Figure 2). Additional studies are needed in order to identify the 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 rubber parts are used here as they consider that ν = 0.5 [37]. The tread model is comprised of 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, consisting of a finite number of parameters, for investigation of stability. The required tire parameters were determined and revised using the theoretical and experimental methods. 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 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 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 stability, and the development of the recommendations on safe movement on the roads using the available characteristics.
Based on the topics discussed above, the main contributions of this paper are as follows: (1) we examine the vehicle stability dynamic models and develop a refined methodology for estimating the key parameters of these models; (2) we produce a revised tire smoothing function model that takes into account tread deflections; (3) we investigate vertical dynamics models and develop recommendations for their use in-vehicle stability studies; (4) we verify the validity of the model of the vertical dynamics using experimental studies; and (5) we establish the vehicle stability model, allowing evaluation of vehicles' technical performance and safe driving conditions of vehicles on roads with different characteristics.

Investigation of the Tire Properties
When dealing with the vehicle stability tasks, the issue of tire-road interaction is one of the most complex tasks. Analysis of the vehicle wheel-to-road contact requires taking into account the deformation properties of the tire. The complexity of the task stems from the uncertainty of the operating conditions of the tire and the complex structure of the tire. The complexity of the problem is defined by the mere fact that no consensus has been reached on the shape of the tire.
The performed studies of the tire-road interaction [36] confirmed that it would be necessary to consider the deformation properties of the band in the investigation of the compensating function of the band. For this purpose, an additional element is added to the quarter-car model ( Figure 2). Additional studies are needed in order to identify the deformation properties of the element.
Appl. Sci. 2021, 11, 6608 4 of 35 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 damping ratio of the tire tread part; c2, k2-tire stiffness and damping ratio; c3, k3-suspension stiffness and damping ratio; m2-unsprung mass; m3-sprung mass.
This paper analyzes radial tires of three types that differ considerably in terms of structure: passenger car (175/70R13; 195/50R15), light-duty truck (185/75R14C), and heavy-duty truck (12.00R20) tires. In the assessment of the vehicle tire-road interaction, it is first necessary to describe the deformation properties of the tire under the vertical load. The static load consists of the vehicle and cargo mass. The loaded tire is subject to deformation when forming the contact area with the road pavement. The tire-road contact is depicted in Figure 3, and tire element characteristics are presented in Table 1. 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 (ld ≈ l0) ( Figure 3a); • there is the linear dependency between the radial load of the tire and radial deformation Δh; • the pressure is distributed evenly across the contact surface.
The second and third assumptions enable determining the length of the tire-road contact: 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' (Figure 4a,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. Figure 2. Quarter-car model: L k -tire-road contact length; m 1 , c 1 and k 1 -mass, stiffness, and damping ratio of the tire tread part; c 2 , k 2 -tire stiffness and damping ratio; c 3 , k 3 -suspension stiffness and damping ratio; m 2 -unsprung mass; m 3 -sprung mass.
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 rubber parts are used here as they consider that ν = 0.5 [37]. The tread model is comprised of 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. This paper analyzes radial tires of three types that differ considerably in terms of structure: passenger car (175/70R13; 195/50R15), light-duty truck (185/75R14C), and heavy-duty truck (12.00R20) tires. In the assessment of the vehicle tire-road interaction, it is first necessary to describe the deformation properties of the tire under the vertical load. The static load consists of the vehicle and cargo mass. The loaded tire is subject to deformation when forming the contact area with the road pavement. The tire-road contact is depicted in Figure 3, and tire element characteristics are presented in Table 1. 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 (l d ≈ l 0 ) ( Figure 3a); • there is the linear dependency between the radial load of the tire and radial deformation ∆h; • the pressure is distributed evenly across the contact surface.
The second and third assumptions enable determining the length of the tire-road contact: where C R -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' (Figure 4a,b) [38]. The ring is considered to be the elastic band; the tensile force N 0 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 35     Where the impact of the cord threads is assessed, the radial deformation model is used, and the tire reaction is assessed using single pressure = . The sidewalls react only to the displacements in the radial direction, and the reaction per length unit equals = . An elastic outer layer with the thickness equal to H before deformation deforms by ΔH in the contact area. Inelastic forces are not subject to assessment. Tread reaction: Arc length x (Figure 4b) is measured from the center of contact. The main measures sought are [38]: • model deflection hz assuming that the band does not deform in the radial direction: load Fz: As a result of integration and rearrangement of Equation (4), the following is obtained: Where the impact of the cord threads is assessed, the radial deformation model is used, and the tire reaction is assessed using single pressure q z = C z w. The sidewalls react only to the displacements in the radial direction, and the reaction per length unit equals q s = C zs w. An elastic outer layer with the thickness equal to H before deformation deforms by ∆H in the contact area. Inelastic forces are not subject to assessment. Tread reaction: Arc length x (Figure 4b) is measured from the center of contact. The main measures sought are [38]: • model deflection h z assuming that the band does not deform in the radial direction: • contact zone length-2x 0 ; • contact pressure distribution-q z (x); • load F z : Appl. Sci. 2021, 11, 6608 As a result of integration and rearrangement of Equation (4), the following is obtained: The unknown measures h z and x 0 are determined using the following dependencies [38]: where c = cosθ and s = sinθ; x-direction forming an angle θ with the reinforcement direction ( Figure 5). The modulus of elasticity of the components and the Poisson's coefficients were taken from the works [36,43]. The methodology used for this purpose is presented in [39][40][41][42]. Direct modulus of elasticity calculations does not provide an answer in relation to the characteristics of the tread and unreinforced ring section that are required for the model provided in Figure 4. To determine the stiffness of the tread Cz, tread deformations are analyzed additionally. The obtained data are presented in Table 2.  The modulus of elasticity of the components and the Poisson's coefficients were taken from the works [36,43]. The methodology used for this purpose is presented in [39][40][41][42]. Direct modulus of elasticity calculations does not provide an answer in relation to the characteristics of the tread and unreinforced ring section that are required for the model provided in Figure 4. To determine the stiffness of the tread C z , tread deformations are analyzed additionally. The obtained data are presented in Table 2.

Investigation of Tread Deformation
Tread deformations were investigated using several models: porous composite and rubber prism models. The research studies analyzing the pneumatic tires specify fairly wide limits of the modulus of elasticity of rubber E r = 2 . . . 20 MPa. In the study [36], a simplified model of a passenger car tire was used, E r = 18 MPa, in order to obtain adequate deformation characteristics of the whole tire. In order to specify the initial characteristics, an experiment was performed to determine the 313 mm wide section of the tread 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 R znom 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 ( Figure 6) is very close to the linear (correlation coefficient r xy = 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, E r = 6.16 MPa, which was much lower than in the study [36]. Note: E-in all directions, E l -the rolling direction of the tire, E t -perpendicular to the rolling direction of the tire, and E r -in the radial direction (in relation to the tire).
Appl. Sci. 2021, 11, 6608 10 of 35 Whereas the tread consists of the parts with and without a pattern and the deformation 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 Whereas the tread consists of the parts with and without a pattern and the deformation 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]: where k u -tread filling factor, where k u = A i /A p , A i -the area of the continuous pattern in the part of the tread part selected according to the repetition of the pattern, A p -area of the selected part of the tread; and E g -tread rubber modulus of elasticity. Whereas the tread consists of the parts with and without a pattern and the deformation 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]: 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. 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: 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 elasticity of the whole tread was calculated using the equation in [44]: 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: where E g -tread rubber modulus of elasticity; k a = δ 1a /δ 1 ; δ 1 = δ 1a + δ 1b ; k bi -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 elasticity of the whole tread was calculated using the equation in [44]: This model does not accurately reflect the operating conditions of the tread. Deformation 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 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: 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]: 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 (Table 4). 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 (Table 2), 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.
Both models offer different characteristics: the modulus of elasticity E 1 or deformation ∆δ corresponding to the radial load R z value. For comparison of the results, the generalized characteristic describing the relationship between load and deformation-stiffness C 1 was used: where C 1 -tread stiffness and R z -radial load of the tire. When determining the measure C 1 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(R z ) dependence was obtained for all the tested tires ( Figure 8). The tire load is calculated using the nominal load of the tire R znom , which is specified in the catalogs, by setting the limits R zmin = 0.2R znom and R zmax = 1.4R znom . The both models were compared to the first model by calculating the following: Appl. Sci. 2021, 11, 6608 13 of 35 Figure 8. Radial loads of the tires (marking according to Table 2).

Figure 9.
Dependence of relative stiffness of the tread part with a pattern on the relative 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  Table 2).
Nonlinear dependence was observed C 1a (R z ) and was related to nonlinear dependence between R z and the length of the contact with the road. This generated the nonlinear q(R z ) 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 C 1a (R z ) has no significant effect on the nature of the dependence ( Figure 9).
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-R z 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.  Table 2).

Figure 9.
Dependence of relative stiffness of the tread part with a pattern on the relative 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. The tape for Nokian passenger car tires 175/70R13 82T was modeled. Its more detailed structure is provided in Table 5. Steel-rubber 0.9 7 Viscose-rubber 0.9 8 Sealing layer 0.9 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][40][41][42] by applying them to the axes presented in Figure 5:  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.
The tape for Nokian passenger car tires 175/70R13 82T was modeled. Its more detailed structure is provided in Table 5. 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][40][41][42] by applying them to the axes presented in Figure 5: 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 anisotropic, 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 identifying the rubber layer between the reinforced layers ( Figure 11a; Table 5) and using simplified separation into two layers (Figure 11b). identifying the rubber layer between the reinforced layers ( Figure 11a; Table 5) and using simplified separation into two layers (Figure 11b). In the investigation of the deformation properties of the whole band, the model used in layered composites was applied initially. The calculated modulus of elasticity of the band as a composite was calculated by taking into consideration the modulus of elasticity of the layers joined in parallel (X and Y directions) or consistently (Z direction).
Since the reinforcement angles of the adjacent breaker layers differed by the sign only (+20° and −20°), it was assumed that the resulting shear forces in the layers were counterbalanced and were not transmitted to other layers.
Considering the large difference between the modulus of elasticity of the reinforcing fiber and the rubber, it was determined that the modulus of elasticity of the layer in the direction of reinforcement was not significantly affected by the modulus of elasticity of the rubber. The method of identification of the layers (a and b) also did not have any effect in this case. At the same time, it was found that the model of elasticity of the layer would be more dependent on the modulus of elasticity of the rubber in the breaker fiber direction X and Z. The modulus of elasticity in X direction was less dependent on the total modulus of elasticity E. By changing the modulus of elasticity of the rubber in the range of 2…20 MPa, Ex differed by 3.7% for model (a) and 3.5% for model (b). Ez dependence was particularly evident: rubber stiffness Er would change 10 times, while Ez would change 9.97 times ( Figure 12). Considering that the deformation properties of the band in Z direction differed considerably from the data provided by other authors, the assumption was that stiffer cord and breaker layers may restrict tread deformations, and in the state of constrained deformation, tread deformation could be significantly influenced by the Poisson'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- In the investigation of the deformation properties of the whole band, the model used in layered composites was applied initially. The calculated modulus of elasticity of the band as a composite was calculated by taking into consideration the modulus of elasticity of the layers joined in parallel (X and Y directions) or consistently (Z direction).
Since the reinforcement angles of the adjacent breaker layers differed by the sign only (+20 • and −20 • ), it was assumed that the resulting shear forces in the layers were counterbalanced and were not transmitted to other layers.
Considering the large difference between the modulus of elasticity of the reinforcing fiber and the rubber, it was determined that the modulus of elasticity of the layer in the direction of reinforcement was not significantly affected by the modulus of elasticity of the rubber. The method of identification of the layers (a and b) also did not have any effect in this case. At the same time, it was found that the model of elasticity of the layer would be more dependent on the modulus of elasticity of the rubber in the breaker fiber direction X and Z. The modulus of elasticity in X direction was less dependent on the total modulus of elasticity E. By changing the modulus of elasticity of the rubber in the range of 2 . . . 20 MPa, E x differed by 3.7% for model (a) and 3.5% for model (b). E z dependence was particularly evident: rubber stiffness E r would change 10 times, while E z would change 9.97 times ( Figure 12). Considering that the deformation properties of the band in Z direction differed considerably from the data provided by other authors, the assumption was that stiffer cord and breaker layers may restrict tread deformations, and in the state of constrained deformation, tread deformation could be significantly influenced by the Poisson'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.
Appl. Sci. 2021, 11, 6608 16 of 35 The numerical experiment showed that the modulus of elasticity of rubber in Z direction would increase greatly and result in an increase in the modulus of elasticity of the entire band in Z direction ( Figure 13).  The numerical experiment confirmed that the models predicting uniform layer deformation 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: 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 direction would increase rapidly. The identified additional rubber layers had a smaller influence in this case ( Figure 14). 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: The sum of the resulting single axial forces in the individual layers in the X and Y axis directions equaled zero: The numerical experiment showed that the modulus of elasticity of rubber in Z direction would increase greatly and result in an increase in the modulus of elasticity of the entire band in Z direction ( Figure 13).
Appl. Sci. 2021, 11, 6608 16 of 35 The numerical experiment showed that the modulus of elasticity of rubber in Z direction would increase greatly and result in an increase in the modulus of elasticity of the entire band in Z direction ( Figure 13).  The numerical experiment confirmed that the models predicting uniform layer deformation 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- The numerical experiment confirmed that the models predicting uniform layer deformation 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: where µ-coefficient of friction between rubber and steel. The numerical experiment confirmed that friction had an impact on the E z measure ( Figure 14). As the coefficient of friction increased, the stiffness of the band in the transverse direction would increase rapidly. The identified additional rubber layers had a smaller influence in this case ( Figure 14). To revise the results obtained, it was attempted to deform the band cut from the tire in the Z direction. The experiment was performed using the 175/70R1382T Nokian tire band used for modeling. The summarized experimental data are presented in Figure 15. The experiment confirmed that the linear relationship between pressure and deformation 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 To revise the results obtained, it was attempted to deform the band cut from the tire in the Z direction. The experiment was performed using the 175/70R1382T Nokian tire band used for modeling. The summarized experimental data are presented in Figure 15. To revise the results obtained, it was attempted to deform the band cut from the tire in the Z direction. The experiment was performed using the 175/70R1382T Nokian tire band used for modeling. The summarized experimental data are presented in Figure 15. The experiment confirmed that the linear relationship between pressure and deformation 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 experiment confirmed that the linear relationship between pressure and deformation 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: where C * i -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 composite with sequentially arranged layers, its stiffness could be estimated by Equation (22). In this case, the member of the stiffness matrix of the tested tire tread E * z would be equal to 27 MPa. This value was used in the research for determination of the C z measure in Equation (22).
For real tread modeling, tread loads and operating conditions need to be adjusted by bringing the numerical model closer to real conditions. The third tread deformation model was used to evaluate the specifics of the tire-road texture interaction. Road texture means unevenness of the road pavement. For asphalt pavement, the height of the unevenness was 1-5 mm. The irregularities were assumed to be hemispherical protrusions of regular spacing ( Figure 16). This assumption enabled modeling of the interaction of an individual unevenness with the tread, assuming that the unevenness acted on the continuous part of the tread. The coefficient of filling with spheres necessary for the conversion of the reference pressure was equal to π/4 = 0.785. Based on the geometrical conditions, it was determined that the gaps between the spheres would be filled at the sphere penetration into the tread equal to: ∆ = 0.455r, where r-radius of the sphere. The known spheres used for the investigation of interaction were the solution for the elastic deformations of the plane contact. The interaction scheme is shown in Figure 17.
Appl. Sci. 2021, 11,6608 For real tread modeling, tread loads and operating conditions need to be adju bringing the numerical model closer to real conditions. The third tread deformation was used to evaluate the specifics of the tire-road texture interaction. Road texture unevenness of the road pavement. For asphalt pavement, the height of the une was 1-5 mm. The irregularities were assumed to be hemispherical protrusions of spacing ( Figure 16). This assumption enabled modeling of the interaction of an ind unevenness with the tread, assuming that the unevenness acted on the continuous the tread. The coefficient of filling with spheres necessary for the conversion of th ence pressure was equal to π/4 = 0.785. Based on the geometrical conditions, it wa mined that the gaps between the spheres would be filled at the sphere penetrat the tread equal to: Δ = 0.455r, where r-radius of the sphere. The known spheres u the investigation of interaction were the solution for the elastic deformations of th contact. The interaction scheme is shown in Figure 17.    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: 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: where Average contact pressure: F measure was to make sure that pressure q acting at the tire-road when the tire was subjected to force R z would be equal to the mean pressure at the sphere-to-tread contact q 0 , taking into account the coefficient of filling with spheres at the contact. Considering that, according to the elastic solution: and after introduction of the condition q 0 = q, the result would be: The tread stiffness values obtained for the individual models are presented in Figure 18. F measure was to make sure that pressure q acting at the tire-road when the tire was subjected to force Rz would be equal to the mean pressure at the sphere-to-tread contact q0, taking into account the coefficient of filling with spheres at the contact. Considering that, according to the elastic solution: (25) and after introduction of the condition q0 = q, the result would be: The tread stiffness values obtained for the individual models are presented in Figure 18. The calculations confirmed that in the case of 5 mm high unevenness, condition Δ < 0.455r (the tread does not fill the entire volume between the unevennesses) would be met for the pressures at the tire-road contact. This result demonstrates that the gaps between the road unevenness elements of this shape and size are not completely filled by the tread and enable better adhesion to the road and removal of moisture from the tire-road contact The calculations confirmed that in the case of 5 mm high unevenness, condition ∆ < 0.455r (the tread does not fill the entire volume between the unevennesses) would be met for the pressures at the tire-road contact. This result demonstrates that the gaps between the road unevenness elements of this shape and size are not completely filled by the tread and enable better adhesion to the road and removal of moisture from the tire-road contact area. Road texture interaction model (r = 5 mm) indicates lower stiffness values; however, dependence CR z remains similar and is related to nonlinear qR z dependence. In addition, it was confirmed that a simplified model estimating the modulus of elasticity of the tread band could be used for the assessment of the tread deformation. The modulus of elasticity of the tread material should be determined by assessing the deformation constraint.

Vehicle Excitation by Road Unevenness
Vehicle stability was assessed using several criteria. One of them was the stability of the linear motion of the steered wheels measured at a critical speed, above which the wheels started to vibrate about the vertical axis. The cause of the vibrations was the suspension vibrations. Therefore, the issue of car stability would be analyzed starting with the formation of the differential equations for a system consisting of a steering wheel mechanism, suspension, and wheels. At the initial moment, the number of elements comprising the system was reduced, i.e., the body vibrations were not taken into account. Body vibrations were not evaluated because its free vibration frequency was significantly lower than the suspension and wheel vibration frequency.
The simplest model used for road-vehicle interaction analysis is the quarter-car model. It consists of three masses: tire element m 1 , unsprung masses (suspensions) m 2 , and sprung masses (body) m 3 (Figure 2). This model enables the calculation of the maximum and characteristic values of the vertical displacements of the wheel and to specify the frequency characteristic of the suspension. More accurate models were found to require a more accurate assessment of tire performance, such as damping. The dimensions of the tire tread element were determined according to Figure 3, while the length of the contact was according to Equation (1) expression. The tread stiffness coefficient was determined using Equation (2). Quarter-car model behavior under kinematic excitation by road unevenness q was described using the following system of equations: where I x -moment of inertia of the sprung mass of the vehicle around the longitudinal axis of the tilt and t-wheel track of the respective axis. Damping matrix: Stiffness matrix: The pavement influence matrix must take into account that the excitation on the left and right sides may be unequal and was therefore expressed as: where q ∆1 = 1 t 1 (q 1k − q 2d ), q ∆2 = 1 t 2 (q 2k − q 2d ), with q ∆1 -and q ∆2 -angular differences in road unevenness at the front and rear axes of the vehicle, respectively.
The quarter-car and 3D car models (the characteristics of the model elements are presented in Table 6, and modeling results are presented in Figure 19) were compared numerically by modeling the compact car movement on a city-type road at the speed of 40 km/h and 80 km/h, taking into account the nonlinearity of the characteristics of the tire and the elastic members. Figure 19 shows that in the case of the 3D model, larger vertical displacements of the sprung masses (up to 20%) were obtained, but the overall curve characteristic remained almost without any changes.
where ∆ = ( − ), ∆ = ( − ), with qΔ1-and qΔ2-angular differences in road unevenness at the front and rear axes of the vehicle, respectively. The quarter-car and 3D car models (the characteristics of the model elements are presented in Table 6, and modeling results are presented in Figure 19) were compared numerically by modeling the compact car movement on a city-type road at the speed of 40 km/h and 80 km/h, taking into account the nonlinearity of the characteristics of the tire and the elastic members.   Figure 19 shows that in the case of the 3D model, larger vertical displacements of the sprung masses (up to 20%) were obtained, but the overall curve characteristic remained almost without any changes.

Excitation during Movement on the Asphalt Pavement
The influence of nonlinear characteristics of the elastic elements of suspension on the vertical displacements of the body and wheels caused by the excitation frequency was investigated during modeling of the stability of the compact car by combining different quality road profile and car speed, as well as assessing the influence of the tire tread. The results were compared by assessing the dependence of the coefficient of dynamism (ratio of sprung and unsprung mass displacements to the height of unevenness) on the road excitation. Figure 20a,b shows the influence of nonlinear damping of shock absorbers and nonlinear stiffness of suspension on the sprung and unsprung mass displacements.

Excitation during Movement on the Asphalt Pavement
The influence of nonlinear characteristics of the elastic elements of suspension on the vertical displacements of the body and wheels caused by the excitation frequency was investigated during modeling of the stability of the compact car by combining different quality road profile and car speed, as well as assessing the influence of the tire tread. The results were compared by assessing the dependence of the coefficient of dynamism (ratio of sprung and unsprung mass displacements to the height of unevenness) on the road excitation. Figure 20a,b shows the influence of nonlinear damping of shock absorbers and nonlinear stiffness of suspension on the sprung and unsprung mass displacements. With nonlinear elements, the frequency dependences shifted towards higher frequency because the damping of the shock absorbers due to nonlinearity increased suddenly, and the displacement (coefficient of dynamism) increased as well. The damping of the vibrations of the unsprung masses was fairly effective, but at a certain excitation frequency, the second peak would occur, causing certain damping deficiencies. Figure 21a,b shows the influence of road unevenness on the tread displacements of the sprung and unsprung masses and the tire tread. As the unevenness increased, the coefficient of dynamism could be observed to increase as well. As long as the amplitude was small, the suspension and shock absorber would operate in a linear section. As the amplitude increased, the influence of the nonlinearity of the elements would be observed (Figure 21b). With nonlinear elements, the frequency dependences shifted towards higher frequency because the damping of the shock absorbers due to nonlinearity increased suddenly, and the displacement (coefficient of dynamism) increased as well. The damping of the vibrations of the unsprung masses was fairly effective, but at a certain excitation frequency, the second peak would occur, causing certain damping deficiencies. Figure 21a,b shows the influence of road unevenness on the tread displacements of the sprung and unsprung masses and the tire tread. As the unevenness increased, the coefficient of dynamism could be observed to increase as well. As long as the amplitude was small, the suspension and shock absorber would operate in a linear section. As the amplitude increased, the influence of the nonlinearity of the elements would be observed (Figure 21b).    Damping of the sprung mass was effective at different speeds. The displacement reached its peak at low frequencies and afterward would decrease considerably. Tire deformations were the most sensitive to the effect of speed.

Experimental Investigation
In an investigation of a car stability, the characteristics of the road microprofile must be available. To measure the road microprofile, a measuring device consisting of the Volkswagen Transporter with DYNATEST road profile measuring equipment and laser sensors mounted on the front of the car were used (Figure 23). Damping of the sprung mass was effective at different speeds. The displacement reached its peak at low frequencies and afterward would decrease considerably. Tire deformations were the most sensitive to the effect of speed.

Experimental Investigation
In an investigation of a car stability, the characteristics of the road microprofile must be available. To measure the road microprofile, a measuring device consisting of the Volkswagen Transporter with DYNATEST road profile measuring equipment and laser sensors mounted on the front of the car were used (Figure 23). Measurement equipment: VAS-21 vibration analysis and processing system [59], PicoScope 3424 signal converter (Table 7), two low-frequency (WILCOXON Model 793L SN7688) and two medium-frequency (WILCOXON Model 784A SN12040) sensors ( Table  8). The data were registered real-time on a PC, recorded, and processed using PICO-SCOPE 5.13.7 software.   Measurement equipment: VAS-21 vibration analysis and processing system [59], PicoScope 3424 signal converter (Table 7), two low-frequency (WILCOXON Model 793L SN7688) and two medium-frequency (WILCOXON Model 784A SN12040) sensors ( Table 8). The data were registered real-time on a PC, recorded, and processed using PICOSCOPE 5.13.7 software. Measurements were performed in the city and on the highway. Road unevenness was measured using the Lithuanian Transport Safety Administration and Transport Competence Agency a VW Transporter car equipped with the measuring devices. The driving speeds in the city were and on the highway: 50, 60, 70, and 80 km/h. Measurement duration: 50 s.
To determine the frequency characteristics of the sprung and unsprung mass, compact car was selected for the study (Table 9). The measurement sensors for compact city car unsprung mass vibrations were mounted on the suspension arms (Figure 24a,b), while the sensors for measuring sprung mass vibrations were mounted on the front fender and car roof (Figure 24c,d).
Appl. Sci. 2021, 11, 6608 27 of 35 Measurements were performed in the city and on the highway. Road unevenness was measured using the Lithuanian Transport Safety Administration and Transport Competence Agency a VW Transporter car equipped with the measuring devices. The driving speeds in the city were and on the highway: 50, 60, 70, and 80 km/h. Measurement duration: 50 s.
To determine the frequency characteristics of the sprung and unsprung mass, compact car was selected for the study (Table 9). The measurement sensors for compact city car unsprung mass vibrations were mounted on the suspension arms (Figure 24a,b), while the sensors for measuring sprung mass vibrations were mounted on the front fender and car roof (Figure 24c,d).  Frequencies close to the frequency of own oscillations of the suspension (for passenger cars 0.8-1.3 Hz) and own oscillation frequency of unsprung masses (for passenger cars 8-13 Hz) were of particular interest in terms of the stability characteristic. It was the travel of the suspension that was defined as the difference between the displacements of the sprung and unsprung masses that had the greatest influence on the changes in the spatial position of the wheel. The equipment registered the peaks of the lowest frequencies in the range of 5-7 Hz (Figure 27), i.e., the frequencies characteristic of own oscillation frequency of the suspensions. In addition, the vibration sensors registered additional vibrations caused by the vibrations of the suspension elements (Figure 27a). Their frequencies exceeded the characteristic vibration frequencies of the tire (80-100 Hz). The oscillations were transmitted to the body ( Figure 27b); however, comparing the frequency spectra of sprung and unsprung masses, the oscillations in the body were observed to be influenced by the structure of elements selected as sensor mounting points, namely, the wing at the suspension fixing point was less sensitive to high-frequency vibrations than the roof panel. Frequencies close to the frequency of own oscillations of the suspension (for passenger cars 0.8-1.3 Hz) and own oscillation frequency of unsprung masses (for passenger cars 8-13 Hz) were of particular interest in terms of the stability characteristic. It was the travel of the suspension that was defined as the difference between the displacements of the sprung and unsprung masses that had the greatest influence on the changes in the spatial position of the wheel. The equipment registered the peaks of the lowest frequencies in the range of 5-7 Hz (Figure 27), i.e., the frequencies characteristic of own oscillation frequency of the suspensions. In addition, the vibration sensors registered additional vibrations caused by the vibrations of the suspension elements (Figure 27a). Their frequencies exceeded the characteristic vibration frequencies of the tire (80-100 Hz). The oscillations were transmitted to the body ( Figure 27b); however, comparing the frequency spectra of sprung and unsprung masses, the oscillations in the body were observed to be influenced by the structure of elements selected as sensor mounting points, namely, the wing at the suspension fixing point was less sensitive to high-frequency vibrations than the roof panel.  The aim of the experiment was to verify the modeling results and investigate their applicability to stability analysis. A Fourier series was used for frequency analysis of the results of direct measurements used for further studies (Figure 27). The Fourier analysis provides a wider spectrum but also does not highlight the characteristic body -and suspension-specific frequencies. This may be related to the specifics of the road. Further research identified how the results of simulation of the suspension operation could be used to define stability. The aim of the experiment was to verify the modeling results and investigate their applicability to stability analysis. A Fourier series was used for frequency analysis of the results of direct measurements used for further studies (Figure 27). The Fourier analysis provides a wider spectrum but also does not highlight the characteristic body-and suspension-specific frequencies. This may be related to the specifics of the road. Further research identified how the results of simulation of the suspension operation could be used to define stability. The aim of the experiment was to verify the modeling results and investigate their applicability to stability analysis. A Fourier series was used for frequency analysis of the results of direct measurements used for further studies (Figure 27). The Fourier analysis provides a wider spectrum but also does not highlight the characteristic body-and suspension-specific frequencies. This may be related to the specifics of the road. Further research identified how the results of simulation of the suspension operation could be used to define stability.

Comparison of Excitation from Road Unevenness and Experimental Study
During the modeling, it was assumed that the compact car was moving on a road of a known profile in the city and on the highway at the speed of 50, 60, 70, and 80 km/h. To compare the results of modeling and experiment, a 100 m long road section was used. Road unevenness of the examined section was measured using profilograph DYNATEST 5051 RSP. The vehicle and suspension data used for the modeling are presented in Table 7.
Matlab program (3D model with nonlinear suspension elements) was used for the modeling. The modeling results are presented in Figures 28-30

Comparison of Excitation from Road Unevenness and Experimental Study
During the modeling, it was assumed that the compact car was moving on a road of a known profile in the city and on the highway at the speed of 50, 60, 70, and 80 km/h. To compare the results of modeling and experiment, a 100 m long road section was used. Road unevenness of the examined section was measured using profilograph DYNATEST 5051 RSP. The vehicle and suspension data used for the modeling are presented in Table  7. Matlab program (3D model with nonlinear suspension elements) was used for the modeling. The modeling results are presented in Figures 28-30. Comparison of the curves of experimental and numerical modeling results suggested that in both cases very similar acceleration values were obtained (acceleration peak values differed by 5-7% for sprung masses and 20-25% for unnsprung masses). Comparison of the dependences of accelerations by means of Fourier analysis ( Figure 29) with the experimental results suggested that a share of the vibrations caused by the resonance and stability of the suspension and body parts may be omitted from the assessment in the research findings. Comparison of the curves of experimental and numerical modeling results suggested that in both cases very similar acceleration values were obtained (acceleration peak values differed by 5-7% for sprung masses and 20-25% for unnsprung masses). Comparison of the dependences of accelerations by means of Fourier analysis ( Figure 29) with the experimental results suggested that a share of the vibrations caused by the resonance and stability of the suspension and body parts may be omitted from the assessment in the research findings. The obtained results show that the recording and analysis of signals required for active suspensions is an area of separate studies, as the additional noise and signal sources were recorded during the experiment, which also interfered with the analysis of directional stability. Therefore, the differences between the displacements of the sprung and The obtained results show that the recording and analysis of signals required for active suspensions is an area of separate studies, as the additional noise and signal sources were recorded during the experiment, which also interfered with the analysis of directional stability. Therefore, the differences between the displacements of the sprung and unsprung masses were examined for further analysis (Figure 30). The maximum amplitude of the difference between the suspension and body displacements was 47 mm for the experiment and 40 mm for the modeling. The obtained results show that the recording and analysis of signals required for active suspensions is an area of separate studies, as the additional noise and signal sources were recorded during the experiment, which also interfered with the analysis of directional stability. Therefore, the differences between the displacements of the sprung and unsprung masses were examined for further analysis (Figure 30). The maximum amplitude of the difference between the suspension and body displacements was 47 mm for the experiment and 40 mm for the modeling.

Discussion
The analysis of the tire tread as the elastic band model using the numerical models has shown that due to the specificity of rubber deformability, the cord and breaker layers have effect on its stiffness in the radial direction by restricting the tread deformations.
The experiment confirmed the theoretical assumptions and showed that the reduced modulus of elasticity of the tread was 4.5 times higher than the modulus of elasticity of the tread material.
It was found that when assessing the tread deformation for different types of tires, summative expressions could be used for the relative coordinates as they use the relative load determined on the basis of the nominal tire load, thus greatly simplifying the assessment of tire deformation.

Discussion
The analysis of the tire tread as the elastic band model using the numerical models has shown that due to the specificity of rubber deformability, the cord and breaker layers have effect on its stiffness in the radial direction by restricting the tread deformations.
The experiment confirmed the theoretical assumptions and showed that the reduced modulus of elasticity of the tread was 4.5 times higher than the modulus of elasticity of the tread material.
It was found that when assessing the tread deformation for different types of tires, summative expressions could be used for the relative coordinates as they use the relative load determined on the basis of the nominal tire load, thus greatly simplifying the assessment of tire deformation.
The analysis of the interaction between the tread and the road texture showed that the tread would even out the unevenness in the texture level, but the depressions would not be completely filled, thereby ensuring better adhesion between the tire and the road.
The specifics of the tire in the road-to-car interaction were determined. The deformation characteristics of the tread in the radial direction did not depend much on the tire or, consequently, on the type of vehicle.
Comparison of the quarter-car and 3D models showed that more detailed information necessary for stability assessment was provided by the 3D model, which predicted the sprung mass displacements higher by up to 20%.
The results of the experiment on the asphalt pavement were compared with the model. Certain issues were observed in the analysis of experimental results due to additional oscillations and vibrations occurring in the body and suspension elements, which could be avoided by modeling the movement after recording the road microprofile with a profilograph. The parameters that have the greatest impact on the stability of the direction-the difference between the displacement of the sprung and unsprung masses in the model versus the field experiment-did not exceed 15%.

Further Investigations
Based on the tire deformation properties obtained in this paper, numerically and experimentally investigated vertical excitation results, the changes in wheel position due to road roughness during suspension excitation can be determined and an analysis of the vehicle's yaw process can be performed. The modeling of the vehicle's directional stability needs to consider the specifics of the suspensions, such as the redistribution of radial loads under lateral force and the changes in wheel geometry due to the kinematic properties of the suspension. The obtained results justify further refinements of the directional stability numerical models: • development of a dynamic model of a comprehensive integrated assessment of the suspension characteristics, suspension kinematics, steering mechanism, and pavement condition. The model would provide the development of a methodology for determination of the vehicles handling characteristics on roads of known quality using the numerical models. • investigation of the sensitivity of the individual suspension types to vertical excitation in terms of directional stability criteria and to evaluate the technical condition of the vehicle if the tendency of changes in suspension parameters during known operation.
To evaluate the technical condition of the vehicle by modifying the stiffness of the suspension components. • analysis of vehicle stability in cornering, taking into account the changing position of the wheels as the vehicle is steered, which may have an effect on the vehicle stability. The stability of a vehicle in a turning situation requires the consideration of two criteria: the loss of stability of the vehicle when skidding, and the loss of stability of the vehicle when overturning. In addition, it is necessary to specify the effect of changes in the vehicle's trajectory in a corner due to the deformation of the suspension and the changes in the position of the wheels as the vehicle is steered. • integrating driver reactions into the vehicle stability model. Without a driver model, such a simulation is not effective, as the vehicle model changes direction due to random effects and requires frequent correction to simulate steering course corrections.

1.
The revised model of the tire equalizing function was developed. The coordinates of relative measures were found to cause no essential difference in assessment of tire deformations of different vehicle tires. This means that the number of tire groups may be reduced for assessment of their impact on stability.

2.
The tread deformation characteristics of the tire were investigated by numerical modeling and experiment and showed that the tread deformations could be assessed by simplified models instead of the modulus of elasticity of the tread material, E r = 6 MPa using the calculated modulus of elasticity E 1z = 27 MPa.

3.
Studies on different types of cars (passenger cars, light-duty trucks and SUVs) have shown that the 3D suspension models more accurately provide reflection of the real situation in terms of stability assessment, while simpler quarter-car models can be used for comfort assessment. Comparison of the quarter-car model with the 3D model has shown that, in case of the latter, the resulting wheel travel is higher by up to 20%.

4.
The results of the experiment on the asphalt pavement were compared with the model. Certain issues were observed in the analysis of experimental results due to additional oscillations and vibrations occurring in the body and suspension elements, which could be avoided by modeling the movement after recording the road microprofile with a profilograph. The parameters that have the greatest impact on the stability of the direction-the difference between the displacement of the sprung and unsprung masses in the model versus the field experiment-did not exceed 15%.

5.
Changes in the car behavior at changing speed are related not only to the slip effect but also to changes in the wheel spatial position due to the suspension kinematics and suspension incompatibility with the steering mechanism. This effect can be particularly pronounced when road defects (pits, ruts) lead to an increase in the suspension travel.