- #1
tricha122
- 18
- 0
Hi,
I am trying to gain insight into using stationary vs. rotating coordinate frames for natural frequency calculations. I have seen many FE codes suggest that critical frequencies can be calculated differently in rotating and inertial frames, so i wanted to do a 1D calc to see for myself how things could be different.
I am trying to model a 1D spring-mass-damper with an rotating unbalance input force in both frames, and solve them numerically for a sweep of input frequencies to view the output displacement vs. frequency. And i expect to see a sharp increase where w = sqrt(K/M)
The equations of motion i am solving are as follows:
Stationary Frame:
M*x'' + C*x' + K*x = m*e*w^2*cos(w*t)
M*y'' + C*y' + K*y = m*e*w^2*sin(w*t)
these equations can be solved independently
For the rotating frame, i apply the transformation:
x = x1*cos(w*t) - y1*sin(w*t)
y = x1*sin(w*t) + y1*cos(w*t)
where x1, y1 are in the rotating frame, x1 is along the radius, y1 is "circumferential"
when i plug in this transformation, and simplify i get the following equations:
Mx1'' = -Kx1 - C(x1'-y1*w) + 2*M*y1'*w+m*(x1+e)*w^2
My1'' = -Ky1 - C(y1'+x1*w) - 2*M*x1'*w+m*y1*w^2
where the equations cannot be solved independently due to the coriolis force. (CF force shows up too).
I solved both stationary & rotating frame equations using M = 200, K = 104000, C = 10, m = 1, w = variable over a long "time" so that the damping could get rid of the "particular solution". After a long "time" i output the max amplitude. This resulted in an amplitude vs. input frequency chart
The stationary frame showed a sharp increase at w = ~22 = (sqrt(104000/200)) as expected.
The rotating frame showed no sharp increase, just steady rising with the forcing function. The strange part however was that at low frequencies, the results were identical to the stationary frame.
Can anyone shed some light on this? i have attached an image of my output. I also uploaded text files of the MATLAB code i was running.
I am trying to gain insight into using stationary vs. rotating coordinate frames for natural frequency calculations. I have seen many FE codes suggest that critical frequencies can be calculated differently in rotating and inertial frames, so i wanted to do a 1D calc to see for myself how things could be different.
I am trying to model a 1D spring-mass-damper with an rotating unbalance input force in both frames, and solve them numerically for a sweep of input frequencies to view the output displacement vs. frequency. And i expect to see a sharp increase where w = sqrt(K/M)
The equations of motion i am solving are as follows:
Stationary Frame:
M*x'' + C*x' + K*x = m*e*w^2*cos(w*t)
M*y'' + C*y' + K*y = m*e*w^2*sin(w*t)
these equations can be solved independently
For the rotating frame, i apply the transformation:
x = x1*cos(w*t) - y1*sin(w*t)
y = x1*sin(w*t) + y1*cos(w*t)
where x1, y1 are in the rotating frame, x1 is along the radius, y1 is "circumferential"
when i plug in this transformation, and simplify i get the following equations:
Mx1'' = -Kx1 - C(x1'-y1*w) + 2*M*y1'*w+m*(x1+e)*w^2
My1'' = -Ky1 - C(y1'+x1*w) - 2*M*x1'*w+m*y1*w^2
where the equations cannot be solved independently due to the coriolis force. (CF force shows up too).
I solved both stationary & rotating frame equations using M = 200, K = 104000, C = 10, m = 1, w = variable over a long "time" so that the damping could get rid of the "particular solution". After a long "time" i output the max amplitude. This resulted in an amplitude vs. input frequency chart
The stationary frame showed a sharp increase at w = ~22 = (sqrt(104000/200)) as expected.
The rotating frame showed no sharp increase, just steady rising with the forcing function. The strange part however was that at low frequencies, the results were identical to the stationary frame.
Can anyone shed some light on this? i have attached an image of my output. I also uploaded text files of the MATLAB code i was running.
