In applied mechanics, bending (also known as flexure) characterizes the behavior of a slender structural element subjected to an external Structural load applied perpendicularly to a longitudinal axis of the element.
The structural element is assumed to be such that at least one of its dimensions is a small fraction, typically 1/10 or less, of the other two.Boresi, A. P. and Schmidt, R. J. and Sidebottom, O. M., 1993, Advanced mechanics of materials, John Wiley and Sons, New York. When the length is considerably longer than the width and the thickness, the element is called a beam. For example, a closet rod sagging under the weight of clothes on is an example of a beam experiencing bending. On the other hand, a shell is a structure of any geometric form where the length and the width are of the same order of magnitude but the thickness of the structure (known as the 'wall') is considerably smaller. A large diameter, but thin-walled, short tube supported at its ends and loaded laterally is an example of a shell experiencing bending.
In the absence of a qualifier, the term bending is ambiguous because bending can occur locally in all objects. Therefore, to make the usage of the term more precise, engineers refer to a specific object such as; the bending of rods, the bending of beams, the bending of plates,Timoshenko, S. and Woinowsky-Krieger, S., 1959, Theory of plates and shells, McGraw-Hill. the Plate theoryLibai, A. and Simmonds, J. G., 1998, The nonlinear theory of elastic shells, Cambridge University Press. and so on.
These last two forces form a couple or moment as they are equal in magnitude and opposite in direction. This bending moment resists the sagging deformation characteristic of a beam experiencing bending. The stress distribution in a beam can be predicted quite accurately when some simplifying assumptions are used.
Consider beams where the following are true:
In this case, the equation describing beam deflection () can be approximated as:
where the second derivative of its deflected shape with respect to is interpreted as its curvature, is the Young's modulus, is the area moment of inertia of the cross-section, and is the internal bending moment in the beam.
If, in addition, the beam is homogeneous along its length as well, and not tapered (i.e. constant cross section), and deflects under an applied transverse load , it can be shown that:
After a solution for the displacement of the beam has been obtained, the bending moment () and shear force () in the beam can be calculated using the relations
Simple beam bending is often analyzed with the Euler–Bernoulli beam equation. The conditions for using simple bending theory are:Shigley J, "Mechanical Engineering Design", p44, International Edition, pub McGraw Hill, 1986,
Compressive and tensile forces develop in the direction of the beam axis under bending loads. These forces induce stresses on the beam. The maximum compressive stress is found at the uppermost edge of the beam while the maximum tensile stress is located at the lower edge of the beam. Since the stresses between these two opposing vary , there therefore exists a point on the linear path between them where there is no bending stress. The locus of these points is the neutral axis. Because of this area with no stress and the adjacent areas with low stress, using uniform cross section beams in bending is not a particularly efficient means of supporting a load as it does not use the full capacity of the beam until it is on the brink of collapse. Wide-flange beams (i-Beam) and truss effectively address this inefficiency as they minimize the amount of material in this under-stressed region.
The classic formula for determining the bending stress in a beam under simple bending is:Gere, J. M. and Timoshenko, S.P., 1997, Mechanics of Materials, PWS Publishing Company.
where are the coordinates of a point on the cross section at which the stress is to be determined as shown to the right, and are the bending moments about the y and z centroid axes, and are the second moments of area (distinct from moments of inertia) about the y and z axes, and is the product of moments of area. Using this equation it is possible to calculate the bending stress at any point on the beam cross section regardless of moment orientation or cross-sectional shape. Note that do not change from one point to another on the cross section.
Large bending considerations should be implemented when the bending radius is smaller than ten section heights h:
With those assumptions the stress in large bending is calculated as:
where
The equation for the bending of a linear elastic, isotropic, homogeneous beam of constant cross-section under these assumptions isRosinger, H. E. and Ritchie, I. G., 1977, On Timoshenko's correction for shear in vibrating isotropic beams, J. Phys. D: Appl. Phys., vol. 10, pp. 1461–1466.
where is the polar moment of inertia of the cross-section, is the mass per unit length of the beam, is the density of the beam, is the cross-sectional area, is the shear modulus, and is a shear correction factor. For materials with Poisson's ratios () close to 0.3, the shear correction factor are approximately
These assumptions imply that
The strain-displacement relations are
The equilibrium equations are
In terms of displacements, the equilibrium equations for an isotropic, linear elastic plate in the absence of external load can be written as
The strain-displacement relations that result from these assumptions are
The equilibrium equations are
The figures below show some vibrational modes of a circular plate.
Quasi-static bending of beams
Euler–Bernoulli bending theory
EI~\cfrac{\mathrm{d}^4 w(x)}{\mathrm{d} x^4} = q(x)
This is the Euler–Bernoulli equation for beam bending.
M(x) = -EI~\cfrac{\mathrm{d}^2 w}{\mathrm{d} x^2} ~;~~ Q(x) = \cfrac{\mathrm{d}M}{\mathrm{d}x}.
where
Extensions of Euler-Bernoulli beam bending theory
Plastic bending
Complex or asymmetrical bending
Large bending deformation
The general solution of the above equation is
\hat{w} = A_1\cosh(\beta x) + A_2\sinh(\beta x) + A_3\cos(\beta x) + A_4\sin(\beta x)
where are constants and
Timoshenko–Rayleigh theory
\begin{align}
& EI~\frac{\partial^4 w}{\partial x^4} + m~\frac{\partial^2 w}{\partial t^2} - \left(J + \frac{E I m}{k A G}\right)\frac{\partial^4 w}{\partial x^2~\partial t^2} + \frac{J m}{k A G}~\frac{\partial^4 w}{\partial t^4} \\6pt
= {} & q(x,t) + \frac{J}{k A G}~\frac{\partial^2 q}{\partial t^2} - \frac{EI}{k A G}~\frac{\partial^2 q}{\partial x^2}
\end{align}
\begin{align}
k &= \frac{5 + 5\nu}{6 + 5\nu} \quad \text{rectangular cross-section}\\[6pt]
&= \frac{6 + 12\nu + 6\nu^2}{7 + 12\nu + 4\nu^2} \quad \text{circular cross-section}
\end{align}
Free vibrations
EI~\cfrac{\mathrm{d}^4 \hat{w}}{\mathrm{d} x^4} + m\omega^2\left(\cfrac{J}{m} + \cfrac{E I}{k A G}\right)\cfrac{\mathrm{d}^2 \hat{w}}{\mathrm{d} x^2} + m\omega^2\left(\cfrac{\omega^2 J}{k A G}-1\right)~\hat{w} = 0
This equation can be solved by noting that all the derivatives of must have the same form to cancel out and hence as solution of the form may be expected. This observation leads to the characteristic equation
\alpha~k^4 + \beta~k^2 + \gamma = 0 ~;~~ \alpha := EI ~,~~ \beta := m\omega^2\left(\cfrac{J}{m} + \cfrac{E I}{k A G}\right) ~,~~ \gamma := m\omega^2\left(\cfrac{\omega^2 J}{k A G}-1\right)
The solutions of this quartic equation are
k_1 = +\sqrt{z_+} ~,~~ k_2 = -\sqrt{z_+} ~,~~ k_3 = +\sqrt{z_-} ~,~~ k_4 = -\sqrt{z_-}
where
z_+ := \cfrac{-\beta + \sqrt{\beta^2 - 4\alpha\gamma}}{2\alpha} ~,~~
z_-:= \cfrac{-\beta - \sqrt{\beta^2 - 4\alpha\gamma}}{2\alpha}
The general solution of the Timoshenko-Rayleigh beam equation for free vibrations can then be written as
\hat{w} = A_1~e^{k_1 x} + A_2~e^{-k_1 x} + A_3~e^{k_3 x} + A_4~e^{-k_3 x}
Quasistatic bending of plates
Kirchhoff–Love theory of plates
\begin{align}
u_\alpha(\mathbf{x}) & = - x_3~\frac{\partial w^0}{\partial x_\alpha}
= - x_3~w^0_{,\alpha} ~;~~\alpha=1,2 \\
u_3(\mathbf{x}) & = w^0(x_1, x_2)
\end{align}
where is the displacement of a point in the plate and is the displacement of the mid-surface.
\begin{align}
\varepsilon_{\alpha\beta} & =
- x_3~w^0_{,\alpha\beta} \\
\varepsilon_{\alpha 3} & = 0 \\
\varepsilon_{33} & = 0
\end{align}
M_{\alpha\beta,\alpha\beta} + q(x) = 0 ~;~~ M_{\alpha\beta} := \int_{-h}^h x_3~\sigma_{\alpha\beta}~dx_3
where is an applied load normal to the surface of the plate.
w^0_{,1111} + 2~w^0_{,1212} + w^0_{,2222} = 0
In direct tensor notation,
\nabla^2\nabla^2 w = 0
Mindlin–Reissner theory of plates
\begin{align}
u_\alpha(\mathbf{x}) & = - x_3~\varphi_\alpha ~;~~\alpha=1,2 \\
u_3(\mathbf{x}) & = w^0(x_1, x_2)
\end{align}
where are the rotations of the normal.
\begin{align}
\varepsilon_{\alpha\beta} & =
- x_3~\varphi_{\alpha,\beta} \\
\varepsilon_{\alpha 3} & = \cfrac{1}{2}~\kappa\left(w^0_{,\alpha}- \varphi_\alpha\right) \\
\varepsilon_{33} & = 0
\end{align}
where is a shear correction factor.
\begin{align}
& M_{\alpha\beta,\beta}-Q_\alpha = 0 \\
& Q_{\alpha,\alpha}+q = 0
\end{align}
where
Q_\alpha := \kappa~\int_{-h}^h \sigma_{\alpha 3}~dx_3
Dynamic bending of plates
Dynamics of thin Kirchhoff plates
M_{\alpha\beta,\alpha\beta} - q(x,t) = J_1~\ddot{w}^0 - J_3~\ddot{w}^0_{,\alpha\alpha}
where, for a plate with density ,
J_1 := \int_{-h}^h \rho~dx_3 ~;~~
J_3 := \int_{-h}^h x_3^2~\rho~dx_3
and
\ddot{w}^0 = \frac{\partial^2 w^0}{\partial t^2} ~;~~
\ddot{w}^0_{,\alpha\beta} = \frac{\partial^2 \ddot{w}^0}{\partial x_\alpha\, \partial x_\beta}
See also
External links
|
|