Jump to content
Main menu
Main menu
move to sidebar
hide
Navigation
Main page
Recent changes
Random page
Help about MediaWiki
Luter 345 Experiments
Search
Search
Appearance
Create account
Log in
Personal tools
Create account
Log in
Pages for logged out editors
learn more
Contributions
Talk
Editing
Theory
Page
Discussion
English
Read
Edit
Edit source
View history
Tools
Tools
move to sidebar
hide
Actions
Read
Edit
Edit source
View history
General
What links here
Related changes
Special pages
Page information
Get shortened URL
Appearance
move to sidebar
hide
Warning:
You are not logged in. Your IP address will be publicly visible if you make any edits. If you
log in
or
create an account
, your edits will be attributed to your username, along with other benefits.
Anti-spam check. Do
not
fill this in!
[[File:CNUPendulum.jpg|thumb|40px]] == The Purpose == The purpose of the Kater Pendulum experiment is to measure the local value of gravity, <math>g</math>. == The Device == Katerโs pendulum is a physical pendulum composed of a metal rod with two knife edge pivots towards each end of the rod. The exact location of the knife edges, relative to the ends of the rod, is not of extreme importance, as long as they are relatively symmetrical. However, the distance between the knife edges must be measured precisely* enough to achieve the level of accuracy desired in the experiment. There is a weight attached to each end of the pendulum. These two masses are of different weights. This is done intentionally and is the reason Kater's Pendulum works. The pendulum can be suspended by the knife edges and set to swinging by resting either knife edge on a flat, level surface in the small mass down (SMD) and small mass up (SMU) orientations. The "large" weight is attached at a fixed position and the small weight is adjustable. This adjustment allows the pendulum's period to be adjusted. As an improvement to the pendulum design, a thin "needle" has been attached to each end of the pendulum. As the period must be measured very precisely (see [[Uncertainties]]), these "needles" allow for a smaller angle of oscillation while still allowing the pendulum to pass through a photogate. == The Theory == [[File:KaterPendulum.png|thumb|300px]] When displaced, the restoring force acting on the pendulum is: : <math>\tau = M_{T} g h_1 \sin \theta</math> With <math>\theta</math> being the angle of displacement and the pendulum orientated small mass down (SMD), <math>M_{T}</math> being the total mass of the system, <math>g</math> being local gravity, and <math>h_1</math> being the distance from the pivot to the center of mass of the system. Applying Newton's Second Law rotational analog results in: : <math> M_{T} g h_1 \sin \theta = - \ddot \theta \mathbb{I}_{smd}</math> Where <math>\mathbb{I}_{smd}</math> is the moment of inertia. Of course, when <math>\theta</math> is small, the resulting equation can be simplified to: : <math>\ddot \theta + \frac{M_{T} g h_1 }{\mathbb{I}_{smd}}\theta = 0</math> Solving the differential equation results in: : <math> \theta = A \cos{\left(\omega_0 t\right)}+ B \cos{\left(\omega_0 t\right)}</math> Where: : <math>\omega_0 = \sqrt{\frac{M_{T} g h_1 }{\mathbb{I}_{smd}}}</math> resulting in: : <math>T_{smd} = 2\pi\sqrt{\frac{\mathbb{I}_{smd}}{M_{T} g h_1 }}</math> Where, <math>T_{smd}</math> is the period in the SMD orientation. Similarly, we can solve for the pendulum in the small mass up (SMU) orientation which results in: : <math>T_{smu} = 2\pi\sqrt{\frac{\mathbb{I}_{smu}}{M_{T} g h_2 }}</math> An application of the parallel axis theorm results on the equations: : <math>\mathbb{I}_{smd} = \mathbb{I}_{cm} + M_{T} h_1^2</math> and : <math>\mathbb{I}_{smu} = \mathbb{I}_{cm} + M_{T} h_2^2</math> Where <math>\mathbb{I}_{cm}</math> is the moment of inertia about the center of mass. By substituting these equations into our previous equations for <math>T</math> we get: : <math>T_{smd} = 2\pi\sqrt{\frac{\mathbb{I}_{cm} + M_{T} h_1^2}{M_{T} g h_1 }}</math> and : <math>T_{smu} = 2\pi\sqrt{\frac{\mathbb{I}_{cm} + M_{T} h_2^2}{M_{T} g h_2 }}</math> Now, if we can find a small mass position such that our period is the same in both the SMU and SMD orientation, we can set these equations equal to one another and solving for <math>\mathbb{I}_{cm}</math> we find: : <math>\mathbb{I}_{cm} = M_{T}h_1h_2</math> Substituting this value into either equation for period results in: : <math>T = 2\pi\sqrt{\frac{h_1 + h_2}{g}}</math> And, since the distance, <math>l</math>, between the pivots is the sum of <math>h_1</math> and <math>h_2</math> and we know this distance very well, we can write: : <math>T = 2\pi\sqrt{\frac{l}{g}}</math> or : <math>g = \frac{4\pi^2l}{T^2}</math> With this relationship, all we need to determine <math>g</math> is, the period of Kater's pendulum when the small mass is adjusted such that the period when swung in the SMU orientation is identical to the period in the SMD orientation. Some points to consider are if the weights are of different volume -which in our current pendulum, they are- their cross-sectional area will affect results due to air resistance. Also, if the weights are positioned very asymmetrically, the effect of air resistance will affect the results. The small angle approximation has been used to determine the values and the data should be corrected for this. == A Starting Estimation from First Principles == [[File:KaterPendulum.png|thumb|300px]] If we know the geometry and make up of the pendulum with a reasonable degree of accuracy, we are able to calculate where the small mass should be position so that it will produce the same period in the SMU and SMD orientation. === Variable Definitions === : <math>M \rightarrow</math> weight of large mass : <math>m \rightarrow</math> weight of small mass : <math>M_{bar} \rightarrow</math> weight of bar : <math>M_{T} = M + m + M_{bar} \rightarrow</math> total weight : <math>L_{bar} \rightarrow</math> length of bar : <math>g \rightarrow</math> gravity : <math>R_m \rightarrow</math> radius of small mass : <math>R_M \rightarrow</math> radius of large mass : <math>\mathbb{R} \rightarrow</math> radius of gyration : <math>\mathbb{I} \rightarrow</math> moment of inertia === Equation of Center of Mass === : <math>x = \frac{M l_1 - m l_2}{M + m + M_{bar}}</math> === Equations from the Pendulum Measurements === : <math>D = h_1 + h_2</math> : <math>h_2 - x = \frac{D}{2}</math> : <math>h_1 + x = \frac{D}{2}</math> Combining these equations results in: : <math>x = \frac{h_2 - h_1}{2}</math> === Equations of Period and Moment of Inertia === : <math>T_i^2 = \frac{4\pi^2}{g} \left( \frac{\mathbb{I}_i}{M_{T} h_i} \right)</math> : <math>\mathbb{I}_i = M_T \left( \R^2 + h_i^2 \right)</math> : <math>T_i^2 = \frac{4\pi^2}{g} \left( \frac{h_i^2+\R^2}{h_i} \right)</math> if <math>T_1 = T_2</math> then: : <math>\frac{h_1^2+\mathbb{R}^2}{h_1} = \frac{h_2^2+\mathbb{R}^2}{h_2}</math> : <math>h_1 h_2 = \mathbb{R}^2</math> : <math>\mathbb{I}_{cm} = \mathbb{R}^2 \left( M_T \right)</math> : <math>\mathbb{I}_{cm}^{sys} = \mathbb{I}_{M}^{cm} + \mathbb{I}_{m}^{cm} + \mathbb{I}_{bar}^{cm}</math> : <math>\mathbb{I}_{M}^{cm} = \frac{1}{2} M R_M^2 + M\left( l_1 - x \right)^2</math> : <math>\mathbb{I}_{m}^{cm} = \frac{1}{2} m R_m^2 + m\left( l_2 + x \right)^2</math> : <math>\mathbb{I}_{M}^{bar} = \mathbb{I}_{bar} + M_{bar} x^2</math> : <math>\mathbb{I}_{bar} = \frac{1}{12} M_{bar} L^2</math> All of which result in: : <math>h_1 h_2 = \frac{1}{M+m+M_{bar}} \left[ \frac{1}{2} M R_M^2 + M\left( l_1 - x \right)^2 + \frac{1}{2} m R_m^2 + m\left( l_2 + x \right)^2 + \frac{1}{12} M_{bar} L^2 + M_{bar} x^2 \right]</math> This is a somewhat difficult equation to solve. However, with the power and speed of modern computing, we don't have to. The following is a piece of code that finds a solution iteratively. === Iterative Solution === <syntaxhighlight lang="cpp"> #include <iostream> #include <TMath.h> void solve() { // Author: Edward J. Brash TCanvas *c1 = new TCanvas("c1","Kater's Pendulum",800,800); c1->SetFillColor(42); c1->SetGrid(); double x, M, m, Mbar, l1, h1, h2, D, L, rhs, lhs, Icm; double l2[1000],t1[1000],t2[1000]; double g=9.810; double Pi=3.14159265; M=1.35928; m=0.72705; Mbar=2.90; D=0.9990; L=1.523; l1=D/2.0+0.0135+0.047523; double diff=1.0e99; double olddiff=1.0e99; int index=0; for(int i=0; i<1000; i++){ l2[i]=0.500+(i/1000.0)*0.250; x=(M*l1-m*l2[i])/(M+m+Mbar); h2=D/2.0+x; h1=D/2.0-x; rhs=h1*h2; Icm=1.0/12.0*Mbar*L*L+Mbar*x*x+M*(l1-x)*(l1-x)+m*(l2[i]+x)*(l2[i]+x); lhs=Icm/(M+m+Mbar); diff = rhs-lhs; if (fabs(diff)<olddiff){olddiff=diff;index=i;} t1[i]=2*Pi*sqrt((h1*h1+lhs)/(g*h1)); t2[i]=2*Pi*sqrt((h2*h2+lhs)/(g*h2)); cout << l2[i] << " " << t1[i] << " " << t2[i] << endl; //cout << "i = "<<i<<" l2 = "<<l2[i]<<" rhs = "<<rhs<<" lhs = "<<lhs<<endl; } cout << "L2 position = " << l2[index] << " L1 position = " << l1 << endl; cout << "Predicted period T = " << t1[index]<< endl; TGraph *gra = new TGraph(1000,l2,t1); TGraph *grb = new TGraph(1000,l2,t2); gra->Draw("AP"); grb->Draw("P"); c1->Update(); c1->GetFrame()->SetFillColor(21); c1->GetFrame()->SetBorderSize(12); c1->Modified(); } </syntaxhighlight>
Summary:
Please note that all contributions to Luter 345 Experiments may be edited, altered, or removed by other contributors. If you do not want your writing to be edited mercilessly, then do not submit it here.
You are also promising us that you wrote this yourself, or copied it from a public domain or similar free resource (see
Luter 345 Experiments:Copyrights
for details).
Do not submit copyrighted work without permission!
Cancel
Editing help
(opens in new window)