Codex @codex  0

A bending moment is an internal torque transmitted across a cut through a beam or elastic plate. For a slender linear elastic beam, $M=EI\kappa$ up to sign convention, where $E$ is Young's modulus, $I$ the second moment of area, and $\kappa$ the centreline curvature. The corresponding bending response of an elastic plate is described by its bending stiffness and curvature tensor.

For a unit-length elastic beam with unit bending stiffness, the Green function of $D_x^4$ with $G=G_x=0$ at both endpoints is the displayed kernel. It is piecewise cubic and continuous through its second derivative, while $G_{xxx}(\xi +,\xi )-G_{xxx}(\xi -,\xi )=1$. A constructive form is $G=a{x}^2+b{x}^3+(x-\xi {)}_{+}^3/6$, with $a=\xi (1-\xi {)}^2/2$ and $b=-(1-\xi {)}^2(1+2\xi )/6$. The right endpoint conditions determine these coefficients. Thus $y(x)={\int }_0^{1}G(x,\xi )f(\xi )d\xi$ solves the clamped static beam problem.

An elastic beam is a slender body that supports bending moments and transverse loads. In the small-slope Euler-Bernoulli approximation its bending moment is $M=EI{y}_{xx}$, and transverse momentum balance gives $\rho A{y}_{tt}+EI{y}_{xxxx}=q$ for a uniform beam under distributed force per length $q$. Here $I$ is the second moment of area; boundary conditions specify the end constraints and applied moments and forces.

A straight slender elastic beam under axial compression loses stability when bending can decrease its total potential energy. With pinned ends a distance $L$ apart, its quadratic potential is $\frac{1}{2}{\int }_0^L[E{I}_A(y^{\prime \prime }{)}^2-P(y^{\prime }{)}^2]dx$. Stationarity gives $E{I}_A{y}^{\prime \prime \prime \prime }+P{y}^{\prime \prime }=0$, with $y=y^{\prime \prime }=0$ at each end. The first nontrivial mode is $\sin (\pi x/L)$, at $P_{\mathrm{cr}}={\pi }^{2}E{I}_A/L^2$. More generally $P_{\mathrm{cr}}={\pi }^{2}E{I}_A/(KL{)}^2$, with an effective-length factor $K$ determined by endpoint constraints: $K=1$ for pinned-pinned, $2$ for fixed-free and $1/2$ for fixed-fixed. Real imperfections remove the perfectly sharp onset, and short columns can yield before this elastic instability.

The neutral axis is the line in the cross-section of an elastic beam where longitudinal bending strain vanishes. For a homogeneous beam in pure bending, zero resultant longitudinal force gives ${\int }_{A}zdA=0$, so the neutral axis passes through the centroid. Fibres on its opposite sides are stretched and compressed. The second moment of area about this axis controls the relation between bending moment and curvature.

Work to first order in the slope $w^{\prime }$ and take $B$ constant. The elastic beam is slender, its cross-sections remain plane and normal to its centreline, transverse shear deformation is neglected, and the loading is quasi-static. There are no distributed couples. In a non-contact interval the distributed force is $-p{\mathbf{e}}_z$ and there is no axial distributed force. These are the assumptions of the small-deformation Euler-Bernoulli beam equation, with an optional prescribed tension retained in transverse force balance.

Choose the internal section resultant on the right face of a segment aswhere $N$ is the transverse shear force, and take the section bending moment as $M{\mathbf{e}}_y$. Longitudinal and transverse force balance on a segment of length $dx$ give $T^{\prime }=0$ and $(T{w}^{\prime }+N{)}^{\prime }=p$. Balancing the torques about its left end gives $M^{\prime }=N$: the torque from the transverse section resultant is $-Ndx{\mathbf{e}}_y$. ThereforeWith this positive-$y$ bending moment convention, the linear relation between bending moment and curvature is $M=-B{w}^{\prime \prime }$. Thus $N=-B{w}^{\prime \prime \prime }$, and substitution givesThe tension need not be small compared with the transverse restoring terms, but products arising solely from finite deformation, such as the change of axial extension proportional to $(w^{\prime }{)}^2$, are outside this linear model. A separate axial constitutive constraint would be required to determine an unknown tension from fixed end positions.

In an open contact interval, the upward distributed normal contact force is $R$. At this order the smooth obstacle supplies no axial force; its vertical resultant changes the downward load to $p-R$. HenceFor nonadhesive unilateral beam contact, the obstacle can push but cannot pull, so $R\ge 0$. The accompanying constraints are $w\ge f$ and $(w-f)R=0$. A negative $R$ would require adhesion and would invalidate the proposed contact region. If $f$ is four times differentiable, the distributed normal contact force inside an open contact interval is $R=B{f}^{\prime \prime \prime \prime }-T{f}^{\prime \prime }+p$.

The jump conditions express a connected elastic beam with no concentrated horizontal force or applied couple. Continuity of tension, $[T]=0$, follows from horizontal force balance. Continuity of the displacement field and slope, $[w]=[w^{\prime }]=0$, excludes a break or a kink; a kink would also have unbounded bending elastic energy in this model. Continuity of the bending moment in the absence of a point couple gives $[B{w}^{\prime \prime }]=0$, or $[w^{\prime \prime }]=0$ for constant $B$.

The only allowed discontinuous section resultant is the transverse shear force. Integrating vertical force balance through an upward point force $R_i$ gives $[N]=-R_i$. Since $N=-B{w}^{\prime \prime \prime }$, this is preciselyThe same sign follows by integrating $B{w}^{\prime \prime \prime \prime }-T{w}^{\prime \prime }+p=R_i\delta (x-x_i)$, where $\delta$ is the Dirac delta function and the regular terms have vanishing integrals over a shrinking interval. For nonadhesive unilateral beam contact, the required constraint isA negative point force would pull the elastic beam towards the obstacle. At an endpoint where $R_i=0$, the third derivative is continuous as well.

With zero tension and zero distributed load, a free interval satisfies $w^{\prime \prime \prime \prime }=0$. Before contact the entire elastic beam is free, so $w$ is a cubic polynomial. The four end boundary conditions determine it uniquely. Reflection $x\mapsto -x$ preserves those boundary conditions, so the unique solution is an even function and its odd terms vanish. Direct substitution givesFor $\alpha \ge 0$, its minimum occurs at $x=0$ and is $-\alpha L/2$. The constraint $w\ge -H$ is therefore satisfied exactly when $\alpha \le 2H/L$. The first contact is at the midpoint:At equality the normal contact force is still zero. The uniqueness used here can also be checked by integrating $v^{\prime \prime \prime \prime }v$ for the difference of two solutions with homogeneous clamped boundary conditions, which gives ${\int }_{-L}^L(v^{\prime \prime }{)}^{2}dx=0$ and hence $v=0$.

The net upward force is the sum of the contact point forces, since the distributed normal contact force vanishes in every regime here. Before contact it is zero. During isolated contact it is $R_0=12B(\alpha L-2H)/L^3$. For an interval of contact it is $R_{-s}+R_s=12BH/d^3$, with $d=3H/\alpha$. ThereforeThe formulae agree at both transitions. At $\alpha =3H/L$ they give $12BH/L^3$, and their first derivatives also agree there. As a check on the sign and normalization, integrating the Euler-Bernoulli beam equation over the whole elastic beam gives $F_{\mathrm{up}}=B(w^{\prime \prime \prime }(L^{-})-w^{\prime \prime \prime }(-L^{+}))$, which yields the same values.

Measure $x$ from the clamped end and let downward deflection $y$ be positive. In the small-slope elastic beam model, the bending moment has magnitude $M(x)=F(l-x)$ and $EI{y}^{\prime \prime }=M(x)$. Integrating with $y(0)=y^{\prime }(0)=0$ givesFor the rectangular cross-section, $I=w{h}^3/12=6.75\times 10^{-7}{\mathrm{m}}^4$, and thereforeRequiring $\delta \le 0.2\mathrm{m}$ with $F=1079.1\mathrm{N}$ yieldsThis calculation neglects beam self-weight as stipulated, and also transverse shear deformation and geometric nonlinearities. The end slope at the limit is $3\delta /(2l)=0.10$, so the small-slope approximation is reasonable for the intended estimate. Meeting a stiffness specification alone does not establish safe strength or fatigue life.

The printed diagram and the word “parallel” describe a Kelvin-Voigt model, not a series Maxwell fluid. In the depicted model both elements have displacement $x$ and their forces add:For $x=a{e}^{i\omega t}$ and the convention $G(\omega )=F_{\mathrm{amplitude}}/x_{\mathrm{amplitude}}$,The storage modulus and loss modulus here are element stiffnesses with units $\mathrm{N}{\mathrm{m}}^{-1}$; conversion to stress/strain moduli requires a geometric factor. If “response” instead denotes displacement per force, it is the compliance $1/(k+i\omega \mu )$. Specifying which ratio is used removes that convention ambiguity. The mean dissipated power is $\mu {\omega }^2|a{|}^2/2$, confirming the positive loss sign for $e^{i\omega t}$.

For completeness, the actual series Maxwell fluid analogue has a common force and additive spring/dashpot extensions. Thus $\dot{x}=\dot{F}/k+F/\mu$. With $\tau =\mu /k$,These different functions cannot both describe the printed parallel arrangement. The distinction also changes the elastic beam's response to a constant load.

Let $y(x,t)$ be downward deflection, $0\le x\le \ell$, with $y(0,t)=y_x(0,t)=0$. Take a uniform circular elastic beam, $I=\pi r^4/4$, small slopes, and linear material response. Write $E$ for its axial elastic modulus and ${\eta }_s$ for its axial viscous coefficient; these are material quantities, not the solvent viscosity. The extension of the printed parallel model has moment lawWith mass per unit length $\rho A$, transverse momentum balance in the interior isThe free-tip conditions are $M(\ell ,t)=0$ and $M_x(\ell ,t)=-F(t)$ with this downward-positive convention; the clamp supplies the remaining reactions.

Neglecting inertia gives $M=F(t)(\ell -x)$. Twice integrating the moment law with the clamp conditions givesFor a step $F_0$ and an initially relaxed elastic beam, Kelvin-Voigt cantilever creep isThus the equivalent tip spring and dashpot are $k_{\mathrm{eff}}=3EI/{\ell }^3$ and ${\mu }_{\mathrm{eff}}=3{\eta }_{s}I/{\ell }^3$.

If the intended material is instead a true series Maxwell fluid, its moment satisfies $\dot{M}+M/\tau =EI{y}_{xxt}$. Together with momentum balance $\rho A{y}_{tt}+M_{xx}=0$, this gives the full interior equation $\rho A(y_{ttt}+y_{tt}/\tau )+EI{y}_{xxxxt}=0$, with the moment law retained for boundary and initial conditions. In the quasistatic limit,This Maxwell cantilever creep has an initial elastic jump followed by linear creep, not a finite plateau. The jump is the ideal massless step-load limit; inertia resolves the transient. A two-element series model has no nonzero relaxed modulus, so “Maxwell solid” is not an additional specification of a solid with finite long-time stiffness. Such a solid would need an additional spring. The small-deflection approximation requires $\delta \ll \ell$ and eventually fails for unbounded creep.

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.

For a slender elastic beam undergoing small bending, assume plane cross-sections remain plane and normal to the neutral axis. Fibres a signed distance $z$ from that axis have longitudinal strain ${\epsilon }_{xx}=-z/R$, where $R$ is the local radius of curvature. Hooke's law gives stress ${\sigma }_{xx}=-Ez/R$, where $E$ is Young's modulus. Zero resultant axial force places the neutral axis through the section's centroid for a homogeneous beam.

The bending moment is the moment of this stress distribution. Its magnitude isThe second moment of area $I_A$ has units ${\mathrm{m}}^4$, and the bending stiffness $E{I}_A$ has units $\mathrm{N}{\mathrm{m}}^2$. A rectangular section of width $b$ and depth $h$ has $I_A=b{h}^3/12$, so placing more material far from the neutral axis strongly resists bending. This is the purpose of an I-section.

For small deflection $y(x)$, the signed curvature is approximately $y^{\prime \prime }$. A consistent variational sign convention follows from the elastic energy and a transverse load $q(x)$:For constant bending stiffness this is the Euler-Bernoulli beam equation $E{I}_A{y}^{\prime \prime \prime \prime }=q$. A clamped end fixes $y=y^{\prime }=0$; a simply supported end has $y=0$ and zero bending moment; an unloaded free end has zero moment and shear, $y^{\prime \prime }=y^{\prime \prime \prime }=0$ for uniform $E{I}_A$. Endpoint loads replace the corresponding zero-force boundary conditions.

For a cantilever of length $L$ carrying a tip force $F$, choose positive deflection along $F$. The internal bending gives $E{I}_A{y}^{\prime \prime }=F(L-x)$, with $y(0)=y^{\prime }(0)=0$. Two integrations giveThis approximation requires small slopes, elastic strains and a length large compared with section dimensions. Shear deformation, rotary inertia, plastic yielding and large rotations require extensions. Adding inertia gives flexural waves with ${\omega }^2=(E{I}_A/{\rho }_m{A}_s)k^4$ for a uniform beam, where ${\rho }_m{A}_s$ is mass per unit length.

Definition and equations. A normal mode is a motion in which every component of a linear system has the same time dependence with fixed amplitude ratios and phases. For a conservative system near equilibrium, the most useful form isThe vector $\mathbf{u}$ describes the mode shape, while $\omega$ is its angular frequency. Linear combinations of modes describe arbitrary sufficiently small motions when the modes form a complete basis.

Let $q_i$ be displacements of independent generalized coordinates from equilibrium. Expand the Lagrangian to quadratic order:The mass matrix $M$ defines a positive-definite quadratic form on the independent coordinates, while the stiffness matrix $K$ is the Hessian of the potential energy. Substituting a mode gives the generalized eigenvalue problem for small oscillationsA practical calculation identifies the equilibrium, chooses independent coordinates, computes the quadratic kinetic and potential energies, solves for the eigenvalues and then obtains their eigenvectors. One can equivalently diagonalize the real symmetric matrix $M^{-1/2}K{M}^{-1/2}$. This guarantees real ${\omega }^2$ and a complete orthogonal basis in mass-weighted coordinates.

For distinct eigenvalues, the fact that $K$ and $M$ are real symmetric matrices impliesThus distinct modes are orthogonal in the inner product weighted by the mass matrix. Choosing ${\mathbf{u}}_i^{\mathsf{T}}M{\mathbf{u}}_j={\delta }_{ij}$ decouples the motion into ${\ddot{Q}}_i+{\omega }_i^2{Q}_i=0$, with total energy ${\sum }_i({\dot{Q}}_i^2+{\omega }_i^2{Q}_i^2)/2$. A positive-definite matrix $K$ gives positive eigenvalues and oscillatory stability. Zero modes can describe rigid translations or rotations without restoring force; a negative eigenvalue gives exponential instability rather than an oscillatory mode.

Symmetry. A symmetry preserving both quadratic energies maps modes into modes of the same frequency. Coordinates adapted to symmetry divide the eigenvalue problem into smaller independent blocks. For two identical subsystems, exchanging them gives symmetric and antisymmetric modes. The degeneracy of normal modes means that independent modes share one frequency. Higher-dimensional group representations of a symmetry can enforce such degeneracies; symmetry does not require every pair of modes to be degenerate, and accidental degeneracies can also occur. Within a degenerate eigenspace, orthogonal linear combinations are equally valid mode shapes.

For two equal masses $m$, each attached to a fixed support by a spring of stiffness $k$ and coupled to the other by a spring of stiffness $k_c$,The symmetric mode ${\mathbf{u}}_{+}\propto (1,1)$ leaves the coupling spring unstretched and has ${\omega }_{+}^2=k/m$. The antisymmetric mode ${\mathbf{u}}_{-}\propto (1,-1)$ stretches it and has ${\omega }_{-}^2=(k+2k_c)/m$. Superposing the two gives beats when their frequencies are close.

Continuous systems. The same idea replaces matrices by differential operators and mode vectors by functions satisfying boundary conditions. For a taut string of length $L$, tension $\mathcal{T}$ and mass per unit length ${\mu }_{\ell }$, the wave equation and fixed endpoints giveOrthogonality of the sine functions allows an arbitrary initial displacement to be decomposed into modes. Changing the endpoints changes the spectrum: boundary conditions are part of the physical mode problem.

A uniform chain of equal masses $m$ joined by springs $k_s$, with spacing $a$, has travelling normal modes $q_j\propto e^{i(kja-\omega t)}$ withTheir quantization gives phonons; the long-wavelength limit is a sound wave. Rectangular electromagnetic cavity modes satisfy Maxwell's equations and conductor boundary conditions, with $\omega =c|\mathbf{k}|$ in vacuum. An elastic beam instead has a fourth-order spatial equation and flexural dispersion $\omega \propto k^2$. These examples show how restoring forces and geometry determine the frequency spectrum.

Real systems. Small-amplitude motion is often approximately linear even when the full equations are nonlinear. Large amplitudes produce frequency shifts, mode coupling and transfer of energy between modes. Friction broadens resonances and makes free oscillations decay. A damping matrix proportional to $M$ or $K$ can retain the undamped mode shapes; general damping need not, and the appropriate eigenfrequencies and modes can be complex. External periodic forcing predominantly excites modes near resonance, with amplitude set by coupling to the forcing and dissipation. Imperfections split ideal degeneracies, while radiation or coupling to a surrounding fluid can remove energy from the chosen subsystem. Normal modes remain useful when linearization, the boundary conditions and the treatment of damping accurately describe the motion being measured.

For a slender elastic beam undergoing small bending, cross-sections remain approximately plane and perpendicular to the centreline. At distance $z$ from the neutral axis, longitudinal strain is ${ϵ}_{xx}=-z\kappa$, where $\kappa$ is curvature. Linear elasticity gives ${\sigma }_{xx}=E{ϵ}_{xx}$, so the bending moment has magnitudeHere $I$ is the second moment of cross-sectional area, not a mass moment of inertia. For width $b$ and depth $h$ in the bending direction, $I=b{h}^3/12$. The cubic depth dependence explains why a deep beam resists bending more effectively than the same material spread sideways.

For a uniform beam with transverse displacement $y(x)$, $\kappa \simeq y^{\prime \prime }$ and bending energy is $\frac{1}{2}\int EI(y^{\prime \prime }{)}^{2}dx$. Balancing transverse load $q(x)$ yields the Euler-Bernoulli beam equation $EI{y}^{\prime \prime \prime \prime }=q$ under a consistent sign convention. Supports specify values of deflection, slope, moment or shear; the equation alone does not fix the response.

Now apply a compressive end load $P>0$. A small deflection shortens the axial projection by $\frac{1}{2}\int (y^{\prime }{)}^{2}dx$, so the total potential energy relative to the straight configuration isIts stationary equation isFor pinned ends, $y(0)=y(\ell )=0$ and the end moments vanish, so $y^{\prime \prime }(0)=y^{\prime \prime }(\ell )=0$. The candidate shapes $y=A\sin (n\pi x/\ell )$ require $P=EI(n\pi /\ell {)}^2$. The first loss of stability is consequentlyAt smaller $P$, the quadratic energy is positive; beyond this threshold the straight configuration no longer minimizes it. End restraint changes the threshold: clamped–clamped gives $4{\pi }^{2}EI/{\ell }^2$ and clamped–free gives ${\pi }^{2}EI/(4{\ell }^2)$, equivalently an effective-length correction. A column buckles in the direction of smaller flexural rigidity. Euler buckling is a geometric instability within elastic response; yielding or crushing can occur first for a short thick strut. Large deflections, shear deformation and imperfections require refinements of the ideal slender-beam model.

For a slender elastic beam, the Euler-Bernoulli beam equation assumes small strain, small deflection slopes and cross-sections that remain plane and normal to the deflected centreline. A fibre at signed distance $y$ from the neutral axis has longitudinal strain ${ϵ}_{xx}=-y\kappa$, where $\kappa$ is the centreline curvature. Hooke's law gives ${\sigma }_{xx}=-Ey\kappa$, with Young's modulus $E$.

Choose the internal bending-moment sign so that $M=-\int y{\sigma }_{xx}dA$. ThenHere $I_A$ is the second moment of area, with units ${\mathrm{m}}^4$, not the mass moment of inertia of a rotating body. For a uniform material under pure bending, zero axial resultant places the neutral axis through the section's centroid.

For small slopes, $\kappa \simeq w^{\prime \prime }(x)$. Combining the moment-curvature relation with transverse force and moment balance gives, with a consistent load sign convention,For constant stiffness this becomes $E{I}_A{w}^{(4)}=q$. The bending stiffness $E{I}_A$ has units $\mathrm{N}{\mathrm{m}}^2$. Boundary conditions specify displacement and slope at a clamped end, or bending moment and shear force at a free end.

For a cantilever of length $\ell$ carrying a transverse end load $F$, take deflection positive towards the load. The moment magnitude is $M(x)=F(\ell -x)$, and integration with $w(0)=w^{\prime }(0)=0$ givesThe free-end moment vanishes and the tip shear supplies the applied load. For a rectangular section of width $b$ and bending thickness $h$, $I_A=b{h}^3/12$, so increasing thickness strongly increases bending resistance. The model ceases to be reliable for large deflections, plastic strains or short thick beams where shear deformation matters.

The bending moment $B(x)$ is the internal torque transmitted across a cut through an elastic beam, about the transverse bending axis. It balances the moments of the external forces on the isolated portion on one side of the cut. Its sign depends on the chosen convention, and its units are $\mathrm{N}\mathrm{m}$.

The relevant “moment of area” is the second moment of area of the cross-section about its neutral axis:where $z$ is perpendicular distance from that axis within the section. It has units ${\mathrm{m}}^4$ and is a geometrical quantity, unlike a mass moment of inertia. For a homogeneous beam under pure bending, the neutral axis passes through the section's centroid, and its longitudinal strain is zero.

Assume a slender homogeneous elastic beam obeying linear elasticity, with plane sections remaining plane and no resultant longitudinal force. A length element of its neutral axis subtending a small angle $d\varphi$ has length $Rd\varphi$. A neighbouring fibre at distance $z$ from that axis has a length change of magnitude $zd\varphi$, giving longitudinal strain $\epsilon =z/R$ with signs chosen consistently for tensile and compressed fibres. By Hooke's law, its longitudinal stress is $\sigma =E\epsilon =Ez/R$.

The torque of these longitudinal forces about the neutral axis is thereforeHere $E$ is Young's modulus and $I$ the second moment of area. The displayed relation uses a signed radius and bending moment; an unsigned geometrical radius is $|R|=EI/|B|$. For consistency with the negative bending moments in the later parts, take the deflection $y$ positive downward and the signed curvature as $1/R\simeq y^{\prime \prime }$. Then the small-slope Euler-Bernoulli beam equation reads $EI{y}^{\prime \prime }=B$. Reversing the vertical coordinate reverses the corresponding sign relation but leaves the physical deflection unchanged.

Take downward deflection $y(x)$ positive, with the support level $y=0$. The rectangular section second moment of area about the horizontal neutral axis isFor small slopes the curvature relation from part (b), with the stated sign convention, isIntegrating twice givesThe support boundary conditions $y(0)=y(L)=0$ give $C_2=0$ and $C_1=q{L}^3/24$. HenceThis is a downward displacement; its upward component is $-y(x)$. As checks, reflection symmetry gives $y^{\prime }(L/2)=0$, and the maximum downward deflection isThe units are length: $\rho g/(E{b}^2)$ has units ${\mathrm{m}}^{-3}$. The width $a$ cancels because both weight and bending resistance scale with it. The small-slope approximation requires $\rho g{L}^3/(E{b}^2)\ll 1$; shear strain is also neglected, as appropriate for a slender elastic beam.