Semi-Active System of Vehicle Vibration Damping

: This article discusses the vibration of a vehicle equipped with four semi-active dampers. The friction forces generated in these dampers have crucial influence on the intensity of the vehicle vibration. To estimate the effect of the vehicle vibration damping, a criterial function depending on the values of these forces has been selected. Solving the problem has the purpose of determining four values of the forces at which the criteria function reaches the minimum at any moment. On the basis of the above, signals for the devices controlling the operation of the semi-active dampers are determined. In this work, the analysis of the discussed problem was conducted in regard to the method for determining the optimal force values. Moreover, the results of the simulation of the vibration of the vehicle equipped with semi-active magneto-rheological dampers are presented. In addition, the analysis of the influence of the form of the criterial function on the vehicle vibration damping efficiency is proposed. Based on the tests carried out on the vehicle with MR dampers, the comfort index improved by 51.6%, compared to the classic suspension.


Introduction
This article discusses the vehicle vibration induced by uneven road surfaces and also by the inertial forces affecting the body. As a result of such a vibration, the driving comfort deteriorates and the dynamic load on the vehicle structure increases. Moreover, wheel pressure on the road surface changes (change in the dynamic loads in the suspension), which can contribute to slippage during acceleration, braking, or driving on a curve. In order to reduce vehicle body vibration, various kinds of devices are commonly used, dispersing the energy, known as vibration dampers.
The problem area related to applying the semi-active dampers to reduce the vibration of vehicles, machines or buildings has long been considered. However, in recent years, technical conditions have emerged that enable practical implementation of such a way of reducing vehicle vibration [1,2]. There has been development in the field of controlled structures based on smart materials. These materials include magnetorheological fluids [3,4], electrorheological fluids [5,6] and vacuum-packed particles [7].
At the end of the 20th century, magnetorheological fluid technology was developed. That laid the groundwork for the development of magnetorheological dampers (MR dampers). In the beginning, the technology was developed by the LORD Corporation, Cary, North Carolina, USA. MR dampers have been used in human prostheses and vehicle seats [8,9]. In the following years, further research related to the development of MR dampers began, which resulted in new solutions presented in publications [10][11][12]. Currently, controlled MR dampers are used in different fields: damping vibrations of buildings [13,14], reducing vibrations of washing machines [15,16], in suspensions of highspeed trains [17,18] and in automotive suspensions [19][20][21]. These conditions were created due to the possibility of using on-board computers with high computing power and semiactive dampers.
An important problem is the determination of the control signal for semi-active dampers. For this purpose, different methods were developed where the control signal is determined on the basis of the measurement signals. In multi-body systems, Zero Residual Vibration-Derivative methods [22,23] are used to determine the control signal. Another solution to the problem is to use Hybrid Prediction Control [24,25]. Sky-hook and ground-hook methods [26,27] are common in semi-active suppression of vehicle vibration. The authors of the work developed original methods of minimizing acceleration in a selected point of the car body [11].
In the presented system, the vibration control of the vehicle is carried out by means of semi-active dampers based on the control signals from the LVDT displacement sensors. The original mechatronic system, based on a dedicated algorithm for determining the control signals, was developed for control. Controlling the frictional force in vibration dampers was adopted to reduce the vertical accelerations affecting the vehicle passengers. In the discussed solution, a controlled magnetorheological damper (MR) was used.

MR Damper
A magnetorheological damper (MR) has a varying damping force that changes due to the changing current values. In such a damper, an enclosed coil generates the magnetic field. Then, the change in the friction force is induced by means of controlling the intensity of the magnetic field of the damper coil. One feature of magnetorheological fluid is its changing viscosity ensued by the magnetic field. Resulting from variations in the viscosity of the MR fluid, a change in the damping forces takes place. Thus, changes in the current intensity in the electronic system enable controlling changes in the dissipated energy. The developed MR damper was subjected to tests. Figure 1 shows the original design of the magnetorheological damper (MR-IP damper), developed at the Institute of Vehicles, Warsaw University of Technology (Instytut Pojazdów Politechniki Warszawskiej). The damper was filled with an MR fluid (3), the viscosity of which varies according to the magnetic field (5), generated by the coil (1), and enclosed in the piston (6). During the piston motion, the fluid is pressed through the orifices (2) in the piston. In the discussed solution, the housing (4) does not move; only the piston gets displaced and (6) mounted onto the piston rod (7). At that point, the change in the current intensity affects the change in the damping force acting on the damper's piston rod. The developed MR-IP damper design served the purpose of conducting the property tests. As a result of the experimental research, the dissipative characteristic of the MR damper were determined. Figure 2 illustrates the conducted tests on the plane: force-displacement and force-velocity, where the coil of the MR damper was supplied with the 2 A current, and with no supply (0 A). Experimental research was conducted with the kinematic excitation at the frequency of 0.8 Hz with 0.18 m amplitude. The diagrams show the change in the frictional force value in the case of tests conducted under no-supply conditions (0 A). The frictional force value at the level of 500 N was obtained, whereas in the case of the tests under the condition of the coil supply 2 A, the frictional force value was at the level of 1500 N. As far as the conducted tests are concerned, it has been postulated that controlling the frictional force in dampers is possible in the course of changing the supply current. Then, during the varying current intensity in the range of 0-2 A, controlling the frictional force value will be possible. In the shown force-velocity characteristics, a hysteresis loop can be observed with regard to the tests conducted under no-supply conditions, as well as with the 2 A current. Figure 2 shows only the limit values of the frictional forces. The change in the friction forces in the range of 0-2 A is not linearly dependent on the current changes, which should be taken into account in the algorithm for selecting the friction forces.
The test results served the purpose of developing the model of the controlled MR damper. In the subject literature, the Bouc-Wen model [28][29][30] is most frequently used to model this type of damper. In this work, a model in the form of a rheological structure was used, presented in Figure 3, originally proposed by Grzesikiewicz. This structure shows the viscoplastic features of the MR damper. The discussed model is characterized by four parameters that are relatively easy to determine on the basis of the experimental research results.
τ ∈ sign y , −1, 1 , when when where T , c , c , k-model parameters, which are positive numbers characterizing the viscous-elastic features of the structure, x, y-coordinates of the model. The values of the parameters (T , c , c , k) are related to the current of the MR damper. In the presented MR damper model, the friction force depends on the velocity of deformation x .
On the basis of the experimental research results, the identification of the damper parameters was conducted. In the course of the identification of the model parameters, it was observed that, along with the change in the current within the range of 0-2 A, the T parameter changes the most. For the sake of numerical tests, the assumption was adopted that the changing current will ensue the variation of only one parameter with the remaining parameters being unchanged. Figure 4 shows the results of the numerical research that corresponds to the experimental research at 0 A and 2 A.  On the basis of the conducted numerical research, the dissipative characteristics of the MR damper model were obtained. Similar to the characteristics regarding the experimental research, ranges of approx. 500 N of the frictional force were obtained at the current intensity 0 A, whereas in the case of the 2 A supply, the value of frictional forces of approx. 1500 N was obtained. The hysteresis loop can also be discerned, which could be observed in the results of the experimental research. Figure 5 shows the comparison of results of the experimental and numerical research on the MR damper with 2A current. As a result of the conducted tests of the MR damper and identification of the damper model, a numerical model was developed. The performed simulation tests validated the assumption that the proposed model may be utilized in numerical research of the MR damper. The determined damper parameters can be used in a model to simulate the tests of the damper supplied with the variable current value in the range of 0-2 A. The proposed damper model was used during further numerical research connected with reduction of the vehicle vibrations. Considering the inconclusive solutions that can be obtained on the basis of the dissipative characteristics of the damper, regularization (simplification of the characteristic) was carried out. This characteristic was used to control the value of damping forces.
The energy-related features of the semi-active damper are characterized by means of the parametric relationship (f) between the force generated in the damper (friction force) (T) and the velocity of its deformation (v), whereas the form of this relationship depends on parameter τ ∈ [0,1], which can be written using the following formula: Schematic plot of such a characteristic is shown in Figure 6. In the relationship under consideration, parameter represents a physical quantity by which the semi-active damper's dissipative features can be controlled, which is shown in Figure 4. The description above ensures that friction force T in the semi-active damper can assume values from the range determined by the following relationship: if Relationships (2) and (3) can be written in an altered form: By means of Formula (8), delineation Ω is defined as follows: if T ∈ R -vector of the friction forces in four dampers, V ∈ R -vector of the deformation velocity of four dampers.
Delineation Ω determines the set of permissible vectors of the friction forces T ∈ R depending on the velocity vector V ∈ R defining velocity values of damper deformation.
The optimization problem discussed subsequently is connected with determining the optimal vector T ∈ Ω(V) for which the criterial function : R → R reaches the minimal value. Thus, the vector is determined, fitting into the following relationship: A detailed description of this relationship will be given further on in this paper.

Vehicle Model
The considered mathematical vehicle model serves the purpose of analyzing the vibration in the vertical plane. Figure 7 shows the scheme of the mechanical system, which was adopted for the equations describing the vehicle vibration to be formulated.
In order to simplify this description, the following were adopted: − Elastic suspension elements have linear characteristics; − Vehicle wheels are in constant contact with the road surface; − Elastic-dissipative features of the suspension are mapped with the help of elements located directly above the wheels; − The vehicle mass is distributed symmetrically relative to the longitudinal axis.
To describe the vehicle vibration, coordinates were selected that determine the vehicle displacements relative to the equilibrium position. Figure 7 shows the selected marked coordinates and the basic model parameters. The equations of the vehicle vibration are defined with the following formulae: where S , T , i = 1, … ,4-elastic and dissipative forces acting between the i-th wheel and the suspension, so that forces resulting from the tire deformation M , M -moments of the force resulting from twisting of the front and rear stabilizers.
if κ , κ -rigidity of the front and rear stabilizers; The stabilizers' forces acting on the wheel The remaining symbols in the figure denote the following: V -vehicle velocity. M , Mmoments of the inertial force acting on the vehicle body during transient Let it be assumed, that functions ξ , i = 1, … ,4 describing kinematic vibration forcing, as well as the function of velocity V , are set. The above-formulated task regarding vehicle vibrations has a solution only when the vector of parameters τ ∈ R ( Figure 6) is known at any moment of time; then, on the basis of Formula (4), the values of forces in the dampers can be determined as follows: Parameters τ (i = 1, … ,4) are set in the system of damper control. The description of this systems are presented in the next section.
In the considerations therein, the following denotations will be used: h ∶= An observation can be made that the equality condition is satisfied: Using the above denotations in Equation (11) enabled the obtaining of the set of three equations describing the acceleration of the vehicle body: In accordance with the adopted assumptions, the vector of force S ∈ R (Formula (9)), determines changes in the load of the elastic elements of the suspension relative to their static load in the state of equilibrium. Figure 7. Scheme of the mechanical system adopted as a vehicle model; S-center of mass; Z-displacement of the bodywork in the Z direction; Φx-rotation about the X-axis; Φy-rotation about the Y-axis; m-mass of the body; Jx-moment of inertia about the X-axis; Jy-moment of inertia about the Y-axis; T1, T2, T3, T4-friction forces in suspension; k1, k2, k3, k4-spring coefficients in suspension; m0-mass of wheel; z01, z02, z03, z04-displacement of the wheels in the Z direction; k0-spring coefficients in wheel; c0-damping coefficients in wheel; , , , -kinematic excitations; Vvehicle-vehicle velocity; a1, a2, b1, b2-position of the suspension in relation to the center of mass; 1, 2, 3, 4-acceleration measurement points.

Model of the Control System
In the system of the damper control, signals are processed that represent the physical quantities characterizing the vehicle state. In the subsequent discussion, denotations of physical quantities and corresponding signals are the same.
The central part of the considered control system is a programmable controller, where the signal vector τ ∈ R is determined, corresponding to the optimal vector of the friction forces T defined by the relationship in Formula (9). Fundamentals applied to formulate this relationship are provided by criterial function , defining the influence of the friction forces on the intensity of the vehicle vibration. In this work, the intensity of the vehicle vibration will be estimated on the basis of vector A ∈ R defining the acceleration of the vehicle body: where vector of forces P ∈ R described with the Formula (14), M-matrix of the body inertia.
Formula (30) determines the relationship between the vectors of friction forces T ∈ R and accelerations of the vehicle body A ∈ R .
On the basis of the above, the first criterial function is determined, which sets the norm for the acceleration vector A ∈ R according to the following formula: Hence, after taking into account Formula (30), the following is obtained: ( ) (T) = (P + HT) M (P + HT), and after transformations, there are: The optimization task is formulated in the following way: vector T ∈ Ω(V) needs to be determined for which the value of function ( ) is the smallest in the set of vectors Ω(V) described with Formula (8).
Taking into account that function ( ) and set Ω(V) are convex allows for concluding that the above task has a solution; because matrix D is not strictly convex-which ensues from Formulae (25)  Thus, the optimizing issue considered herein comes down to determining a subset of vectors T, as mentioned before, and selecting one of those T . This task can be rendered as follows: For the vector of forces T ∈ Ω(V) so determined, the norm of the vector of body accelerations, denoted with Formula (32), reaches the minimum. It needs to be noted that vector T , selected in such a way, changes along the vector U, V ∈ R because the force vector P ∈ R (Formula (26)) depends on them and the set Ω (Formula (8)). This means, that relationship (37) is non-stationary because vectors U, V change over time.
The method of determining solution (37) can be found further in this work. Here, the second criterial function is presented, defined on the basis of the Karnopp principle [26,27] according to which effective damping of the vehicle vibration can be achieved by means of the forces that depend on the vector of the body velocity V : The force determined in such a way is known as sky hook. Figure 8 shows a schematic model of the vehicle with the damper control according to the sky-hook strategy.
The relationship of such a force (F ) by means of four semi-active dampers is possible when the following relationships are fulfilled: The relationships (41), however, can be fulfilled only in particular states of the vehicle body motion. For this reason, selecting of force T ∈ R is considered, which is the one that best maps the action of force F ∈ R . With the task determined in such a way, the second criterial function can be defined: that sets the measure for the vector deviation T ∈ R from the set of solutions to equation HT − F = 0. After the appropriate transformation of Formula (42), the following is obtained: The observation can be made that the form of criterial function is analogous to the form of criterial function . Moreover, criterial function is also convex but not strictly convex and also non-stationary because force F depends on the body velocity. Thus, vector T ∈ Ω(V), which optimally translates the action of force F ∈ R , is determined by means of the relationship from Formula (37): In further considerations, a method for determination of optimal vectors of friction forces, defined by relationships (37) and (45), is presented. Figure 8. Scheme of the mechanical system adopted as a vehicle model-sky hook; Z-displacement of the bodywork in the Z direction; Φx-rotation about the X axis; Φy-rotation about the Y axis; cz-damping coefficient relative to the sky in the Z direction; cx-angular damping coefficient relative to the sky in about the X axis; cy-angular damping coefficient relative to the sky in about the Y axis; T1, T2, T3, T4-friction forces in suspension; k1, k2, k3, k4-spring coefficients in suspension; m0-mass of wheel; z01, z02, z03, z04-displacement of the wheels in the Z direction; k0-spring coefficients in wheel; c0-damping coefficients in wheel; , , , -kinematic excitations; Vvehicle-vehicle velocity.

Solution to the Optimization Task
A homogenous description of both abovementioned tasks is considered here. Let the following criterial function be the focus: the matrix of which D ∈ R is semi-positive and its nucleus is determined by the onedimensional space: where vector h ∈ R is defined in Formula (24). The optimal solution is sought in convex set Ω(V) ⊂ R and is determined by the following relationship: The set of optimal vectors so determined is not empty, and each pair of vectors satisfies the condition T ( ) − T ( ) ∈ kerD.
On the basis of the Kuhn-Tucker theorem, and after taking into account the specific form of the set Ω (8), the conditions are obtained, which should be satisfied by relationship solution (48): λ ∈ R -of such vectors, that λ ≥ 0, i = 1, … ,4; ρ ∈ R -any positive number; π : R → R -mapping of the projection on the set Ω; The schematic plot illustrating relationship (51) is shown in Figure 9. The conditions from Formula (46) can be used to determine the iterative method of solving the optimization problem (45). In the first step, vector T * ∈ R is determined, which belongs to the set of stationary points of the function (43): In the case of criterial function ( ) (16), this set has the following form: and for criterial function (22), it is the following: Vector T * with the smallest norm in this set is defined by the following formula: where D ∈ R is the Moore-Penrose inverse of a matrix. For criterial functions ( ) and there are, respectively, the following: if and vector P is defined by Formula (12). If T * ∈ Ω, then vector T * determines the sought solution to the relationship (27).
If T * ∉ Ω(V), then the solution to relationship (45) or (46) is determined by means of the approximate iterative method according to the following scheme: and the condition of termination of the iteration has the form: where ε > 0 the number determining accuracy of the solution.
Equations (59) and (60) can be replaced with the one following equation: The solution of the relationship (45) shown earlier is determined by vector T ∈ Ω(V), for which criterial function reaches the smallest value. After determining vector Formula (9) enables determining the corresponding vector of parameters τ ∈ [0,1] , that is, where f function inverse relative to the second argument. Vector τ defines the signals controlling the working system with the help of which a physical quantity can be altered that affects the dissipative features of the semi-active damper. The assumption is that operation of the effective element controlled with signal τ is described by the following equation: where θ-time constant characterizing the working element.

Processing of Signals in the Control System
This section presents the principle of signal processing in the discussed control system. First, the signals related to the first criterial function ( ) are considered. In this case, the input signal to the control system are vectors U ∈ R , V ∈ R , representing displacements and velocities described by Formulae (8). The signals enumerated are processed in the programmable controller, whose microprocessor determines the solution to the optimization task in the form of signals τ ∈ R . Next, the effective element alters the parametric values τ , i = 1, … ,4 in accordance with Equation (65); as a result, the forces defined by Formula (9) are operative in the dampers.
It should be noted that the form of criterial function ( ) depends on the inertial parameters of the vehicle body, determined by the matrix of inertia M ∈ R , which is dependent upon the number of passengers and the mass of the load carried. Therefore, in order to correct matrix M, signal U can be used, determined at the initial moment.
In the case of control performed in accordance with criterial function , signal processing is analogous to the one described above. However, in such a case, the input signals are vector V ∈ R , described by Formula (9), and signal V ∈ R , determining the velocity of the vehicle body, described by Formula (34).

Simulation Research of Vehicle Vibration
Vibration simulations for a medium-sized passenger car are considered here, the model of which was mentioned in Section 3. The following parameter values regarding this vehicle have been adopted: m = 1250 kg-mass of the vehicle body; m = 25 kgreduced mass of the wheel and suspension elements; J = mρ , ρ = 0,6 m-moment of the body inertia relative to the longitudinal axis; J = mρ , ρ = 1,15 m-of the body inertia relative to the transverse axis; a = 1,4 m, a = 1,45 m, b = b = 0,725 m-measurements defining location of the body mass center relative to the wheels; h = 0,5 mheight of the vehicle mass center above the road surface ( Figure 8); k = k = 14,5 kN/m-rigidity of front and rear wheel springs; k = 200 kN/m-stiffness of the wheel tire; c = 2,5 • 10 Ns/m-coefficient of tire damping; κ = κ = 1,5 • 10 Ns/radrigidity of front and rear stabilizers. The vehicle considered is equipped with the magnetorheological dampers (MR dampers) whose characteristic is shown in Figure 6. The following parameters were adopted, shown in the Figure  In the MR damper, the parametric variable represents the current flowing through the winding of the coil located at the damper's side [11,31]. On the basis of the presented data, the damper characteristics have been formulated (1) as follows: where I-intensity of the current flowing through the coil winding. Moreover, the assumption was made that the vehicle moves along the road, the roughness of which is described by the following function: ξ(s) = ξ sin2π , where ξ = 5 mm, L = 12 m, s-the stretch of road covered. The vehicle is moving at variable velocity, the plot of which is shown in Figure 11. Another assumption was that the road's uneven surface under the left-side and right-side wheels is not exactly the same but shifted in phase by ∆s = a . In the presented paper, harmonic excitation was used. In this type of research, excitation based on the road profile in accordance with ISO 8608 [32,33] is often used. The assumption was made that parameters of the virtual damper, included in Formula (22), are determined on the basis of non-dimensional damping coefficients amounting to γ = 0,413; γ = 0,814; γ = 0,470. Due to the comparative character of the performed tests, it was assumed that the reference vehicle is equipped with classic shock absorbers (dampers) whose bi-directional characteristic is determined by two constants: c = 500 Ns/m when the shock absorber is squeezed, and c = 2,5 • c when the shock absorber is stretched. While using the presented model, the vehicle motion on the described road was simulated with the velocity given in Figure 11, and the drive time amounted to 30 s. It was assumed that at the initial moment (t = 0), all coordinates and corresponding velocities (see: Equation (10)) were equal to zero. Figure 12 shows a scheme of the supply system for the MR damper coil, i.e., of the system featuring the parametric variable τ, which is the current flowing through the coil.
The waveforms for the currents I , i = 1, … ,4 flowing through coils in the dampers, determined according to the system of equations. The following parameter values of the coil supply systems were adopted R = 2 Ω, L = 6 mH, I = 2 A.
where I is computed on the basis of the damper characteristic described by Formula (69): To evaluate the body vibrations, an indicator was adopted whose values determine the sum of the acceleration squares in the four points of the body above the wheels ( Figure 6): where a -accelerations of the enumerated points, i = 1, … ,4. Additionally, the value of an indicator characterizing change in the wheels' pressure on the road surface was also determined: if N (t) = (S (t) + T (t)) , i = 1, … , 4, where forces are described by Formula (16), and Q -static wheel pressure. Three vibration simulations were conducted for the vehicle moving under the conditions described above. During the first simulation, the vehicle was equipped with classic dampers whose characteristic is given earlier in point 4. The Figure 13 shows the change in the wheels' pressure on the road surface on the front axle, and Figure 14 shows the change in vertical accelerations over the front axle.  The second simulation was conducted with regard to the vehicle equipped with the MR dampers whose characteristic is given in Figure 11; optimal values of forces in the dampers were determined in this simulation in accordance with the first criterial function ( ) . The changes in pressure on the road surface on the front axle is shown in Figure 15, and the vertical acceleration on the front axle is shown in Figure 16.  The third simulation was analogous to the second, but the forces in the dampers were determined in accordance with the first criterial function . Figure 17 shows the change in the wheels' pressure on the road surface on the front axle, and Figure 18 shows the change in vertical acceleration over the front axle.  The purpose of this simulation research was the effectiveness comparison of reduction in the vehicle vibration in the three situations described above. Evaluation of this effectiveness was performed based on indicators W and W , which characterize the intensity of the vibration. At the same time, the smaller the value of the indicator, the greater the efficiency of the vibration reduction. Table 1 presents the values of the considered indicators. The results obtained show that the values of indicators W -characterizing the vibration intensity with regard to comfort-are significantly smaller for a vehicle equipped with the MR dampers than for a vehicle with classic dampers. In the ( ) system, the W indicator has shown an improvement of 51.6%, compared to the classic suspension, and in the ( ) system, one of 42.5%.
However, the influence of the form of the criterial function used in the control system on the value W is insignificant.
In the case of indicator W -determining changes in the wheel vertical forces-the impact of the MR dampers controlled in accordance with the criterion ( ) , is little (4.9% reduction). However, as a result of the control according to the criterion, there was a 90.8% increase in the W value, which indicates that the conditions of impact of the wheels on the road surface deteriorated.

Conclusions
This paper focuses on the optimization issue related to the semi-active system for damping the vehicle vibration. The problem was formulated and the method for determining the optimal values of the forces damping the vehicle vibration proposed. In the next step, two vibration simulations of the vehicle equipped with the semi-active vibration damping system were conducted. In the course of computations, there were optimal values of the four forces determined at all moments of time. Additionally, the simulation of the vibration of the vehicle equipped with the classic dampers was performed. On the basis of the simulation results, the values of the indicators used to evaluate the vibration intensity of a vehicle moving along the established route were determined. Analysis of these results showed that due to the use of the semi-active system for vibration damping, a substantial decrease in the body accelerations occurred, which contributes to improvement in ride comfort. In a vehicle with controlled MR dampers by the criterial function ( ) , the W1 indicator was improved by 51.6%, compared to the classic suspension, while maintaining a constant W2 dynamic load index. The semi-active system of reduction did not, however, significantly affect the minimizing of the change in the wheel pressure on the road surface. The lack of this effect isin the authors' opinions-caused by the active inertia forces arising during the acceleration or braking of the vehicle. Such slow-changing vehicle loads cannot be reduced by the action of semi-active vibration dampers.
The work presented in the article was carried out to develop a determining algorithm for the damper control signal that will be applied to the semi-active vehicle's vibration damping systems. Further testing will be performed, using a Honda del Sol vehicle equipped with a controlled magnetorheological damper. Figure 19 shows the suspension of the vehicle where the controlled MR dampers and the displacement sensors have been mounted to measure the suspension deflection. Further simulation studies of a vehicle model with controlled MR dampers are planned. This paper presents studies where the vehicle vibrations w ere excited by a harmonic function. In order to relate the simulation results to road conditions, tests with random function excitation will be conducted [32,33].
In some vehicles moving at high speeds, there are special aerodynamic air deflectors (spoilers) generating forces that push the wheels against the road surface [34]. As a result, wheel slippage is prevented, which was mentioned in the introduction section of this paper. In their oncoming works, the authors plan to develop a model of the vehicle equipped with controlled aerodynamic plates for the purpose of reducing changes in vertical forces acting on the wheels during transient movement, which would create greater driving safety. Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.