System of Dynamics
System of Dynamics
Part 1:Design
1. We have the Free Body Diagram(FBD) as follows. By force balance, we get the equation of motion as
x
x, x,
kx
m
cx
Figure 1: FBD for Q1
m
x = kx cx m
x + cx + kx = 0
Now, assuming trial solution x = Aet and substituting above, we get on simplifying,
Aet (m 2 + c + k) = 0 Aet = x(t) = 0 or m 2 + c + k = 0
Case 1(x(t) = 0) is trivial and uninteresting/negligible solution. Case 2 gives solution to this quadratic equation as
c c2 4mk
c + c2 4mk
, 2 =
1 =
2m
2m
This gives on substitution to x(t), the two non-trivial basis solutions as follows.
x1 (t) = Ae1 t , x2 (t) = Ae2 t
And the general solution is a linear combination of these two basis solutions where the initial conditions decide the
combination coefficients below.
x(t) = c1 x1 (t) + c2 x2 (t) = C1 e1 t + C2 e2 t
Now, given the problem data for parameters and initial conditions as follows.
Mass m (kg) x(0) = x0 (m) x(0)(initial
4mk c2
c
1 = a + b, , 2 = a b
where b =
, a=
2m
2m
Now, using Euler theorem,
x(t) = C1 eat (cos bt + sin bt)) + C2 ea (cos bt sin bt))
= eat ((C1 + C2 ) cos bt + (C1 C2 ) sin bt)
Assuming C1 + C2 = cos and C1 C2 = sin and simplifying, we get solution of the form as follows.
x(t) = Keat cos(bt ), where K & are arbitrary constants
dx(t)
dv(t)
v(t) =
, a(t) =
dt
dt
2. We can see maximum displacement is defined by the amplitude term in solution namely Keat and maximum accelaration also occurs at the same point but scaled by b2 ie
Keat dmax & b2 Keat amax Keat Keat 0.0075 & b2 Keat 0.75
The coding for the same is done in MATLAB and script is attached. Please run the file Q1.m to get plot.