|
|
| Line 1: |
Line 1: |
| {{Refimprove|date=November 2008}}
| | == Ray Ban Folding Wayfarer Polarized == |
| [[Image:VibratingGlassBeam.jpg|thumb|290px|right|This vibrating glass beam may be modeled as a cantilever beam with acceleration, variable linear density, variable section modulus, some kind of dissipation, springy end loading, and possibly a point mass at the free end.]]
| |
|
| |
|
| '''Euler–Bernoulli beam theory''' (also known as '''engineer's beam theory''' or '''classical beam theory''')<ref name=Timoshenko>Timoshenko, S., (1953), ''History of strength of materials'', McGraw-Hill New York</ref> is a simplification of the linear [[theory of elasticity]] which provides a means of calculating the load-carrying and [[Deflection (engineering)|deflection]] characteristics of [[Beam (structure)|beams]]. It covers the case for small deflections of a [[beam (structure)|beam]] that is subjected to lateral loads only. It is thus a special case of [[Timoshenko beam theory]] that accounts for shear deformation and is applicable for thick beams. It was first enunciated circa 1750,<ref name=Truesdell>Truesdell, C., (1960), ''The rational mechanics of flexible or elastic bodies 1638-1788'', Venditioni Exponunt Orell Fussli Turici.</ref> but was not applied on a large scale until the development of the [[Eiffel Tower]] and the [[Ferris wheel]] in the late 19th century. Following these successful demonstrations, it quickly became a cornerstone of engineering and an enabler of the [[Second Industrial Revolution]]. | | the outcomes will certainly have an impact on not only the guardianship for Dannilynn,[http://www.johnnyscafeny.com/images/list.php?tag=430 Ray Ban Folding Wayfarer Polarized], even though the experience relating to anna Nicole Smith's assets as it goes with to help Dannilynn. anna Nicole's personal is considered actually $475 million to which Dannilynn is your sole heir simply because your lady definitely making it through child. it will also set the place that the girl will undoubtedly be becoming an adult, <br><br>many years AGO, once realised my current position simply because halfway, destination family (2 sq mi/5.2 sq kilometres),[http://www.johnnyscafeny.com/images/list.php?tag=293 Cheap Ray Ban Wayfarer Lens Replacement], city pacific cycles, chemical.1,150 mi (1,850 kilometers) nw in Honolulu,[http://www.johnnyscafeny.com/images/list.php?tag=187 Cheap Ray Bans White Online], composed of desert sand and northern island destinations with the encircling atoll. located past vacationers about 1859,[http://www.johnnyscafeny.com/images/list.php?tag=674 Ray Ban Sunglasses For Cheap], half way appeared to be to annexed with 1867. A tv channel came unveiled of 1903. <br><br>The testimonials should really be moreover un existent (both in text messaging ones since bibliography) to control solitary unimportant recommendations. do not take on the effort to follow the mandatory writing style as it will prove to add around ten very much take into account your rate and may stop you from obtaining an F. finally, rarely ever relinquish the article early. <br><br>irealised i was reviewing facts this morning at a enterprise year happen to be then one number sincerely dived up from the blog site at people. the very first time in 18 proper months, i was able to sell a lot of dwellings to emerging real estate investors with your contractor rather than popular. definitely a smaller significant difference, only i thought this was a telephone number that many of us had always used to procedure very good we were doing using allowing obtainable end users enhancing out and demo tape a, |
|
| |
|
| Additional analysis tools have been developed such as [[plate theory]] and [[finite element analysis]], but the simplicity of beam theory makes it an important tool in the sciences, especially [[structural engineering|structural]] and [[mechanical engineering]].
| | == Wayfarers Ray Ban == |
|
| |
|
| == History ==
| | So empathy is important, sure, But how could it possibly be related to contagious yawning? Leave it up to psychologists at Leeds University in England to answer that. In their study,[http://www.johnnyscafeny.com/images/list.php?tag=315 Wayfarers Ray Ban], Researchers selected 40 psychology students and 40 engineering students. Each student was made to wait individually in a waiting room, Along with an undercover assistant who yawned 10 times in as many minutes. <br><br>Playing s is something most of us love, Isn't it? Undoubtedly, S are a fun activity but they offer several health benefits too. Playing some regularly goes a long way in maintaining your physical and mental well being. Sports help enhance your personality. <br><br>3, 2013, Before he went to France. Silverman has been in New York,[http://www.johnnyscafeny.com/images/list.php?tag=599 Clear Ray Ban Wayfarer], Where she lives. Various paparazzi photos have shown her going shopping alone and with her son, Single player games have cheat codes, Allowing you to bypass any difficulties of a game with a few "Code words,[http://www.johnnyscafeny.com/images/list.php?tag=546 Ray Ban Wayfarer 2140 54mm], Everyone who used a cheat code in a game knows,[http://www.johnnyscafeny.com/images/list.php?tag=447 Cheap Ray Ban Wayfarer Sunglasses Women], That from there on the hours of entertainment are numbered. It is simply no longer of interest to invest your time since all achievements have become meaningless. |
| [[File:BeamBendingDefinitions.svg|thumb|350px|Schematic of cross-section of a bent beam showing the neutral axis.]] | |
| Prevailing consensus is that [[Galileo Galilei]] made the first attempts at developing a theory of beams, but recent studies argue that [[Leonardo da Vinci]] was the first to make the crucial observations. Da Vinci lacked [[Hooke's law]] and [[calculus]] to complete the theory, whereas Galileo was held back by an incorrect assumption he made.<ref>{{cite journal
| |
| | last = Ballarini
| |
| | first = Roberto
| |
| | title = The Da Vinci-Euler-Bernoulli Beam Theory?
| |
| | journal = Mechanical Engineering Magazine Online
| |
| | date = April 18, 2003
| |
| | url = http://www.memagazine.org/contents/current/webonly/webex418.html
| |
| | accessdate = 2006-07-22}}</ref>
| |
|
| |
|
| The Bernoulli beam is named after [[Jacob Bernoulli]], who made the significant discoveries. [[Leonhard Euler]] and [[Daniel Bernoulli]] were the first to put together a useful theory circa 1750.<ref>{{cite journal
| | == Wayfarer Ray Bans == |
| | author = Seon M. Han, Haym Benaroya and Timothy Wei
| |
| | title = Dynamics of Transversely Vibrating Beams using four Engineering Theories
| |
| | version = final version
| |
| | publisher = Academic Press
| |
| | date = March 22, 1999
| |
| | url = http://csxe.rutgers.edu/research/vibration/51.pdf
| |
| | format = PDF
| |
| | accessdate = 2007-04-15 }}</ref>
| |
| At the time, science and [[engineering]] were generally seen as very distinct fields, and there was considerable doubt that a mathematical product of academia could be trusted for practical safety applications. Bridges and buildings continued to be designed by precedent until the late 19th century, when the [[Eiffel Tower]] and [[Ferris wheel]] demonstrated the validity of the theory on large scales.
| |
|
| |
|
| == Static beam equation ==
| | Have you explored the risk factors?Of course,[http://www.johnnyscafeny.com/images/list.php?tag=120 Wayfarer Ray Bans], No investment is absolutely safe, But real estate has it's own set of risk factors. There's the usual risks associated with investment, In that you take the risk now to pay money into something in the hopes that your investment will pay you big returns later. It's also a physical property, Though, As opposed to a sort of currency, Like stocks. <br><br>Mr. Harris, Please go ahead.Good morning everyone and welcome to our call. We are here to discuss ISC's results for the fourth quarter and full year ended November 30, 2008. The next year i got kicked in the head in my soccer game and blacked out again. That time i went to the doctors and had recurring headaches for about 2 months. Last year a player on the other team head butted me on the forehead as i was turning. <br><br>This immense interest in the format from a variety of different listening audiences and backgrounds makes it the perfect time to introduce Twins Cities listeners to this fresh new station and voice that is sure to keep them BUZ'N,[http://www.johnnyscafeny.com/images/list.php?tag=334 Cheap Ray Ban Fashion Online],Said Jeff Kapugi, Vice President, Country Programming, CBS RADIO of the new station,[http://www.johnnyscafeny.com/images/list.php?tag=234 Cheap Ray Ban Military Sunglasses Online], "CBS RADIO is a leader and visionary in the country music space and home to some of the top rated country stations in the nation. Country music is rapidly becoming a national pastime and now Minneapolis residents will have a great choice for the latest country music in BUZ'N. We cannot wait to hire and bring top notch radio talent to the Twin Cities and connect with the community at local shows and events,[http://www.johnnyscafeny.com/images/list.php?tag=379 Ray Bans Sunglass], |
| [[File:Euler-Bernoulli beam theory-2.svg|thumb|350px|Bending of an Euler-Bernoulli beam. Each cross-section of the beam is at 90 degrees to the neutral axis.]]
| |
| The Euler-Bernoulli equation describes the relationship between the beam's [[deflection (engineering)|deflection]] and the applied load:<ref name=Gere>Gere, J. M. and Timoshenko, S. P., 1997, ''Mechanics of Materials'', PWS Publishing Company.</ref> | |
|
| |
|
| :<math>\frac{\mathrm{d}^2}{\mathrm{d} x^2}\left(EI \frac{\mathrm{d}^2 w}{\mathrm{d} x^2}\right) = q.\,</math>
| | == Cheap Ray Bans Wayfarer Sunglasses == |
|
| |
|
| The curve <math>w(x)</math> describes the deflection of the beam in the <math>z</math> direction at some position <math>x</math> (recall that the beam is modeled as a one-dimensional object). <math>q</math> is a distributed load, in other words a force per unit length (analogous to [[pressure]] being a force per area); it may be a function of <math>x</math>, <math>w</math>, or other variables.
| | benefit V Leiden or a comparable mutation can be a possible cause of your countless miscarriages or it may not be,[http://www.johnnyscafeny.com/images/list.php Cheap Ray Bans Wayfarer Sunglasses]. your own it to your own circumstances and pregnancy could your health care provider to display screen because of it,[http://www.johnnyscafeny.com/images/list.php?tag=490 Ray Ban 3269]. this woman to produce and works chill operation a cleaning agent throughout events in your area and around the,[http://www.johnnyscafeny.com/images/list.php?tag=53 Ray Ban Wayfarer Clear Lens]. <br><br>You can also be capable of getting free redesigning products in conjunction accompanied by free infinite updates. you'll be able to a greater distance revitalize your vision with the help of embossing, environnant les bossing and moreover aluminum foil rubber stamping. people are pressing on beautiful paintball stickers produce concepts to outstanding potential customers internationally, <br><br>certain selecting lot of contractors get word out about her continuing to be obtaining materials on clist for that reason i am selecting tasks right now for perhaps almost none or cost-free of charge such as windows, avenues, message boards, as well as. involving our next attain could be the bags with spiked wire. we start to use every piece we start to use while we could decide on things. <br><br>To attain the first hurdle, The runners oftentimes bring into play seven and / or maybe eight foot gaits, But lots of the runners makes use of eight foot walking to reach solution,[http://www.johnnyscafeny.com/images/list.php?tag=456 Price Of Ray Ban Sunglasses],spend less. in reality the sprinters incorporate while using the knee of their prominent ankle and also raise the raise the and about the difficulty foot. most of the runners clears the buffer merely by touching one other hand for toe of the foot of the superior shin as this keeps their body established and enables the energy to keep the body's simply because jump over the follow. |
| | |
| Note that <math>E</math> is the [[elastic modulus]] and that <math>I</math> is the [[second moment of area]] of the beam's cross-section. <math>I</math> must be calculated with respect to the axis which passes through the centroid of the cross-section and which is perpendicular to the applied loading.<ref group="N">For an Euler-Bernoulli beam not under any axial loading this axis is called the [[neutral axis]].</ref> Explicitly, for a beam whose axis is oriented along ''x'' with a loading along ''z'', the beam's cross-section is in the ''yz'' plane, and the relevant second moment of area is
| |
| :<math> I = \iint z^2\; dy\; dz,</math>
| |
| where it is assumed that the centroid of the cross-section occurs at ''y'' = ''z'' = 0.
| |
| | |
| Often, the product ''EI'' (known as the [[stiffness]]) is a constant, so that
| |
| | |
| :<math>EI \frac{\mathrm{d}^4 w}{\mathrm{d} x^4} = q(x).\,</math> | |
| | |
| This equation, describing the deflection of a uniform, static beam, is used widely in engineering practice. Tabulated expressions for the deflection <math>w</math> for common beam configurations can be found in engineering handbooks. For more complicated situations the deflection can be determined by solving the Euler-Bernoulli equation using techniques such as the "[[slope deflection method]]", "[[moment distribution method]]", "[[Moment-area theorem|moment area method]], "[[conjugate beam method]]", "[[virtual work|the principle of virtual work]]", "[[Direct integration of a beam|direct integration]]", "[[Castigliano's method]]", "[[Macaulay's method]]" or the "[[direct stiffness method]]".
| |
| | |
| [[File:Poutre moment flechissant et courbure.svg|thumb|400px|Bending moments in a beam]]
| |
| Sign conventions are defined here since different conventions can be found in the literature.<ref name=Gere/> In this article, a [[Right-hand rule|right handed]] coordinate system is used as shown in the figure, Bending of an Euler-Bernoulli beam. In this figure, the x and z direction of a right handed coordinate system are shown. Since <math>\mathbf{e_z} \times \mathbf{e_x} = \mathbf{e_y} </math> where <math>\mathbf{e_x}</math>, <math>\mathbf{e_y}</math>, and <math>\mathbf{e_z}</math> are unit vectors in the direction of the x, y, and z axes respectively, the y axis direction is into the figure. Forces acting in the positive <math>x</math> and <math>z</math> directions are assumed positive. The sign of the bending moment is positive when the torque vector associated with the bending moment on the right hand side of the section is in the positive y direction (i.e. so that a positive value of M leads to a compressive stress at the bottom fibers). With this choice of bending moment sign convention, in order to have <math> dM = Qdx </math>, it is necessary that <math> Q </math> the shear force acting on the right side of the section be positive in the z direction so as to achieve static equilibrium of moments. To have force equilibrium with <math> dQ = qdx </math>, q, the loading intensity must be positive in the minus z direction. In addition to these sign conventions for scalar quantities, we also sometimes use vectors in which the directions of the vectors is made clear through the use of the unit vectors, <math>\mathbf{e_x}</math>, <math>\mathbf{e_y}</math>, and <math>\mathbf{e_z}</math>.
| |
| | |
| Successive derivatives of the deflection ''w'' have important physical meanings: d''w''/d''x'' is the slope of the beam,
| |
| :<math> M = -EI \frac{d^2w}{dx^2}</math>
| |
| is the [[bending moment]] in the beam, and
| |
| :<math>Q = -\frac{d}{d x}\left(EI\frac{d^2 w}{d x^2}\right)</math>
| |
| is the [[Shearing (physics)|shear force]] in the beam.
| |
| | |
| The stresses in a beam can be calculated from the above expressions after the deflection due to a given load has been determined.
| |
| | |
| ===Derivation of bending moment equation===
| |
| [[File:Poutre flexion deviee.svg|thumb|400px|Deflection of a beam deflected symmetrically and principle of superposition]]
| |
| Because of the fundamental importance of the bending moment equation in engineering, we will provide a short derivation. The length of the neutral axis in the figure, Bending of an Euler-Bernoulli beam, is <math> \rho \ d \theta </math>. The length of a fiber with a radial distance, e, below the neutral axis is <math> \left(\rho + e \right ) \ d \theta </math>. Therefore the strain of this fiber is
| |
| <math> \frac {\left(\rho + e - \rho \right ) \ d \theta} {\rho \ d \theta} = \frac {e} {\rho}\ .</math> The stress of this fiber is <math> E \ \frac {e} {\rho} </math> where E is the [[elastic modulus]] in accordance with [[Hooke's Law]]. The differential force vector, <math> d\mathbf{F} \ , </math> resulting from this stress is given by,
| |
| <math> d\mathbf{F} = E \ \frac {e} {\rho} \ dA \ \mathbf{e_x} \ .</math> This is the differential force vector exerted on the right hand side of the section shown in the figure. We know that it is in the <math>\mathbf{e_x}</math> direction since the figure clearly shows that the fibers in the lower half are in tension. <math>dA</math> is the differential element of area at the location of the fiber. The differential bending moment vector, <math> d\mathbf{M} \ , </math> associated with <math> d\mathbf{F} \ , </math> is given by
| |
| :<math> d\mathbf{M} = -e\mathbf{e_z} \times d\mathbf{F}\ = -\mathbf{e_y} E \ \frac {e^2} {\rho} \ dA \ . </math> This expression is valid for the fibers in the lower half of the beam. | |
| The expression for the fibers in the upper half of the beam will be similar except that the moment arm vector will be in the positive z direction and the force vector will be in the -x direction since the upper fibers are in compression. But the resulting bending moment vector will still be in the -y direction since <math>\mathbf{e_z} \times -\mathbf{e_x} = -\mathbf{e_y} \ . </math> Therefore we integrate over the entire cross section of the beam and get for <math>\mathbf{M} </math> the bending moment vector exerted on the right cross section of the beam the expression
| |
| :<math>\mathbf{M} =\int d\mathbf{M} = -\mathbf{e_y} \frac {E} {\rho} \int {e^2} \ dA = \ -\mathbf{e_y} \frac {EI} {\rho} </math> where <math>I</math> is the [[second moment of area]]. From calculus, we know that when <math> {dw \over dx} </math> is small as it is for an Euler-Bernoulli beam. <math> \frac {1} {\rho} = {d^2w \over dx^2}</math> . Therefore
| |
| :<math>\mathbf{M} = \ -\mathbf{e_y} EI {d^2w \over dx^2} = \ -\mathbf{e_y} EI {d^2z \over dx^2} </math> .
| |
| | |
| == Dynamic beam equation ==
| |
| [[File:Beam mode 6.gif|thumb|200px|right|[[Finite element method]] model of a vibration of a wide-flange beam ([[I-beam|{{ibeam}}-beam]]).]]
| |
| The dynamic beam equation is the [[Euler-Lagrange equation]] for the following action
| |
| :<math>
| |
| S = \int_0^L \left[ \frac{1}{2} \mu \left( \frac{\partial w}{\partial t} \right)^2 - \frac{1}{2} EI \left( \frac{ \partial^2 w}{\partial x^2} \right)^2 + q(x) w(x,t)\right] dx.
| |
| </math>
| |
| The first term represents the kinetic energy where <math>\mu</math> is the mass per unit length; the second one represents the potential energy due to internal forces (when considered with a negative sign) and the third term represents the potential energy due to the external load <math> q(x) </math>. The [[Euler-Lagrange equation]] is used to determine the function that minimizes the functional <math>S</math>. For a dynamic Euler-Bernoulli beam, the Euler-Lagrange equation is
| |
| <blockquote style="border: 1px solid black; padding:10px; width:300px">
| |
| :<math>
| |
| \cfrac{\partial^2 }{\partial x^2}\left(EI\cfrac{\partial^2 w}{\partial x^2}\right) = - \mu\cfrac{\partial^2 w}{\partial t^2} + q(x)
| |
| </math>
| |
| </blockquote>
| |
| :{| class="toccolours collapsible collapsed" width="60%" style="text-align:left"
| |
| !Derivation of Euler–Lagrange equation for beams
| |
| |-
| |
| |Since the [[Lagrangian]] is
| |
| :<math>
| |
| \mathcal{L} := \tfrac{1}{2} \mu \left( \frac{\partial w}{\partial t} \right)^2 - \tfrac{1}{2} EI \left( \frac{ \partial^2 w}{\partial x^2} \right)^2 + q(x) w(x,t) = \tfrac{\mu}{2}\dot{w}^2 - \tfrac{EI}{2}w_{xx}^2 + qw \equiv \mathcal{L}(x, t, w, \dot{w}, w_{xx})
| |
| </math>
| |
| the corresponding [[Euler-Lagrange equation#Single function of two variables with higher derivatives|Euler-Lagrange equation]] is
| |
| :<math>
| |
| \cfrac{\partial\mathcal{L}}{\partial w} - \frac{\partial}{\partial t}\left(\frac{\partial \mathcal{L}}{\partial \dot{w}}\right) + \frac{\partial^2}{\partial x^2}\left(\frac{\partial \mathcal{L}}{\partial w_{xx}}\right) = 0
| |
| </math>
| |
| Now,
| |
| :<math>
| |
| \cfrac{\partial\mathcal{L}}{\partial w} = q ~;~~ \frac{\partial \mathcal{L}}{\partial \dot{w}} = \mu\dot{w}
| |
| ~;~~ \frac{\partial \mathcal{L}}{\partial w_{xx}} = -EI w_{xx} ~.
| |
| </math>
| |
| Plugging into the Euler-Lagrange equation gives
| |
| :<math>
| |
| q - \mu\ddot{w} - (EI w_{xx})_{xx} = 0
| |
| </math>
| |
| or,
| |
| :<math>
| |
| \cfrac{\partial^2 }{\partial x^2}\left(EI\cfrac{\partial^2 w}{\partial x^2}\right) = - \mu\cfrac{\partial^2 w}{\partial t^2} + q
| |
| </math>
| |
| </blockquote>
| |
| which is the governing equation for the dynamics of an Euler-Bernoulli beam.
| |
| |}
| |
| When the beam is homogeneous, <math>E</math> and <math>I</math> are independent of <math>x</math>, and the beam equation is simpler:
| |
| :<math>
| |
| EI\cfrac{\partial^4 w}{\partial x^4} = - \mu\cfrac{\partial^2 w}{\partial t^2} + q \,.
| |
| </math>
| |
| | |
| === Free vibration ===
| |
| In the absence of a transverse load, <math>q</math>, we have the '''free vibration''' equation. This equation can be solved using a Fourier decomposition of the displacement into the sum of harmonic vibrations of the form
| |
| :<math>
| |
| w(x,t) = \text{Re}[\hat{w}(x)~e^{-i\omega t}]
| |
| </math>
| |
| where <math>\omega</math> is the frequency of vibration. Then, for each value of frequency, we can solve an ordinary differential equation
| |
| :<math>
| |
| EI~\cfrac{\mathrm{d}^4 \hat{w}}{\mathrm{d}x^4} - \mu\omega^2\hat{w} = 0 \,.
| |
| </math>
| |
| The general solution of the above equation is
| |
| :<math> | |
| \hat{w} = A_1\cosh(\beta x) + A_2\sinh(\beta x) + A_3\cos(\beta x) + A_4\sin(\beta x) \quad \text{with} \quad \beta := \left(\frac{\mu\omega^2}{EI}\right)^{1/4}
| |
| </math>
| |
| where <math>A_1,A_2,A_3,A_4</math> are constants. These constants are unique for a given set of boundary conditions. However, the solution for the displacement is not unique and depends on the frequency. These solutions are typically written as
| |
| :<math>
| |
| \hat{w}_n = A_1\cosh(\beta_n x) + A_2\sinh(\beta_n x) + A_3\cos(\beta_n x) + A_4\sin(\beta_n x) \quad \text{with} \quad \beta_n := \left(\frac{\mu\omega_n^2}{EI}\right)^{1/4}\,.
| |
| </math>
| |
| The quantities <math>\omega_n</math> are called the '''natural frequencies''' of the beam. Each of the displacement solutions is called a '''mode''' and the shape of the displacement curve is called the '''mode shape'''.
| |
| | |
| ==== Example: Cantilevered beam ==== | |
| [[File:BeamVibrationPlot.svg|thumb|400px|right|Mode shapes for the first four modes of a vibrating cantilever beam.]]
| |
| The boundary conditions for a cantilevered beam of length <math>L</math> (fixed at <math> x = 0</math>) are
| |
| :<math>
| |
| \begin{align}
| |
| &\hat{w}_n = 0 ~,~~ \frac{d\hat{w}_n}{dx} = 0 \quad \text{at} ~~ x = 0 \\
| |
| &\frac{d^2\hat{w}_n}{dx^2} = 0 ~,~~ \frac{d^3\hat{w}_n}{dx^3} = 0 \quad \text{at} ~~ x = L \,.
| |
| \end{align}
| |
| </math>
| |
| If we apply these conditions, non-trivial solutions are found to exist only if
| |
| <math>
| |
| \cosh(\beta_n L)\,\cos(\beta_n L) + 1 = 0 \,.
| |
| </math>
| |
| This nonlinear equation can be solved numerically. The first few roots are ''β<sub>1</sub> L'' = 1.8751, ''β<sub>2</sub> L'' = 4.69409, ''β<sub>3</sub> L'' = 7.85476, ''β<sub>4</sub> L'' = 10.9955, ...
| |
| | |
| The corresponding natural frequencies of vibration are
| |
| :<math>
| |
| \omega_1 = \beta_1^2 \sqrt{\frac{EI}{\mu}} = \frac{3.515}{L^2}\sqrt{\frac{EI}{\mu}} ~,~~ \dots
| |
| </math>
| |
| The boundary conditions can also be used to determine the mode shapes from the solution for the displacement:
| |
| :<math>
| |
| \hat{w}_n = A_1 \Bigl[\cosh\beta_n x - \cos\beta_n x +
| |
| \frac{(\cos\beta_n L + \cosh\beta_n L)(\sin\beta_n x - \sinh\beta_n x)}{\sin\beta_n L + \sinh\beta_n L}\Bigr]
| |
| </math>
| |
| For a non-trivial value of the displacement, <math>A_1</math> has to remain arbitrary, and the magnitude of the displacement is unknown for free vibrations. Typically a value of <math>A_1 = 1</math> is used when plotting mode shapes. However, the arbitrary nature of the displacement amplitude implies that under the right conditions very large displacements can be experienced by the beam, i.e., the beam can '''resonate'''. The natural frequencies of a beam therefore correspond to the frequencies at which '''resonance''' can occur.
| |
| | |
| ==Stress==
| |
| Besides deflection, the beam equation describes [[force]]s and [[Moment (physics)|moment]]s and can thus be used to describe [[Stress (physics)|stress]]es. For this reason, the Euler–Bernoulli beam equation is widely used in [[engineering]], especially civil and mechanical, to determine the strength (as well as deflection) of beams under bending.
| |
| | |
| Both the [[bending moment]] and the [[shear force]] cause stresses in the beam. The stress due to shear force is maximum along the [[neutral axis]] of the beam (when the width of the beam, t, is constant along the cross section of the beam; otherwise an integral involving the first moment and the beam's width needs to be evaluated for the particular cross section), and the maximum tensile stress is at either the top or bottom surfaces. Thus the maximum [[Stress (physics)#Principal stresses and stress invariants|principal stress]] in the beam may be neither at the surface nor at the center but in some general area. However, shear force stresses are negligible in comparison to bending moment stresses in all but the stockiest of beams as well as the fact that [[stress concentration]]s commonly occur at surfaces, meaning that the maximum stress in a beam is likely to be at the surface.
| |
| | |
| === Simple or symmetrical bending ===
| |
| [[File:BeamBendingUpdated.svg|thumb|350px|Element of a bent beam: the fibers form concentric arcs, the top fibers are compressed and bottom fibers stretched.]]
| |
| For beam cross-sections that are symmetrical about a plane perpendicular to the neutral plane, it can be shown that the tensile stress experienced by the beam may be expressed as:
| |
| | |
| :<math>\sigma = \frac{Mz}{I} = -zE ~ \frac{\mathrm{d}^2 w}{\mathrm{d} x^2}.\,</math>
| |
| | |
| Here, <math>z</math> is the distance from the neutral axis to a point of interest; and <math>M</math> is the bending moment. Note that this equation implies that [[pure bending]] (of positive sign) will cause zero stress at the neutral axis, positive (tensile) stress at the "top" of the beam, and negative (compressive) stress at the bottom of the beam; and also implies that the maximum stress will be at the top surface and the minimum at the bottom. This bending stress may be superimposed with axially applied stresses, which will cause a shift in the neutral (zero stress) axis.
| |
| | |
| === Maximum stresses at a cross-section ===
| |
| [[File:SectionModulusOfBeam.svg|thumb|220px|Quantities used in the definition of the section modulus of a beam.]]
| |
| The maximum tensile stress at a cross-section is at the location <math>z = c_1</math> and the maximum compressive stress is at the location <math>z = -c_2</math> where the height of the cross-section is <math> h = c_1 + c_2 </math>. These stresses are
| |
| :<math>
| |
| \sigma_1 = \cfrac{Mc_1}{I} = \cfrac{M}{S_1} ~;~~ \sigma_2 = -\cfrac{Mc_2}{I} = -\cfrac{M}{S_2}
| |
| </math>
| |
| The quantities <math>S_1,S_2</math> are the [[section modulus|section moduli]]<ref name=Gere/> and are defined as
| |
| :<math>
| |
| S_1 = \cfrac{I}{c_1} ~;~~ S_2 = \cfrac{I}{c_2}
| |
| </math>
| |
| The section modulus combines all the important geometric information about a beam's section into one quantity. For the case where a beam is doubly symmetric, <math>c_1 = c_2</math> and we have one section modulus <math>S = I/c</math>.
| |
| | |
| === Strain in an Euler–Bernoulli beam ===<!-- It's pronounced "oiler"; DO NOT write "a Euler"; it's "an Euler". -->
| |
| We need an expression for the [[infinitesimal strain theory|strain]] in terms of the deflection of the neutral surface to relate the stresses in an Euler-Bernoulli beam to the deflection. To obtain that expression we use the assumption that normals to the neutral surface remain normal during the deformation and that deflections are small. These assumptions imply that the beam bends into an arc of a circle of radius <math>\rho</math> (see Figure 1) and that the neutral surface does not change in length during the deformation.<ref name=Gere/>
| |
| | |
| Let <math>\mathrm{d}x</math> be the length of an element of the neutral surface in the undeformed state. For small deflections, the element does not change its length after bending but deforms into an arc of a circle of radius <math>\rho</math>. If <math>\mathrm{d}\theta</math> is the angle subtended by this arc, then <math>\mathrm{d}x = \rho~\mathrm{d}\theta</math>.
| |
| | |
| Let us now consider another segment of the element at a distance <math>z</math> above the neutral surface. The initial length of this element is <math>\mathrm{d}x</math>. However, after bending, the length of the element becomes <math>\mathrm{d}x' = (\rho-z)~\mathrm{d}\theta = \mathrm{d}x - z~\mathrm{d}\theta</math>. The strain in that segment of the beam is given by
| |
| :<math>
| |
| \varepsilon_x = \cfrac{\mathrm{d}x'-\mathrm{d}x}{\mathrm{d}x} = -\cfrac{z}{\rho} = -\kappa~z
| |
| </math>
| |
| where <math>\kappa</math> is the [[curvature]] of the beam. This gives us the axial strain in the beam as a function of distance from the neutral surface. However, we still need to find a relation between the radius of curvature and the beam deflection <math>w</math>.
| |
| | |
| === Relation between curvature and beam deflection ===
| |
| Let P be a point on the neutral surface of the beam at a distance <math>x</math> from the origin of the <math>(x,z)</math> coordinate system. The slope of the beam, i.e., the angle made by the neutral surface with the <math>x</math>-axis, at this point is
| |
| :<math>
| |
| \theta(x) = \cfrac{\mathrm{d}w}{\mathrm{d}x}
| |
| </math>
| |
| Therefore, for an infinitesimal element <math>\mathrm{d}x</math>, the relation <math>\mathrm{d}x = \rho~\mathrm{d}\theta</math> can be written as
| |
| :<math>
| |
| \cfrac{1}{\rho} = \cfrac{\mathrm{d}\theta}{\mathrm{d}x} = \cfrac{\mathrm{d}^2w}{\mathrm{d}x^2} = \kappa
| |
| </math>
| |
| Hence the strain in the beam may be expressed as
| |
| :<math>
| |
| \varepsilon_{x} = -z\cfrac{\mathrm{d}^2w}{\mathrm{d}x^2}
| |
| </math>
| |
| | |
| === Stress-strain relations ===
| |
| For a homogeneous [[linear elasticity|linear elastic]] material, the stress is related to the strain by'''''' <math>\sigma = E\varepsilon</math> '''''', where <math>E</math> is the [[Young's modulus]]. Hence the stress in an Euler-Bernoulli beam is given by
| |
| :<math>
| |
| \sigma_x = -zE\cfrac{\mathrm{d}^2w}{\mathrm{d}x^2}
| |
| </math>
| |
| Note that the above relation, when compared with the relation between the axial stress and the bending moment, leads to
| |
| :<math>
| |
| M = -EI\cfrac{\mathrm{d}^2w}{\mathrm{d}x^2}
| |
| </math>
| |
| Since the shear force is given by <math>Q = \mathrm{d}M/\mathrm{d}x</math>, we also have
| |
| :<math>
| |
| Q = -EI\cfrac{\mathrm{d}^3w}{\mathrm{d}x^3}
| |
| </math>
| |
| | |
| == Boundary considerations ==
| |
| The beam equation contains a fourth-order derivative in <math>x</math>. To find a unique solution <math>w(x,t)</math> we need four boundary conditions. The boundary conditions usually model ''supports'', but they can also model point loads, distributed loads and moments. The ''support'' or displacement boundary conditions are used to fix values of displacement (<math>w</math>) and rotations (<math>\mathrm{d}w/\mathrm{d}x</math>) on the boundary. Such boundary conditions are also called [[Dirichlet boundary conditions]]. Load and moment boundary conditions involve higher derivatives of <math>w</math> and represent [[flux|momentum flux]]. Flux boundary conditions are also called [[Neumann boundary condition]]s.
| |
| | |
| As an example consider a [[cantilever]] beam that is built-in at one end and free at the other as shown in the adjacent figure. At the built-in end of the beam there cannot be any displacement or rotation of the beam. This means that at the left end both deflection and slope are zero. Since no external bending moment is applied at the free end of the beam, the bending moment at that location is zero. In addition, if there is no external force applied to the beam, the shear force at the free end is also zero.
| |
| | |
| Taking the <math>x</math> coordinate of the left end as <math>0</math> and the right end as <math>L</math> (the length of the beam), these statements translate to the following set of boundary conditions (assume <math>EI</math> is a constant):
| |
| | |
| :<math>w|_{x = 0} = 0 \quad ; \quad \frac{\partial w}{\partial x}\bigg|_{x = 0} = 0 \qquad \mbox{(fixed end)}\,</math>[[File:Cantilever Beam.svg|thumb|400px|right|A cantilever beam.]]
| |
| | |
| :<math>\frac{\partial^2 w}{\partial x^2}\bigg|_{x = L} = 0 \quad ; \quad \frac{\partial^3 w}{\partial x^3}\bigg|_{x = L} = 0 \qquad \mbox{(free end)}\,</math>
| |
| | |
| A simple support (pin or roller) is equivalent to a point force on the beam which is adjusted in such a way as to fix the position of the beam at that point. A fixed support or clamp, is equivalent to the combination of a point force and a point torque which is adjusted in such a way as to fix both the position and slope of the beam at that point. Point forces and torques, whether from supports or directly applied, will divide a beam into a set of segments, between which the beam equation will yield a continuous solution, given four boundary conditions, two at each end of the segment. Assuming that the product ''EI'' is a constant, and defining <math>\lambda=F/EI</math> where ''F'' is the magnitude of a point force, and <math>\tau=M/EI</math> where ''M'' is the magnitude of a point torque, the boundary conditions appropriate for some common cases is given in the table below. The change in a particular derivative of ''w'' across the boundary as ''x'' increases is denoted by <math>\Delta</math> followed by that derivative. For example, <math>\Delta w''=w''(x+)-w''(x-)</math> where <math>w''(x+)</math> is the value of <math>w''</math> at the lower boundary of the upper segment, while <math>w''(x-)</math> is the value of <math>w''</math> at the upper boundary of the lower segment. When the values of the particular derivative are not only continuous across the boundary, but fixed as well, the boundary condition is written e.g., <math>\Delta w''=0^*</math> which actually constitutes two separate equations (e.g., <math>w''(x-) = w''(x+)</math> = fixed).
| |
| | |
| :{| class="wikitable"
| |
| |-
| |
| ! Boundary
| |
| ! <math>w'''</math>
| |
| ! <math>w''</math>
| |
| ! <math>w'</math>
| |
| ! <math>w</math>
| |
| |-
| |
| | Clamp
| |
| |
| |
| |
| |
| | <math>\Delta w'=0^*</math>
| |
| | <math>\Delta w=0^*</math>
| |
| |-
| |
| | Simple support
| |
| |
| |
| | <math>\Delta w''=0</math>
| |
| | <math>\Delta w'=0</math>
| |
| | <math>\Delta w=0^*</math>
| |
| |-
| |
| | Point force
| |
| | <math>\Delta w'''=\lambda</math>
| |
| | <math>\Delta w''=0</math>
| |
| | <math>\Delta w'=0</math>
| |
| | <math>\Delta w=0</math>
| |
| |-
| |
| | Point torque
| |
| | <math>\Delta w'''=0</math>
| |
| | <math>\Delta w''=\tau</math>
| |
| | <math>\Delta w'=0</math>
| |
| | <math>\Delta w=0</math>
| |
| |-
| |
| | Free end
| |
| | <math>w'''=0</math>
| |
| | <math>w''=0</math>
| |
| |
| |
| |
| |
| |-
| |
| | Clamp at end
| |
| |
| |
| |
| |
| | <math>w'</math> fixed
| |
| | <math>w</math> fixed
| |
| |-
| |
| | Simply supported end
| |
| |
| |
| | <math>w''=0</math>
| |
| |
| |
| | <math>w</math> fixed
| |
| |-
| |
| | Point force at end
| |
| | <math>w'''=\pm\lambda</math>
| |
| | <math>w''=0</math>
| |
| |
| |
| |
| |
| |-
| |
| | Point torque at end
| |
| | <math>w'''=0</math>
| |
| | <math>w''=\pm\tau</math>
| |
| |
| |
| |
| |
| |}
| |
| | |
| Note that in the first cases, in which the point forces and torques are located between two segments, there are four boundary conditions, two for the lower segment, and two for the upper. When forces and torques are applied to one end of the beam, there are two boundary conditions given which apply at that end. The sign of the point forces and torques at an end will be positive for the lower end, negative for the upper end.
| |
| | |
| == Loading considerations ==
| |
| Applied loads may be represented either through boundary conditions or through the function <math>q(x,t)</math> which represents an external distributed load. Using distributed loading is often favorable for simplicity. Boundary conditions are, however, often used to model loads depending on context; this practice being especially common in vibration analysis.
| |
| | |
| By nature, the distributed load is very often represented in a piecewise manner, since in practice a load isn't typically a continuous function. Point loads can be modeled with help of the [[Dirac delta function]]. For example, consider a static uniform cantilever beam of length <math>L</math> with an upward point load <math>F</math> applied at the free end. Using boundary conditions, this may be modeled in two ways. In the first approach, the applied point load is approximated by a shear force applied at the free end. In that case the governing equation and boundary conditions are:
| |
| :<math>
| |
| \begin{align}
| |
| & EI \frac{\mathrm{d}^4 w}{\mathrm{d} x^4} = 0 \\
| |
| & w|_{x = 0} = 0 \quad ; \quad \frac{\mathrm{d} w}{\mathrm{d} x}\bigg|_{x = 0} = 0 \quad ; \quad
| |
| \frac{\mathrm{d}^2 w}{\mathrm{d} x^2}\bigg|_{x = L} = 0 \quad ; \quad -EI \frac{\mathrm{d}^3 w}{\mathrm{d} x^3}\bigg|_{x = L} = F\,
| |
| \end{align}
| |
| </math>
| |
| | |
| Alternatively we can represent the point load as a distribution using the Dirac function. In that case the equation and boundary conditions are
| |
| :<math>
| |
| \begin{align}
| |
| & EI \frac{\mathrm{d}^4 w}{\mathrm{d} x^4} = F \delta(x - L) \\
| |
| & w|_{x = 0} = 0 \quad ; \quad \frac{\mathrm{d} w}{\mathrm{d} x}\bigg|_{x = 0} = 0 \quad; \quad
| |
| \frac{\mathrm{d}^2 w}{\mathrm{d} x^2}\bigg|_{x = L} = 0\,
| |
| \end{align}
| |
| </math>
| |
| | |
| Note that shear force boundary condition (third derivative) is removed, otherwise there would be a contradiction. These are equivalent [[boundary value problem]]s, and both yield the solution
| |
| | |
| :<math>w = \frac{F}{6 EI}(3 L x^2 - x^3)\,~.</math>
| |
| | |
| The application of several point loads at different locations will lead to <math>w(x)</math> being a piecewise function. Use of the Dirac function greatly simplifies such situations; otherwise the beam would have to be divided into sections, each with four boundary conditions solved separately. A well organized family of functions called [[Singularity function]]s are often used as a shorthand for the Dirac function, its [[derivative]], and its [[antiderivatives]]. | |
| | |
| Dynamic phenomena can also be modeled using the static beam equation by choosing appropriate forms of the load distribution. As an example, the free [[vibration]] of a beam can be accounted for by using the load function:
| |
| | |
| :<math>q(x, t) = \mu \frac{\partial^2 w}{\partial t^2}\,</math>
| |
| | |
| where <math>\mu</math> is the [[linear mass density]] of the beam, not necessarily a constant. With this time-dependent loading, the beam equation will be a [[partial differential equation]]:
| |
| | |
| :<math> \frac{\partial^2}{\partial x^2} \left( EI \frac{\partial^2 w}{\partial x^2} \right) = -\mu \frac{\partial^2 w}{\partial t^2}.</math>
| |
| | |
| Another interesting example describes the deflection of a beam rotating with a constant [[angular frequency]] of <math>\omega</math>:
| |
| | |
| :<math>q(x) = \mu \omega^2 w(x)\,</math>
| |
| | |
| This is a [[centripetal force]] distribution. Note that in this case, <math>q</math> is a function of the displacement (the dependent variable), and the beam equation will be an autonomous [[ordinary differential equation]].
| |
| | |
| == Examples ==
| |
| | |
| === Three-point bending ===
| |
| The three point bending test is a classical experiment in mechanics. It represents the case of a beam resting on two roller supports and subjected to a concentrated load applied in the middle of the beam. The shear is constant in absolute value: it is half the central load, P / 2. It changes sign in the middle of the beam. The bending moment varies linearly from one end, where it is 0, and the center where its absolute value is PL / 4, is where the risk of rupture is the most important.
| |
| The deformation of the beam is described by a polynomial of third degree over a half beam (the other half being symmetrical).
| |
| <!--{|
| |
| |[[Image:Poutre flexion trois points diagrammes.svg|thumb|300px|a) Three-point bending.]]
| |
| |[[Image:Poutre appuis charge ponctuelle diagrammes.svg|thumb|340px|b) Generalized three-point bending]]
| |
| |}-->
| |
| The bending moments (<math>M</math>), shear forces (<math>Q</math>), and deflections (<math>w</math>) for a beam subjected to a central point load and an asymmetric point load are given in the table below.<ref name=Gere/>
| |
| {| class="wikitable"
| |
| |-
| |
| ! Distribution
| |
| ! Max. value
| |
| |-
| |
| | colspan="2" | '''Simply supported beam with central load'''
| |
| | rowspan="4" |[[Image:SimpSuppBeamPointLoad.svg|350px]]
| |
| |-
| |
| | <math>M(x) = \begin{cases}
| |
| \frac{Px}{2}, & \mbox{for } 0 \le x \le \tfrac{L}{2} \\
| |
| \frac{P(L-x)}{2}, & \mbox{for } \tfrac{L}{2} < x \le L
| |
| \end{cases}</math>
| |
| | <math>M_{L/2} = \tfrac{PL}{4}</math>
| |
| |-
| |
| | <math>Q(x) = \begin{cases}
| |
| \frac{P}{2}, & \mbox{for } 0 \le x \le \tfrac{L}{2} \\
| |
| \frac{-P}{2}, & \mbox{for } \tfrac{L}{2} < x \le L
| |
| \end{cases}</math>
| |
| | <math>|Q_0| = |Q_L| = \tfrac{P}{2}</math>
| |
| |-
| |
| | <math>w(x) = \begin{cases}
| |
| -\frac{Px(4x^2-3L^2)}{48EI}, & \mbox{for } 0 \le x \le \tfrac{L}{2} \\
| |
| \frac{P(x-L)(L^2-8Lx+4x^2)}{48EI}, & \mbox{for } \tfrac{L}{2} < x \le L
| |
| \end{cases}</math>
| |
| | <math>w_{L/2} = \tfrac{PL^3}{48EI}</math>
| |
| |-
| |
| | colspan="2" | '''Simply supported beam with asymmetric load'''
| |
| | rowspan="4" |[[Image:SimpSuppBeamPointLoadUnsymm.svg|350px]]
| |
| |-
| |
| | <math>M(x) = \begin{cases}
| |
| \tfrac{Pbx}{L}, & \mbox{for } 0 \le x \le a \\
| |
| \tfrac{Pbx}{L}-P(x-a)=\tfrac{Pa(L-x)}{L}, & \mbox{for } a < x \le L
| |
| \end{cases}</math>
| |
| | <math>M_B = \tfrac{Pab}{L}</math>
| |
| |-
| |
| | <math>Q(x) = \begin{cases}
| |
| \tfrac{Pb}{L}, & \mbox{for } 0 \le x \le a \\
| |
| \tfrac{Pb}{L}-P, & \mbox{for } a < x \le L
| |
| \end{cases}</math>
| |
| | <math>Q_A = \tfrac{Pb}{L}</math>
| |
| <math>Q_C = \tfrac{Pa}{L}</math>
| |
| |-
| |
| | <math>w(x) = \begin{cases}
| |
| \tfrac{Pbx(L^2-b^2-x^2)}{6LEI}, & 0 \le x \le a \\
| |
| \tfrac{Pbx(L^2-b^2-x^2)}{6LEI}+\tfrac{P(x-a)^3}{6EI}, & a < x \le L
| |
| \end{cases}</math>
| |
| | <math>w_{\mathrm{max}} = \tfrac{\sqrt{3}Pb(L^2-b^2)^{\frac{3}{2}}}{27LEI}</math>
| |
| at <math>x = \sqrt{\tfrac{L^2-b^2}{3}}</math>
| |
| |}
| |
| | |
| === Cantilever beams ===
| |
| Another important class of problems involves [[cantilever]] beams. The bending moments (<math>M</math>), shear forces (<math>Q</math>), and deflections (<math>w</math>) for a cantilever beam subjected to a point load at the free end and a uniformly distributed load are given in the table below.<ref name=Gere/>
| |
| {| class="wikitable" style="margin:auto;"
| |
| |-
| |
| ! Distribution
| |
| ! Max. value
| |
| |-
| |
| | colspan="2" | '''Cantilever beam with end load'''
| |
| | rowspan="4" |[[Image:CantBeamPointLoad.svg|350px]]
| |
| |-
| |
| | <math>M(x) = P(x-L)</math>
| |
| | <math>M_A = PL</math>
| |
| |-
| |
| | <math>Q(x) = P </math>
| |
| | <math>Q_{\mathrm{max}} = P</math>
| |
| |-
| |
| | <math>w(x) = \tfrac{Px^2(3L-x)}{6EI} </math>
| |
| | <math>w_C = \tfrac{PL^3}{3EI}</math>
| |
| |-
| |
| | colspan="2" | '''Cantilever beam with uniformly distributed load'''
| |
| | rowspan="4" |[[Image:CantBeamDistLoad.svg|350px]]
| |
| |-
| |
| | <math>M(x) =-\tfrac{q(L^2-2Lx+x^2)}{2}</math>
| |
| | <math>M_A = \tfrac{qL^2}{2}</math>
| |
| |-
| |
| | <math>Q(x) = q(L-x), </math>
| |
| | <math>Q_A = qL</math>
| |
| |-
| |
| | <math>w(x) = \tfrac{qx^2(6L^2-4Lx+x^2)}{24EI} </math>
| |
| | <math>w_C = \tfrac{qL^4}{8EI}</math>
| |
| |}
| |
| | |
| Solutions for several other commonly encountered configurations are readily available in textbooks on mechanics of materials and engineering handbooks.
| |
| | |
| === Statically indeterminate beams ===
| |
| The [[bending moment]]s and [[shear force]]s in Euler-Bernoulli beams can often be determined directly using static balance of [[force]]s and [[moment (physics)|moment]]s. However, for certain boundary conditions, the number of reactions can exceed the number of independent equilibrium equations.<ref name=Gere/> Such beams are called ''[[statically indeterminate]]''.
| |
| | |
| The built-in beams shown in the figure below are statically indeterminate. To determine the stresses and deflections of such beams, the most direct method is to solve the Euler–Bernoulli beam equation with appropriate boundary conditions. But direct analytical solutions of the beam equation are possible only for the simplest cases. Therefore, additional techniques such as linear superposition are often used to solve statically indeterminate beam problems.
| |
| | |
| The superposition method involves adding the solutions of a number of statically determinate problems which are chosen such that the boundary conditions for the sum of the individual problems add up to those of the original problem.
| |
| | |
| {| class="wikitable" style="background:white; text-align:center;"
| |
| |-
| |
| |[[Image:Fem1.png|center|400px]]<br>(a) Uniformly distributed load ''q''. ||[[Image:Fem3.png|center|400px]]<br>
| |
| (b) Linearly distributed load with maximum ''q''<sub>0</sub>
| |
| |-
| |
| |<math>M_{\mathrm{max}} = \cfrac{qL^2}{12} ~;~~ w_{\mathrm{max}} = \cfrac{qL^4}{384EI}</math>||
| |
| |-
| |
| |[[Image:Fem2.png|center|400px]]<br>(c) Concentrated load P||[[Image:Fem4.png|center|400px]]<br>(d) Moment M<sub>0</sub>
| |
| |}
| |
| | |
| Another commonly encountered statically indeterminate beam problem is the [[cantilever|cantilevered beam]] with the free end supported on a roller.<ref name=Gere/> The bending moments, shear forces, and deflections of such a beam are listed below.
| |
| {| class="wikitable"
| |
| |-
| |
| ! Distribution
| |
| ! Max. value
| |
| | rowspan="4" |[[Image:Beam1svg.svg|350px]]
| |
| |-
| |
| | <math>M(x) = -\tfrac{q}{8}(L^2-5Lx+4x^2)</math>
| |
| | <math>\begin{align} M_B & = -\tfrac{9qL^2}{128} \mbox{ at } x=\tfrac{5L}{8}\\
| |
| M_A & = \tfrac{qL^2}{8} \end{align}</math>
| |
| |-
| |
| | <math>Q(x) = -\tfrac{q}{8}(8x -5L)</math>
| |
| | <math>Q_A = -\tfrac{5qL}{8}</math>
| |
| |-
| |
| | <math>w(x) = \tfrac{qx^2}{48EI}(3L^2-5Lx+2x^2)</math>
| |
| | <math>w_{\mathrm{max}} = \tfrac{qL^4}{185EI} \mbox{ at } x=0.5785L</math>
| |
| |}
| |
| | |
| == Extensions ==
| |
| | |
| The kinematic assumptions upon which the Euler–Bernoulli beam theory is founded allow it to be extended to more advanced analysis. Simple superposition allows for three-dimensional transverse loading. Using alternative [[constitutive equation]]s can allow for [[Viscoelasticity|viscoelastic]] or [[plastic bending|plastic]] beam deformation. Euler–Bernoulli beam theory can also be extended to the analysis of curved beams, [[Buckling|beam buckling]], composite beams, and geometrically nonlinear beam deflection.
| |
| | |
| Euler–Bernoulli beam theory does not account for the effects of transverse [[Shear stress|shear]] strain. As a result it underpredicts deflections and overpredicts natural frequencies. For thin beams (beam length to thickness ratios of the order 20 or more) these effects are of minor importance. For thick beams, however, these effects can be significant. More advanced beam theories such as the [[Timoshenko beam theory]] (developed by the Russian-born scientist [[Stephen Timoshenko]]) have been developed to account for these effects.
| |
| | |
| === Large deflections ===
| |
| [[File:EulerBernoulliBeam.png|right|thumb|400px|Euler-Bernoulli beam]]
| |
| The original Euler-Bernoulli theory is valid only for [[infinitesimal strain theory|infinitesimal strains]] and small rotations. The theory can be extended in a straightforward manner to problems involving moderately large rotations provided that the strain remains small by using the [[Theodore von Kármán|von Kármán]] strains.<ref name=Reddy>Reddy, J. N., (2007), ''Nonlinear finite element analysis'', Oxford University Press.</ref>
| |
| | |
| The Euler-Bernoulli hypotheses that plane sections remain plane and normal to the axis of the beam lead to displacements of the form
| |
| :<math>
| |
| v_1 = v_0(x) - z \cfrac{\mathrm{d}w_0}{\mathrm{d}x} ~;~~v_2 = 0 ~;~~ v_3 = w_0(x)
| |
| </math>
| |
| Using the definition of the Lagrangian Green strain from [[finite strain theory]], we can find the ''von Karman strains'' for the beam that are valid for large rotations but small strains. These strains have the form
| |
| :<math>
| |
| \begin{align}
| |
| \varepsilon_{11} & = \cfrac{\mathrm{d}u_0}{dx_1} - x_3\cfrac{\mathrm{d}^2w_0}{\mathrm{d}x_1^2} +
| |
| \frac{1}{2}\left[
| |
| \left(\cfrac{\mathrm{d}u_0}{\mathrm{d}x_1}-x_3\cfrac{\mathrm{d}^2w_0}{\mathrm{d}x_1^2}\right)^2 +
| |
| \left(\cfrac{\mathrm{d}w_0}{\mathrm{d}x_1}\right)^2\right] \\
| |
| \varepsilon_{22} & = 0 \\
| |
| \varepsilon_{33} & = \frac{1}{2}\left(\cfrac{\mathrm{d}w_0}{\mathrm{d}x_1}\right)^2 \\
| |
| \varepsilon_{23} & = 0 \\
| |
| \varepsilon_{31} & =
| |
| \frac{1}{2}\left(\cfrac{\mathrm{d}w_0}{\mathrm{d}x_1}-\cfrac{\mathrm{d}w_0}{\mathrm{d}x_1}\right) -
| |
| \frac{1}{2}\left[\left(\cfrac{\mathrm{d}u_0}{\mathrm{d}x_1}-x_3\cfrac{\mathrm{d}^2w_0}{\mathrm{d}x_1^2}\right)
| |
| \left(\cfrac{\mathrm{d}w_0}{\mathrm{d}x_1}\right)\right] \\
| |
| \varepsilon_{12} & = 0
| |
| \end{align}
| |
| </math>
| |
| From the [[principle of virtual work]], the balance of forces and moments in the beams gives us the equilibrium equations
| |
| :<math>
| |
| \begin{align}
| |
| \cfrac{\mathrm{d}N_{xx}}{\mathrm{d}x} + f(x) & = 0 \\
| |
| \cfrac{\mathrm{d}^2M_{xx}}{\mathrm{d}x^2} + q(x) +
| |
| \cfrac{\mathrm{d}}{\mathrm{d}x}\left(N_{xx}\cfrac{\mathrm{d}w_0}{\mathrm{d}x}\right) & = 0
| |
| \end{align}
| |
| </math>
| |
| where <math>f(x)</math> is the axial load, <math>q(x)</math> is the transverse load, and
| |
| :<math>
| |
| N_{xx} = \int_A \sigma_{xx}~ \mathrm{d}A ~;~~ M_{xx} = \int_A z\sigma_{xx}~ \mathrm{d}A
| |
| </math>
| |
| To close the system of equations we need the [[constitutive equations]] that relate stresses to strains (and hence stresses to displacements). For large rotations and small strains these relations are
| |
| :<math>
| |
| \begin{align}
| |
| N_{xx} & =
| |
| A_{xx}\left[\cfrac{\mathrm{d}u_0}{dx} + \frac{1}{2}\left(\cfrac{\mathrm{d}w_0}{\mathrm{d}x}\right)^2 \right] -
| |
| B_{xx}\cfrac{\mathrm{d}^2w_0}{\mathrm{d}x^2} \\
| |
| M_{xx} & =
| |
| B_{xx}\left[\cfrac{du_0}{\mathrm{d}x} + \frac{1}{2}\left(\cfrac{\mathrm{d}w_0}{\mathrm{d}x}\right)^2 \right] -
| |
| D_{xx}\cfrac{\mathrm{d}^2w_0}{\mathrm{d}x^2}
| |
| \end{align}
| |
| </math>
| |
| where
| |
| :<math>
| |
| A_{xx} = \int_A E~\mathrm{d}A ~;~~ B_{xx} = \int_A zE~\mathrm{d}A ~;~~ D_{xx} = \int_A z^2E~\mathrm{d}A ~.
| |
| </math>
| |
| The quantity <math>A_{xx}</math> is the ''[[stiffness|extensional stiffness]]'',<math>B_{xx}</math> is the coupled ''extensional-bending stiffness'', and <math>D_{xx}</math> is the ''[[bending stiffness]]''.
| |
| | |
| For the situation where the beam has a uniform cross-section and no axial load, the governing equation for a large-rotation Euler-Bernoulli beam is
| |
| <blockquote style="border: 1px solid black; padding:10px; width:400px">
| |
| :<math>
| |
| EI~\cfrac{\mathrm{d}^4 w}{\mathrm{d}x^4} - \frac{3}{2}~EA~\left(\cfrac{\mathrm{d} w}{\mathrm{d}x}\right)^2\left(\cfrac{\mathrm{d}^2 w}{\mathrm{d}x^2}\right) = q(x)
| |
| </math>
| |
| </blockquote>
| |
| | |
| == See also ==
| |
| * [[Applied mechanics]]
| |
| * [[Bending]]
| |
| * [[Bending moment]]
| |
| * [[Buckling]]
| |
| * [[Flexural rigidity]]
| |
| * [[Plate theory]]
| |
| * [[Sandwich theory]]
| |
| * [[Shear and moment diagram]]
| |
| * [[Singularity function]]
| |
| * [[Strain (materials science)]]
| |
| * [[Timoshenko beam theory]]
| |
| * [[Theorem of three moments]] (Clapeyron's theorem)
| |
| * [[Three point flexural test]]
| |
| | |
| == Notes ==
| |
| {{Reflist|group="N"}}
| |
| | |
| ==References==
| |
| {{Reflist}}
| |
| | |
| * {{cite conference | author=E.A. Witmer | title=Elementary Bernoulli-Euler Beam Theory | booktitle=MIT Unified Engineering Course Notes | year=1991-1992 | pages=5–114 to 5–164}}
| |
| | |
| {{DEFAULTSORT:Euler-Bernoulli Beam Equation}}
| |
| [[Category:Elasticity (physics)]]
| |
| [[Category:Solid mechanics]]
| |
| [[Category:Structural analysis]]
| |
| [[Category:Mechanical engineering]]
| |
| [[Category:Equations]]
| |