 Solution

Use the Euler-Bernoulli beam equation for a uniform slender elastic beam, small deflection and small slope, neglecting axial tension, shear deformation, rotational inertia of a cross-section and damping. Let $E$ be Young's modulus, $I={\int }_A{z}^{2}dA$ the second moment of area about the neutral axis, and $\rho$ the mass per unit length. If density instead denotes mass per volume, replace this $\rho$ by ${\rho }_{\mathrm{bulk}}A$. The bending stiffness is $EI$; the small-slope curvature is $y_{xx}$. The bending and kinetic energies areTwice integrating the variation of $U$ by parts gives the restoring load $-EI{y}_{xxxx}$ per unit length. Newton's second law therefore givesAt a pivot, displacement is fixed but slope is free; zero end bending moment gives $EI{y}_{xx}=0$. This also follows from the natural boundary term $EI{y}_{xx}\delta y_x$ in the variation.

With separation of variables, $y=XT$ givesChoose the time normalization $T(0)=1$. Initial rest removes the sine time term, so $T(t)=\cos \omega t$. The spatial characteristic roots are $\pm ik,\pm k$, givingThe two conditions at $x=0$ give $B+D=0$ and $-B+D=0$, hence $B=D=0$. At $L$, addition and subtraction of $X(L)=X^{\prime \prime }(L)=0$ give $C\sinh kL=0$ and $A\sin kL=0$. Nontrivial positive-frequency modes therefore have $C=0$ and $k_n=n\pi /L$. Thus the simply supported beam normal modes areThe zero-frequency cubic solution also vanishes under these four endpoint conditions, so $n=0$ is not an additional mode.

For the specified release, small slope requires $|a|L\ll 1$. Expand the initial deflection in the Fourier sine series. Its coefficients areWriting ${\omega }_1={\pi }^2\sqrt{EI/\rho }/L^2$, the resulting motion isOnly odd normal modes occur because the initial profile is symmetric about the midpoint, and their amplitudes decrease as $n^{-3}$. The frequencies grow as $n^2$, so the motion is dispersive and is not a single shape multiplied by one cosine. Nevertheless the ideal undamped solution is exactly periodic with period $2\pi /{\omega }_1$. At $t=\pi /(2{\omega }_1)$ all odd-mode cosines vanish, so the bar is instantaneously straight; at $t=\pi /{\omega }_1$ it has the opposite initial deflection. The first mode dominates displacement, with third-mode amplitude only $1/27$ of the first.

There is initial curvature incompatibility at a moment-free beam end: the prescribed parabola has $y_{xx}(x,0)=2a$, whereas a moment-free pivot requires zero curvature at an endpoint. Therefore the data do not define a classical solution with curvature continuous through the two end-time corners. They do define a finite-energy weak solution after release of the forces or end moments that prepared the initial shape. The sine expansion represents that solution: its displacement reproduces the initial parabola, each mode satisfies the released endpoint conditions, and ${\sum }_n{n}^4|A_n{|}^2<\infty$ gives finite bending energy. One must not differentiate the initial series twice at its endpoints and assign the interior parabola's curvature there. This qualification leaves the modal displacement prediction intact.