6 Engineering Example 4

6.1 Deflection of a uniformly loaded beam

Introduction

A uniformly loaded beam of length L is supported at both ends as shown in Figure 9. The deflection y ( x ) is a function of horizontal position x and obeys the ordinary differential equation (ODE)

d 4 y d x 4 ( x ) = 1 E I q ( x ) (1)

where E is Young’s modulus, I is the moment of inertia and q ( x ) is the load per unit length at point x . We assume in this problem that q ( x ) = q a constant. The boundary conditions are (i) no deflection at x = 0 and x = L (ii) no curvature of the beam at x = 0 and x = L .

Figure 9 :

{ The bending beam, parameters involved in the mathematical formulation}

Problem in words

Find the deflection of a beam, supported so that that there is no deflection and no curvature of the beam at its ends, subject to a uniformly distributed load, as a function of position along the beam.

Mathematical statement of problem

Find the equation of the curve y ( x ) assumed by the bending beam that satisfies the ODE (1). Use the coordinate system shown in Figure 9 where the origin is at the left extremity of the beam. In this coordinate system, the boundary conditions, which require that there is no deflection at x = 0 and x = L , and that there is no curvature of the beam at x = 0 and x = L , are

  1. y ( 0 ) = 0
  2. y ( L ) = 0
  3. d 2 y d x 2 x = 0 = 0
  4. d 2 y d x 2 x = L = 0
  5. d 2 y d x 2 x = L = 0

Note that d y ( x ) d x and d 2 y ( x ) d x 2 are respectively the slope and the radius of curvature of the curve at point ( x , y ) .

Mathematical analysis

Integrating Equation (1) leads to:

E I d 3 y d x 3 ( x ) = q x + A (2)

Integrating a second time:

E I d 2 y d x 2 ( x ) = q x 2 2 + A x + B (3)

Integrating a third time:

E I d y d x ( x ) = q x 3 6 + A x 2 2 + B x + C (4)

Integrating a fourth time:

E I y ( x ) = q x 4 24 + A x 3 6 + B x 2 2 + C x + D . (5)

The boundary conditions 1. and 2. enable determination of the constants of integration A , B , C , D . Indeed, the boundary condition 1., y ( 0 ) = 0 , and Equation (5) give

E I y ( 0 ) = q × ( 0 ) 4 24 + A × ( 0 ) 3 6 + B × ( 0 ) 2 2 + C × ( 0 ) + D = 0

which yields D = 0 .

The boundary condition 2., y ( L ) = 0 , and Equation (5) give

E I y ( L ) = q L 4 24 + A L 3 6 + B L 2 2 + C L + D .

Using the newly found value for D one writes

q L 4 24 + A L 3 6 + B L 2 2 + C L = 0 (6)

The boundary condition 3. obtained from the definition of the radius of curvature, d 2 y d x 2 ( 0 ) = 0 , and Equation (3) give

I d 2 y d x 2 ( 0 ) = q × ( 0 ) 2 2 + A × ( 0 ) + B

which yields B = 0 . The boundary condition 4., d 2 y d x 2 ( L ) = 0 , and Equation (3) give

E I d 2 y d x 2 ( L ) = q L 2 2 + A L = 0

which yields A = q L 2 . The expressions for A , B , D are introduced in Equation (6) to find the last unknown constant C . This leads to q L 4 24 q L 4 12 + C L = 0 or C = q L 3 24. Finally, Equation (5) and the values of constants lead to the solution

y ( x ) = [ q x 4 24 q L x 3 12 + q L 3 x 24 ] E I . (7)

Interpretation

The predicted deflection is zero at both ends as required, and you may check that it is symmetrical about the centre of the beam by switching to the coordinate system ( X , Y ) with L 2 x = X and y = Y and verifying that the deflection Y ( X ) is symmetrical about the vertical axis, i.e. Y ( X ) = Y ( X ) .

Exercises
  1. In an RC circuit (a resistor and a capacitor in series) the applied emf is a constant E . Given that d q d t = i where q is the charge in the capacitor, i the current in the circuit, R the resistance

    and C the capacitance the equation for the circuit is

    R i + q C = E .

    If the initial charge is zero find the charge subsequently.

  2. If the voltage in the RC circuit is E = E 0 cos ω t find the charge and the current at time t .
  3. An object is projected from the Earth’s surface. What is the least velocity (the escape velocity) of projection in order to escape the gravitational field, ignoring air resistance.

    The equation of motion is

    m v d v d x = m g R 2 x 2

    where the mass of the object is m , its distance from the centre of the Earth is x and the radius of the Earth is R .

  4. The radial stress p at distance r from the axis of a thick cylinder subjected to internal pressure is given by p + r d p d r = A p where A is a constant. If p = p 0 at the inner wall ( r = r 1 ) and is negligible ( p = 0 ) at the outer wall ( r = r 2 ) find an expression for p .
  5. The equation for an LCR circuit with applied voltage E is

    L d i d t + R i + 1 C q = E .

    By differentiating this equation find the solution for q ( t ) and i ( t ) if L = 1 , R = 100 , C = 1 0 4 and E = 1000 given that q = 0 and i = 0 at t = 0 .

  6. Consider the free vibration problem in Section 19.4 subsection 2 (page 57) when m = 1 , n = 1 and k = 2 (critical damping).

    Find the solution for x ( t ) .

  7. Repeat Exercise 6 for the case    m = 1 , n = 1 and k = 1.5  (light damping)
  8. Consider the forced vibration problem in Section 19.4 subsection 2 with m = 1 , n = 25 , k = 8 , E = sin 3 t , x 0 = 0 with an initial velocity of 3.
  9. This refers to the Task on page 55 concerning modelling the dissolving of a pill in the stomach.

    An alternative model supposes that the pill is very rapidly permeated by stomach acids and the small granules contained in the capsule dissolve individually. In this case, the rate of change of volume is assumed to be directly proportional to the volume. Using the experimental data given in the Task, estimate the time for 95% of the pill to dissolve, based on this alternative model, and compare results.

  1. Use the equation in the form R d q d t + q C = E or d q d t + 1 R C q = E R .

    The integrating factor is e t R C and the general solution is

    q = E C ( 1 e t R C ) and as t q E C .

  2. q = E 0 C 1 + ω 2 R 2 C 2 cos ω t e t R C + ω R C sin ω t

    i = d q d t = E 0 C 1 + ω 2 R 2 C 2 ω sin ω t + 1 R C e t R C + ω 2 R C cos ω t .

  3. v min = 2 g R . If R = 6378 km and g = 9.81 m s 2 then v min = 11.2 km s 1 .
  4. p = p 0 r 1 2 r 1 2 r 2 2 1 r 2 2 r 2
  5. q = 0.1 1 10 3 e 50 t ( sin 50 3 t + 3 cos 50 3 t ) i = 20 3 e 50 t sin 50 3 t .
  6. x = x 0 ( 1 + t ) e t
  7. x = x 0 e 0.75 t ( cos 7 4 t + 3 7 sin 7 4 t )
  8. x = 1 104 e 4 t 3 cos 3 t + 106 sin 3 t 3 cos 3 t + 2 sin 3 t
  9. This leads to d V d t = k V and V = V 0 e k t where k = 1 3 ln 2 . The time taken is about 4 hr 19 min. This is much less than the other model, as should be expected.