What Frequency Range Keeps the Spring Safe in Driven Oscillations?

Click For Summary
SUMMARY

The discussion focuses on determining the safe frequency range for a driven oscillation of a 2 kg mass suspended by a spring with a strength of 2000 N/m. The driving force is given as 36 cos(pt) N, and the spring will yield if its extension exceeds 4 cm. The correct formulation of the differential equation is established as mx'' = 36 cos(pt) - kx - mg, where k is the spring constant and g is the acceleration due to gravity. The participants emphasize the importance of correctly accounting for forces acting on the mass, including gravitational and spring forces, to derive the appropriate frequency range.

PREREQUISITES
  • Understanding of Newton's second law of motion
  • Familiarity with differential equations, specifically second-order linear ODEs
  • Knowledge of harmonic motion and spring dynamics
  • Basic principles of driven oscillations and resonance
NEXT STEPS
  • Study the derivation of the equation of motion for driven harmonic oscillators
  • Learn about the effects of damping in oscillatory systems
  • Explore the concept of resonance and its implications in mechanical systems
  • Investigate the role of initial conditions in solving second-order differential equations
USEFUL FOR

Students and professionals in physics, mechanical engineering, and applied mathematics who are involved in the analysis of oscillatory systems and their stability under external forces.

  • #31
Hootenanny said:
Your solution looks okay, except that you're missing a factor i from your complementary function. In the first line of your solution you say assume that

x=e^ipt

And then you determine the values of p,

p_{1,2}=\pm10\sqrt{10}

Hence, your complementary function should be

x= c1eip1*t + c2eip2*t
Do you follow? It is the fact that these exponentials are complex, that we allow us to write the complimentary function in a form that is little easier to deal with.
yes but what about the function for the driven response ? I still have p^2 in my equation; I will right down what I mean:

x''-1000x= 36eipt -10
let x= ceipt. then x'= cipeipt and x''= -cp^2eipt;
Therefore ,-cp^2eipt-1000( ceipt)=36eipt -10 ; canceling out eipt , I am left with:

-cp^2-1000=36-10e-ipt; I still have that pesky e-iptterm in my equation , and therefore cannot have an amplitude composed of only real and imaginary parts and therefore cannot calculate xD.


To answer your question regarding extension, x is the extension of the spring beyond it's equilibrium position.

not sure what you mean. Should I subtract 4 cm from the equilibrium point?
 
Physics news on Phys.org
  • #32
Why are you choosing e^ipt as a particular solution?

The nonhomogeneous terms in your ODE (mx''+kx=36cos(pt)-mg) are the 36cos(pt) and -mg terms. Are you familiar with the method of undetermined coefficients? If so, what are the forms of the suggested trial particular solutions for each of these two terms?
 
  • #33
gabbagabbahey said:
Why are you choosing e^ipt as a particular solution?

The nonhomogeneous terms in your ODE (mx''+kx=36cos(pt)-mg) are the 36cos(pt) and -mg terms. Are you familiar with the method of undetermined coefficients? If so, what are the forms of the suggested trial particular solutions for each of these two terms?

There is a problem like this in my textbook where F(t) has a cosine term and Euler's formula is applied.

the solutions to this problem say p is safe if p < 20 rd/s and p>40 rd/s; I don't know what that means since I found that p1= +10\sqrt{10} and p2= -10\sqrt{10}
 
  • #34
Benzoate said:
There is a problem like this in my textbook where F(t) has a cosine term and Euler's formula is applied.

the solutions to this problem say p is safe if p < 20 rd/s and p>40 rd/s; I don't know what that means since I found that p1= +10\sqrt{10} and p2= -10\sqrt{10}

But p doesn't need to be one of those two values (p1,p2). p1 and p2 are the values needed to obey the homogeneous part of the ODE, not the inhomogeneous part.

The general form of the particular solution should be Ae^ipt +Be^-ipt +C not just e^ipt. The first two terms are necessary to account for the 36cos(pt) term on the RHS, while the constant term is needed to account for the -mg constant term.

What do you get when you plug this particular solution into the ODE?
 
  • #35
gabbagabbahey said:
But p doesn't need to be one of those two values (p1,p2). p1 and p2 are the values needed to obey the homogeneous part of the ODE, not the inhomogeneous part.

The general form of the particular solution should be Ae^ipt +Be^-ipt +C not just e^ipt. The first two terms are necessary to account for the 36cos(pt) term on the RHS, while the constant term is needed to account for the -mg constant term.

What do you get when you plug this particular solution into the ODE?

So I wouldn't write out use Euler term because of the extra constant term?

x= Ae^ipt +Be^-ipt +C

x'= ipAe^ipt -ipBe^-ipt + 0

x''=-p^2Ae^ipt -B^2e^-ipt+0


x''-1000x= 36eipt -10 , but F(t)=36 cos(t), so trig terms will not completely go away .
 
  • #36
I'm going to jump out here, two people helping in a thread isn't constructive.

Let me know if you need me to jump back in.
 
  • #37
Benzoate said:
So I wouldn't write out use Euler term because of the extra constant term?

x= -p^2Ae^ipt -B^2e^-ipt

x'= ipAe^ipt -ipBe^-ipt + 0

x''=-p^2Ae^ipt -B^2e^-ipt+0


x''-1000x= 36eipt -10 , but F(t)=36 cos(t), so trig terms will not completely go away .

You're missing a 1/m and cos(pt)≠e^ipt:

x''-1000x= (36/m)cos(pt) -10=18cos(pt)-10

=>-p^2Ae^ipt -Bp^2e^-ipt-1000(Ae^ipt -B^2e^-ipt)=18cos(pt)-10

But what is cos(pt) in terms of complex exponentials?
 
  • #38
gabbagabbahey said:
You're missing a 1/m and cos(pt)≠e^ipt:

x''-1000x= (36/m)cos(pt) -10=18cos(pt)-10

=>-p^2Ae^ipt -Bp^2e^-ipt-1000(Ae^ipt -B^2e^-ipt)=18cos(pt)-10

But what is cos(pt) in terms of complex exponentials?

Euler Formula:

ei\theta= cos(\theta)+sin(\theta)
therefore

18*eipt-10= 18*(cos(pt) + sin (pt))-10
 
Last edited:
  • #39
But you don't have 18 e^ipt! You have 18 cos(pt).

What is cos(pt) in term of complex exponentials?

Hint: look under the section "Relationship to Trignometry" here
 
Last edited:
  • #40
gabbagabbahey said:
But you don't have 18 e^ipt! You have 18 cos(pt).

What is cos(pt) in term of complex exponentials?

Hint: look under the section "Relationship to Trignometry" here


sorry if you are becoming frustrated. I have a copy of Gregory Douglass's Classical mechanics books and there is an example like this in that book on p. 10 cos(t) and they say the comple counter part is 10e^it

18 cos(pt)= 18e^ipt - 10(e^0)
 
  • #41
cos(pt)≠e^ipt

Euler's formula gives cos(pt)=Real[e^ipt]=(e^ipt+e^-ipt)/2 (Real[z] is the real part of the complex number z)

so, 18cos(pt)-10=9e^ipt+9e^-ipt-10

=>-p^2Ae^ipt -Bp^2e^-ipt+1000(Ae^ipt +Be^-ipt+C)=9e^ipt+9e^-ipt-10

Compare the coefficients in front of each of the e^ipt,e^-ipt, and constant terms...what must A,B and C be (in terms of p)?
 
Last edited:
  • #42
gabbagabbahey said:
cos(pt)≠e^ipt

Euler's formula gives cos(pt)=Real[e^ipt]=(e^ipt+e^-ipt)/2 (Real[z] is the real part of the complex number z)

so, 18cos(pt)-10=9e^ipt+9e^-ipt-10

=>-p^2Ae^ipt -Bp^2e^-ipt-1000(Ae^ipt -B^2e^-ipt+C)=9e^ipt+9e^-ipt-10

Compare the coefficients in front of each of the e^ipt,e^-ipt, and constant terms...what must A,B and C be (in terms of p)?

A(p^2-1000)=9 ==> A=9/(p^2-1000)
B(1000-p^2)=9 ==>B=9/(1000-p^2)
C=-10/1000?
 
  • #43
Close, there were some typos in my last equation; you should get A=9/(1000-p^2), your B and C are correct though.

Now, since A and B are equal you have: x(t)=A(e^ipt+e^-ipt)+C=2Acos(pt)+C

Now, if the maximum extension of the spring is 4cm, what must the Value of A be?
 
  • #44
gabbagabbahey said:
Close, there were some typos in my last equation; you should get A=9/(1000-p^2), your B and C are correct though.

Now, since A and B are equal you have: x(t)=A(e^ipt+e^-ipt)+C=2Acos(pt)+C

Now, if the maximum extension of the spring is 4cm, what must the Value of A be?
4=9/(1000-p^2)*cos(pt)-1/100 am I trying to find p?
 
  • #45
Yes, since x(t)=Acos(pt)+C, it should be clear that the angular frequency of oscillation is p. So you want to find p.

You need to be careful of your units though; 4cm =0.04m so you should have:

0.04=9/(1000-p^2)*cos(pt)-1/100

since the rest of the quantities in the equation are in meters.
 
  • #46
gabbagabbahey said:
Yes, since x(t)=Acos(pt)+C, it should be clear that the angular frequency of oscillation is p. So you want to find p.

You need to be careful of your units though; 4cm =0.04m so you should have:

0.04=9/(1000-p^2)*cos(pt)-1/100

since the rest of the quantities in the equation are in meters.

won't there be two solutions to p? The p's I finding for the xD are entirely different from the p's in the complementary function
 
  • #47
Sorry to bumped this thread again, even though its been two days since its been active.
for my final solution I get c=36/(1000-p^2) and c=.04==> p1=10 and p2=-10; I'm not sure what my text means when it says 'spring is safe if p=<20 rad/s and p>=40. rad/s
 

Similar threads

Replies
4
Views
2K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
9
Views
2K
Replies
8
Views
3K
Replies
2
Views
2K
  • · Replies 17 ·
Replies
17
Views
3K
Replies
1
Views
1K
  • · Replies 20 ·
Replies
20
Views
4K