Chapter One 1
Chapter One 1
Chapter One 1
𝐹 E, A, L 𝐹
𝐹 𝑥 𝐴𝐸
𝜎 = 𝐸𝜖 → = 𝐸 → 𝐹 = (𝐴𝐸/𝐿)𝑥 → 𝑘 =
𝐴 𝐿 𝐿
J, G, L, c T
𝑇𝑐 𝜙𝑐 𝐽𝐺 𝐽𝐺
𝜏 = 𝐺𝛾 → =𝐺 →𝑇= 𝜙→𝑘=
𝐽 𝐿 𝐿 𝐿
𝐶𝑙𝑎𝑚𝑝𝑒𝑑 − 𝑐𝑙𝑎𝑚𝑝𝑒𝑑:
𝐸𝐼 𝐸𝐼
𝑎 𝐹 = 192 3 𝑥 → 𝑘 = 192 3
𝑙 𝑙 𝑎
𝐸, 𝐹
𝐼,
𝑃𝑖𝑛𝑛𝑒𝑑 − 𝑝𝑖𝑛𝑛𝑒𝑑 𝐸, 𝐹
3𝐸𝐼𝑙 𝐼,
𝑙 𝐹= 2 𝑥→ 𝑙
𝑙−𝑎 𝑎 (𝑙 − 𝑎)2
3𝐸𝐼𝑙 𝑙−𝑎
𝑘= 2
𝑎 (𝑙 − 𝑎)2
Springs: In vibration modeling springs are mass less elastic
elements, the ratio of the applied force to the corresponding
deflection or the slope of the force deflection plot is called
spring stiffness.
ΔF
Oscillation around 𝑥 = 𝑥𝑜 = 0 → 𝑘 = |𝑥𝑜 =0 F, [N]
Δx
1 2
𝐹 = 𝑘 𝑥, 𝑉𝑒 = 𝑘𝑥
2 x, [m]
k
k, [N/m]=slope
Linearization:.
ΔF F, [N]
Oscillation around 𝑥 = 𝑥𝑜 ≠ 0 →𝑘= |𝑥𝑜
Δx
𝑥 = 𝑥𝑜 x, [m]
Spring or elastic force 𝐹 is a function of the deflection 𝑥
measured from the equilibrium point 𝑥 = 𝑥𝑜 .
Springs in series and parallel connections:
𝑘1
𝑘2 𝐹1 𝐹2 𝐹𝑛
F 𝑘3 F 𝐹 𝐹
𝑘𝑛 𝑘1
𝑘2, 𝑘𝑛
Springs in parallel connections. Springs in series connections.
𝐹 = 𝑘1 𝑥 + 𝑘2 𝑥+. . .𝑘𝑛 𝑥= 𝐹 = 𝑘1 𝑥1 = 𝑘2 𝑥2 =, … =𝑘𝑛 𝑥𝑛
𝑘1 + 𝑘2 +. . .𝑘𝑛 𝑥 = 𝑘eq 𝑥 𝐹
𝑥 = 𝑥1 + 𝑥2 +. . . +𝑥𝑛 =
𝑘eq = 𝑘1 + 𝑘2 +. . .𝑘𝑛 𝑘eq
𝐹 𝐹 𝐹 𝐹
𝑥= + +. . . + =
𝑘1 𝑘2 𝑘𝑁 𝑘eq
𝑘1 𝑘2 𝑘3 . . .𝑘𝑛
𝑘eq =
𝑘1 + 𝑘2 + 𝑘3 +. . . +𝑘𝑛
1
𝐹 = 𝑘eq 𝑥; 𝑉𝑒 = 𝑘eq 𝑥 2
2
Mass: The mass used here represent equivalent lumped mass of a distributed one
in the mechanical system, it is an effective mass, or effective polar mass moment
of inertia. The mass is the only mean to store kinetic energy.
𝑚 𝑚
𝑉
. .
𝑂 .
𝑚 = 𝑚𝑒𝑞 𝜃 𝜃
1 𝑚𝑒𝑞 =𝐽𝑜
𝑇 = 𝑚𝑒𝑞 𝑉 2
2 1 2
1 .2
𝑇 = 𝑚𝑒𝑞 𝑉 = 𝐽𝑜 𝜃
2 2
𝑂
.
𝑚 𝐽𝑜 𝜃, 𝜃
1 2
1 .2
𝑇 = 𝑚𝑒𝑞 𝑉 = 𝐽𝑜 𝜃
2 2
Viscous Damper: Is an Energy dissipation element. The force act in
a viscous damper is proportional to the velocity. The
constant of proportionality is the damping coefficient c,
[N.s/m].
Equivalent damper for a group of dampers in series or parallel follow the same
rule as springs.
Examples: Model the following cantilever beam by an equivalent lumped mass at its
,
tip
,
and an equivalent spring (single degree of freedom system).
𝑦 𝑥 , and k have expressions you found in the material mechanics textbook
Equivalent mass of actual spring (Coil).
Find the equivalent mass of the spring at point A,
where y l , if the mass of the spring per unit length
is m . First find the Kinetic. Energy of the spring
1 1
dT dm[ x( y )]2 (mdy )[ x( y )]2 , (1-9)
2 2
Undamped Free Vibration Model:
No damping exists, and free since there is no external excitation
𝐹, 𝑁
Equilibrium position:
𝑚𝑔 + 𝑘𝛿 = 0
Undamped Free Vibration Model (cont’ed): The system vibrates because of initial
disturbance called initial conditions. These are initial displacement , velocity or
both displacement and velocity . If we select the position reference at the
equilibrium point, the static deflection force 𝑘 Δ𝑥 will cancel the weight 𝑊, (𝑊
= 𝑚𝑔), So whenever the reference datum we use to measure the motion is located
at the equilibrium position do not include the weight force .
Steps to write the equation of motion
1) Draw free body diagram,
Isolate the mass, chose
coordinate system and indicate
the positive direction,
determine whether gravity is
included or not.
Indicate (Draw) all forces act on the
body. Draw the inertia force in the
opposite positive direction.
2. Sum all forces: The summation of all
Equilibrium position 𝑘(𝑥 + 𝛿)
forces including the inertia is zero.
𝐹 = 0: 𝑚𝑔 − 𝑘 𝑥 + 𝛿 − 𝑚𝑥ሷ = 0 𝑚
𝑚𝑔 + 𝑘𝛿 = 0: eq. position
x, 𝑥,ሶ 𝑥ሷ
(𝑚𝑔 + 𝑘𝛿) − 𝑘 𝑥 + 𝛿 − 𝑚𝑥ሷ = 0 𝑚𝑥ሷ + 𝑘𝑥 = 0 𝑚𝑔 𝑚𝑥ሷ
Draw the inertia force in the opposite positive direction.
2) Sum all forces: The summation of all forces including the inertia is zero.
x 0 A cos() v0 1 ( v 0 ) tan 1 ( v 0 )
n tan( ), tan
v 0 n A sin() x0 n x 0 n x 0
(1-2d)
v0 2 2 2x 2
A sin() v v n 0
n 2 0 2
A 2 x0 A
0
2
x 0 A cos() n n
v0 2 2 2x 2
A sin() v v n 0
n 2 0 2
A 2 x0 A
0
2
x 0 A cos() n n
The solution is, v02 n2 x02
cos(n t tan 1(
v0
x(t)= ))
n2 n x0
v max
v 0 2m /s
x (t ), m A 1.2 m
x 0 1, m
t ,s
1 2
T ,s
f 3
7) Draw the velocity v(t), and displacement, x(t) and on the same p
v (t)
x(t)
t, s
9) Notice the amplitude of v(t), x(t), and a(t)and notice the phase shift between them. what is
the value of A for which the x (t ) v (t ) a(t )
Free Damped Vibrations:
Besides the spring- mass element model used
to model the un-damped system, for damped
systems we need to add some kind of energy
dissipation. The easiest way is to consider
the viscous damper, the force in the viscous
damper at any moment is proportional to the
velocity of the moving mass.
Find the solution for the displacement of
the mass m, x(t).
Where
A2
1
tan ( ) and A Ai2 A 22
A1
1 A1
or tan ( ) , D D12 D 22 , (1-5b)
A2
Apply Initial Conditions: x (0) x 0 , x (0) v 0
Substitute these initial conditions into the last
solution to get, two algebraic equations in terms of A1 and A2.
n t
x (t ) n e (A1 cos(d t ) A 2 sin(d t ))
e n t (d A1 sin(d t ) d A 2 cos(d t ))
x (0) v 0 n A1 d A2
v 0 n x 0
A2 , (1-5c)
d
So the solution is,
v 0 n x 0
x (t ) e n t (x 0 cos(d t ) sin(d t ))
d
v n x 0 2
x (t ) e n t x 02 ( 0 ) sin(d t )
d
n t 1 x0d
1
x(t ) e x02d2 (v0 n x0 ) sin{d t tan (
2
)} , (1-5d)
d v0 n x0
Draw the solution for
0.1,0.3,0.6, n 4 rad / s, ,
x0 1 m, v0 1m / s
You need to find A and
x(t)
3
2.0 11
0.5
2 2
n t ( x0 wd ) ( v0 z wn x0 )
2 e
wd
x(t)
1
n( wd t ph)
2.0 11
Set up your conclusions
x(t), [m] about the effect of on the response x(t ) .
Initial conditions on the starting point location and
slop.
t, [s]
0.6
0.3
0.1
Case 2: Over damped motion, 1,
1 n ( 2 1), 2 n ( 2 1) ,
c2 (3r )
c1(r )
Eq.
Oy Position
J
Ox
k (r ) x
Mo J , , (1-8a)
m 9m
J Jdisk Jrod (1/ 2)(m1)r2 (1/ 3)m2 (3r )2 r2 ( 1 2 ) , (1-8b)
2 3
m 9m
c1(r )r k (r )r c2 (3r )3r r2 ( 1 2 ) , , , (1-8c)
2 3
m 9m
r 2 ( 1 2 ) r 2 (c1 c2 9) r 2 k 0 , (1-8d)
2 3
1 1 , (1-8e)
m 9m 1 x (c
1 2c 9) x ( k ) x 0 ,
( 1 2) 3r 3r
2 3 3r
ceq keq
meq
(c1 c2 9) k
x x ( ) x 0 , (1-8f)
m1 9m2 m1 9m2
( ) ( )
2 3 2 3
2eq n n2
This is exactly the same equation of motion of the equivalent
linear system with mass m meq , k keq and c ceq .
Notice we found the equivalent spring stiffness keq ,
the equivalent mass meq and the equivalent damping ceq
as well as eq are based on the linear displacement coordinate x.
Energy Methods:
𝑑𝐸 𝑑𝑇 𝑑𝑉
𝐸 = 𝑇 + 𝑉, =0= +
𝑑𝑡 𝑑𝑡 𝑎𝑡
Datum
𝑳=𝑻−𝑼
𝒅 𝝏𝑳 𝝏𝑳
− =𝟎
𝒅𝒕 𝝏𝒒ሶ 𝝏𝒒
Datum
𝒒 = 𝜽ሶ
𝟏 𝟏 𝟏 𝟏 𝑳 𝟐 𝟏
𝑻= 𝑱𝑮 𝜽ሶ 𝟐 + 𝒎 𝒗𝟐 = 𝑱𝑮 𝜽ሶ 𝟐 + 𝒎 ( 𝜽)
ሶ = 𝑱𝒐 𝜽ሶ 𝟐
𝟐 𝟐 𝟐 𝟐 𝟐 𝟐
𝟏 𝟏 𝟏 𝑳
𝑻 = 𝑱𝒐 𝜽ሶ 𝒐𝒓 𝑻 = 𝑱𝑮 𝜽ሶ + 𝒎 𝒗
𝟐 𝟐 𝟐 𝑼 = 𝒎𝒈( (𝟏 − 𝒄𝒐𝒔𝜽)
𝟐
𝟐 𝟐 𝟐
𝟏 𝑳 𝝏𝑳 𝒅 𝝏𝑳
𝑳=𝑻−𝑼= 𝑱𝒐 𝜽ሶ 𝟐 − 𝒎𝒈( 𝟏 − 𝒄𝒐𝒔𝜽 , = 𝑱𝒐 𝜽,ሶ = 𝑱𝒐 𝜽ሷ
𝟐 𝟐 𝝏𝒒ሶ 𝒅𝒕 𝝏𝒒ሶ
𝝏𝑳 𝑳 𝒅 𝝏𝑳 𝝏𝑳 𝑳 𝑳
= 𝒎 𝒈 𝜽, − = 𝟎, 𝑱𝒐 𝜽ሷ − −𝒎𝒈 𝜽 = 𝟎 ሶ
𝑱𝒐 𝜽 + 𝒎 𝒈 𝜽 = 𝟎
𝝏𝒒 𝟐 𝒅𝒕 𝝏𝒒ሶ 𝝏𝒒 𝟐 𝟐
Example: Assume the wheel rolls without slipping, write the equation
of motion using the energy method.
A
Nonlinear system;Coulomb friction: Bilinear system
1 2
0
x(t)
-1.5
-0.5
0.5
1.5
-2
-1
2
Free body Diagram 0
i=1
2
1 n
2 2
0
n
3 2
i=2
1 n
4 2
t, [s]
n 0
5
i=3
1 n
6
n
2 mgn
Slope
k
stable system
If for other system the value of k is negative and m is positive, the solution
would be
Which no longer bounded for any time, instead increases without bound as t increases. |x(t)|
is no longer remain finite, and such solution is called unstable or divergent
Divergent
instability
Flutter
instability
Example: Inverted pendulum
Linearize,
Rearrange,
Spring Stiffness: can be found
knowing the geometry and 3
0 . 01
x(t1)
0 . 00 5 x(t2 )
x(t3 )
z t wn 0 . 00 2 5 x(t4 )
e xx sin( wd t p h)
z t wn x ,[m] t2
e xx
0 t1 0 . 05 0 . 1t3 0 . 15t4 0.2 0 . 25 0.3 0 . 35 0.4
z t wn
e xx
0 . 00 2 5 t ,[s]
0 . 00 5
0 . 00 7 5
Damping Coefficient:
7 . 98 41 0
3
0 . 01
0 t .4
0 . 01
3
7 . 98 41 0
0 . 00 7 5 x(t1)
0 . 00 5 x(t2 )
x(t3 )
z t wn 0 . 00 2 5 x(t4 )
e xx sin( wd t p h)
z t wn x ,[m] t2
e xx
0 t1 0 . 05 0 . 1t3 0 . 15t4 0.2 0 . 25 0.3 0 . 35 0.4
z t wn
e xx
0 . 00 2 5 t ,[s]
0 . 00 5
0 . 00 7 5
3
7 . 98 41 0
0 . 01
0 . 00 5 x(t2 )
x(t3 )
z t wn 0 . 00 2 5 x(t4 )
e xx sin( wd t p h)
z t wn x ,[m] t2
e xx
0 t1 0 . 05 0 . 1t3 0 . 15t4 0.2 0 . 25 0.3 0 . 35 0.4
z t wn
e xx
0 . 00 2 5 t ,[s]
0 . 00 5
0 . 00 7 5
3
7 . 98 41 0
0 . 01
2 2
n ,
0 t .4
n 1 2 1 2 4 2 2
0 . 00 5 x(t2 )
x(t3 )
z t wn 0 . 00 2 5 x(t4 )
e xx sin( wd t p h)
z t wn x ,[m] t2
e xx
0 t1 0 . 05 0 . 1t3 0 . 15t4 0.2 0 . 25 0.3 0 . 35 0.4
z t wn
e xx
0 . 00 2 5 t ,[s]
0 . 00 5
0 . 00 7 5
3
7 . 98 41 0
0 . 01
0 t .4
Consider the system of Figure 2.40(a). Let x
describe the downward displacement of m1
from the system’s equilibrium position.
(a) Derive the differential equation governing
x(t).
(b) Determine the reaction at the center of the
disk at the pin support in terms x(t), dx/dt and
dx^2/dt^2
At equilibrium: Summing moments about the
pin support to zero with positive moments Static Equilibrium
counterclockwise leads to
Use of Equation (a) in Equation (c) gives