Predicting frequencies, decays, levels, and phases of modes of a circular membrane (drum head)?

I have been working on some acoustic synthesis software for modelling several acoustic instruments. I am using the principle of modal synthesis, where every mode of vibration is individually synthesized, generally using a resonant bandpass filter (or sine wave). I have used this approach to get great results on stringed instruments like guitar/cello where the relationship between modes is simple and reasonably predictable. However, I am struggling with drums.

As a 2D circular membrane, the vibrational properties of a drumhead (and even more so, a two headed drum with a drum shell between) are much more complicated than a simple string. However, I still believe this should be possible to synthesize well given modern computing power. With the correct data, I can easily synthesize 500+ modes simultaneously to give a likely good sound.

Imagine a simple 2D circular membrane excited by a strike of energy at it's dead center. The information I ideally need for each mode is:

1) Frequency ratio relative to the fundamental (0,1) - this is easily given in an ideal circular membrane by the Bessel zeros. Things would be different in a "nonideal" membrane, but this is easy to get from any program like ANSYS either way, and I would prefer an "ideal" simulation first.

2) Maximum amplitude of sound of each mode relative to the fundamental (0,1). If a point of reference is needed, let's say we are measuring sound at 2-3" from the center of the excited membrane.

3) Decay rate of each mode, where decay rate is given by time to reach 1/e amplitude relative to that mode's highest initial amplitude, or alternatively in dB/s.

4) Time delay from excitation to beginning of a mode's oscillation in ms or in radians/degrees for that mode's frequency.

If I have a table of that data for the first 500 significant resonances of a circular membrane (or full drum model), I can easily put that into synthesis to see what I get.

ANSYS is okay even in simple usage for providing basic modal frequencies (#1). But I am uncertain if/how it or another program can possibly provide #2-4 on that list.

Is this a very simple or challenging set of data to try to get? How would you approach it, ie. with what program or modelling technique?

Alternatively, the biggest component of the "sound" besides the frequency ratios is the decay rates. If you are aware of any equation that can, using an arbitrary damping coefficient "c", express the theoretical decay rate of any (m,n) mode from a given decay rate of (0,1), that would likely work well enough too. It was easy enough in strings to work in this way, but I don't think it will be so easy in 2D. I'm sure such an equation can be derived, but I don't know it, and I'm not sure anyone else does either. I am hoping to be able to work that relationship out from modeled data if I can get modeling working.

Thanks.

• How do you account for tension - is a real drumskin tensioned perfectly evenly in all directions? Is the thickness constant? Is it likely that those natural variations are what makes the natural sound « warmer » compared to synthesized ? – Solar Mike Jun 17 '18 at 6:02
• I want to synthesize a perfectly ideal drumskin with perfectly ideal tension. I am not interested in "warmth" but rather capturing the ideal relationship between the frequencies, levels, decays, and phases of the modes. The Bessels zeros give the ideal frequency ratios, but not the rest. I am sure there must be equations that can predict the rest, but they are beyond my knowledge. Any ideas, particularly with respect to the decay rates? The math was easy for a string but I don't know where to start working in two dimensions. That is why I proposed modelling first, and working back from that. – mike Jun 17 '18 at 6:35
• You are missing a very important part of the model: the elasticity of the air inside the drum couples the vibration modes of the heads together, and the combined mode frequencies may be completely different from a single membrane in a vacuum, which is what the Bessel function model will give you. Try doing a fluid-solid coupled vibration analysis in Ansys. – alephzero Jun 17 '18 at 8:12
• To model the damping, you also need a coupled fluid-solid model that models the rate at which energy is radiated from the drum head into the air. Acoustic simulation software can do that. You can't easily do it with a standard FE program because the volume of air outside the drum that you need to model is (theoretically) infinite - if you try to model the drum in a "room", i.e. a "large box with rigid walls" the results will be confusing because of all the room resonances! Alternatively, estimate the decay time for the sound of some real drums, and set your damping factor(s) to match that. – alephzero Jun 17 '18 at 8:18
• Hi @alephzero. I am not worried about coupling of drum heads to start, although in a subsequent model I will want this. Single headed drums, even without a proper drum shell definitely sound drum-like! eg. youtube.com/watch?v=z4P7FTDqqww So to begin with I just need to model the properties of an ideal circular membrane and get data from that. From what you say, a fluid-solid coupled vibration analysis in ANSYS will have this capability? Will it be able to give me decay rate in dB/s of the modes? If so, any tips on how to get this from a basic membrane model? Thanks a bunch! – mike Jun 17 '18 at 17:13

I got a reply from someone on another site suggesting the following:

If the air damps it linearly enough, you can probably solve it analytically. Use plate theory to generate a PDE, then work out all the eigenmodes. The decay rate will be determined by the real components of the eigenvalues, and can be converted into dBs-1 using a few logs. https://en.wikipedia.org/wiki/Vibration_of_plates#Isotropic_Kirchhoff%E2%80%93Love_plates

Does this sound like a reasonable approach, and if so, any further tips on how to go about this? Should the air "damp it linearly enough"?

An analytical solution would be more ideal than a modelling solution.

I see partial differential equations for a circular membrane explained here: http://ramanujan.math.trinity.edu/rdaileda/teach/s12/m3357/lectures/lecture_3_29.pdf

But I have no formal physics/math past the 100 level courses I took in undergrad so this is going to take some learning curve for me to sort out. Any help if this makes sense?

Thanks

• a couple of things to consider, the drum skin material properties will almost certainly be non-linear whether the skin is animal (in which case definitely non-linear) or a polymer. If the vibrations are acting in the linear region of deformation then no problem but I wouldn't be surprised if they aren't. Also there will be damping from the drum skin material and boundary. I think you mentioned you weren't that bothered about modelling the boundary? It might be worth investigating what visco-elastic properties the skin material has tho. This could be significant & thickness dependent. – DrBwts Jun 19 '18 at 16:10