Mechanism of time-delay feedback control of suspension damping with an annular vibration-absorbing structure

. With the aim of enhancing both the ride comfort and the safety of the vehicle, we propose a new type of suspension with an annular vibration-absorbing structure, and establish a 3-DOF 1/4 vehicle model. The structure parameters and time-delay feedback control parameters are determined by particle swarm optimization algorithms, which take the root mean values of body acceleration, suspension dynamic deflection, and tire dynamic displacement as their optimization objectives. We analyze the stability of the suspension control system to ensure the stability of the time-delay control system through the Routh-Hurwitz stability criterion, characteristic root method, and stability switching method. Then, we compare and analyze the response characteristics of conventional suspension, new suspension without time-delay feedback control, and new suspension with time-delay feedback control under simple harmonic excitation and random excitation. The results show that the new suspension with time-delay feedback control has a significant damping effect on the body under the premise of ensuring the stability of the system.


Introduction
The vehicle is a complex multi-body system, and the vibration of the engine, the collision between different parts, and the impact of the unevenness of the road can cause the vibration of the vehicle.Severe vibration can have a negative impact on the health of passengers and vehicle components.The suspension system, as the most important vibration-damping component of the vehicle, functions to transmit the forces and moments exerted between the wheels and the body, to cushion the impact on the ground, and to reduce the vibration caused by various dynamic loads, thus ensuring that the vehicle can drive normally [1].However, suspension is an automotive component that is difficult to achieve the ideal state because it has to take into account both ride comfort and handling stability, but these two properties are contradictory to each other.Nowadays, the research on improving the overall performance of the car and reconciling the contradiction between smoothness and handling stability is mainly focused on two aspects: firstly, designing a good suspension system; secondly, finding the optimal control method.
Rabinow [2], the American National Bureau of Standards, invented the development of magnetorheological fluid materials in 1948, the appearance of new materials provided new ideas and ways for the development of the suspension.American Lord company installed magnetorheological fluid dampers into the car seat suspension system and carried out road tests, and the results showed that the developed car seat damper has a good effect in improving vehicle ride comfort.Researchers at Carnegie Mellon University, USA, applied the parallel structure of magnetorheological dampers and Stewart as a semi-active vibration isolation platform to the suspension system of an off-road vehicle and achieved the desired vibration isolation effect in a real vehicle test study [3].Hu et al. [4] studied a magnetorheological damper for a vehicle suspension system, as well as designed a switch-fuzzy hybrid control applied to the damping control of a semi-active suspension, and the simulation results showed that the vibration performance of the vehicle under impact excitation was significantly improved.In 1909, Frahm [5] installed an anti-sway tank on the ship by using the power-absorbing technique, which effectively suppressed the ship's sway, and the prototype of the power-absorbing vibration absorber came into being.Zhang et al. [6] compared the effects of suspensions with and without power absorbers on vehicle performance to conclude that suspensions with power absorbers can effectively control body and wheel vibrations in the high-frequency range.In 2001, University of Cambridge scholar Smith [7] creatively proposed the concept of inertia vessel, through the study of the characteristics of inertia vessel components, Smith found that inertia vessels can meet the vibration isolation requirements of suspensions, and then carried out the research work of applying them to vehicle suspensions, exploring a new way to improve the performance of conventional suspensions.In 2018, Li [8] proposed a three-element integrated ISD suspension, improving the cab and body vibration isolation performance effectively.However, these methods to solve vehicle vibration are using vibration damping and vibration isolation devices, without solving the problem from the suspension structure.
In the control process of a suspension system, the execution of the controller, the issuance and acquisition of signals, the response of the actuator, and the movement of the mechanical structure all produce a certain time-delay phenomenon, which is known as time-delay.Inherent time-delay is an unfavorable factor, which affects the stability of the system, deteriorates the system's performance, and even causes chaos and bifurcation.Therefore, it has received a lot of attention from scholars and has been studied in depth.
At the end of the last century, Olgac's [9], [10] research team proposed to use the amount of time-delay as an active control parameter to exploit the inherent damping property of time-delay for damping control, and by selecting appropriate time-delay feedback control parameters, the vibration of the main system under simple harmonic excitation can be attenuated to zero, thus transforming the time-delay into a favorable factor.Jalili et al. [11] used multiple identical time-delay dynamic absorbers to suppress the vibration of a multi-degree-of-freedom system, simulated and analyzed under a single simple harmonic excitation, with results demonstrating that time-delay feedback control is still desirable.Zhao et al. [12] investigated the damping performance of time-delay nonlinear dynamical vibration absorbers on the main system based on the multiscale method, and the study showed that the acceleration amplitude of the main system was reduced by 90 % with suitable feedback gain coefficients and time-delay values compared to that of the non-linear dynamical vibration absorbers without time-delay.Li et al. [13] studied the damping of LQR control and time-delay feedback control for the whole vehicle, and it was proven that the reasonable selection of time-delay feedback control parameters is more effective than LQR control for the whole vehicle damping.Zhang et al. [14] applied time-delay feedback control to automotive seat damping, obtained the critical time-delay and stability interval of the system by frequency domain scanning method, and simulated the system under simple harmonic excitation and random excitation, which demonstrated that the time-delay power absorber can greatly suppress the vibration response of the seat under either road excitation, effectively improving the ride comfort.
For these issues, this paper proposes a new type of suspension with annular vibration absorbing structure and applies the time-delay feedback control to the new type of suspension, using particle swarm optimization algorithm (PSO) to obtain the parameters of the annular vibration absorbing structure and the time-delay feedback control parameters under different road excitations, furthermore, applying the Routh-Hurwitz stability criterion, the characteristic root method and the stability switching method to obtain the stability interval of the system, as well as perform numerical validation.Finally, we establish the models of conventional suspension (CS), new suspension without time-delay feedback control (PVASS), and new suspension with time-delay feedback control (TDVASS), to compare and analyze the effect of time-delay feedback control on system damping using the above models.

The mechanical model
In order to research the effect of the new suspension with the annular vibration-absorbing structure on the body vertical vibration better intuitively, the structure schematic shown in Fig. 1 is constructed, in which the annular vibration-absorbing mass block is connected to the wheel hub and the swing arm support bracket through the spring, damping, and mass block support bracket.For the problem to be researched in this paper, the structure is simplified to a 3-DOF system vibration model including the body, the annular vibration-absorbing structure, and the wheel, as shown in Fig. 2(a).The tire displacement time-delay feedback control is introduced to this model to form the main research object of this paper, as shown in Fig. 2(b), and is analyzed in comparison with the conventional passive suspension, as shown in Fig. 2(c).The physical significance and values of the parameters in the model are shown in Table 1, and due to the small damping of the tires, all the researches in this paper do not consider the tire damping.
where,    −  denotes the tire displacement time-delay feedback control force,  is the time-delay feedback gain coefficient, and  is the time-delay amount.
According to the vehicle suspension performance index and passenger ride comfort performance index, for this purpose, selected body displacement  , body velocity  , annular vibration-absorbing structure displacement  , annular vibration-absorbing structure velocity  , tire displacement  , tire velocity  as state space variables, the differential equation of motion of the system Eq.( 1) can be written in the form of the following state space equation: where

Suspension inherent characteristics analysis
In order to analyze the vibration characteristics of the new suspension system controlled by tire displacement time-delay feedback in the frequency domain, the Fourier transform of Eq. ( 1) is performed separately to obtain:  .Thus, the frequency response function of the body droop displacement for the pavement excitation displacement input and the amplitude function of the body droop acceleration for the pavement excitation displacement input can be obtained as:

Annular vibration-absorbing structure parameters determination
The design of new suspension parameters with annular vibration-absorbing structure includes  ,  and  .The three parameters have a direct impact on the damping performance of the body.As for the design requirements of the suspension system, one is to achieve compact structure and occupy small space size; the other is to ensure the reliability of all force and moment transmission between the body and the wheels, and to ensure sufficient strength while guaranteeing the mass of components to be small.Therefore, the mass  of the annular vibration-absorbing structure is chosen to be 0.1 times the wheel mass  , i.e.,  = 0.1 [15].The stiffness coefficient  and damping coefficient  of the annular vibration-absorbing structure were obtained using a particle swarm optimization algorithm [16].In order to prevent the optimized parameters from not conforming to the actual characteristics of the annular vibration-absorbing structure, the optimization ranges of  and  were set to (0,  ) and (0,  ), respectively.

Establishment of the objective function
In this paper, the optimization objective of the new suspension with annular vibration-absorbing structure is to suppress the body vibration to the maximum extent, improve the ride comfort, handling stability and tire grounding of the vehicle, so the RMS values of BA, SDD and TDD are taken as the optimization indexes, and the optimization objective function  is established: where,  = 0.4 ,  = 0.2 , and  = 0.4 are weighting coefficients.

Particle swarm optimization algorithm (PSO)
PSO is a classical swarm intelligence algorithm, which mainly uses the sharing of information of individuals in the group, so that the motion of the whole group can produce the evolution process from disorder to order in the problem solution space, and thus obtain the optimal solution of the problem, so the PSO is used to optimize the solution of the stiffness coefficient  and damping coefficient  of the annular vibration-absorbing structure.The specific steps of the particle swarm optimization algorithm are shown below, Step 1: Initialize the parameters of the optimization algorithm and generate  particles randomly within the parameter optimization range, giving each particle a random initial position  and an initial velocity  .
Step 2: Design of the fitness function, and in order to optimize the damping performance of the system, the fitness function is therefore designed as shown in Eq. ( 6).
Step 3: Calculate the fitness value of the generated particles and select the parameter that minimizes the fitness value as the individual best value and the objective function that minimizes the fitness value as the global best value.
Step 4: Update the position and velocity of each particle as shown in Eq. ( 7): where  is the number of current iterations,  is the inertia weight,  and  are learning factors,  and  are random values of [0,1],  is the individual optimal particle position,  is the global optimal particle position,  is the particle flight speed, and  is the current position of the particle.
Step 5: Optimize the second set of parameters, get the individual best and global best according to steps 2 and 3, and compare them with the first set, and update the individual best and the global best.
Step 6: According to the above steps, the individual optimum and the global optimum are updated continuously until the set number of iterations  is satisfied, and the system will output the optimized optimum parameters.
The flow chart of the algorithm steps is shown in Fig. 3. Step 1 Step 2 Step 4 Step 6 Step 3 Step 5 By the above particle swarm optimization algorithm, the stiffness and damping coefficients of the annular vibration-absorbing structure are optimized, and the optimal parameters of the annular vibration absorbing structure are obtained as  = 6444 N/m,  = 348 N⋅s/m.Thus, the vehicle model parameters are shown in Table 1.

Stability analysis
The selection of time-delay feedback control parameters is decisive for the system to be in a stable state.If the time-delay feedback control parameters are not selected properly, chaos and bifurcation will occur, and in serious cases, the system will be destabilized.In order to ensure that the system is in a stable state, this paper uses the Routh-Hurwitz stability criterion, the characteristic root method and the stability switching method to analyze the stability of the system.
By applying Laplace transform to Eq. ( 1), the differential equation of motion of the system is deformed into the characteristic equation: where,  is the characteristic root of the equation, () and () are the coefficients of the polynomial, respectively: which:

Time-delay independent stability
Whether the time-delay feedback control can ensure that the system is stable depends on the value of the feedback gain  and the time-delay amount .When any value of the time-delay amount is taken, the stability of the system is not changed with the time-delay amount.The stability of the system does not change with the time-delay amount when the time-delay amount takes any value, but only changes with the feedback gain, which is called the full time-delay stable state, i.e., time-delay independent stability.For the time-delay independent stability of the system, it is necessary to satisfy the following two conditions, (1) When  = 0, the characteristic Eq. (8) of the system satisfies Routh-Hurwitz stability.
(2) When  > 0, the characteristic Eq. (8) of the system has no positive real roots  .Condition (1): when  = 0, the characteristic equation of the system polynomial Eq. ( 8) is expressed as: It is concluded from the Routh-Hurwitz stability criterion that for the system to be in a stable state, firstly, all coefficients of the characteristic Eq. (10) of the system are required to be greater than zero; secondly, all terms in the first column of the Routh-Hurwitz table are required to be positive, and the Routh-Hurwitz discriminant table constructed from the polynomial coefficients can be obtained according to Eq. ( 10), as shown in Fig. 4, which, Substituting the parameters of Table 1 into Fig. 4, the stability conditions  ,  ,  ,  ,  ,  , and  can be found to be constant greater than zero.Other stability regions drawn according to the Routh-Hurwitz stability conditions are shown in Fig. 5. From Fig. 5, it is clear that when  ∈ (−∞, 171870.8], the polynomial coefficients of the characteristic equation of the system and the first column coefficients in the Routh-Hurwitz discriminant table are greater than zero, which satisfies the Routh-Hurwitz stability criterion and indicates that the no time-delay system is stable in this interval. Condition (2): When  > 0, it is known from stability theory that the asymptotic stability of the time-delay system needs to satisfy that the characteristic roots of the polynomial of the characteristic equation don't contain positive real roots.When the characteristic roots contain pure imaginary numbers, the system is at the boundary of instability and there exists stability switching as well as critical time-delay stability interval.Assuming the existence of imaginary roots  =  in the characteristic equation, the characteristic polynomial of the system Eq.( 8) can be expressed as (, ) = ( , ) = 0 and simplified by using Euler's formula  = cos( ) − isin( ) for ( , ) = 0, which can be obtained: where: Separating the real and imaginary parts shows that if the polynomial Eq. ( 11) is zero, it must satisfy that both the real and imaginary parts are zero, i.e.: Further derivation leads to: With sin ( ) + cos ( ) = 1, the trigonometric term is eliminated to obtain a polynomial equation for  , i.e.: Such that ( ) =  ( ) +  ( ) − [ ( ) +  ( ) ], the characteristic root discriminant polynomial of the system can be obtained, i.e.: where: After a series of derivations, ( , ) = 0 can finally be reduced to ( ) = 0.By solving the value and number of positive real roots of the polynomial equation ( ) = 0, the range of its feedback gain  can be determined when the time-delay independent stability is achieved.Fig. 6 shows the time-delay independent stability region of the new suspension under the tire displacement time-delay feedback control.
According to Fig. 6, when  ∈ [−24427.7,24427.7],the characteristic equation of the system has no positive real roots, and the range of values of this gain coefficient  is in the stability interval of condition Eq. ( 1), which means that the system is time-delay independent stable in this interval range.When the value range of this feedback gain is exceeded, the polynomial equation has positive real roots, which means that the stability of the system will switch a certain number of times and needs further discussion and analysis.

Stability switching
When the time-delay feedback control parameter  is outside the region of time-delay independent stability, the stability of the system is not fixed, and the sign of the real part of the characteristic root changes with time-delay, which means that the stability of the system changes with time-delay, and this situation is called stability switching.Therefore, finding the critical time-delay is crucial to analyze whether the system is stable or not.
Assuming that the polynomial Eq. ( 15) has positive real roots  , all the critical time-delays corresponding to the positive real roots  can be solved by the following equation: where,  =  ,  ∈ [0,2].After all the critical time-delays are determined by Eq. ( 16), the change of the real part of the characteristic root corresponding to each critical time-delay is then judged to further determine whether the system is stable or not.The characteristic equation polynomial Eq. ( 8) can be regarded as a function of  with respect to .Its derivative can be obtained from the stability switching discriminant formula: Then, the change of the real part of the characteristic root can be determined according to the positive or negative sign of .In the case of  = +1, as the time-delay  increases, it means that the real part of the characteristic root crosses the imaginary axis from left to right in the complex plane, thus switching the system from stable to unstable state; in the case of  = −1, it means that the real part of the characteristic root crosses the imaginary axis from right to left in the complex plane, and the system switches from unstable to stable state.The following is a concrete example to further illustrate the switching principle of stability.
When  = 50000 Ν/m, the number of positive real roots of the characteristic equation is two,  * = 62.3487,  * = 75.1147,according to Eq. ( 15), as shown in Fig. 6.Substituting these two positive real roots into Eq.( 17), we can get  ( * ) < 0 and  ( * ) > 0, and thus the critical timedelay corresponding to these two positive real roots can be found according to Eq. ( 16) as:  * = 0.0501,  * = 0.1338,  * = 0.2174,  * = 0.3011,  * = 0.3875, ⋯ ,  * = 0.0945,  * = 0.1952,  * = 0.2960,  * = 0.3968,  * = 0.4976, ⋯ . ( The obtained critical time-delays are sorted to obtain: Since  ( ) < 0, i.e.,  = −1, when the time-delay  crosses the time-delay  , * ( = 1,2,3, ⋯ ), the real part of the characteristic root changes from positive to negative, and the system switches from unstable to stable state; similarly, when  ( ) > 0, i.e.,  = +1, the timedelay  crosses  , * ( = 1,2,3, ⋯ ), and the real part of the characteristic root changes from negative to positive, then the system switches from stable to unstable state.From the above analysis, it can be seen that when  ∈ [0,  * ), the system is asymptotically stable; when  ∈ [ * ,  * ], the system is unstable.From Eq. ( 19), it can be seen that the critical time-delay  * is followed by  * , no switching occurs, for two consecutive unstable critical time-delays, then it shows that after  * , the system is always in an unstable state.The final time-delay interval that puts the system in a stable state when When  = -50000 Ν/m, the positive real roots of the characteristic equation are obtained as  * = 62.3487,  * = 75.1147,respectively, as shown in Fig. 6.Further solving for  ( ) > 0 and  ( ) < 0, then the critical time-delays corresponding to these two positive real roots are obtained as:  * = 0.0083,  * = 0.0919,  * = 0.1756,  * = 0.2592,  * = 0.3429, ⋯ ,  * = 0.0441,  * = 0.1449,  * = 0.2456,  * = 0.3464,  * = 0.4472, ⋯ . (20) The obtained critical time-delays are sorted to obtain: It can be concluded that no switching occurs after the critical time-delay  * , which directly becomes  * , for two consecutive unstable critical time-delays, then it is derived that after  * the system are in an unstable state.Therefore, when  = -50000 N/m, the critical time-delay interval that keeps the system stable is Fig. 7 shows the critical time-delay curves of (, ) under tire displacement time-delay feedback control.The intersection of the purple and light blue curves on the right side of Fig. 7 is before the critical time-delay point, and the system is stable, while the intersection of the green and purple curves is after the critical time-delay point, and the system is unstable, and the system state has changed, i.e., the stability switching occurs between [ * ,  * ] shown in the figure .Similarly, the intersection of the purple and light blue curves on the left side of Fig. 7 is on the outside of the critical time-delay point, and the system is stable, while the intersection of the green and purple curves is on the inside of the critical time-delay point, and the system is unstable, which indicates that the system has switched stability between [ * ,  * ].In the case that the gain coefficient is determined, the polynomial solution of the characteristic equation is used to get the positive real roots.If there are no positive real roots in the equation, the system is stable regardless of the value of ; if there are positive real roots in the equation, the system will undergo a certain stability transformation and finally become unstable.According to the change of the number of positive real roots, we can get the full time-delay stability region and stability switching region of the system, and then judge whether the system corresponding to the critical time-delay point is in stable state according to the stability switching discriminant formula, and finally we can get the critical stability region in the (, ) plane under the tire displacement time-delay feedback control, as shown in Fig. 8.

Stability verification
In order to verify the correction of the region in Fig. 8, three points are randomly selected for verification in each stability interval and instability interval of the system, where points  * ( = 40000,  = 0.03),  * ( = -40000,  = 0.25), and  * ( = 40000,  = 0.3) correspond to the stability interval of the system; points  * ( = -40000,  = 0.11),  * ( = 40000,  = 0.15), and  * ( = -40000,  = 0.4) correspond to the instability interval of the system.The coordinates of each point are substituted into Eq.( 1), and the initial condition is taken as From the numerical verification results in Fig. 9, it can be seen that when the initial excitation of the system is given, the BA vibration amplitude gradually converges to zero, which indicates that the system is asymptotically stable; the BA vibration amplitude gradually diverges to infinity, which indicates that the system is unstable.The verification results in Fig. 9 are consistent with the corresponding system stability cases in Fig. 8, and the numerical simulation results verify the correctness of the critical stability time-delay stability region map.

Method of solving the time-delay equation
The solution is difficult because the characteristic equation contains a transcendental function term.Therefore, in this paper, the fine integration method of time-delay differential equations is used to solve the numerical solution of each response quantity of the 1/4 vehicle model vibration.
At first, treating the time-delay feedback control term () in Eq. ( 2) as the excitation term, the solution of the equation can be written as: Next, the numerical discretization of Eq. ( 22) is performed by assuming that the simulation time is .Taking the time step Δ =  −  ( ∈ [0, ]), then  = /Δ, the solution of the equation is discretized into the following stepwise integral equation in recursive form: The above equation can be further derived as: where,  = /Δ, is expressed as the number of nodal steps that the time-delay quantity  is divided into during the solution process.Further solving the solution of Eq. ( 24), assuming the initial condition () = [0 0 0 0 0 0] of the equation, the system vibration response of each mass block (body, annular vibration-absorbing structure, wheels) at each time point  is found as: where,  ,  represents the vibration displacement response and vibration velocity response of the body mass  at time  ;  ,  represents the vibration displacement response and vibration velocity response of the annular vibration-absorbing structure mass  at time  ;  ,  represents the vibration displacement response and vibration velocity response of the tire mass  at time  , respectively.Finally, by substituting Eq. (25) into Eq.( 1), the vibration acceleration of the body, the annular vibration-absorbing structure and the tire in the time domain are obtained respectively: (26)

Optimization of feedback control parameters
For the optimization of the time-delay feedback control parameters the PSO described above is still used, and its objective function is shown in Eq. ( 6).With 100 randomly generated particles and 200 iterations, the fine integration method for solving the time-delay and the pavement displacement are introduced into the PSO to optimize the objective function in the time domain and solve the appropriate time-delay feedback control parameters for different pavements.The simple harmonic excitation and random excitation are selected for the pavement displacement, and the variation curves of the adaptation degree values with the number of iterations for different pavement excitations are shown in Fig. 10.From Fig. 10, it can be seen that the optimal time-delay feedback control parameters are basically obtained at the later stage of the iteration number, and the optimal time-delay feedback control parameters under simple harmonic excitation and random excitation are finally obtained as  = 5386 N/m,  = 2×10 -4 s and  = 2528 N/m,  = 1.28×10 -4 s, respectively.The optimized two sets of parameters are in the stability region in Fig. 8, which can make the system in a stable state, indicating that the optimized time-delay feedback control parameters are reasonable.

Simulation analysis under simple harmonic excitation
The above optimized tire displacement time-delay feedback control parameters  and  are substituted into Eq.( 1), and the simple harmonic excitation  = 0.05sin(5) is simulated as the tire road input.The effects of the CS, the PVASS, the TDVASS on the system damping performance in the time domain are analyzed.The simulation curves of each response of the system were obtained as shown in Fig. 11, and the RMS values of each response quantity of the system were shown in Table 2.
As can be seen from Fig. 11, under the optimal feedback control parameters, the BA, SDD and TDD of the TDVASS are reduced compared with the PVASS, and the amplitude is lower in the full frequency domain than that of the PVASS; while compared with the CS, the BA and TDD are significantly reduced, and the SDD compared with the CS, the BA and TDD are significantly reduced, while the SDD is slightly deteriorated, but the increase is still within the design requirement (the suspension dynamic travel of ±100 mm is selected in this paper [17]), which doesn't exceed the limit travel of the suspension dynamic deflection, i.e., it doesn't increase the probability of hitting the limit block.In Table 2, it can be obtained that compared with the CS, the RMS values of BA and TDD of the TDVASS are reduced from 1.6081, 3.1109 to 0.3068, 0.6251, respectively, and the optimization ratios are 80.92 % and 79.91 %; compared with the PVASS, the RMS values of BA, SDD, TDD of the TDVASS are reduced from 0.9715, 0.0693, 1.6606 to 0.3068, 0.0388, 0.6251, respectively, and the optimized ratios are 68.42 %, 44.01 %, 62.36 %, respectively.The above analysis shows that the tire displacement time-delay damping control method can effectively alleviate the impact of road excitation on the body vibration, so that its comfort can be significantly improved.At the same time, the grounding of the tires is also improved, which guarantees the safety of the vehicle.

Simulation analysis under random excitation
The simple harmonic excitation is only for a certain fixed frequency pavement, which is more desirable for a road, while the random excitation is more in line with the actual road conditions.According to the filtered white noise method described in ISO 8608:1995 as shown in Eq. ( 27), a standard C-level random road surface [18] is constructed as shown in Fig. 12(a), which is used to further verify the damping effect of the TDVASS on the body.To illustrate the effectiveness of the time domain random pavement excitation model generated based on the filtered white noise method, a comparative power spectral density plot of the C-level pavement is established, as shown in Fig. 12(b), which shows that the simulated power spectral density of the pavement unevenness is in good agreement with the standard power spectral density, indicating the correctness of the generated pavement.This pavement can be used for simulation analysis of the model:    According to Fig. 13 and Table 3, it can be concluded that the TDVASS has greater advantages compared to the CS and the PVASS.Compared with the CS, the RMS values of BA and TDD are reduced from 0.3184, 0.6399 to 0.1014, 0.5289, respectively, with the optimized ratio of 81.72 % and 19.41 %, and the SDD is slightly deteriorated, but its vibration amplitude is still within the dynamic travel of the suspension and does not exceed its limit travel.Compared with the PVASS, the RMS values of BA, SDD, and TDD were reduced from 0.1512, 0.0105, and 0.5607 to 0.1014, 0.0073, and 0.5289, respectively, with optimization ratios of 61.51 %, 28.57%, and 8.03 %.As a result, it can be derived that the time-delay feedback control can also effectively reduce the body vibration and improve the comfort of the vehicle and the grounding of the tires under the random excitation with the selection of suitable feedback parameters.

Amplitude and frequency characteristics analysis
Through frequency domain analysis of CS, PVASS and TDVASS, the response quantities of each vibration characteristic of the suspension system at different frequencies in the studied frequency band can be obtained, and the effect of different time-delay feedback control on vehicle damping can be further investigated.With Eq. ( 4) and ( 5), the amplitude frequency characteristic functions of the system can be obtained.The optimal time-delay feedback control parameters  and  obtained under random road excitation are substituted into the system to obtain the amplitude frequency characteristic curves of BD and BA under tire displacement time-delay feedback, as shown in Fig. 14  It can be seen that due to the addition of the under-spring mass annular vibration-absorbing structure to form the PVASS on the basis of the CS , and the introduction of the tire displacement time-delay feedback control to form the PVASS, which results in the low frequency resonance frequencies of the PVASS and the TDVASS being smaller than those of the CS, and the high frequency resonance frequencies are almost unchanged.It can also be concluded that compared with the CS, the TDVASS only has a certain degree of deterioration in the low frequency resonance region, and the optimization effects in the rest of the frequency bands are substantially improved.With PVASS, the amplitude frequency characteristic curves of BD and BA in TDVASS show substantial improvement in both the low frequency region and the medium frequency region, and the curves almost overlap in the high frequency region with no obvious optimization effect.The above analysis shows that TDVASS has obvious effect on reducing body vibration in the low and medium frequency bands, and the optimization effect is not obvious in the high frequency band.

Conclusions
In this paper, a new suspension with an annular vibration absorbing structure is designed to reduce body vibration, and a time-delay feedback control strategy based on tire displacement is introduced to further suppress body droop vibration and improve comfort.The main conclusions are as follows: 1) Designed a new suspension with an annular vibration-absorbing structure, and derived the parameters of the annular vibration-absorbing structure by PSO.And the time-delay feedback control of the tire displacement is introduced, and the parameter values of the time-delay feedback control determine whether the system is stable or not.
2) The time-delay independent stability region of the system is obtained by using the characteristic root method, in which the critical time-delay stability region is obtained by using the stability switching method, meanwhile, the correctness of the stability switching method is verified by numerical simulation.
3) The vibration differential equations of the time-delay system are solved by the fine integration method, and the time-delay feedback control parameters are optimized by introducing the fine integration method and different road excitations into the PSO, and finally the optimal time-delay feedback control parameters for different road surfaces are obtained, so that the system can obtain the best vibration damping control effect.
4) Simulation analysis of the optimal feedback control parameters for different road surface excitation in time domain and frequency domain shows that the tire displacement time-delay feedback control has a significant effect on the body vibration damping.It shows that the timedelay feedback control can significantly improve the ride comfort and tire grounding of the vehicle while ensuring the driving safety.

10 .
a) Simple harmonic excitation b) Random excitation Fig. Adaptation function change curves a) Body droop acceleration b) Suspension dynamic deflection c) Tire dynamic displacement Fig. 11.Simulation comparison under simple harmonic excitation where,  is the under cut-off frequency, generally take  = 0.1;  ( ) is the A-H level pavement unevenness coefficient;  is the reference spatial frequency,  = 0.1 m -1 ;  is the car driving speed, this paper takes  = 20 m/s; () is the Gaussian white noise signal with zero mean value.a) C-level pavement unevenness b) Power spectrum density comparison chart Fig. 12. C-level pavementUsing Fig.12(a) as the random excitation of the pavement and substituting the optimal feedback control parameters  and  into Eq.(1), the simulation curves of each response quantity of the system are shown in Fig.13, and the root mean square values of each response quantity of the system are shown in Table3 .

14 .
a) Body droop displacement b) Body droop acceleration Fig. Amplitude frequency characteristic curves

Yang
Nan studied in Shandong University of Technology, School of transportation and vehicle engineering, Master of Science in Mechanical.His current main research interests are energy analysis, vibration analysis, and vehicle vibration reduction.Pengcheng Zhou studied in Shandong University of Technology, School of transportation and vehicle engineering, Master of Science in Mechanical.His current main research direction is four-wheel steering.

Table 1 .
Vehicle model parameters meaning and values

Table 2 .
RMS values of vehicle performance indicators under simple harmonic excitation

Table 3 .
RMS values of vehicle performance indicators under random excitation