buckling_spring


buckling_spring, an Octave code which illustrates solutions of the buckling spring equations.

We are given three points A, B, and C:

Springs:

Forces:

  • A vertical load MU is applied at point B (downward is positive);
  • A horizontal load LAMBDA is applied at point C (leftware is positive);
  • The spring force is applied perpendicularly to the axes of the two springs.
  • If we compute F(1), the force along the axis of one spring, and F(2), the force perpendicular to the axis of one spring, we have that F(L,THETA,LAMBDA,MU) is given by:

    To explore these equations, we can first solve for MU and LAMBDA in terms of L and THETA:

    We can then study the behavior of solutions by choosing a fixed value of THETA (say pi/8), and plotting LAMBDA(L,THETA) versus MU(L,THETA) over a range of values of L, say from 0.25 to 1.75. Recall that L = 1 when the springs are at their "natural" length. This approach is taken by the function BUCKLING_SPRING_L.

    We can also plot the same data, but instead fix a value of L, typically in the range of 0.25 to 1.75, and the plot LAMBDA(L,THETA) versus MU(L,THETA) over a range of values of THETA, say from -3pi/8 to +3pi/8. This approach is taken by the function BUCKLING_SPRING_THETA.

    Usage:

    To draw lines of constant L, with THETA varying along a line:

    buckling_spring_l ( l_num, theta_num )
    where

    To draw lines of constant THETA, with L varying along a line:

    buckling_spring_theta ( l_num, theta_num )
    where

    Licensing:

    The computer code and data files described and made available on this web page are distributed under the MIT license

    Languages:

    buckling_spring is available in a Mathematica version and a MATLAB version and an Octave version.

    Related Data and Programs:

    buckling_spring_test

    Reference:

    1. Tim Poston, Ian Stewart,
      Catastrophe Theory and its Applications,
      Dover, 1996,
      ISBN13: 978-0486692715,
      LC: QA614.58.P66.

    Source Code:


    Last revised on 12 June 2023.