Beam README: the near 2:1 at h = 0.2 is not a resonance (mode 2 antisymmetric)
1 changed file+1−1
src/beam/README.mdmodified+1−1View file
@@ -15,7 +15,7 @@ Live site: https://chaos-apps--listening-for-chaos.vault1-sites.magland.org/beam
1515 - Potential: springs between neighbours (and to each support), k d^2 / 2 with d = u_j - u_{j-1}; bending κ c^2 / 2 at each mass with both neighbours, c = u_{j-1} - 2 u_j + u_{j+1} (supports at u = 0). Pinned support: no extra term (ghost u_{-2} = -u_0 gives c = 0 there). Clamped: a joint at the support with ghost u_{-2} = u_0, c = 2 u_0, energy 2 κ u_0^2. Free end: no spring or joint beyond the last mass.
1616 - Stretching: "local", k3 d^4 / 4 per spring (FPUT); or "Kirchhoff", k3 S^2 / (4 nS) with S the sum of d^2 over the nS springs, which acts as one tension k3 S / nS shared by all springs. They agree when all springs are stretched alike.
1717 - Curvature (toward gongs and cymbals): rise h (m, default 0 = straight) of a stress-free rest shape w0, a half sine between two supports, 1 - cos(pi x / 2L) for a cantilever (flat at the clamp). Stretching is measured from it: d^2 becomes e = d^2 + 2 D0 d (D0 the spring's rest offset) in both forms; bending is unchanged (linear). This adds linear stiffness and a quadratic nonlinearity (softening) to the quartic one (hardening). The animation draws the beam at w0 + u, with the rest shape as the faint line through the rails, both scaled to fit the largest |w0 + u| (so at rest the arch fills the upper half; small vibrations on a large arch look small). The nonlinear scale in the derived values counts only the quartic part. Checked (Node): h = 0 bit-identical to before (all ends, both forms, k = 0 and 10, 50000 damped steps); forces = -grad V to 5e-11 and the linear stiffness = the Hessian at rest to 7e-9 for h = 0.05 and 0.3; energy error <= 7e-5 over 200 T1 up to 100 J.
18-- Pitch glide (mode-1 frequency from zero crossings, clamped, local, τ 540 s, cents from the linear mode 1): straight beam hardening, pitch falls as the note decays (mode-1 strike 2.51 J: +1237 c at 2 J, +352 at 0.1 J, +25 at 4 mJ). Curvature flips it: h = 0.1 m hardening at high energy, softening below ~0.1 J (pitch falls, then rises slightly); h = 0.2 m +302 c at 2 J, -232 at 0.5 J, -27 at 20 mJ (mode 2 at 1.98 x mode 1, near 2:1); h = 0.3 m softening throughout (-358 c at 2 J, -49 at 0.1 J), pitch rises as the note decays. Kirchhoff is similar. Default mass-2 strike, h = 0.3: lowest line rises from ~78 to 100 Hz over ~3 s of sound (17 s note, x219: mode 1 is 0.457 Hz, so the same τ covers fewer periods per second of sound), with steady difference tones (e.g. mode 2 - mode 1 at 79 Hz) from the quadratic term.
18+- Pitch glide (mode-1 frequency from zero crossings, clamped, local, τ 540 s, cents from the linear mode 1): straight beam hardening, pitch falls as the note decays (mode-1 strike 2.51 J: +1237 c at 2 J, +352 at 0.1 J, +25 at 4 mJ). Curvature flips it: h = 0.1 m hardening at high energy, softening below ~0.1 J (pitch falls, then rises slightly); h = 0.2 m +302 c at 2 J, -232 at 0.5 J, -27 at 20 mJ (mode 2 at 1.98 x mode 1, but antisymmetric, so by symmetry it has no quadratic coupling to mode 1: not a 1:2 resonance); h = 0.3 m softening throughout (-358 c at 2 J, -49 at 0.1 J), pitch rises as the note decays. Kirchhoff is similar. Default mass-2 strike, h = 0.3: lowest line rises from ~78 to 100 Hz over ~3 s of sound (17 s note, x219: mode 1 is 0.457 Hz, so the same τ covers fewer periods per second of sound), with steady difference tones (e.g. mode 2 - mode 1 at 79 Hz) from the quadratic term.
1919 - Normal modes of the linear beam from the stiffness matrix (assembled from the force routine) by Jacobi; mode energies as in fput-chain. Checked: pinned frequencies equal the closed form omega_n^2 = (4 k s^2 + 16 κ s^4) / m, s = sin(pi n / (2 (N + 1))), to 1e-14; N = 32 ratios clamped 2.75, 5.35, 8.78 (beam theory 2.76, 5.40, 8.93), cantilever 6.26, 17.48, 34.09 (6.27, 17.55, 34.39). At N = 8: clamped 1, 2.59, 4.68, 6.97, top 13.2; pinned 1, 3.88, 8.29, top 32.2; cantilever 1, 6.13, 16.4, top 75.
2020 - Nonlinear scale: the energy of mode 1 alone at the amplitude where its stretching energy equals its linear energy (N = 8, local: clamped 176 mJ, pinned 8.7 mJ, cantilever 1.1 mJ).
2121 - Damping: drag proportional to mass, every mode's amplitude decays as exp(-t / τ); applied exactly in half steps around each Verlet step, dissipated energy in W.