Application of wave based method for predicting the response of coupled vibro-acoustic system with unconstrained damping layer

Xiaojun Xia1 , Zhongming Xu2 , Shiyang Lai3 , Zhifei Zhang4 , Yansong He5

1, 2State Key Laboratory of Mechanical Transmission, Chongqing, China

2, 3, 4, 5College of Automotive Engineering of Chongqing University, Chongqing, China

2Corresponding author

Vibroengineering PROCEDIA, Vol. 5, 2015, p. 27-32.
Accepted 21 August 2015; published 18 September 2015

Copyright © 2015 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Abstract.

The Wave Based Method (WBM) is a deterministic prediction method that is computational efficiency as compared to other deterministic prediction techniques in mid-frequency problems. This paper discusses the application of WBM for predicting the dynamic displacement of plate with an unconstrained damping layer based on Kirchhoff theory. Further, the prediction of acoustic response of the coupled vibro-acoustic system with unconstrained damping is realized on the use of WBM. A numerical example is introduced, and the comparison of numerical result obtained by WBM and FEM is acquired. It is seen that the WBM is applicable for vibro-acoustic system with unconstrained damping and is expected to yield faster and more accurate predictions. The limitation of the method caused by simplify hypothesis is described in combination with modelling ways and numerical results.

Keywords: wave based method, vibro-acoustic, unconstrained damping, Kirchhoff theory.

1. Introduction

The element based methods such as FEM and BEM are limited for low-frequency problems because of its high computational cost and increasing error as frequency increasing [1]. While the statistical methods like SEA [2] are merely suit for high frequency problems. So an efficient numerical technique for simulating and predicting the mid-frequency problems is desired. Wave based method (WBM) as a deterministic method this paper focuses on is a better choice. At present, WBM has shown its high accuracy and low compute cost in contrast with existing deterministic method when predicting the vibro-acoustic response in mid-frequency [3]. Meanwhile, as one of the most common and convenient methods for vibration and noise reduction, adding an unconstrained damping layer on plate like structure is wildly adopted in engineering. The analytic method to analysis the dynamic displacement of composite is proposed based on the Kirchhoff theory [4, 5].The application of such efficient numerical predictive technique for vibro-acoustic system with unconstrained damping is indispensable and valuable.

This paper starts with the description of dynamic vibration of a plate with unconstrained damping based on the Kirchhoff theory and the description of government equations of coupled vibro-acoustic system. The following section describes the basic principle of WBM for coupled vibro-acoustic system with unconstrained damping layer. A numerical example demonstrates the validity and the efficiency of the presented method finally.

2. Basic methodology of coupled vibro-acoustic system

2.1. Vibration of plate with an unconstrained damping layer

Fig. 1 shows the section of composite plate. t1, t2 is the thickness of steel layer and damping layer respectively. The dynamic displacement government equation based on Kirchhoff theory is:

(1)
B * 4 w t x s , y s + m ~ 2 w t x s , y s t 2 = Q t .

With m~=ρ1t1+ρ2t2 the unit mass of plate, Qt the external load. B* represents the complex bending stiffness under the Kirchhoff hypothesis, which defined as:

(2)
B * = E 1 * t 1 3 12 1 - u 1 2 + E 2 * 1 - u 2 2 t 2 3 12 + t 2 + t 21 2 ,

where t21=t1+t2/2, E1*=E11+jη1, E2*=E21+jη2, e2=E2/E1, h2=t2/t1. Ei, ηi, i=1, 2 is elasticity modulus and material loss factor of each layer correspondingly.

Fig. 1. Structure of a plate with unconstrained damping

 Structure of a plate with unconstrained damping

2.2. Government equation of vibro-acoustic system

For a steady-state harmonic excitation Qt=Qejωt, displacement and the excitation can be written as wzxs,ysejωt, subsequently, the steady government equation of the composite plate of vibro-acoustic system is redefined as:

(3)
4 w z x s , y s - k b * 4 w z x s , y s = Q B * δ x F s , y F s + p x a s , y a s , z a s B * ,

where 4=4/x4+4/x2y2+4/y4,kb*=m~ω2/B*4 is the plate bending wave number, xFs,yFs is the location of external load. pxas,yas,zas/B* describes the acoustic load acting on the coupled interface xas,yas,zas.

The Helmholtz equation of the steady pressure of associated cavity is:

(4)
2 p r + k 2 p r = 0 ,

with 22/x2+2/y2+2/z2 the Laplace operator, r the coordinate in cavity, k the acoustic wavenumbers.

3. Wave based method

The WBM belongs to the category of indirect Trefftz methods. As the Trefftz principle, the field variables within convex domain are approximated by an expansion of basis functions, which exactly satisfy the governing dynamic equations. The errors on the boundary is forced to zero through weighted residual formulation, then the contribution factor of each basis function are determined. Subsequently, the displacement and acoustic response of vibro-acoustic system are obtained.

3.1. Field variable expansion

The variables, the out of plane displacement wz of composite plate and the pressure p of cavity, are approximated by the following field variable expansion:

(5)
w z i = 1 n s w s i ψ s i + i = 1 n a p a i w s i + w ^ F = ψ s w s + p a φ a s + w ^ F ,
(6)
p i = 1 n a p a i ϕ a i = ϕ a p a ,

With ψs, ϕa the wave functions vectors, ws, pa the corresponding contribution factors vectors of displacement and pressure respectively. And the wave function defined as:

(7)
ψ s i = e - j k x s i x + k y s i y ,       ϕ a i = e - j k x a i x + k y a i y + k z a i z .

With kxsi2+kysi22=kb*4, and kxai2+kyai2+kzai2=k2.

w ^ F serves as the particular solution for external point force excitation:

(8)
w ^ F = - j F 8 k b * 2 B * H 0 2 k b * r F - H 0 2 - j k b * r F ,

where rF is the distance to the location of exciting point xF,yF.

ϕ a s is the pressure load results from the acoustic wave functions ϕa, defined as:

(9)
φ a s ( x s , y s ) = - j 8 k 2 * 2 B * Ω s   ϕ a x a s , y a s , z a s H 0 2 ( k b * r s   ) d s
            + j 8 k 2 * 2 B * Ω s   ϕ a x a s , y a s , z a s H 0 2 ( - j k b * r s   ) d s .

With rs the distance to xs,ys.

3.2. Coupled vibro-acoustic wave model

The wave functions of plate and acoustics satisfy the corresponding government equations. The boundary residual is enforced to zero through weighted formula, like the Galerkin weighting procedure used in FEM. This yields a matrix equation consisting of ns+na algebraic equations in the ns+na unknown wave function contribution factors:

(10)
A a a + C a a C a s C s a A s s p a w s = f a f s ,
(11)
A a a = Γ v   ϕ a T L v ϕ a d Γ ,
(12)
C a a = - Γ s   j ω ϕ a T ϕ a s d Γ ,
(13)
C a s = - Γ s   j ω ϕ a T ψ s d Γ ,
(14)
C s a = Γ w θ   L Q n ψ s T w a s d Γ + Γ w θ   L m n ψ s T L θ n w a s d Γ ,
(15)
A s s = Γ w θ   L Q n ψ s T ψ s d Γ + Γ w θ   L m n ψ s T L θ n ψ s d Γ ,
(16)
f a = Γ v   ϕ a T v - n d Γ - Γ s   ϕ a T j ω w F d Γ ,
(17)
f s = Γ w θ   L Q n ψ s T w - z - w ^ F d Γ + Γ w θ   L m n ψ S T θ - n - L θ n w ^ F d Γ ,

where Γv, Γwθ, Γs is the acoustic boundary, plate boundary and coupled interface respectively. Lv, LQn, Lθn and Lmn is the differential operators for acoustic velocity, the normal rotation of plate, the bending moment of plate and the generalized shear force of plate. After solving the matrix equations, the obtained contributions pa and ws are substituted back to Eqs. (5) and (6) to acquire the variables p and wz.

4. Numerical example

4.1. Problem description

In order to validate the described methodology and show its capabilities a numerical example is given. The considered geometry of coupled vibro-acoustic system is shown in Fig. 2. The parameters of composite plate are presented in Table 1. A unit point force is applied at (0.2 m, 0.2 m, 0.5 m) and a reference point R(0.4 m, 0.3 m, 0.5 m) in cavity is selected. WBM model is built in MATLAB R2015a, and the FEM model is built with MSC/Nastran.

Fig. 2. Problem geometry

 Problem geometry

Table 1. Parameters of composite plate

Material
Thickness (mm)
Density (kg/m3)
Elasticity modulus (pa)
Poisson coefficient
Loss factor
Aluminum
1.0
2700
7.2e10
0.33
0
Damper
2.0
1300
1e8
0.45
0.8

4.2. Result

As we know, the accuracy of FEM is improving with decreasing the size of element to some degree, which, however, increases the computational cost and limits its capacity for higher frequency problems. To illustrate the influence of element size, coarse FE model with 20 element size and fine FE model with 5 element size are built. Fig. 3 plots the acoustic response of reference point calculated by FE and WBM at 50-500 Hz. The comparison shows that the result obtained by WBM agrees with the result obtained by fine FE model better than the coarse one. The influence of damping layer is revealed through the comparison with the undamping model, and the average amplitude of acoustic response is reduced dramatically compared with the damping model.

Fig. 4 shows the displacement contour of composite plate calculated by WBM with 262 wave functions and FEM with fine mesh at 100 Hz. It is seen that the displacement is zero at the clamped edge that satisfy the boundary condition. And the result obtained by WBM can agree with the result obtained by fine FE model basically.

The error generated by this methodology is mainly caused by: (1) WBM makes use of the Kirchhoff plate theory to model plate bending problems, whereas the FEM in MSC/Nastran is based on the Reissner-Mindlin theory that considers the rotatory inertia and shear deformation. (2) The composite plate is treated as a single plate in WBM, while the damping layer is built by solid element in FEM. So the thicker damping layer, the influence more serious, which restricts that this method is applied only to the thin plate that according with the Kirchhoff hypothesis.

Fig. 3. The acoustic response curves of system: a) with and b) without unconstrained damping

The acoustic response curves of system: a) with and b) without unconstrained damping

a)

The acoustic response curves of system: a) with and b) without unconstrained damping

b)

In contrast with the FEM, the advantages of the application using WBM for predicting coupled system with unconstrained damping is obvious and significant.it has a higher convergence rate and convenient modeling produce. Table 2 shows convergence of WBM and FEM. The comparison demonstrates that WBM can get the acceptable result with much less computational cost than FEM especially for mid-frequency problems.

Fig. 4. Displacement contour of composite plate at 100 Hz: a) WBM and b) FEM

Displacement contour of composite plate at 100 Hz: a) WBM and b) FEM

a)

Displacement contour of composite plate at 100 Hz: a) WBM and b) FEM

b)

Table 2. The DOF and corresponding pressure amplitude of reference point at 60 Hz

WBM
FEM
Pressure (Pa)
Wave functions
Pressure (Pa)
Nodes
2.6053
132
2.7659
10925
2.5088
164
2.5543
20834
2.4368
198
2.4768
45758
2.4217
230
2.3784
148718
2.4217
258
2.3768
1134133

5. Conclusion and future work

This paper describes the method to predict vibration of plate with an unconstrained damping layer through the complex stiffness. Subsequently, the methodology for predicting the coupled vibro-acoustic system with unconstrained damping is presented. The last section gives a numerical example, which compares the result obtained by WBM to the results obtained by FEM with different sizes. In the example, the WBM shows a higher convergence rate and lower computational cost in contrast with FEM. The comparison validates the effectiveness of this method for predicting coupled vibro-acoustic system with unconstrained damping layer.

An interesting next step is to explore the availability for acoustic optimization based on this method. Furthermore, an experimental validation of the obtained numerical results is foreseen.

Acknowledgements

The research work of Xia Xiaojun is financed by Chongqing Graduate Student Research Innovation Project (CYB14036).

References

  1. Hartmann F., Zotemantel R. The direct boundary element method in plate bending. International Journal for Numerical Methods in Engineering, Vol. 21, Issue 11, 1986, p. 2049-2069. [Search CrossRef]
  2. Lyon R., DeJong R. Theory and Application of Statistical Energy Analysis (2nd Edition). Butterworth Heinemann, 1995. [Search CrossRef]
  3. Desmet W. A Wave Based Prediction Technique for Coupled Vibro-Acoustic Analysis. Leuven, 1998. [Search CrossRef]
  4. Zheng Baizhe, Li Yahong, Zhu Yuxiang Theoretical and experimental research of the vibration characteristics of the damped plate composed of an elastic metal layer and an unconstrained viscoelastic layer. Engineering Mechanics, Vol. 8, Issue 3, 1991, p. 129-138. [Search CrossRef]
  5. Leissa A. Vibration of Plates. Acoustical Society of America, Woodbury, New York, 1993. [Search CrossRef]
  6. Vanmaele C., Vandepitte D. An efficient wave based prediction technique for plate bending vibrations. Computer Methods in Applied Mechanics and Engineering, Vol. 196, Issues 33-34, 2007, p. 3178-3189. [Search CrossRef]
  7. luymers B., et al. Trefftz-based methods for time-harmonic acoustics. Archives of Computational Methods in Engineering, Vol. 14, Issue 4, 2007, p. 343-381. [Search CrossRef]
  8. Vergote K., et al. An efficient wave based approach for the time-harmonic vibration analysis of 3D plate assemblies. Journal of Sound and Vibration, Vol. 332, Issue 8, 2013, p. 1930-1946. [Search CrossRef]