R.1.1: Equation of Motion for a Spring-dashpot Parallel System in Series with Mass and Applied Force[ 1]
edit
Given a spring-dashpot system in parallel with an applied force, find the equation of motion.
For this problem, Newton's second law is used,
∑
F
=
m
a
{\displaystyle \sum F=ma}
(1-1 )
and applied to the mass at the end of a spring and damper in parallel.
Assuming no rotation of the mass,
y
=
y
m
=
y
k
=
y
c
{\displaystyle y=y_{m}=y_{k}=y_{c}}
(1-2 )
Therefore, the spring and damper forces can be written as, respectively,
f
k
=
k
y
{\displaystyle f_{k}=ky}
(1-3 )
and
f
c
=
c
y
′
{\displaystyle f_{c}=cy'}
(1-4 )
And the resultant force on the mass can be written as
F
=
m
a
=
m
y
″
{\displaystyle F=ma=my''}
Now, from equations (1-3 ) and (1-4 ) all of the forces can be substituted into Newton's second law, and since each distance is equal,
∑
F
=
m
y
″
=
−
f
k
−
f
c
+
f
(
t
)
{\displaystyle \sum F=my''=-f_{k}-f_{c}+f(t)}
(1-5 )
Substituting these exact force equations,
m
y
″
=
−
k
y
−
c
y
′
+
f
(
t
)
{\displaystyle my''=-ky-cy'+f(t)}
A little algebraic manipulation yields
m
y
″
+
c
y
′
+
k
y
=
f
(
t
)
{\displaystyle my''+cy'+ky=f(t)}
Finally, dividing by the mass to put the equation in standard form gives the final equation of motion for the mass:
y
″
+
c
m
y
′
+
k
m
y
=
1
m
f
(
t
)
)
{\displaystyle y''+{\frac {c}{m}}y'+{\frac {k}{m}}y={\frac {1}{m}}f(t))}
(1-6 )
↑ EGM4313: Section 1 Notes (.djvu)
R.1.2: Equation of Motion for a Spring-mass-dashpot System with Applied Force[ 1]
edit
For this problem, the task was to derive the equation of motion of the spring-mass-dashpot in Fig. 53, in K 2011 p.85 with an applied force r(t).
To solve this problem, note that the ODE of a damped mass-spring system is
m
y
″
+
c
y
′
+
k
y
=
0
{\displaystyle my''+cy'+ky=0}
(2-1 )
When there is a external formal added to the model r(t) on the right. This then gives up the equation
m
y
″
+
c
y
′
+
k
y
=
r
(
t
)
{\displaystyle my''+cy'+ky=r(t)}
(2-2 )
In this problem referencing back to equation (1-2 )
y
=
y
k
=
y
c
{\displaystyle y=y_{k}=y_{c}}
The resultant force for the system can be described as stated in equation (1-5 )
∑
F
=
m
y
″
=
−
f
k
−
f
c
+
f
(
t
)
{\displaystyle \sum F=my''=-f_{k}-f_{c}+f(t)}
Since r(t) is the external force,
f
(
t
)
=
r
(
t
)
{\displaystyle f(t)=r(t)}
The model of a mass-spring system ODE with external force on the right is modeled as
m
y
″
+
f
i
n
t
e
r
n
a
l
=
r
(
t
)
{\displaystyle my''+f_{internal}=r(t)}
The internal force is equal to the force of the spring making this equation
m
y
″
+
k
y
=
r
(
t
)
{\displaystyle my''+ky=r(t)}
When a dashpot is added, the force of the dashpot cy’ is added to the equation making it
m
y
″
+
c
y
′
+
k
y
=
r
(
t
)
{\displaystyle my''+cy'+ky=r(t)}
(2-3 )
For the characteristic equation m is divided throughout the whole equation making it
↑ EGM4313: Section 1 Notes (.djvu)
R.1.3: Free-body Diagram and Equation of Motion for Spring-dashpot-mass System[ 1]
edit
For this problem we were to draw the FBDs and derive the equation of motion for the following system.
To start solving this problem, we have to first look at how the spring, dashpot, and mass interact with each other by analyzing their respective free-body diagrams.
Next we take a look at what we already know.
We know from kinematics that:
y
=
y
k
+
y
c
{\displaystyle y=y{_{k}}+y{_{c}}}
(3-1 )
We also know from kinetics that:
m
y
″
+
f
i
n
t
e
r
n
a
l
=
f
(
t
)
{\displaystyle my''+f{_{internal}}=f(t)}
(3-2 )
f
i
=
f
k
=
f
c
{\displaystyle f{_{i}}=f{_{k}}=f{_{c}}}
(3-3 )
From the constitutive relations we know:
f
k
=
k
y
k
{\displaystyle f_{k}=ky_{k}}
(3-4 )
f
c
=
c
y
c
′
{\displaystyle f{_{c}}=cy{_{c}}'}
(3-5 )
We also know:
y
″
=
y
k
″
+
y
c
″
{\displaystyle y''=y{_{k}}''+y{_{c}}''}
(3-6 )
f
(
k
)
=
f
(
c
)
⇒
k
y
k
=
c
y
c
′
{\displaystyle f(k)=f(c)\Rightarrow ky{_{k}}=cy{_{c}}'}
(3-7 )
by solving for
y
c
′
{\displaystyle y_{c}'}
we get
y
c
′
=
k
c
y
k
{\displaystyle y_{c}'={\frac {k}{c}}y_{k}}
(3-8 )
To derive the equation of motion we have to manipulate and combine a few formulas that we know.
We first take what we found in (3-2 )
and substitute in
y
k
″
+
y
c
″
{\displaystyle y{_{k}}''+y{_{c}}''}
(3-9 )
for
y
″
{\displaystyle y''}
. From equations (3-3 ), (3-4 ), and (3-5 ) we can substitute
f
i
{\displaystyle f{_{i}}}
for
k
y
k
+
c
y
′
{\displaystyle ky{_{k}}+c{_{y}}'}
After the substitutions into equation (3-2 ) we have
m
(
y
k
″
+
y
c
″
)
+
k
y
k
+
c
y
′
=
f
(
t
)
{\displaystyle m(y{_{k}}''+y{_{c}}'')+ky{_{k}}+c{_{y}}'=f(t)}
(3-10 )
From the equation (3-8 ) we know
y
c
′
=
k
c
y
k
{\displaystyle y_{c}'={\frac {k}{c}}y_{k}}
We can then substitute this into equation (3-10 ) getting
m
(
y
k
″
+
k
c
y
k
′
)
+
k
y
k
+
k
c
y
k
=
f
(
t
)
{\displaystyle m(y{_{k}}''+{\frac {k}{c}}y{_{k}}')+ky{_{k}}+{\frac {k}{c}}y{_{k}}=f(t)}
To get the answer in its final form we divide both sides of the equation by m, and our final answer is
y
k
″
+
k
m
c
y
k
′
+
k
m
(
1
+
1
c
)
y
k
=
1
m
f
(
t
)
{\displaystyle y{_{k}}''+{\frac {k}{mc}}y{_{k}}'+{\frac {k}{m}}(1+{\frac {1}{c}})y{_{k}}={\frac {1}{m}}f(t)}
(3-11 )
↑ EGM4313: Section 1 Notes (.djvu)
R.1.4: Derivation of Voltage-Charge and Voltage-Current Relationships[ 1]
edit
For this problem, the goal is to derive two different equations from the circuit equation
V
=
L
C
d
2
v
c
d
t
2
+
R
C
d
v
c
d
t
+
v
c
{\displaystyle V=LC{\frac {d^{2}v_{c}}{dt^{2}}}+RC{\frac {dv_{c}}{dt}}+v_{c}}
(4-1 )
Each separate derivation is presented in Voltage-Charge Derivation and Voltage-Current Derivation
To solve these derivation, it is wise to note that capacitance is
Q
=
C
v
c
{\displaystyle Q=Cv_{c}}
(4-2 )
It can also be noted that
Q
=
∫
i
d
t
{\displaystyle Q=\int idt}
This means that
∫
i
d
t
=
C
v
C
{\displaystyle \int idt=Cv_{C}}
Completing the integration results in
I
=
C
d
v
C
d
t
{\displaystyle I=C{\frac {dv_{C}}{dt}}}
(4-3 )
It should also be noted that (4-2 ) can be written in the form
v
c
=
1
C
Q
{\displaystyle v_{c}={\frac {1}{C}}Q}
(4-4 )
It should also be noted that (4-3 ) can be written in the form
d
v
C
d
t
=
1
C
I
{\displaystyle {\frac {dv_{C}}{dt}}={\frac {1}{C}}I}
(4-5 )
Voltage-Charge Derivation
edit
For this problem, the Voltage-Charge equation
V
=
L
Q
″
+
R
Q
′
+
1
C
Q
{\displaystyle V=LQ''+RQ'+{\frac {1}{C}}Q}
(4-6 )
are derived from (4-1 )
Taking the derivative of (4-2 ) are taken with respect to time
Q
′
=
C
d
v
c
d
t
{\displaystyle Q'=C{\frac {dv_{c}}{dt}}}
(4-7 )
Q
″
=
C
d
2
v
c
d
t
2
{\displaystyle Q''=C{\frac {d^{2}v_{c}}{dt^{2}}}}
(4-8 )
Substituting (4-4 ), (4-7 ), and (4-8 ) into (4-1 ) results in (4-6 ) such that
L
C
d
2
v
c
d
t
2
+
R
C
d
v
c
d
t
+
v
c
=
L
Q
″
+
R
Q
′
+
1
C
Q
{\displaystyle LC{\frac {d^{2}v_{c}}{dt^{2}}}+RC{\frac {dv_{c}}{dt}}+v_{c}=LQ''+RQ'+{\frac {1}{C}}Q}
(4-9 )
Voltage-Current Derivation
edit
For this problem, the following equation
V
′
=
L
I
″
+
R
I
′
+
1
C
I
{\displaystyle V'=LI''+RI'+{\frac {1}{C}}I}
(4-10 )
are derived from (4-1 )
The first step is to take the derivative of (4-1 ) resulting in
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 V'=LC{\frac {d^{3}v_{c}}{dt^{3}}}+RC{\frac {d^{2}v_{c}}{dt^{2}}}+{\frac {dv_{c}}{dt}}}
(4-11 )
Taking the derivative of (4-3 ) are taken with respect to time
I
′
=
C
d
2
v
c
d
t
2
{\displaystyle I'=C{\frac {d^{2}v_{c}}{dt^{2}}}}
(4-12 )
I
″
=
C
d
3
v
c
d
t
3
{\displaystyle I''=C{\frac {d^{3}v_{c}}{dt^{3}}}}
(4-13 )
Substituting (4-3 ), (4-12 ), and (4-13 ) into (4-11 ) results in (4-10 ) such that
L
C
d
3
v
c
d
t
3
+
R
C
d
2
v
c
d
t
2
+
d
v
c
d
t
=
L
I
″
+
R
I
′
+
1
C
I
{\displaystyle LC{\frac {d^{3}v_{c}}{dt^{3}}}+RC{\frac {d^{2}v_{c}}{dt^{2}}}+{\frac {dv_{c}}{dt}}=LI''+RI'+{\frac {1}{C}}I}
(4-16 )
↑ EGM4313: Section 2 Notes (.djvu)
R.1.5: Solutions of General 2nd Order ODEs[ 1]
edit
Consider a second-order homogeneous linear ODE with constant coefficients a and b.
y
″
+
a
y
′
+
b
y
=
0
{\displaystyle y''+ay'+by=0}
(5-1 )
To solve this problem, note that the solution to a first-order linear ODE of the form:
y
′
+
k
y
=
0
{\displaystyle y'+ky=0}
is an exponential function, yielding a solution of the form:
y
=
e
λ
x
{\displaystyle y=e^{\lambda x}}
(5-2 )
Based upon this, we will expand the solution to apply to second-linear homogeneous ODEs. To do so, we must first find equations describing
y
′
{\displaystyle y'}
and
y
″
{\displaystyle y''}
, and so we will take the derivatives (with respect to x) of (5-2 ).
y
′
=
λ
e
λ
x
{\displaystyle y'=\lambda e^{\lambda x}}
(5-3 )
y
″
=
λ
2
e
λ
x
{\displaystyle y''=\lambda ^{2}e^{\lambda x}}
(5-4 )
We will now substitute (5-2 ), (5-3 ), and (5-4 ) into (5-1 ) to obtain the relationship:
λ
2
e
λ
x
+
a
λ
e
λ
x
+
b
e
λ
x
=
0
{\displaystyle \lambda ^{2}e^{\lambda x}+a\lambda e^{\lambda x}+be^{\lambda x}=0}
Simplifying, we have:
(
λ
2
+
a
λ
+
b
)
e
λ
x
=
0
{\displaystyle (\lambda ^{2}+a\lambda +b)e^{\lambda x}=0}
λ
2
+
a
λ
+
b
=
0
{\displaystyle \lambda ^{2}+a\lambda +b=0}
(5-5 )
Since (5-5 ) follows the same form as the quadratic equation, we can solve for
λ
1
{\displaystyle \lambda _{1}}
and
λ
2
{\displaystyle \lambda _{2}}
as follows:
λ
1
=
1
2
(
−
a
+
(
a
2
−
4
b
)
{\displaystyle \lambda _{1}={\frac {1}{2}}(-a+{\sqrt {(a^{2}-4b)}}}
λ
2
=
1
2
(
−
a
−
(
a
2
−
4
b
)
{\displaystyle \lambda _{2}={\frac {1}{2}}(-a-{\sqrt {(a^{2}-4b)}}}
Referring back to algebra, we know that the solution to these two equations can be one of three cases:
Two real roots if...
a
2
−
4
b
>
0
{\displaystyle a^{2}-4b>0}
These two roots give us two solutions:
y
1
=
e
λ
1
x
{\displaystyle y_{1}=e^{\lambda _{1}x}}
and
y
2
=
e
λ
2
x
{\displaystyle y_{2}=e^{\lambda _{2}x}}
The corresponding general solution then takes the form of the following:
y
=
c
1
e
λ
1
x
+
c
2
e
λ
2
x
{\displaystyle y=c_{1}e^{\lambda _{1}x}+c_{2}e^{\lambda _{2}x}}
(5-6 )
A real double root if...
a
2
−
4
b
=
0
{\displaystyle a^{2}-4b=0}
This yields only one solution:
y
1
=
e
λ
x
{\displaystyle y_{1}=e^{\lambda x}}
In order to form a basis, a set of two linearly independent solutions, we must find another solution. To do so, we will use the method of reduction of order, where:
y
2
=
u
y
1
{\displaystyle y_{2}=uy_{1}}
y
2
′
=
u
′
y
1
+
u
y
1
′
{\displaystyle y'_{2}=u'y_{1}+uy'_{1}}
and
y
2
″
=
u
″
y
1
+
2
u
′
y
1
′
+
u
y
1
″
{\displaystyle y''_{2}=u''y_{1}+2u'y'_{1}+uy''_{1}}
yeilding
(
u
″
y
1
+
2
u
′
y
1
′
+
u
y
1
″
)
+
a
(
u
′
y
1
+
u
y
1
′
)
+
b
u
y
1
=
0
{\displaystyle (u''y_{1}+2u'y'_{1}+uy''_{1})+a(u'y_{1}+uy'_{1})+buy_{1}=0}
u
″
y
1
+
u
′
(
2
y
1
′
+
a
y
1
)
+
u
(
y
1
″
+
a
′
y
1
+
b
y
1
)
=
0
{\displaystyle u''y_{1}+u'(2y'_{1}+ay_{1})+u(y''_{1}+a'y_{1}+by_{1})=0}
Since
y
1
{\displaystyle y_{1}}
is a solution of (5-1), the last set of parentheses is zero.
Similarly, the first parentheses is zero as well, because
2
y
1
′
=
−
a
e
a
x
/
2
=
−
a
y
1
{\displaystyle 2y'_{1}=-ae^{ax/2}=-ay_{1}}
From integration, we get the solution:
u
=
c
1
x
+
c
2
{\displaystyle u=c_{1}x+c_{2}}
If we set
c
1
=
1
{\displaystyle c_{1}=1}
and
c
2
=
0
{\displaystyle c_{2}=0}
y
2
=
x
y
1
{\displaystyle y_{2}=xy_{1}}
Thus, the general solution is:
y
=
(
c
1
x
+
c
2
)
e
λ
2
x
{\displaystyle y=(c_{1}x+c_{2})e^{\lambda _{2}x}}
(5-7 )
Complex conjugate roots if...
a
2
−
4
b
<
0
{\displaystyle a^{2}-4b<0}
In this case, the roots of (5-5 ) are complex:
y
1
=
e
−
a
x
/
2
c
o
s
(
ω
x
)
{\displaystyle y_{1}=e^{-ax/2}cos(\omega x)}
and
y
2
=
e
−
a
x
/
2
s
i
n
(
ω
x
)
{\displaystyle y_{2}=e^{-ax/2}sin(\omega x)}
Thus, the corresponding general solution is of the form:
y
=
e
−
a
x
/
2
(
A
c
o
s
(
ω
x
)
+
B
s
i
n
(
ω
x
)
)
{\displaystyle y=e^{-ax/2}(Acos(\omega x)+Bsin(\omega x))}
(5-8 )
For this problem, we were tasked with finding the general solution to the homogeneous linear ODE given below:
y
″
+
4
y
′
+
(
π
2
+
4
)
y
=
0
{\displaystyle y''+4y'+(\pi ^{2}+4)y=0}
(5a-1 )
Following the process that yields (5-5 ), we find the equation:
λ
2
e
λ
x
+
4
λ
e
λ
x
+
(
π
2
+
4
)
e
λ
x
=
0
{\displaystyle \lambda ^{2}e^{\lambda x}+4\lambda e^{\lambda x}+(\pi ^{2}+4)e^{\lambda x}=0}
(5a-2 )
To find the form of the general solution, we must solve for the values of
λ
{\displaystyle \lambda }
using the quadratic formula:
Since the solutions to
λ
{\displaystyle \lambda }
are complex numbers, our solution takes the form of (5-8 ):
y
=
e
−
2
x
(
A
c
o
s
(
π
x
)
+
B
s
i
n
(
π
x
)
)
,
{\displaystyle y=e^{-2x}(Acos(\pi x)+Bsin(\pi x)),}
(5a-3 )
A and B are constant coefficients of unknown value. Had we been given initial values of
y
(
x
)
{\displaystyle y(x)}
and
y
′
(
x
)
{\displaystyle y'(x)}
, we would solve for those values as well.
Confirmation of Solution
edit
To confirm that the solution is true, we will substitute
y
{\displaystyle y}
,
y
′
{\displaystyle y'}
, and
y
″
{\displaystyle y''}
into (5a-1 ). If the final result is a tautology[ 3] then the solution is confirmed. If the result is a contradiction, a statement that is never true, then our solution is disproved.
y
=
e
−
2
x
(
A
c
o
s
(
π
x
)
+
B
s
i
n
(
π
x
)
)
{\displaystyle y=e^{-2x}(Acos(\pi x)+Bsin(\pi x))}
y
′
=
−
2
e
−
2
x
(
A
c
o
s
(
π
x
)
+
B
s
i
n
(
π
x
)
)
+
e
−
2
x
(
−
A
π
s
i
n
(
π
x
)
+
B
π
c
o
s
(
π
x
)
)
{\displaystyle y'=-2e^{-2x}(Acos(\pi x)+Bsin(\pi x))+e^{-2x}(-A\pi sin(\pi x)+B\pi cos(\pi x))}
(5a-4 )
y
″
=
4
e
−
2
x
(
A
c
o
s
(
π
x
)
+
B
s
i
n
(
π
x
)
)
−
2
e
−
2
x
(
−
A
π
s
i
n
(
π
x
)
+
B
π
c
o
s
(
π
x
)
)
−
2
e
−
2
x
(
−
A
π
s
i
n
(
π
x
)
{\displaystyle y''=4e^{-2x}(Acos(\pi x)+Bsin(\pi x))-2e^{-2x}(-A\pi sin(\pi x)+B\pi cos(\pi x))-2e^{-2x}(-A\pi sin(\pi x)}
+
B
π
c
o
s
(
π
x
)
)
+
e
−
2
x
(
−
A
π
2
c
o
s
(
π
x
)
−
B
π
2
s
i
n
(
π
x
)
)
{\displaystyle +B\pi cos(\pi x))+e^{-2x}(-A\pi ^{2}cos(\pi x)-B\pi ^{2}sin(\pi x))}
(5a-5 )
Substituting (5a-3 ), (5a-4 ), and (5a-5 ) into (5a-1 ), we are left with:
[
4
e
−
2
x
(
A
c
o
s
(
π
x
)
+
B
s
i
n
(
π
x
)
)
−
2
e
−
2
x
(
−
A
π
s
i
n
(
π
x
)
+
B
π
c
o
s
(
π
x
)
)
−
2
e
−
2
x
(
−
A
π
s
i
n
(
π
x
)
+
B
π
c
o
s
(
π
x
)
)
+
e
−
2
x
(
−
A
π
2
c
o
s
(
π
x
)
−
B
π
2
s
i
n
(
π
x
)
)
]
{\displaystyle [4e^{-2x}(Acos(\pi x)+Bsin(\pi x))-2e^{-2x}(-A\pi sin(\pi x)+B\pi cos(\pi x))-2e^{-2x}(-A\pi sin(\pi x)+B\pi cos(\pi x))+e^{-2x}(-A\pi ^{2}cos(\pi x)-B\pi ^{2}sin(\pi x))]}
+
4
[
−
2
e
−
2
x
(
A
c
o
s
(
π
x
)
+
B
s
i
n
(
π
x
)
)
+
e
−
2
x
(
−
A
π
s
i
n
(
π
x
)
+
B
π
c
o
s
(
π
x
)
)
]
+
(
π
2
+
4
)
[
e
−
2
x
(
A
c
o
s
(
π
x
)
+
B
s
i
n
(
π
x
)
)
]
=
0
{\displaystyle +4[-2e^{-2x}(Acos(\pi x)+Bsin(\pi x))+e^{-2x}(-A\pi sin(\pi x)+B\pi cos(\pi x))]+(\pi ^{2}+4)[e^{-2x}(Acos(\pi x)+Bsin(\pi x))]=0}
Combining similar terms allows us to clean up this solution and check our answer:
e
−
2
x
[
c
o
s
(
π
x
)
[
(
4
A
−
2
B
π
−
2
B
π
−
A
π
2
)
+
(
−
8
A
+
4
B
π
)
+
(
A
π
2
+
4
A
)
]
+
s
i
n
(
π
x
)
[
(
4
B
+
2
A
π
+
2
A
π
−
B
π
2
)
+
(
−
8
B
−
4
A
π
)
+
(
B
π
2
+
4
B
)
]
]
=
0
{\displaystyle e^{-2x}[cos(\pi x)[(4A-2B\pi -2B\pi -A\pi ^{2})+(-8A+4B\pi )+(A\pi ^{2}+4A)]+sin(\pi x)[(4B+2A\pi +2A\pi -B\pi ^{2})+(-8B-4A\pi )+(B\pi ^{2}+4B)]]=0}
e
−
2
x
[
c
o
s
(
π
x
)
[
A
(
4
+
4
−
8
−
π
2
+
π
2
)
+
B
(
−
2
π
−
2
π
+
4
π
)
]
+
s
i
n
(
π
x
)
[
A
(
2
π
+
2
π
−
4
π
)
+
B
(
4
+
4
−
8
−
π
2
+
π
2
)
]
]
=
0
{\displaystyle e^{-2x}[cos(\pi x)[A(4+4-8-\pi ^{2}+\pi ^{2})+B(-2\pi -2\pi +4\pi )]+sin(\pi x)[A(2\pi +2\pi -4\pi )+B(4+4-8-\pi ^{2}+\pi ^{2})]]=0}
e
−
2
x
[
c
o
s
(
π
x
)
(
0
)
+
s
i
n
(
π
x
)
(
0
)
]
=
0
{\displaystyle e^{-2x}[cos(\pi x)(0)+sin(\pi x)(0)]=0}
0
=
0
{\displaystyle 0=0}
Since the outcome of our check is a tautology, the solution is confirmed for all values of A and B and we affirm that the final solution is:
y
=
e
−
2
x
(
A
c
o
s
(
π
x
)
+
B
s
i
n
(
π
x
)
)
{\displaystyle y=e^{-2x}(Acos(\pi x)+Bsin(\pi x))}
For this problem, we were tasked with finding the general solution to the homogeneous linear ODE given below:
y
″
+
2
π
y
′
+
π
2
y
=
0
{\displaystyle {y}''+2\pi {y}'+\pi ^{2}y=0}
(5b-1 )
Following the process that yields (5-5 ), we find the equation:
λ
2
e
λ
x
+
2
π
λ
e
λ
x
+
π
2
e
λ
x
=
0
{\displaystyle \lambda ^{2}e^{\lambda x}+2\pi \lambda e^{\lambda x}+\pi ^{2}e^{\lambda x}=0}
(5b-2 )
To find the form of the general solution, we must solve for the values of
λ
{\displaystyle \lambda }
. Instead of going through the quadratic formula, we can analyze the discriminant.
(
2
π
)
2
−
4
π
2
=
0
{\displaystyle {\sqrt {(2\pi )^{2}-4\pi ^{2}}}=0}
Thus we see that the discriminant follows Case 2 and the general solution must follow the form of (5-7 ), and we have the following:
y
=
(
c
1
+
c
2
x
)
e
−
π
x
{\displaystyle y=(c_{1}+c_{2}x)e^{-\pi x}}
(5b-3 )
where
c
1
{\displaystyle c_{1}}
and
c
2
{\displaystyle c_{2}}
are unknown constant coefficients. If this were an initial value problem with values of
y
(
x
)
{\displaystyle y(x)}
and
y
′
(
x
)
{\displaystyle y'(x)}
, we would solve for those values as well.
Confirmation of Solution
edit
To confirm that the solution is true, we will substitute
y
{\displaystyle y}
,
y
′
{\displaystyle y'}
, and
y
″
{\displaystyle y''}
into (5a-1 ) and check to see if the final result is a tautology or a contradiction.
y
′
=
−
π
(
c
1
+
c
2
x
)
e
−
π
x
+
c
2
e
−
π
x
{\displaystyle y'=-\pi (c_{1}+c_{2}x)e^{-\pi x}+c_{2}e^{-\pi x}}
(5b-4 )
y
″
=
π
2
(
c
1
+
c
2
x
)
e
−
π
x
+
−
2
π
c
2
e
−
π
x
{\displaystyle y''=\pi ^{2}(c_{1}+c_{2}x)e^{-\pi x}+-2\pi c_{2}e^{-\pi x}}
(5b-5 )
Substituting (5b-3 ), (5b-4 ), and (5b-5 ) into (5b-1 ), we are left with:
[
π
2
(
c
1
+
c
2
x
)
e
−
π
x
+
−
2
π
c
2
e
−
π
x
]
+
2
π
[
−
π
(
c
1
+
c
2
x
)
e
−
π
x
+
c
2
e
−
π
x
]
+
π
2
[
(
c
1
+
c
2
x
)
e
−
π
x
]
=
0
{\displaystyle [\pi ^{2}(c_{1}+c_{2}x)e^{-\pi x}+-2\pi c_{2}e^{-\pi x}]+2\pi [-\pi (c_{1}+c_{2}x)e^{-\pi x}+c_{2}e^{-\pi x}]+\pi ^{2}[(c_{1}+c_{2}x)e^{-\pi x}]=0}
After combining similar terms, we have the form:
π
2
[
(
c
1
+
c
2
x
)
e
−
π
x
−
2
(
c
1
+
c
2
x
)
e
−
π
x
+
(
c
1
+
c
2
x
)
e
−
π
x
]
+
2
π
[
c
2
e
−
π
x
−
c
2
e
−
π
x
]
=
0
{\displaystyle \pi ^{2}[(c_{1}+c_{2}x)e^{-\pi x}-2(c_{1}+c_{2}x)e^{-\pi x}+(c_{1}+c_{2}x)e^{-\pi x}]+2\pi [c_{2}e^{-\pi x}-c_{2}e^{-\pi x}]=0}
0
=
0
{\displaystyle 0=0}
Since the above result is a tautology for all values of
c
1
{\displaystyle c_{1}}
and
c
2
{\displaystyle c_{2}}
, (5-3 ) is confirmed and the general solution to (5-1 ) is:
y
=
(
c
1
+
c
2
x
)
e
−
π
x
{\displaystyle y=(c_{1}+c_{2}x)e^{-\pi x}}
↑ EGM4313: Section 2 Notes (.djvu)
↑ Kreyszig, Erwin; Herbert Kreyszig; Edward J. Norminton (2011). Advanced Engineering Mathematics . John Wily & Sons, Inc . ISBN 978-0-470-45836-5 . Pg 52-58
↑ Introduction to Propositional Logic
R.1.6: Determination of ODE Order, Linearity, and Application of the Superposition Principle[ 1]
edit
Given the eight Ordinary Differential Equations shown below:
(1) Determine the order
(2) Determine the linearity (or lack of)
(3) Show whether the principle of superposition can be applied for each equation
Consider a function
y
¯
(
x
)
{\displaystyle {\bar {y}}(x)}
, defined as the sum of the homogeneous solution,
y
h
(
x
)
{\displaystyle y_{h}(x)}
, and the particular solution,
y
p
(
x
)
{\displaystyle y_{p}(x)}
:
y
¯
(
x
)
=
y
p
(
x
)
+
y
h
(
x
)
{\displaystyle {\bar {y}}(x)=y_{p}(x)+y_{h}(x)}
(6-1 )
Example:
Consider the Differential Equation in standard form:
y
″
+
p
(
x
)
y
′
+
q
(
x
)
y
=
r
(
x
)
{\displaystyle y''+p(x)y'+q(x)y=r(x)}
(6-2 )
The homogeneous solution is then:
y
h
″
+
p
(
x
)
y
h
′
+
q
(
x
)
y
h
=
0
{\displaystyle y_{h}''+p(x)y_{h}'+q(x)y_{h}=0}
(6-3 )
and the particular solution is:
y
p
″
+
p
(
x
)
y
p
′
+
q
(
x
)
y
p
=
r
(
x
)
{\displaystyle y_{p}''+p(x)y_{p}'+q(x)y_{p}=r(x)}
(6-4 )
Now we sum up equations (6-3 ) and (6-4 ) to get:
(
y
h
″
+
y
p
″
)
+
p
(
x
)
(
y
h
′
+
y
p
′
)
+
q
(
x
)
(
y
h
+
y
p
)
=
r
(
x
)
{\displaystyle (y_{h}''+y_{p}'')+p(x)(y_{h}'+y_{p}')+q(x)(y_{h}+y_{p})=r(x)}
(6-5 )
If the Equation is Linear then some reductions can take place:
y
h
″
+
y
p
″
=
(
y
h
+
y
p
)
″
=
y
¯
″
{\displaystyle y_{h}''+y_{p}''=(y_{h}+y_{p})''={\bar {y}}''}
and
y
h
′
+
y
p
′
=
(
y
h
+
y
p
)
′
=
y
¯
′
{\displaystyle y_{h}'+y_{p}'=(y_{h}+y_{p})'={\bar {y}}'}
Now we can make some substitutions into (6-5 )using (6-3 )and (6-4 ):
y
¯
″
+
p
(
x
)
y
¯
′
+
q
(
x
)
y
¯
=
r
(
x
)
{\displaystyle {\bar {y}}''+p(x){\bar {y}}'+q(x){\bar {y}}=r(x)}
(6-6 )
Which is the same as equation (6-2 ):
y
″
+
p
(
x
)
y
′
+
q
(
x
)
y
=
r
(
x
)
{\displaystyle y''+p(x)y'+q(x)y=r(x)}
y
″
=
g
=
c
o
n
s
t
a
n
t
{\displaystyle y''=g=constant}
(6a-1 )
Order:
2
n
d
{\displaystyle 2^{nd}}
Order ODE Linearity: Linear
Superposition: Yes
This can be shown using the method from the background theory.
Recall:
y
¯
(
x
)
=
y
p
(
x
)
+
y
h
(
x
)
{\displaystyle {\bar {y}}(x)=y_{p}(x)+y_{h}(x)}
The homogeneous solution is:
y
h
″
=
0
{\displaystyle y_{h}''=0}
(6a-2 )
The particular solution is:
y
p
″
=
g
{\displaystyle y_{p}''=g}
(6a-3 )
Now we add Equations (6a-2 ) and (6a-3 ):
y
h
″
+
y
p
″
=
g
{\displaystyle y_{h}''+y_{p}''=g}
(6a-4 )
We can use (6-1 ) to simplify (6a-4 ):
(
y
h
+
y
p
)
″
=
g
{\displaystyle (y_{h}+y_{p})''=g}
(6a-5 )
y
¯
″
=
g
{\displaystyle {\bar {y}}''=g}
(6a-6 )
Since (6a-1 ) and (6a-6 ) are the same then we can apply Superposition.
m
v
′
=
m
g
−
b
v
2
{\displaystyle mv'=mg-bv^{2}}
(6b-1 )
Order:
1
s
t
{\displaystyle 1^{st}}
Order ODE Linearity: Non-Linear
Superposition: No
This can be shown using the method from the background theory.
Before we start it is easier to do the problem if the equation is rearranged as such:
m
v
′
+
b
v
2
=
m
g
{\displaystyle mv'+bv^{2}=mg}
(6b-2 )
Recall:
v
¯
(
x
)
=
v
p
(
x
)
+
v
h
(
x
)
{\displaystyle {\bar {v}}(x)=v_{p}(x)+v_{h}(x)}
The homogeneous solution is:
m
v
h
′
+
m
v
h
2
=
0
{\displaystyle mv_{h}'+mv_{h}^{2}=0}
(6b-3 )
The particular solution is:
m
v
p
′
+
m
v
p
2
=
m
g
{\displaystyle mv_{p}'+mv_{p}^{2}=mg}
(6b-4 )
Now we add Equations (6b-3 ) and (6b-4 ):
m
v
h
′
+
m
v
h
2
+
m
v
p
′
+
m
v
p
2
=
m
g
{\displaystyle mv_{h}'+mv_{h}^{2}+mv_{p}'+mv_{p}^{2}=mg}
(6b-5 )
We can use (6-1 ) to simplify (6b-5 ):
m
(
y
h
+
y
p
)
′
+
m
(
v
h
2
+
v
p
2
)
=
m
g
{\displaystyle m(y_{h}+y_{p})'+m(v_{h}^{2}+v_{p}^{2})=mg}
(6b-6 )
Note that:
(
v
h
2
+
v
p
2
)
≠
v
¯
2
{\displaystyle (v_{h}^{2}+v_{p}^{2})\neq {\bar {v}}^{2}}
So we cannot apply Super position.
h
′
=
−
k
h
{\displaystyle {h}'=-k{\sqrt {h}}}
(6c-1 )
Order:
1
s
t
{\displaystyle 1^{st}}
Order ODE Linearity: Non-Linear
Superposition: No
This can be shown using the method from the background theory.
Before we start it is easier to do the problem if the equation is rearranged like this:
h
′
+
k
h
=
0
{\displaystyle {h}'+k{\sqrt {h}}=0}
(6c-2 )
Recall:
h
¯
(
x
)
=
h
h
(
x
)
+
h
p
(
x
)
{\displaystyle {\bar {h}}(x)=h_{h}(x)+h_{p}(x)}
The homogeneous solution is:
h
h
′
+
k
h
h
=
0
{\displaystyle {h}_{h}'+k{\sqrt {h_{h}}}=0}
(6c-3 )
The particular solution is:
h
p
′
+
k
h
p
=
0
{\displaystyle {h_{p}}'+k{\sqrt {h_{p}}}=0}
(6c-4 )
Now we add Equations (6c-3 ) and (6c-4 ):
h
h
′
+
k
h
h
+
h
p
′
+
k
h
p
=
0
{\displaystyle {h_{h}}'+k{\sqrt {h_{h}}}+{h_{p}}'+k{\sqrt {h_{p}}}=0}
(6c-5 )
We can use (6-1 ) to simplify (6c-5 ):
(
h
h
+
h
p
)
′
+
k
(
h
h
+
h
p
)
=
0
{\displaystyle ({h_{h}}+{h_{p}})'+k({\sqrt {h_{h}}}+{\sqrt {h_{p}}})=0}
(6c-6 )
Note that:
(
h
h
+
h
p
)
≠
h
¯
{\displaystyle ({\sqrt {h_{h}}}+{\sqrt {h_{p}}})\neq {\sqrt {\bar {h}}}}
So we cannot apply Super position.
m
y
″
+
k
y
=
0
{\displaystyle m{y}''+ky=0}
(6d-1 )
Order:
2
n
d
{\displaystyle 2^{nd}}
Order ODE Linearity: Linear
Superposition: Yes
This can be shown using the method from the background theory.
Recall:
y
¯
(
x
)
=
y
p
(
x
)
+
y
h
(
x
)
{\displaystyle {\bar {y}}(x)=y_{p}(x)+y_{h}(x)}
The homogeneous solution is:
m
y
h
″
+
k
y
h
=
0
{\displaystyle m{y}_{h}''+ky_{h}=0}
(6d-2 )
The particular solution is:
m
y
p
″
+
k
y
p
=
0
{\displaystyle m{y}_{p}''+ky_{p}=0}
(6d-3 )
Now we add Equations (6d-2 ) and (6d-3 ):
m
y
h
″
+
k
y
h
+
m
y
p
″
+
k
y
p
=
0
{\displaystyle m{y}_{h}''+ky_{h}+m{y}_{p}''+ky_{p}=0}
(6d-4 )
We can use (6-1 ) to simplify (6d-4 ):
m
(
y
h
+
y
p
)
″
+
k
(
y
h
+
y
p
)
=
0
{\displaystyle m({y}_{h}+{y}_{p})''+k(y_{h}+y_{p})=0}
(6d-5 )
m
y
¯
″
+
k
y
¯
=
0
{\displaystyle m{\bar {y}}''+k{\bar {y}}=0}
(6d-6 )
Since (6d-1 ) and (6d-6 ) are the same then we can apply Superposition.
y
″
+
ω
0
2
y
=
cos
ω
t
,
ω
0
≈
ω
{\displaystyle {y}''+\omega _{0}^{2}y=\cos \omega t,\omega _{0}\approx \omega }
(6e-1 )
Order:
2
n
d
{\displaystyle 2^{nd}}
Order ODE Linearity: Linear
Superposition: Yes
This can be shown using the method from the background theory.
Recall:
y
¯
(
x
)
=
y
p
(
x
)
+
y
h
(
x
)
{\displaystyle {\bar {y}}(x)=y_{p}(x)+y_{h}(x)}
The homogeneous solution is:
y
h
″
+
ω
0
2
y
h
=
0
{\displaystyle {y}_{h}''+\omega _{0}^{2}y_{h}=0}
(6e-2 )
The particular solution is:
y
p
″
+
ω
0
2
y
p
=
cos
ω
t
{\displaystyle {y}_{p}''+\omega _{0}^{2}y_{p}=\cos \omega t}
(6e-3 )
Now we add Equations (6e-2 ) and (6e-3 ):
y
h
″
+
ω
0
2
y
h
+
y
p
″
+
ω
0
2
y
p
=
cos
ω
t
{\displaystyle {y}_{h}''+\omega _{0}^{2}y_{h}+{y}_{p}''+\omega _{0}^{2}y_{p}=\cos \omega t}
(6e-4 )
We can use (6-1 ) to simplify (6e-4 ):
(
y
h
+
y
p
)
″
+
ω
0
2
(
y
h
+
y
p
)
=
cos
ω
t
{\displaystyle ({y}_{h}+{y}_{p})''+\omega _{0}^{2}(y_{h}+y_{p})=\cos \omega t}
(6e-5 )
y
¯
″
+
ω
0
2
y
¯
=
cos
ω
t
{\displaystyle {\bar {y}}''+\omega _{0}^{2}{\bar {y}}=\cos \omega t}
(6e-6 )
Since (6e-1 ) and (6e-6 ) are the same then we can apply Superposition.
L
I
″
+
R
I
′
+
1
C
I
=
E
′
{\displaystyle LI''+RI'+{\frac {1}{C}}I=E'}
(6f-1 )
Order:
2
n
d
{\displaystyle 2^{nd}}
Order ODE Linearity: Linear
Superposition: Yes
This can be shown using the method from the background theory.
Recall:
I
¯
(
x
)
=
I
p
(
x
)
+
I
h
(
x
)
{\displaystyle {\bar {I}}(x)=I_{p}(x)+I_{h}(x)}
The homogeneous solution is:
L
I
h
″
+
R
I
h
′
+
1
C
I
h
=
0
{\displaystyle LI_{h}''+RI_{h}'+{\frac {1}{C}}I_{h}=0}
(6f-2 )
The particular solution is:
L
I
p
″
+
R
I
p
′
+
1
C
I
p
=
E
′
{\displaystyle LI_{p}''+RI_{p}'+{\frac {1}{C}}I_{p}=E'}
(6f-3 )
Now we add Equations (6f-2 ) and (6f-3 ):
L
I
h
″
+
R
I
h
′
+
1
C
I
h
+
L
I
p
″
+
R
I
p
′
+
1
C
I
p
=
E
′
{\displaystyle LI_{h}''+RI_{h}'+{\frac {1}{C}}I_{h}+LI_{p}''+RI_{p}'+{\frac {1}{C}}I_{p}=E'}
(6f-4 )
We can use (6-1 ) to simplify (6f-4 ):
L
(
I
h
+
I
p
)
″
+
R
(
I
h
+
I
p
)
′
+
1
C
(
I
h
+
I
p
)
=
E
′
{\displaystyle L(I_{h}+I_{p})''+R(I_{h}+I_{p})'+{\frac {1}{C}}(I_{h}+I_{p})=E'}
(6f-5 )
L
I
¯
″
+
R
I
¯
′
+
1
C
I
¯
=
E
′
{\displaystyle L{\bar {I}}''+R{\bar {I}}'+{\frac {1}{C}}{\bar {I}}=E'}
(6f-6 )
Since (6f-1 ) and (6f-6 ) are the same then we can apply Superposition.
E
I
y
i
v
=
f
(
x
)
{\displaystyle EIy^{iv}=f(x)}
(6g-1 )
Order:
4
t
h
{\displaystyle 4^{th}}
Order ODE Linearity: Linear
Superposition: Yes
This can be shown using the method from the background theory.
Recall:
y
¯
(
x
)
=
y
p
(
x
)
+
y
h
(
x
)
{\displaystyle {\bar {y}}(x)=y_{p}(x)+y_{h}(x)}
The homogeneous solution is:
E
I
y
h
i
v
=
0
{\displaystyle EIy_{h}^{iv}=0}
(6g-2 )
The particular solution is:
E
I
y
p
i
v
=
f
(
x
)
{\displaystyle EIy_{p}^{iv}=f(x)}
(6g-3 )
Now we add Equations (6g-2 ) and (6g-3 ):
E
I
y
h
i
v
+
E
I
y
p
i
v
=
f
(
x
)
{\displaystyle EIy_{h}^{iv}+EIy_{p}^{iv}=f(x)}
(6g-4 )
We can use (6-1 ) to simplify (6g-4 ):
E
I
(
y
h
+
y
p
)
i
v
=
f
(
x
)
{\displaystyle EI(y_{h}+y_{p})^{iv}=f(x)}
(6g-5 )
E
I
(
y
¯
)
i
v
=
f
(
x
)
{\displaystyle EI({\bar {y}})^{iv}=f(x)}
(6g-6 )
Since (6g-1 ) and (6g-6 ) are the same then we can apply Superposition.
L
θ
″
+
g
sin
θ
=
0
{\displaystyle L\theta ''+g\sin \theta =0}
(6h-1 )
Order:
2
n
d
{\displaystyle 2^{nd}}
Order ODE Linearity: Non-Linear
Superposition: No
This can be shown using the method from the background theory.
Recall:
θ
¯
(
x
)
=
θ
p
(
x
)
+
θ
h
(
x
)
{\displaystyle {\bar {\theta }}(x)=\theta _{p}(x)+\theta _{h}(x)}
The homogeneous solution is:
L
θ
h
″
+
g
sin
θ
h
=
0
{\displaystyle L\theta _{h}''+g\sin \theta _{h}=0}
(6h-2 )
The particular solution is:
L
θ
p
″
+
g
sin
θ
p
=
0
{\displaystyle L\theta _{p}''+g\sin \theta _{p}=0}
(6h-2 )
Now we add Equations (6h-2 ) and (6h-3 ):
L
θ
h
″
+
g
sin
θ
h
+
L
θ
p
″
+
g
sin
θ
p
=
0
{\displaystyle L\theta _{h}''+g\sin \theta _{h}+L\theta _{p}''+g\sin \theta _{p}=0}
(6h-4 )
We can use (6-1 ) to simplify (6h-4 ):
L
(
θ
h
+
θ
p
)
″
+
g
(
sin
θ
h
+
sin
θ
p
)
=
0
{\displaystyle L(\theta _{h}+\theta _{p})''+g(\sin \theta _{h}+\sin \theta _{p})=0}
(6h-5 )
Note that:
(
sin
θ
h
+
sin
θ
p
)
≠
sin
θ
¯
{\displaystyle (\sin \theta _{h}+\sin \theta _{p})\neq \sin {\bar {\theta }}}
So we cannot apply Super position.
↑ 1.0 1.1 EGM4313: Section 2 Notes (.djvu)
↑ Kreyszig, Erwin; Herbert Kreyszig; Edward J. Norminton (2011). Advanced Engineering Mathematics . John Wily & Sons, Inc . ISBN 978-0-470-45836-5 . Pg 46-48
Contributing Members
edit
Team Contribution Table
Problem Number
Solved and Typed By
Proofread By
1.1
Maxwell Shuman
Yamil Herrera
1.2
Dalwin Marrero
Maxwell Shuman
1.3
Ron D'Amico
Dalwin Marrero
1.4
Jennifer Melroy
Ron D'Amico
1.5
Mark James and Avery Cornell
Jennifer Melroy
1.6
Yamil Herrera
Avery Cornell
Diagrams and Wiki Set-up
Mark James
Team