Published: 31 December 2013

Research on active sound absorption based on minimality of quadratic sum of reflected sound pressure in low frequency

Congyun Zhu1
Jianru Shi2
Shufeng Yang3
Qibai Huang4
1, 2, 3School of Mechanical Science & Engineering, Zhongyuan University of Technology, Zhengzhou, Henan 450007, P.R.China
4School of Mechanical Science & Engineering Hua Zhong University of Science and Technology, Wuhan, Huibei 430074, P.R.China
Corresponding Author:
Congyun Zhu
Views 29
Reads 13
Downloads 1127

Abstract

In this paper, the minimality of the quadratic sum of the reflected sound pressure is used to outline the method of active sound absorption. A piezoelectric ceramic is attached to a simply supported plate, and two microphones are placed in the front of the simply supported plate. The reflected sound pressure is measured. According to the measured reflected sound pressure, the voltage is added to the surface of the piezoelectric ceramic. Then, the simply supported plate is vibrated. Thus, the total sound pressure is expressed as a combination of the reflected sound pressure and the radiated sound pressure. A condition of the quadratic sum being minimal is presented and the voltage (added to the surface of the piezoelectric ceramic) is calculated. At last, the numerical calculation and the experiment are carried out, which demonstrate that the method outlined in this paper is possible.

1. Introduction

In 1960s, active sound elimination raised a great interest amongst scholars. Active sound elimination is a method that uses primary sound source sound waves in order to produce equal and opposite in phase waves from a secondary sound source (destructive interference between the two sources results in sound elimination).

The expected control goals are achieved and the purpose of the elimination is attained. The method was proposed by Paul Lueg, who patented it in 1934 and explained its basic theory in 1936 [1-3]. The electronic control system could not be manufactured due to electronic technology imitations. Therefore, the technology had not been developed until the late 1960s. With the development of electronic technology, scholars posed the method of sound elumination again. The present study hotspot is using structural acoustic radiation to offset the incidental acoustic energy in order improve the absorption effect. There have been a lot of reseach achievements and some research results have been used in practice.

The basic principle of the elimination is shown in Fig. 1, when the sound wave is incident to the flat plate. The sound signal is detected by the sensor placed in front of the plate. The detected sound signal is transmitted to the control filters, which give the controlled signal to the actuator attached to the plate (the plate is assumed to be elastic). When the elastic plate has no secondary force source effect, the acoustic field above the plate includes: incidental sound pressure pi(r,t), the reflected sound pressure pr(r,t) by the elastic plate and the acoustic radiation pressure pps(r,t). These elements form the following formula [4-6].

1
ppr,t=pir,t+prr,t+ppsr,t.

Fig. 1Active noise elimination is using a plate structure

Active noise elimination is using a plate structure

When the secondary force source effects the elastic plate, the radiated sound pressure is pcs(r,t). So the total sound pressure is:

2
pr,t=ppr,t+pcsr,t=pir,t+prr,t+ppsr,t+pcsr,t.

Similarly, when the elastic plate has no secondary force source, the velocity in the space above the elastic plate includes three parts: (1) he velocity Vi(r,t) generated by the incident sound wave; (2) the reflected velocity Vr(r,t) generated by the elastic plate; (3) the radiated velocity Vps(r,t), which is shown as the following formula:

3
Vpr,t=Vir,t+Vrr,t+Vpsr,t.

When the secondary force source effects the elastic plate, the radiated velocity Vcs(r,t) is generated and the total velocity is shown as the following formula:

4
Vr,t=Vpr,t+Vcsr,t=Vir,t+Vrr,t+Vpsr,t+Vcsr,t.

Therefore, when the secondary force source effects on the elastic plate, the overall sound power of the elastic plate W is:

5
W=pr,tVr,t.

In order to achieve the purpose of the activity eliminated, the overall sound power W should be minimal and the control target is:

6
Wmin.

The incident sound pressure is constant, which has no affection on the derivation. the minimality of the quadratic sum of the reflected sound pressure is used to outline the method of active sound absorption. A piezoelectric ceramic is attached to a simply supported plate, and two microphones are placed in the front of the simply supported plate. The reflected sound pressure is measured. According to the measured reflected sound pressure, the voltage is added to the surface of the piezoelectric ceramic. Then, the simply supported plate is vibrated. Thus, the total sound pressure is expressed as a combination of the reflected sound pressure and the radiated sound pressure. A condition of the quadratic sum being minimal is presented and the voltage (added to the surface of the piezoelectric ceramic) is calculated.

2. The theory calculation method of the quadratic sum of the reflected sound pressure with monolithic piezoelectric ceramic chips

The arrangement of the monolithic piezoelectric ceramic chip is shown in the Fig. 2.

Fig. 2Active sound absorption is based on the minimality of quadratic sum of the reflected sound pressure

Active sound absorption is based on the minimality of quadratic sum of the reflected sound pressure

For the external voltage Vi, the strain δi of the ith piezoelectric ceramic chips is:

7
δi=ΔLiLi=d31Vitia,

where tia is the thickness of the ith piezoelectric ceramic chip and d31 is the strain constant of the piezoelectric ceramic chip. In the Fig. 2, the piezoelectric ceramic chips are attached to the simply supported plate and the double direction longitudinal force reads [7-15]:

8
Mix=Miy=ρia(2+ρia)4(1+βρia(3+3ρia+ρia2))h2γδi,
9
ρia=2tiah.

According to the formula (8), the motion equation of the simply-supported plate is:

10
D4wi+ρhw¨i=Mixδx-xi1x-δx-xi2xuy-yi1-uy-yi2
+ δy-yi1y-δy-yi2yux-xi1-ux-xi2.

For the formula (10), the displacement wi(x,y,t) of the simply-supported plate at any location and time can be expressed using the vibration modal method:

11
wix,y,t=ejωtm=1n=1Aimnψimnx,y,

where Aimn is (m,n) modal amplitude caused by the i piece of piezoelectric ceramic chip and ψimn(x,y) is the vibration modal shape function of the position (x,y):

12
ψimnx,y=sinmπxLxsinnπyLy.

Putting the formula (11) and (12) into the formula (10) and then multiplying ψimn(x,y) on both sides of formula (10) yields:

13
0Lx0LysinmπxLxsinnπyLym=1n=1Aimnx,ysinmπxLxsinnπyLyDmπLx2+nπLy2-ρhω2dxdy.

When the model is mutually orthogonal, the integral of any two different vibration models along the surface of simply supported plate is equal to zero. Therefore, the equation (13) reads:

14
LxLy2AimnDmπLx2+nπLy2-ρhω2.

According to the formula (14), the inherent frequency ωmn is:

15
ωmn2=DmπLx2+nπLy2ρh.

ψimn(x,y) is multiplied by the right of formula (9) and yields:

16
Mix0Lx0LysinmπxLxsinnπyLyδx-xi1x-δx-xi2xuy-yi1-uy-yi2+δy-yi1y-δy-yi2yux-xi1-ux-xi2dxdy.

According to the formula (16), the following formula yields:

17
-2MixLxLymπLx2+nπLy2mnπ2cosmπxi1Lx-cosmπxi2Lxcosnπyi1Ly-cosnπyi2Ly.

To make formula (14) and formula (17) equal, the modal amplitude Aimn is expressed as:

18
Aimn=4MixmπLx2+nπLy2ρhmnπ2ωmn2-ωi2cosmπxi2Lx-cosmπxi1Lxcosnπyi1Ly-cosnπyi2Ly.

According to the above formula, the vibration displacement wi(x,y,t) of the simply supported plate can be yield (which is differential). The vibration velocity u(x,y,t) of the simply supported plate’s surface is:

19
ux,y,t=i=1Njωiejωitm=1n=1Aimnψimnx,y,

where N is the number of piezoelectric ceramic chips that are attached to the simply supported plate. According to the above analysis and the calculation method, the sound pressure pi(xQ,yQ) is expressed as matrix:

20
Pi=pk=1,l=1i,,pk=1,l=ai,pk=2,l=1i,pk=2,l=ai,,pk=b,l=1i,pk=b,l=2i,,pk=b,l=ai.

The radiated sound pressure of the simply supported plate caused by N piezoelectric ceramic chip is:

21
P=i=1NPi=i=1Npk=1,l=1i,,i=1Npk=1,l=ai,i=1Npk=2,l=1i,
i=1Npk=2,l=ai,,i=1Npk=b,l=1i,,i=1Npk=b,l=ai.

The simple supported plate is divided into a×b units as showed in the Fig. 3.

Fig. 3Element division of the supported plate

Element division of the supported plate

According to the Fig. 3, the center coordinations of the No. Mthelement and the No. Qthelement are:

22
x:Coordinatei-1aLx,iaLx,y:Coordinatej-1bLy,jbLy,x:Coordinatek-1aLx,kaLx,y:Coordinatel-1bLy,lbLy,

where i and j are the orders directing x and yof the No. Mth rectangular element in the Fig. 3. k and l is the order directing x and y of the No. Qth rectangular element. Due to the elements being tiny, the coordinate of the No. Mth rectangular element is i-12aLx,j-12bLy and the coordinate of the No. Qth rectangular element is k-12aLx,l-12bLy. The distance between the No. Mth rectangular element and No. Qth rectangular element is:

23
rM,Q=i-kaLx2+j-lbLy2.

Due to the elements being tiny, the surface vibration velocity of every element is equal, thus the velocity of the No. Mth rectangular element is:

24
uMx,y,t=jωejωtm=1n=1Amnψmni-12aLx,j-12bLy.

Accordingly, formula (23) and formula (24) are put into the formula (23):

25
pxQ,yQ=jωρ0ΔS2πk=1al=1buxM,yMe-jkrr,

where ρ0 is the density of the air, ΔS is the area of each element, ΔS=LxLyab. The formula (25) can be denoted as the following formula:

26
P=pk=1,l=1,pk=1,l=2,,pk=1,l=a,pk=2,l=1,pk=2,l=2,,pk=2,l=a,
, pk=b,l=1,pk=b,l=2,,pk=b,l=a.

According to the above mentioned, when there are no secondary force effects on the simple supported plate, the reflected sound consists of two parts: (1) geometric reflected waves pr(r,t) (in the condition that the elastic is rigid plate) and (2) the acoustic radiation wave pps(r,t):

27
pr,t=prr,t+ppsr,t.

When the piezoelectric ceramic chip attached to the simply supported plate has voltage, the radiation sound pressure is p. The total reflected sound pressure is:

28
pr,t=prr,t+ppsr,t+p.

Assuming the incident sound wave doesn't work, (the incident sound pressure pps(r,t) is zero) the formula (28) is:

29
pr,t=prr,t+p.

The incident sound wave is a plane wave, therefore the reflected sound pressure can be shown as the following formula:

30
Pr=[pr,pr,pr]1×(a×b).

Put the formula (26) and the formula (30) into the formula (29):

31
Pr,t=pr+pk=1,l=1,pr+pk=1,l=2,,pr+pk=b,l=1,pr+pk=b,l=2,,pr+pk=b,l=a.

According to the formula (31), the standards of the method mentioned in this paper are attained. The corresponding target function is:

32
Jp=PHr,tPr,t.

Jp is the function of the voltage and the angular frequency of the piezoelectric ceramic chip. If the value of the formula (32) is minimum, the following formula is:

33
JpV1=PH(r,t)P(r,t)V1=0,Jpω1=PH(r,t)P(r,t)ω1=0,JpVi=PH(r,t)P(r,t)Vi=0,Jpωi=PH(r,t)P(r,t)ωi=0,JpVN=PH(r,t)P(r,t)VN=0,JpωN=PH(r,t)P(r,t)ωN=0.

The calculated voltage and frequency (according to the above formula) is the voltage and frequency of the piezoelectric ceramic chip.

3. Measurement theory of the incident and reflection wave

If the wave is a plane sound wave, the incidence and reflection sound waves can be expressed as:

34
pix,t=piexpjωt-kx,
35
prx,t=prexpjωt+kx,

where k is the wave number and ω is the frequency of the sound wave.

The signals of microphone 1 and microphone 2 are:

36
p1t=pi0,t+pr0,t,
37
p2t=pid,t+prd,t=pi0,t+τ+pr0,t-τ,

where τ is the propagation time between microphone 1 and microphone 2, τ=dc. c is the propagation velocity of the wave in the air. If the sound pressure p1(t) is delayed by τ, then the following equation yields:

38
p1τt=p1t-τ=pi0,t-τ+pr0,t-τ.

Signal p3(t) is defined as:

39
p3t=p1τ-p2=pi0,t-τ-pi0,t+τ.

If the signal p3(t) is delayed by τ, then the signal p5(t) is:

40
p5t=p3t-τ=pi0,t-2τ-pi0,t=2pi0,te-jπ2+ωτsinωτ.

If the signal p2(t) is delayed by τ, then the following equation yields:

41
p2τt=p2t-τ=pi0,t+pr0,t-2τ.

Signal p4(t) is defined as:

42
p4t=p2τ-p1=pr0,t-2τ-pr0,t=2pr0,te-jπ2+ωτsinωτ.

According to the Eq. 41 and Eq. 42, the following equation can be attained:

43
p4'=-p4+p4'e-2jωτ=pr0,t,
44
p5'=-p5+p5'e-2jωτ=pi0,t.

The signal pr(0,t) is the reflective sound pressure of microphone 1. The reflective sound pressure in the place of the front piezoelectric ceramic is:

45
pr=pr2d,t=expj2kdpr0,t.

The principle of measuring the incidence and reflection sound pressure by two microphones is shown in Fig. 4.

Fig. 4The principle of measuring the incidence and reflection sound pressure

The principle of measuring the incidence and reflection sound pressure

4. Numerical calculation and experiment

The experimental diagram is shown as the Fig. 5.

The rectangular tube for the experiment is 400 mm×400 mm×3000 mm. The primary source is a speaker, which is located at the end of the tube. The distance between the two microphone sensors is 5 cm and the distance between microphone 2 and the simply supported plate is 2.5 cm. The density of the air is ρ=1.21 kg/m3. The propagation velocity in the air is c=343 ms-1. The density of the simply supported plate is ρp=7860 Kg/m3; the thickness is h=0.5 mm; elastic modulus is Ep=1011 Pa; the Poisson’s ratio is νp=0.3. The thickness of piezoelectric ceramic chip is ta=0.5 mm; the piezoelectric constant of the piezoelectric composites is d31=1.66×10-10 m/V; elastic modulus is Ea=6.3×1010 Pa; Poisson’s ratio is νp=0.35.

Fig. 5The experiment arrangement

The experiment arrangement

Fig. 6Layout diagram of piezoelectric ceramic chip

Layout diagram of piezoelectric ceramic chip

a)

Layout diagram of piezoelectric ceramic chip

b)

Layout diagram of piezoelectric ceramic chip

c)

In this experiment, piezoelectric ceramic chips are divided into: single chip piezoelectric ceramic chips; three piezoelectric ceramic chips; five piezoelectric ceramic chips. Fig. 6 is the various layout diagram.

The sound wave pressure figure is following one piece of the piezoelectric ceramic chip. In order to simplify, in this experiment, the incident sound wave is a sine wave. The amplitude is 1.0 Pa and the frequency is 1000 Hz. The incident sound wave and the uncontrolled sound waves of microphone 1 and microphone 2 are shown in Fig. 7-9.

Fig. 7The incident sound wave of amplitude is 1.0 Pa and frequency is 1000 Hz

The incident sound wave of amplitude is 1.0 Pa and frequency is 1000 Hz

Fig. 8The uncontrolled sound waves of microphone 1 and microphone 2

The uncontrolled sound waves of microphone 1 and microphone 2

a)

The uncontrolled sound waves of microphone 1 and microphone 2

b)

Fig. 9The controlled sound waves of microphone 1 and microphone 2

The controlled sound waves of microphone 1 and microphone 2

a)

The controlled sound waves of microphone 1 and microphone 2

b)

The amplitude is 1.0 Pa and the frequency is 2000 Hz. The incident sound wave, the uncontrolled sound waves of microphone 1 and microphone 2 are shown in Fig. 10-12.

Fig. 10The incident sound wave of amplitude is 1.0 Pa and frequency is 2000 Hz

The incident sound wave of amplitude is 1.0 Pa and frequency is 2000 Hz

Fig. 11The uncontrolled sound wave of microphone 1 and microphone 2

The uncontrolled sound wave of microphone 1 and microphone 2

a)

The uncontrolled sound wave of microphone 1 and microphone 2

b)

Fig. 12The controlled sound wave of microphone 1 and microphone 2

The controlled sound wave of microphone 1 and microphone 2

a)

The controlled sound wave of microphone 1 and microphone 2

b)

The amplitude is 1.0 Pa and the frequency is 3000 Hz, The incident sound wave, uncontrolled sound wave of microphone 1 and microphone 2 are shown as Fig. 13-15.

Fig. 13The incident sound wave of amplitude is 1.0 Pa and frequency is 3000 Hz

The incident sound wave of amplitude is 1.0 Pa and frequency is 3000 Hz

Fig. 14The uncontrolled sound wave of microphone 1 and microphone 2

The uncontrolled sound wave of microphone 1 and microphone 2

a)

The uncontrolled sound wave of microphone 1 and microphone 2

b)

Fig. 15The controlled sound wave of microphone 1 and microphone 2

The controlled sound wave of microphone 1 and microphone 2

a)

The controlled sound wave of microphone 1 and microphone 2

b)

The amplitude is 1.0 Pa and the frequency is 4000 Hz. The incident sound wave, uncontrolled sound waves of microphone 1 and microphone 2 are shown in Fig. 16-18.

Fig. 16The incident sound wave of amplitude is 1.0 Pa and frequency is 4000 Hz

The incident sound wave of amplitude is 1.0 Pa and frequency is 4000 Hz

Fig. 17The uncontrolled sound wave of microphone 1 and microphone 2

The uncontrolled sound wave of microphone 1 and microphone 2

a)

The uncontrolled sound wave of microphone 1 and microphone 2

b)

Fig. 18The controlled sound wave of microphone 1 and microphone 2

The controlled sound wave of microphone 1 and microphone 2

a)

The controlled sound wave of microphone 1 and microphone 2

b)

The amplitude is 1.0 Pa and the frequency is 5000 Hz. The incident sound wave, uncontrolled sound waves of microphone 1 and microphone 2 are shown in Fig. 19-21

Fig. 19The incident sound wave of amplitude is 1.0 Pa and frequency is 5000 Hz

The incident sound wave of amplitude is 1.0 Pa and frequency is 5000 Hz

Fig. 20The uncontrolled sound wave of microphone 1 and microphone 2

The uncontrolled sound wave of microphone 1 and microphone 2

a)

The uncontrolled sound wave of microphone 1 and microphone 2

b)

Fig. 21The controlled sound wave of microphone 1 and microphone 2

The controlled sound wave of microphone 1 and microphone 2

a)

The controlled sound wave of microphone 1 and microphone 2

b)

For one piece of the piezoelectric ceramic chip, the sqaure sums of sound before controlling (voltage not applied) and after controlling (voltage applied) are shown in the Fig. 22.

Fig. 22The sound pressure squares of a piezoelectric ceramic chip before and after controlling

The sound pressure squares of a piezoelectric ceramic chip before and after controlling

For three pieces of the piezoelectric ceramic chips, the sqaure sums of sound before controlling (voltage not applied) and after controlling (voltage applied) are shown in the Fig. 23.

Fig. 23The sound pressure squares of three piezoelectric ceramic chips before and after controlling

The sound pressure squares of three piezoelectric ceramic chips before and after controlling

For five pieces of the piezoelectric ceramic chips, the sqaure sums of sound before controlling (voltage not applied) and after controlling (voltage applied) are shown in the Fig. 24.

Fig. 24The sound pressure squares of five piezoelectric ceramic chips before and after controlling

The sound pressure squares of five piezoelectric ceramic chips before and after controlling

From Fig. 7 to Fig. 24, the experimental results obtained the absorption effect. And yet the absorption effect is different from the amount of the piezoelectric ceramic chips. When the amount of piezoelectric ceramic chips is three pieces, the sound absorption effect is superb. If the amount is increased, the absorption is not obvious.

5. Conclusions

Active sound elimination raised a great interest amongst the scholars. Active sound elimination is a method that uses primary sound source sound waves in order to produce equal and opposite in phase waves from a secondary sound source (destructive interference between the two sources results in sound elimination). The expected control goals are achieved and the purpose of the elimination is attained. This effect is unsatisfactory, so this paper puts forward an absorption method based on the minimalism of the reflections power. The minimality of quadratic sum of the reflected sound pressure is used to outline the method of active sound absorption. A piezoelectric ceramic is attached to a simply supported plate, and two microphones are placed in the front of the simply supported plate. The reflected sound pressure is measured. According to the measured reflected sound pressure, the voltage is added to the surface of the piezoelectric ceramic. Then, the simply supported plate is vibrated. Thus, the total sound pressure is expressed as a combination of the reflected sound pressure and the radiated sound pressure. A condition of the quadratic sum being minimal is presented and the voltage (added to the surface of the piezoelectric ceramic) is calculated. In this paper, all measured signals are harmonic. In order to conveniently calculate the equations and finish the experiment, the excited force is a harmonic sound wave. Since the excited force is harmonic sound wave, the microphones signals are harmonic. According to the related equations, the sound signals can be transformed to the sound pressure signals and then sound pressure squares can be calculated. Of course, if the excited force is random, the equation and the conclusion are correct too.

References

  • Kenneth A. Cunefare, Gary H. Koopmann. Global optimum active noise control: Surface and far field effects. Journal of the Acoustical Society of America, Vol. 90, Issue 1, 1991, p. 365-373.
  • Kenneth A. Cunefare, Gary. H. Koopmann. Acoustic design sensitivity for structural radiators. ASME Journal of Vibration and Acoustics, Vol. 114, Issue 2, 1992, p. 178-186.
  • Kenneth A. Cunefare. The effect of modal interaction on sound radiation from vibrating structures. AIAA Journal, Vol. 30, Issue 12, 1992, p. 2819-2828.
  • Borgiotti G. V. The power radiated by a vibrating body in an acoustic fluid and its determination from boundary measurements. Journal of the Acoustical Society of America, Vol. 88, Issue 4, 1990, p. 1884-1893.
  • Photiais D. M. The relationships of singular value decomposition to wave-vector by three-dimensional structures. Journal of the Acoustical Society of America, Vol. 88, Issue 4,1990, p. 1152-1159.
  • Cunfeare K. A. On the exterior acoustic radiation modes of structure. Journal of the Acoustical Society of America, Vol. 96, Issue 4, 1994, p. 2302-2312.
  • Curnfare K. A. The radiation modes of baffled finite plates. Journal of the Acoustical Society of America, Vol. 98, Issue 3, 1995, p. 1570-1580.
  • Elliott S. J. Radiation modes and the active control of sound power. Journal of the Acoustical Society of America, Vol. 94, Issue 4, 1993, p. 2194-2204.
  • D. Guicking, K. Karcher. Active impedance control for one-dimensional sound. Journal Vibration Acoustic Stress Reliability, Vol. 106, Issue 1, 1984, p. 393-396.
  • D. Guicking, E. Lorenz. An Active sound absorber with porous plate. Journal Vibration Acoustic Stress Reliability, Vol. 106, Issue 2, 1984, p. 389-392.
  • Guicking D., Karcher K., Roll wage M. Active control of the acoustic reflection coefficient at low frequencies. In Proceeding of Inter-Noise’83, 1983, p. 419-422.
  • Guicking D., Albrecht M. An electric loudspeaker for active acoustic systems. Journal of vibration, Acoustics, Stress, and Reliability in design, Vol. 106, Issue 1, 1984, p. 397-398.
  • C. R. Fuller, M. J. Bronzel, C. H. Gentry, D. E. Whittington. Control of sound radiation/reflection with adaptive foams. In Proceedings of Noise-Con 1994 Fort Lauderdale, Florida.
  • Thenail D., Lacour O., Galland M. A., Furtoss M. The active control of wall impedance. Acoustic, Vol. 83, Issue 1, 1983, p. 1039-1044.
  • M. Furstoss, D. Thenail, M. A. Galland. Surface impedance control for sound absorption: Direct and hybrid passive/active strategies. Journal of sound and vibration, Vol. 203, Issue 2, 1997, p. 219-236.

About this article

Received
11 September 2013
Accepted
12 November 2013
Published
31 December 2013
Keywords
reflected sound pressure
simply supported plate
piezoelectric ceramic
active sound absorption