Engr/Math/Physics 25. Accelerating Pendulum. Bruce Mayer, PE Licensed Electrical & Mechanical Engineer [email protected]. Recall 3 rd order Transformation. A 3 rd order Transformation (2). A 3 rd order Transformation (3). Thus the 3-Eqn 1 st Order ODE System. - PowerPoint PPT Presentation
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Bruce Mayer, PE Engineering/Math/Physics 25: Computational Methods
Problem 9.34 Accelerating
Pendulum For an Arbitrary Lateral-Acceleration Function, a(t), the ANGULAR Position, θ, is described by the (nastily) NONlinear 2nd Order, Homogeneous ODE
Bruce Mayer, PE Engineering/Math/Physics 25: Computational Methods
All Done for Today
FoucaultPendulum
While our clocks are set by an average 24 hour day for the passage of the Sun from noon to noon, the Earth rotates on its axis in 23 hours 56 minutes and 4.1 seconds with respect to the rest of the universe. From our perspective here on Earth, it appears that the entire universe circles us in this time. It is possible to do some rather simple experiments that demonstrate that it is really the rotation of the Earth that makes this daily motion occur.
In 1851 Leon Foucault (1819-1868) was made famous when he devised an experiment with a pendulum that demonstrated the rotation of the Earth.. Inside the dome of the Pantheon of Paris he suspended an iron ball about 1 foot in diameter from a wire more than 200 feet long. The ball could easily swing back and forth more than 12 feet. Just under it he built a circular ring on which he placed a ridge of sand. A pin attached to the ball would scrape sand away each time the ball passed by. The ball was drawn to the side and held in place by a cord until it was absolutely still. The cord was burned to start the pendulum swinging in a perfect plane. Swing after swing the plane of the pendulum turned slowly because the floor of the Pantheon was moving under the pendulum.
Bruce Mayer, PE Engineering/Math/Physics 25: Computational Methods
Prob 9.34 Script File
% Bruce Mayer, PE * 05Nov11% ENGR25 * problem 9.34% file = Demo_Prob9_34.m% %This script file calls FUNCTION pendacc%clear % clears memoryglobal m b; % globalize accel calc constants% Acceleration, a(t) = m*t + b% ask user for max time; suggest starting at 25tmax = input('tmax = '); %%set the case consts, and IC's y(0) & dy(0)/dt%=> remove the leading "%" to toggle between casesm = 0, b = 5, y0 = [0.5 0]; % case-a%m = 0, b = 5, y0 = [3 0]; % case-b%m = 0.5, b = 0, y0 = [3 0]; % case-c% m = 0.4, b = -4, y0 = [1.7 2.3]; % case-d => EXTRA%%Call the ode45 routine with the above data inputs[t,x]=ode45('pendacc', [0, tmax], y0);%%PLot theta(t)subplot(1,1,1)plot(t,x(:,1)), xlabel('t (sec)'), ylabel('theta (rads)'),...
title('P9.34 - Accelerating Pendulum'), grid;disp('Plotting ONLY theta - Hit Any Key to continue')pause%Plot the FIRST column of the solution “matrix” %giving x1 or y.subplot(2,1,1)plot(t,x(:,1)), xlabel('t (sec)'), ylabel('theta (rads)'),...
title('P9.34 - Accelerating Pendulum'), grid;%Plot the SECOND column of the solution “matrix” %giving x2 or dy/dt.subplot(2,1,2)plot(t,x(:,2)), xlabel('t (sec)'), ylabel('dtheta/dt (r/s)'), grid;disp('Plotting Both theta and dtheta/dt; hit any key to continue')
Bruce Mayer, PE Engineering/Math/Physics 25: Computational Methods
Prob 9.34 Function Filefunction dxdt = pendacc(t_val,z);% Bruce Mayer, PE * 05Nov05% ENGR25 * Prob 8-30% %This is the function that makes up the system %of differential equations solved by ode45%% the Vector z contains yk & [dy/dt]k%%Globalize the Constants used to calc the Accelglobal m b% set the physical constantsL = 1; % in mg = 9.81; % in m/sq-Sec%%DEBUG § => remove semicolons to reveal t_val & zt_val; z;%% Calc the Cauchy (State) valuesdxdt(1)= z(2); % at t=0, dxdt(1) = dy(0)/dtdxdt(2)= ((m*t_val + b)*cos(z(1)) - g*sin(z(1)))/L;% at t = 0, dxdt(2) =((m*t_val + b)*cos(y(0)) - g*sin(y(0)))/L; %% make the dxdt into a COLUMN vectordxdt = [dxdt(1); dxdt(2)];
Bruce Mayer, PE Engineering/Math/Physics 25: Computational Methods
Θ with Torsional Damping The Angular Position, θ, of a linearly
accelerating pendulum with a Journal Bearing mount that produces torsional friction-damping can be described by this second-order, non-linear Ordinary Differential Equation (ODE) and Initial Conditions (IC’s) for θ(t):
L
m W = mg
taD
0cossin btngDL
rads 8.20 secrads 9.100
tdt
d
L = 1.6 meters D = 0.07 meters/sec g = 9.8 meters/sec2