Classical mechanics of particles and fields
|
|
- Brooke Watkins
- 5 years ago
- Views:
Transcription
1 Classical mechanics of particles and fields D.V. Skryabin Department of Physics, University of Bath PACS numbers: The concise and transparent exposition of many topics covered in this unit can be found in: L. Landau and E. Lifshitz, Mechanics. Course of theoretical physics vol. 1. (library has many copies) 1. Introduction to variational calculus Main idea of the Lagrangian mechanics is that the true trajectory from the poin in space-time to the point 2 gives the minimal value of the integral A = A = (kin.energy pot.energy) ( ) mẋ 2 U(x) 2 It means: if we calculate A for some x(t) we always will find a number, which is greater than A calculated for the true, i.e. physically realizable, x(t). A is called action. To formulate this idea more rigorously we need to learn how to find functions which give stationary points (e.g. maxima or minima) for this type of integrals. We need to learn Variational calculus Functional is a mathematical operation, which takes a function x(t) and produces a number: A = A[x(t)] FIG. 1: Possible trajectories. The true one gives minimum of the action integral.
2 2 Typical example of a functional is a definite integral: A[x(t)] = L(x(t), t) Here L is a GIVEN function of its arguments (x and t) and x is itself a function of t. A is a number which obviously depends on our choice of function x(t) for a given L. One shell remember that our functional operates on an infinite set of functions x(t). Useful analogy is a function of a single independent variable operating on the entire set of real numbers from to +. To find extrema of a functional we, first, take a small deviation δx from a chosen (out of infinite number of possibilities) function x(t). δx is called variation of x. Then we calculate A[x + δx] A[x] upto the first order in δx and the result is called variation (or first variation) of A Let us consider A[x(t)] = t 2 using the outlined method. L(x(t), t) and calculate its variation and variational derivative A[x(t) + δx(t)] = L(x(t) + δx(t), t) but L is a usual function so we can expand into Taylor series. and Thus A[x(t) + δx(t)] A[x(t)] = L(x(t) + δx(t), t) = L(x(t), t) + L x δx(t) +... [ L(x(t), t) + L ] x δx(t) +... L(x(t), t) L δa = x δx. Function between the integral sign and δx is called variational derivative of a functional A: δa/δx. The above definition of the variational derivative is analogous to the following fact known from the usual calculus: if F (t) = x(t), then df/ = x(t). δa δx = L x Function x 0 (t) for which our functional reaches its stationary value must be such that δa = 0. x 0 (t) is then called an extremum of the functional A
3 3 FIG. 2: Lagrange (left): and Leonard Euler (right): Because δx is arbitrary, it means that δa is zero only if δa δx = L x = 0. The trivial results δa = L is valid only if L = L(x(t), t), but not if L = L(x, ẋ, t) δx x Let s work out variation of A[x(t)] if L = L(x, ẋ, t) using the same methodology... δa = { δx(t) [ ] [ ]} L L + δẋ x ẋ Using integration by parts we manipulate the second part of the integral above δẋ [ ] L dδx = ẋ [ ] [ ] t2 L L = ẋ ẋ δx t2 δx d t 1 We consider variations with the fixed ends, i.e. δx( ) = δx(t 2 ) = 0. {[ ]} L ẋ This gives δa = [ L x d ] L δx ẋ Thus if x 0 (t) is such that it solves the equation δa δx = L x d L ẋ L x d L ẋ = 0
4 4 FIG. 3: Snell s law then our functional A[x(t)] reaches its stationary value A 0 for x = x 0 (t) and x 0 (t) is called extremum of A Equation above is called Euler-Lagrange s or simply Lagrange s equation. Thus we have demonstrated that the problem of finding stationary values for the functional can be reduced to differential equations. 1.1 Variational formulation of the geometrical optics: Fermat s principle A ray of light follows a path such that the total path time is stationary (minimal) Mathematically it means that along the true path δ = 0. l = dl v = n(l)dl c where n(l) is refractive index varying along the path. Hence δ n(l)dl = 0. l Thus the equivalent formulation of the Fermat principle is that a ray follows a trajectory such that the total optical path (length refractive index) is stationary. In cartesian coordinates we reformulate the problem, as finding such y(x), that gives a stationary value to the functional x2 F = n(x, y) 1 + (dy/dx) 2 dx x 1 Let s consider an example when n is a function of y only and does not explicitly contain x. Then δf = 1 x2 dxδy d [ ] n(y) = 0 2 dx 1 + ẏ 2 Thus x 1 n(y) 1 + ẏ 2 = k = const
5 5 FIG. 4: William Hamilton: k is the constant determined by initial conditions. Take a simple example of an interface (at y = 0) between two homogeneous media with different refractive indices n 1 and n 2. ẏ = tanφ hence Hence which is the Snell s law n cos φ = const = n sin θ n 1 sin θ 1 = n 2 sin θ 2 2. Hamilton s principle of mechanics or principle of stationary action For conservative systems, i.e. when forces can be derived from potentials the following is true Any mechanical system follows a path that yields the stationary value of the action functional: A = t 2 L. L is called Lagrangian and it is the difference between the kinetic and potential energies of the system. δa = 0 for the true path. It is an axiomatic law, which is proved by the comparison with experiment. In many, but NOT ALL, practical cases A has a minimum value, when calculated along the true trajectory. Therefore the above principle is often called principle of least action. However, it can be a maximum sometimes, or it can be a saddle point, i.e. second variations along different coordinates can have different signs. To understand whether a given stationary trajectory corresponds to a minimum or maximum one has to calculate second variations (we are not doing it in this course). Lagrangian of a particle in a potential We do not have time to verify our theory by experiments, therefore we will verify it by comparison with the Newton equations. mẍ = F or equivalently dp = du dx,
6 where U is the potential energy and p = mẋ is the momentum. Thus one can see that the unusually looking Lagrange equation d L ẋ = L x is surprisingly similar to the 2nd law of Newton and even equivalent to it providing that Using that L = m L ẋ p p = L ẋ, L x = du dx and integrating the first of the above equations we find that L = p2 2m + ψ(x), where ψ is an unknown function of x. Integrating the 2nd of the above equations we find that L = U(x) + φ(p), where φ is an unknown function of φ. Requiring compatibility of the 2 expressions for L we find that ψ(x) = U(x) and that φ(p) = p 2 /(2m). Thus Lagrangian of a particle moving in the potential U(x) is L = p2 U(x) = T U, 2m where T is the kinetic energy and V is the potential one. We have demonstrated that the Lagrangian is the kinetic energy minus potential energy. Principle of Stationary-Action also can be formulated as: The path of a particle obeys the Euler-Lagrange equation [ ] d L = L ẋ x One might like to remember it in the form d [ ] L = L v x In 3D the principle of stationary action gives system of 3 Lagrange equations [ ] d L = L [ ] ẋ x, d L = L [ ] ẏ y, d L = L ż z Here L = m 2 (ẋ2 + ẏ 2 + ż 2 ) U(x, y, z) For shortness the 3D Lagrange equations can be written as [ ] d L r = L r where the vector derivatives should be understood as gradients and then corresponding gradients components in the left-hand side should be made equal to the ones in the right-hand side. Examples Ball dropped under the angle. L = mẋ 2 /2 + mẏ 2 /2 + mgx and Lagrange equation is ẍ = g ÿ = 0 (Galileo equation). Mass on a horizontal spring. T = mẋ 2 /2, U = kx 2 /2 and the Lagrange equation is mẍ = kx. 6
7 7 3. Universality of Lagrange equations Generalized coordinates, momenta and forces Let us now consider a system of particles, which interact one with another, but not with other bodies or external forces. Lagrangian of this system is a sum of Lagrangians of the free particles, minus the interaction energy. L = 1 N m a v 2 a U( r 1, r 2,...) 2 a=1 a numbers the particles. The fact that the potential energy depends only on positions of particles, but not on their velocities, shows, that interactions are instantaneous. If the propagation of interactions is not instantaneous, then that velocity is different in different inertial reference frames, which contradicts to the relativity principle. Equations of motion are d L = L dx a, m a ẋ a x a = U x a If L is written in terms of cartesian coordinates x, y, z, than the Euler-Lagrange equations are equivalent to Newtons F = ma equations. But what about the cases where we use polar, spherical, or some other coordinates or dynamical variables to describe motion? Many particle systems is a particularly obvious example, when the coordinates different from cartesian ones can be very convenient If we try to describe a system of N particles, we need in general 3N coordinates to specify the position of its components. If external constraints are imposed on the system (as it often happens in practice), the number of degrees of freedom may be less. If there are m constraints applied, the number of degrees of freedom will be s = 3N m (give some example (pendulum, molecule)). The coordinates do not need to be the coordinates of a coordinate system, but can be any set of quantities that completely specifies the state of the system. These quantities are called generalized coordinates: q n. The coordinate space based on the generalized coordinates is called configuration space. The instantaneous state of the system is fully specified by a point in the s = 3N m dimensional configuration space. The time evolution of the system can be described by a path in the configuration space. One of the great advantages of using Lagrangians is that Euler-Lagrange equation look the same in all coordinate systems. You can use any convenient degrees of freedom and both inertial and non-inertial frames. In any case d L q i = L q i. If we assume x a = f a (q 1, q 2,... q s ), ẋ a = k f a q k q k and substitute these expressions into the Lagrangian L = 1 2 Na=1 m a v a 2 U( r 1, r 2,...), we find L = i,k 1 2 a ik( q) q i q k U( q) The kinetic energy in the generalized coordinates is still a quadratic function of the velocities, but may depend on the coordinates as well
8 8 Whatever is the choice of the coordinates we have q k = q k (x a, y a, z a, t) d [ ] L q k = L q k For example consider 2D motion and choose (r, φ) as the generalized coordinates: x = r cos φ, y = r sin φ. Transformation of the coordinate system for the 2nd law of Newton is not a simple task, but using Lagrangian approach it is straightforward to write T = mẋ2 2 + mẏ2 2 = m 2 [r2 φ2 + ṙ 2 ]. Let us assume that the force is f = Ar α 1 /α = r U, then the potential energy U = Ar α. Thus Lagrangian is The corresponding equations of motion are Definitions: p k = L q k is the generalized momentum F k = L q k is the generalized force Thus EL equations are L = m 2 (ṙ2 + r 2 φ2 ) Ar α m( r r φ 2 ) = A α rα 1, m d (r2 φ) = 0 dp k = F k In polar coordinates L/ φ = mr 2 φ = l is the angular momentum and L/ ṙ is the momentum of radial motion. Generalized angular force L/ φ is the torque: dl/ = τ 3.1 Lagrange s equation and motion in a non-inertial frame of reference The power of the Lagrange s approach is that the Lagrange eq. d L L = r r is valid in any coordinate system including the ones, which are not inertial. Inertial frame of reference is the frame, where free motion of a particle, i.e. motion without forces, happens with a constant velocity. It is known that the standard form of the Newton s 2nd law applies only in the inertial frames. Let s assume that in an inertial frame of reference K 0 a particle has the Lagrangian L 0 = 1 2 m v2 0 U
9 9 Naturally the corresponding equation of motion coincides with the 2nd law of Newton m d v 0 = U r 0 First we consider a translational motion with acceleration and then rotational. We consider a non-inertial frame K 1 moving along a straight line with the velocity V relative to K 0 v 0 = v 1 + V (t) We do NOT assume that V is time independent, hence the system K 1 can accelerate. The Lagrange equation in the new system is d L 1 r = L 1 1 r 1 Now we take L 0 and replace v 0 with v 1 to get L 1 L 1 = 1 2 m( v 1) 2 + m v 1 V m V 2 U Before making further transformations we need to note, that if we consider two Lagrangians L a(x, ẋ, t) and L b (x, ẋ, t), such that L b = L a + df/, where f = f(x, t) (does not depend on ẋ!). Then one can easily show that the action A b = A a + f(x 2, t 2) f(x 1, ). Hence δa a = δa b (δx() = δx(t 2) = 0). Hence the resulting equations of motion are the same in the a and b cases. Thus the Lagrangian is defined up to an additive total time derivative of an arbitrary function of co-ordinates and time (but not velocity). Omitting V 2 term as not giving any contribution into the Lagrange equation and replacing v 1 = d r 1 / we find v 1 V = d( V r 1 )/ r 1 d V / Omitting the total time derivative (see above) from the Lagrangian we finally have L 1 = 1 2 m( v 1) 2 m W r 1 U, W = d V / W is the acceleration of the K 1 The equation of motion one can derive is m d v 1 = U r 1 m W Thus an accelerated translational motion of a reference frame is equivalent to an additional force m W. This type of forces resulting from the motion of a reference frame are called inertial forces. Now we want to study effects of rotation and we assume that there is a 3rd reference frame K. Its origin coincides with the origin of K 1, but it rotates with respect to it with the variable angular velocity ω(t). Thus K executes both translational and rotational motions relative to the inertial frame K 0. v 0, r 0 is the velocity and position in K 0, v 1, r 1 in K 1 and v, r in K. Note that r 1 = r, since the origin of these two systems is the same and the vectors go to the same point. Thus we have and v 1 = v + ω r L = 1 2 m( v)2 + m v ( ω r) m( ω r)2 m W r U
10 [note (see the intro to the relativistic chapter), that the rotation of the frame leads to the term depending on the direction of v in the Lagrangian.] Calculation of derivatives needed for the Lagrange eq. is difficult and we first carefully write dl, where L( r, v, t) dl = m v d v + md v ω r + m v ω d r+ 10 +m( ω r) ( ω d r) m W d r ( U/ r) d r = If we organize the equation above in a way that d r and d v enter via scalar products only, then it is easy to get partial derivatives = m v d v + md v ω r + md r v ω+ +m( ω r) ω d r m W d r ( U/ r) d r Thus The Lagrange equation L v = m v + m ω r L r = m v ω + m ω r ω m W ( U/ r) md v/ = ( U/ r) m W + +m r ω + 2m v ω + m ω r ω Thus there are 3 inertial forces due to rotation of a reference frame. The first term is non-zero only if the rotation is non-uniform. The 2nd term is the Coriolis force. The Coriolis force is linear in ω and is the only one which depends on v. The last term is the centrifugal force. Note that acceleration of the frame is equivalent to destruction of the isotropy and homogeneity of space, through the violation of the symmetries of the Lagrangian. 4. Symmetries and conservation laws Noether s theorem (flexible formulation) All fundamental conservation laws of nature can be derived from symmetries. In particular, from the symmetries of space and time, expressed as symmetries of a Lagrangian. 4.1 Conservation of momentum Elementary case. Let s consider 1D Lagrangian L = mẋ 2 /2 + U, d L ẋ = L x
11 11 FIG. 5: Emmy Noether: p = L/ ẋ = p x is the momentum. If U does not depend on x, then dp x / = 0, p x = const which means that momentum is conserved. The generalization of this observation is that if Lagrangian L(q i, q i, t) does not depend explicitly a generalized coordinate q n, so that L/ q n = 0, then the associated momentum L/ q n = p n, is a conserved quantity or an integral of motion so that p n = const. If homogeneity of space is not destroyed by presence of a position dependent potential, then all three components of momentum of a mechanical system are conserved. Homogeneity of space and homogeneity of Lagrangian are the same things. Mathematically it is expressed as: if we take the Lagrangian L( r, r, t) and shift the coordinate system by an arbitrary vector ɛ, i.e. replace r with r + ɛ, ( r r + ɛ), then the Lagrangian has exactly the same appearance in the new coordinate system. Example L = m/2(ẋ 2 + ẏ 2 ). Introducing new coordinates x 0 = x + a, y 0 = y + b ( ɛ = (a, b)) we find that Lagrangian L 0 = m/2(ẋ ẏ0) 2 is the same. If L = m/2(ẋ 2 + ẏ 2 ) mgy, then clearly new Lagrangian will be different. This is because the homogeneity of space along y direction is destroyed by the field of gravity. Proof for the space homogeneous in all 3 directions: Let s consider a small change of the Lagrangian dl caused by a small (time independent) differential change of the radius vector d r (using vector notations and that d r = 0) dl( r, r, t) = L r d r + L r d r = ( i x L + j y L + k z L) d r If Lagrangian is known to be invariant (symmetric) with respect to an arbitrary shift of the reference frame, it means that change of the Lagrangian caused by the shift is zero, i.e., dl = 0, which is possible only if x L = y L = z L = 0 = L r This, through the Lagrange equation, implies, that all 3 components of the momentum are conserved. In a less trivial situation, there can be only a specific direction in space shifts along which preserve the Lagrangian. Let s associate the vector r 0 with this direction, and calculate change of the Lagrangian caused by the shift ɛ r 0, where ɛ 1 dl = L r r 0ɛ
12 Since in this case only the scalar ɛ is arbitrary, we have that L r r 0 = 0, which through the Lagrange equations is equivalent to hence d [ ] L r r 0 = 0 p r 0 = const The above equation implies that the projection of p on r 0 is an integral of motion. Directing a coordinate axis along r 0 will obviously make the momentum component associated with this coordinate to be a conserved quantity. Thus conservation of momentum follows from the homogeneity of space or invariance of a Lagrangian with respect to the coordinate shifts. Other words translational invariance, implies conservation of momentum. All above is equally valid for the many particle case, when we need to sum the momenta of the individual particles, and for example the last equation will be replaced with p a r 0 = const. a Also all above is valid for any generalized coordinates q i different from the Cartesian ones. A generalized coordinate q i, such that L/ q i = 0, is called a cyclic coordinate. Momenta associated with cyclic coordinates do not change with time. Below is the brief summary of the lecture presentation, which is fine, but in future it will be replaced with the more transparent and detailed explanations presented above. Corresponding change in velocities is Q i Q i + ɛk i(q 1, Q 2, Q 3) V i V i + ɛ d Ki(Q1, Q2, Q3) Symmetry implies, that when we make the above substitution in the Lagrangian, we should have that Conservation law is then Worked example: i dl dɛ = 0 dl ( ) L dɛ = Q i Q i ɛ + L V i = V i ɛ i ( ) L K i + L dk i = d L K i = 0 Q i V i V i i L V i K i = const i 12
13 13 6. Lagrangian of a mechanical system is given by Are x and y cyclic coordinates?? Verify that x x + ɛ, y y + 2ɛ is a symmetry. It is obvious by a direct substitution. Work out a conserved momentum: ṗ = 0 L = m 2 (5ẋ2 2ẋẏ + 2ẏ 2 ) + C(2x y) It can be done by couple of equivalent methods. First, lets calculate a change of L, i.e. dl, caused by the nonzero ɛ in the above substitution. It is clear that symmetry implies that dl = 0, simply because the original Lagrangian and Lagrangian with x replaced by x + ɛ and y replaced by y + 2ɛ are the same. On the other hand dl = L(x + ɛ, y + 2ɛ, ẋ, ẏ) L(x, y, ẋ, ẏ) = ɛ( xl + 2 yl). However Lagrange s eqs. are dpx = xl, dp y = yl. Hence dl = ɛ d (px + 2py) = 0. Thus px + 2py = const Another method is to apply r 0 p = const directly. In this example r 0 = (1, 2) and p = (p x, p y). 4.2 Time invariant Lagrangians and conservation of energy If Lagrangian does not explicitly depend on time, i.e. it varies in time only by means of coordinates and velocities, then L L = L(q i, v i ), t = 0, dl 0 Thus in this case the differential change of the Lagrangian dl caused by is given by dl = dl = [ L dq i i q i + L ] dv i = d L v i v i i v i Thus [ ] [ ] d L v i L = 0, H = p i v i L i v i i H is called Hamiltonian or generalized energy. Hamiltonian is a conserved quantity only if Lagrangian does not change after the substitution: t t + ɛ In general H T + U, but in some cases it is. For example, in Cartesian basis i p i v i = 2T and hence H = T + U. 4.3 Conservation of angular momentum Isotropy of space means that mechanical properties of a closed system do not vary, when it is rotated as a whole in space. Formally it means that the Lagrangian is rotationally invariant. If a particle with a radius vector r undergoes rotation by dφ around the vertical axis (specifying direction of d φ and making angle θ with r), then the resulting change of the radius vector is The resulting velocity change is d r = d φ r d v = d φ v
14 14 FIG. 6: The Lagrangian is unchanged on rotation hence dl resulting from d r and d v is 0. dl = a ( L r a d r a + L v a d v a ) = 0 From the Euler-Lagrange equations d L v x = L we have x ( ṗ a d r a + p a dv a ) = ( ṗ a δφ r a + p a dφ v a ) = 0 a a Using permutation of the mixed product A ( B C) = B ( C A) = C ( A B), we find d φ a ( r a ṗ a + v a p a ) = d φ d r a p a = 0 a dφ is arbitrary and the conservation of the total angular momentum around the chosen axis of symmetry is obvious r a p a = const a 5. 1D motion in a potential Consider Lagrangian L = 1 2 mẋ2 U(x) then hamiltonian is H = 1 2 mẋ2 + U(x) ẋ 2 = 2 (H U(x)) m
15 15 FIG. 7: Potential of the harmonic oscillator. Phase plane of the harmonic oscillator. E = H(t = 0) is a constant which can be found from initial conditions. It has meaning of the energy stored in the system at all moments of time including the initial time t = 0. E U(x) for any x. This is simply because there is also kinetic energy If U(x) = E then the corresponding x is a turning point, i.e. the point where velocity of a particle changes its sign going through zero (kinetic energy=0). ± dx 2 = (E U(x)) m t t 0 = t t 0 = ± x 2 x 0 dx m (E U(x )) If there are 2 turning points, then period of oscillations is twice of the time required to go from one turning point to the other T = x2 2m x 1 dx E U(x) the simplest example is the harmonic oscillator with U(x) = kx 2 /2. Classically forbidden regions: E < U(x) Allowed regions: E U(x) 5.1 Pendulum More sophisticated example is the pendulum x = l sin φ, y = l cos φ (we count y from the point of support), U = mgl cos φ, T = ml 2 φ2 /2 EL equation is l φ = g sin φ
16 16 FIG. 8: Small amplitude of oscillations sin φ φ 1 3! φ Then problem is reduced to the standard harmonic oscillator. l φ gφ. Hence ω 2 g/l and period 2π/ω 2π l/g. H = ml φ 2 /2 mgl cos φ 0 < ml 2 φ2 /2 = E + mgl cos φ = E U(φ) Kinetic energy is zero at the turning points φ = φ 0, i.e. E = U(φ 0 ) Therefore, at the turning points, we have E = mgl 2 cos φ 0 Since energy is conserved along the trajectory, we derive ml 2 φ2 /2 = mgl(cos φ cos φ 0 ) If we supply enough energy to the pendulum then oscillatory motion changes to rotation. In the critical case φ 0 = π and E cr = U(π) = mgl. In this case the exact integration of equations of motion is possible. General trigonometric expression tells us that cos φ + 1 = 2 cos 2 φ. Therefore for φ 2 0 = π our equations of motion are reduced to φ 2 = 4g cos 2 φ g l 2, φ = ±2 l cos φ 2 ± corresponds to clockwise or anti-clockwise motion. The ode in question is now 1st order and separable. For convenience we take ψ = φ/2 dψ g cos ψ = l t New integration variable u = tan ψ 1 u2, cos ψ = u, du 2 dψ = 1 2 (1 + u2 )
17 17 FIG. 9: this gives us dψ = 2du/(1 + u 2 ) and 2 g du g 1 u = t 2 l, arctanh(u) = t l (t 2, u = tanh 1 ) g 2 l φ = 4arctan tanh t g l 2 t = 0, then φ = 0, t =, then φ = π It means if we push pendulum with just the correct energy, it will move towards its upright position over the infinite time. If we just below the required energy it will be an oscillatory motion with very long period, if we are just above it will be rotational motion with very long revolution time. 5.2 Stability of equilibria Dynamical system is said to be in an equilibrium if time derivatives of all generalized coordinates and velocities are equal to zero. Equilibrium is stable if small deviations do not drive generalized coordinates far from their values at equilibrium. (these definitions are not very rigorous but sufficient for our purposes) Consider a pendulum l φ = g sin φ There are two equilibrium values of φ: φ 1 = 0, φ 2 = π. φ 1 is stable and φ 2 is unstable. Let s prove it mathematically. Consider a small deviations x from φ 1, i.e. make substitution φ = φ 1 + x lẍ = g sin(φ 1 + x) g[sin φ 1 + (cos φ 1 )x +...] lẍ = gx x = C cos ωt + D sin ωt = Ae iωt + A e iωt = 2 A cos(ωt + ψ 0 ), lω 2 = g
18 18 x is bounded and φ stays close to φ 1 providing C, D are small. Consider small deviations x from φ 2, i.e. make substitution φ = φ 2 + x lẍ = g sin(φ 2 + x) g[sin φ 2 + (cos φ 2 )x +...] lẍ = gx x = Ae λt + Be λt, lλ 2 = g x is increasing with time to very large values and φ does not stay close to φ 2 even if A, B are small. Pendulum falls down towards its stable equilibrium. Using small deviations from equilibria as new generalized coordinates we can work out kinetic and potential energies upto terms which are quadratic in x and ẋ. Write down a quadratic Lagrangian and quadratic Hamiltonian. Derive EL equation which will coincide with above equations. If Lagrangian is quadratic in generalized coordinates and velocities, it means that EL equations derived from this Lagrangian will be linear T = ml2 2 ( φ) 2 ml2 2 (ẋ)2 Taylor expansion for cos, cos φ = cos φ 0 (sin φ 0 )(φ φ 0 ) 1/2(cos φ 0 )(φ φ 0 ) U = mgl cos φ mgl(cos φ 1,2 1 2 cos φ 1,2x 2 ) L(x, ẋ) = ml2 2 ẋ2 mgl 2 cos φ 1,2x 2 generalized coordinate here is x Important Note: terms in the Lagrangian which do not depend on the generalized coordinates and velocities can be disregarded. This is because they do not make any contribution into the EL equation Quadratic Hamiltonian is H(x, ẋ) = ml2 2 ẋ2 + mgl 2 cos φ 1,2x 2 Using this Hamiltonian we can calculate energy of oscillations near the stable equilibrium φ 1 = 0. In this case we have H = ml2 2 ẋ2 + mgl 2 x2 and known solution to the equations of motion is x = 2 A cos(ωt + ψ 0 ). Thus H = 2 A 2 mgl One can ask what is happening with H-quadratic calculated for unstable solution. Can H depend on time for this solution?? Breaking conservation of energy?? and making our approximation not
19 physical?? No! Everything is Ok. H near unstable point φ 2 = π is H = ml2 2 ẋ2 mgl 2 x2 substituting in there our solution x = Ae λt + Be λt we find that H = 2mgl AB Note if initial conditions are such that either A = 0 or B = 0, then H = Problems: One-star problems are for everyone. 2- and 3-star problems are only for enthusiasts. 1. Consider functional A[x(t)] = t 2 L(x, ẋ, ẍ, t). Assume that δx(,2 ) = δẋ(,2 ) = 0 and work out δa δx =?? 2. Find extremum of A[x(t)] = 1 0 (x2 + ẋ 2 + 2xe t ) for x(0) = 0, x(1) = e Answer is x 0 (t) = e 2 2(e 2 1) (et e t ) + tet 2 3. Mass on a vertical spring. Consider mass m oscillating on a vertical spring in the field of gravity. Sketch plot showing your coordinate system. Write down expressions for kinetic energy, potential energy, Lagrangian and Euler-Lagrange equations. Make sure that the EL equations coincide with the corresponding 2nd law of Newton. 4. Lagrangian of a particle is L = m 2 (ẋ2 +ẏ 2 +ż 2 )+mgy. Specify cyclic coordinates and conserved momenta 5. Work out 3D Lagrangian of a free particle in polar coordinates (r, φ, z), x = r cos φ, y = r sin φ. Specify cyclic coordinates and write down expression for the corresponding conserved momenta. 6. Lagrangian of a mechanical system is given by L = m 2 (5ẋ2 2ẋẏ + 2ẏ 2 ) + C(2x y) Are x and y cyclic coordinates?? Work out two generalized momenta: p x and p y. Are they conserved?? Verify that x x + ɛ, y y + 2ɛ is a symmetry. Workout a conserved momentum and verify its conservation by a direct calculation: ṗ = 0 7. A particle of mass m moves in one dimension under the influence of a force F : F = k x 2 e t/τ where k and τ are positive constants. Compute the Lagrangian and Hamiltonian functions. Compare the Hamiltonian and the total energy and discuss the conservation of energy for the system.
20 8. Consider particle with kinetic energy T = ẋ 2 and potential energy U = x 2 x 4. Specify intervals of E values which allow oscillatory motion of the particle. E is the total (initial) energy. Specify stable and unstable equilibria. 9. Work out an approximate expression for dependence of the pendulum frequency on amplitude of oscillations. Plan: Start from l φ = g sin φ and assume that φ is small. Replace sin φ with its Taylor expansion, but retain the first term nonlinear in φ. To solve the resulting equation assume φ = a(t)e iω0t + a(t) e iω0t, where ω 0 is the frequency of the pendulum found in the approximation sin φ φ. Assume that a is a slow function of time: ä ω 0 ȧ. 2nd part. Higher order harmonics. Consider substitution x = a(t)e iω0t + a(t) e iω0t + b(t)e i3ω0t + b(t) e i3ω0t and using the same approximations as above derive the system of first order differential equations for a and b. Using amplitude-phase substitutions a = r 1 e iθ and b = r 3 e i3θ find solutions of this system Write Lagrangian for the pendulum if its hanging point oscillates horizontally as f cos Ωt. Solve the corresponding Lagrange equation in the limit of small oscillations. 11. Write Lagrangian for the pendulum if its hanging point oscillates vertically as f cos Ωt. Write the corresponding Lagrange equation in the limit of small oscillations around φ = 0 and φ = π.
Classical Mechanics Comprehensive Exam Solution
Classical Mechanics Comprehensive Exam Solution January 31, 011, 1:00 pm 5:pm Solve the following six problems. In the following problems, e x, e y, and e z are unit vectors in the x, y, and z directions,
More informationPhysics 5153 Classical Mechanics. Canonical Transformations-1
1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant
More informationSymmetries 2 - Rotations in Space
Symmetries 2 - Rotations in Space This symmetry is about the isotropy of space, i.e. space is the same in all orientations. Thus, if we continuously rotated an entire system in space, we expect the system
More information28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)
28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.
More informationPhysics 235 Chapter 7. Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics
Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics Many interesting physics systems describe systems of particles on which many forces are acting. Some of these forces are immediately
More informationCP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017
CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS Prof. N. Harnew University of Oxford TT 2017 1 OUTLINE : CP1 REVISION LECTURE 3 : INTRODUCTION TO CLASSICAL MECHANICS 1. Angular velocity and
More informationLecture 4. Alexey Boyarsky. October 6, 2015
Lecture 4 Alexey Boyarsky October 6, 2015 1 Conservation laws and symmetries 1.1 Ignorable Coordinates During the motion of a mechanical system, the 2s quantities q i and q i, (i = 1, 2,..., s) which specify
More informationVariation Principle in Mechanics
Section 2 Variation Principle in Mechanics Hamilton s Principle: Every mechanical system is characterized by a Lagrangian, L(q i, q i, t) or L(q, q, t) in brief, and the motion of he system is such that
More informationLagrangian and Hamiltonian Mechanics (Symon Chapter Nine)
Lagrangian and Hamiltonian Mechanics (Symon Chapter Nine Physics A301 Spring 2005 Contents 1 Lagrangian Mechanics 3 1.1 Derivation of the Lagrange Equations...................... 3 1.1.1 Newton s Second
More informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms
More informationPhysics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top
Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More information7 Pendulum. Part II: More complicated situations
MATH 35, by T. Lakoba, University of Vermont 60 7 Pendulum. Part II: More complicated situations In this Lecture, we will pursue two main goals. First, we will take a glimpse at a method of Classical Mechanics
More informationP321(b), Assignement 1
P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates
More informationLecture Notes for PHY 405 Classical Mechanics
Lecture Notes for PHY 405 Classical Mechanics From Thorton & Marion s Classical Mechanics Prepared by Dr. Joseph M. Hahn Saint Mary s University Department of Astronomy & Physics September 1, 2005 Chapter
More informationConsider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be.
Chapter 4 Energy and Stability 4.1 Energy in 1D Consider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be T = 1 2 mẋ2 and the potential energy
More informationis conserved, calculating E both at θ = 0 and θ = π/2 we find that this happens for a value ω = ω given by: 2g
UNIVERSITETET I STAVANGER Institutt for matematikk og naturvitenskap Suggested solutions, FYS 500 Classical Mechanics Theory 2016 fall Set 5 for 23. September 2016 Problem 27: A string can only support
More informationPhysics 351, Spring 2015, Homework #6. Due at start of class, Friday, February 27, 2015
Physics 351, Spring 2015, Homework #6. Due at start of class, Friday, February 27, 2015 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page
More informationM2A2 Problem Sheet 3 - Hamiltonian Mechanics
MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,
More informationTwo-Body Problem. Central Potential. 1D Motion
Two-Body Problem. Central Potential. D Motion The simplest non-trivial dynamical problem is the problem of two particles. The equations of motion read. m r = F 2, () We already know that the center of
More informationKinematics
Classical Mechanics Kinematics Velocities Idea 1 Choose the appropriate frame of reference. Useful frames: - Bodies are at rest - Projections of velocities vanish - Symmetric motion Change back to the
More informationSolution Set Two. 1 Problem #1: Projectile Motion Cartesian Coordinates Polar Coordinates... 3
: Solution Set Two Northwestern University, Classical Mechanics Classical Mechanics, Third Ed.- Goldstein October 7, 2015 Contents 1 Problem #1: Projectile Motion. 2 1.1 Cartesian Coordinates....................................
More information1 Simple Harmonic Oscillator
Massachusetts Institute of Technology MITES 2017 Physics III Lectures 02 and 03: Simple Harmonic Oscillator, Classical Pendulum, and General Oscillations In these notes, we introduce simple harmonic oscillator
More informationPhysics 351, Spring 2015, Final Exam.
Physics 351, Spring 2015, Final Exam. This closed-book exam has (only) 25% weight in your course grade. You can use one sheet of your own hand-written notes. Please show your work on these pages. The back
More informationHamilton-Jacobi theory
Hamilton-Jacobi theory November 9, 04 We conclude with the crowning theorem of Hamiltonian dynamics: a proof that for any Hamiltonian dynamical system there exists a canonical transformation to a set of
More information06. Lagrangian Mechanics II
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 06. Lagrangian Mechanics II Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationPhysical Dynamics (PHY-304)
Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.
More informationHAMILTON S PRINCIPLE
HAMILTON S PRINCIPLE In our previous derivation of Lagrange s equations we started from the Newtonian vector equations of motion and via D Alembert s Principle changed coordinates to generalised coordinates
More informationSketchy Notes on Lagrangian and Hamiltonian Mechanics
Sketchy Notes on Lagrangian and Hamiltonian Mechanics Robert Jones Generalized Coordinates Suppose we have some physical system, like a free particle, a pendulum suspended from another pendulum, or a field
More informationMultibody simulation
Multibody simulation Dynamics of a multibody system (Euler-Lagrange formulation) Dimitar Dimitrov Örebro University June 16, 2012 Main points covered Euler-Lagrange formulation manipulator inertia matrix
More informationthe EL equation for the x coordinate is easily seen to be (exercise)
Physics 6010, Fall 2016 Relevant Sections in Text: 1.3 1.6 Examples After all this formalism it is a good idea to spend some time developing a number of illustrative examples. These examples represent
More informationPHYS2330 Intermediate Mechanics Fall Final Exam Tuesday, 21 Dec 2010
Name: PHYS2330 Intermediate Mechanics Fall 2010 Final Exam Tuesday, 21 Dec 2010 This exam has two parts. Part I has 20 multiple choice questions, worth two points each. Part II consists of six relatively
More informationAnalytical Mechanics - Extra Problems
Analytical Mechanics - Extra Problems Physics 105, F17 (R) are review problems. Review problems are those that have already been covered in prior courses, mostly Intro to Physics I and II. Some are math
More informationLegendre Transforms, Calculus of Varations, and Mechanics Principles
page 437 Appendix C Legendre Transforms, Calculus of Varations, and Mechanics Principles C.1 Legendre Transforms Legendre transforms map functions in a vector space to functions in the dual space. From
More informationPhysics 351, Spring 2017, Homework #12. Due at start of class, Friday, April 14, 2017
Physics 351, Spring 2017, Homework #12. Due at start of class, Friday, April 14, 2017 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page at
More informationLecture 41: Highlights
Lecture 41: Highlights The goal of this lecture is to remind you of some of the key points that we ve covered this semester Note that this is not the complete set of topics that may appear on the final
More informationSolutions 2: Simple Harmonic Oscillator and General Oscillations
Massachusetts Institute of Technology MITES 2017 Physics III Solutions 2: Simple Harmonic Oscillator and General Oscillations Due Wednesday June 21, at 9AM under Rene García s door Preface: This problem
More informationPhysical Dynamics (SPA5304) Lecture Plan 2018
Physical Dynamics (SPA5304) Lecture Plan 2018 The numbers on the left margin are approximate lecture numbers. Items in gray are not covered this year 1 Advanced Review of Newtonian Mechanics 1.1 One Particle
More informationPhys 7221 Homework # 8
Phys 71 Homework # 8 Gabriela González November 15, 6 Derivation 5-6: Torque free symmetric top In a torque free, symmetric top, with I x = I y = I, the angular velocity vector ω in body coordinates with
More informationHamiltonian. March 30, 2013
Hamiltonian March 3, 213 Contents 1 Variational problem as a constrained problem 1 1.1 Differential constaint......................... 1 1.2 Canonic form............................. 2 1.3 Hamiltonian..............................
More informationPhysics 106a, Caltech 13 November, Lecture 13: Action, Hamilton-Jacobi Theory. Action-Angle Variables
Physics 06a, Caltech 3 November, 08 Lecture 3: Action, Hamilton-Jacobi Theory Starred sections are advanced topics for interest and future reference. The unstarred material will not be tested on the final
More informationThe Principle of Least Action
The Principle of Least Action In their never-ending search for general principles, from which various laws of Physics could be derived, physicists, and most notably theoretical physicists, have often made
More informationAnalytical Dynamics: Lagrange s Equation and its Application A Brief Introduction
Analytical Dynamics: Lagrange s Equation and its Application A Brief Introduction D. S. Stutts, Ph.D. Associate Professor of Mechanical Engineering Missouri University of Science and Technology Rolla,
More informationForces of Constraint & Lagrange Multipliers
Lectures 30 April 21, 2006 Written or last updated: April 21, 2006 P442 Analytical Mechanics - II Forces of Constraint & Lagrange Multipliers c Alex R. Dzierba Generalized Coordinates Revisited Consider
More informationPart II. Classical Dynamics. Year
Part II Year 28 27 26 25 24 23 22 21 20 2009 2008 2007 2006 2005 28 Paper 1, Section I 8B Derive Hamilton s equations from an action principle. 22 Consider a two-dimensional phase space with the Hamiltonian
More informationProblem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension
105A Practice Final Solutions March 13, 01 William Kelly Problem 1: Lagrangians and Conserved Quantities Consider the following action for a particle of mass m moving in one dimension S = dtl = mc dt 1
More informationAnalytical Mechanics ( AM )
Analytical Mechanics ( AM ) Olaf Scholten KVI, kamer v8; tel nr 6-55; email: scholten@kvinl Web page: http://wwwkvinl/ scholten Book: Classical Dynamics of Particles and Systems, Stephen T Thornton & Jerry
More informationPhysics 351 Wednesday, February 14, 2018
Physics 351 Wednesday, February 14, 2018 HW4 due Friday. For HW help, Bill is in DRL 3N6 Wed 4 7pm. Grace is in DRL 2C2 Thu 5:30 8:30pm. Respond at pollev.com/phys351 or text PHYS351 to 37607 once to join,
More information4.1 Important Notes on Notation
Chapter 4. Lagrangian Dynamics (Most of the material presented in this chapter is taken from Thornton and Marion, Chap. 7) 4.1 Important Notes on Notation In this chapter, unless otherwise stated, the
More informationPHYSICS 311: Classical Mechanics Final Exam Solution Key (2017)
PHYSICS 311: Classical Mechanics Final Exam Solution Key (017) 1. [5 points] Short Answers (5 points each) (a) In a sentence or two, explain why bicycle wheels are large, with all of the mass at the edge,
More informationFrom quantum to classical statistical mechanics. Polyatomic ideal gas.
From quantum to classical statistical mechanics. Polyatomic ideal gas. Peter Košovan peter.kosovan@natur.cuni.cz Dept. of Physical and Macromolecular Chemistry Lecture 5, Statistical Thermodynamics, MC260P105,
More informationIntroductory Physics. Week 2015/05/29
2015/05/29 Part I Summary of week 6 Summary of week 6 We studied the motion of a projectile under uniform gravity, and constrained rectilinear motion, introducing the concept of constraint force. Then
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationRotational motion problems
Rotational motion problems. (Massive pulley) Masses m and m 2 are connected by a string that runs over a pulley of radius R and moment of inertia I. Find the acceleration of the two masses, as well as
More informationCoordinate systems and vectors in three spatial dimensions
PHYS2796 Introduction to Modern Physics (Spring 2015) Notes on Mathematics Prerequisites Jim Napolitano, Department of Physics, Temple University January 7, 2015 This is a brief summary of material on
More informationResonance and response
Chapter 2 Resonance and response Last updated September 20, 2008 In this section of the course we begin with a very simple system a mass hanging from a spring and see how some remarkable ideas emerge.
More informationMechanics IV: Oscillations
Mechanics IV: Oscillations Chapter 4 of Morin covers oscillations, including damped and driven oscillators in detail. Also see chapter 10 of Kleppner and Kolenkow. For more on normal modes, see any book
More informationPHYSICS 110A : CLASSICAL MECHANICS
PHYSICS 110A : CLASSICAL MECHANICS 1. Introduction to Dynamics motion of a mechanical system equations of motion : Newton s second law ordinary differential equations (ODEs) dynamical systems simple 2.
More informationThe Particle-Field Hamiltonian
The Particle-Field Hamiltonian For a fundamental understanding of the interaction of a particle with the electromagnetic field we need to know the total energy of the system consisting of particle and
More informationLecture 38: Equations of Rigid-Body Motion
Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can
More informationUse conserved quantities to reduce number of variables and the equation of motion (EOM)
Physics 106a, Caltech 5 October, 018 Lecture 8: Central Forces Bound States Today we discuss the Kepler problem of the orbital motion of planets and other objects in the gravitational field of the sun.
More informationPhysics 106b/196b Problem Set 9 Due Jan 19, 2007
Physics 06b/96b Problem Set 9 Due Jan 9, 2007 Version 3: January 8, 2007 This problem set focuses on dynamics in rotating coordinate systems (Section 5.2), with some additional early material on dynamics
More informationCircular motion. Aug. 22, 2017
Circular motion Aug. 22, 2017 Until now, we have been observers to Newtonian physics through inertial reference frames. From our discussion of Newton s laws, these are frames which obey Newton s first
More informationLecture 38: Equations of Rigid-Body Motion
Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can
More informationL = 1 2 a(q) q2 V (q).
Physics 3550 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6, 11.1 11.3 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More information3 Space curvilinear motion, motion in non-inertial frames
3 Space curvilinear motion, motion in non-inertial frames 3.1 In-class problem A rocket of initial mass m i is fired vertically up from earth and accelerates until its fuel is exhausted. The residual mass
More informationPhysics 351, Spring 2017, Homework #3. Due at start of class, Friday, February 3, 2017
Physics 351, Spring 2017, Homework #3. Due at start of class, Friday, February 3, 2017 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page at
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamics in High Energy Accelerators Part 2: Basic tools and concepts Nonlinear Single-Particle Dynamics in High Energy Accelerators This course consists of eight lectures: 1.
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationPhysics 351, Spring 2017, Homework #4. Due at start of class, Friday, February 10, 2017
Physics 351, Spring 2017, Homework #4. Due at start of class, Friday, February 10, 2017 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page
More informationOscillating Inverted Pendulum and Applications
Oscillating Inverted Pendulum and Applications Ramon Driesen, Jaden Johnston, Massimo Pascale, and Evan Ridley Department of Mathematics, University of Arizona, Tucson, Arizona, 85719 (Dated: March 20,
More informationLecture 19: Calculus of Variations II - Lagrangian
Lecture 19: Calculus of Variations II - Lagrangian 1. Key points Lagrangian Euler-Lagrange equation Canonical momentum Variable transformation Maple VariationalCalculus package EulerLagrange 2. Newton's
More informationMathematical Methods of Physics I ChaosBook.org/ predrag/courses/phys Homework 1
PHYS 6124 Handout 6 23 August 2012 Mathematical Methods of Physics I ChaosBook.org/ predrag/courses/phys-6124-12 Homework 1 Prof. P. Goldbart School of Physics Georgia Tech Homework assignments are posted
More informationNoether s Theorem. 4.1 Ignorable Coordinates
4 Noether s Theorem 4.1 Ignorable Coordinates A central recurring theme in mathematical physics is the connection between symmetries and conservation laws, in particular the connection between the symmetries
More informationNumerical Methods for ODEs. Lectures for PSU Summer Programs Xiantao Li
Numerical Methods for ODEs Lectures for PSU Summer Programs Xiantao Li Outline Introduction Some Challenges Numerical methods for ODEs Stiff ODEs Accuracy Constrained dynamics Stability Coarse-graining
More information15. Hamiltonian Mechanics
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 15. Hamiltonian Mechanics Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationChapter 15 Periodic Motion
Chapter 15 Periodic Motion Slide 1-1 Chapter 15 Periodic Motion Concepts Slide 1-2 Section 15.1: Periodic motion and energy Section Goals You will learn to Define the concepts of periodic motion, vibration,
More informationProblem 1 : Solution/marking scheme Two Problems in Mechanics (10 points)
Problem 1 : Solution/marking scheme Two Problems in Mechanics (10 points) Part A. The Hidden Disk (3.5 points) A1 (0.8 pt) Find an expression for b as a function of the quantities (1), the angle φ and
More informationQuestion 1: Spherical Pendulum
Question 1: Spherical Pendulum Consider a two-dimensional pendulum of length l with mass M at its end. It is easiest to use spherical coordinates centered at the pivot since the magnitude of the position
More informationL(q, q) = m 2 q2 V (q) 2 m + V (q)
Lecture 7 Phase Space, Part 1 MATH-GA 71.1 Mechanics 1 Phase portraits 1.1 One dimensional system Consider the generic one dimensional case of a point mass m described by a generalized coordinate q and
More informationPhysics 141, Lecture 7. Outline. Course Information. Course information: Homework set # 3 Exam # 1. Quiz. Continuation of the discussion of Chapter 4.
Physics 141, Lecture 7. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 07, Page 1 Outline. Course information: Homework set # 3 Exam # 1 Quiz. Continuation of the
More informationCurves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,
Physics 6010, Fall 2010 Hamiltonian Formalism: Hamilton s equations. Conservation laws. Reduction. Poisson Brackets. Relevant Sections in Text: 8.1 8.3, 9.5 The Hamiltonian Formalism We now return to formal
More informationNewton s laws. Chapter 1. Not: Quantum Mechanics / Relativistic Mechanics
PHYB54 Revision Chapter 1 Newton s laws Not: Quantum Mechanics / Relativistic Mechanics Isaac Newton 1642-1727 Classical mechanics breaks down if: 1) high speed, v ~ c 2) microscopic/elementary particles
More informationEnergy and Angular Momentum
Chapter 3 Energy and Angular Momentum In this chapter, we generalize the discussion of Chapter 2 to the case of motion in two or three dimensions. Throughout this chapter, we shall be concerned with the
More informationHW 3 Solution Key. Classical Mechanics I HW # 3 Solution Key
Classical Mechanics I HW # 3 Solution Key HW 3 Solution Key 1. (10 points) Suppose a system has a potential energy: U = A(x 2 R 2 ) 2 where A is a constant and R is me fixed distance from the origin (on
More informationLecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More informationElectric and Magnetic Forces in Lagrangian and Hamiltonian Formalism
Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Benjamin Hornberger 1/26/1 Phy 55, Classical Electrodynamics, Prof. Goldhaber Lecture notes from Oct. 26, 21 Lecture held by Prof. Weisberger
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,
More information13. Rigid Body Dynamics II
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 13. Rigid Body Dynamics II Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationNov : Lecture 22: Differential Operators, Harmonic Oscillators
14 Nov. 3 005: Lecture : Differential Operators, Harmonic Oscillators Reading: Kreyszig Sections:.4 (pp:81 83),.5 (pp:83 89),.8 (pp:101 03) Differential Operators The idea of a function as something that
More informationKEELE UNIVERSITY PHYSICS/ASTROPHYSICS MODULE PHY OSCILLATIONS AND WAVES PRACTICE EXAM
KEELE UNIVERSITY PHYSICS/ASTROPHYSICS MODULE PHY-10012 OSCILLATIONS AND WAVES PRACTICE EXAM Candidates should attempt ALL of PARTS A and B, and TWO questions from PART C. PARTS A and B should be answered
More informationPHY 5246: Theoretical Dynamics, Fall September 28 th, 2015 Midterm Exam # 1
Name: SOLUTIONS PHY 5246: Theoretical Dynamics, Fall 2015 September 28 th, 2015 Mierm Exam # 1 Always remember to write full work for what you do. This will help your grade in case of incomplete or wrong
More informationFirst Year Physics: Prelims CP1 Classical Mechanics: DR. Ghassan Yassin
First Year Physics: Prelims CP1 Classical Mechanics: DR. Ghassan Yassin MT 2007 Problems I The problems are divided into two sections: (A) Standard and (B) Harder. The topics are covered in lectures 1
More information2.3 Calculus of variations
2.3 Calculus of variations 2.3.1 Euler-Lagrange equation The action functional S[x(t)] = L(x, ẋ, t) dt (2.3.1) which maps a curve x(t) to a number, can be expanded in a Taylor series { S[x(t) + δx(t)]
More informationSOLUTIONS, PROBLEM SET 11
SOLUTIONS, PROBLEM SET 11 1 In this problem we investigate the Lagrangian formulation of dynamics in a rotating frame. Consider a frame of reference which we will consider to be inertial. Suppose that
More informationNIU PHYS 500, Fall 2006 Classical Mechanics Solutions for HW6. Solutions
NIU PHYS 500, Fall 006 Classical Mechanics Solutions for HW6 Assignment: HW6 [40 points] Assigned: 006/11/10 Due: 006/11/17 Solutions P6.1 [4 + 3 + 3 = 10 points] Consider a particle of mass m moving in
More informationPHY 3221 Fall Homework Problems. Instructor: Yoonseok Lee. Submit only HW s. EX s are additional problems that I encourage you to work on.
PHY 3221 Fall 2012 Homework Problems Instructor: Yoonseok Lee Submit only HW s. EX s are additional problems that I encourage you to work on. Week 1: August 22-24, Due August 27 (nothing to submit) EX:
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture 2.1 Dynamic Behavior Richard M. Murray 6 October 28 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
More informationLecture 6, September 1, 2017
Engineering Mathematics Fall 07 Lecture 6, September, 07 Escape Velocity Suppose we have a planet (or any large near to spherical heavenly body) of radius R and acceleration of gravity at the surface of
More information