# Loudspeaker Analysis and Experiments: Part II

If you haven’t read Part I of this series, it is strongly recommended that the user starts there.

__Brüel and Kjær__published a quality report a while ago on the transfer function for a loudspeaker in free air:

*ω*is the frequency input to the system,

*C*is the mechanical compliance,

_{mr}*M*is the mass of the acoustic and mechanical diaphragm and coil loading,

_{mr}*R*is the mechanical resistance, and

_{mr}*Bl*is a product of magnetic field strength and coil length which determines the amount of force imposed on the system in accordance with

__Faraday's law of induction.__

The constants above can be found for specific changes in impedance due to frequency, which we will continue to explore in the sections to follow.

If we take the derivative of the transfer function above:

The frequency at which the derivative is equal to zero is the value which we call the resonance frequency:

*ω*:

_{r}and rewriting in a more familiar format:

*C*is the mechanical compliance, and

_{mr}*M*is the mass of the acoustic and mechanical diaphragm and coil loading.

_{mr}One method for finding the mechanical compliance and mass is to add known masses to the mechanical mass and observe the change in resonance frequency. This allows us to setup a relationship between resonance frequency and added mass in order to solve for mechanical compliance and mass in a linear function. We can start by writing the added mass and new resonance as a function of one another:

*m*is the added mass amount and

_{i}*f*is the resonance at that added mass. We can rewrite the above equation in linear form:

_{r,i}where:

*C*and

_{mr}*M*is we use multiple masses and resonance calculations to fit a line and solve for slope and intercept. I used resonances for added masses at: 0g (actual resonance), 2g, 5g, 10g, 15g, 20g. The linear fit based on the square relationship above is shown below:

_{mr}#### Linear fit for finding the mechanical compliance and mass based on the added mass technique

The plot above gives us the approximate values for the compliance and mass:

The Python code for implementing the added mass and resonance fit is shown below. All the user needs to do is add the mass and measure the resonance frequency, then input the values into the code and it prints out the plot shown above.

## linear fitting the mass loading of loudspeaker import numpy as np from scipy import stats import matplotlib.pyplot as plt plt.rcParams["font.family"] = 'Times New Roman' plt.rcParams["mathtext.fontset"] = 'stix' plt.style.use('ggplot') m = np.array([0.0,2.0,5.0,10.0,15.0,20.0])*0.001 f = np.array([86.4,74.0,62.6,52.1,45.0,40.0]) y = 1.0/((2.0*np.pi*f)**2.0) slope,intercept,r_val,p_val,std_err = stats.linregress(m,y) fig = plt.figure(figsize=(12,9)) ax = fig.add_subplot(111) plt.plot(m,y,label='data',linewidth=3.0) plt.plot(m,(m*slope)+intercept,label='fit',linewidth=3.0) plt.legend(fontsize=20) diff_val = (y-((m*slope)+intercept)) MAPE = 100.0*np.mean(np.abs(diff_val/y)) print('MAPE: {0:2.2f}%'.format(MAPE)) stats_str = '% Error = {0:2.2f}%\n$R^2$ = {1:2.3f}'.format(MAPE,r_val**2.0) plt.text(0.0, .000012, stats_str, size=20, ha="left", va="center", bbox=dict(boxstyle="round", ec=(0.2, 0.2, 0.2), fc=(0.95, 0.95, 0.95), ) ) eqn_str = '$(2\pi f_{r,i})^{-2} = $'+'{0:.2E}'.format(slope)+'$m_i$+'+'{1:.2E}'.format(slope,intercept) eqn_str = eqn_str.replace('E','e') plt.text(np.mean(m), .000005, eqn_str, size=20, ha="left", va="center", ) plt.xlabel('$m_i$ [grams]',fontsize=20) plt.ylabel('$(2 \pi f_{r,i})^{-2}$ [1/Hz$^2$]',fontsize=20) plt.savefig('mass_loading.png',dpi=300,facecolor=np.array([252,252,252])/255.0) plt.show()

*V*, which is a measure of the amount of air that the suspension system is capable of moving. It is measured in liters or m

_{ar}^{3}, and generally is a measure of what size enclosure the loudspeaker should be housed in. It can be calculated as follows:

*ρ*is the density of air,

*c*is the speed of sound in air,

*S*is the effective area of the moving speaker cone, and

_{d}*C*is the compliance calculated above. For the speaker cone area I measure across the width of the speaker until about halfway past the speaker surround. This gives an effective diameter of about 90-95% of the cited diameter.

_{mr}For the speaker used in this series, the equivalent compliance volume is calculated as follows:

For comparison with other drivers of similar size and specifications, this value seems reasonable. Now we will continue with the analysis of impedance by investigating the significance of impedance magnitudes in relation to damping and the rest of the parameters associated with loudspeaker characterization.

*C*and

_{mr}*M*is to quanitfy the quality factor (inverse of damping),

_{mr}*Q*, for both mechanical and electrical domains. The quality factor for an oscillating system

__can be defined as__:

*Q*is the quality factor, which in this case is defined by two resonance sidelobe frequencies. Typically, the sidelobe frequencies are chosen at -3dB from the peak value (

*Z*) at resonance, which is also one half of the value. Therefore, if we wanted to define the

_{max}*Q*, we would have to solve for the frequencies where the impedance is half the maximum. In my case, the sidelobe frequencies are roughly 75 Hz and 99 Hz, resulting in:

_{mr}*Q*, can be approximated using the mechanical with the following relation:

_{er}It is customary to combine the quality factors in series to produce a total quality factor:

The *Bl *product, which is a measure of the magnetic field strength with respect to the voice coil, comes quite easily from the quality factor in the electrical domain:

which gives us the result:

*Bl*product above is in Tesla-meters - a unit of magnetic field strenght and length.

which appears to be quite poor, however, in comparison with other loudspeakers - it’s not far from the average for its size and design.

This concludes the investigation into the Thiele-Small parameters for a loudspeaker. A lot of the analysis and calculation was based on decades-old methods and parameters. In the next entry into this series I will focus on experiments with loudspeakers, such as an application with microphones and even an experiment with a non-newtonian fluid.

*See More in Acoustics and Engineering:*