## Physics-based methods for modeling nonlinear vibrating strings,” Acta Acustica united with (2005)

Venue: | Acustica |

Citations: | 3 - 0 self |

### BibTeX

@ARTICLE{Pakarinen05physics-basedmethods,

author = {Jyri Pakarinen and Vesa Välimäki and Matti Karjalainen},

title = {Physics-based methods for modeling nonlinear vibrating strings,” Acta Acustica united with},

journal = {Acustica},

year = {2005},

pages = {312--325}

}

### OpenURL

### Abstract

Nonlinearity in the vibration of a string is responsible for interesting acoustical features in many plucked string instruments, resulting inacharacteristic and easily recognizable tone. For this reason, synthesis models have to be capable of modeling this nonlinear behavior, when high quality results are desired. This study presents two novel physical modeling algorithms for simulating the tension modulation nonlinearity in vibrating strings in a spatially distributed manner. The first method uses fractional delay filters within a digital waveguide structure, allowing the length of the string to be modulated during run time. The second method uses a nonlinear finite difference approach, where the string state is approximated between sampling instants using similar fractional delay elements, thus allowing run-time modulation of the temporal sampling location. The magnitude of the tension modulation is evaluated from the elongation of the string at every time step in both cases. Simulation results of the two models are presented and compared. Real-time sound synthesis of the kantele, a traditional Finnish plucked-string instrument with strong effect of tension modulation, has been implemented using the nonlinear digital waveguide algorithm. PACS no. 43.75.Gh, 43.75.Wx 1.

### Citations

171 | Physical modeling using digital waveguides - Smith - 1992 |

165 |
The physics of musical instruments
- Fletcher, Rossing
- 1991
(Show Context)
Citation Context ... � (2) K K x K z where K is the string tension and is the linear mass density of the string. Obviously, every real string vibration decays with time. This results mainly from three damping mechanisms =-=[5]-=-: (1) air damping, (2) internal damping in the string medium, and (3) mechanical energy transfer through string terminations. If the losses are simply divided into frequency-dependent and frequency in... |

61 |
Principles of digital waveguide models of musical instruments
- Smith
- 1998
(Show Context)
Citation Context ...single delay-loop with two consecutive samples averaged, as is done in the classic KarplusStrong algorithm [21]. An excellent introduction to DWGs used in modeling musical instruments can be found in =-=[22]-=-. The entire digital waveguide methodology is based on the traveling-wave solution of the wave equation. This means that the solution to equation (1) can be seen as a superposition of two waveforms tr... |

58 |
Splitting the unit delay: Tools for fractional delay filter design
- Laakso, Valimaki, et al.
- 1996
(Show Context)
Citation Context ...l string is shown in Figure 3. Correct tuning of the waveguide can be enabled by adding a fractional delay filter (i.e. a filter capable of producing also noninteger delay values) inside the DWG loop =-=[24]-=-. Frequency-independent losses can easily be modeled in a DWG structure by inserting simple scaling coefficients in the ideal DWG string structure [22]. Frequency� -1 314Pakarinen et al.: Modeling no... |

56 |
Digital synthesis of plucked-string and drum timbres
- Karplus, Strong
- 1983
(Show Context)
Citation Context ...ies, and returns to its initial position. At its simplest form, this can be implemented as a single delay-loop with two consecutive samples averaged, as is done in the classic KarplusStrong algorithm =-=[21]-=-. An excellent introduction to DWGs used in modeling musical instruments can be found in [22]. The entire digital waveguide methodology is based on the traveling-wave solution of the wave equation. Th... |

48 |
Pluckedstring models: from the Karplus–Strong algorithm to digital waveguides and beyond
- Karjalainen, Välimäki, et al.
- 1998
(Show Context)
Citation Context ...e coupling to being one-directional in order to avoid stability problems which would otherwise rise due to the feedback loop formed from the interconnected strings, as suggested by Karjalainen et al. =-=[32]-=-. Clearly, the actual physical coupling is two-directional. The elongation approximations of the strings and the resulting allpass filter coefficient values are evaluated separately for the two DWG mo... |

47 |
Discrete-Time Modeling of Acoustic Tubes Using Fractional Delay Filters
- Välimäki
- 1995
(Show Context)
Citation Context ...hase response error caused by the allpass filter is not considered to pose a problem since its effect is negligible in the audio frequency range, assuming the sampling frequency to be reasonably high =-=[26]-=-. 3.2. Distributed nonlinear DWG string Previous works [27], [28], [29] use a single fractional delay element in a single-polarization string model, or a DWG string terminated with a nonlinear double-... |

34 |
Synthesizing musical sounds by solving the wave equation for vibrating objects. Part I
- Hiller, Ruiz
- 1971
(Show Context)
Citation Context ...nd (3) mechanical energy transfer through string terminations. If the losses are simply divided into frequency-dependent and frequency independent terms, the lossy 1-D wave equation can be written as =-=[6]-=-, [7]: Kyxx(t� x) = ytt(t� x) +d1yt(t� x) ; d2yxxx(t� x): (3) Here, d1 and d2 are coefficients that simulate the frequency-independent and frequency-dependent damping, respectively. Real strings also ... |

33 |
Numerical simulations of piano strings. I. A physical model for a struck string using finite difference methods
- Chaigne, Askenfelt
- 1994
(Show Context)
Citation Context ...) mechanical energy transfer through string terminations. If the losses are simply divided into frequency-dependent and frequency independent terms, the lossy 1-D wave equation can be written as [6], =-=[7]-=-: Kyxx(t� x) = ytt(t� x) +d1yt(t� x) ; d2yxxx(t� x): (3) Here, d1 and d2 are coefficients that simulate the frequency-independent and frequency-dependent damping, respectively. Real strings also exper... |

30 |
The physics of the violin
- Cremer
- 1985
(Show Context)
Citation Context ... [11], [12], and [13], to name few. The bow-string nonlinearity is mainly caused by the stick-slip contact between the string and the bow. Also this interaction is coveredinmanystudies,[14],[15], and =-=[16]-=- present some of them. In the following, we will study more thoroughly the tension modulation nonlinearity. When a real string is displaced, its length, and therefore also its tension is increased. Wh... |

19 | Modeling of tension modulation nonlinearity in plucked strings
- Tolonen, Välimäki, et al.
- 2000
(Show Context)
Citation Context ...d to pose a problem since its effect is negligible in the audio frequency range, assuming the sampling frequency to be reasonably high [26]. 3.2. Distributed nonlinear DWG string Previous works [27], =-=[28]-=-, [29] use a single fractional delay element in a single-polarization string model, or a DWG string terminated with a nonlinear double-spring [30] to simulate the nonlinear string. This is done in ord... |

19 |
A passive nonlinear digital filter design which facilitates physics-based sound synthesis of highly nonlinear musical instruments
- PIERCE, DUYNE
- 1997
(Show Context)
Citation Context ...Distributed nonlinear DWG string Previous works [27], [28], [29] use a single fractional delay element in a single-polarization string model, or a DWG string terminated with a nonlinear double-spring =-=[30]-=- to simulate the nonlinear string. This is done in order to reduce the computational complexity of the model, butithas also some shortcomings. Since the system is nonlinear, the FD elements cannot be ... |

19 |
On the use of finite differences for musical synthesis. Application to plucked stringed instruments
- Chaigne
- 1992
(Show Context)
Citation Context ...ferences had already earlier been used in mathematics for numerical solving of partial differential equations. Afine introduction to FDM in the synthesis of plucked string instruments can be found in =-=[35]-=-. Below, we follow the guidelines provided in [22] in deriving the FDM recurrence equation. 5.1. Ideal finite difference string The partial derivatives in the 1-D wave equation (1) can be replaced by ... |

19 |
Numerical simulations of xylophones. i. time-domain modelling of the vibrating bars
- Chaigne
- 1997
(Show Context)
Citation Context ...f the two algorithms were also compared. As stated in section 5.1, the explicit finite difference scheme was chosen for simplicity. Another option would be to use an implicit scheme, such as a scheme =-=[40]-=-, where the temporal and spatial derivatives of the wave equation (equation 1) are averaged in space and time, re323ACTA ACUSTICA UNITED WITH ACUSTICA Vol. 91 (2005) Pakarinen et al.: Modeling nonlin... |

14 |
Nonlinear generation of missing modes on a vibrating string
- Legge, Fletcher
- 1984
(Show Context)
Citation Context ...the bridge is able to move in the z-direction, this vibration will be coupled with the string. With a rigid bridge, no coupling will take place, and the missing harmonics will not be generated (after =-=[19]-=-). to the instrument body, and therefore cannot be heard directly. In instruments where this is not the case, partials created by TM can be found at twice the frequencies of the transversal modes, and... |

14 | Pluckedstring synthesis algorithms with tension modulation nonlinearity
- Välimäki, Tolonen, et al.
- 1992
(Show Context)
Citation Context ...sidered to pose a problem since its effect is negligible in the audio frequency range, assuming the sampling frequency to be reasonably high [26]. 3.2. Distributed nonlinear DWG string Previous works =-=[27]-=-, [28], [29] use a single fractional delay element in a single-polarization string model, or a DWG string terminated with a nonlinear double-spring [30] to simulate the nonlinear string. This is done ... |

12 |
BlockCompiler: efficient simulation of acoustic and audio systems
- Karjalainen
- 2003
(Show Context)
Citation Context ...ithm used is efficient enough to process a five-string kantele model on an ordinary laptop computer at a 44.1 kHz samplingrate. A detailed description of the BlockCompiler is presented by Karjalainen =-=[33]-=-. 5. Finite difference approach In the previous section we discussed string synthesis via discretizing the d’Alembert’s solution to the 1-D wave equation. Another approach is to discretize the wave eq... |

11 | Välimäki: Acoustical analysis and model-based sound synthesis of the kantele - Erkut, Karjalainen, et al. |

10 |
1-D digital waveguide modeling for improved sound synthesis
- Karjalainen
- 2002
(Show Context)
Citation Context ...=0.Thisresults in a phase-inverting termination, which suits perfectly the case of an ideal string. For other types of string termination, several models have been introduced (see e.g. [6], [35], and =-=[36]-=-). Generally, the nonrigid string terminations lead to frequency-dependent losses in the string model. For the FDM string excitation, a useful method has been proposed in [36]. It is conceptually simp... |

9 |
Digital waveguides vs. finite difference schemes: Equivalence and mixed modeling
- Karjalainen, Erkut
- 2004
(Show Context)
Citation Context ...f the synthesis engine. There has also been a trend towards comparing and unifying the existing physical modeling methods with a focus on generating more flexible and efficient sound synthesis models =-=[1]-=-, [2]. This paper discusses and compares two physical models for simulating the behavior of the nonlinear string: a spatially distributed nonlinear digital waveguide string recently developed by the a... |

8 |
Dispersion of Waves in Piano Strings
- Podlesak, Lee
- 1988
(Show Context)
Citation Context ...n the sound synthesis of strings with relatively high elasticity and small diameters, the effects of stiffness are not discussed here. More in-depth study about this topic can be found for example in =-=[8]-=-. 2.2. Nonlinear string Although the linearity assumption usually simplifies calculations, the vibration of real strings can be considered linear in the most coarse approximations only. More realistic... |

7 |
Finite difference method vs. digital waveguide method in string instrument modeling and synthesis
- Erkut, Karjalainen
(Show Context)
Citation Context ... synthesis engine. There has also been a trend towards comparing and unifying the existing physical modeling methods with a focus on generating more flexible and efficient sound synthesis models [1], =-=[2]-=-. This paper discusses and compares two physical models for simulating the behavior of the nonlinear string: a spatially distributed nonlinear digital waveguide string recently developed by the author... |

7 | Principle of fractional delay filters
- Valimaki, Laakso
- 2000
(Show Context)
Citation Context ...er approximating a given fractional delay [24], and partly because the first-order allpass filter is the best choice for the fractional delay element when delay values around unity are to be obtained =-=[25]-=-. The phase response error caused by the allpass filter is not considered to pose a problem since its effect is negligible in the audio frequency range, assuming the sampling frequency to be reasonabl... |

6 |
Piano string excitation in the case of a small hammer mass, J Acoust Soc Am 79
- Hall
- 1986
(Show Context)
Citation Context ... string instruments. These phenomena fall out of the scope of this paper. The hammer-string nonlinearity is mainly due to the compression of the hammer felt, and it is covered in many earlier studies =-=[9]-=-, [10], [11], [12], and [13], to name few. The bow-string nonlinearity is mainly caused by the stick-slip contact between the string and the bow. Also this interaction is coveredinmanystudies,[14],[15... |

6 |
Model analysis of a hammer-string interaction
- Suzuki
- 1987
(Show Context)
Citation Context ...ts. These phenomena fall out of the scope of this paper. The hammer-string nonlinearity is mainly due to the compression of the hammer felt, and it is covered in many earlier studies [9], [10], [11], =-=[12]-=-, and [13], to name few. The bow-string nonlinearity is mainly caused by the stick-slip contact between the string and the bow. Also this interaction is coveredinmanystudies,[14],[15], and [16] presen... |

6 |
Modeling the longitudinal vibration of piano strings
- Bank, Sujbert
- 2003
(Show Context)
Citation Context ...artials created by TM can be found at twice the frequencies of the transversal modes, and they are called phantom partials. Thegeneration of these partials is discussed in more detail in [17], and in =-=[18]-=-. Due to decaying vibration, the string displacement is clearly greater just after the pluck than at the end of the vibration. Therefore also the amplitude envelope of TM has the same form as the tran... |

5 |
On the fundamentals of bowed string dynamics Acustica 43 93–108
- McIntyre, Woodhouse
- 1979
(Show Context)
Citation Context ...dies [9], [10], [11], [12], and [13], to name few. The bow-string nonlinearity is mainly caused by the stick-slip contact between the string and the bow. Also this interaction is coveredinmanystudies,=-=[14]-=-,[15], and [16] present some of them. In the following, we will study more thoroughly the tension modulation nonlinearity. When a real string is displaced, its length, and therefore also its tension i... |

5 | Karjalainen: Simulation of plucked strings exhibiting tension modulation drivingforce
- Tolonen, Erkut, et al.
(Show Context)
Citation Context ...he two one-polarizational vibrations are denoted by fy(n) and fz(n). Theconnection from the elongation approximation block to the output simulates the direct coupling of the TM to the instrument body =-=[20]-=-. A scaling coefficient, denoted by C, isused to control the amount of this coupling. The output of the whole two-polarization string model is finally presented as a force signal fout(n) excerted to t... |

4 |
Piano string excitation II: General solution for a hard narrow hammer
- Hall
- 1987
(Show Context)
Citation Context ...ng instruments. These phenomena fall out of the scope of this paper. The hammer-string nonlinearity is mainly due to the compression of the hammer felt, and it is covered in many earlier studies [9], =-=[10]-=-, [11], [12], and [13], to name few. The bow-string nonlinearity is mainly caused by the stick-slip contact between the string and the bow. Also this interaction is coveredinmanystudies,[14],[15], and... |

4 |
Piano string excitation III: General solution for a soft narrow hammer
- Hall
- 1987
(Show Context)
Citation Context ...truments. These phenomena fall out of the scope of this paper. The hammer-string nonlinearity is mainly due to the compression of the hammer felt, and it is covered in many earlier studies [9], [10], =-=[11]-=-, [12], and [13], to name few. The bow-string nonlinearity is mainly caused by the stick-slip contact between the string and the bow. Also this interaction is coveredinmanystudies,[14],[15], and [16] ... |

4 |
Virtual strings based on a 1-D FDTD waveguide model: Stability, losses, and traveling waves
- Erkut, Karjalainen
- 2002
(Show Context)
Citation Context ...x m,butsince homogeneous strings are considered here, this dependencyisomitted. Values q = p 2 � jpj 1 (24) ensure the stability of a linear finite difference string with frequency-independent losses =-=[37]-=-. Note that the sign difference of p and q in [37] has already been taken care of in equation (23). Modeling of frequency-dependent losses by discretizing the lossy wave equation leads to an implicit ... |

4 | Audibility of initial pitch glides in string instrument sounds
- Järveläinen, Välimäki
- 2001
(Show Context)
Citation Context ...ntal dotted line approximates the mean value of perceptual detection threshold of an initial pitch glide. The psychoacoustic detection threshold in the frequency region of these tones is about 5.4 Hz =-=[38]-=-. This shows that the fundamental frequency glide is an audible phenomenon in plucked string instruments such as the kantele even at modest plucking amplitudes, and thus it must be included in a synth... |

2 |
Välimäki: Modeling and real-time synthesis of the kantele using distributed tension modulation
- Pakarinen, Karjalainen, et al.
(Show Context)
Citation Context ...This paper discusses and compares two physical models for simulating the behavior of the nonlinear string: a spatially distributed nonlinear digital waveguide string recently developed by the authors =-=[3]-=-, and a new spatially distributed nonlinear finite difference string model, created as a part of a Master’s Thesis work at Helsinki University of Technology [4]. Since the purpose of physics-based mod... |

2 |
Spatially distributed computational modeling of a nonlinear vibrating string. Diploma Thesis
- Pakarinen
- 2004
(Show Context)
Citation Context ...tring recently developed by the authors [3], and a new spatially distributed nonlinear finite difference string model, created as a part of a Master’s Thesis work at Helsinki University of Technology =-=[4]-=-. Since the purpose of physics-based modeling is to simulate the physical phenomena in the system of interest, it is of paramount importance to be familiar with the behavior of the real-world case bef... |

2 |
Idealised models of a bowed string. Acustica 79
- Woodhouse
- 1993
(Show Context)
Citation Context ...[9], [10], [11], [12], and [13], to name few. The bow-string nonlinearity is mainly caused by the stick-slip contact between the string and the bow. Also this interaction is coveredinmanystudies,[14],=-=[15]-=-, and [16] present some of them. In the following, we will study more thoroughly the tension modulation nonlinearity. When a real string is displaced, its length, and therefore also its tension is inc... |

2 |
Välimäki: New designs for the kantele with improved sound radiation
- Pölkki, Erkut, et al.
(Show Context)
Citation Context ...oducing beating [29]. A more detailed structure and acoustical analysis of the kantele can be found in [29]. A study of the history of kantele and an acoustically improved new design are presented in =-=[31]-=-. 4.2. A novel kantele string model The novel synthesis model of a single kantele string is constructed using two single-polarization time-varying DWG models, illustrated in Figure 5, and connecting t... |

2 |
Strikwerda: Finite difference schemes and partial differential equations
- C
- 1989
(Show Context)
Citation Context ...e nonlinear DWG waveguide has noreal stability problems when synthesis of natural plucked-instrument sounds are desired. We studied the stability of the NFDTD algorithm using the Von Neumann analysis =-=[39]-=- in the time-invariant case, i.e. parameter a of equation (26) was kept constant. The basic idea of this method is to calculate the spatial Fourier spectrum of the system under discussion at two conse... |

1 |
Boutillon: Model for piano hammers: Experimental determination and digital simulation
- unknown authors
- 1988
(Show Context)
Citation Context ...phenomena fall out of the scope of this paper. The hammer-string nonlinearity is mainly due to the compression of the hammer felt, and it is covered in many earlier studies [9], [10], [11], [12], and =-=[13]-=-, to name few. The bow-string nonlinearity is mainly caused by the stick-slip contact between the string and the bow. Also this interaction is coveredinmanystudies,[14],[15], and [16] present some of ... |

1 |
Conklin: Generation of partials due to nonlinear mixing in a stringed instrument
- A
- 1999
(Show Context)
Citation Context ...t the case, partials created by TM can be found at twice the frequencies of the transversal modes, and they are called phantom partials. Thegeneration of these partials is discussed in more detail in =-=[17]-=-, and in [18]. Due to decaying vibration, the string displacement is clearly greater just after the pluck than at the end of the vibration. Therefore also the amplitude envelope of TM has the same for... |