Skip to main content

Calculus III - 3-Dimensional Space: Equations of Lines

In this section we need to take a look at the equation of a line in  R 3 R 3 . As we saw in the previous section the equation  y = m x + b y = m x + b  does not describe a line in  R 3 R 3 , instead it describes a plane. This doesn’t mean however that we can’t write down an equation for a line in 3-D space. We’re just going to need a new way of writing down the equation of a curve. So, before we get into the equations of lines we first need to briefly look at vector functions. We’re going to take a more in depth look at vector functions later. At this point all that we need to worry about is notational issues and how they can be used to give the equation of a curve. The best way to get an idea of what a vector function is and what its graph looks like is to look at an example. So, consider the following vector function. → r ( t ) = ⟨ t , 1 ⟩ r → ( t ) = ⟨ t , 1 ⟩ A vector function is a function that takes one or more variables, one in this case, and returns a vector. Note as we

Differential Equations - Second Order: Mechanical Vibrations - i

It’s now time to take a look at an application of second order differential equations. We’re going to take a look at mechanical vibrations. In particular we are going to look at a mass that is hanging from a spring.

Vibrations can occur in pretty much all branches of engineering and so what we’re going to be doing here can be easily adapted to other situations, usually with just a change in notation.

Let’s get the situation setup. We are going to start with a spring of length l, called the natural length, and we’re going to hook an object with mass m up to it. When the object is attached to the spring the spring will stretch a length of L. We will call the equilibrium position the position of the center of gravity for the object as it hangs on the spring with no movement.
Below is sketch of the spring with and without the object attached to it.

As denoted in the sketch we are going to assume that all forces, velocities, and displacements in the downward direction will be positive. All forces, velocities, and displacements in the upward direction will be negative.

Also, as shown in the sketch above, we will measure all displacement of the mass from its equilibrium position. Therefore, the u=0 position will correspond to the center of gravity for the mass as it hangs on the spring and is at rest (i.e. no movement).

Now, we need to develop a differential equation that will give the displacement of the object at any time t. First, recall Newton’s Second Law of Motion.
ma=F
In this case we will use the second derivative of the displacement, u, for the acceleration and so Newton<’s Second Law becomes,
mu=F(t,u,u)
We now need to determine all the forces that will act upon the object. There are four forces that we will assume act upon the object. Two that will always act on the object and two that may or may not act upon the object.
Here is a list of the forces that will act upon the object.
  1. Gravity, FgThe force due to gravity will always act upon the object of course. This force is
    Fg=mg
  2. Spring, FsWe are going to assume that Hooke’s Law will govern the force that the spring exerts on the object. This force will always be present as well and is
    Fs=k(L+u)Hooke’s Law tells us that the force exerted by a spring will be the spring constant, k>0, times the displacement of the spring from its natural length. For our set up the displacement from the spring’s natural length is L+u and the minus sign is in there to make sure that the force always has the correct direction.
    Let’s make sure that this force does what we expect it to. If the object is at rest in its equilibrium position the displacement is L and the force is simply Fs=kL which will act in the upward position as it should since the spring has been stretched from its natural length.
    If the spring has been stretched further down from the equilibrium position then L+u will be positive and Fs will be negative acting to pull the object back up as it should be.
    Next, if the object has been moved up past its equilibrium point, but not yet to its natural length then u will be negative, but still less than L and so L+u will be positive and once again Fs will be negative acting to pull the object up.
    Finally, if the object has been moved upwards so that the spring is now compressed, then u will be negative and greater than L. Therefore, L+u will be negative and now Fs will be positive acting to push the object down.
    So, it looks like this force will act as we expect that it should.
  3. Damping, FdThe next force that we need to consider is damping. This force may or may not be present for any given problem.
    Dampers work to counteract any movement. There are several ways to define a damping force. The one that we’ll use is the following.
    Fd=γuwhere, γ>0 is the damping coefficient. Let’s think for a minute about how this force will act. If the object is moving downward, then the velocity (u) will be positive and so Fd will be negative and acting to pull the object back up. Likewise, if the object is moving upward, the velocity (u) will be negative and so Fd will be positive and acting to push the object back down.
    In other words, the damping force as we’ve defined it will always act to counter the current motion of the object and so will act to damp out any motion in the object.
  4. External Forces, F(t)This is the catch all force. If there are any other forces that we decide we want to act on our object we lump them in here and call it good. We typically call F(t) the forcing function.
Putting all of these together gives us the following for Newton’s Second Law.
mu=mgk(L+u)γu+F(t)
Or, upon rewriting, we get,
mu+γu+ku=mgkL+F(t)
Now, when the object is at rest in its equilibrium position there are exactly two forces acting on the object, the force due to gravity and the force due to the spring. Also, since the object is at rest (i.e. not moving) these two forces must be canceling each other out. This means that we must have,
(1)mg=kL
Using this in Newton’s Second Law gives us the final version of the differential equation that we’ll work with.
(2)mu+γu+ku=F(t)
Along with this differential equation we will have the following initial conditions.
(3)u(0)=u0Initial displacement from the equilibrium position.u(0)=u0Initial velocity.
Note that we’ll also be using (1) to determine the spring constant, k.
Okay. Let’s start looking at some specific cases.

Free, Undamped Vibrations

This is the simplest case that we can consider. Free or unforced vibrations means that F(t)=0 and undamped vibrations means that γ=0. In this case the differential equation becomes,
mu+ku=0
This is easy enough to solve in general. The characteristic equation has the roots,
r=±ikm
This is usually reduced to,
r=±ω0i
where,
ω0=km
and ω0 is called the natural frequency. Recall as well that m>0 and k>0 and so we can guarantee that this quantity will not be complex. The solution in this case is then
(4)u(t)=c1cos(ω0t)+c2sin(ω0t)
We can write (4) in the following form,
(5)u(t)=Rcos(ω0tδ)
where R is the amplitude of the displacement and δ is the phase shift or phase angle of the displacement.
When the displacement is in the form of (5) it is usually easier to work with. However, it’s easier to find the constants in (4) from the initial conditions than it is to find the amplitude and phase shift in (5) from the initial conditions. So, in order to get the equation into the form in (5)we will first put the equation in the form in (4), find the constants, c1 and c2 and then convert this into the form in (5).
So, assuming that we have c1 and c2 how do we determine R and δ? Let’s start with (5) and use a trig identity to write it as
(6)u(t)=Rcos(δ)cos(ω0t)+Rsin(δ)sin(ω0t)
Now, R and δ are constants and so if we compare (6) to (4) we can see that
c1=Rcosδc2=Rsinδ
We can find R in the following way.
c12+c22=R2cos2δ+R2sin2δ=R2
Taking the square root of both sides and assuming that R is positive will give
(7)R=c12+c22
Finding δ is just as easy. We’ll start with
c2c1=RsinδRcosδ=tanδ
Taking the inverse tangent of both sides gives,
(8)δ=tan1(c2c1)
Before we work any examples let’s talk a little bit about units of mass and the Imperial vs. metric system differences.
Recall that the weight of the object is given by
W=mg
where m is the mass of the object and g is the gravitational acceleration. For the examples in this problem we’ll be using the following values for g.
Imperial : g=32ft/s2Metric : g=9.8m/s2
This is not the standard 32.2 ft/s2 or 9.81 m/s2, but using these will make some of the numbers come out a little nicer.
In the metric system the mass of objects is given in kilograms (kg) and there is nothing for us to do. However, in the British system we tend to be given the weight of an object in pounds (yes, pounds are the units of weight not mass…) and so we’ll need to compute the mass for these problems.
At this point we should probably work an example of all this to see how this stuff works.


Example 1 A 16 lb object stretches a spring 89 ft by itself. There is no damping and no external forces acting on the system. The spring is initially displaced 6 inches upwards from its equilibrium position and given an initial velocity of 1 ft/sec downward. Find the displacement at any time tu(t).

We first need to set up the IVP for the problem. This requires us to get our hands on m and k.
This is the Imperial system so we’ll need to compute the mass.
m=Wg=1632=12Now, let’s get k. We can use the fact that mg=kL to find k. Don’t forget that we’ll need all of our length units the same. We’ll use feet for the unit of measurement for this problem.
k=mgL=168/9=18We can now set up the IVP.
12u+18u=0u(0)=12u(0)=1For the initial conditions recall that upward displacement/motion is negative while downward displacement/motion is positive. Also, since we decided to do everything in feet we had to convert the initial displacement to feet.
Now, to solve this we can either go through the characteristic equation or we can just jump straight to the formula that we derived above. We’ll do it that way. First, we need the natural frequency,
ω0=181/2=36=6The general solution, along with its derivative, is then,
u(t)=c1cos(6t)+c2sin(6t)u(t)=6c1sin(6t)+6c2cos(6t)Applying the initial conditions gives
12=u(0)=c1c1=121=u(0)=6c2cos(6t)c2=16The displacement at any time t is then
u(t)=12cos(6t)+16sin(6t)Now, let’s convert this to a single cosine. First let’s get the amplitude, R.
R=(12)2+(16)2=106=0.52705You can use either the exact value here or a decimal approximation. Often the decimal approximation will be easier.
Now let’s get the phase shift.
δ=tan1(1/61/2)=0.32175We need to be careful with this part. The phase angle found above is in Quadrant IV, but there is also an angle in Quadrant II that would work as well. We get this second angle by adding π onto the first angle. So, we actually have two angles. They are
δ1=0.32175δ2=δ1+π=2.81984We need to decide which of these phase shifts is correct, because only one will be correct. To do this recall that
c1=Rcosδc2=RsinδNow, since we are assuming that R is positive this means that the sign of cosδ will be the same as the sign of c1 and the sign of sinδ will be the same as the sign of c2. So, for this particular case we must have cosδ<0 and sinδ>0. This means that the phase shift must be in Quadrant II and so the second angle is the one that we need.
So, after all of this the displacement at any time t is.
u(t)=0.52705cos(6t2.81984)Here is a sketch of the displacement for the first 5 seconds.

Now, let’s take a look at a slightly more realistic situation. No vibration will go on forever. So, let’s add in a damper and see what happens now.


Free, Damped Vibrations

We are still going to assume that there will be no external forces acting on the system, with the exception of damping of course. In this case the differential equation will be.
mu+γu+ku=0
where mδ, and k are all positive constants. Upon solving for the roots of the characteristic equation we get the following.
r1,2=γ±γ24mk2m
We will have three cases here.
  1. γ24mk=0In this case we will get a double root out of the characteristic equation and the displacement at any time t will be.
    u(t)=c1eγt2m+c2teγt2mNotice that as t the displacement will approach zero and so the damping in this case will do what it’s supposed to do.

    This case is called critical damping and will happen when the damping coefficient is,
    γ24mk=0γ2=4mkγ=2mk=γCRThe value of the damping coefficient that gives critical damping is called the critical damping coefficient and denoted by γCR.
  2. γ24mk>0In this case let’s rewrite the roots a little.
    r1,2=γ±γ24mk2m=γ±γ14mkγ22m=γ2m(1±14mkγ2)Also notice that from our initial assumption that we have,
    γ2>4mk1>4mkγ2Using this we can see that the fraction under the square root above is less than one. Then if the quantity under the square root is less than one, this means that the square root of this quantity is also going to be less than one. In other words,
    14mkγ2<1Why is this important? Well, the quantity in the parenthesis is now one plus/minus a number that is less than one. This means that the quantity in the parenthesis is guaranteed to be positive and so the two roots in this case are guaranteed to be negative. Therefore, the displacement at any time t is,
    u(t)=c1er1t+c2er2tand will approach zero as t. So, once again the damper does what it is supposed to do.
    This case will occur when
    γ2>4mkγ>2mkγ>γCRand is called over damping.
  3. γ24mk<0In this case we will get complex roots out of the characteristic equation.
    r1,2=γ2m±γ24mk2m=λ±μiwhere the real part is guaranteed to be negative and so the displacement is
    u(t)=c1eλtcos(μt)+c2eλtsin(μt)=eλt(c1cos(μt)+c2sin(μt))=Reλtcos(μtδ)Notice that we reduced the sine and cosine down to a single cosine in this case as we did in the undamped case. Also, since λ<0 the displacement will approach zero as t and the damper will also work as it’s supposed to in this case.
    We will get this case will occur when
    γ2<4mkγ<2mkγ<γCRand is called under damping.
Let’s take a look at a couple of examples here with damping.


Example 2 Take the spring and mass system from the first example and attach a damper to it that will exert a force of 12 lbs when the velocity is 2 ft/s. Find the displacement at any time tu(t).

The mass and spring constant were already found in the first example so we won’t do the work here. We do need to find the damping coefficient however. To do this we will use the formula for the damping force given above with one modification. The original damping force formula is,
Fd=γuHowever, remember that the force and the velocity are always acting in opposite directions. So, if the velocity is upward (i.e. negative) the force will be downward (i.e. positive) and so the minus in the formula will cancel against the minus in the velocity. Likewise, if the velocity is downward (i.e. positive) the force will be upwards (i.e. negative) and in this case the minus sign in the formula will cancel against the minus in the force. In other words, we can drop the minus sign in the formula and use
Fd=γuand then just ignore any signs for the force and velocity.
Doing this gives us the following for the damping coefficient
12=γ(2)γ=6The IVP for this example is then,
12u+6u+18u=0u(0)=12u(0)=1Before solving let’s check to see what kind of damping we’ve got. To do this all we need is the critical damping coefficient.
γCR=2km=2(18)(12)=29=6So, it looks like we’ve got critical damping. Note that this means that when we go to solve the differential equation we should get a double root.
Speaking of solving, let’s do that. I’ll leave the details to you to check that the displacement at any time t is.
u(t)=12e6t2te6tHere is a sketch of the displacement during the first 3 seconds.

Notice that the “vibration” in the system is not really a true vibration as we tend to think of them. In the critical damping case there isn’t going to be a real oscillation about the equilibrium point that we tend to associate with vibrations. The damping in this system is strong enough to force the “vibration” to die out before it ever really gets a chance to do much in the way of oscillation.

Comments

Popular posts from this blog

Digital Signal Processing - Basic Continuous Time Signals

To test a system, generally, standard or basic signals are used. These signals are the basic building blocks for many complex signals. Hence, they play a very important role in the study of signals and systems. Unit Impulse or Delta Function A signal, which satisfies the condition,   δ ( t ) = lim ϵ → ∞ x ( t ) δ ( t ) = lim ϵ → ∞ x ( t )   is known as unit impulse signal. This signal tends to infinity when t = 0 and tends to zero when t ≠ 0 such that the area under its curve is always equals to one. The delta function has zero amplitude everywhere except at t = 0. Properties of Unit Impulse Signal δ(t) is an even signal. δ(t) is an example of neither energy nor power (NENP) signal. Area of unit impulse signal can be written as; A = ∫ ∞ − ∞ δ ( t ) d t = ∫ ∞ − ∞ lim ϵ → 0 x ( t ) d t = lim ϵ → 0 ∫ ∞ − ∞ [ x ( t ) d t ] = 1 Weight or strength of the signal can be written as; y ( t ) = A δ ( t ) y ( t ) = A δ ( t ) Area of the weighted impulse signal can

Differential Equations - First Order: Bernoulli

In this section we are going to take a look at differential equations in the form, y ′ + p ( x ) y = q ( x ) y n y ′ + p ( x ) y = q ( x ) y n where  p ( x ) p ( x )  and  q ( x ) q ( x )  are continuous functions on the interval we’re working on and  n n  is a real number. Differential equations in this form are called  Bernoulli Equations . First notice that if  n = 0 n = 0  or  n = 1 n = 1  then the equation is linear and we already know how to solve it in these cases. Therefore, in this section we’re going to be looking at solutions for values of  n n  other than these two. In order to solve these we’ll first divide the differential equation by  y n y n  to get, y − n y ′ + p ( x ) y 1 − n = q ( x ) y − n y ′ + p ( x ) y 1 − n = q ( x ) We are now going to use the substitution  v = y 1 − n v = y 1 − n  to convert this into a differential equation in terms of  v v . As we’ll see this will lead to a differential equation that we can solve. We are going to have to be c

Differential Equations - Systems: Solutions

Now that we’ve got some of the basics out of the way for systems of differential equations it’s time to start thinking about how to solve a system of differential equations. We will start with the homogeneous system written in matrix form, → x ′ = A → x (1) (1) x → ′ = A x → where,  A A  is an  n × n n × n  matrix and  → x x →  is a vector whose components are the unknown functions in the system. Now, if we start with  n = 1 n = 1 then the system reduces to a fairly simple linear (or separable) first order differential equation. x ′ = a x x ′ = a x and this has the following solution, x ( t ) = c e a t x ( t ) = c e a t So, let’s use this as a guide and for a general  n n  let’s see if → x ( t ) = → η e r t (2) (2) x → ( t ) = η → e r t will be a solution. Note that the only real difference here is that we let the constant in front of the exponential be a vector. All we need to do then is plug this into the differential equation and see what we get. First notice that

Digital Signal Processing - Miscellaneous Signals

There are other signals, which are a result of operation performed on them. Some common type of signals are discussed below. Conjugate Signals Signals, which satisfies the condition  x ( t ) = x ∗ ( − t ) are called conjugate signals. Let  x ( t ) = a ( t ) + j b ( t ) So,  x ( − t ) = a ( − t ) + j b ( − t ) And  x ∗ ( − t ) = a ( − t ) − j b ( − t ) By Condition,  x ( t ) = x ∗ ( − t ) If we compare both the derived equations 1 and 2, we can see that the real part is even, whereas the imaginary part is odd. This is the condition for a signal to be a conjugate type. Conjugate Anti-Symmetric Signals Signals, which satisfy the condition  x ( t ) = − x ∗ ( − t ) are called conjugate anti-symmetric signal Let  x ( t ) = a ( t ) + j b ( t ) So  x ( − t ) = a ( − t ) + j b ( − t ) And  x ∗ ( − t ) = a ( − t ) − j b ( − t ) − x ∗ ( − t ) = − a ( − t ) + j b ( − t ) By Condition  x ( t ) = − x ∗ ( − t ) Now, again compare, both the equations just as w

Calculus III - 3-Dimensional Space: Equations of Lines

In this section we need to take a look at the equation of a line in  R 3 R 3 . As we saw in the previous section the equation  y = m x + b y = m x + b  does not describe a line in  R 3 R 3 , instead it describes a plane. This doesn’t mean however that we can’t write down an equation for a line in 3-D space. We’re just going to need a new way of writing down the equation of a curve. So, before we get into the equations of lines we first need to briefly look at vector functions. We’re going to take a more in depth look at vector functions later. At this point all that we need to worry about is notational issues and how they can be used to give the equation of a curve. The best way to get an idea of what a vector function is and what its graph looks like is to look at an example. So, consider the following vector function. → r ( t ) = ⟨ t , 1 ⟩ r → ( t ) = ⟨ t , 1 ⟩ A vector function is a function that takes one or more variables, one in this case, and returns a vector. Note as we

Differential Equations - First Order: Modeling - i

We now move into one of the main applications of differential equations both in this class and in general. Modeling is the process of writing a differential equation to describe a physical situation. Almost all of the differential equations that you will use in your job (for the engineers out there in the audience) are there because somebody, at some time, modeled a situation to come up with the differential equation that you are using. This section is not intended to completely teach you how to go about modeling all physical situations. A whole course could be devoted to the subject of modeling and still not cover everything! This section is designed to introduce you to the process of modeling and show you what is involved in modeling. We will look at three different situations in this section : Mixing Problems, Population Problems, and Falling Objects. In all of these situations we will be forced to make assumptions that do not accurately depict reality in most cases, but wi

Differential Equations - Basic Concepts: Definitions

Differential Equation The first definition that we should cover should be that of  differential equation . A differential equation is any equation which contains derivatives, either ordinary derivatives or partial derivatives. There is one differential equation that everybody probably knows, that is Newton’s Second Law of Motion. If an object of mass  m m  is moving with acceleration  a a  and being acted on with force  F F  then Newton’s Second Law tells us. F = m a (1) (1) F = m a To see that this is in fact a differential equation we need to rewrite it a little. First, remember that we can rewrite the acceleration,  a a , in one of two ways. a = d v d t OR a = d 2 u d t 2 (2) (2) a = d v d t OR a = d 2 u d t 2 Where  v v  is the velocity of the object and  u u  is the position function of the object at any time  t t . We should also remember at this point that the force,  F F  may also be a function of time, velocity, and/or position. So, with all these things in

Differential Equations - Partial: Summary of Separation of Variables

Throughout this chapter we’ve been talking about and solving partial differential equations using the method of separation of variables. However, the one thing that we’ve not really done is completely work an example from start to finish showing each and every step. Each partial differential equation that we solved made use somewhere of the fact that we’d done at least part of the problem in another section and so it makes some sense to have a quick summary of the method here. Also note that each of the partial differential equations only involved two variables. The method can often be extended out to more than two variables, but the work in those problems can be quite involved and so we didn’t cover any of that here. So with all of that out of the way here is a quick summary of the method of separation of variables for partial differential equations in two variables. Verify that the partial differential equation is linear and homogeneous. Verify that the boundary condi

Differential Equations - Systems: Repeated Eigenvalues - i

This is the final case that we need to take a look at. In this section we are going to look at solutions to the system, → x ′ = A → x x → ′ = A x → where the eigenvalues are repeated eigenvalues. Since we are going to be working with systems in which  A A  is a  2 × 2 2 × 2  matrix we will make that assumption from the start. So, the system will have a double eigenvalue,  λ λ . This presents us with a problem. We want two linearly independent solutions so that we can form a general solution. However, with a double eigenvalue we will have only one, → x 1 = → η e λ t x → 1 = η → e λ t So, we need to come up with a second solution. Recall that when we looked at the double root case with the second order differential equations we ran into a similar problem. In that section we simply added a  t t  to the solution and were able to get a second solution. Let’s see if the same thing will work in this case as well. We’ll see if → x = t e λ t → η x → = t e λ t η → will also be a

Differential Equations - First Order: Modeling - ii

Example 4  A 50 kg object is shot from a cannon straight up with an initial velocity of 10m/s off a bridge that is 100 meters above the ground. If air resistance is given by 5 v v  determine the velocity of the mass when it hits the ground. First, notice that when we say straight up, we really mean straight up, but in such a way that it will miss the bridge on the way back down. Here is a sketch of the situation. Notice the conventions that we set up for this problem. Since the vast majority of the motion will be in the downward direction we decided to assume that everything acting in the downward direction should be positive. Note that we also defined the “zero position” as the bridge, which makes the ground have a “position” of 100. Okay, if you think about it we actually have two situations here. The initial phase in which the mass is rising in the air and the second phase when the mass is on its way down. We will need to examine both situations and set up an IVP for