Next  |  Prev  |  Up  |  Top  |  Index  |  JOS Index  |  JOS Pubs  |  JOS Home  |  Search

Mass and Dashpot in Series

This is our first example illustrating a series connection of wave digital elements. Figure F.28 gives the physical scenario of a simple mass-dashpot system, and Fig.F.29 shows the equivalent circuit. Replacing element voltages and currents in the equivalent circuit by wave variables in an infinitesimal waveguides produces Fig.F.30.

Figure F.28: External force driving a mass which in turn drives a dashpot terminated on the other end by a rigid wall.
\includegraphics{eps/massdash}

Figure F.29: Electrical equivalent circuit of the mass and dashpot system of Fig.F.28.
\includegraphics{eps/massdashec}

Figure F.30: Intermediate wave-variable model of the mass and dashpot of Fig.F.29.
\includegraphics{eps/massdashdt}

Figure F.31: Wave digital filter for an ideal force source in parallel with the series combination of a mass $ m$ and dashpot $ \mu $ . The parallel and series adaptors are joined at an impedance $ R$ which is calculated to suppress reflection from port 1 of the series adaptor.
\includegraphics{eps/massdashjunc}

The system can be described as an ideal force source $ f(t)$ connected in parallel with the series connection of mass $ m$ and dashpot $ \mu $ . Figure F.31 illustrates the resulting wave digital filter. Note that the ports are now numbered for reference. Two more symbols are introduced in this figure: (1) the horizontal line with a dot in the middle indicating a series adaptor, and (2) the indication of a reflection-free port on input 1 of the series adaptor (signal $ f^{{+}}_1(n)$ ). Recall that a reflection-free port is always necessary when connecting two adaptors together, to avoid creating a delay-free loop.

Let's first calculate the impedance $ R$ necessary to make input 1 of the series adaptor reflection free. From Eq.(F.44), we require

$\displaystyle R = m + \mu
$

That is, the impedance of the reflection-free port must equal the series combination of all other port impedances meeting at the junction.

The parallel adaptor, viewed alone, is equivalent to a force source driving impedance $ R=m+\mu$ . It is therefore realizable as in Fig.F.22 with the wave digital spring replaced by the mass-dashpot assembly in Fig.F.31. However, we can also carry out a quick analysis to verify this: The alpha parameters are

\begin{eqnarray*}
\alpha_1 &\isdef & \frac{2\Gamma _1}{\Gamma _1+\Gamma _2}
= \frac{2\cdot\infty}{\infty+\frac{1}{m}+\frac{1}{\mu}}
= 2\\
\alpha_2 &\isdef & \frac{2\Gamma _2}{\Gamma _1+\Gamma _2}
= \frac{2\left(\frac{1}{m}+\frac{1}{\mu}\right)}{\infty+\left(\frac{1}{m}+\frac{1}{\mu}\right)}
= 0
\end{eqnarray*}

Therefore, the reflection coefficient seen at port 1 of the parallel adaptor is $ \rho = \alpha_1 - 1 = 1$ , and the Kelly-Lochbaum scattering junction depicted in Fig.F.22 is verified.

Let's now calculate the internals of the series adaptor in Fig.F.31. From Eq.(F.33), the beta parameters are

\begin{eqnarray*}
\beta_1 &\isdef & \frac{2R_1}{R_1+R_2+R_3}
= \frac{2(m+\mu)}{(m+\mu)+m+\mu}
= 1\\
\beta_2 &\isdef & \frac{2R_2}{R_1+R_2+R_3}
= \frac{2\mu}{(m+\mu)+m+\mu}
= \frac{\mu}{m+\mu}\\
\beta_3 &\isdef & \frac{2R_3}{R_1+R_2+R_3}
= \frac{2m}{(m+\mu)+m+\mu}
= \frac{m}{m+\mu}
\end{eqnarray*}

Following Eq.(F.37), the series adaptor computes

\begin{eqnarray*}
f^{{+}}_J(n) &=& f^{{+}}_1(n)+f^{{+}}_2(n)+f^{{+}}_3(n)
= f(n) + f^{{+}}_3(n)
\\ [5pt]
f^{{-}}_1(n) &=& f^{{+}}_1(n) - \beta_1f^{{+}}_J(n)
= f^{{+}}_3(n)\\ [5pt]
f^{{-}}_3(n) &=& f^{{+}}_3(n) - \beta_3f^{{+}}_J(n)\\
&=& (1-\beta_3)f^{{+}}_3(n) - \beta_3 f(n)\\
&=& \frac{\mu}{m+\mu}f^{{+}}_3(n) - \frac{m}{m+\mu} f(n)
\end{eqnarray*}

We do not need to explicitly compute $ f^{{-}}_2(n)$ because it goes into a purely resistive impedance $ \mu $ and produces no return wave. For the same reason, $ f^{{+}}_2(n)\equiv\message{CHANGE eqv TO equiv IN SOURCE}0$ . Figure F.32 shows a wave flow diagram of the computations derived, together with the result of elementary simplifications.

Figure F.32: Wave flow diagram for the WDF modeling an ideal force source in parallel with the series combination of a mass $ m$ and dashpot $ \mu $ .
\includegraphics{eps/massdashwdf}

Because the difference of the two coefficients in Fig.F.32 is 1, we can easily derive the one-multiply form in Fig.F.33.

Figure F.33: One-multiply form of the WDF in Fig.F.32.
\includegraphics{eps/massdashwdfom}



Subsections
Next  |  Prev  |  Up  |  Top  |  Index  |  JOS Index  |  JOS Pubs  |  JOS Home  |  Search

[How to cite this work]  [Order a printed hardcopy]  [Comment on this page via email]

``Physical Audio Signal Processing'', by Julius O. Smith III, W3K Publishing, 2010, ISBN 978-0-9745607-2-4
Copyright © 2023-08-20 by Julius O. Smith III
Center for Computer Research in Music and Acoustics (CCRMA),   Stanford University
CCRMA