Singularity Functions

Tim Lahey http://www.cgl.uwaterloo.ca/tjlahey/ February 28, 2000

Description of Singularity Functions

To handle the discontinuities in V ( x) and M( x) curves we introduce a family of functions called singularity functions. We dene the function as: ( x a)n x a, f ( x) < x a > n = 0 x < a.

(1)

The basic singularity functions are shown in Figure 1. They are (from top to bottom): unit doublet, unit impulse, unit step, unit ramp, and unit acceleration. The unit impulse is sometimes referred to as the Dirac delta function. The integration rule for singularity functions is: < x a >1 Z x < x a >n dx = < x a >0 < xa>n+1
n+1

n = 2, n = 1, n 0. (2)

The loading of beams can be determined from a superposition of singularity functions for the load distribution function q( x). The unit doublet is the distribution function representation for the applied moment and the unit impulse is the representation for an applied load. For example, an applied torque of M0 at point x = a is q( x) = M0 < x a >2 . Once we have the

<x - a>-2

a
-1 <x - a>

<x - a>0

<x - a>1

<x - a>2

Figure 1: Basic Singularity Functions 2

distribution function q( x), we can integrate to get the shear V ( x) and the moment M( x) functions. V ( x) = M( x) =
Z Z

q( x) + C1

(3) (4)

V ( x) + C1 < x >0 +C2

Integration of the singularity functions in the load distribution and the resulting internal shear force distribution and internal bending moment distribution is summarised in Table 1.

Loading M0

Distribution q(x)

Shear = q( x)dx
R

Moment = q( x)dx
R

x
Table 1: Summary of Singularity Functions

a P

q = M0 < x a >2

V = M0 < x a >1

M = M0 < x a >0

x w

q = P < x a > 1

V = P < x a >0

M = P < x a >1

a m

q = w < x a > 0

V = w < x a >1

2 M = w 2 < xa >

q = m < x a > 1

V=

m 2

< x a >2

3 M = m 6 < xa >

Example 1

To illustrate the use of singularity functions in getting the shear and moment equations, consider this simple example. Solving for the reactions:

500 lb

200 lb/ft

x
4 6

Free Body Diagram 500 lb 200 lb/ft

M V
Figure 2: Simple Loading Case

VL = 500 + 200(6) = 1700lb. ML = 500(10) 1200(3) = 8600lb. f t. 5

(5) (6)

First, nd q(x): q( x) = 500 < x >1 200 < x 4 >0 +1700 < x 10 >1 +8600 < x 10 >2 (7) To get the shear, we integrate the loading according as shown below: V ( x) =
Z x
0

q( x)dx

(8)

substituting the loading and integrating: V ( x) = 500 < x >0 +200 < x 4 >1 1700 < x 10 >0 8600 < x 10 >1 (9) We can now integrate the shear to get the bending moment equation as shown below: M( x) = M( x) = 500 < x >1
Z x
0

V ( x)dx

(10)

200 < x 4 >2 +1700 < x 10 >1 +8600 < x 10 >0 2 (11) (12)

M( x) = 500 x 100 < x 4 >2 +8600 < x 10 >0

The shear and bending moment equations can be represented in graphical form:

1800

Shear Force Diagram

1600

1400

Net Curve
1200

1000 V(x)

800

1 200<x - 4>

600

400

500<x>0

200

5 Position (ft)

10

Figure 3: Shear Force Diagram for Example 1

Bending Moment Diagram

-1000

2 -100<x - 4>

-2000

-3000

-500x

-4000 M(x) -5000

Net Curve
-6000

-7000

-8000

-9000

5 Position (ft)

10

Figure 4: Bending Moment Diagram for Example 1

Example 2

To compare and contrast the two techniques for getting the shear and moment equations we will look at the following example.

8kN 4kN/m

4kNm B 3m

D 2m

0.5m 1.5m

Free Body Diagram

8kN 4kN/m

4kNm B R1= 10 kN

R 2 = 26 kN

Figure 5: Multiple Loadings

Total Length: L=3+

1 1 + 1 + 2 = 7m 2 2

(13)

3.1 Method of Sections

VAB + R1 4 x = 0 VAB = 4 x 10 M AB 10 x +
Z x
0

(14) (15) (16) (17) (18)

4 xdx = 0

1 M AB = 4 x2 + 10 x 2 M AB = 2 x2 + 10 x Sections from B to C: 3 x < 3.5

VBC = VAB = 4 x 10 M BC 4 10 x +
Z x
0

(19) (20) (21)

4 xdx = 0

M BC = 2 x2 + 10 x + 4 Sections from C to D: 3.5 x < 5

VCD + 10 4 x 8 = 0 VCD = 4 x 2 MCD 4 10 x + 8( x 3.5) +

(22) (23) 4 xdx = 0 (24) (25) (26)

MCD = 4 + 10 x 8 x + 28 2 x2 MCD = 2 x2 + 2 x + 32

10

Sections from D to E: 5 x 7

VDE + 10 + 26 8 4 x = 0 VDE = 4 x 28 M DE 4 10 x + 8( x 3.5) 26( x 5) + M DE = 2 x2 + 28 x 98

(27) (28) 4 xdx = 0 (29) (30)

3.2 Singularity Functions

q( x)dx

(32)

= 10 < x >0 4 < x 3 >1 +8 < x 3.5 >0 26 < x 5 >0 +4 < x >1 +V0 (33)
Since we have no unknown shear forces, V0 is zero, so our shear force equation is: V ( x) = 10 < x >0 4 < x 3 >1 +8 < x 3.5 >0 26 < x 5 >0 +4 < x >1 (34) Plotting the shear force diagram:

11

30

Shear Force Diagram

20

Net Curve
10

4x 8<x - 3.5>0

V(x)

-10

-10<x>0 0 -26<x - 5>

-20

-30

Position (ft)

Figure 6: Shear Force Diagram for Example 2

12

Integrating the shear force equation: M( x) =

V ( x)dx

(35)

4 = 10 < x >1 +4 < x 3 >0 8 < x 3.5 >1 +26 < x 5 >1 < x >2 + M0 2 (36) We have no unknown moments so, M0 is zero. Our bending moment equation becomes: M( x) = 10 < x >1 +4 < x 3 >0 8 < x 3.5 >1 +26 < x 5 >1 2 < x >2 (37) Finally, plotting the bending moment diagram:

13

80

Bending Moment Diagram

60

10<x>1
40

26<x-5>1 Net Curve 4<x-3>0

20

0 M(x) -20

-8<x-3.5>1

-40

-60

-80

-2<x>2

-100

Position (ft)

Figure 7: Bending Moment Diagram for Example 2

14