Research on Low-Frequency Noise Control of Automobiles Based on Acoustic Metamaterial

With the transformation of the trend of vehicle electrification, the overall noise level in the vehicle is gradually reduced. The problem of low-frequency noise in the vehicle, which was previously ignored, is becoming more and more prominent. To solve the vehicle low-frequency noise problem, a combination of real-vehicle tests and simulation analysis is carried out. During the test, the driver and passengers feel that there is a relatively obvious low-frequency roar noise in the car, which results from the structural radiation noise of the trunk door vibration. Therefore, to solve this problem, we design an acoustic metamaterial with lightweight and miniaturized features based on the local resonance principle of phononic crystals. Firstly, the selection of the resonant unit configuration and the design of the band gap are implemented. Then, the layout planning of the whole vehicle, the layout of the resonance unit and the design of the base frame are implemented. The actual vehicle test results show that: after attaching the designed acoustic metamaterial, the low-frequency noise sound pressure levels in the front and rear of the vehicle were reduced by 2.0 dB (A) and 2.3 dB (A), respectively, meanwhile, the interior noise sound quality was improved. The sound pressure level at the driver’s right ear in the car has an abnormal peak of around 35Hz. In addition, the driver and passengers feel that there is a relatively obvious low-frequency roar noise in the car, and through low-pass filtering of the collected signals, it is confirmed that the peak frequency is the main cause of the low-frequency roar in the car. The low-frequency steady-state noise of the car is generally considered to be the low-frequency vibration of the body panel and the radiation occurs. Through the finite element simulation analysis (Grid Participation Analysis) of the abnormal peak frequency, the results show that the low-frequency roar is caused by the low-frequency vibration of the tailgate sheet metal, and the problem peak frequency is not coupled with the acoustic cavity mode. Facing the problem of the low-frequency roar radiated into the car by the vibration of the tailgate sheet metal parts, based on the local resonance band gap theory, we developed a design to suppress the 35 Hz vibration of the tailgate sheet metal parts and meet the characteristics of lightweight and miniaturization. By attaching the acoustic metamaterial to the tailgate and performing CAE simulation of the whole vehicle, it is determined that the structure can indeed reduce the 35 Hz noise in the car and the peak value of the tailgate sheet metal vibration.


Introduction
In the field of automotive NVH (Noise, Vibration, Harshness), high-frequency noise and vibration control can achieve ideal results through acoustic packaging or attaching

Introduction
In the field of automotive NVH (Noise, Vibration, Harshness), high-frequency noise and vibration control can achieve ideal results through acoustic packaging or attaching vibration damping, but controlling the radiation caused by the vibration of automotive sheet metal parts low-frequency noise has always been a major and difficult problem in the field of automotive noise and vibration control [1][2][3]. For automobiles, the low-frequency noise and vibration problems radiated into the car caused by the vibration of sheet metal parts will not only seriously affect the ride comfort of the driver and passengers, but also have a great impact on the driving safety and service life of the car [4]. In order to obtain a more comfortable and safe driving and riding environment, and to improve the market competitiveness of vehicles, it is extremely important to control the low-frequency noise radiated into the vehicle by the vibration of the automobile sheet metal parts [5]. The generation path of low-frequency noise radiated by automobile (tailgate) sheet metal parts is shown in Figure 1. The car is driving on rough roads, and the displacement excitation of the road is transmitted to the body sheet metal parts through the tire-suspension system [6]. Excitation generates low-frequency vibration and radiates low-frequency structural noise into the car, causing serious discomfort to passengers and drivers [7,8]. At present, there are three main ways to reduce the vibration of sheet metal parts in response to the problem of low-frequency noise radiated by the vibration of automobile sheet metal parts at home and abroad [9,10]. First, using spring-damping vibration-reducing materials or rubber vibration-reducing materials. With the development of current industrial technology, it is found that the use of spring-damping vibration-reducing materials or rubber damping vibration-reducing materials has a more obvious effect on highfrequency vibration, but it is difficult to achieve an ideal suppression effect of the middle and low-frequency noise [11]. Second, applying reinforcing ribs to the sheet metal parts to increase the natural frequency of the low-frequency band of the sheet metal structure by increasing the rigidity, but the structure of the automobile sheet metal parts is greatly changed which will cause the production cost of the model to increase sharply, and seriously affect the economic benefits of the enterprise [12]. Third, installing the dynamic vibration absorber, but its cost is higher and the mass is higher [13]. In view of the above problems, the local resonance type vibration-damping acoustic metamaterial with negative dynamic equivalent parameter characteristics provides a new idea for noise control of low-frequency vibration radiation of automobile sheet metal parts [14][15][16][17][18][19].
In recent years, studies on the locally resonant phonon crystal formed by the periodic arrangement of elastomers and mass scatterers have been studied with great vigor [20]. Since the band gap produced by it is not limited by the size of the unit cell, it has a broad application market for low-frequency vibration and noise reduction in automobiles. In At present, there are three main ways to reduce the vibration of sheet metal parts in response to the problem of low-frequency noise radiated by the vibration of automobile sheet metal parts at home and abroad [9,10]. First, using spring-damping vibration-reducing materials or rubber vibration-reducing materials. With the development of current industrial technology, it is found that the use of spring-damping vibration-reducing materials or rubber damping vibration-reducing materials has a more obvious effect on high-frequency vibration, but it is difficult to achieve an ideal suppression effect of the middle and lowfrequency noise [11]. Second, applying reinforcing ribs to the sheet metal parts to increase the natural frequency of the low-frequency band of the sheet metal structure by increasing the rigidity, but the structure of the automobile sheet metal parts is greatly changed which will cause the production cost of the model to increase sharply, and seriously affect the economic benefits of the enterprise [12]. Third, installing the dynamic vibration absorber, but its cost is higher and the mass is higher [13]. In view of the above problems, the local resonance type vibration-damping acoustic metamaterial with negative dynamic equivalent parameter characteristics provides a new idea for noise control of low-frequency vibration radiation of automobile sheet metal parts [14][15][16][17][18][19].
In recent years, studies on the locally resonant phonon crystal formed by the periodic arrangement of elastomers and mass scatterers have been studied with great vigor [20]. Since the band gap produced by it is not limited by the size of the unit cell, it has a broad application market for low-frequency vibration and noise reduction in automobiles. In 2000, when Liu Zhengyou et al. studied the three-dimensional phononic crystal formed by the cubic lattice structure of lead balls coated with viscoelastic silicone rubber soft material and embedded in epoxy resin, they found that the wavelength of the phononic crystal forbidden band is much larger [21]. Due to the size of the lattice, it breaks through the limitation of the Bragg scattering mechanism and can control elastic waves two orders of magnitude lower than the resonant unit. Moreover, the scatterers are not strictly periodic or even randomly distributed, the composite structure also has a band gap effect [22]. This puts forward the local resonance mechanism of the elastic wave band gap, which provides a new theoretical basis for the follow-up automotive low-frequency noise and vibration control field [23]. In 2001, BAZ introduced active control to study the transmission characteristics of elastic waves in periodic beams with two materials alternately arranged, gave the corresponding simplified theoretical analysis model of the spring-mass system, explained the active adjustment method of passband and forbidden band, and proposed the directional design method of local resonance type acoustic metamaterial [24]. In 2008, Hong Kong University of Science and Technology Yang Zhiyu and others proposed a thin-film acoustic metamaterial based on the local resonance theory [25][26][27]. The lightweight structure with a thickness as low as millimeters successfully achieved outstanding sound insulation capabilities in the low-frequency range of 500~800 Hz, which provided a new solution for low-frequency noise reduction. The elastic film needs external tension to propagate the vibration [28]. As the stiffness of the elastic film is very low and the performance is unstable, a small change in tension will cause a frequency shift of tens or even hundreds of hertz. Therefore, no matter from a theoretical or technical point of view, effective control of film tension is very necessary. In 2015, Ma Fuyin and others at Xi'an Jiaotong University proposed a square lattice thin-film acoustic sound insulation metamaterial that requires almost no tension by increasing the thickness of the film [29]. In 2019, Niu Jiamin and others from Xi'an Jiaotong University proposed that under the asymmetric mode, the overall sound absorption performance of the structure has been greatly improved, and the sound absorption bandwidth has been expanded from the original narrow frequency to a wide frequency [30].
In summary, although local resonance acoustic metamaterials seem to be relatively new products now, there are a lot of data on the research of its mechanism, and there are also certain physical samples to verify its vibration suppression effect. However, the current acoustic metamaterials developed by technical institutions have poor stability and the vibration suppression frequency is difficult to reduce to less than 100 Hz. The acoustic metamaterial research is almost limited to simulation design and theoretical research.
Focusing on the problem of the low-frequency roar radiated into the car by the vibration of the tailgate sheet metal of a certain vehicle model, based on the theory of local resonance bandgap, a design is designed to suppress the 35 Hz vibration of the tailgate sheet metal with stable and small equivalent stiffness. Through the vehicle finite element simulation, it is verified that the structure can indeed reduce the 35 Hz noise peak in the vehicle and reduce the vibration of the tailgate sheet metal parts, thereby providing a new idea and new solution for the field of vehicle low-frequency noise and vibration control.

Elastic Wave Equation in Acoustic Metamaterial
An acoustic metamaterial is a kind of man-made periodic composite structure material, which is developed from crystals in solid mechanics [31,32]. Therefore, the theory of the energy band of crystals in solid mechanics is also applicable to periodic composite structures. The content of acoustic metamaterial exploration is the propagation of elastic waves in the medium, and the band gap is due to the special effects of elastic waves when they are transmitted in the structure [33]. According to the knowledge of elastic dynamics theory, under the condition of ensuring an ideal elastic medium, small displacement and no initial stress, consider any particle in the homogeneous acoustic metamaterial medium, and use the displacement method used in solving dynamic problems to accurately establish three types of equations describing the relationship between the mass point force, displacement, stress and strain in the acoustic metamaterial are as follows: Differential Equation of Motion: Geometric Equation: Physical Equation: The above equation involves the six components of the stress tensor (σ x , σ y , σ z , τ yz , τ zx , τ xy ), six components of strain tensor (γ xy , γ xz , γ yz , ε x , ε y , ε z ) and three displacement components (u, v, w), a total of 15 unknowns of spatial coordinates x, y, z and time t. σ x is the normal stress, τ xy is the shear stress, ρ is mass point density, γ x is shear strain, ε x is the normal strain, λ, µ are Lame coefficients, θ t = ε x + ε y + ε z .
Not every equation contains the unknowns in the above equation, by using a displacement method often used to solve the equation of elastic dynamics, considering some of these unknown function hypotheses into the mathematic expression of the "basic unknown functions", it can be brought into the rest of the equation, and can obtain the expressed in the displacement of stress components, the final three displacements as unknown functions derived the elastic dynamics equation: In the equation, → u is the displacement vector; λ, µ are the Lame coefficient; ρ is the density of the medium; ∇ is the Laplace operator; the relation between them and the wave velocity is expressed as: Longitudinal wave velocity: C l = (λ + 2µ)/ρ, Shear wave velocity: C t = µ/ρ.

Acoustic Metamaterial Band Gap Mechanism
When the vibrational elastic wave propagates in the local resonant phonon crystal, it will be affected by the periodic elastic scatterers, and the elastic wave within a certain intrinsic frequency range cannot continue to propagate through the elastic scatterers, which is called the forbidden band of the acoustic metamaterial. Other eigenfrequency ranges that can pass through these elastic scatterers without obstruction are called passbands. As shown in Figure 2a, it is a typical phonon crystal structure with local resonance. The structure is composed of mass scatterers, an elastic coating layer and matrix. The structure can be simplified to a typical spring-mass system, as shown in Figure 2b. The mass scatterers are composed of a dense metal material, which provides the mass component (expressed as m) for the cell structure. The elastic coating layer is composed of a thin film with certain elasticity, which provides the elastic component of the crystal cell structure (expressed as k). The matrix is rigid material (expressed as M), which is the carrier of vibration elastic wave transmission [34][35][36][37].
The structure is composed of mass scatterers, an elastic coating layer and matrix. The structure can be simplified to a typical spring-mass system, as shown in Figure 2b. The mass scatterers are composed of a dense metal material, which provides the mass component (expressed as m) for the cell structure. The elastic coating layer is composed of a thin film with certain elasticity, which provides the elastic component of the crystal cell structure (expressed as k). The matrix is rigid material (expressed as M), which is the carrier of vibration elastic wave transmission [34][35][36][37].
Consider the base body and spring-mass system as a whole: In the equation: is the system equivalent mass, ( ) is the system displacement frequency response function, ω is the natural circular frequency, 0 = √ ⁄ is the natural frequency of the internal spring oscillator. It can be seen from the above formula that: with the change of external excitation frequency, the system shows different characteristics due to its different equivalent mass, as shown in Figure 3. Suppose the displacement of matrix M and mass component m is X and x, the matrix excitation is F, and the reaction force of the mass component is F 1 . By Newton's second law and Hooke's law: Consider the base body and spring-mass system as a whole: In the equation: m e is the system equivalent mass, H(w) is the system displacement frequency response function, ω is the natural circular frequency, w 0 = √ k /m is the natural frequency of the internal spring oscillator.
It can be seen from the above formula that: with the change of external excitation frequency, the system shows different characteristics due to its different equivalent mass, as shown in Figure 3.  According to Newton's second law and Hooke's law theory and the external excitation frequency and displacement frequency response function, the relationship between the dynamic equivalent quality of the system and the excitation frequency can be inferred, as shown in Table 1.   According to Newton's second law and Hooke's law theory and the external excitation frequency and displacement frequency response function, the relationship between the dynamic equivalent quality of the system and the excitation frequency can be inferred, as shown in Table 1. Table 1. The relationship between system dynamic equivalent mass m e and excitation frequency w.

Excitation Frequency Range Dynamic Equivalent Mass System Movement Status
M and m resonate, and the system vibration tends to be static The equivalent mass is negative, M and m move in opposite directions In conclusion, the band gap is caused by the resonance of the internal mass component. The band gap range is determined by the natural frequency of the whole structure and internal mass component. The band gap effect occurs when the system is in the negative equivalent mass phase, the excitation frequency is in the range of w 0 , w 0 √ (M + m) /M . Therefore, the frequency range generated by bandgap can be adjusted directionally by changing the mass or spring of the equivalent simplified model of the cell unit. Based on the above theory, according to the vibration frequency of auto sheet metal parts (equivalent to excitation frequency), an acoustic metamaterial can be designed to suppress the specific frequency.

Analysis of Low-Frequency Roaring in Car
Aiming at the relatively obvious low-frequency roar (20-100 Hz) generated by a certain hatch car when it is running at a constant speed of 30 km/h on a rough road surface, the LMS equipment is used to collect the sound pressure level at the driver's right ear and at the middle back passenger seat. Sound pressure sensors are arranged according to relevant national standards GB/T18697-2002, and the location of measuring points is shown in Figure 4. During the signal acquisition process, the adjustable seat is back as far as possible to put it in a vertical position. The test is empty, with only the driver and measurement of personnel in the car. Before the test we ensure that the car is in a closed state, the four-wheel positioning is normal and in the test environment, the influence of background noise is less than the allowable error. When the vehicle speed is adjusted to a constant speed of 30 km/h, LMS onboard equipment is used to collect the sound pressure signal inside the vehicle. The noise test results are shown in Figure 5. metal parts (equivalent to excitation frequency), an acoustic metamaterial can be designed to suppress the specific frequency.

Analysis of Low-Frequency Roaring in Car
Aiming at the relatively obvious low-frequency roar (20-100 Hz) generated by a certain hatch car when it is running at a constant speed of 30 km/h on a rough road surface, the LMS equipment is used to collect the sound pressure level at the driver's right ear and at the middle back passenger seat. Sound pressure sensors are arranged according to relevant national standards GB/T18697-2002, and the location of measuring points is shown in Figure 4. During the signal acquisition process, the adjustable seat is back as far as possible to put it in a vertical position. The test is empty, with only the driver and measurement of personnel in the car. Before the test we ensure that the car is in a closed state, the four-wheel positioning is normal and in the test environment, the influence of background noise is less than the allowable error. When the vehicle speed is adjusted to a constant speed of 30 km/h, LMS onboard equipment is used to collect the sound pressure signal inside the vehicle. The noise test results are shown in Figure 5.    The test results show that the car has an obvious peak value at 35 Hz in the frequency range of 20-100 Hz. The peak value at 35 Hz in the right ear of the driver is 50.39 dB (A), and the peak value at 35 Hz in the middle of the passenger rear seat is 50.78 dB (A), as shown in Figure 4 (ordinate of Figure 4a is dB (A) and ordinate of Figure 4b is Pa). In the LMS test software, the collected signals were filtered to filter the peak value at 35 Hz, and the filtered noise signals were listened to by earphones, and it was found that the lowfrequency roar had disappeared. Therefore, it can be determined that the low-frequency roar of this vehicle is caused by the peak noise at 35 Hz. As the peak frequency is very low, the continuous steady-state noise of low frequency is generally considered structural noise [38,39].

Analysis of Interior Noise Problem
Generally speaking, there are two ways to generate structural noise: the sound cavity mode of the whole vehicle and the coupled resonance of the body parts amplify the noise to form a low-frequency roar; the low-frequency vibration of body sheet metal radiates to the car, forming a low-frequency roar. In view of the above two approaches, this paper solves the acoustic cavity mode of the vehicle and analyzes the joint contribution of the panel parts based on the peak noise frequency of 35 Hz in the simulation analysis of road noise. The analysis process is shown in Figure 6.   Figure 4b is Pa). In the LMS test software, the collected signals were filtered to filter the peak value at 35 Hz, and the filtered noise signals were listened to by earphones, and it was found that the low-frequency roar had disappeared. Therefore, it can be determined that the low-frequency roar of this vehicle is caused by the peak noise at 35 Hz. As the peak frequency is very low, the continuous steady-state noise of low frequency is generally considered structural noise [38,39].

Analysis of Interior Noise Problem
Generally speaking, there are two ways to generate structural noise: the sound cavity mode of the whole vehicle and the coupled resonance of the body parts amplify the noise to form a low-frequency roar; the low-frequency vibration of body sheet metal radiates to the car, forming a low-frequency roar. In view of the above two approaches, this paper solves the acoustic cavity mode of the vehicle and analyzes the joint contribution of the panel parts based on the peak noise frequency of 35 Hz in the simulation analysis of road noise. The analysis process is shown in Figure 6. Through CAE simulation analysis, the sound mode of the vehicle has no resonance with the body parts at 35 Hz. Therefore, the peak noise frequency of 35 Hz is mainly caused by the vibration radiation of the automobile panel.

The Establishment of Finite Element Model
In order to carry out acoustic cavity modal analysis, road noise simulation and body plate joint contribution analysis of the whole vehicle, it is necessary to establish the sound cavity finite element model and acoustic-solid coupled finite element model.
The TB (Trim Body) body is further constructed on the basis of the finite element model of body-in-white that has been constructed and modified through tests, which includes all parts except the engine and chassis system, many of which are replaced by centralized unit mass and one-dimensional unit. Determine the modeling standards and guidelines between the connections of various components: (1) Modeling of spot welding structural adhesives. Spot welding structural adhesives are mainly used in key parts of the car body, which can greatly improve the stiffness of the car body when used in conjunction with welding spots. Modeling of RBE3-hexa-RBE3 for spot welding of structural adhesive. (2) Modeling of shock-absorbing expansion adhesive. Shock-absorbing expansion adhesive mainly connects the outer plate of the top cover with the beam of the top cover, the outer plate of the open and closed parts with the support plate or the inner plate, and mainly plays a shock-absorbing role. (3) Modeling of reinforcement plate and damping plate. A reinforcement glue plate is generally used for side circumference and rear door plates with a thickness of 1.5 mm. It mainly plays the role of stiffness reinforcement and has a certain damping dissipation function. (4) Seal strip modeling. The seal strip is simulated using BUSH (used to define the nominal property values for a generalized spring and damper structural element) or ELAST (used to define the stiffness and stress coefficient of a scalar elastic element) in three translational directions. (5) Limit buffer block modeling. A limit buffer block is mainly used in the engine room cover and Through CAE simulation analysis, the sound mode of the vehicle has no resonance with the body parts at 35 Hz. Therefore, the peak noise frequency of 35 Hz is mainly caused by the vibration radiation of the automobile panel.

The Establishment of Finite Element Model
In order to carry out acoustic cavity modal analysis, road noise simulation and body plate joint contribution analysis of the whole vehicle, it is necessary to establish the sound cavity finite element model and acoustic-solid coupled finite element model.
The TB (Trim Body) body is further constructed on the basis of the finite element model of body-in-white that has been constructed and modified through tests, which includes all parts except the engine and chassis system, many of which are replaced by centralized unit mass and one-dimensional unit. Determine the modeling standards and guidelines between the connections of various components: (1) Modeling of spot welding structural adhesives. Spot welding structural adhesives are mainly used in key parts of the car body, which can greatly improve the stiffness of the car body when used in conjunction with welding spots. Modeling of RBE3-hexa-RBE3 for spot welding of structural adhesive.
(2) Modeling of shock-absorbing expansion adhesive. Shock-absorbing expansion adhesive mainly connects the outer plate of the top cover with the beam of the top cover, the outer plate of the open and closed parts with the support plate or the inner plate, and mainly plays a shock-absorbing role. (3) Modeling of reinforcement plate and damping plate. A reinforcement glue plate is generally used for side circumference and rear door plates with a thickness of 1.5 mm. It mainly plays the role of stiffness reinforcement and has a certain damping dissipation function. (4) Seal strip modeling. The seal strip is simulated using BUSH (used to define the nominal property values for a generalized spring and damper structural element) or ELAST (used to define the stiffness and stress coefficient of a scalar elastic element) in three translational directions. (5) Limit buffer block modeling. A limit buffer block is mainly used in the engine room cover and back door, to play a shock absorption and support role. The BUSH element is used to simulate and only three translational stiffnesses are set. (6) Connection between car door lock and lock. The BUSH unit can be used to simulate the model car door lock and lock. (7) Hinge modeling for open and closed parts. The hinge needs solid grid modeling, the grid size can be set to 3-5 mm, and the solid grid requires at least three layers of body elements in the direction of thickness. If the hinge is modeled using a solid grid or a shell element, the stiffness of the hinge will be overestimated. Doors generally contain two hinges connected by RBE2 and CBEAM, and the RBE2 center of the two hinges should be set to be collinear to ensure doors can rotate freely around the hinge. The hinge shaft is simulated using three CBEAM units and releases rotational degrees of freedom at both ends. The TB finite element model established is shown in Figure 7.  7) Hinge modeling for open and closed parts. The hinge needs solid grid modeling, the grid size can be set to 3-5 mm, and the solid grid requires at least three layers of body elements in the direction of thickness. If the hinge is modeled using a solid grid or a shell element, the stiffness of the hinge will be overestimated. Doors generally contain two hinges connected by RBE2 and CBEAM, and the RBE2 center of the two hinges should be set to be collinear to ensure doors can rotate freely around the hinge. The hinge shaft is simulated using three CBEAM units and releases rotational degrees of freedom at both ends. The TB finite element model established is shown in Figure 7. In Hypermesh 2017.01 software, Altair Engineering Inc.,we extract the inner surface of the body's internal contact with air, to form a completely closed acoustic cavity model, which defines the spoke boundary as a fully rigid wall. When dividing the finite element mesh, it should be noted that the size of the acoustic element and the structural element should be consistent, and the acoustic element should have at least six elements in each wavelength range. As we are focusing on frequency within 100 Hz, the minimum wavelength λ = c/f = 344/100 = 3.44 m, and the maximum length of the sound cavity grid element is 0.52 m. In this paper, a hexahedral element with a side length of 50 mm is used to automatically divide the cavity mesh. The generated results of the finite element model of the body cavity are shown in Figure 8.  In Hypermesh 2017.01 software, Altair Engineering Inc. (Troy, MI, USA), we extract the inner surface of the body's internal contact with air, to form a completely closed acoustic cavity model, which defines the spoke boundary as a fully rigid wall. When dividing the finite element mesh, it should be noted that the size of the acoustic element and the structural element should be consistent, and the acoustic element should have at least six elements in each wavelength range. As we are focusing on frequency within 100 Hz, the minimum wavelength λ = c/f = 344/100 = 3.44 m, and the maximum length of the sound cavity grid element is 0.52 m. In this paper, a hexahedral element with a side length of 50 mm is used to automatically divide the cavity mesh. The generated results of the finite element model of the body cavity are shown in Figure 8.  7) Hinge modeling for open and closed parts. The hinge needs solid grid modeling, the grid size can be set to 3-5 mm, and the solid grid requires at least three layers of body elements in the direction of thickness. If the hinge is modeled using a solid grid or a shell element, the stiffness of the hinge will be overestimated. Doors generally contain two hinges connected by RBE2 and CBEAM, and the RBE2 center of the two hinges should be set to be collinear to ensure doors can rotate freely around the hinge. The hinge shaft is simulated using three CBEAM units and releases rotational degrees of freedom at both ends. The TB finite element model established is shown in Figure 7. In Hypermesh 2017.01 software, Altair Engineering Inc.,we extract the inner surface of the body's internal contact with air, to form a completely closed acoustic cavity model, which defines the spoke boundary as a fully rigid wall. When dividing the finite element mesh, it should be noted that the size of the acoustic element and the structural element should be consistent, and the acoustic element should have at least six elements in each wavelength range. As we are focusing on frequency within 100 Hz, the minimum wavelength λ = c/f = 344/100 = 3.44 m, and the maximum length of the sound cavity grid element is 0.52 m. In this paper, a hexahedral element with a side length of 50 mm is used to automatically divide the cavity mesh. The generated results of the finite element model of the body cavity are shown in Figure 8.  The body structure and the sound cavity system are coupled through the vibration of the boundary nodes. In Hypermesh2017.01 software, the coupling between the structure and the sound cavity is realized by using the ACMODF card to search the structure nodes automatically through the nodes on the sound cavity boundary. The acoustic-solid coupling model of the body system is established as shown in Figure 9 (the left door and part of the side circumference are hidden to display the acoustic cavity model). The body structure and the sound cavity system are coupled through the vibration of the boundary nodes. In Hypermesh2017.01 software, the coupling between the structure and the sound cavity is realized by using the ACMODF card to search the structure nodes automatically through the nodes on the sound cavity boundary. The acoustic-solid coupling model of the body system is established as shown in Figure 9 (the left door and part of the side circumference are hidden to display the acoustic cavity model).

Simulation Results of Road Noise
Based on the acoustic-solid coupled finite element model, the time-domain signals of acceleration measured by engineers on the actual road surface at the steering knuckle of the front and rear suspension are converted into PSD frequency-domain signal files to stimulate the four wheel cores of the vehicle, and the solution response points are the sound pressure level at the driver's right ear and the middle sound pressure level of the rear passenger's seat. The solution frequency range is set as 20-100 Hz. The comparison between the test data and the simulation data of road noise is shown in Figures 10 and 11.

Simulation Results of Road Noise
Based on the acoustic-solid coupled finite element model, the time-domain signals of acceleration measured by engineers on the actual road surface at the steering knuckle of the front and rear suspension are converted into PSD frequency-domain signal files to stimulate the four wheel cores of the vehicle, and the solution response points are the sound pressure level at the driver's right ear and the middle sound pressure level of the rear passenger's seat. The solution frequency range is set as 20-100 Hz. The comparison between the test data and the simulation data of road noise is shown in Figures 10 and 11. The body structure and the sound cavity system are coupled through the vibration of the boundary nodes. In Hypermesh2017.01 software, the coupling between the structure and the sound cavity is realized by using the ACMODF card to search the structure nodes automatically through the nodes on the sound cavity boundary. The acoustic-solid coupling model of the body system is established as shown in Figure 9 (the left door and part of the side circumference are hidden to display the acoustic cavity model).

Simulation Results of Road Noise
Based on the acoustic-solid coupled finite element model, the time-domain signals of acceleration measured by engineers on the actual road surface at the steering knuckle of the front and rear suspension are converted into PSD frequency-domain signal files to stimulate the four wheel cores of the vehicle, and the solution response points are the sound pressure level at the driver's right ear and the middle sound pressure level of the rear passenger's seat. The solution frequency range is set as 20-100 Hz. The comparison between the test data and the simulation data of road noise is shown in Figures 10 and 11.   The results show that the sound pressure level of the driver's right ear and the middle of the passenger's rear seat has an obvious peak value at 35 Hz in the low-frequency band, which is corresponding to the peak frequency of the test results of the real car and the overall trend is basically corresponding. The sound pressure level in the right ear of the real car test driver is 50.39 dB (A), and the simulation data is 49.30 dB (A). The test sound pressure level in the middle of the rear passenger row is 50.79 dB (A), and the simulation data is 50.50 dB (A). In the simulation and test, the difference in the sound pressure level between the driver's right ear and the middle of the rear passenger's row is relatively small at the peak point, indicating that the previous modified finite element model has high accuracy. In view of this kind of steady-state structural noise, the main contribution of the peak frequency can be effectively determined by analyzing the joint contribution of the body panel. The results show that the sound pressure level of the driver's right ear and the middle of the passenger's rear seat has an obvious peak value at 35 Hz in the low-frequency band, which is corresponding to the peak frequency of the test results of the real car and the overall trend is basically corresponding. The sound pressure level in the right ear of the real car test driver is 50.39 dB (A), and the simulation data is 49.30 dB (A). The test sound pressure level in the middle of the rear passenger row is 50.79 dB (A), and the simulation data is 50.50 dB (A). In the simulation and test, the difference in the sound pressure level between the driver's right ear and the middle of the rear passenger's row is relatively small at the peak point, indicating that the previous modified finite element model has high accuracy. In view of this kind of steady-state structural noise, the main contribution of the peak frequency can be effectively determined by analyzing the joint contribution of the body panel.

Analysis of Joint Contribution of Body Parts
The control path of vehicle interior noise is generally: excitation source-transfer path-response point. For mass-produced models, it is very difficult to control the excitation source (chassis-suspension system). Modifying the chassis-suspension system will have a great impact on the smoothness, operation stability and safety performance of the whole vehicle, and the modification cost is huge. Therefore, for the noise problem of body vibration radiation, general enterprises tend to consider it from the perspective of transmission path (body sheet metal). The car body and the car sound cavity can be regarded as a closed cavity. The closed cavity sheet metal S low-frequency vibration radiates into the cavity to form the car's low-frequency roar. Take a point inside the space O as the origin of coordinates, and the rest of the points in the space use it as a reference point. Assume that point J is a point on the closed cavity sheet metal part S, and its radiation sound pressure on a point i in the car is: where ( 0 , ) is the radiation sound pressure from point J to point I, is the volume sound source at point i, ( , ) is the sound pressure at point j, ( , ) is the speed that the sound source j transmits to point i, 0 is the distance from point i to point o, r is the distance from point j to point o, w is the frequency. Many points on the body closed cavity sheet metal parts S radiation noise to point i, and the proportion of point j is:

Analysis of Joint Contribution of Body Parts
The control path of vehicle interior noise is generally: excitation source-transfer pathresponse point. For mass-produced models, it is very difficult to control the excitation source (chassis-suspension system). Modifying the chassis-suspension system will have a great impact on the smoothness, operation stability and safety performance of the whole vehicle, and the modification cost is huge. Therefore, for the noise problem of body vibration radiation, general enterprises tend to consider it from the perspective of transmission path (body sheet metal). The car body and the car sound cavity can be regarded as a closed cavity. The closed cavity sheet metal S low-frequency vibration radiates into the cavity to form the car's low-frequency roar. Take a point inside the space O as the origin of coordinates, and the rest of the points in the space use it as a reference point. Assume that point J is a point on the closed cavity sheet metal part S, and its radiation sound pressure on a point i in the car is: where P ij (r 0 , w) is the radiation sound pressure from point J to point I, Q i is the volume sound source at point i, P j (r, w) is the sound pressure at point j, v i (r, w) is the speed that the sound source j transmits to point i, r 0 is the distance from point i to point o, r is the distance from point j to point o, w is the frequency. Many points on the body closed cavity sheet metal parts S radiation noise to point i, and the proportion of point j is: P(r 0 , w) is the total sound pressure at point i, and t j (w) is the contribution of point j on the body closed cavity sheet metal part S to the interior noise. When all the point contributions are plotted together, the contribution of each area to the interior acoustic radiation can be determined. Through the analysis of the contribution source of the nodes, the main contribution node area can be found, so as to determine the source of the lowfrequency roar contribution. Analysis results of 35 Hz body panel joint contribution were extracted, as shown in Figure 12. The results show that the rear door has the largest contribution to 35 Hz noise.
( 0 , ) is the total sound pressure at point i, and ( ) is the contribution of point j on the body closed cavity sheet metal part S to the interior noise. When all the point contributions are plotted together, the contribution of each area to the interior acoustic radiation can be determined. Through the analysis of the contribution source of the nodes, the main contribution node area can be found, so as to determine the source of the lowfrequency roar contribution. Analysis results of 35 Hz body panel joint contribution were extracted, as shown in Figure 12. The results show that the rear door has the largest contribution to 35 Hz noise.

Simulation Analysis of Vehicle Acoustic Cavity Mode
In the process of driving a car, the sheet metal parts of the car body and other sealing components form a closed cavity, and the acoustic cavity mode of the car reflects the inherent properties of the interior air fluid. When the structure resonates with the acoustic cavity mode, the interior noise will be amplified to form the low-frequency roaring sound (20 Hz-100 Hz). By comparing the result of free mode analysis of the sound cavity with the noise in the car, it can be judged whether the 35 Hz noise in the car is caused by the resonance between the sound cavity and the body panel.
The modal solution card was set up in the cavity finite element model established earlier, and the modal solution file was imported into Nastran14.0 solver for solving. Through calculation, the modal frequency and mode shape of the acoustic cavity of the structure are obtained in Table 2, Figures 13 and 14:

Simulation Analysis of Vehicle Acoustic Cavity Mode
In the process of driving a car, the sheet metal parts of the car body and other sealing components form a closed cavity, and the acoustic cavity mode of the car reflects the inherent properties of the interior air fluid. When the structure resonates with the acoustic cavity mode, the interior noise will be amplified to form the low-frequency roaring sound (20 Hz-100 Hz). By comparing the result of free mode analysis of the sound cavity with the noise in the car, it can be judged whether the 35 Hz noise in the car is caused by the resonance between the sound cavity and the body panel.
The modal solution card was set up in the cavity finite element model established earlier, and the modal solution file was imported into Nastran14.0 solver for solving. Through calculation, the modal frequency and mode shape of the acoustic cavity of the structure are obtained in Table 2, Figures 13 and 14:     Through simulation, the front section has determined that the peak noise of 35 Hz in Figure 14. Model vibration shape of the second-order acoustic cavity.

Establishment of Mathematical Model of Acoustic Metamaterial
Through simulation, the front section has determined that the peak noise of 35 Hz in the car is mainly caused by low-frequency vibration radiation of the tail door. Therefore, according to the corresponding size structure of the plate of the tail door, an acoustic metamaterial with specific frequency suppression can be designed. According to the mathematical model of the tailgate, the installation position of the acoustic metamaterial is determined corresponding to the position of the red box shown in Figure 15 (the rectangular shape of the acoustic metamaterial is 60 × 204 × 3 mm, 1 is the inner panel of the tail door, 2 is the reinforced floor of the tail door, 3 is the outer panel of the tail door, 4 is the car tail door lock, and the red dotted box is the attachment position of the acoustic metamaterial of the tail door). The acoustic metamaterial is composed of three parts: Acoustic metamaterial matrix plate 5, adjustable mass block 6 and acoustic metamaterial damping layer 7, as is shown in Figure 16. The designed acoustic metamaterial is directly attached to the sheet metal parts of the tailgate through the magnetic damping layer, eliminating the need for a complicated installation process. The adjustable mass block and the cantilever beam in the acoustic metamaterial layer constitute a typical acoustic metamaterial spring-mass system, which directionally regulates the band gap frequency of the acoustic metamaterial by adjusting the weight of the mass block.

Frequency Design of Acoustic Metamaterial Resonance Element
The mathematical model of acoustic metamaterial was imported into permesh2017.01, which was processed by solid cutting and placed in the mapping s A hexahedral mesh of 1mm was established for each cut entity. The total number of ac tic metamaterials grid is 40,324. Assign material attributes to each component, as sho in Table 3. After designing the modal solution card, it is imported into Nastran solver solving the calculation. Based on the principle of local resonance of acoustic metamate the resonance element has to produce local resonance with the sheet metal. The nat

Frequency Design of Acoustic Metamaterial Resonance Element
The mathematical model of acoustic metamaterial was imported into Hy-permesh2017.01, which was processed by solid cutting and placed in the mapping state. A hexahedral mesh of 1mm was established for each cut entity. The total number of acoustic metamaterials grid is 40,324. Assign material attributes to each component, as shown in Table 3. After designing the modal solution card, it is imported into Nastran solver for solving the calculation. Based on the principle of local resonance of acoustic metamaterial, the resonance element has to produce local resonance with the sheet metal. The natural

Frequency Design of Acoustic Metamaterial Resonance Element
The mathematical model of acoustic metamaterial was imported into Hypermesh2017.01, which was processed by solid cutting and placed in the mapping state. A hexahedral mesh of 1mm was established for each cut entity. The total number of acoustic metamaterials grid is 40,324. Assign material attributes to each component, as shown in Table 3. After designing the modal solution card, it is imported into Nastran solver for solving the calculation. Based on the principle of local resonance of acoustic metamaterial, the resonance element has to produce local resonance with the sheet metal. The natural frequency of the resonant element is calculated as f = 1 2π k m (k is the equivalent stiffness, m is the equivalent mass). By adjusting the thickness of attached mass block 6 (which is equivalent to changing the value of f by changing m in the calculation formula of natural frequency), the first-order natural frequency of each resonance element in the acoustic metamaterial is within the range of 35 ± 0.5 Hz. The adjusted calculation results are shown in Figure 17.

Results and Discussion
After the corresponding acoustic metamaterial is designed, the acoustic metamaterial finite element model is attached to the design part, and the whole vehicle finite element model with the acoustic metamaterial attached to the tailgate is established, as shown in Figure 18 (in order to see the installation position of the acoustic metamaterial, the attached sheet metal parts is displayed transparently). The response points are the sound pressure level in the driver's right ear, the middle pressure level in the rear passenger's seat and the vibration acceleration in the middle of the tailgate panel.

Results and Discussion
After the corresponding acoustic metamaterial is designed, the acoustic metamaterial finite element model is attached to the design part, and the whole vehicle finite element model with the acoustic metamaterial attached to the tailgate is established, as shown in Figure 18 (in order to see the installation position of the acoustic metamaterial, the attached sheet metal parts is displayed transparently). The response points are the sound pressure level in the driver's right ear, the middle pressure level in the rear passenger's seat and the vibration acceleration in the middle of the tailgate panel. The same road noise excitation signal and excitation point positions as in Section 3 were used to conduct road noise simulation analysis and vibration response analysis of the plate metal parts of the tailgate before and after the attaching of the acoustic-solid coupled finite element model. The simulation results are shown in Figures 19-22. By comparing the data before and after the acoustic metamaterial, the vibration and noise reduction effect of the acoustic metamaterial was verified.  The same road noise excitation signal and excitation point positions as in Section 3 were used to conduct road noise simulation analysis and vibration response analysis of the plate metal parts of the tailgate before and after the attaching of the acoustic-solid coupled finite element model. The simulation results are shown in Figures 19-22. By comparing the data before and after the acoustic metamaterial, the vibration and noise reduction effect of the acoustic metamaterial was verified.    In the original state without an acoustic metamaterial attached, the sound pressure level at the driver's right ear and the middle seat of the rear passenger's seat has an obvious peak value at 35 Hz; the amplitude of the sound pressure level of the driver's right      By comparing the data verified by the above simulation, it is found that the local resonant acoustic metamaterial has an ideal effect on the control of low-frequency noise radiated from the vibration of sheet metal parts to the vehicle.
Regarding the durability of the device, components of acoustic metamaterial: aluminum plate and the damping layer have good durability and have been widely used in the automotive industry, its durability can meet the needs of mass production applications. In addition to that, the three-interval method of random vibration fatigue life calculation based on Gaussian distribution and Miner's linear cumulative damage criterion shows that the structure can carry out infinite cycles under the random load, which meets the engineering requirements.

Conclusions
Focusing on the problem of low-frequency noise in a certain car, this paper determines the frequency point of an abnormal peak value of noise in the car through test and simulation analysis. Through CAE simulation analysis, it is concluded that the low-frequency noise is caused by the low-frequency vibration of sheet metal parts of the tailgate. To solve this problem, an acoustic metamaterial is designed based on the band gap theory of local resonant phonon crystal. Through simulation verification, it is judged that the structure can indeed reduce the amplitude of interior noise and vibration acceleration, which verifies the effectiveness of the acoustic metamaterial. It provides a new way of thinking about the problem that the low-frequency noise of vibration and radiation of sheet metal parts is difficult to be controlled.

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