Next Article in Journal
Design and Optimization of Low Impact Shift Control Strategy for Aviation Transmission Power System Based on Response Surface Methodology
Previous Article in Journal
Sampling-Based Path Planning Algorithm for a Plug & Produce Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Propagation Characteristics of Tire Cavity Resonance Noise in the Automobile Suspension

1
School of Mechanical, Electronic and Control Engineering, Beijing Jiaotong University, Beijing 100044, China
2
School of Transportation Science and Engineering, Beihang University, Beijing 100191, China
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2023, 13(22), 12113; https://doi.org/10.3390/app132212113
Submission received: 11 August 2023 / Revised: 31 October 2023 / Accepted: 2 November 2023 / Published: 7 November 2023
(This article belongs to the Section Acoustics and Vibrations)

Abstract

:

Featured Application

The research of propagation characteristics can provide a significance and value for McPherson suspension system design.

Abstract

Tire cavity resonance noise (TCRN) is one of main contributors to vehicle interior noise, which has long been a concern in the automotive industry and academia. As suspension is crucial for the propagation of TCRN energy into the vehicle, the propagation characteristics of energy in automobile suspension is studied in this research. Firstly, the finite element model of a McPherson suspension system connected to an aluminum alloy wheel with a Y-shaped spoke was established. Then, a modal analysis and response calculations of the McPherson suspension system connected to the aluminum alloy wheel with a Y-shaped spoke were carried out. Finally, the propagation characteristics of TCRN in the McPherson suspension system connected to the aluminum alloy wheel with a Y-shaped spoke were studied and analyzed by the power flow method under different working conditions. The power flow output via the lower arm front bushing was the largest, while the output via the rear bushing was the smallest in the Y-spoke aluminum alloy wheel and suspension system. The areas in the suspension system with high stress are located at the steering knuckle, lower swing arm, and shock absorber. Therefore, study of the propagation characteristics can provide a basis for a McPherson suspension system design.

1. Introduction

With growing requirements for comfort in vehicles, TCRN has attracted the attention of researchers and has become an important area of tire–noise research around the world in recent years [1,2]. Although many results related to research on TCRN have been published, research of the propagation characteristics is still ongoing, and there is still no simple and easier method to suppress TCRN which could guide the design of automobiles [3,4].
The energy of TCRN propagates in the wheel and suspension system; an appropriate method must be chosen to reduce the propagation of energy in the wheel and suspension system. The power flow method is a widely applicable method to clarify the propagation characteristics of an elastic wave. In 1962, the power flow method was first put forward [5]. It was used to study the energy transfer to the seat through a vibration isolator, and a method was proposed to calculate the power flow of a structure [6]. The power flow method was used to compute the energy transfer of a thin-walled structure and to calculate power flow in simple two-dimensional structures [7]. Related research on surface energy and surface effects has proven to be helpful in power flow research [8,9,10]. The method of combining finite element analysis and vibration power flow calculations was used to obtain complex elastic foundation modal information using finite element software and to substitute it into the elastic foundation admittance [11]. The power flow of the structural sound intensity method was introduced to study the propagation law of vibration energy in a thin-walled rectangular beam, and the difference in power flow at response points before and after applying a concentrated mass was examined [12]. Based on our power flow analysis, a power flow vibration and shock-fusion controller was proposed. This fusion controller covered the vibration control module that attenuated vibration excitation and the shock control module that buffered shock excitation and vibration [13]. The theory of vibration power flow was used to optimize the selection of the damping element bushing. The dynamic response of the system before and after parameter optimization and the influence of the bushing parameters on the power transmission were discussed in ref. [14].
The above works show that the power flow method has wide and mature applications in many fields, but so far, it has not been applied to research on the propagation characteristics of TCRN in a McPherson suspension system connected to an aluminum alloy wheel. Therefore, we applied the power flow method to research the propagation characteristics of TCRN in a McPherson suspension system connected to an aluminum alloy wheel. The laws of propagation characteristics of TCRN in a McPherson suspension system under different working conditions (loads, inflation pressures, or speeds) can provide the necessary basis for the design of a McPherson suspension system.
To this end, a finite element model of a McPherson suspension system connected to an aluminum alloy wheel was established, and modal analyses and response calculations of the McPherson suspension system were carried out. Finally, the power flow method was used to research the propagation characteristics of TCRN in a McPherson suspension system connected to an aluminum alloy wheel with a Y-shaped spoke, and the propagation characteristics of TCRN in the McPherson suspension system under different working conditions were obtained.

2. An Introduction to the Power Flow Method

To research the propagation of vibration energy in a structure, the power flow method is widely used. With it, the transmission path of the vibration energy in a structure can be identified [15]. Vibration power flow describes the energy consumption of structural vibrations or the ability of an external force to do work per unit time. The average vibration power flow in a period is generally used as a reference, and input power flow p in a structure is determined as follows:
p t = F t V ( t )
where F(t) is the external force acting somewhere in the structure and V(t) is the velocity at that location produced by F(t).
Average vibration power flow P corresponding to a certain period T is calculated as follows:
P = 1 T 0 T F t V t d t
To calculate the power flow in a structure based on the finite element method, the nodal force and velocity vectors need to be calculated. The power flow per unit area of a certain section of the structure is its structural sound intensity, so stress can be introduced into the calculation. For an elastic body, the stress state at any point inside it can be expressed with six stress components: shear stress, i.e., τ x y , τ y z , and τ z x , and normal stress, i.e., σ x , σ y , and σ z . As an example, the stress distribution of a micro-element is shown in Figure 1.
The formula of structure sound intensity p n is as follows:
p n = 1 2 R e τ n 1 v 1 + τ n 2 v 2 + σ n v n
where τ n 1 and τ n 2 are the shear stress in the 1-direction and 2-direction, respectively, σ n is the normal stress in the normal n-direction, and v 1 , v 2 , and v n are the velocity complex conjugates in the 1-direction, 2-direction, and n-direction, respectively.
It can be seen from the literature [16] that stress and displacement can be evaluated by calculating the harmonic response results. A solid element has three directions, x , y , and z , so the structural sound intensity in each direction is expressed as:
p x = ω 2 I m τ x y u y + τ x z u z + σ x u x
p y = ω 2 I m τ y x u x + τ y z u z + σ y u y
p z = ω 2 I m τ z x u x + τ z y u y + σ z u z
where ω is the angular velocity, τ x y , τ x z , τ y x , τ y z , τ z x , and τ z y are the shear stress in the x , y , and z directions, respectively, u x , u y , and u z are the displacement conjugates of the node in the x , y , and z directions, respectively, and p x , p y , and p z are the structural sound intensity in the x , y , and z directions, respectively.

3. Establishment and Response Calculation of a McPherson Suspension System Model

3.1. Estabalishment of a McPherson Suspension System Connected to an Aluminum Alloy Wheel

The McPherson suspension system has a simple structure, low cost, and good performance, so it is widely used in the front suspension of automobiles. It consists of a coil spring, knuckle assembly, lower transverse arm, shock absorber, etc. With reasonable simplification, a McPherson suspension system model is shown in Figure 2a. With a connected 15-inch Y-spoke aluminum alloy wheel, the system model is shown in Figure 2b.
According to a tire cavity sound pressure distribution under the condition of 4000 N load, 0.22 MPa inflation pressure, and 60 km/h speed, the excitation acting on the rim is as shown in Figure 3, as described in previous literature [17,18]. The resonant sound pressure distribution function of the tire cavity is shown in Equation (7).
P θ , t = ( 80 30 cos θ + 30 sin θ + 32 cos 2 θ 24 sin 2 θ + 3 cos 3 θ + 14 sin 3 θ ) cos ( 2 π f t )
where θ is the included angle between the radius of a certain point of the wheel rim and the radius of the contact position between the tire and the drum, and f is the tire cavity resonance frequency.
For our simulation, some parts of the McPherson suspension system were simplified. For example, the outer end of the lower transverse arm was in bound contact with the knuckle assembly, and the inner end of the lower transverse arm was connected to the subframe through a bushing. The upper end of the shock absorber was fixed, and a spring unit was used to simulate the stiffness of the coil spring and the damping of the shock absorber. The spring unit was connected to the knuckle assembly through a cylindrical joint and to the upper end of the shock absorber at the other end. The knuckle assembly was in bound contact with the hub, and there was also bound contact between the hub and the aluminum alloy wheel. The stiffness and damping [19,20,21,22,23] of the associated bushings and spring units are shown in Table 1, and the simplified McPherson suspension system and transverse arm bushings are shown in Figure 4. A point was selected at the center of the bushing, and the bushing was linked to that point, with the connection between this point and the bushing representing the stiffness and damping. The input surface of the TCRN is the wheel rim surface. The output surface is the inner surface of the cylindrical joint support and the front and rear bushing surfaces of the lower transverse arm.

3.2. Modal Analysis and Response Calculation of a McPherson Suspension System

As a key link of TCRN propagation to the vehicle, the natural frequency of a suspension system is close to the tire cavity resonance frequency, which directly affects the propagation of this noise energy. We performed a modal analysis of a suspension system model using the Lanczos method in Abaqus 6.14. We focused only on a low-order modal analysis in the vast majority of cases and generally selected the first six orders of modes. Then, the first six orders of natural frequencies and vibration modes of the suspension system were taken, as shown in Figure 5. The finite element type was C3D10HS, which is a 10-node, general purpose quadratic tetrahedron with improved surface stress visualization. The DOF of C3D10HS is ten, and the behavior of the element type is in the form of a Lagrange description. The number of elements was 393,719.
It can be concluded that the steering knuckle bends in the longitudinal plane from the third order to the fourth order, and the lower swing arm bends in the longitudinal plane from the first to the sixth order. The natural frequencies of the suspension system are shown in Table 2. It can be seen that the sixth natural frequency was close to the resonance frequency of the tire cavity [17]. So, the FE suspension results were verified.
The Y-spoke aluminum alloy wheel and suspension system were simulated and analyzed. The suspension was made of steel and the type of elastic was isotropic. Our analysis was linear, using the implicit direct method of integration, where the steady state response of the system was calculated by directly integrating the original equations of the model. The dynamic response of tire cavity resonance noise was calculated, and the propagation characteristics of tire cavity resonance noise in suspension were obtained, providing a basis for the suppression of this noise. In our direct steady-state dynamics analysis of the suspension system, the response at the first frequency of tire cavity resonance was relevant. The defined sweep frequency range was 210–260 Hz, and 101 frequency points were defined, with one frequency excitation point every 0.5 Hz. Using this method, the stress and strain of the suspension system at the first natural frequency of the tire cavity resonance noise (238 Hz) were evaluated. The frequency (238 Hz) came from an experiment described in the literature [17], i.e., this value was not the result of other simulations. The stress and strain contours of the suspension system, obtained in the frequency domain, are shown in Figure 6. The unit of stress is MPa.
It can be seen from Figure 6 that the stress and strain of the suspension system were relatively large at the lower swing arm, steering knuckle, and shock absorber. Therefore, it can be speculated that the energy of the tire cavity resonance noise propagating to this area would also be large.

4. Power Flow Analysis of Tire Cavity Resonance Noise Propagation in a Suspension System

According to the steps of calculating the power flow in the wheel [18], the Y-spoke aluminum alloy wheel and suspension system were simulated. The distribution of structural sound intensity in the system under the conditions of 4000 N load, 0.22 MPa inflation pressure, and 60 km/h speed could be obtained, as shown in Figure 7. It can be seen that there was a large degree of structural sound intensity at the local spokes.
After the suspension system with a wheel had been studied, the real and imaginary values of stress and displacement at each point, the excitation frequency of the system response, and the point set data of the input and output surfaces were derived. The above txt format data were read by MATLAB R2014a, and the structural stress at each point in the wheel and the structural sound intensity of the input and output surfaces could be calculated according to the structural sound intensity formula. After obtaining the structural sound intensity and power flow, a vector diagram could be obtained by further programming. The contour of the structural sound intensity could be imported into the simulation software through a program compiled in the Python language to realize the secondary development. Then, the contour of the structural sound intensity could be obtained. According to the steps described in ref. [18], the structural sound intensity vector distribution figures and a partial enlargement system in different planes could be obtained, as shown in Figure 8.
As shown in Figure 8, the structural sound intensity distribution of the suspension system was uneven. In the partial enlargement of Figure 8, the length of the arrows represents the amplitude of structural sound intensity, and the direction of the arrows represents the direction of structural sound intensity. During the transmission of structural sound intensity vectors, structural sound intensities of different directions and sizes meet and merge into a whole vector toward the hub center. Structural sound intensity vectors with the same frequency and different phases propagate from the spokes to the hub and continuously converge and dissipate.
For the Y-spoke aluminum alloy wheel and suspension system, data such as the sum of structural sound intensity and cross-sectional area at each point of the wheel input surface and output surface could be calculated. The power flow could be determined by formulas presented in ref. [18], and the results are shown in Table 3. Power flow is the product of the average structural sound intensity at the node and the corresponding cross-sectional area, as shown in Equation (8).
P f = p s i N S
where Pf is the power flow of cross-section S, psi is the structural sound intensity of point i, S is the area of the corresponding section, and N is the number of points.
The first output, second output, and third output in the table are the inner surface of the cylindrical auxiliary support and the front and rear bushing surfaces of the lower swing arm, respectively.
It can be seen from Table 3 that for the suspension system under 4000 N load, the output power flow of the front bushing of the lower swing arm was the largest and that of the rear bushing was the smallest. Therefore, special attention should be paid to the front bushing of the lower swing arm in order to suppress the transmission of the resonant noise energy of the tire cavity into the vehicle.
For different vehicle speeds and load conditions, the sound pressure amplitude function of the tire cavity resonance noise experienced by the wheels is shown in Table 4.
According to the amplitude distribution of resonant sound pressure in the tire cavity under different load conditions at 0.22 MPa inflation pressure and 60 km/h speed, as shown in Table 4 (the meanings of symbols in the function are the same as before), the model of a Y-spoke aluminum alloy wheel and suspension system after loading is shown in Figure 9.
The above five suspension system models under different loads were analyzed using the finite element method, and the response of an aluminum alloy wheel under the action of tire cavity resonance noise was calculated. The stress contour of the aluminum alloy wheel and suspension system under different loads is shown in Figure 10, and the corresponding strain contour is shown in Figure 11. The unit of stress is MPa.
It can be seen from the stress and strain contour of the suspension system under the above five different load conditions that the areas with high stress are located at the steering knuckle, lower swing arm, and shock absorber. The resonance noise energy of the tire cavity acted behind the rim, spread to the inside of the wheel, and finally spread to the steering knuckle, lower swing arm, and shock absorber through the spokes, resulting in high stress and strain in these three places. Therefore, it can be concluded that the energy transmitted to this area was concentrated, and the structural sound intensity was correspondingly large. With an increase of load, the stress and strain of our suspension system increased.
Similar to the steps in the previous section, the structural sound intensity of a suspension system under different loads could be evaluated, as shown in Figure 12. It can be seen that there was a large degree of structural sound intensity near the rim edge and hub. With an increase of load, the structural sound intensity of suspension system increased.
For the suspension system under different loads, we calculated data such as the sum of the structural sound intensity and cross-sectional area at each point of the wheel input surface and three output surfaces. The power flow was obtained as shown in Table 5.
It can be seen from Table 5 that the total power flow of the three output surfaces of the suspension system under different loads was quite different, and different loads had obvious influences on the transmission of power flow. Additionally, the output power flow of the suspension system decreased slightly with an increase in load. It was this increase in load that led to an increase of input energy. In the process of energy transmission, damping consumes a great deal of energy, but the output energy is lower.
Based on the sound pressure amplitude distribution of the tire cavity resonance noise under different inflation pressure conditions under a load of 3500 N and at a speed of 50 km/h, as shown in Table 4, the power flow results of the input surface and three output surfaces of the suspension system were calculated, as shown in Table 6.
From Table 6, it can be seen that the difference of the total power flow of the three output surfaces of the suspension system under different inflation pressures was small. Additionally, the influence of different inflation pressures on the transmission of power flow was not significant. The inflation pressure was 0.22 MPa, and the output power flow of the suspension system was the largest. Under the other four inflation pressures, the output power flow of the suspension system was relatively low.
Based on the sound pressure amplitude distribution of the tire cavity resonance noise under different speed conditions with 0.22 MPa inflation pressure and 3000 N load, the power flow results of the input surface and three output surfaces of the suspension system were calculated, as shown in Table 7.
From Table 7, it can be seen that for the suspension system at different speeds, the total power flow of the three output surfaces was quite different. Additionally, different speeds had an obvious influence on the transmission of power flow. As the speed increased, the output power flow of the suspension system decreased slightly.

5. Conclusions

In this research, the finite element model of a McPherson suspension system was established and the correctness of this model was verified by modal analysis. Finally, the power flow method was used to research the propagation characteristics of TCRN in the suspension system, and the propagation characteristics of TCRN at different running speeds, inflation pressures, and loads were obtained. The related conclusions are as follows:
(1)
During the transmission of tire cavity resonance noise in a Y-spoke aluminum alloy wheel and suspension system, the power flow output by the lower arm front bushing is the largest and the power flow output by the rear bushing is the smallest.
(2)
The areas in the suspension system with high stress are located at the steering knuckle, lower swing arm, and shock absorber. With an increase in load, the stress and strain of a suspension system increase, and the structural sound intensity of the suspension system increases.
(3)
The output power flow of a suspension system decreases slightly with an increase in load. In our study, the output power flow of the suspension system was the largest when the inflation pressure was 0.22 MPa; the output power flow was relatively low under the other four inflation pressures. As the speed increased, the output power flow of the suspension system decreased slightly.
Studying the propagation characteristics of tire cavity resonance noise in suspension systems can provide a basis for the design of wheels and suspension systems which could suppress the propagation of tire cavity resonance noise.

Author Contributions

Conceptualization, X.L. and X.H.; methodology, X.H.; software, X.H.; validation, X.L., X.H. and Y.S.; investigation, T.H.; data curation, X.H.; writing—original draft preparation, X.H.; writing—review and editing, X.L.; visualization, X.H.; supervision, X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of China grant number 51675021.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the requirements of the university.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arpanahi, R.A.; Eskandari, A.; Mohammadi, B.; Hashemi, S.H. Study on the effect of viscosity and fluid flow on buckling behavior of nanoplate with surface energy. Results Eng. 2023, 18, 101078. [Google Scholar] [CrossRef]
  2. Arpanahi, R.A.; Hashemi, K.H.; Mohammadi, B.; Hashemi, S.H. Investigation of the vibration behavior of nano piezoelectric rod using surface effects and non-local elasticity theory. Eng. Res. Express 2023, 5, 035029. [Google Scholar] [CrossRef]
  3. Arpanahi, R.A.; Mohammadi, B.; Ahmadian, M.T.; Hashemi, S.H. Study on the buckling behavior of nonlocal nanoplate submerged in viscous moving fluid. Int. J. Dyn. Control. 2023, 11, 2820–2830. [Google Scholar] [CrossRef]
  4. Deng, S.T. Investigation on Vibration Energy Transfer Characteristics of Beam Structures in Ladder Frame Chassis. Master’s Thesis, Hunan University, Changsha, China, 2021. [Google Scholar]
  5. Dzierzek, S. Experiment-based modeling of cylindrical rubber bushings for the simulation of wheel suspension dynamic behavior. SAE Trans. 2000, 109, 78–85. [Google Scholar]
  6. Feng, B. Methods for Modeling and Parameters Identification of Rubber Bushing. Master’s Thesis, Jilin University, Changchun, China, 2015. [Google Scholar]
  7. Gavrić, L.; Pavić, G. A finite element method for computation of structural intensity by the normal mode approach. J. Sound Vib. 1993, 164, 29–43. [Google Scholar] [CrossRef]
  8. Hu, X.; Liu, X.; Wan, X.; Shan, Y.; Yi, J. Experimental analysis of sound field in the tire cavity arising from the acoustic cavity resonance. Appl. Acoust. 2019, 161, 107172. [Google Scholar] [CrossRef]
  9. Hu, X.; Liu, X. Research on Tire Cavity Resonance Mechanism and Its Noise Characteristics. Doctoral Thesis, Beihang University, Beijing, China, 2021. [Google Scholar]
  10. Kim, G.J.; Holland, K.R.; Lalor, N. Identification of the airborne component of tyre-induced vehicle interior noise. Appl. Acoust. 1997, 51, 141–156. [Google Scholar] [CrossRef]
  11. Lyon, R.H.; Maidanik, G. Power flow between linearly coupled oscillators. J. Acoust. Soc. Am. 1962, 34, 623–639. [Google Scholar] [CrossRef]
  12. Li, Z. Analysis of the Performance of Rear Drive Axle Suspension with Rubber Bushing. Master’s Thesis, Nanjing University of Aeronautics and Astronautics, Nanjing, China, 2016. [Google Scholar]
  13. Li, X. Research on Fusion Control Strategy of Semi-Active Seat Suspension System for Both Vibration and Shock Mitigation. Master’s Thesis, Hefei University of Technology, Hefei, China, 2021. [Google Scholar]
  14. Li, X.P.; Li, F.J.; Yang, L.X.; Liu, X.L. Optimization design and dynamics analysis of vehicle suspension system. J. Northeast. Univ. (Nat. Sci.) 2020, 41, 1097–1102. [Google Scholar]
  15. Mohamed, Z.; Wang, X.; Jazar, R. A survey of wheel tyre cavity resonance noise. Int. J. Veh. Noise Vib. 2013, 9, 276–293. [Google Scholar] [CrossRef]
  16. Mohamed, Z.; Egab, L.; Wang, X. Tyre cavity coupling resonance and countermeasures. Appl. Mech. Mater. 2014, 471, 3–8. [Google Scholar] [CrossRef]
  17. Mohamed, Z. A Study of Tyre Cavity Resonance Noise Mechanism and Countermeasures Using Vibroacoustic Analysis. Doctoral Thesis, RMIT University, Melbourne, VIC, Australia, 2014. [Google Scholar]
  18. Pinnington, R.J.; White, R.G. Power flow through machine isolators to resonant and non-resonant beams. J. Sound Vib. 1981, 75, 179–197. [Google Scholar] [CrossRef]
  19. Qiu, S. Simulation and Analysis of Vehicle Ride Comfort Affected by Suspension Rubber Bushing. Master’s Thesis, Beijing Institute of Technology, Beijing, China, 2016. [Google Scholar]
  20. Stephen, A.H.; Richard, P.S. Predictions of structural intensity fields using solid finite elements. Noise Control Eng. J. 1999, 47, 209–217. [Google Scholar]
  21. Wu, Y.Y. Analysis and Optimization Design of Vibration Characteristics of Electric Vehicle Powertrain Mounting System. Master’s Thesis, Shandong University, Jinan, China, 2020. [Google Scholar]
  22. Weisser, T.; Foltête, E.; Bouhaddi, N.; Gonidou, L.O. A power flow mode approach dedicated to structural interface dynamic characterization. J. Sound Vib. 2015, 334, 202–218. [Google Scholar] [CrossRef]
  23. Zhao, P. Analysis and Optimization of the Rubber Bushing Influence on the K&C Characteristic of Suspension. Master’s Thesis, Chongqing University of Technology, Chongqing, China, 2013. [Google Scholar]
Figure 1. Micro-stress distribution.
Figure 1. Micro-stress distribution.
Applsci 13 12113 g001
Figure 2. McPherson suspension system model. (a) McPherson suspension system. (b) Aluminum alloy wheel and McPherson suspension system.
Figure 2. McPherson suspension system model. (a) McPherson suspension system. (b) Aluminum alloy wheel and McPherson suspension system.
Applsci 13 12113 g002
Figure 3. Load state of a McPherson suspension system.
Figure 3. Load state of a McPherson suspension system.
Applsci 13 12113 g003
Figure 4. Simplified McPherson suspension system and transverse arm bushings. (a) Simplified McPherson suspension. (b) Transverse arm front bushing. (c) Transverse arm rear bushing.
Figure 4. Simplified McPherson suspension system and transverse arm bushings. (a) Simplified McPherson suspension. (b) Transverse arm front bushing. (c) Transverse arm rear bushing.
Applsci 13 12113 g004
Figure 5. The first six orders of vibration mode shapes of the suspension system.
Figure 5. The first six orders of vibration mode shapes of the suspension system.
Applsci 13 12113 g005
Figure 6. Stress and strain contour of a suspension system. (a) Stress contour of a suspension system. (b) Strain contour of a suspension system.
Figure 6. Stress and strain contour of a suspension system. (a) Stress contour of a suspension system. (b) Strain contour of a suspension system.
Applsci 13 12113 g006
Figure 7. Power flow contour of aluminum alloy wheel and suspension system.
Figure 7. Power flow contour of aluminum alloy wheel and suspension system.
Applsci 13 12113 g007
Figure 8. Distribution of the sound intensity vector of a suspension system in different planes and with partial enlargement. (a) XZ plane sound intensity vector distribution diagram and partial enlargement. (b) XY plane sound intensity vector distribution diagram and partial enlargement. (c) YZ plane sound intensity vector distribution diagram and partial enlargement.
Figure 8. Distribution of the sound intensity vector of a suspension system in different planes and with partial enlargement. (a) XZ plane sound intensity vector distribution diagram and partial enlargement. (b) XY plane sound intensity vector distribution diagram and partial enlargement. (c) YZ plane sound intensity vector distribution diagram and partial enlargement.
Applsci 13 12113 g008aApplsci 13 12113 g008b
Figure 9. Loading state of a suspension system under different loads.
Figure 9. Loading state of a suspension system under different loads.
Applsci 13 12113 g009
Figure 10. Stress contour of a suspension system under different loads. (a) Stress contour under 2500 N load. (b) Stress contouar under 3000 N load. (c) Stress contour under 3500 N load. (d) Stress contour under 4000 N load. (e) Stress contour under 4500 N load.
Figure 10. Stress contour of a suspension system under different loads. (a) Stress contour under 2500 N load. (b) Stress contouar under 3000 N load. (c) Stress contour under 3500 N load. (d) Stress contour under 4000 N load. (e) Stress contour under 4500 N load.
Applsci 13 12113 g010
Figure 11. Strain contour of a suspension system under different loads. (a) Strain contour under 2500 N load. (b) Strain contour under 3000 N load. (c) Strain contour under 3500 N load. (d) Strain contour under 4000 N load. (e) Strain contour under 4500 N load.
Figure 11. Strain contour of a suspension system under different loads. (a) Strain contour under 2500 N load. (b) Strain contour under 3000 N load. (c) Strain contour under 3500 N load. (d) Strain contour under 4000 N load. (e) Strain contour under 4500 N load.
Applsci 13 12113 g011
Figure 12. Power flow contour of a suspension system under different loads. (a) Power flow contour under 2500 N load. (b) Power flow contour under 3000 N load. (c) Power flow contour under 3500 N load. (d) Power flow contour under 4000 N load. (e) Power flow contour under 4500 N load.
Figure 12. Power flow contour of a suspension system under different loads. (a) Power flow contour under 2500 N load. (b) Power flow contour under 3000 N load. (c) Power flow contour under 3500 N load. (d) Power flow contour under 4000 N load. (e) Power flow contour under 4500 N load.
Applsci 13 12113 g012
Table 1. Stiffness and damping of the bushing and spring unit.
Table 1. Stiffness and damping of the bushing and spring unit.
LocationDirectionStiffness (N/Mm)Damping (N·s/mm)
Shock absorber upper endX10,0001.5
Y10,000
Z10,000
Spring and shock absorber 251
Lower transverse arm front bushingX10,0002
Y10,000
Z1000
Lower transverse arm rear bushingX5000.8
Y500
Z100
Table 2. Natural frequency of suspension system in each order.
Table 2. Natural frequency of suspension system in each order.
OrderFirstSecondThirdFourthFifthSixth
Frequency (Hz)28.4460.3081.56131.62201.82245.92
Table 3. Finite element calculation results of input and output sectional power flow.
Table 3. Finite element calculation results of input and output sectional power flow.
Load ConditionSum of Structural Sound Intensity (10−3 W/mm2)Number of PointsCross Sectional Area (mm2)Power Flow of Cross Section (10−3 W)
4000 NInput2.897 × 10−215,06481,3671.565 × 10−1
First output1.068 × 10−4184070314.081 × 10−4
Second output1.153 × 10−446227976.980 × 10−4
Third output1.808 × 10−535013520.698 × 10−4
Table 4. Sound pressure amplitude function of wheels under various working conditions.
Table 4. Sound pressure amplitude function of wheels under various working conditions.
Tire Working ConditionSound Pressure Amplitude Function
0.22 MPa 60 km/h2500 N P θ , t = ( 43 18 cos θ + 19 sin θ + 12 cos 2 θ 5 sin 2 θ cos 3 θ 4 sin 3 θ ) cos ( 2 π f t )
3000 N P θ , t = ( 55 9 cos θ + 45 sin θ 2 cos 2 θ 27 sin 2 θ + 8 cos 3 θ + sin 3 θ ) cos ( 2 π f t )
3500 N P θ , t = ( 88 34 cos θ + 6 sin θ + 18 cos 2 θ 6 sin 2 θ ) cos ( 2 π f t )
4000 N P θ , t = ( 80 30 cos θ + 30 sin θ + 32 cos 2 θ 24 sin 2 θ + 3 cos 3 θ + 14 sin 3 θ ) cos ( 2 π f t )
4500 N P θ , t = ( 92 37 cos θ + 34 sin θ + 29 cos 2 θ 35 sin 2 θ ) cos ( 2 π f t )
3500 N
50 km/h
0.16 MPa P θ , t = ( 33 8 cos θ + 7 sin θ + 13 cos 2 θ 11 sin 2 θ + cos 3 θ + 6 sin 3 θ ) cos ( 2 π f t )
0.19 MPa P θ , t = ( 53 11 cos θ + 7 sin θ + 5 cos 2 θ 11 sin 2 θ + cos 3 θ 6 sin 3 θ ) cos ( 2 π f t )
0.22 MPa P θ , t = ( 49 20 cos θ + 6 sin θ + 8 cos 2 θ 8 sin 2 θ ) cos ( 2 π f t )
0.25 MPa P θ , t = ( 39 cos θ + 11 sin θ + 12 cos 2 θ 14 sin 2 θ 0.1 cos 3 θ + 7 sin 3 θ ) cos ( 2 π f t )
0.28 MPa P θ , t = ( 45 12 cos θ + 5 sin θ + 9 cos 2 θ 11 sin 2 θ cos 3 θ 5 sin 3 θ ) cos ( 2 π f t )
0.22 MPa
3000 N
20 km/h P θ , t = ( 28 + 5 cos θ + 11 sin θ + 2 cos 2 θ 12 sin 2 θ ) cos ( 2 π f t )
30 km/h P θ , t = ( 39 + 2 cos θ + 12 sin θ + 3 cos 2 θ 13 sin 2 θ ) cos ( 2 π f t )
40 km/h P θ , t = ( 48 + cos θ + 13 sin θ + 2 cos 2 θ 12 sin 2 θ ) cos ( 2 π f t )
50 km/h P θ , t = ( 48 9 cos θ + 14 sin θ + 2 cos 2 θ 10 sin 2 θ ) cos ( 2 π f t )
60 km/h P θ , t = ( 55 9 cos θ + 45 sin θ 2 cos 2 θ 27 sin 2 θ + 8 cos 3 θ + sin 3 θ ) cos ( 2 π f t )
Table 5. Finite element calculation results of total power flow of input and three output sections.
Table 5. Finite element calculation results of total power flow of input and three output sections.
Load ConditionSum of Structural Sound Intensity (10−3 W/mm2)Number of PointsCross Sectional Area (mm2)Power Flow of Cross Section (10−3 W)
2500 NInput2.831 × 10−215,06481,3671.529 × 10−1
Output2.793 × 10−4265211,1801.177 × 10−3
3000 NInput2.867 × 10−215,06481,3671.549 × 10−1
Output2.627 × 10−4265211,1801.107 × 10−3
3500 NInput2.848 × 10−215,06481,3671.538 × 10−1
Output2.503 × 10−4265211,1801.055 × 10−3
4000 NInput2.897 × 10−215,06481,3671.565 × 10−1
Output2.402 × 10−4265211,1801.013 × 10−3
4500 NInput2.975 × 10−215,06481,3671.607 × 10−1
Output2.290 × 10−4265211,1809.652 × 10−4
Table 6. Simulation results of sum power flow of input and three output sections.
Table 6. Simulation results of sum power flow of input and three output sections.
Inflation Pressure ConditionSum of Structural Sound Intensity (10−3 W/mm2)Number of PointsCross Sectional Area (mm2)Power Flow of Cross Section (10−3 W)
0.16 MPaInput2.849 × 10−215,06481,3671.539 × 10−1
Output2.603 × 10−4265211,1801.097 × 10−3
0.19 MPaInput2.830 × 10−215,06481,3671.529 × 10−1
Output2.648 × 10−4 265211,1801.116 × 10−3
0.22 MPaInput2.811 × 10−215,06481,3671.518 × 10−1
Output2.702 × 10−4265211,1801.139 × 10−3
0.25 MPaInput2.821 × 10−215,06481,3671.523 × 10−1
Output2.662 × 10−4265211,1801.122 × 10−3
0.28 MPaInput2.836 × 10−215,06481,3671.532 × 10−1
Output2.615 × 10−4 265211,1801.102 × 10−3
Table 7. Simulation results of total power flow of input and three output sections.
Table 7. Simulation results of total power flow of input and three output sections.
Speed ConditionSum of Structural Sound Intensity (10−3 W/mm2)Number of PointsCross Sectional Area (mm2)Power Flow of Cross Section (10−3 W)
20 km/hInput2.566 × 10−215,06481,3671.386 × 10−1
Output2.958 × 10−4 265211,1801.247 × 10−3
30 km/hInput2.635 × 10−215,06481,3671.423 × 10−1
Output2.894 × 10−4 265211,1801.220 × 10−3
40 km/hInput2.701 × 10−215,06481,3671.459 × 10−1
Output2.815 × 10−4 265211,1801.187 × 10−3
50 km/hInput2.782 × 10−215,06481,3671.503 × 10−1
Output2.722 × 10−4265211,1801.148 × 10−3
60 km/hInput2.867 × 10−215,06481,3671.549 × 10−1
Output2.627 × 10−4265211,1801.107 × 10−3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hu, X.; Liu, X.; Shan, Y.; He, T. Research on Propagation Characteristics of Tire Cavity Resonance Noise in the Automobile Suspension. Appl. Sci. 2023, 13, 12113. https://doi.org/10.3390/app132212113

AMA Style

Hu X, Liu X, Shan Y, He T. Research on Propagation Characteristics of Tire Cavity Resonance Noise in the Automobile Suspension. Applied Sciences. 2023; 13(22):12113. https://doi.org/10.3390/app132212113

Chicago/Turabian Style

Hu, Xiaojun, Xiandong Liu, Yingchun Shan, and Tian He. 2023. "Research on Propagation Characteristics of Tire Cavity Resonance Noise in the Automobile Suspension" Applied Sciences 13, no. 22: 12113. https://doi.org/10.3390/app132212113

APA Style

Hu, X., Liu, X., Shan, Y., & He, T. (2023). Research on Propagation Characteristics of Tire Cavity Resonance Noise in the Automobile Suspension. Applied Sciences, 13(22), 12113. https://doi.org/10.3390/app132212113

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop