Dynamic Analysis
Dynamic Analysis
Dynamic Analysis
DYNAMIC ANALYSIS
Prof. A. S. Moharir
Piping systems in process plants are subjected to a variety of loads. It is
customary to categorize these loads as sustained loads, thermal loads and occasional
loads. Dynamic analysis pertains to the behaviour of piping systems under occasional
loads. While the origin of these loads could be in process related or layout related or
ambient related causes, these loads are typically oscillatory or impulsive in nature. These
can be magnified due to resonance with the natural frequency of the piping systems or its
specific segments. These can result into vibrations of large enough amplitude to cause
severe bending moments and stress cyclicity, which eventually could lead to fatigue
failure. Dynamic analysis is done to safeguard piping systems against such failures.
This articles attempts to prepare the reader with basic background to appreciate
the causes of vibrations, their quantification, basic mathematical description of the
dynamic phenomena, possible vibration arresters in real life and the code stipulation as
regards the stresses that are allowable under combination of static, thermal and
occasional loads.
equilibrium
position
m
x0
Consider a massless spring attached at one end to a rigid support and at another
end to a mass which can slide over a horizontal surface. Friction between the sliding mass
and the surface is assumed to be zero.
While in this condition, the spring will be relaxed. The mass is not moving and
hence the system has no kinetic energy as well. If some work is done on the system and it
is pulled a distance x 0 in the horizontal plane, the work would be stored as energy of the
system. The spring acquires potential energy ( = kx0 , k being spring constant) which is
equal to this work done. If the mass held stationary at x 0 is released now, it starts
oscillating. The spring-mass systems energy is conserved during this oscillations. Only
the potential energy of the spring and the kinetic energy of the mass interchange forms.
mv + kx = m(v + dv ) + k ( x + dx )
or
mvdv + kxdx = 0
As
v=
dx
dv
=0
+ kx
dt
dt
dv d 2 x
dx
and
=
, we get
dt dt 2
dt
d 2x
dx
m. 2 + kx
=0
dt
dt
dx
or velocity cannot be zero at every time (except
dt
when the mass is at its two extremes),
d 2x
m 2 + kx = 0
dt
or
d 2x k
+ x =0
dt 2 m
(i.e. at x = x 0 or x 0 )
x = c1 cos t + c 2 sin t
dx
= c1 sin t + c 2 cos t
dt
Using the initial conditions,
c1 = x 0 and c2 = 0
The solution describing the oscillation of the ideal mass-spring system is thus
x = x0 cos t
dx
and
= wx 0 sin t .
dt
Substituting the solution in the differential equation, one gets the circular frequency ( )
= k m radians/s
The mass exhibits a harmonic motion with the following frequency and period.
1
Frequency f =
cycles/s or Hz
k m
2
Period
= 2 m k
seconds.
The x vs t and dx dt = v vs t can be plotted. The ideal spring mass system thus
undergoes undamped oscillations over infinite period, once it is disturbed/perturbed.
What happens if the mass is not at rest (v = 0) when it is at x 0 at time t = 0 but
has a velocity v0 ?
What happens if the mass was hung to the spring vertically (see figure below) and
then perturbed by a vertical push (up or down)?
k
m
F0
mv + kx = m(v + dv ) + k ( x + dx ) F0 dx
or
mvdv + kxdx = F0 dx
As
v=
dv
dx
dx
+ kx
= F0
dt
dt
dt
dv d 2 x
dx
and cannot be zero all the time and with
=
, one gets
dx dt 2
dt
m
d 2x
+ kx = F0
dt 2
or
d 2x k
k F0
+ x=
2
m
m k
dt
or
d 2x
+ 2 x = 2 (F0 k )
2
dt
If the system was not oscillating, the force F 0 would have caused a deflection
F0 k . Let us call this static deflection ( x stat ).
Therefore, the equation governing the undamped oscillations of a spring-mass
system subjected to a steady force F 0 is
d 2x
+ 2 x = 2 x stat
2
dt
s.t.
x max
= 2.
x stat
Dynamic response in this case is twice the static response. This is an important
observation.
Plot the deflection and velocity vs t. What are the amplitude, frequency and period of
vibrations?
Spring-Mass System with Friction
k
m
Frictional force
f (moving left)
f (moving right)
Sliding of a mass in reality will encounter friction with the horizontal surface.
Similarly, movement (oscillations) of a pipe would encounter friction with ambient air
and end connections or supports. Energy will thus get dissipated with each cycle and one
would expect the oscillating system to come to a halt after a finite number of cycles after
initial perturbation.
This situation is similar to the previous case, but with a difference. The friction
force acts on the mass in the direction opposing its motion. It thus acts from left to right
as the mass swings from extreme right position to its extreme left position (first half
cycle). It acts from right to left as the mass swings from leftmost to rightmost position
(second half cycle)
Carry out the derivation and solution to arrive at the displacement vs time plot of
these damped oscillations. Show that amplitude reduces by (2 f k ) in magnitude every
v
v
f D = cv
Coulomb
Viscous
The governing equations and the solutions can be obtained for a spring-mass system with
viscous damping.
mv + kx = v m + kx = f D v
2
dt 2
dt
f D = cv
dx
d 2x
m 2 + c + kx = 0
dt
dt
k
m
c
c
c
=
=
2m 2 km cc
d = 1 2 .
x max = x0
2
v
0
1 +
+
d x0 d
= tan 1 [ d x0 (v 0 + x 0 )]
The value of is very important and decides the behaviour of the damped oscillator.
>1
<1
=1
k
m
F0
F
F
d 2x
dx
+ 2
+2x = 0 = 2 0
2
dt
m
k
dt
General solution
F
x(t ) = e t ( A cos d t + B sin d t ) + 0
k
Particular solution here is the static displacement and complimentary solution is the
vibration solution with external force.
For x = 0 at t = 0 and v = 0 at t = 0 ,
x(t ) =
F0
1 e t cos d t +
sin d t
k
1 2
k
m
F( t )
F( t )
t
Let the external force be
F (t ) = F0 sin t
The equations and solutions are:
d 2x
dx
+ 2
+ 2 x = (F0 m )sin t
2
dt
dt
10
x max =
(F0 k )
2 2
1 + (2 )2
= tan 1 2 1
1
2 2
1 + (2 )2
11
Oscillations
k (lateral stiffness)
The spring is ideal and mass-less. The mass is perturbed in a lateral direction and
oscillates. If a force F causes a displacement x at static equilibrium, F x is called the
lateral stiffness. Find the frequency and period of undamped oscillation.
Also analyse the following system with steady force F0 acting on the mass in
lateral direction.
Also analyse a spring-mass system with viscous dampner.
m
F0
k
c
12
k
c
Z (t ) = Z 0 sin t
Consider a system with two masses and 3 springs. For simplification let the
masses be same and springs identical.
x1 (t )
k
x 2 (t )
k
k
m
dt 2
2
2
2
2
d 2 x1
d 2 x2
v1 m 2 + kx1 k ( x 2 x1 ) + v 2 m 2 + kx2 + k ( x 2 x1 ) = 0
dt
dt
13
Each bracketed expression must vanish independently for this to be true. Thus,
d 2 x1
+ 2 x1 2 ( x 2 x1 ) = 0
2
dt
d 2 x2
+ 2 x 2 + 2 ( x 2 x1 ) = 0
2
dt
or
d 2 x1
+ 2 2 x1 2 x 2 = 0
dt 2
d 2 x2
2 x1 + 2 2 x 2 = 0
2
dt
Visual inspection of the governing equations indicate that two modes of oscillations of
the system can satisfy the system of equations.
Mode 1 : x1 (t ) = x2 (t )
Mode 2 : x1 (t ) = x2 (t )
Let us interpret these modes. In mode 1, we hold both the masses and move them a
distance x 0 to the right (or left) and leave them. Symmetry of the situation tells us that
the masses would oscillate in phase with same frequency and amplitude.
In mode 2, we hold both the masses, move them towards each other by a distance x 0
each and then leave them to oscillate. Symmetry again tells us that the two masses would
oscillate with same amplitude and frequency, but out of phase.
Let in mode 1, the two masses oscillate with a frequency 1 . The expected solution is
then
x1 1
x = 1 x0 cos 1t
2
(i.e. x1 = x 2 = x0 cos 1t )
Substituting this in any of the two equations of motion, we get
1 2 =
k
1 = k m
m
14
f1 =
1
1
=
k m
2 2
x1 1
x = 1 x0 cos 2 t
2
Substituting in the equation of motion, one gets
22 =
3k
3k
2 =
m
m
f2 =
2
1
=
2 2
3k
m
The system thus appears much stiffer while vibrating in mode 2 than in mode 1. Can one
offer an explanation for this?
Visualization of the modes for this spring-mass system was quite easy. For more complex
system, this may not be. We therefore look for some mathematics behind identification
and definition of the modes.
The two equations of motion can be written in the matrix form as
d 2 x1
2 k 2 1 x 0
1
dt +
=
2
d x 2 m 1 2 x 2 0
2
dt
x1 1
1
The solution vectors for mode 1 and were = x0 cos 1t and x0 cos 2 t
1
x 2 1
1
1
respectively. The vectors and are eigen vectors of the stiffness matrix
1
1
2 1
1 2 . Each eigen vector thus indeed defines a mode.
In a more complex system, the stiffness matrix will have larger dimensions and larger
number of eigen vectors and modes of vibration. Frequency of each mode can be found
out by using the solution vector defined by eigen vector for that mode. Natural
frequencies of vibration can then be defined, arranged in ascending order of stiffness.
15
We have now all the ingredients that go in the dynamic analysis of a complex
piping system.
The piping system is defined for the stress analysis s/w to carry out dead weight
analysis and static flexibility analysis. For dynamic analysis, this is then converted to an
equivalent spring-mass network. The mass and spring system with appropriate values
attached to mass and spring stiffness simulate the vibration characteristics of a piping
section between any two modes. Coulomb damping, viscous damping can be added
wherever necessary. The system under analysis is thus transformed to a 3-D spring-mass
network with spring straining possible in axial/lateral direction. Modes of vibration are
then identified through eigen vector approach. Behaviour of the system under steady
force, sinusoidal force or base excitation can then be studied. These external
forces/displacements could be for variety of reasons such as seismicity, wind load, water
hammer, relief, blowdown, machine vibrations, etc. Appropriate modes would resonate
and cause displacements. Maximum displacement of a section is used to calculate
bending moments and stresses at nodes. These must be within allowable limits stipulated
by codes.
If stress anywhere crosses allowable limits, frequency of vibration is altered to
remove or reduce resonance possibilities. These measures could be by additional
supports, snubbers, etc.
A system safe from dead weight, thermal flexibility and occasional load
considerations is then passed for construction.
16