2011-04-01, 03:08 AM
Suppose that there is a mass hanging from a spring. Then the differential equation regarding this motion is:
m*x'' = -k*x + m*g
m*x'' + k*x = m*g
The convention is that downwards is positive and upwards is negative.
If we guess that x1 = e^(r*t) is a solution to the homogeneous equation:
m*x'' + k*x = 0
Then:
m*r^2*e^(r*t)+k*e^(r*t) = 0
e^(r*t)*(m*r^2+k)=0
r^2 = -k/m
r = +/- sqrt(-k/m)
r = +/- i*sqrt(k/m)
x1 = C1*e^(i*t*sqrt(k/m))+C2*e^(-i*t*sqrt(k/m))
x1 = (C1+C2)*cos(sqrt(k/m)*t)+(C1-C2)*i*sin(sqrt(k/m)*t)
C1=(C3-C4*i)/2, C2=(C3+C4*i)/2
x1 = C3*cos(sqrt(k/m)*t)+C4*i*sin(sqrt(k/m)*t)
If we guess that the solution to the particular equation:
m*x''+k*x=m*g
is in the form x2=A, where A is a constant:
k*A=m*g
A=m*g/k
x2=m*g/k
Then the general solution is:
x=C3*cos(sqrt(k/m)*t)+C4*sin(sqrt(k/m)*t)+m*g/k
Let's suppose that we pulled the spring some distance x(0) and released it. The initial condition is therefore:
x(0)=x0
x'(0)=0
x(0)=x0=C3-m*g/k
C3=x0+m*g/k
x'(t)=-sqrt(k/m)*C3*cos(sqrt(k/m)*t)+sqrt(k/m)*C4*cos(sqrt(k/m)*t)
x'(0)=0=sqrt(k/m)*C4
C4=0
Therefore, the complete solution to the initial value problem is:
x=(x0+m*g/k)*cos(sqrt(k/m)*t)+m*g/k
Mathematically, this works. Or at least, it should work out. However, logically, this doesn't make any sense. The mass is hanging in the presence of gravity, so the time it takes for the mass to fall should be shorter than the time it takes to rise. In other words, the mass should fall faster than it rises because it's in a gravitational field. Yet the equation implies that it doesn't. The first question is: wtp?
Also, a completely mathematical problem that I encountered is the conversion of this step:
x1 = C1*e^(i*t*sqrt(k/m))+C2*e^(-i*t*sqrt(k/m))
to sines and cosines. The first way is shown above, but I'll copy and paste it down here:
x1 = C1*e^(i*t*sqrt(k/m))+C2*e^(-i*t*sqrt(k/m))
x1 = (C1+C2)*cos(sqrt(k/m)*t)+(C1-C2)*i*sin(sqrt(k/m)*t)
C1=(C3-C4*i)/2, C2=(C3+C4*i)/2
x1 = C3*cos(sqrt(k/m)*t)+C4*i*sin(sqrt(k/m)*t)
The second way is:
x1 = C1*e^(i*t*sqrt(k/m))+C2*e^(-i*t*sqrt(k/m))
C1 = C5*^(C6*i), C2 = C7*e^(C8*i)
x1 = C5*e^(C6*i)*e^(i*t*sqrt(k/m))+C7*e^(C8*i)*e^(-i*t*sqrt(k/m))
x1 = C5*e^(i*t*sqrt(k/m)+i*C6)+C7*e^(-i*t*sqrt(k/m)+i*C8)
x1 = (C5+C7)*cos(sqrt(k/m)*t+C6)+(C5-C7)*i*sin(sqrt(k/m)*t+C8)
C3=(C9-C10*i)/2, C5=(C9+C10*i)/2
x1 = C9*cos(sqrt(k/m)*t+C6)+C10*sin(sqrt(k/m)*t+C7)
Both of them are mathematically correct, so hypothetically, they should also be mathematically equivalent. Yet I don't see how that's possible. Using the first simplification, the range is -(C3+C4)<x1<C3+C4. Using the second simplification, the range is -(C9+C10)< x1< C9+C10. |C9| and |C10| must be less than |C3| and |C4|, because a part of them went into shifting the phase of sine and cosine. The second question is: wtp?
tl;dr: What am I doing wrong?
m*x'' = -k*x + m*g
m*x'' + k*x = m*g
The convention is that downwards is positive and upwards is negative.
If we guess that x1 = e^(r*t) is a solution to the homogeneous equation:
m*x'' + k*x = 0
Then:
m*r^2*e^(r*t)+k*e^(r*t) = 0
e^(r*t)*(m*r^2+k)=0
r^2 = -k/m
r = +/- sqrt(-k/m)
r = +/- i*sqrt(k/m)
x1 = C1*e^(i*t*sqrt(k/m))+C2*e^(-i*t*sqrt(k/m))
x1 = (C1+C2)*cos(sqrt(k/m)*t)+(C1-C2)*i*sin(sqrt(k/m)*t)
C1=(C3-C4*i)/2, C2=(C3+C4*i)/2
x1 = C3*cos(sqrt(k/m)*t)+C4*i*sin(sqrt(k/m)*t)
If we guess that the solution to the particular equation:
m*x''+k*x=m*g
is in the form x2=A, where A is a constant:
k*A=m*g
A=m*g/k
x2=m*g/k
Then the general solution is:
x=C3*cos(sqrt(k/m)*t)+C4*sin(sqrt(k/m)*t)+m*g/k
Let's suppose that we pulled the spring some distance x(0) and released it. The initial condition is therefore:
x(0)=x0
x'(0)=0
x(0)=x0=C3-m*g/k
C3=x0+m*g/k
x'(t)=-sqrt(k/m)*C3*cos(sqrt(k/m)*t)+sqrt(k/m)*C4*cos(sqrt(k/m)*t)
x'(0)=0=sqrt(k/m)*C4
C4=0
Therefore, the complete solution to the initial value problem is:
x=(x0+m*g/k)*cos(sqrt(k/m)*t)+m*g/k
Mathematically, this works. Or at least, it should work out. However, logically, this doesn't make any sense. The mass is hanging in the presence of gravity, so the time it takes for the mass to fall should be shorter than the time it takes to rise. In other words, the mass should fall faster than it rises because it's in a gravitational field. Yet the equation implies that it doesn't. The first question is: wtp?
Also, a completely mathematical problem that I encountered is the conversion of this step:
x1 = C1*e^(i*t*sqrt(k/m))+C2*e^(-i*t*sqrt(k/m))
to sines and cosines. The first way is shown above, but I'll copy and paste it down here:
x1 = C1*e^(i*t*sqrt(k/m))+C2*e^(-i*t*sqrt(k/m))
x1 = (C1+C2)*cos(sqrt(k/m)*t)+(C1-C2)*i*sin(sqrt(k/m)*t)
C1=(C3-C4*i)/2, C2=(C3+C4*i)/2
x1 = C3*cos(sqrt(k/m)*t)+C4*i*sin(sqrt(k/m)*t)
The second way is:
x1 = C1*e^(i*t*sqrt(k/m))+C2*e^(-i*t*sqrt(k/m))
C1 = C5*^(C6*i), C2 = C7*e^(C8*i)
x1 = C5*e^(C6*i)*e^(i*t*sqrt(k/m))+C7*e^(C8*i)*e^(-i*t*sqrt(k/m))
x1 = C5*e^(i*t*sqrt(k/m)+i*C6)+C7*e^(-i*t*sqrt(k/m)+i*C8)
x1 = (C5+C7)*cos(sqrt(k/m)*t+C6)+(C5-C7)*i*sin(sqrt(k/m)*t+C8)
C3=(C9-C10*i)/2, C5=(C9+C10*i)/2
x1 = C9*cos(sqrt(k/m)*t+C6)+C10*sin(sqrt(k/m)*t+C7)
Both of them are mathematically correct, so hypothetically, they should also be mathematically equivalent. Yet I don't see how that's possible. Using the first simplification, the range is -(C3+C4)<x1<C3+C4. Using the second simplification, the range is -(C9+C10)< x1< C9+C10. |C9| and |C10| must be less than |C3| and |C4|, because a part of them went into shifting the phase of sine and cosine. The second question is: wtp?
tl;dr: What am I doing wrong?


![[Image: church_mass_kneel.jpg]](http://www.stclementshrine.org/images/church_mass_kneel.jpg)
![[Image: water%20spring.jpg]](http://www.filterclean.co.uk/images/infospringimages/water%20spring.jpg)