<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">WJM</journal-id><journal-title-group><journal-title>World Journal of Mechanics</journal-title></journal-title-group><issn pub-type="epub">2160-049X</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/wjm.2018.89024</article-id><article-id pub-id-type="publisher-id">WJM-87003</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Reaction Forces on a Fixed Ladder in Static Equilibrium: Analysis and Definitive Experimental Test of the Ladder Problem
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>M.</surname><given-names>P. Silverman</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Department of Physics, Trinity College, Hartford, USA</addr-line></aff><pub-date pub-type="epub"><day>29</day><month>08</month><year>2018</year></pub-date><volume>08</volume><issue>09</issue><fpage>311</fpage><lpage>342</lpage><history><date date-type="received"><day>2,</day>	<month>August</month>	<year>2018</year></date><date date-type="rev-recd"><day>27,</day>	<month>August</month>	<year>2018</year>	</date><date date-type="accepted"><day>30,</day>	<month>August</month>	<year>2018</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The development of a theoretical model to predict the four equilibrium forces of reaction on a simple ladder of non-adjustable length leaning against a wall has long remained an unresolved matter. The difficulty is that the problem is statically indeterminate and therefore requires complementary information to obtain a unique solution. This paper reports 1) a comprehensive theoretical analysis of the three fundamental models based on treating the ladder as a single Euler-Bernoulli beam, and 2) a detailed experimental investigation of the forces of reaction as a function of applied load and location of load. In contrast to previous untested proposals that the solution to the ladder problem lay in the axial constraint on compression or the transverse constraint on flexure, the experimental outcome of the present work showed unambiguously that 1) the ladder could be modeled the best by a pinned support at the base (on the ground) and a roller support at the top (at the wall), and 2) the only complementary relation needed to resolve the static indeterminacy is the force of friction at the wall. Measurements were also made on the impact loading of a ladder by rapid ascent and descent of a climber. The results obtained were consistent with a simple dynamical model of the ladder as a linear elastic medium subject to a pulse perturbation. The solution to the ladder problem herein presented provides a basis for theoretical extension to other types of ladders. Of particular importance, given that accidents involving ladders in the workplace comprise a significant fraction of all industrial accidents, the theoretical relations reported here can help determine whether a collapsed structure, against which a ladder was applied, met regulatory safety limits or not.
 
</p></abstract><kwd-group><kwd>Reaction Forces on a Ladder</kwd><kwd> Reaction Forces on a Beam</kwd><kwd> Impulse Load on a Ladder</kwd><kwd> Ladder in Static Equilibrium</kwd><kwd> Statically Indeterminate Forces</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction: The Ladder Problem</title><p>Although determination of the reaction forces on a beam is a standard part of the mechanics of continuous media, the exact conditions under which an actual physical structure comprised of beams is stressed are not always apparent. Models based on different assumptions can predict hugely different reaction forces, and careful experimentation is required to determine which model may most accurately characterize a structure. One of the most important mechanical structures for which reaction forces need to be determined is that of a ladder. Although ladders have been in use since antiquity, to date no definitive model of the reaction forces on a ladder has been experimentally tested and confirmed. This paper reports what the author believes to be the most comprehensive theoretical analysis and definitive experimental test of the reaction forces on a fixed ladder.</p><sec id="s1_1"><title>1.1. The Fixed Ladder as an Euler-Bernoulli Beam</title><p>A fixed ladder is defined by the American Ladder Institute as a non-self-supporting ladder of fixed (i.e. non-adjustable) length [<xref ref-type="bibr" rid="scirp.87003-ref1">1</xref>] . It comprises two identical rails of uniform density connected by short rungs. As a suitable first approximation to modeling such a ladder, the rungs are ordinarily disregarded and the two rails are merged into a single uniform beam, such as shown in panel A of <xref ref-type="fig" rid="fig1">Figure 1</xref>. The single-beam ladder of length L is in static equilibrium supported obliquely by the flat horizontal ground and a flat vertical wall. The weight W of the ladder is regarded as a uniformly distributed load of linear density</p><p>w = W / L (1)</p><p>that acts vertically downward at the center of mass of any isolated segment of the beam. A climber is modeled as a point load P acting vertically downward at a distance β L ( 1 ≥ β ≥ 0 ) along the beam from the contact point at the ground, which is taken to be the origin.</p><p>As shown in panel A of the figure,there are four reaction forces:1) normal force at the ground ( R 1 ) , 2) parallel force at the ground ( R 2 ) , 3) normal force at the wall ( R 3 ) , and 4) parallel force at the wall ( R 4 ) . The prediction of these four forces constitutes what is called the ladder problem. Before discussing the theoretical solutions and experimental test of the ladder problem, it is worthwhile to consider why the problem is a significant one. This significance is both conceptual and practical.</p><p>First, conceptually, the solution to any model serving as an archetype for a class of real-world problems within some branch of physics is significant by virtue of its fundamentality. The single-beam ladder is an archetype within the mechanics of continuous media. Moreover, understanding the forces on a fixed</p><p>ladder would provide useful modeling guidelines for predicting the reaction forces on more complex ladder types, such as extension ladders, step ladders, trestle ladders, and others. Second, from a practical standpoint, the ability to predict the forces on ladders has direct impact on matters relating to safety in the use of ladders and to legal issues that arise when use of a ladder in the workplace results in injury or death. For example, with a valid fixed ladder model, the analyst can determine the acceptable range of inclination angles to avoid slippage, or determine whether the net force exerted by a ladder on a wall or railing that collapsed was within or exceeded regulatory limits [<xref ref-type="bibr" rid="scirp.87003-ref2">2</xref>] .</p><p>One might think that prediction of the reaction forces on a structure as simple as a single-beam ladder is a relatively straightforward matter. Indeed, the problem is trivially solvable if the parallel reaction at the wall is ignored, which is the case in countless physics textbooks (see, for example, [<xref ref-type="bibr" rid="scirp.87003-ref3">3</xref>] ) and expository internet articles. The fixed ladder leaning against a frictionless wall has been an iconic example of static equilibrium in elementary mechanics books for at least a century. This problem is statically determinate; i.e. the three reaction forces ( R 1 , R 2 , R 3 ) are uniquely solvable from the three equations of static equilibrium in which the net horizontal force, net vertical force, and rotational moment about any point at rest all vanish. The resulting solution, however, is of academic interest only and does not describe correctly the interactions of real ladders with real surfaces.</p><p>When account is taken of the parallel reaction at the wall, the problem becomes statically indeterminate, and therefore more complicated, since there are now 4 unknowns and 3 equations. One proposed solution assumed that the ladder undergoes an axial deformation which is constrained by the supports [<xref ref-type="bibr" rid="scirp.87003-ref4">4</xref>] . This model cannot be correct as it stands because the transverse (bending) deformation―as calculated in Section 2―is approximately two orders of magnitude greater than the axial deformation and should not have been neglected. Another proposal [<xref ref-type="bibr" rid="scirp.87003-ref5">5</xref>] took account of both axial compression and transverse flexion, and concluded that the problem could not be solved by static analysis because it was not possible to determine which deformation dominated. This conclusion is not substantiated by the findings of the present paper. Neither of these proposals provided experimental measurements to test their validity.</p><p>Experimentally, there have been many investigations, of which some representative examples are [<xref ref-type="bibr" rid="scirp.87003-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.87003-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.87003-ref8">8</xref>] , to elucidate the conditions under which ladders can be used safely. In general, these investigations focused on matters relating to friction, angle of inclination, climbing speed, weight of climber, and other empirically accessible variables, but did not attempt to provide a model that would succinctly incorporate these data in a mathematical theory.</p><p>This paper reports the theoretical solution of three basic models by which to analyze the reaction forces of a fixed ladder represented as a single Euler-Bernoulli (E-B) beam. The characteristics of an E-B beam pertinent to the present study are as follows [<xref ref-type="bibr" rid="scirp.87003-ref9">9</xref>] :</p><p>1) the length of the beam is much greater than the linear dimension of the cross-section;</p><p>2) the long axis of the beam lies within the neutral surface;</p><p>3) the cross section of the beam remains plane and perpendicular to the long axis during deflection;</p><p>4) deformation of the cross section within its own plane is neglected;</p><p>5) the beam is a linear elastic material subject to Hooke’s law; i.e. normal stress within a cross section varies linearly with perpendicular distance from the neutral axis.</p><p>The neutral surface of a beam is the interface that separates the fibers under compression from the fibers under tension when the beam is deflected transversely. The neutral axis within any cross section is the line of intersection of the neutral surface with the cross section.</p></sec><sec id="s1_2"><title>1.2. The Analytical Models</title><p>The three basic ladder models analyzed in this paper are distinguished by the nature of the supports at the points of contact with the surfaces of the ground and wall. In the mechanics of beams, the three common types of supports are 1) roller, 2) pin, and 3) fixed [<xref ref-type="bibr" rid="scirp.87003-ref10">10</xref>] . The names refer to the constraints on the motion of the affected segment of the beam, and not necessarily to the actual presence of a roller, pin, or clamping mechanism. Thus, if one end of a beam is supported by a roller, it can displace laterally, but not vertically. If the end of the beam is pinned, it cannot displace either laterally or vertically, but it can rotate about the pin. And, finally, if the end is fixed, it cannot displace or rotate at all. [There is a fourth kind of support usually designated as simple, which is idealized as a frictionless surface. This would apply to the support of the single-beam ladder at a wall as treated in the elementary physics books.]</p><p>If a right-handed Cartesian coordinate system ( x , y , z ) is chosen such that an initially horizontal beam lies along the x axis, and a load is applied that can deflect the beam vertically downward along the negative y axis, then the boundary and continuity conditions at the three kinds of supports can be summarized as follows</p><p>Roller :       Δ y = 0 ,       d y / d x ≠ 0 ,     M z = 0 Pin :                 Δ y = 0 ,       d y / d x ≠ 0 ,     M z = 0 Fixed :         Δ y = 0 ,       d y / d x = 0 ,     M z ≠ 0 (2)</p><p>where the symbol Δ signifies displacement, and M z is an internal bending moment (torque) about the z axis. The salient features of relations (2) are that a) the displacement and slope of the deflected beam are zero at a fixed support but the support creates an internal bending moment, whereas b) at roller and pin supports there is no internal moment, and the slope of the bending curve there need not be zero.</p><p>The models analyzed in this paper are based on the supports described by Equation (2).</p><p>・ In Model 1, the ground support of the single-beam ladder is fixed and the wall support is equivalent to a pin. Under a vertical load due to the distributed weight of the ladder and the point weight of a hypothetical climber, the ladder is subject to both axial and transverse stresses, but is constrained from compression, elongation, or rotation at the supports. This is a case of static indeterminateness in two variables: one axial and one transverse. Although friction is presumably the source of the constraints, no specific assumptions are made in regard to the mathematical form of the frictional force.</p><p>・ In Model 2, the ground support of the single-beam ladder is fixed, and the wall support is equivalent to a roller. The axial component of the load can compress the ladder, and thereby remove the static indeterminateness along the length. The reaction at the wall is assumed to be governed by the phenomenological relation for friction taken to be proportional to the normal force at the wall [<xref ref-type="bibr" rid="scirp.87003-ref11">11</xref>] . There remains a static indeterminateness in one of the transverse reactions because of the unknown internal moment at the ground support.</p><p>・ In Model 3, the ground support of the single-beam ladder is equivalent to a pin and the wall support is equivalent to a roller. The ladder is free to compress under an axial force and flex about the ground support, thereby relieving all static constraints. The reaction at the wall is assumed governed by friction as in Model 2. This is a case of static indeterminateness in one variable in which the equations of static equilibrium, complemented by the frictional force law, suffice to determine a unique solution.</p><p>The three models are schematically illustrated in <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p></sec><sec id="s1_3"><title>1.3. Outline of Paper</title><p>In Section 2 the reaction forces on a ladder characterized as a single E-B beam are calculated for each model and examined as a function of load, location of load along the ladder, and angle of inclination of the ladder with respect to the ground. In the case of Models 1 and 3, the full bending curve is calculated within the Euler-Bernoulli approximation (neglect of the square of the slope d y / d x ), and the magnitudes of transverse and axial deformations are compared. (The bending curve of Models 2 and 1 are the same).</p><p>In Section 3 are reported experiments using horizontal and vertical force platforms to measure the reactions on a fixed metal ladder inclined against a wall. It is noted here briefly that the model that was found to account for the observed results was Model 3.</p><p>In Section 4 the impact loading on a ladder by a climber is examined theoretically and experimentally.</p><p>In Section 5 are reported measurements of the reactions on a wood beam inclined obliquely against a wall. The intent of the experiment was to ascertain whether the same theoretical model that best applied to a ladder modeled as a beam applied as well to an actual single beam.</p><p>In Section 6 conclusions drawn from this research are summarized.</p></sec></sec><sec id="s2"><title>2. Analysis of a Single-Beam Fixed Ladder</title><p>Panel A of <xref ref-type="fig" rid="fig1">Figure 1</xref></p><p>A 1 = R 1 sin θ + R 2 cos θ B 1 = R 1 cos θ − R 2 sin θ A 2 = − R 3 cos θ + R 4 sin θ B 2 = R 3 sin θ + R 4 cos θ (3)</p><p>with inverse relations</p><p>R 1 = B 1 cos θ + A 1 sin θ R 2 = − B 1 sin θ + A 2 cos θ R 3 = − B 2 sin θ − A 2 cos θ R 4 = B 2 cos θ + A 2 sin θ (4)</p><p>The equations of static equilibrium then take either of the forms</p><p>(I) R 1 + R 4 = W + P R 2 − R 3 = 0 (II) B 1 + B 2 = W n + P n A 1 + A 2 = W a + P a (5)</p><p>in which the normal and axial components of the ladder weight and a point load representative of a climber are respectively</p><p>W n = W cos θ               W a = W sin θ P n = P cos θ                   P a = P sin θ . (6)</p><p>The sets of equations in (5) do not contain an equation for the vanishing of a rotational moment about a designated point at rest. In the absence of information regarding the supports of the ladder at the ground, it is not possible to write a torque equation based on static equilibrium alone. This is a consequence of the fact that the system may contain an initially unknown internal moment created by the ground support in order to maintain static equilibrium.</p><p>Since there are four unknown reactions in Equation (5) and fewer than four equations of static equilibrium, the set of equations is statically indeterminate and must be supplemented by what are termed complementary conditions. Under the circumstances stated above, these conditions can be applied to the axial and normal directions independently. The complementary conditions take different forms depending on the specific model, as discussed in the following subsections.</p><sec id="s2_1"><title>2.1. Model 1: Fixed Ground and Pinned Wall Supports</title><p>With respect to the schematic diagram in panel B of <xref ref-type="fig" rid="fig1">Figure 1</xref>, Model 1 entails a fixed support at the left end x = 0 (with initially unknown internal moment M z ( 0 ) ) and a pin support at the right end x = L . The total strain energy U of the beam</p><p>U = U a + U b (7)</p><p>is the sum of the axial strain energy</p><p>U a = 1 2 E A ∫ 0 L F a ( x ) 2 d x (8)</p><p>due to the compression or tension by force F a ( x ) and the energy of flexure (bending)</p><p>U b = 1 2 E I z ∫ 0 L M z ( x ) 2 d x (9)</p><p>due to the moment M z ( x ) about the z axis (see Ref. [<xref ref-type="bibr" rid="scirp.87003-ref15">15</xref>] , pp 188-196). The material parameters in Equations (8) and (9) are the beam cross section area A, beam area moment of inertia I z about the transverse z axis, and the elastic modulus (Young’s modulus) E. For a beam of width w and height h as portrayed in panel C of <xref ref-type="fig" rid="fig1">Figure 1</xref>, the cross section and area moment of inertia take the forms (see [<xref ref-type="bibr" rid="scirp.87003-ref15">15</xref>] , pp 448-449)</p><p>A = w h (10)</p><p>I z = 1 12 w h 3 . (11)</p><p>Because the overall load comprises both distributed and point contributions, the force F a ( x ) and moment M z ( x ) are piecewise continuous functions. These functions can be represented economically by use of the Heaviside unit step function [<xref ref-type="bibr" rid="scirp.87003-ref16">16</xref>] defined here as</p><p>H ( x − x 0 ) = { 0           x &lt; x 0 1           x ≥ x 0 . (12)</p><p>Mathematically, the Heaviside function is usually defined by pure inequalities in both partitions of the real axis and left undefined at the point where the argument is zero. In the mechanics of continuous media, however, the point load occurs at the partition between segments and therefore it is more useful to define the Heaviside function as unity there. One can then express the axial force and bending moment in the following ways</p><p>F a ( x ) = A 2 − W a L ( L − x ) − P a H ( β L − x ) (13)</p><p>M z ( x ) = B 2 ( L − x ) − W n 2 L ( L − x ) 2 − P n ( β L − x ) H ( β L − x ) . (14)</p><p>From the complementary energy theorem (Castigliano’s theorem) of continuum mechanics (see [<xref ref-type="bibr" rid="scirp.87003-ref15">15</xref>] pp 201-217) for systems where the force-deformation relation is linear, the deflection of a point i on a beam in the direction of a load F i applied at that point is given by</p><p>δ i = ∂ U ∂ F i . (15)</p><p>Since a physical condition of Model 1 is that the right endpoint of the beam in panel B of <xref ref-type="fig" rid="fig1">Figure 1</xref> suffers no deflection, it follows that</p><p>δ L ( a ) = ∂ U ∂ A 2 = 0 (16)</p><p>and</p><p>δ L ( b ) = ∂ U ∂ B 2 = 0 (17)</p><p>where superscript a signifies axial and superscript b signifies bending. Equations (16) and (17) are the complementary relations required to supplement the set of relations (5) and resolve the static indeterminacy of Model 1.</p><p>Substitution of relations (13) and (14) respectively into strain energies (8) and (9), taking derivatives (16) and (17), and setting the resulting integrals to zero lead to expressions for the axial reaction A 2 and transverse reaction B 2 . The reactions A 1 and B 1 are then obtained from Equation (5) of static equilibrium. The resulting set of four axial and transverse reactions are found to be</p><p>A 1 = 1 2 W sin θ + ( 1 − β ) P sin θ A 2 = 1 2 W sin θ + β P sin θ B 1 = 5 8 W cos θ + ( 1 − 3 2 β 2 + 1 2 β 3 ) P cos θ B 2 = 3 8 W cos θ + ( 3 2 β 2 − 1 2 β 3 ) P cos θ (18)</p><p>from which follow by transformation (4) the reactions normal and parallel with respect to the ground and wall</p><p>R 1 = 1 2 W ( 1 + 1 4 cos 2 θ ) + P [ 1 − ( 3 2 β 2 − 1 2 β 3 ) cos 2 θ − β sin 2 θ ] R 2 = R 3 = [ 1 16 W + ( 1 2 β − 3 4 β 2 + 1 4 β 4 ) P ] sin 2 θ R 4 = 1 2 W ( 1 − 1 4 cos 2 θ ) + P [ ( 3 2 β 2 − 1 2 β 3 ) cos 2 θ + β sin 2 θ ] . (19)</p><p>The directions of R 2 and R 3 are shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>, and no algebraic sign denoting direction is needed since these two forces are not superposed in any calculation in this paper.</p><p><xref ref-type="fig" rid="fig3">Figure 3</xref> shows the variation in normal force at the ground (panel A) and at the wall (panel B) as a function of the fractional location ( 1 ≥ β ≥ 0 ) of the applied load along the beam for different load magnitudes ranging from 0 to 100 pounds (lbs). The ladder weight is taken to be 25 lbs (111.2 N), in accordance with the use of English units in ladder specifications given in the US, but calculations are performed and recorded in metric units (newtons N). The angle of inclination was taken to be 60˚. Plots of the parallel force at the ground ( R 2 ) and at the wall ( R 4 ) are not shown since by relations (18) the former are the same as the plots in panel B, and the latter correspond to reversing the horizontal scale (abscissa) in panel A.</p><p>A complementary perspective is given in <xref ref-type="fig" rid="fig4">Figure 4</xref>, which shows the variation in normal force at the ground (panel A) and at the wall (panel B) as a function of the angle of inclination for different fractional locations of the applied load along the beam. According to Model 1, the location of a hypothetical climber (of</p><p>fixed weight) influences the normal reaction at the wall strongly and the normal reaction at the ground relatively weakly as the inclination of the ladder is changed.</p><p>Applying Equation (15) to the total strain energy U, one can determine the deformation at any point in the elastic curve of the single-beam ladder and thereby compare the magnitudes of the axial compression δ ( a ) and transverse</p><p>(bending) deflection δ ( b ) . Taking β = 1 2 for illustration yields the deformations</p><p>δ 1 2 L ( a ) = ∂ U ∂ P a | β = 1 2 = L 2 A E ( W + P ) sin θ (20)</p><p>δ 1 2 L ( b ) = ∂ U ∂ P n | β = 1 2 = L 3 768 E I z ( 4 W + 7 P ) cos θ . (21)</p><p>Given typical parameters of beam length L = 2.5 m, width w = 2 cm, height h = 5 cm, inclination angle θ = 45 ∘ , ladder weight W = 25 lbs (111.2 N), climber weight P = 150 lbs (667.2 N), and elastic modulus of aluminum</p><p>E = 69 GPa, Equations (20) and (21) lead to deformations δ 1 2 L ( a ) ≈ 10 &#181;m, δ 1 2 L ( b ) ≈ 5.1 mm. The deformation due to bending is approximately 500 times the</p><p>deformation due to axial compression. Therefore, any model of a ladder subject to ordinary loads that takes account of axial deformation cannot neglect deformation due to flexure.</p><p>As a final consideration, it is useful to derive the deflection curve of Model 1 as a check on the calculated deformations and to ascertain that the reactions of a single-beam ladder derived from Model 1 correctly fulfill the complementary conditions. The exact flexure formula for an E-B beam relates the bending moment to the radius of curvature and takes the form of a nonlinear second-order differential equation (see [<xref ref-type="bibr" rid="scirp.87003-ref10">10</xref>] pp. 576-583)</p><p>d 2 y / d x 2 [ 1 + ( d y / d x ) 2 ] 3 / 2 = 1 E I z M z ( x ) (22)</p><p>whose solution is called the elastica [<xref ref-type="bibr" rid="scirp.87003-ref17">17</xref>] . Although an exact solution would be required for a highly flexible beam such as a fishing rod [<xref ref-type="bibr" rid="scirp.87003-ref9">9</xref>] , the deflection of the rail of a ladder is sufficiently low that one can neglect the square of the slope in the denominator of Equation (22), as is conventionally done in engineering practice. The equation then becomes linear and can be readily solved subject to the imposed boundary conditions</p><p>y ( 0 ) = y ( L ) = 0 y ( x ) i s c o n t i n u o u s a t x = β L d y / d x iscontinuousat x = β L (23)</p><p>to yield the piecewise continuous function</p><p>y ( x ) = y 1 ( x ) H ( β L − x ) + y 2 ( x ) H ( x − β L ) (24)</p><p>with</p><p>y 1 ( x ) Model1 = 1 E I z [ − L ( 1 16 W n + ( 1 2 β − 3 4 β 2 + 1 4 β 3 ) P n ) x 2                                                           + ( 5 48 W n + 1 6 ( 1 − 3 2 β 2 + 1 2 β 3 ) P n ) x 3 − W n 24 L x 4 ] y 2 ( x ) Model1 = 1 E I z [ 1 6 β 3 L 3 P n − 1 2 ( β 2 L 2 P n ) x + L ( − 1 16 W n + ( 3 4 β 2 − 1 4 β 3 ) P n ) x 2 ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​                                                         + ( 5 48 W n − 1 6 ( 3 2 β 2 − 1 2 β 3 ) P n ) x 3 − W n 24 L x 4 ] . (25)</p><p><xref ref-type="fig" rid="fig5">Figure 5</xref> shows plots of the deflection curve y ( x ) for different locations of a 150 lb load on a ladder 2.5 m in length and 25 lbs in weight inclined at an angle of 45˚ against a vertical wall. The cross section of the single beam is 2   cm &#215; 5   cm and the modulus of elasticity is E = 69 GPa. As seen in the figure, the deflection curve is the same for a load directly at the base x = 0 or at the top of the</p><p>ladder x = L . The deflection δ 1 2 L ( b ) ∼ 5.1 mm occurs at the mid-position β = 1 2 ,</p><p>which is the result predicted by Equation (21) obtained by Castigliano’s theorem. The deflection curve is not symmetric for any location of the load, even at mid-position, because the slope is required to be 0 at the ground support (left end) but not at the wall support (right end).</p><p>Calculation of deflections (20), (21) and the elastica curve (24) involved</p><p>integration of expressions of the form x n H ( x − x 0 ) for non-negative integer n. The integration of Heaviside functions is summarized in the Appendix.</p></sec><sec id="s2_2"><title>2.2. Model 2: Fixed Ground and Roller Wall Supports</title><p>With respect to the schematic diagram in panel B of <xref ref-type="fig" rid="fig1">Figure 1</xref>, Model 2 entails a fixed support at the left end and a roller support at the right end. Resolution of the static indeterminateness of the bending reactions is accomplished in exactly the same way as for Model 1, and therefore forces B 1 and B 2 are given as shown in Equations (18). To determine the other forces, one makes use of the transformation (3) for B 2 and the phenomenological relation for friction</p><p>R 4 = μ R 3 (26)</p><p>withcoefficient of friction μ to obtain</p><p>B 2 = R 3 sin θ + R 4 cos θ = R 3 ( sin θ + μ cos θ ) , (27)</p><p>fromwhich follows</p><p>R 3 = R 2 = R 4 μ = 3 8 W + ( 3 2 β 2 − 1 2 β 3 ) P tan θ + μ (28)</p><p>andtherefore from Equation (5)</p><p>R 1 = W + P − R 4 = W + P − μ R 3 = W + P − μ   [ 3 8 W + ( 3 2 β 2 − 1 2 β 3 ) P tan θ + μ ] . (29)</p><p>The axial forces A 1 and A 2 are given by the transformation (3)</p><p>A 1 = ( W + P ) sin θ + ( 3 8 W + ( 3 2 β 2 − 1 2 β 3 ) P cos θ ) [ 1 − μ tan θ tan θ + μ ] A 2 = − ( 3 8 W + ( 3 2 β 2 − 1 2 β 3 ) P cos θ ) [ 1 − μ tan θ tan θ + μ ] (30)</p><p>ratherthan by a complementary condition constraining axial compression.</p><p><xref ref-type="fig" rid="fig6">Figure 6</xref>shows the variation in normal force at the ground (panel A) and at the wall (panel B) as a function of load location for different loads under the same conditions as for <xref ref-type="fig" rid="fig3">Figure 3</xref>. The coefficient of friction was taken to be 0.4. Correspondingly, <xref ref-type="fig" rid="fig7">Figure 7</xref>shows the variation in these two reactions as a function of angle of inclination with the ground under the same conditions as for <xref ref-type="fig" rid="fig4">Figure 4</xref>. Models 1 and 2 are seen to predict reactions that depend very differently on load location.</p><p>Since the reaction forces normal to the beam ( B 1 , B 2 ) are the same as for Model 1, the bending deflection curve, which does not depend on the axial reactions, is also the same.</p></sec><sec id="s2_3"><title>2.3. Model 3: Pinned Ground and Roller Wall Supports</title><p>With respect to the schematic diagram in panel B of <xref ref-type="fig" rid="fig1">Figure 1</xref>, Model 3 entails a</p><p>pinned support at the left end and a roller support at the right end. Because there is no unknown internal moment at pinned and roller supports, the set of equations (5) can be supplemented by an equation for the vanishing of the external moment M z about any fixed point. For example, if the pivot is taken to be the point of contact with the ground, then the vanishing of M z ( 0 ) can be expressed by either of the equations</p><p>R 3 sin θ + R 4 cos θ = 1 2 W cos θ + β P cos θ (31)</p><p>or</p><p>B 2 = 1 2 W n + β P n . (32)</p><p>Together with the complementary relation (26) for force of friction and the other equations of static equilibrium, one obtains</p><p>R 3 = R 2 = R 4 μ = 1 2 W + β P tan θ + μ (33)</p><p>B 1 = 1 2 W n + ( 1 − β ) P n (34)</p><p>and</p><p>R 1 = W + P − μ R 3 = W ( tan θ + 1 2 μ ) + P ( tan θ + ( 1 − β ) μ ) tan θ + μ . (35)</p><p>From the set of transformations (3) follow the axial reactions</p><p>A 1 = ( W + P ) sin θ + ( 3 8 W + ( 3 2 β 2 − 1 2 β 3 ) P cos θ ) [ 1 − μ tan θ tan θ + μ ] A 2 = − ( 3 8 W + ( 3 2 β 2 − 1 2 β 3 ) P cos θ ) [ 1 − μ tan θ tan θ + μ ] ​ ​ (36)</p><p><xref ref-type="fig" rid="fig8">Figure 8</xref> and <xref ref-type="fig" rid="fig9">Figure 9</xref> respectively show the variation in normal force at the ground (panel A) and at the wall (panel B) as a function of load location along the ladder and as a function of inclination angle with the ground. The plots are similar to, but not numerically identical with, the plots of <xref ref-type="fig" rid="fig6">Figure 6</xref> and <xref ref-type="fig" rid="fig7">Figure 7</xref>.</p><p>The bending curve (elastica) of Model 3 takes a different shape than that of Models 1 and 2 because of the difference in the nature of the support at the ground. Solving the Euler-Bernoulli approximation to the flexure equation (22)</p><p>d 2 y / d x 2 = M z ( x ) / E I z (37)</p><p>with moment (see panel B of <xref ref-type="fig" rid="fig1">Figure 1</xref>)</p><p>M z ( x ) = B 1 x − W n 2 L x 2 − P n ( x − β L ) H ( x − β L ) (38)</p><p>and boundary conditions (23) leads to the piecewise continuous function (24) with components</p><p>y 1 ( x ) Model3 = 1 E I z [ − L 2 ( 1 24 W n + ( 1 3 β − 1 2 β 2 + 1 6 β 3 ) P n ) x                                                           + ( 1 12 W n + 1 6 ( 1 − β ) P n ) x 3 − W n 24 L x 4 ] y 2 ( x ) Model3 = 1 E I z [ β 3 L 3 6 P n − L 2 ( 1 24 W n + ( 1 3 β + 1 6 β 3 ) P n ) x + ( 1 2 β L P n ) x 2 ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​                                                       + ( 1 12 W n − 1 6 β P n ) x 3 − W n 24 L x 4 ] (39)</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>0 shows plots of the deflection curve y ( x ) for different locations of a 150 lb load on a ladder 2.5 m in length and 25 lbs in weight inclined at an angle of 45˚ against a vertical wall. As seen in the figure, the deflection curve is again the same for a load directly at the base x = 0 or at the top of the ladder x = L . In contrast to the deflection curves of Models 1 and 2, the deflection curve of Model 3 is symmetric for a load applied at the middle of the beam; the slopes are nonzero and equal at both the ground and wall contact points. For other locations of the load, the deflection curve is no longer symmetric, but at no location are the slopes zero at the end points.</p><p>The maximum deflection δ 1 2 L ( b ) ∼ 11.8 mm occurs at the mid-position</p><p>β = 1 / 2 . This is the same result one would get from applying Castigliano’s theorem, as implemented in Section 2.1, which leads to the deformation</p><p>δ 1 2 L ( b ) = L 3 384 E I z ( 5 W + 8 P ) cos θ . (40)</p><p>The deflection produced by the axial forces of Equation (36) is given by</p><p>δ β L ( a ) = L A E ( A 2 − ( 1 2 W + β P ) sin θ ) ( 1 − μ tan θ tan θ + μ ) . (41)</p><p>For the load at mid-position and the other parameters of <xref ref-type="fig" rid="fig1">Figure 1</xref>0, Equation (41) leads to an axial compression 13.3 &#181;m. Thus, as was the case with Model 1, the bending deflection is again seen to be at least two orders of magnitude greater than the corresponding deformation due to axial compression.</p></sec></sec><sec id="s3"><title>3. Measured Reaction Forces on a Fixed Ladder</title><p>To determine which, if any, of the three basic models describes most closely the reaction forces on a fixed ladder, an experiment was performed to measure forces R 1 , R 2 , and R 3 . The experiment was performed as follows. A ladder of length 244 cm and weight 18.9 lbs (84.07 N) was placed against a vertical wall at an angle of 68˚ to the ground. A two-axis force platform (Pasco Model PS-2142) was situated horizontally on the ground beneath each of the two rails of the ladder. The top of the ladder rested on a third force platform attached vertically to the wall. Each force platform, used in conjunction with a computer and the appropriate interface, was designed to measure the applied force normal to its surface (up to a maximum of 4400 N) and along an axis parallel to its surface (up to 1100 N).</p><p>The ladder comprised seven rungs located at fractional lengths along the rails of</p><p>β = ( 32 244 , 62 244 , 92 244 , 122 244 , 152 244 , 182 244 , 212 244 ) ≈ ( 0.131 , 0.254 , 0.377 , 0.500 , 0.623 , 0.746 , 0.869 )</p><p>of which the first six rungs were used in the experiment. External loads were applied to the rungs in the form of lead bricks, each brick weighing 25 lbs (111.21 N). Measurements of the normal and parallel forces were made for discrete loads of 1, 2, 3, or 4 bricks placed on each rung in progressive order from rung 1 to rung 6. The sum of the forces normal to the horizontal surface of force platforms 1 and 2 (supporting the rails) constituted reaction R 1 . The sum of the forces parallel to the surface of platforms 1 and 2 constituted reaction R 2 . The force normal to the vertical surface of force platform 3 (against which the top of the ladder rested) constituted reaction R 3 . Although the initial intention was to measure reaction R 4 as well, this measurement could not be made because the force platforms did not measure shear forces accurately when the platform was oriented vertically<sup>1</sup>. The omitted reaction, however, was not needed to establish the conclusions drawn from this experiment.</p><p><xref ref-type="table" rid="table1">Table 1</xref> records the measurements of the 3 reactions ( R 1 , R 2 , R 3 ) at each of 6 rungs for each of 4 loads (25, 50, 75, 100 lbs). The experimental forces are respectively compared with the corresponding theoretical predictions of the three models in Figures 11-13.</p><p>In order to apply the theoretical relations of Models 2 and 3, it is necessary to know the coefficient of friction μ of the ladder with the wall―or, more accurately, with the surface of the force platform attached to the wall. This coefficient cannot be calculated theoretically, but must be obtained empirically from the data. The method employed here was trial and inspection; i.e. to assume a value for μ , calculate and plot the resulting reactions, and adjust as necessary the value of μ until a satisfactory match with data was achieved, provided that the model itself is a valid representation of the underlying mechanical processes. The entire procedure can be executed in a matter of seconds with a fast computer and available mathematical software such as Maple or Mathcad. It is to be emphasized that the totality of the 72 measurements (3 reactions &#215; 4 loads &#215; 6 load locations) tightly constrained the value of μ . For Model 2 or 3 to be valid it is necessary that a single value of μ suffice to account for all the data. If different statistically significant values of μ for the various reactions are needed to match the data, or if no value of μ matches the data, then Model 2 or 3 would not describe the static equilibrium of a ladder against a wall.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Experimental Normal and Parallel Reactions on a Ladder</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >θ = 68.0 ∘</th><th align="center" valign="middle" >LOAD (lbs)</th><th align="center" valign="middle" >25</th><th align="center" valign="middle" >50</th><th align="center" valign="middle" >75</th><th align="center" valign="middle" >100</th></tr></thead><tr><td align="center" valign="middle" >REACTIONS (N)</td><td align="center" valign="middle" >RUNG</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >Ground Normal R 1</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >347.81</td><td align="center" valign="middle" >303.74</td><td align="center" valign="middle" >412.54</td><td align="center" valign="middle" >526.04</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >188.54</td><td align="center" valign="middle" >302.98</td><td align="center" valign="middle" >413.19</td><td align="center" valign="middle" >526.67</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >193.31</td><td align="center" valign="middle" >304.92</td><td align="center" valign="middle" >409.82</td><td align="center" valign="middle" >522.06</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >189.95</td><td align="center" valign="middle" >303.37</td><td align="center" valign="middle" >413.52</td><td align="center" valign="middle" >520.52</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >187.31</td><td align="center" valign="middle" >300.74</td><td align="center" valign="middle" >409.77</td><td align="center" valign="middle" >516.84</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >6</td><td align="center" valign="middle" >188.50</td><td align="center" valign="middle" >298.47</td><td align="center" valign="middle" >403.50</td><td align="center" valign="middle" >508.11</td></tr><tr><td align="center" valign="middle" >Ground Parallel R 2</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >20.19</td><td align="center" valign="middle" >23.87</td><td align="center" valign="middle" >31.56</td><td align="center" valign="middle" >37.00</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >22.92</td><td align="center" valign="middle" >35.01</td><td align="center" valign="middle" >47.29</td><td align="center" valign="middle" >58.86</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >32.47</td><td align="center" valign="middle" >48.67</td><td align="center" valign="middle" >63.63</td><td align="center" valign="middle" >80.60</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >37.15</td><td align="center" valign="middle" >60.90</td><td align="center" valign="middle" >86.58</td><td align="center" valign="middle" >106.12</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >37.50</td><td align="center" valign="middle" >67.05</td><td align="center" valign="middle" >95.88</td><td align="center" valign="middle" >121.87</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >6</td><td align="center" valign="middle" >47.62</td><td align="center" valign="middle" >81.22</td><td align="center" valign="middle" >113.40</td><td align="center" valign="middle" >145.17</td></tr><tr><td align="center" valign="middle" >Wall Normal R 3</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >10.01</td><td align="center" valign="middle" >14.09</td><td align="center" valign="middle" >17.70</td><td align="center" valign="middle" >21.40</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >16.45</td><td align="center" valign="middle" >28.49</td><td align="center" valign="middle" >36.74</td><td align="center" valign="middle" >46.61</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >24.14</td><td align="center" valign="middle" >39.44</td><td align="center" valign="middle" >53.67</td><td align="center" valign="middle" >69.85</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >28.82</td><td align="center" valign="middle" >51.80</td><td align="center" valign="middle" >74.75</td><td align="center" valign="middle" >95.94</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >32.76</td><td align="center" valign="middle" >60.15</td><td align="center" valign="middle" >87.46</td><td align="center" valign="middle" >112.59</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >6</td><td align="center" valign="middle" >40.07</td><td align="center" valign="middle" >73.52</td><td align="center" valign="middle" >104.20</td><td align="center" valign="middle" >138.55</td></tr></tbody></table></table-wrap><p>From examination of Figures 11-13 one can conclude the following:</p><p>・ Model 1 does not account for any of the measured reactions.</p><p>・ Models 2 and 3 account more or less equally well for the normal reaction R 1 at the ground support (<xref ref-type="fig" rid="fig1">Figure 1</xref>1).</p><p>・ Model 2 does not account for the measurements of the parallel reaction R 2 at the ground (<xref ref-type="fig" rid="fig1">Figure 1</xref>2) or normal reaction R 3 at the wall (<xref ref-type="fig" rid="fig1">Figure 1</xref>3).</p><p>・ Model 3 accounts satisfactorily for all three reaction forces for μ ≈ 0.25 estimated from the data by trial and inspection.</p><p>To test the consistency of Model 3, the reaction parallel to the wall R 4 was obtained from the data by the two independent relations</p><p>( a )       R 4 = μ R 3 ( b )       R 1 + R 4 = W + P (42)</p><p>where relation (a) is one of the defining conditions of the model, and relation (b) is required by Newton’s 2nd law applied to static equilibrium. Panel A of <xref ref-type="fig" rid="fig1">Figure 1</xref>4</p>Interpretation of Results<p>From the foregoing comprehensive experimental test, one can conclude that, of the various mechanical conditions presumed to describe the static equilibrium of a ladder against a wall, those upon which Model 3 is based best represent the actual physical situation. It is useful to summarize explicitly the nature of these conditions.</p><p>・ The elimination of Model 1 indicates that the length of the ladder is not constrained by axial forces at the ground or wall. Recall that it was this constraint that provided a complementary relation for resolving static indeterminateness in the axial direction. However, the reported experiment is consistent with no unrelieved axial strain. This means that the ladder, in response to its own weight and/or to the applied load, has relaxed to its (very slightly) compressed length in the course of reaching static equilibrium.</p><p>・ The elimination of Model 2 (and Model 1) likewise indicates that the bending of the ladder does not engender a torque at the ground support. Recall that it was this initially unknown moment that necessitated a complementary relation to resolve the static indeterminateness in the transverse direction. The reported experiment is consistent with no internal rotational moment at the base. This means that the ladder, initially presumed to be straight, has relaxed (very slightly) to the shape of an elastic deformation curve like those shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>0 in the course of reaching static equilibrium.</p><p>・ The satisfactory accord between Model 3 and the experimental results tends to confirm that, for a ladder stationary, but not clamped, at its base, it is the friction at the wall that provides the complementary relation needed to resolve static indeterminateness so that all four forces of reaction can be calculated. The theoretical expressions for the reactions do not depend on the elastic modulus of the ladder. They do depend, however, on the coefficient of friction of the ladder at the wall. The experimental results are consistent with an empirical relation whereby the friction at a surface is proportional to the normal force on that surface.</p><p>The theoretical models tested in this paper each embodied one or more distinct physical attributes conceivably responsible for the forces of reaction on a fixed ladder modeled by a single beam. In this regard, the models were never intended to provide an exact theory of the stresses and strains within an actual ladder. To achieve that would require a detailed mechanical model, outside the scope of this paper, taking account of the exact shape, weight distribution, and area moments of inertia of the components of the ladder, including rungs. This caveat must be borne in mind when comparing theory and experiment in Figures 11-14. The dashed lines signifying Model 3 in the figures do not pass through all the experimental points, although they are sufficiently close for the purposes of distinguishing among the models. Small adjustments in the value of the coefficient of friction μ could achieve better agreement in one reaction or another, but there is little reason to make such refinements. Nor is there any reason to attempt to obtain an optimum value of μ by a statistical fitting process such as the method of maximum likelihood, or least square’s analysis, or by Bayes’ theorem [<xref ref-type="bibr" rid="scirp.87003-ref18">18</xref>] . A nonlinear fit of that kind would be computationally intensive and, assuming it even converged, would yield values of μ not much better than the estimate obtained simply and quickly by trial and inspection.</p><p>The great practical utility of Model 3, now that it is validated by experiment, is that it provides a simple, but satisfactory, way to estimate the forces exerted by ladders on structures as a function of the climber’s weight and position of ascent and the ladder’s angle of inclination. Moreover, it is readily generalizable in special situations of practical interest that may arise in engineering design or matters of litigation, such as cases where the ground is flat but sloped, or the wall is flat but not vertical, or the ladder inclines against a railing, not a wall, and overhangs it.</p></sec><sec id="s4"><title>4. Effects of Impact Loads</title><p>Although the primary focus of this paper has been on the static equilibrium of a ladder against a wall, the manner in which ladders are used in practice makes it necessary to discuss, albeit briefly, the question of impact loads. The experimental data of <xref ref-type="table" rid="table1">Table 1</xref> were obtained by carefully and sequentially placing calibrated masses at designated rungs along the ladder so as to avoid imparting momentum to the ladder. <xref ref-type="fig" rid="fig1">Figure 1</xref>5 shows the normal and parallel reactions (in arbitrary units) at the wall recorded as a function of time as a climber slowly ascends and then slowly descends the ladder, resting for a few seconds at each rung. The horizontal red bars mark the approximate values of the equilibrium reactions when the climber is stationary at each rung.</p><p>One sees, however, significant fluctuations about the equilibrium values as the climber transitions from one rung to the next. In particular, on the ascent, the reactions at the wall drop suddenly as the climber lifts his foot from a rung, and then rises suddenly as he places his foot on the next higher rung. On the descent, the reverse occurs for the normal reaction at the wall: the reaction increases dramatically as the climber raises his foot to step on the rung below.</p><p>The conveyance of momentum to a structure over a time interval is referred to as impact or impulse loading [<xref ref-type="bibr" rid="scirp.87003-ref19">19</xref>] . According to Newton’s laws of motion, when the climber accelerates upward, the ladder must recoil downward, thereby exerting an additional force on whatever the ladder is leaning against. In physics, the integral of a force F ( t ) over time t is designated impulse. An essential</p><p>matter in the safe use of a ladder is to determine how large a force the impulsive loading by a climber creates in comparison to the forces of static equilibrium.</p><p>One way to do this is to model the ladder as linear elastic structure governed by Hooke’s law</p><p>P ( t ) = − k ​ y ( t ) (43)</p><p>in which k is the force constant and y ( t ) is the deformation under the load P ( t ) . It is well known from the physics of oscillators [<xref ref-type="bibr" rid="scirp.87003-ref20">20</xref>] that the force constant k and the oscillation angular frequency ω are related by</p><p>k = M ω 2 (44)</p><p>where M is the mass of the object. Newton’s second law, in the absence of dissipation of mechanical energy, then takes the form</p><p>d 2 y d t 2 + ω 2 y = P ( t ) M (45)</p><p>for which the general solution, derivable by means of Green’s functions [<xref ref-type="bibr" rid="scirp.87003-ref20">20</xref>] , is</p><p>y ( t ) = 1 M ω ∫ 0 t sin ω ( t − t ′ ) P ( t ′ ) d t ′ (46)</p><p>with initial conditions</p><p>y ( 0 ) = d y d t | t = 0 = 0 . (47)</p><p>The impulse load P ( t ) can be modeled as a pulse of amplitude P 0 and width τ by the difference of two Heaviside functions</p><p>P ( t ) = P 0 [ H ( t ) − H ( t − τ ) ] . (48)</p><p>Substitution of Equation (48) into (46) leads to the solution</p><p>y ( t ) = P 0 M ω 2 [ cos ω ( t − τ ) − cos ω t ] (49)</p><p>from which one obtains by standard calculus the maximum deflection</p><p>y max = P 0 M ω 2 2 ( 1 − cos ω τ ) (50)</p><p>and therefore from Equations (43), (44) and (50) the maximum force of impact</p><p>P max = | M ω 2 y max | ω ​ τ = π = 2 P 0 . (51)</p><p>P 0 is interpretable as the time-averaged load or, equivalently, the load at static equilibrium―i.e. the weight of the climber in the context of a ladder.</p><p>To the extent that the conditions of the pulse model realistically apply to a ladder, Equation (51) indicates that a climber can generate a maximum impact load that is twice his weight. <xref ref-type="fig" rid="fig1">Figure 1</xref>6 shows a plot of the normal reaction at the ground (i.e. normal to a horizontal force platform) as a climber rapidly ascends and descends a ladder―in contrast to the slow ascent and descent of <xref ref-type="fig" rid="fig1">Figure 1</xref>5 (note that the time base in <xref ref-type="fig" rid="fig1">Figure 1</xref>6 is 1/3 that of <xref ref-type="fig" rid="fig1">Figure 1</xref>5). The blue horizontal bars mark the approximate reactions at static equilibrium at each rung. Of</p><p>particular interest is the region within the black rectangle, which corresponds to the transition of the climber from rung 3 to rung 4 within a time interval Δ t of about 0.2 s. Since the vertical separation between the rungs is Δ s = 30 cm, the vertical acceleration a of the climber can be estimated from the kinematics of uniform acceleration [<xref ref-type="bibr" rid="scirp.87003-ref21">21</xref>] to be</p><p>a = 2 Δ s ( Δ t ) 2 ≈ 2 &#215; 0.3   m ( 0.2   s ) 2 = 15   m ⋅ s − 2 (52)</p><p>or, relative to the acceleration of gravity g = 9.81   m ⋅ s − 2 ,</p><p>a g = P climber P 0 ≈ 15 10 = 1.5 . (53)</p><p>Equation (53) is consistent with Equation (51) of the dynamical pulse model. Indeed, in a more comprehensive model, the dissipation of mechanical energy must be taken into account and would lead to a ratio ( P max / P 0 ) &lt; 2 . Nevertheless, the import of Equation (51) is that the numerical difference between impact forces and forces at static equilibrium can be consequential for those who use ladders in their work or who are responsible for setting safety guidelines and regulations with respect to ladders.</p></sec><sec id="s5"><title>5. Reaction Forces on a Single Wood Beam Inclined against a Wall</title><p>The single-beam models developed in Section 2 to describe the reaction forces on a ladder inclined against a wall should, in principle, apply to an actual single beam inclined against a wall. The author is unaware of any previously reported experiment to measure these forces and determine which, if any, of the theoretical models best accounts for these forces. This section reports the results of such a test.</p><p>Measurements using one horizontal force platform at the base and one vertically mounted force platform at the top were made on a wood beam of mass 6.92 kg, length L = 2.04 m, width w = 7.5 cm, height h = 8.5 cm (symbols correspond to panel C of <xref ref-type="fig" rid="fig1">Figure 1</xref>) inclined against a wall at an angle of 57.7˚ to the ground. Starting at the center of the beam, round eyelets were screwed into the narrower side at intervals of 30 cm for a total of six eyelets. Loads were applied to the beam in units of 2 kg masses stacked on a 444 g weight hanger hooked to one of the eyelets. The eyelets, which served the equivalent function of rungs of the ladder, permitted load placements at fractional lengths of the beam</p><p>β = ( 42 204 , 72 204 , 102 204 , 132 204 , 162 204 , 192 204 ) ≈ ( 0.206 , 0.353 , 0.500 , 0.647 , 0.794 )</p><p>of which the first five eyelets were used in the experiment. Total loads (apart from the weight of the beam) ranged from about 4.4 N to 102.4 N.</p><p><xref ref-type="table" rid="table2">Table 2</xref> records the measurements of the four forces of reaction. With a single beam, in contrast to a ladder, it was possible to measure directly the parallel force at the wall because the beam contacted the center of the vertical force platform where the shear-sensitive strain gauge was located.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>7 shows a comparison of the measured reactions with the predictions of the three theoretical models for four different loads. One sees again, as in the case of the ladder measurements, that Model 1 (dashdot line) does not fit any of the data. For a coefficient of friction at the wall μ = 0.35 , Models 2 and 3 give close predictions for the normal force at the ground (blue circles and curves) and the parallel force at the wall (maroon circles and curves), although the curve due to Model 3 (dashed line) passes through the data points (solid circles) and the curve due to Model 2 (solid line) does not. However, the parallel force at the ground (red circles and curves) and the normal force at the wall (green circles and curves) suffice to eliminate Model 2. The linear curve due to Model 3 (dashed line) very closely matches the measured points, whereas the curve due to Model 2 (solid line) has a marked curvature and passes well below the red solid circles. Only Model 3 accounts satisfactorily for all the data for a single value of the coefficient of friction.</p><p>This experiment lends strong support to the conclusion that the reaction forces on a single beam are described by the same theoretical model that describes the reaction forces on a fixed ladder, and that, as long as the base of either does not slip, the complementary relation needed to resolve the static indeterminateness of the problem is the force of friction at the wall.</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Experimental normal and parallel reactions on a single wood beam</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >θ = 57.7 ∘</th><th align="center" valign="middle" >LOAD (N)</th><th align="center" valign="middle" >43.55</th><th align="center" valign="middle" >63.15</th><th align="center" valign="middle" >82.75</th><th align="center" valign="middle" >102.35</th></tr></thead><tr><td align="center" valign="middle" >REACTIONS (N)</td><td align="center" valign="middle" >EYELET</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >Ground Normal R 1</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >100.3</td><td align="center" valign="middle" >118.8</td><td align="center" valign="middle" >137.7</td><td align="center" valign="middle" >156.6</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >99.6</td><td align="center" valign="middle" >118.6</td><td align="center" valign="middle" >137.3</td><td align="center" valign="middle" >154.6</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >100.0</td><td align="center" valign="middle" >119.2</td><td align="center" valign="middle" >134.1</td><td align="center" valign="middle" >150.9</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >98.0</td><td align="center" valign="middle" >115.8</td><td align="center" valign="middle" >131.8</td><td align="center" valign="middle" >149.0</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >99.0</td><td align="center" valign="middle" >113.8</td><td align="center" valign="middle" >131.9</td><td align="center" valign="middle" >148.3</td></tr><tr><td align="center" valign="middle" >Ground Parallel R 2</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >21.8</td><td align="center" valign="middle" >23.1</td><td align="center" valign="middle" >15.9</td><td align="center" valign="middle" >29.5</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >25.8</td><td align="center" valign="middle" >31.0</td><td align="center" valign="middle" >35.3</td><td align="center" valign="middle" >38.5</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >31.9</td><td align="center" valign="middle" >36.8</td><td align="center" valign="middle" >46.9</td><td align="center" valign="middle" >49.4</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >34.6</td><td align="center" valign="middle" >40.4</td><td align="center" valign="middle" >45.6</td><td align="center" valign="middle" >53.1</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >36.2</td><td align="center" valign="middle" >44.3</td><td align="center" valign="middle" >53.6</td><td align="center" valign="middle" >62.6</td></tr><tr><td align="center" valign="middle" >Wall Normal R 3</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >22.7</td><td align="center" valign="middle" >25.0</td><td align="center" valign="middle" >28.2</td><td align="center" valign="middle" >31.1</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >27.5</td><td align="center" valign="middle" >32.8</td><td align="center" valign="middle" >38.0</td><td align="center" valign="middle" >42.1</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >33.3</td><td align="center" valign="middle" >40.1</td><td align="center" valign="middle" >45.7</td><td align="center" valign="middle" >51.5</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >36.7</td><td align="center" valign="middle" >45.8</td><td align="center" valign="middle" >53.5</td><td align="center" valign="middle" >62.2</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >42.4</td><td align="center" valign="middle" >44.3</td><td align="center" valign="middle" >63.0</td><td align="center" valign="middle" >73.0</td></tr><tr><td align="center" valign="middle" >Wall Parallel R 4</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >9.0</td><td align="center" valign="middle" >9.7</td><td align="center" valign="middle" >10.6</td><td align="center" valign="middle" >10.9</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >9.2</td><td align="center" valign="middle" >10.3</td><td align="center" valign="middle" >10.7</td><td align="center" valign="middle" >13.0</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >9.1</td><td align="center" valign="middle" >9.7</td><td align="center" valign="middle" >13.5</td><td align="center" valign="middle" >16.2</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >10.7</td><td align="center" valign="middle" >12.8</td><td align="center" valign="middle" >16.1</td><td align="center" valign="middle" >18.4</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >10.4</td><td align="center" valign="middle" >14.3</td><td align="center" valign="middle" >15.9</td><td align="center" valign="middle" >18.7</td></tr></tbody></table></table-wrap></sec><sec id="s6"><title>6. Conclusions</title><p>The question of how to calculate the four forces of reaction on a ladder in static equilibrium against a wall has been an unresolved issue for many years. The core of the difficulty is that the problem is statically indeterminate and requires complementary relations beyond the equations of static equilibrium. Three fundamental types of models, differing in the assumptions with regard to supports at the ground and wall, were examined. The comprehensive theoretical analysis and experimental measurements of this paper showed unambiguously that the forces of reaction were best accounted for by the model employing a single complementary relation characterizing friction at the wall.</p><p>The force law for friction that resolved the static indeterminateness in the third model is a single-parameter empirical relation. This parameter, the coefficient</p><p>of friction μ , was estimated from, and tightly constrained by, the experimental data. A critical test of the third model was, in fact, to show that a single value of this coefficient sufficed to yield the normal and parallel reactions at the ground and at the wall for all values and locations of the loads.</p><p>As part of the analyses reported here, the axial and transverse deformations of the single-beam ladder were calculated by means of Castigliano’s theorem, and then confirmed by a full calculation of the deflection curve (elastica) of the beam. The analyses demonstrated that the deformation due to flexure is about two to three orders of magnitude larger than the deformation due to compression. Thus, any model that takes into account only axial compression, but neglects flexure, to resolve the static indeterminateness of reactions on an inclined beam is likely to be inaccurate.</p><p>Measurements were also made of the impact loading on a ladder by the ascent and descent of a climber. In a typical ascent without hesitation, the impact loading</p><p>was found to be about 1 1 2 times the climber’s weight. This result is consistent</p><p>with the predicted maximum value of 2 for the ratio of impact load to static load by a dynamic model of the ladder as a linear elastic medium subject to a pulse perturbation.</p><p>The successful description of the reaction forces on a ladder modeled as a single beam prompted the question of whether any of the three theoretical models satisfactorily accounted for the reaction forces on an actual single beam. Measurements analogous to those made on the ladder were also made on a wood beam of comparable length. The results showed that the same model that accounted for the reaction forces on the ladder also best predicted the reaction forces on the single wood beam. Whether this result is completely general, or whether a different outcome might arise for a beam of different composition, length, and/or cross section can only be determined by further experimentation.</p><p>A final point of importance concerns the inclusiveness of the three fundamental models. In other words, is there another independent way to analyze an Euler-Bernoulli (E-B) beam ladder that might produce better agreement with experiment than Model 3 of this paper? In the author’s opinion, the answer is probably “no”. The reason is that other models within the single E-B beam framework are either irrelevant, equivalent, or inconsistent.</p><p>For example, two models not analyzed in this paper include a ladder with simple (i.e. frictionless) support at the wall and a ladder with fixed support (i.e. infinite friction) at the wall. Neither model is relevant because the contact at the wall is neither frictionless nor immovable. (The models of this paper all assume non-slippage at the ground, for otherwise there would be no static equilibrium.)</p><p>Consider, however, a variation of Model 1. Model 1 resolved static indeterminateness by applying two boundary conditions to the wall support in the configuration of Panel B of <xref ref-type="fig" rid="fig1">Figure 1</xref>: no displacement axially (in the direction of force A 2 ) or transversely (in the direction of force B 2 ). Although this model was in marked disagreement with the experiment, a question might be raised as to whether the outcome could be different if the boundary conditions were applied to the wall support in the configuration of Panel A of <xref ref-type="fig" rid="fig1">Figure 1</xref>, which shows the usual orientation of a ladder. The modified boundary conditions would then be: no displacement vertically (in the direction of force R 4 ) or horizontally (in the direction of force R 3 ). To implement these boundary conditions, one would express the strain energy U in terms of R 3 and R 4 (by transforming A 2 and B 2 ), then set to zero the derivatives of U with respect to R 3 and R 4 . This procedure, however, leads to exactly the same reactions as obtained in the original Model 1.</p><p>Another conceivable model might be to express U in terms of R 3 and R 4 but impose only one boundary condition at the wall, e.g. no displacement vertically. The resulting equation, which couples the vertical and horizontal reactions at the wall, is then used together with Equation (31) which expresses zero torque in static equilibrium. This procedure is inconsistent, however, because the first equation leads to a rotational moment at the ground support, whereas the second equation describes the absence of such a moment. The resulting set of reactions depends on the length (L) and width (h) of the beam (in conflict with experimental results) and leads to large reaction forces that violate physical laws.</p><p>In conclusion, for any model of a ladder (or other physical structure) it is paramount that the equations of static equilibrium and the complementary relations be consistent and compatible with the boundary conditions.</p></sec><sec id="s7"><title>Acknowledgements</title><p>E. Donovan, W. Emery, J. Hallquist, W. Scully, L. Sewanan, and W. Strange each helped take data at some point during the initial phase of this project. The author also thanks Trinity College for partial support through the research fund associated with the George A. Jarvis Chair of Physics.</p></sec><sec id="s8"><title>Conflicts of Interest</title><p>The author declares no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s9"><title>Cite this paper</title><p>Silverman, M.P. (2018) Reaction Forces on a Fixed Ladder in Static Equilibrium: Analysis and Definitive Experimental Test of the Ladder Problem. World Journal of Mechanics, 8, 311-342. https://doi.org/10.4236/wjm.2018.89024</p></sec><sec id="s10"><title>Appendix: Integration of the Heaviside Function</title><p>The Heaviside function defined in Equation (12) and repeated below</p><p>H ( x − x 0 ) = { 0       x &lt; x 0 1           x ≥ x 0 (54)</p><p>provides an economic way to represent piecewise continuous functions, such as frequently occur in the mechanics of continuous media. In the representation of static forces, impact forces, rotational moments, and strain energy densities it is often necessary to integrate over expressions of the form z n H ( z − a ) where a is a non-negative constant and n is a non-negative integer. The general result is</p><p>∫ 0 x z n H ( z − a ) d z = ( x n + 1 − a n + 1 n + 1 ) H ( x − a )                 ( n = 0 , 1 , 2 , ⋯ ) . (55)</p><p>A general proof, which will not be given here, can be constructed inductively from the examples below for n = 0 , 1 .</p><p>1) Case n = 0 :</p><p>∫ 0 x H ( z − a ) d z = ∫ 0 x H ( z − a ) d ( z − a ) = ∫ − a x − a H ( u ) d u = ( x − a ) H ( x − a ) (56)</p><p>Steps: Transform the integration variable from z to u = z − a and use the properties of H ( u ) in Equation (54).</p><p>2) Case n = 1 :</p><p>∫ 0 x z H ( z − a ) d z = ∫ 0 x ( z − a ) H ( z − a ) d ( z − a ) + a ∫ 0 x H ( z − a ) d ( z − a ) = ∫ − a x − a u H ( u ) d u + a ∫ − a x − a H ( u ) d u = ( x − a ) 2 2 H ( x − a ) + a [ ( x − a ) H ( x − a ) ] = x 2 − a 2 2 H ( x − a ) (57)</p><p>Steps:</p><p>・ Line 1: Transform the integrand into a function of u = z − a by appropriate subtraction and addition of terms.</p><p>・ Line 2: Perform the first integral and import the result obtained in the case n = 0 for the second integral.</p><p>・ Lines 3 and 4: Expand and simplify the algebraic expressions to obtain the form of Equation (55).</p></sec><sec id="s11"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.87003-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Ladders 101. https://www.americanladderinstitute.org/page/Ladders101</mixed-citation></ref><ref id="scirp.87003-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Silverman, M.P. (2009) Final Report: Force of a Ladder on a Railing under Static and Dynamic Conditions. Case of Shattuck v. Wynfield, Hartford Superior Court, Hartford, CT.</mixed-citation></ref><ref id="scirp.87003-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Knight, R.D. (2017) Physics for Scientists and Engineers: A Strategic Approach. 4th Edition, Pearson Education Inc., Boston, MA, 113.</mixed-citation></ref><ref id="scirp.87003-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Mendelson, K.S. (1994) Statics of a Ladder Leaning against a Rough Wall. American Journal of Physics, 63, 148-150. https://doi.org/10.1119/1.17972</mixed-citation></ref><ref id="scirp.87003-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Gonzalez, A.G. and Gratton, J. (1996) Reaction Forces on a Ladder Leaning against a Rough Wall. American Journal of Physics, 64, 1001-1005. https://doi.org/10.1119/1.18317</mixed-citation></ref><ref id="scirp.87003-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Hakkinen, K.K., Pesonen, J. and Rajamaki, E. (1988) Experiments on Safety in the Use of Portable Ladders. Journal of Occupational Accidents, 10, 1-19. https://doi.org/10.1016/0376-6349(88)90002-8</mixed-citation></ref><ref id="scirp.87003-ref7"><label>7</label><mixed-citation publication-type="book" xlink:type="simple">Lee, Y.H. and Tung, E.K. (1992) Body and Ladder Mechanical Stresses Analysis in a Climbing Strike. In: Kumar, S., Ed., Advances in Industrial Ergonomics and Safety IV, Taylor and French, Montreal, 1007-1014.</mixed-citation></ref><ref id="scirp.87003-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Chang, W.-R., Chang, C.-C., Matz, S. and Son, D.H. (2004) Friction Requirements for Different Climbing Conditions in Straight Ladder Ascending. Safety Science, 42, 791-805. https://doi.org/10.1016/j.ssci.2004.02.002</mixed-citation></ref><ref id="scirp.87003-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Silverman, M.P. and Farrah, J. (2018) Bending of a Tapered Rod: Modern Application and Experimental Test of Elastica Theory. World Journal of Mechanics, 8, 272-300. https://doi.org/10.4236/wjm.2018.87022</mixed-citation></ref><ref id="scirp.87003-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Hibbeler, R.C. (1997) Mechanics of Materials. 3rd Edition, Prentice Hall, Upper Saddle River, NJ, 582-583.</mixed-citation></ref><ref id="scirp.87003-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Fowles, G.R. and Cassiday, G.L. (1993) Analytical Mechanics. 5th Edition, Harcourt Brace, New York, 54-55.</mixed-citation></ref><ref id="scirp.87003-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Elastic Modulus of Aluminum. https://www.engineeringtoolbox.com/young-modulus-d_417.html</mixed-citation></ref><ref id="scirp.87003-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Material Property Data. http://www.matweb.com/Search/MaterialGroupSearch.aspx?GroupID=316</mixed-citation></ref><ref id="scirp.87003-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Elastic Modulus of Wood. https://www.engineersedge.com/lumber.htm</mixed-citation></ref><ref id="scirp.87003-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Budynas, R.G. (1977) Advanced Strength and Applied Stress Analysis. McGraw-Hill, New York, 240-247.</mixed-citation></ref><ref id="scirp.87003-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Matthews, J. and Walker, R.L. (1973) Mathematical Methods of Physics. W.A. Benjamin, Menlo Park, 107-109.</mixed-citation></ref><ref id="scirp.87003-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Levien, R. (2008) The Elastica: A Mathematical History. Technical Report No. UCB/EECS 2008-103, Electrical Engineering and Computer Sciences University of California, Berkeley. http://www.eecs.berkeley.edu/Pubs/TechRpts/2008/EECS-2008-103.html</mixed-citation></ref><ref id="scirp.87003-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Silverman, M.P. (2014) A Certain Uncertainty: Nature’s Random Ways. Cambridge University Press, Cambridge, 45-83. https://doi.org/10.1017/CBO9781139507370</mixed-citation></ref><ref id="scirp.87003-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Paz, M. (1985) Structural Dynamics: Theory &amp; Computation. Van Nostrand Reinhold Co., New York, 63-65.</mixed-citation></ref><ref id="scirp.87003-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Marion, J.B. and Thornton, S.T. (1988) Classical Dynamics of Particles and Systems. 3rd Edition, Harcourt Brace Jovanovich, New York, 98-135.</mixed-citation></ref><ref id="scirp.87003-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Symon, K.R. (1971) Mechanics. 3rd Edition, Addison-Wesley, Reading, 11-13.</mixed-citation></ref></ref-list></back></article>