Rods - United States Naval Academy

MECHANICAL ENGINEERING DEPARTMENT
UNITED STATES NAVAL ACADEMY
EM423 - INTRODUCTION TO MECHANICAL VIBRATIONS
CONTINUOUS SYSTEMS – LONGITUDINAL VIBRATION IN RODS
SYMBOLS
ρ
mass density
AX
Cross-sectional area
u
Longitudinal deflection from equilibrium position
x
Distance along the rod
P
Internal longitudinal (axial) force. Varies with position along the rod.
INTRODUCTION
This theory is applicable to longitudinal vibration in slender rods. With this type of vibration, very small
amplitudes of motion can produce very large forces. Typical problems include damage to thrust bearings
in ship propulsion systems, noise and damage caused by engine support props in helicopters, and
radiated noise problems in submarines.
ASSUMPTIONS
1.
The rod is homogeneous and isotropic.
2.
The material is within the elastic limit, and obeys Hooke's Law.
3.
The rod is thin compared with its length.
4.
Plane sections remain plane.
5.
The lateral deflection caused by changes in the length of the rod is small.
Rods - 1
THEORY
Consider a small element of a rod that is vibrating axially. The element is at position x along the rod,
where we assume that the left end of the rod is at x = 0, with positive x to the right. At a certain instant
in time the element is displaced an amount u from the equilibrium position.
P
P
x
P+
dx
x
u+
dx
u
x
dx
u
We first apply strength of materials relationships to the element.
The element changes length by an amount
πœ•πœ•πœ•πœ•
𝑑𝑑𝑑𝑑
πœ•πœ•πœ•πœ•
strain =
and so:
(change in length) πœ•πœ•πœ•πœ•
=
(original length)
πœ•πœ•πœ•πœ•
From Hooke's Law, E = (stress) / (strain), and (stress) = (force) / (area), so
Hence
Rearrange and differentiate
𝐸𝐸 =
𝑃𝑃
(stress) οΏ½ �𝐴𝐴π‘₯π‘₯ οΏ½
=
πœ•πœ•πœ•πœ•
(strain)
οΏ½ οΏ½
πœ•πœ•πœ•πœ•
𝑃𝑃
πœ•πœ•πœ•πœ•
=
πœ•πœ•πœ•πœ• 𝐴𝐴π‘₯π‘₯ 𝐸𝐸
πœ•πœ• 2 𝑒𝑒
πœ•πœ•πœ•πœ•
= 𝐴𝐴π‘₯π‘₯ 𝐸𝐸 2
πœ•πœ•π‘₯π‘₯
πœ•πœ•πœ•πœ•
We now apply Newton’s 2nd Law (f=ma) to the element.
Resolve forces along the rod.
Rods - 2
�𝑃𝑃 +
Combining the two results yields:
πœ•πœ•πœ•πœ•
𝑑𝑑2 𝑒𝑒
𝛿𝛿𝛿𝛿� βˆ’ 𝑃𝑃 = (ρ𝐴𝐴π‘₯π‘₯ 𝛿𝛿𝛿𝛿) 2
πœ•πœ•πœ•πœ•
𝑑𝑑𝑑𝑑
1 𝑑𝑑2 𝑒𝑒
πœ•πœ• 2 𝑒𝑒
=
πœ•πœ•π‘₯π‘₯ 2 𝑐𝑐 2 𝑑𝑑𝑑𝑑 2
With
𝐸𝐸
𝑐𝑐 = οΏ½ = wave velocity
𝜌𝜌
This is the same wave equation as derived for string vibrations, but with a different wave velocity. Using
the same method to solve the equation, firstly separate the variables.
Hence:
𝑒𝑒(π‘₯π‘₯, 𝑑𝑑) = π‘ˆπ‘ˆ(π‘₯π‘₯). 𝐺𝐺(𝑑𝑑)
πœ”πœ”πœ”πœ”
πœ”πœ”πœ”πœ”
with
π‘ˆπ‘ˆ(π‘₯π‘₯) = �𝐴𝐴1 𝑠𝑠𝑠𝑠𝑠𝑠 οΏ½ οΏ½ + 𝐴𝐴2 𝑐𝑐𝑐𝑐𝑐𝑐 οΏ½ οΏ½οΏ½
𝑐𝑐
𝑐𝑐
and
𝐺𝐺(𝑑𝑑) = 𝑠𝑠𝑠𝑠𝑠𝑠(πœ”πœ”πœ”πœ”)
πœ”πœ”πœ”πœ”
πœ”πœ”πœ”πœ”
𝑒𝑒(π‘₯π‘₯, 𝑑𝑑) = π‘ˆπ‘ˆ(π‘₯π‘₯). 𝐺𝐺(𝑑𝑑) = �𝐴𝐴1 𝑠𝑠𝑠𝑠𝑠𝑠 οΏ½ οΏ½ + 𝐴𝐴2 𝑐𝑐𝑐𝑐𝑐𝑐 οΏ½ οΏ½οΏ½ . 𝑠𝑠𝑠𝑠𝑠𝑠(πœ”πœ”πœ”πœ”)
𝑐𝑐
𝑐𝑐
COMMON BOUNDARY CONDITIONS
For a FREE end the internal stress is zero, since there is nothing for the end of the rod to react against,
and thus no axial force can be generated.
therefore:
For a FIXED end the displacement is zero.
stress = 𝐸𝐸
πœ•πœ•πœ•πœ•
=0
πœ•πœ•πœ•πœ•
πœ•πœ•πœ•πœ•
πœ•πœ•πœ•πœ•
𝑒𝑒 = 0
Rods - 3
NATURAL FREQUENCIES OF A FREE-FREE ROD, length = L
πœ”πœ”πœ”πœ”
πœ”πœ”πœ”πœ”
𝑒𝑒(π‘₯π‘₯, 𝑑𝑑) = π‘ˆπ‘ˆ(π‘₯π‘₯). 𝐺𝐺(𝑑𝑑) = �𝐴𝐴1 𝑠𝑠𝑠𝑠𝑠𝑠 οΏ½ οΏ½ + 𝐴𝐴2 𝑐𝑐𝑐𝑐𝑐𝑐 οΏ½ οΏ½οΏ½ . 𝑠𝑠𝑠𝑠𝑠𝑠(πœ”πœ”πœ”πœ”)
𝑐𝑐
𝑐𝑐
At x = 0,
πœ•πœ•πœ•πœ•
πœ•πœ•πœ•πœ•
= 0 for all time. Therefore
πœ”πœ”πœ”πœ”
πœ”πœ”πœ”πœ”
π‘ˆπ‘ˆ(π‘₯π‘₯) = �𝐴𝐴1 𝑠𝑠𝑠𝑠𝑠𝑠 οΏ½ οΏ½ + 𝐴𝐴2 𝑐𝑐𝑐𝑐𝑐𝑐 οΏ½ οΏ½οΏ½
𝑐𝑐
𝑐𝑐
πœ•πœ•πœ•πœ•
πœ”πœ”
πœ”πœ”πœ”πœ”
πœ”πœ”
πœ”πœ”πœ”πœ”
= �𝐴𝐴1 οΏ½ οΏ½ 𝑐𝑐𝑐𝑐𝑐𝑐 οΏ½ οΏ½ βˆ’ 𝐴𝐴2 οΏ½ οΏ½ 𝑠𝑠𝑠𝑠𝑠𝑠 οΏ½ οΏ½οΏ½ = 0
πœ•πœ•πœ•πœ•
𝑐𝑐
𝑐𝑐
𝑐𝑐
𝑐𝑐
Therefore
At x = L,
πœ•πœ•πœ•πœ•
πœ•πœ•πœ•πœ•
𝐴𝐴1 = 0
= 0 for all time. Therefore
πœ•πœ•πœ•πœ•
πœ”πœ”
πœ”πœ”πœ”πœ”
= οΏ½0 βˆ’ 𝐴𝐴2 οΏ½ οΏ½ 𝑠𝑠𝑠𝑠𝑠𝑠 οΏ½ οΏ½οΏ½ = 0
πœ•πœ•πœ•πœ•
𝑐𝑐
𝑐𝑐
πœ”πœ”πœ”πœ”
οΏ½
𝑐𝑐
This is only true for all time if 𝑠𝑠𝑠𝑠𝑠𝑠 οΏ½
= 0, so
πœ”πœ”π‘›π‘› 𝐿𝐿
= 𝑛𝑛𝑛𝑛
𝑐𝑐
with n = 1, 2, 3, …
Solving for the natural frequency in Hz:
𝑓𝑓𝑛𝑛 =
πœ”πœ”π‘›π‘›
𝑛𝑛 𝐸𝐸
=
οΏ½
2πœ‹πœ‹ 2𝐿𝐿 𝜌𝜌
with n = 1, 2, 3, …
Q: What are the first 3 natural frequencies of a steel rod that is 170 m long? Both ends are unrestrained.
A: From the datasheet:
𝐸𝐸 = 206 kN/mm2 = 206 × 109 N/m2 = 206 × 109 Pa
𝜌𝜌 = 1738 kg/m3
𝑓𝑓𝑛𝑛 = 𝑓𝑓𝑛𝑛 =
206 × 109
πœ”πœ”π‘›π‘›
𝑛𝑛 𝐸𝐸
𝑛𝑛
οΏ½
=
οΏ½ =
= 𝑛𝑛 × 32.02 Hz
2πœ‹πœ‹ 2𝐿𝐿 𝜌𝜌 2 × 170
1738
n =1 gives f1 = 16 Hz
n =2 gives f1 =48 Hz
n =3 gives f1 = 80 Hz
Rods - 4
NATURAL FREQUENCIES OF A FIXED-FREE ROD, length = L
πœ”πœ”πœ”πœ”
πœ”πœ”πœ”πœ”
𝑒𝑒(π‘₯π‘₯, 𝑑𝑑) = π‘ˆπ‘ˆ(π‘₯π‘₯). 𝐺𝐺(𝑑𝑑) = �𝐴𝐴1 𝑠𝑠𝑠𝑠𝑠𝑠 οΏ½ οΏ½ + 𝐴𝐴2 𝑐𝑐𝑐𝑐𝑐𝑐 οΏ½ οΏ½οΏ½ . 𝑠𝑠𝑠𝑠𝑠𝑠(πœ”πœ”πœ”πœ”)
𝑐𝑐
𝑐𝑐
At x = 0, 𝑒𝑒 = 0 for all time. Therefore
πœ”πœ”πœ”πœ”
πœ”πœ”πœ”πœ”
π‘ˆπ‘ˆ(π‘₯π‘₯) = �𝐴𝐴1 𝑠𝑠𝑠𝑠𝑠𝑠 οΏ½ οΏ½ + 𝐴𝐴2 𝑐𝑐𝑐𝑐𝑐𝑐 οΏ½ οΏ½οΏ½ = 0
𝑐𝑐
𝑐𝑐
Therefore
At x = L,
πœ•πœ•πœ•πœ•
πœ•πœ•πœ•πœ•
𝐴𝐴2 = 0
= 0 for all time. Therefore
πœ•πœ•πœ•πœ•
πœ”πœ”
πœ”πœ”πœ”πœ”
= �𝐴𝐴1 οΏ½ οΏ½ 𝑐𝑐𝑐𝑐𝑐𝑐 οΏ½ οΏ½ βˆ’ 0οΏ½ = 0
πœ•πœ•πœ•πœ•
𝑐𝑐
𝑐𝑐
πœ”πœ”πœ”πœ”
οΏ½
𝑐𝑐
This is only true for all time if 𝑐𝑐𝑐𝑐𝑐𝑐 οΏ½
= 0, so
1
πœ”πœ”π‘›π‘› 𝐿𝐿
= �𝑛𝑛 βˆ’ οΏ½ πœ‹πœ‹
2
𝑐𝑐
Hence
𝑓𝑓𝑛𝑛 =
πœ”πœ”π‘›π‘› �𝑛𝑛 βˆ’ 1οΏ½2οΏ½ 𝐸𝐸
=
οΏ½
2πœ‹πœ‹
2𝐿𝐿
𝜌𝜌
with n = 1, 2, 3, …
with n = 1, 2, 3, …
Q: What are the first 3 natural frequencies of a steel rod that is 3 m long? One end is unrestrained, the
other is fixed.
A: From the datasheet:
𝐸𝐸 = 206 kN/mm2 = 206 × 109 N/m2 = 206 × 109 Pa
𝜌𝜌 = 1738 kg/m3
πœ”πœ”π‘›π‘› �𝑛𝑛 βˆ’ 1οΏ½2οΏ½ 𝐸𝐸 �𝑛𝑛 βˆ’ 1οΏ½2οΏ½ 206 × 109
οΏ½
𝑓𝑓𝑛𝑛 =
=
οΏ½ =
= �𝑛𝑛 βˆ’ 1οΏ½2οΏ½ × 1814.5 Hz
2πœ‹πœ‹
2𝐿𝐿
𝜌𝜌
2×3
1738
n =1 gives f1 = 907 Hz
n =2 gives f1 =2722 Hz
n =3 gives f1 = 4536 Hz
Rods - 5
DISCUSSION
1.
The solution for longitudinal vibrations in a rod is similar to flexural vibrations of a taught string,
but with a different wave speed.
2.
This means the same discussion applies to rod vibrations. The only difference is that, for
graphical presentation, longitudinal displacements are usually drawn as if they were transverse.
3.
The wave speed for longitudinal waves travelling in many common metals is almost the same.
Calculate the wave speeds and fundamental natural frequencies for the following materials and
complete the table.
E (kN/mm2)
Density (kg/m3)
Longitudinal
wave speed
(m/s)
Fundamental
natural
frequency for a
3m long rod (Hz)
Steel
206
7843
Aluminum
70
2720
Brass
105
8410
Nickel Alloy
207
8580
Magnesium
45
1738
Rods - 6
EXAMPLES
1.
A steel propulsion shaft on a ship is 25 m long. Its outside and inside diameters are 32 cm and
15 cm respectively. What are the first 2 natural frequencies of longitudinal vibration? Solve the problem
2 times. First, assume the gearbox and propeller are very light. Second, assume the gearbox is light, but
the propeller is very heavy.
L = 25 m; ρ = 7843 kg/m3; E = 206 kN/mm2
Solve for free-free boundary conditions.
𝑓𝑓𝑛𝑛 =
So
𝑓𝑓1 =
𝑓𝑓2 =
𝑛𝑛 𝐸𝐸
οΏ½
2𝐿𝐿 𝜌𝜌
206 × 109
1
οΏ½
= 102.5𝐻𝐻𝐻𝐻
7843
2 × 25
206 × 109
2
οΏ½
= 205.0𝐻𝐻𝐻𝐻
7843
2 × 25
Solve for fixed-free boundary conditions.
𝑓𝑓𝑛𝑛 =
So
𝑓𝑓1 =
�𝑛𝑛 βˆ’ 1οΏ½2οΏ½ 𝐸𝐸
οΏ½
2𝐿𝐿
𝜌𝜌
οΏ½1 βˆ’ 1οΏ½2οΏ½ 206 × 109
οΏ½
= 51.2𝐻𝐻𝐻𝐻
7843
2 × 25
οΏ½2 βˆ’ 1οΏ½2οΏ½ 206 × 109
οΏ½
𝑓𝑓2 =
= 153.7𝐻𝐻𝐻𝐻
7843
2 × 25
2.
What difference does it make to the previous question if you assume both the propeller and
gearbox are very heavy?
The solution is identical to that for a free-free rod. Why?
Rods - 7
ASSIGNMENTS
1.
Find, from first principles, the velocity of longitudinal waves along a thin steel rod.
2.
The figure shows a schematic of a magnetostriction
oscillator. The frequency of the oscillator is determined by
L
the length of the nickel alloy rod, which generates an
alternating voltage in the surrounding coils equal to the
frequency of the fundamental longitudinal vibration of the
rod. Determine the full length, L, of the rod for a frequency
of 20 kHz.
3.
(extra credit) A uniform rod of length L and cross-sectional area Ax is fixed at the upper end and
is loaded with a weight W=mg on the other end. Show that the natural frequencies of longitudinal
vibration are determined from the equation:
𝐸𝐸
𝐴𝐴π‘₯π‘₯ 𝜌𝜌𝜌𝜌𝜌𝜌
𝐸𝐸
πœ”πœ”πœ”πœ”οΏ½ 𝑑𝑑𝑑𝑑𝑑𝑑 οΏ½πœ”πœ”πœ”πœ”οΏ½ οΏ½ =
𝜌𝜌
π‘Šπ‘Š
𝜌𝜌
SOLUTIONS
1.
Find, from first principles, the velocity of longitudinal waves along a thin steel rod.
See the Course Handout.
2.
The figure shows a schematic of a magnetostriction oscillator. The frequency of the oscillator is
determined by the length of the nickel alloy rod, which generates an alternating voltage in the
surrounding coils equal to the frequency of the fundamental longitudinal vibration of the rod.
Determine the full length, L, of the rod for a frequency of 20 kHz.
We consider half of the rod, with boundary conditions of fixed at one end (x = 0), and free at the other
end (x = L/2). Note that in this situation we include only half of the total length of the rod.
πœ”πœ”πœ”πœ”
πœ”πœ”πœ”πœ”
𝑒𝑒(π‘₯π‘₯, 𝑑𝑑) = �𝐴𝐴1 𝑠𝑠𝑠𝑠𝑠𝑠 οΏ½ οΏ½ + 𝐴𝐴2 𝑐𝑐𝑐𝑐𝑐𝑐 οΏ½ οΏ½οΏ½ 𝑠𝑠𝑠𝑠𝑠𝑠(πœ”πœ”πœ”πœ”)
𝑐𝑐
𝑐𝑐
πœ•πœ•πœ•πœ•
πœ”πœ”
πœ”πœ”πœ”πœ”
πœ”πœ”
πœ”πœ”πœ”πœ”
= �𝐴𝐴1 οΏ½ οΏ½ 𝑐𝑐𝑐𝑐𝑐𝑐 οΏ½ οΏ½ βˆ’ 𝐴𝐴2 οΏ½ οΏ½ 𝑠𝑠𝑠𝑠𝑠𝑠 οΏ½ οΏ½οΏ½ 𝑠𝑠𝑠𝑠𝑠𝑠(πœ”πœ”πœ”πœ”)
πœ•πœ•πœ•πœ•
𝑐𝑐
𝑐𝑐
𝑐𝑐
𝑐𝑐
Rods - 8
Boundary condition at x = 0, u = 0. The A2 = 0.
At x = L/2,
πœ•πœ•πœ•πœ•
πœ•πœ•πœ•πœ•
= 0 for all time. Therefore
πœ”πœ”πœ”πœ”
οΏ½
2𝑐𝑐
This is only true if 𝑐𝑐𝑐𝑐𝑐𝑐 οΏ½
From which
= 0, so
πœ•πœ•πœ•πœ•
πœ”πœ”
πœ”πœ”πœ”πœ”
= 0 = �𝐴𝐴1 οΏ½ οΏ½ 𝑐𝑐𝑐𝑐𝑐𝑐 οΏ½ οΏ½ βˆ’ 0οΏ½
πœ•πœ•πœ•πœ•
𝑐𝑐
2𝑐𝑐
πœ”πœ”π‘›π‘› 𝐿𝐿
1
= �𝑛𝑛 βˆ’ οΏ½ πœ‹πœ‹
2
2𝑐𝑐
Rearrange and solve for L with n =1.
𝐿𝐿 =
𝑓𝑓𝑛𝑛 =
with n = 1, 2, 3, …
πœ”πœ”π‘›π‘› �𝑛𝑛 βˆ’ 1οΏ½2οΏ½ 𝐸𝐸
=
οΏ½
2πœ‹πœ‹
𝐿𝐿
𝜌𝜌
οΏ½1 βˆ’ 1οΏ½2οΏ½ 𝐸𝐸 οΏ½1 βˆ’ 1οΏ½2οΏ½ 206 × 109
οΏ½
οΏ½ =
= 0.128 m = 128 mm
7843
𝜌𝜌
20,000
𝑓𝑓1
The total length of the rod must be 128.2 mm.
3.
(extra credit) A uniform rod of length L and cross-sectional area Ax is fixed at the upper end and
is loaded with a weight W=mg on the other end. Show that the natural frequencies are determined from
the equation:
𝐸𝐸
𝐴𝐴π‘₯π‘₯ 𝜌𝜌𝜌𝜌𝜌𝜌
𝐸𝐸
πœ”πœ”πœ”πœ”οΏ½ 𝑑𝑑𝑑𝑑𝑑𝑑 οΏ½πœ”πœ”πœ”πœ”οΏ½ οΏ½ =
𝜌𝜌
π‘Šπ‘Š
𝜌𝜌
The solutions to extra credit problems are available from your instructor.
Rods - 9