By: Team 16
File:R1.1.jpg
From Figure: Mass, Spring, Damper system with spring and damper in parallel
First, we take the + y direction to be downward and z to be out of the page.
Within this coordinate system, the displacement is described as a vector quantity y with the velocity and acceleration of the system denoted by
y
˙
{\displaystyle {\dot {y}}\!}
and
y
¨
{\displaystyle {\ddot {y}}\!}
respectively.
Multiplying by the mass denoted
m
{\displaystyle m\!}
we derive the kinematics of the system to be
m
y
¨
+
F
1
=
F
y
{\displaystyle m{\ddot {y}}+F1=F_{y}\!}
F
1
=
F
k
+
F
c
{\displaystyle F1=F_{k}+F_{c}\!}
To develop the kinetics of the system, we consider the two forces:
For the force limiting displacement in the y direction, we look at Hooke's law
F
k
=
−
k
y
{\displaystyle F_{k}=-ky\!}
For the force opposing the motion, the damping is related to the velocity of the system by
F
c
=
−
c
y
˙
{\displaystyle F_{c}=-c{\dot {y}}\!}
where
c
{\displaystyle c\!}
is the coefficient of viscous friction
Equating the kinematics to the kinetics we preliminarily have
m
y
¨
+
k
y
+
c
y
¨
=
F
y
¨
{\displaystyle m{\ddot {y}}+ky+c{\ddot {y}}={\ddot {F_{y}}}\!}
and simplifying with:
y
=
y
k
=
y
c
{\displaystyle y=y_{k}=y_{c}\!}
y
˙
=
y
c
˙
=
y
k
˙
{\displaystyle {\dot {y}}={\dot {y_{c}}}={\dot {y_{k}}}\!}
Final Equation:
F
(
t
)
=
m
y
¨
+
c
y
˙
+
k
y
{\displaystyle F(t)=m{\ddot {y}}+c{\dot {y}}+ky\!}
Author: Nyal Jennings
Last Edit: 01/28/2012 - 11:01 PM
Reviewed By: Martin Dolgiej
Example alt text Figure 1
From Figure 1 -
Spring:
−
k
y
=
F
k
{\displaystyle \displaystyle -ky=F_{k}}
(
E
q
.1
)
{\displaystyle \displaystyle (Eq.1)}
Mass:
m
y
″
=
F
m
{\displaystyle \displaystyle my''=F_{m}}
(
E
q
.2
)
{\displaystyle \displaystyle (Eq.2)}
Damp:
c
y
′
=
F
c
{\displaystyle \displaystyle cy'=F_{c}}
(
E
q
.3
)
{\displaystyle \displaystyle (Eq.3)}
Forcing Funtion:
r
(
t
)
=
F
r
{\displaystyle \displaystyle r(t)=F_{r}}
(
E
q
.4
)
{\displaystyle \displaystyle (Eq.4)}
Balancing the forces in the y direction:
∑
F
y
=
F
m
−
F
r
=
F
k
−
F
c
{\displaystyle \sum F_{y}=F_{m}-F_{r}=F_{k}-F_{c}\!}
Substituting equations (1),(2),(3), and (4):
m
y
″
−
r
(
t
)
=
−
k
y
−
c
y
′
{\displaystyle my''-r(t)=-ky-cy'\!}
m
y
″
+
c
y
′
+
k
y
=
r
(
t
)
{\displaystyle my''+cy'+ky=r(t)\!}
Author: Martin Dolgiej
Last Edit: 01/28/2012 - 11:32 PM
Reviewed By: Matt Hilsenrath
FBD of Spring
f
s
=
k
∗
y
s
{\displaystyle f_{s}=k*y_{s}\!}
FBD of Dashpot
C
y
D
′
=
f
c
{\displaystyle Cy_{D}'=f_{c}\!}
FBD of Mass
(Newton's 2nd Law)
f
=
m
y
″
{\displaystyle f=my''\!}
∑
F
(
t
)
=
m
y
″
+
f
s
−
f
c
{\displaystyle \sum F(t)=my''+f_{s}-f_{c}\!}
f
1
=
f
s
−
F
c
{\displaystyle f_{1}=f_{s}-F_{c}\!}
F
(
t
)
=
m
y
″
+
f
1
{\displaystyle F(t)=my''+f_{1}\!}
Author: Ryan McDaniel
Reviewed by: Erin Standish
Given :
V
=
L
C
d
2
v
C
d
t
2
+
R
C
d
v
C
d
t
+
v
C
{\displaystyle \displaystyle V=LC{\frac {d^{2}v_{C}}{dt^{2}}}+RC{\frac {dv_{C}}{dt}}+v_{C}}
(
E
q
.2
)
{\displaystyle \displaystyle (Eq.2)}
Solution :
Deriving equation (3):
Taking the derivative of (Eq. 2) gives:
V
′
=
L
C
d
3
v
C
d
t
3
+
R
C
d
2
v
c
d
t
2
+
d
v
C
d
t
{\displaystyle \displaystyle V'=LC{\frac {d^{3}v_{C}}{dt^{3}}}+RC{\frac {d^{2}v_{c}}{dt^{2}}}+{\frac {dv_{C}}{dt}}}
(
E
q
.
A
)
{\displaystyle \displaystyle (Eq.A)}
By definition:
i
=
C
d
v
C
d
t
{\displaystyle i=C{\frac {dv_{C}}{dt}}}
i
′
=
C
d
2
v
C
d
t
2
{\displaystyle i'=C{\frac {d^{2}v_{C}}{dt^{2}}}}
i
″
=
C
d
3
v
C
d
t
3
{\displaystyle i''=C{\frac {d^{3}v_{C}}{dt^{3}}}}
Plugging theses definitions into (Eq. A):
V
′
=
L
I
″
+
R
I
′
+
I
C
{\displaystyle \displaystyle V'=LI''+RI'+{\frac {I}{C}}}
(
E
q
.3
)
{\displaystyle \displaystyle (Eq.3)}
Deriving equation (4):
By definition:
Q
=
C
v
c
{\displaystyle Q=Cv_{c}\!}
Q
′
=
C
d
v
C
d
t
{\displaystyle Q'=C{\frac {dv_{C}}{dt}}}
Q
″
=
C
d
2
v
C
d
t
2
{\displaystyle Q''=C{\frac {d^{2}v_{C}}{dt^{2}}}}
Using Q and its time derivatives, (2) can be written as
V
=
L
Q
″
+
R
Q
′
+
Q
C
{\displaystyle \displaystyle V=LQ''+RQ'+{\frac {Q}{C}}}
(
E
q
.4
)
{\displaystyle \displaystyle (Eq.4)}
Authors: Erin Standish , James Welch
Reviewed by: Brandon Wright Ryan McDaniel
For question 1.5a we have to find the general solution to the equation:
y
″
+
4
y
′
+
(
π
2
+
4
)
y
=
0
{\displaystyle \displaystyle y''+4y'+(\pi ^{2}+4)y=0}
{\displaystyle \displaystyle }
(1)
and we have the relation
y
=
e
λ
x
{\displaystyle \displaystyle y=e^{\lambda x}}
{\displaystyle \displaystyle }
(2)
Its first two derivatives are
y
′
=
λ
e
λ
x
{\displaystyle \displaystyle y'=\lambda e^{\lambda x}}
{\displaystyle \displaystyle }
(3)
and
y
″
=
λ
2
e
λ
x
{\displaystyle \displaystyle y''=\lambda ^{2}e^{\lambda x}}
{\displaystyle \displaystyle }
(4)
Incorporating (2), (3), and (4) into equation (1),
it can be shown that
λ
2
e
λ
x
+
4
λ
e
λ
x
+
(
π
2
+
4
)
e
λ
x
=
0
{\displaystyle \displaystyle \lambda ^{2}e^{\lambda x}+4\lambda e^{\lambda x}+(\pi ^{2}+4)e^{\lambda x}=0}
{\displaystyle \displaystyle }
Factoring out the exponential term,
e
λ
x
(
λ
2
+
4
λ
+
π
2
+
4
)
=
0
{\displaystyle \displaystyle e^{\lambda x}(\lambda ^{2}+4\lambda +\pi ^{2}+4)=0}
{\displaystyle \displaystyle }
We recognize that an exponential term cannot equal zero, so it can be divided out.
Now we are left with the characteristic equation
λ
2
+
4
λ
+
π
2
+
4
=
0
{\displaystyle \displaystyle \lambda ^{2}+4\lambda +\pi ^{2}+4=0}
{\displaystyle \displaystyle }
Recognizing the relation, (7) can be factored to
(
λ
+
2
)
2
+
π
2
=
0
{\displaystyle \displaystyle (\lambda +2)^{2}+\pi ^{2}=0}
{\displaystyle \displaystyle }
to reveal the roots of the characteristic equation to be
λ
=
−
2
+
π
i
{\displaystyle \displaystyle \lambda =-2+\pi i}
{\displaystyle \displaystyle }
and
λ
=
−
2
−
π
i
{\displaystyle \displaystyle \lambda =-2-\pi i}
{\displaystyle \displaystyle }
so the standard form of the equation can be written
y
=
e
−
2
x
(
A
c
o
s
(
π
x
)
+
B
s
i
n
(
π
x
)
)
{\displaystyle \displaystyle y=e^{-2x}(Acos(\pi x)+Bsin(\pi x))}
{\displaystyle \displaystyle }
For question 1.5b we have to find the general solution to the equation:
y
″
+
2
π
y
′
+
π
2
y
=
0
{\displaystyle \displaystyle y''+2\pi y'+\pi ^{2}y=0}
{\displaystyle \displaystyle }
and we have the relation
y
=
e
λ
x
{\displaystyle \displaystyle y=e^{\lambda x}}
{\displaystyle \displaystyle }
its first two derivatives are
y
′
=
λ
e
λ
x
{\displaystyle \displaystyle y'=\lambda e^{\lambda x}}
{\displaystyle \displaystyle }
and
y
″
=
λ
2
e
λ
x
{\displaystyle \displaystyle y''=\lambda ^{2}e^{\lambda x}}
{\displaystyle \displaystyle }
Incorporating these equations it can be shown that
λ
2
e
λ
x
+
2
π
λ
e
λ
x
+
π
2
e
λ
x
=
0
{\displaystyle \displaystyle \lambda ^{2}e^{\lambda x}+2\pi \lambda e^{\lambda x}+\pi ^{2}e^{\lambda x}=0}
{\displaystyle \displaystyle }
Dividing out,
e
λ
x
{\displaystyle \displaystyle e^{\lambda x}}
{\displaystyle \displaystyle }
we have
e
λ
x
(
λ
2
+
2
π
λ
+
π
2
)
=
0
{\displaystyle \displaystyle e^{\lambda x}(\lambda ^{2}+2\pi \lambda +\pi ^{2})=0}
{\displaystyle \displaystyle }
and since
e
λ
x
{\displaystyle \displaystyle e^{\lambda x}}
{\displaystyle \displaystyle }
can't equal 0, we must evaluate the roots of the equation.
(
λ
2
+
2
π
λ
+
π
2
)
=
0
{\displaystyle \displaystyle (\lambda ^{2}+2\pi \lambda +\pi ^{2})=0}
{\displaystyle \displaystyle }
Factoring we have,
(
λ
+
π
)
2
=
0
{\displaystyle \displaystyle (\lambda +\pi )^{2}=0}
{\displaystyle \displaystyle }
Which places a double root at
x
=
−
π
{\displaystyle \displaystyle x=-\pi }
{\displaystyle \displaystyle }
So to describe the system in standard form,
y
=
(
c
1
+
c
2
x
)
e
−
π
x
{\displaystyle \displaystyle y=(c_{1}+c_{2}x)e^{-\pi x}}
{\displaystyle \displaystyle }
Author: Brandon Wright
Reviewed by: Matt Hilsenrath
Example alt text Figure 1
Order: Second Order ODE
Linearity: Linear
Superposition:
y
¯
(
x
)
=
y
p
(
x
)
+
y
h
(
x
)
{\displaystyle {\bar {y}}(x)=y_{p}(x)+y_{h}(x)\!}
From Figure 1 -
y
″
=
g
{\displaystyle y''=g\!}
Splitting differential equation into particular and homogeneous:
y
p
″
=
g
{\displaystyle y_{p}''=g\!}
y
h
″
=
0
{\displaystyle y_{h}''=0\!}
Summing above equations together:
y
p
″
+
y
h
″
=
g
{\displaystyle y_{p}''+y_{h}''=g\!}
(
y
p
+
y
h
)
″
=
g
=
y
¯
(
x
)
″
{\displaystyle (y_{p}+y_{h})''=g={\bar {y}}(x)''\!}
Therefore superposition IS valid.
Example alt text Figure 2
Order: First Order ODE
Linearity: Non-linear
Superposition:
v
¯
(
x
)
=
v
p
(
x
)
+
v
h
(
x
)
{\displaystyle {\bar {v}}(x)=v_{p}(x)+v_{h}(x)\!}
From Figure 2 -
m
v
′
=
m
g
−
b
v
2
{\displaystyle mv'=mg-bv^{2}\!}
Splitting differential equation into particular and homogeneous:
m
v
p
′
+
b
v
p
2
=
m
g
{\displaystyle mv_{p}'+bv_{p}^{2}=mg\!}
m
v
h
′
+
b
v
h
2
=
0
{\displaystyle mv_{h}'+bv_{h}^{2}=0\!}
Summing above equations together:
m
v
p
′
+
m
v
h
′
+
b
v
p
2
+
b
v
h
2
=
m
g
{\displaystyle mv_{p}'+mv_{h}'+bv_{p}^{2}+bv_{h}^{2}=mg\!}
m
(
v
p
+
v
h
)
′
+
b
(
v
p
2
+
b
v
h
2
)
=
m
g
≠
m
(
v
p
+
v
h
)
′
+
b
(
v
p
+
b
v
h
)
2
{\displaystyle m(v_{p}+v_{h})'+b(v_{p}^{2}+bv_{h}^{2})=mg\not =m(v_{p}+v_{h})'+b(v_{p}+bv_{h})^{2}\!}
Therefore, superposition is NOT valid.
Example alt text Figure 3
Order: First Order ODE
Linearity: Non-linear
Superposition:
h
¯
(
x
)
=
h
p
(
x
)
+
h
h
(
x
)
{\displaystyle {\bar {h}}(x)=h_{p}(x)+h_{h}(x)\!}
From Figure 3 -
h
′
=
−
k
h
1
/
2
{\displaystyle h'=-kh^{1/2}\!}
Splitting differential equation into particular and homogeneous:
h
h
′
=
−
k
h
h
1
/
2
{\displaystyle h_{h}'=-kh_{h}^{1/2}\!}
h
p
′
=
−
k
h
p
1
/
2
{\displaystyle h_{p}'=-kh_{p}^{1/2}\!}
Summing above equations together:
h
h
′
+
h
p
′
=
−
k
h
h
1
/
2
−
k
h
p
1
/
2
{\displaystyle h_{h}'+h_{p}'=-kh_{h}^{1/2}-kh_{p}^{1/2}\!}
(
h
h
+
h
p
)
′
=
−
k
h
h
1
/
2
−
k
h
p
1
/
2
≠
−
k
(
h
h
+
h
p
)
1
/
2
{\displaystyle (h_{h}+h_{p})'=-kh_{h}^{1/2}-kh_{p}^{1/2}\not =-k(h_{h}+h_{p})^{1/2}\!}
Therefore superposition is NOT valid.
Example alt text Figure 4
Order: Second Order ODE
Linearity: Linear
Superposition:
y
¯
(
x
)
=
y
p
(
x
)
+
y
h
(
x
)
{\displaystyle {\bar {y}}(x)=y_{p}(x)+y_{h}(x)\!}
From Figure 4 -
m
y
″
+
k
y
=
0
{\displaystyle my''+ky=0\!}
Splitting differential equation into particular and homogeneous:
m
y
h
″
+
k
y
h
=
0
{\displaystyle my_{h}''+ky_{h}=0\!}
m
y
p
″
+
k
y
p
=
0
{\displaystyle my_{p}''+ky_{p}=0\!}
Summing above equations together:
m
y
h
″
+
m
y
p
″
+
k
y
h
+
k
y
p
=
0
{\displaystyle my_{h}''+my_{p}''+ky_{h}+ky_{p}=0\!}
m
(
y
h
+
y
p
)
″
+
k
(
y
h
+
y
p
)
=
0
=
m
y
¯
(
x
)
″
+
k
y
¯
(
x
)
{\displaystyle m(y_{h}+y_{p})''+k(y_{h}+y_{p})=0=m{\bar {y}}(x)''+k{\bar {y}}(x)\!}
Therefore superposition IS valid.
Example alt text Figure 5
Order: Second Order ODE
Linearity: Linear
Superposition:
y
¯
(
x
)
=
y
p
(
x
)
+
y
h
(
x
)
{\displaystyle {\bar {y}}(x)=y_{p}(x)+y_{h}(x)\!}
From Figure 5 -
y
″
+
w
0
2
y
=
cos
w
t
{\displaystyle y''+w_{0}^{2}y=\cos {wt}\!}
Splitting differential equation into particular and homogeneous:
y
h
″
+
w
0
2
y
h
=
0
{\displaystyle y_{h}''+w_{0}^{2}y_{h}=0\!}
y
p
″
+
w
0
2
y
p
=
cos
w
t
{\displaystyle y_{p}''+w_{0}^{2}y_{p}=\cos {wt}\!}
Summing above equations together:
y
p
″
+
y
h
″
+
w
0
2
y
p
+
w
0
2
y
h
=
cos
w
t
{\displaystyle y_{p}''+y_{h}''+w_{0}^{2}y_{p}+w_{0}^{2}y_{h}=\cos {wt}\!}
(
y
p
+
y
h
)
″
+
w
0
2
(
y
p
+
y
h
)
=
cos
w
t
=
y
¯
(
x
)
″
+
w
0
2
y
¯
(
x
)
{\displaystyle (y_{p}+y_{h})''+w_{0}^{2}(y_{p}+y_{h})=\cos {wt}={\bar {y}}(x)''+w_{0}^{2}{\bar {y}}(x)\!}
Therefore superposition IS valid.
Example alt text Figure 6
Order: Second Order ODE
Linearity: Linear
Superposition:
I
¯
(
x
)
=
I
p
+
I
h
{\displaystyle {\bar {I}}(x)=I_{p}+I_{h}\!}
From Figure 6 -
L
I
″
+
R
I
′
+
1
C
I
=
E
′
{\displaystyle LI''+RI'+{\frac {1}{C}}I=E'\!}
Splitting differential equation into particular and homogeneous:
L
I
p
″
+
R
I
p
′
+
1
C
I
p
=
E
′
{\displaystyle LI_{p}''+RI_{p}'+{\frac {1}{C}}I_{p}=E'\!}
L
I
h
″
+
R
I
h
′
+
1
C
I
h
=
0
{\displaystyle LI_{h}''+RI_{h}'+{\frac {1}{C}}I_{h}=0\!}
Summing above equations together:
L
I
p
″
+
L
I
h
″
+
R
I
p
′
+
R
I
h
′
+
1
C
I
p
+
1
C
I
h
=
E
′
{\displaystyle LI_{p}''+LI_{h}''+RI_{p}'+RI_{h}'+{\frac {1}{C}}I_{p}+{\frac {1}{C}}I_{h}=E'\!}
L
(
I
p
+
I
h
)
″
+
R
(
I
p
+
I
h
)
′
+
1
C
(
I
p
+
I
h
)
=
E
′
=
L
I
¯
″
+
R
I
¯
′
+
1
C
I
¯
{\displaystyle L(I_{p}+I_{h})''+R(I_{p}+I_{h})'+{\frac {1}{C}}(I_{p}+I_{h})=E'=L{\bar {I}}''+R{\bar {I}}'+{\frac {1}{C}}{\bar {I}}\!}
Therefore superposition IS valid.
Example alt text Figure 7
Order: iv Order ODE
Linearity: Non-linear
Superposition:
y
¯
(
x
)
=
y
p
+
y
h
{\displaystyle {\bar {y}}(x)=y_{p}+y_{h}\!}
From Figure 7 -
E
I
y
i
v
=
f
(
x
)
{\displaystyle EIy^{iv}=f(x)\!}
Splitting differential equation into particular and homogeneous:
E
I
y
p
i
v
=
f
(
x
)
{\displaystyle EIy_{p}^{iv}=f(x)\!}
E
I
y
h
i
v
=
0
{\displaystyle EIy_{h}^{iv}=0\!}
Summing above equations together:
E
I
y
p
i
v
+
E
I
y
h
i
v
=
f
(
x
)
{\displaystyle EIy_{p}^{iv}+EIy_{h}^{iv}=f(x)\!}
E
I
(
y
p
i
v
+
y
h
i
v
)
=
f
(
x
)
≠
E
I
y
¯
i
v
=
E
I
(
y
p
+
y
h
)
i
v
{\displaystyle EI(y_{p}^{iv}+y_{h}^{iv})=f(x)\not =EI{\bar {y}}^{iv}=EI(y_{p}+y_{h})^{iv}\!}
Therefore superposition is NOT valid.
Example alt text Figure 8
Order: Second Order ODE
Linearity: Non-linear
Superposition:
θ
¯
=
θ
p
+
θ
h
{\displaystyle {\bar {\theta }}=\theta _{p}+\theta _{h}\!}
From Figure 8 -
L
θ
″
+
g
sin
θ
=
0
{\displaystyle L\theta ''+g\sin {\theta }=0\!}
Splitting differential equation into particular and homogeneous:
L
θ
h
″
+
g
sin
θ
h
=
0
{\displaystyle L\theta _{h}''+g\sin {\theta _{h}}=0\!}
L
θ
p
″
+
g
sin
θ
p
=
0
{\displaystyle L\theta _{p}''+g\sin {\theta _{p}}=0\!}
Summing above equations together:
L
θ
h
″
+
L
θ
p
″
+
g
sin
θ
h
+
g
sin
θ
p
=
0
{\displaystyle L\theta _{h}''+L\theta _{p}''+g\sin {\theta _{h}}+g\sin {\theta _{p}}=0\!}
L
(
θ
h
+
θ
p
)
″
+
g
(
sin
θ
h
+
sin
θ
p
)
=
0
≠
L
θ
¯
″
+
g
sin
θ
¯
{\displaystyle L(\theta _{h}+\theta _{p})''+g(\sin {\theta _{h}}+\sin {\theta _{p}})=0\not =L{\bar {\theta }}''+g\sin {\bar {\theta }}\!}
Therefore superposition is NOT valid.
Author: Martin Dolgiej
Last Edited: 01/28/2011 - 4:10 PM
Reviewed by:James Welch , Matt Hilsenrath