As it can be observed from Fig. a) and b), applying Newton s Law for the tangential force component results into: mg sin mat
|
|
- Brent Wilkinson
- 5 years ago
- Views:
Transcription
1 PHY4HF Exercse : Numercal ntegraton methods The Pendulum Startng wth small angles of oscllaton, you wll get expermental data on a smple pendulum and wll wrte a Python program to solve the equaton of moton. You wll plot the graph to vsualze the soluton (poston (angle) vs. tme). You wll also solve the energy equaton and plot t. You wll have to dscuss the output and eventually optmze the code. Background knowledge for Exercses -: Python: lsts, arrays, numercal ntegraton, scpy, pylab, leastsq. Error analyss: ch squared, goodness of the ft. R.Knght: Physcs for Scentsts and Engneers, nd ed., 8, 4.6: The Pendulum. Introducton Physcs of pendulum at small angles s based on applyng Newton s second law to derve the equaton of moton. The angle from the vertcal s θ, the dstance from the pvot pont s L; g = 9.8m/s. As t can be observed from Fg. a) and b), applyng Newton s Law for the tangental force component results nto: mg sn mat d s The tangental acceleraton s gven by: a t where the arc length dsplacement s s related to the angle θ by: s = Lθ. Brngng together all needed quanttes, we can wrte the equaton of moton of a pendulum at small angles of oscllaton as: () d g L
2 The most common method used to fndng solutons to equatons of moton s by settng up a par of coupled ordnary dfferental equatons. Gven: m = mass, q = coordnate, p = momentum and F = force, we can wrte: dp dq p F and : () m Consderng p and q to be the ntal values, we shall try to fnd solutons p( and q(.. Numercal methods To smplfy the equaton of moton, we approxmate the dervatves to: dy y( t y( t (3) Therefore, we can wrte: p( t p( F( q( ) t (4) p( q( t q( t m (4) s a set of update formulae that allow us to determne the numercal soluton (poston and momentum) at Δt, Δt, etc., gven the startng tme t = t. The numercal soluton wll approach the actual soluton as Δt. We can re-wrte (4) n a way closer to our Python code: p p F t p q q t m,,,..., t t.3 Numercal methods and the smple pendulum The equaton of moton for the smple pendulum () can be wrtten n the coupled form: d d and (6) where: Ω = g/l; θ s angle from the vertcal; ω s angular velocty. Our ntal condtons wll be: θ = 5 (small angle approxmaton) and ω =. Usng (6), the numercal approxmaton can be wrtten as: t t,,3... In equaton (7), q and s angular velocty. Ths s called the Forward Euler Method, because the rght-hand sde of (7) s evaluated at the ntal pont of the teraton step (5) (7)
3 .4 Python programmng (prelmnary) The basc steps you have to take are the followng: - defne constants - wrte ntal condtons - use numercal approxmaton (7) to step forward n tme - loop untl done - plot the graph - nterpret the result Remember that comments start wth # and the Python code s case-senstve..5 The lab exercse (All requrements marked by () have to be submtted to your TA). A pendulum consstng of a steel wre and a bob s attached to a rotatonal moton sensor, connected to a Natonal Instruments nterface. The output of the NI nterface s analyzed by a LabVew applcaton. Level the horzontal arm of the stand by usng the level provded and the knobs at the base. Measure the length of the pendulum and wegh the bob. Open the nd Yr Lab Fles folder from the lab computer desktop. Double clck on the RMS.v shortcut to open the LabVew applcaton. Before begnnng data acquston, get famlar wth cursors postons and the use of the Graph Palette to resze the graph. Clck on Acqure to start the program. In order to stop the acquston, clck on Acqure agan. The STOP button exts the program. When takng the pendulum out of equlbrum, rotate slowly the moton sensor wheel (do not touch the bob) to avod wobblng. Take the pendulum out of equlbrum by ~5 o (ths s the upper lmt of the lnear approxmaton leadng to equaton ()). Use the dsplay to setup the ntal angle. Remember that the rotatonal moton sensor sets up ts startng pont where you ntally brng the pendulum, and the LabVew applcaton reads angles n degrees. Start the acquston and let the pendulum swng for ~ seconds. In order to do the Python exercse, you wll need to determne Ω. Take the average of 5-6 oscllatons and use the cursors to obtan the perod of oscllaton T..6 Python programmng (plan) Get through ths program carefully because you may use t as a template for future applcatons. - Import the needed modules: pylab wll be needed for mathematcal lbrares and also for submodules (matplotlb.pyplo used to plottng the graphs: - Defne the needed constants: tme step Δt and Ω. - We need to calculate all the values for the plot. You took the expermental data for a total tme of seconds, so t s not a bad dea to do the calculaton up to t =. s. Gven deltat =. and t =, we ll have values. The best way to manpulate all the calculated values s to place them n an array (see the Compwk tutoral). Arrays are ndexed startng at zero. Therefore an array wth N elements has an ndex runnng from to N-. Our array wll be ntally set to zero. 3
4 - Known ntal condtons for our problem are: ntal angular velocty s zero. We have to setup the angle and the momentum ntal values. Angles have to be expressed n radans - Next step wll be wrtng the loop that ncreases the tme by Δt at each teraton. At step, the correspondng tme wll be Δt. The last tme value wll be. s - The last step wll be to plot the calculated data. We shall plot the angle on y-axs and tme on x-axs. We shall label the axes and put the plot on the screen. Wrte and run the program. Save the program on your memory stck Note: Sometmes, IDLE and Pylab do not work well together. You may need to open fewer wndows n order to keep the program stable. What happens? Are you confused? Can you nterpret the graph? To get more nformaton, try two other plots: ) angular velocty vs. tme (plot tme on x-axs and ω on y-axs) ) angular velocty vs. angle the phase plot (plot θ on x-axs and ω on y-axs). To do ths, you have to change the last lnes of code. Save each verson under a dfferent name. It could be very useful to analyze the total energy of the pendulum. All frst-year textbooks clam that energy of a smple pendulum s conserved: mv K(, U ( ) where: K s knetc energy ( K ml ) and U s E tot potental energy ( U mgh E mgl ). The energy expresson you have to use s: d ml mgl Modfy your program to calculate the energy at each step. Note that energy s not zero at t =. You have to nclude the length of the pendulum (L), the mass (m) and the gravtatonal constant (g) to the constants secton of your program. Plot: ) energy vs. tme What does the energy plot suggest? Does t explan the strange appearance of the ) and ) plots? Q. For a smple pendulum, the phase plot should be an ellpse. Usng energy conservaton, explan why. Now try to gve an explanaton for your phase plot. Q Determne the leadng error n our numercal method: perform a Taylor expanson of y(t + Δ and fnd the terms we have gnored n Equaton (4). Answer all the questons and submt all the code fles and plots to your demonstrator.. Numercal Integrators In the frst part of ths exercse, we analyzed the equaton of moton for a smple pendulum and wrote the Python code needed to numercally ntegrate the equaton. If you strctly used the numercal approxmaton, eq. (7), you must have found that everythng faled catastrophcally: angle ampltudes ncreased and energy was not conserved. The reason was the ntegraton method. We used the most prmtve numercal ntegraton method, called Euler Forward (explc: (8) 4
5 [ ] [ ] [ ] t [ ] [ ] [ ] t (9) The method ncreases the angle through an nterval Δt usng dervatve nformaton from only the begnnng of the nterval. It can be proved (see Reference at the end) that the method s unstable whch means that oscllaton ampltude as well as total energy monotoncally ncreases n tme. The numercal soluton shows a spral orbt n the (θ, ω) phase space. The method s accuracy and stablty can be mproved by decreasng the tme step, whch makes t attractve because of ts smplcty of mplementaton. Smlar to (9), the Euler Backward Method (mplc s gven by: [ ] [ ] [ ] t [ ] [ ] [ ] t The method s mplct because both [ ], [ ] are used on the rght hand sde; t s stable and therefore allows large tme steps to be taken. However, t nvolves some numercal dsspaton (see Reference). (9 ) A very smple remedy s to combne the two Euler methods nto: [ ] [ ] [ ] t () [ ] [ ] [ ] t The resultant s an explct method, stable, wth no numercal dsspaton, called Euler- Cromer or Symplectc Euler Method (SEM). Insert the new code lnes from SEM nto the program you wrote last tme and plot: Energy vs. tme and the phase plot. What happens? Compare and dscuss the phase plots from: Forward Euler, and Symplectc Method. Use the same tme step n all. Keep n mnd that n ntegratng conservatve problems t s essental to use a symplectc method.. Pendulum at large angles Wrte the Python code to ntegrate the equaton of moton of the pendulum at large angle. What you have to change are the followng: - ntal angle (n radans) - equatons of moton become: d d sn and - numercal soluton has to be re-wrtten - energy expresson s now: E ml mgl( cos ) 5
6 Use the modfed code template from wth SEM and change t to account for the large oscllaton angle. Extend the tme to mnutes. Plot Angle vs. tme, Energy vs. tme and the phase plot..3 Addng a dampng term. There s a large dscrepancy between the angle vs. tme Python plot and the real expermental data (open the RMS.v applcaton and let the pendulum swng for a couple of mnutes). It s obvous that the equatons of moton you used so far modeled the pendulum wthout dampng (physcal dsspaton). In general, the dampng force exerted on a body movng n ar or water depends on velocty v of the body relatve to the medum accordng to: C Av v () F d Coeffcents used n (3) are: C = drag coeffcent (dmensonless) A = cross-sectonal area perpendcular to the flow (m ) = densty of the medum (kg/m 3 ) v = lnear velocty of the body relatve to the medum (m/s) The drecton of the dampng (drag) force s always opposte to the drecton of velocty. The drag coeffcent C s not constant: C depends on body velocty, but also on vscosty of the medum, the shape of the body, and the roughness of ts surface. The Reynolds number R e has been found to be a useful dmensonless quantty that characterzes the dependence of the drag coeffcent on velocty. The Reynolds number s the rato of the nertal force of the medum to the vscous force: v R e () R e = Reynolds number (dmensonless) = Characterstc length of the body along the drecton of flow (m) = Dynamc vscosty of the medum (N s/m ) = Densty of the medum (kg/m 3 ) v = Lnear velocty of the body relatve to the medum (m/s) If the flow s lamnar, the Reynolds number takes small values (R e <3) and the drag coeffcent s nversely proportonal to the velocty. Ths makes the drag force drectly proportonal to the velocty: F d v When the flow s turbulent, the Reynolds number s large (4 < R e ) the drag coeffcent s approxmately constant and the drag force s dependent on the square of the velocty. Calculate the Reynolds number for the large angle pendulum. Moton and speed of a typcal pendulum bob n ar at large oscllaton angles correspond to small Reynolds numbers. Therefore, the equaton of moton of our pendulum would be wrtten as: d d sn (3) 6
7 d The dampng term was wrtten to nclude the lnear velocty: v L L. Usng cursors, take 5-6 readngs of ampltude. Assumng an exponental decay of the oscllaton envelope, estmate the decay constant γ (γ/ s the nverse of tme for whch ampltude falls to /e of the ntal value). You do not have to use Python to do ths calculaton..4 The Python applcaton The coupled equatons we need to formulate for our next Python applcaton are: d d sn (4) Use the template from part., modfy t accordng to (4) and plot Angle vs. tme. Dscuss the plot..5 Compare wth expermental data (qualtatvely) Open the RMS.v applcaton. Remember that the ntal poston of the rotatonal moton sensor sets up the orgn of the vertcal axs. Take the pendulum out of equlbrum by 3 o and let t swng for as long as you need n order to see a sgnfcant decay n ampltude. Qualtatvely compare the expermental data wth the output of your program. What do you thnk makes them dfferent? Note. Comparson wth some expermental data was justfed only for qualtatve purposes.we have not attempted to ft the data. Ths wll be the topc of the next exercse. All requrements marked by () have to be submtted to your TA. Reference: W.H. Press, S.A. Teukolsky, W.T. Vetterlng & B.P. Flannery (99), Numercal Recpes n C: the art of scentfc computng, nd ed. Cambrdge Unversty Press. Wrtten by Ruxandra Serbanescu (v.6, 9-5) 7
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationChapter 8. Potential Energy and Conservation of Energy
Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal
More informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationGravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11)
Gravtatonal Acceleraton: A case of constant acceleraton (approx. hr.) (6/7/11) Introducton The gravtatonal force s one of the fundamental forces of nature. Under the nfluence of ths force all objects havng
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationMoments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.
Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these
More information11. Dynamics in Rotating Frames of Reference
Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons
More informationSpin-rotation coupling of the angularly accelerated rigid body
Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationPHYS 705: Classical Mechanics. Newtonian Mechanics
1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]
More informationPlease initial the statement below to show that you have read it
EN40: Dynamcs and Vbratons Mdterm Examnaton Thursday March 5 009 Dvson of Engneerng rown Unversty NME: Isaac Newton General Instructons No collaboraton of any knd s permtted on ths examnaton. You may brng
More information10/23/2003 PHY Lecture 14R 1
Announcements. Remember -- Tuesday, Oct. 8 th, 9:30 AM Second exam (coverng Chapters 9-4 of HRW) Brng the followng: a) equaton sheet b) Calculator c) Pencl d) Clear head e) Note: If you have kept up wth
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationPhysics 111: Mechanics Lecture 11
Physcs 111: Mechancs Lecture 11 Bn Chen NJIT Physcs Department Textbook Chapter 10: Dynamcs of Rotatonal Moton q 10.1 Torque q 10. Torque and Angular Acceleraton for a Rgd Body q 10.3 Rgd-Body Rotaton
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More informationWeek 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product
The Vector Product Week 11: Chapter 11 Angular Momentum There are nstances where the product of two vectors s another vector Earler we saw where the product of two vectors was a scalar Ths was called the
More informationPart C Dynamics and Statics of Rigid Body. Chapter 5 Rotation of a Rigid Body About a Fixed Axis
Part C Dynamcs and Statcs of Rgd Body Chapter 5 Rotaton of a Rgd Body About a Fxed Axs 5.. Rotatonal Varables 5.. Rotaton wth Constant Angular Acceleraton 5.3. Knetc Energy of Rotaton, Rotatonal Inerta
More informationHow Differential Equations Arise. Newton s Second Law of Motion
page 1 CHAPTER 1 Frst-Order Dfferental Equatons Among all of the mathematcal dscplnes the theory of dfferental equatons s the most mportant. It furnshes the explanaton of all those elementary manfestatons
More informationEN40: Dynamics and Vibrations. Homework 4: Work, Energy and Linear Momentum Due Friday March 1 st
EN40: Dynamcs and bratons Homework 4: Work, Energy and Lnear Momentum Due Frday March 1 st School of Engneerng Brown Unversty 1. The fgure (from ths publcaton) shows the energy per unt area requred to
More informationPHYSICS 231 Review problems for midterm 2
PHYSICS 31 Revew problems for mdterm Topc 5: Energy and Work and Power Topc 6: Momentum and Collsons Topc 7: Oscllatons (sprng and pendulum) Topc 8: Rotatonal Moton The nd exam wll be Wednesday October
More informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
More informationStudy Guide For Exam Two
Study Gude For Exam Two Physcs 2210 Albretsen Updated: 08/02/2018 All Other Prevous Study Gudes Modules 01-06 Module 07 Work Work done by a constant force F over a dstance s : Work done by varyng force
More informationPhysics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.
Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
More informationLab 2e Thermal System Response and Effective Heat Transfer Coefficient
58:080 Expermental Engneerng 1 OBJECTIVE Lab 2e Thermal System Response and Effectve Heat Transfer Coeffcent Warnng: though the experment has educatonal objectves (to learn about bolng heat transfer, etc.),
More informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More informationEN40: Dynamics and Vibrations. Homework 7: Rigid Body Kinematics
N40: ynamcs and Vbratons Homewor 7: Rgd Body Knematcs School of ngneerng Brown Unversty 1. In the fgure below, bar AB rotates counterclocwse at 4 rad/s. What are the angular veloctes of bars BC and C?.
More informationThe classical spin-rotation coupling
LOUAI H. ELZEIN 2018 All Rghts Reserved The classcal spn-rotaton couplng Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 louaelzen@gmal.com Abstract Ths paper s prepared to show that a rgd
More informationHandout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system. and write EOM (1) as two first-order Eqs.
Handout # 6 (MEEN 67) Numercal Integraton to Fnd Tme Response of SDOF mechancal system State Space Method The EOM for a lnear system s M X + DX + K X = F() t () t = X = X X = X = V wth ntal condtons, at
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationCelestial Mechanics. Basic Orbits. Why circles? Tycho Brahe. PHY celestial-mechanics - J. Hedberg
PHY 454 - celestal-mechancs - J. Hedberg - 207 Celestal Mechancs. Basc Orbts. Why crcles? 2. Tycho Brahe 3. Kepler 4. 3 laws of orbtng bodes 2. Newtonan Mechancs 3. Newton's Laws. Law of Gravtaton 2. The
More informationFirst Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.
Secton 1. Dynamcs (Newton s Laws of Moton) Two approaches: 1) Gven all the forces actng on a body, predct the subsequent (changes n) moton. 2) Gven the (changes n) moton of a body, nfer what forces act
More informationModeling of Dynamic Systems
Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More information6.3.4 Modified Euler s method of integration
6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from
More informationENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15
NGN 40 ynamcs and Vbratons Homework # 7 ue: Frday, Aprl 15 1. Consder a concal hostng drum used n the mnng ndustry to host a mass up/down. A cable of dameter d has the mass connected at one end and s wound/unwound
More informationThree views of mechanics
Three vews of mechancs John Hubbard, n L. Gross s course February 1, 211 1 Introducton A mechancal system s manfold wth a Remannan metrc K : T M R called knetc energy and a functon V : M R called potental
More informationPhysics 141. Lecture 14. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 14, Page 1
Physcs 141. Lecture 14. Frank L. H. Wolfs Department of Physcs and Astronomy, Unversty of Rochester, Lecture 14, Page 1 Physcs 141. Lecture 14. Course Informaton: Lab report # 3. Exam # 2. Mult-Partcle
More informationSo far: simple (planar) geometries
Physcs 06 ecture 5 Torque and Angular Momentum as Vectors SJ 7thEd.: Chap. to 3 Rotatonal quanttes as vectors Cross product Torque epressed as a vector Angular momentum defned Angular momentum as a vector
More informationMEASUREMENT OF MOMENT OF INERTIA
1. measurement MESUREMENT OF MOMENT OF INERTI The am of ths measurement s to determne the moment of nerta of the rotor of an electrc motor. 1. General relatons Rotatng moton and moment of nerta Let us
More informationChapter 11: Angular Momentum
Chapter 11: ngular Momentum Statc Equlbrum In Chap. 4 we studed the equlbrum of pontobjects (mass m) wth the applcaton of Newton s aws F 0 F x y, 0 Therefore, no lnear (translatonal) acceleraton, a0 For
More informationSpring 2002 Lecture #13
44-50 Sprng 00 ecture # Dr. Jaehoon Yu. Rotatonal Energy. Computaton of oments of nerta. Parallel-as Theorem 4. Torque & Angular Acceleraton 5. Work, Power, & Energy of Rotatonal otons Remember the md-term
More informationAPPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS
Unversty of Oulu Student Laboratory n Physcs Laboratory Exercses n Physcs 1 1 APPEDIX FITTIG A STRAIGHT LIE TO OBSERVATIOS In the physcal measurements we often make a seres of measurements of the dependent
More informationPHYS 1441 Section 002 Lecture #16
PHYS 1441 Secton 00 Lecture #16 Monday, Mar. 4, 008 Potental Energy Conservatve and Non-conservatve Forces Conservaton o Mechancal Energy Power Today s homework s homework #8, due 9pm, Monday, Mar. 31!!
More informationConservation of Angular Momentum = "Spin"
Page 1 of 6 Conservaton of Angular Momentum = "Spn" We can assgn a drecton to the angular velocty: drecton of = drecton of axs + rght hand rule (wth rght hand, curl fngers n drecton of rotaton, thumb ponts
More informationCOMPLEX NUMBERS AND QUADRATIC EQUATIONS
COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationONE-DIMENSIONAL COLLISIONS
Purpose Theory ONE-DIMENSIONAL COLLISIONS a. To very the law o conservaton o lnear momentum n one-dmensonal collsons. b. To study conservaton o energy and lnear momentum n both elastc and nelastc onedmensonal
More informationLecture 2: Numerical Methods for Differentiations and Integrations
Numercal Smulaton of Space Plasmas (I [AP-4036] Lecture 2 by Lng-Hsao Lyu March, 2018 Lecture 2: Numercal Methods for Dfferentatons and Integratons As we have dscussed n Lecture 1 that numercal smulaton
More informationPhysics 207: Lecture 20. Today s Agenda Homework for Monday
Physcs 207: Lecture 20 Today s Agenda Homework for Monday Recap: Systems of Partcles Center of mass Velocty and acceleraton of the center of mass Dynamcs of the center of mass Lnear Momentum Example problems
More informationRotational Dynamics. Physics 1425 Lecture 19. Michael Fowler, UVa
Rotatonal Dynamcs Physcs 1425 Lecture 19 Mchael Fowler, UVa Rotatonal Dynamcs Newton s Frst Law: a rotatng body wll contnue to rotate at constant angular velocty as long as there s no torque actng on t.
More informationPhysics 2A Chapters 6 - Work & Energy Fall 2017
Physcs A Chapters 6 - Work & Energy Fall 017 These notes are eght pages. A quck summary: The work-energy theorem s a combnaton o Chap and Chap 4 equatons. Work s dened as the product o the orce actng on
More informationHomework 2: Kinematics and Dynamics of Particles Due Friday Feb 7, 2014 Max Score 45 Points + 8 Extra Credit
EN40: Dynamcs and Vbratons School of Engneerng Brown Unversty Homework : Knematcs and Dynamcs of Partcles Due Frday Feb 7, 014 Max Score 45 Ponts + 8 Extra Credt 1. An expermental mcro-robot (see a descrpton
More informationWorkshop: Approximating energies and wave functions Quantum aspects of physical chemistry
Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department
More informationLagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013
Lagrange Multplers Monday, 5 September 013 Sometmes t s convenent to use redundant coordnates, and to effect the varaton of the acton consstent wth the constrants va the method of Lagrange undetermned
More informationChapter 11 Angular Momentum
Chapter 11 Angular Momentum Analyss Model: Nonsolated System (Angular Momentum) Angular Momentum of a Rotatng Rgd Object Analyss Model: Isolated System (Angular Momentum) Angular Momentum of a Partcle
More informationPhysics 106a, Caltech 11 October, Lecture 4: Constraints, Virtual Work, etc. Constraints
Physcs 106a, Caltech 11 October, 2018 Lecture 4: Constrants, Vrtual Work, etc. Many, f not all, dynamcal problems we want to solve are constraned: not all of the possble 3 coordnates for M partcles (or
More informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationIrregular vibrations in multi-mass discrete-continuous systems torsionally deformed
(2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected
More informationWeek 8: Chapter 9. Linear Momentum. Newton Law and Momentum. Linear Momentum, cont. Conservation of Linear Momentum. Conservation of Momentum, 2
Lnear omentum Week 8: Chapter 9 Lnear omentum and Collsons The lnear momentum of a partcle, or an object that can be modeled as a partcle, of mass m movng wth a velocty v s defned to be the product of
More informationAdvanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)
Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed
More informationA Tale of Friction Basic Rollercoaster Physics. Fahrenheit Rollercoaster, Hershey, PA max height = 121 ft max speed = 58 mph
A Tale o Frcton Basc Rollercoaster Physcs Fahrenhet Rollercoaster, Hershey, PA max heght = 11 t max speed = 58 mph PLAY PLAY PLAY PLAY Rotatonal Movement Knematcs Smlar to how lnear velocty s dened, angular
More informationKinematics in 2-Dimensions. Projectile Motion
Knematcs n -Dmensons Projectle Moton A medeval trebuchet b Kolderer, c1507 http://members.net.net.au/~rmne/ht/ht0.html#5 Readng Assgnment: Chapter 4, Sectons -6 Introducton: In medeval das, people had
More informationComparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy
Comparatve Studes of Law of Conservaton of Energy and Law Clusters of Conservaton of Generalzed Energy No.3 of Comparatve Physcs Seres Papers Fu Yuhua (CNOOC Research Insttute, E-mal:fuyh1945@sna.com)
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationPHYS 705: Classical Mechanics. Hamilton-Jacobi Equation
1 PHYS 705: Classcal Mechancs Hamlton-Jacob Equaton Hamlton-Jacob Equaton There s also a very elegant relaton between the Hamltonan Formulaton of Mechancs and Quantum Mechancs. To do that, we need to derve
More informationCHAPTER 10 ROTATIONAL MOTION
CHAPTER 0 ROTATONAL MOTON 0. ANGULAR VELOCTY Consder argd body rotates about a fxed axs through pont O n x-y plane as shown. Any partcle at pont P n ths rgd body rotates n a crcle of radus r about O. The
More informationLAB # 4 - Torque. d (1)
LAB # 4 - Torque. Introducton Through the use of Newton's three laws of moton, t s possble (n prncple, f not n fact) to predct the moton of any set of partcles. That s, n order to descrbe the moton of
More information= 1.23 m/s 2 [W] Required: t. Solution:!t = = 17 m/s [W]! m/s [W] (two extra digits carried) = 2.1 m/s [W]
Secton 1.3: Acceleraton Tutoral 1 Practce, page 24 1. Gven: 0 m/s; 15.0 m/s [S]; t 12.5 s Requred: Analyss: a av v t v f v t a v av f v t 15.0 m/s [S] 0 m/s 12.5 s 15.0 m/s [S] 12.5 s 1.20 m/s 2 [S] Statement:
More informationLecture 16. Chapter 11. Energy Dissipation Linear Momentum. Physics I. Department of Physics and Applied Physics
Lecture 16 Chapter 11 Physcs I Energy Dsspaton Lnear Momentum Course webste: http://aculty.uml.edu/andry_danylov/teachng/physcsi Department o Physcs and Appled Physcs IN IN THIS CHAPTER, you wll learn
More informationReview of Taylor Series. Read Section 1.2
Revew of Taylor Seres Read Secton 1.2 1 Power Seres A power seres about c s an nfnte seres of the form k = 0 k a ( x c) = a + a ( x c) + a ( x c) + a ( x c) k 2 3 0 1 2 3 + In many cases, c = 0, and the
More information10/24/2013. PHY 113 C General Physics I 11 AM 12:15 PM TR Olin 101. Plan for Lecture 17: Review of Chapters 9-13, 15-16
0/4/03 PHY 3 C General Physcs I AM :5 PM T Oln 0 Plan or Lecture 7: evew o Chapters 9-3, 5-6. Comment on exam and advce or preparaton. evew 3. Example problems 0/4/03 PHY 3 C Fall 03 -- Lecture 7 0/4/03
More informationTHE CURRENT BALANCE Physics 258/259
DSH 1988, 005 THE CURRENT BALANCE Physcs 58/59 The tme average force between two parallel conductors carryng an alternatng current s measured by balancng ths force aganst the gravtatonal force on a set
More information10/9/2003 PHY Lecture 11 1
Announcements 1. Physc Colloquum today --The Physcs and Analyss of Non-nvasve Optcal Imagng. Today s lecture Bref revew of momentum & collsons Example HW problems Introducton to rotatons Defnton of angular
More information= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.
Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationA particle in a state of uniform motion remain in that state of motion unless acted upon by external force.
The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,
More informationTHERMAL DISTRIBUTION IN THE HCL SPECTRUM OBJECTIVE
ame: THERMAL DISTRIBUTIO I THE HCL SPECTRUM OBJECTIVE To nvestgate a system s thermal dstrbuton n dscrete states; specfcally, determne HCl gas temperature from the relatve occupatons of ts rotatonal states.
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationChapter Newton s Method
Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve
More informationAssessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion
Assessment of Ste Amplfcaton Effect from Input Energy Spectra of Strong Ground Moton M.S. Gong & L.L Xe Key Laboratory of Earthquake Engneerng and Engneerng Vbraton,Insttute of Engneerng Mechancs, CEA,
More informationPhysics 207 Lecture 6
Physcs 207 Lecture 6 Agenda: Physcs 207, Lecture 6, Sept. 25 Chapter 4 Frames of reference Chapter 5 ewton s Law Mass Inerta s (contact and non-contact) Frcton (a external force that opposes moton) Free
More informationAngular Momentum and Fixed Axis Rotation. 8.01t Nov 10, 2004
Angular Momentum and Fxed Axs Rotaton 8.01t Nov 10, 2004 Dynamcs: Translatonal and Rotatonal Moton Translatonal Dynamcs Total Force Torque Angular Momentum about Dynamcs of Rotaton F ext Momentum of a
More informationIterative General Dynamic Model for Serial-Link Manipulators
EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators 1. Introducton Iteratve General Dynamc Model for Seral-Lnk Manpulators In ths set of notes, we are gong to develop a method for computng a general
More informationSome modelling aspects for the Matlab implementation of MMA
Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton
More informationTIME OF COMPLETION NAME SOLUTION DEPARTMENT OF NATURAL SCIENCES. PHYS 2211, Exam 2 Section 1 Version 1 October 18, 2013 Total Weight: 100 points
TIME OF COMPLETION NAME SOLUTION DEPARTMENT OF NATURAL SCIENCES PHYS, Exam Secton Verson October 8, 03 Total Weght: 00 ponts. Check your examnaton or completeness pror to startng. There are a total o nne
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationPhysics 4B. Question and 3 tie (clockwise), then 2 and 5 tie (zero), then 4 and 6 tie (counterclockwise) B i. ( T / s) = 1.74 V.
Physcs 4 Solutons to Chapter 3 HW Chapter 3: Questons:, 4, 1 Problems:, 15, 19, 7, 33, 41, 45, 54, 65 Queston 3-1 and 3 te (clockwse), then and 5 te (zero), then 4 and 6 te (counterclockwse) Queston 3-4
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationAP Physics 1 & 2 Summer Assignment
AP Physcs 1 & 2 Summer Assgnment AP Physcs 1 requres an exceptonal profcency n algebra, trgonometry, and geometry. It was desgned by a select group of college professors and hgh school scence teachers
More informationFTCS Solution to the Heat Equation
FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use
More informationC/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1
C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned
More information( ) = ( ) + ( 0) ) ( )
EETOMAGNETI OMPATIBIITY HANDBOOK 1 hapter 9: Transent Behavor n the Tme Doman 9.1 Desgn a crcut usng reasonable values for the components that s capable of provdng a tme delay of 100 ms to a dgtal sgnal.
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More information