/ challenges-old / cymbal
Sign in
challenges-old / cymbal
cymbal / challenges / cymbal-tapered-v1 / README.md
4.6 KBPreviewCodeBlameHistoryRaw

cymbal-tapered-v1: the hammered cymbal lathed to a tapered thickness#

This is the fifth challenge of the cymbal tier. It is cymbal-hammered-v1 with one change: the thickness is no longer uniform. It is 1.2 mm over the bell and decreases linearly along the bow to 0.7 mm at the rim, the taper that lathing gives a 16 inch crash. Everything else, the bow, the bell, the hole, the 150 dents, the material, the damping, the strike, the listener, the output and the error measures, is unchanged, so the difference between the two challenges is the effect of the taper alone. A thinner rim lowers the modes that live near the edge and makes them radiate and cascade differently from a uniform sheet; a heavier center raises the low modes of the whole body. The deliverable is again six seconds of pressure at the listener.

The machine-readable version of this specification is challenge.json.

Model#

Reference geometry#

The mid-surface SS is that of cymbal-hammered-v1: the graph Z(x,y)Z(x, y) of the bow, the bell and the dents over the annulus rh≤r≤Rr_h \le r \le R, with unit normal n(x,y)\mathbf{n}(x, y). The thickness is a function of the distance from the axis,

h(r)={hc,r≤rb,hc+(he−hc) r−rbR−rb,rb<r≤R,h(r) = \begin{cases} h_c, & r \le r_b, \\[1ex] h_c + (h_e - h_c)\, \dfrac{r - r_b}{R - r_b}, & r_b < r \le R, \end{cases}

constant over the bell and linear on the bow, continuous at the bell edge. The body is the set of points within half the local thickness of the mid-surface along its normal,

Ω={X+s n(X):X∈S, ∣s∣<h(r(X))/2},\Omega = \{ \mathbf{X} + s\, \mathbf{n}(\mathbf{X}) : \mathbf{X} \in S,\ |s| < h(r(\mathbf{X})) / 2 \},

so the top face is S+={X+(h/2) n}S^+ = \{ \mathbf{X} + (h/2)\, \mathbf{n} \} and the bottom face S−={X−(h/2) n}S^- = \{ \mathbf{X} - (h/2)\, \mathbf{n} \} with the local hh, and the rim and the hole edge are the ruled surfaces swept by the normal segment of length h(R)h(R) at r=Rr = R and h(rh)=hch(r_h) = h_c at r=rhr = r_h. Since hh varies, the faces are not parallel to the mid-surface, by a slope of ∣h′∣/2=(hc−he)/(2(R−rb))≈0.0017|h'| / 2 = (h_c - h_e) / (2 (R - r_b)) \approx 0.0017, negligible against the slopes of the bow.

Material, equations of motion, strike, and output#

All as before: the Saint Venant-Kirchhoff solid with Kelvin-Voigt viscosity and mass-proportional loss, traction-free everywhere except the strike, at rest at t=0t = 0, free, with the rigid motion linearized in the output. The strike is the same vertical raised-cosine pulse at 0.85R0.85 R on the xx axis; its normalization constant CC is the surface integral of the footprint bump over the top face, whose position now depends on the local thickness (h≈0.80h \approx 0.80 mm at the strike), and is given in the parameter table. The output is the same integral of the normal acceleration over S+S^+ with the retarded time and the same wav format and scale.

Parameters#

Those of cymbal-hammered-v1 with the uniform thickness hh replaced by the two thickness values.

Symbol Value Unit Meaning
hch_c 0.0012 m thickness over the bell and at its edge
heh_e 0.0007 m thickness at the rim
CC see challenge.json m^2 strike normalization integral

The mass of the body is close to that of the uniform 1 mm disk (the mean thickness over the area is about 0.86 mm), so the drift and the kick of the strike are similar; the bending stiffness, which goes with h3h^3, is 1.7 times the uniform value at the bell and 0.34 times it at the rim. The dents, a millimeter deep, are now deeper than the local thickness near the rim (0.7 mm).

Error measures, reference, field export#

As for cymbal-hammered-v1, including the caution that at this strike the response is chaotic beyond about 0.15 s.

Appendix: shallow-shell reduction#

The Marguerre equations apply with variable coefficients: the mass per unit area ρh\rho h, the bending stiffness D=Eh3/(12(1−ν2))D = E h^3 / (12 (1 - \nu^2)) and the membrane stiffness EhE h all vary with rr. The bending operator becomes ∂xx[D(wxx+νwyy)]+∂yy[D(wyy+νwxx)]+2(1−ν)∂xy[Dwxy]\partial_{xx}[D (w_{xx} + \nu w_{yy})] + \partial_{yy}[D (w_{yy} + \nu w_{xx})] + 2 (1 - \nu) \partial_{xy}[D w_{xy}], and the compatibility equation for the Airy function becomes

∂yyΦyy−νΦxxEh+∂xxΦxx−νΦyyEh−2 ∂xy(1+ν)ΦxyEh=−12L(w,w)−L(w0,w),\partial_{yy} \frac{\Phi_{yy} - \nu \Phi_{xx}}{E h} + \partial_{xx} \frac{\Phi_{xx} - \nu \Phi_{yy}}{E h} - 2\, \partial_{xy} \frac{(1 + \nu) \Phi_{xy}}{E h} = -\tfrac{1}{2} \mathcal{L}(w, w) - \mathcal{L}(w_0, w),

which for constant hh reduces to ∇4Φ/(Eh)\nabla^4 \Phi / (E h) on the left. A modal solver that uses the eigenfunctions of this variable-coefficient operator (with Φ=∂Φ/∂n=0\Phi = \partial \Phi / \partial n = 0 on the rim) in place of the clamped biharmonic modes recovers the same structure of equations as for the uniform shell.

↑↓ move↵openescclose