- +1-530-264-8006
- info@tutorsglobe.com

18,76,764

Questions

Asked

21,311

Experts

9,67,568

Questions

Answered

Start Excelling in your courses, Ask an Expert and get answers for your homework and assignments!!

Submit Assignment2015 © Tutors Globe. All rights reserved.

## Lagrange and Hamilton Formulation of Mechanics, Physics tutorial

:D'Alembert's PrincipleD'Alembert's Principle: The definition of virtual work was δW = Σ

_{ij}F^{→}_{ij}.δr^{→}_{i}where sum comprises all (constraint and non-constraint) forces. Assuming the position coordinates are in the inertial frame (but not essentially our generalized coordinates), Newton's second law states Σ_{ij}F^{→}_{ij}= P_{i}sum of all forces applying on the particle provide rate of change of momentum. We may then rewrite δW:δW = Σ

_{i}Σ_{j}F^{→}_{ij}.δr^{→}_{i}= Σ_{i}dP_{ i}^{ →}.δr_{ i}^{ →}Thus, relation can be written as

Σ

_{i}[Σ_{j}F^{→}_{ij}^{(nc)}-dP^{→}.δr_{ i}^{ →}= 0,...................................Eq.1Above equation is referred to as D'Alembert's principle. Expression illustrates that rate of change of momentum is determined only by non-constraint forces. In this form, it is not much use, but conclusion that rate of change of momentum is determined only by non-constraint forces is the significant physical statement.

:Generalized Equations of MotionFrom Newtonian mechanics that work is associated to kinetic energy, so it is natural to expect

T ≡ Σ

_{i}1/2m_{i}dr_{i}^{→}_{ }. dr_{ i}^{ →}= T({q_{k}}, {dq_{k}}, t)T must be obtained by first writing T in terms of position velocities {r

^{→}_{i}} and then using definition of position coordinates in terms of generalized coordinates to re-write T as the function of generalized coordinates and velocities.Recalling D'Alembert's principle, equation can be written as:

Σ

_{i}F_{i}^{→(nc)}.∂r_{ i}^{ →}/∂qk = f_{k}= d/dt(∂T/∂dq_{k})-(∂T/∂q_{k})This is generalized equation of motion.

The Lagrangian and the Lagrange's Equations:For conservative non-constraint forces, we can obtain a slightly more compact form of the generalized equation of motion, known as the Euler-Lagrange equations.

Generalized Conservative Forces:Now let us specialize to non-constraint forces which are conservative; i.e. F

^{→}_{i}^{(nc)}= -∇_{i}^{→}U({r_{ j}^{ →}}). Where ∇_{I}^{→}indicates gradient with respect to r_{ i}^{ →}. Whether constraint forces are conservative is unrelated; we will only explicitly require potential for non-constraint forces. U is supposed to be the function of coordinate positions only; there is no explicit dependence on time or on velocities, ∂U/∂t = 0 and ∂U/∂r^{→}_{i}. Let us employ this expression in writing generalized force:f

_{k}= Σ_{i}F_{i}^{→(nc)}.∂r_{ i}^{ →}/∂q_{k}= -Σ_{i}∇_{i}^{→}U({r_{ j}^{ →}}).∂r_{ i}^{ →}/∂q_{k}= ∂/∂q_{k}U({q_{1}},t).We made use of holonomic constraints to re-write U as the function of {q

_{l}} and possibly t, and realize that previous line is just partial derivative of U with respect to q_{k}. Therefore, rather than finding the equation of motion by computing generalized force from non-constraint forces and coordinate transformation relations, we can re-write potential energy as the function of generalized coordinates and compute generalized force by gradients thereof.The Euler-Lagrange Equations:An even simpler method exists; we may re-write the generalized equation of motion using the above relationship between generalized force and gradient of the potential energy as

∂U/∂q

_{k}= d/dt(∂T/∂dq_{k}) - ∂T/∂q_{k}State the Lagrangian as

L ≡ T - U

As we have assumed holonomic constraints, we have ∂U/∂q

_{k}= 0. This permits replacement of d/dt(∂T/∂dq_{k}) with d/dt(∂L/∂dq_{k}) providingd/dt(∂L/∂dq

_{k}) - (∂L/∂q_{k}) = 0This is Euler-Lagrange equation; there is one for each generalized coordinate q

_{k}.:Lagrange's Mechanics:Newton's Second LawWe consider N particles moving in three-dimensional space, and we explain location of each particle using Cartesian coordinates. We let m

_{i}be mass of particle i, and we let x_{i}, y_{i}, and z_{i}be respectively x, y and z-coordinates of particle i. For time derivatives of coordinates (and all other physical observables), we use dot notation first introduced by Isaac Newtonx

^{.}_{i }= dx_{i}/dt and x^{..}_{i}= d^{2}x_{i}/dt^{2}We let Fx

_{i}be x-component of force on particle i. Then, Newton's Second Law takes form Fx_{i}= mx_{i}.For instance, if we study motion of single particle of mass m moving in one dimension in the harmonic potential with related force F

_{x}= -k_{x}. Newton's second law takes form- k_{x}= m_{x}.Conservative Systems:In the qualitative sense conservative systems are those for which total energy E is sum of kinetic and potential energies. For any isolated system E is conserved (i.e. dE/dt = 0), and for conservative systems, sum of potential energy and kinetic energy is conserved. Explicitly, we have following definition:

The classical mechanical system is conservative if there exists the function U(x

_{1}, y_{1}, z_{1}, x_{2},...z_{n}) known as potential energy such that for any coordinate x_{i}(or y_{i}or zi_{)}, we can writeF

_{x}i = ∂U/∂x_{i}(or Fy_{i}= ∂U/∂y_{i}or Fz_{i}= ∂U/∂z_{i})Where Fx

_{i}(or Fy_{i}or Fz_{i}) is x (or y or z) component of force on particle i. As the example, we can consider one dimensional particle moving in harmonic well with force F = -k_{x}. For such a system, potential energy exists and is provided by U(x) = 1/2 k_{x}^{2}By differentiating potential energy with respect to x, force is attained. We can then be sure that for harmonic oscillator total energy is conservedE = 1/2 mdx

^{2}+ 1/2 kx^{2}:Lagrange's EquationsNow derive Lagrange's equations for special case of the conservative system in Cartesian coordinates. Start with equation for kinetic energy. We first distinguish expression with respect to one of velocities

∂T/∂x

_{j}= m_{j}dx_{j}Definition: Classical Lagrangian is provided by

L = T - U

Classical Lagrangian is a difference between kinetic and potential energies of system. Using the definition we get

d/dt(∂L/∂x

_{j}) - ∂L/∂x_{j}= 0The equation given above is Lagrange's equations in Cartesian coordinates. We use plural (equations), as Lagrange's equations are the set of equations. We have separate equation for each coordinate xj. A entirely analogous set of equations is attained for other Cartesian directions y and z.

Generalized Coordinates:Much of Lagrange's work was concerned with methods useful for systems subject to external constraints.

Lagrange's Equations in Generalized Coordinates:Lagrange has shown that the form of Lagrange's equations is invariant to the particular set of generalized coordinates chosen. For any set of generalized coordinates, Lagrange's equations take the form

(d/dt)(∂L/∂qi) = 0

exactly the same form that we derived in Cartesian coordinates. We now illustrate how to use Lagrange's equations in generalized coordinates by applying the approach to the free motion of a particle confined to move on the perimeter of a ring as discussed previously. The meaning of the expression of "free particle" is the absence of any external forces. We can arbitrarily set the potential energy U to zero. Then, in Cartesian coordinates, the Lagrangian for any free particle in the xy-plane can be expressed as

L = 1/2m(dx

^{2}+d y^{2})Because of the constraint, Lagrangian is the function of single coordinate f. We finally give Lagrange's equations

∂L/∂Φ = mR

^{2}Φ, and d/dt(∂L/∂Φ) = mR2d2Φ and (∂L/∂Φ)= 0So therefore

d/dt(∂L/∂Φ) = d/dtmR

^{2}dΦ/dt = 0Equation given above means acceleration of coordinates Φ is zero so that particle moves with the constant generalized velocity Φ.

Generalized Momenta:Equation can be interpreted to mean that quantity ∂L/∂Φ = mR

^{2}Φ is conserved. To fully explore meaning of conservation of a quantity like ∂L/∂Φ, consider the Lagrangian in Cartesian coordinates for the particle of mass in one dimension.L = 1/2mdx

^{2}- U(x)By differentiating L with respect to velocity dx, we get linear momentum.

∂L/∂x = mdx

Which is conserved in case of no external forces; that is linear momentum is conserved if U(x) is the constant. Using this simple equation, we are lead to given definition:

Generalized momentum pi conjugate to coordinate q

_{i}is stated byp

_{i }= ∂L/∂q_{i}Lagrangian for Some Physical Systems:i)

1-D motion-the pendulum:One of the simplest nonlinear systems is one-dimensional physical pendulum (so called to differentiate it from linearized harmonic oscillator approximation). The pendulum comprises of the light rigid rod of length l, making the angle Θ with vertical, swinging from the fixed pivot at one end and with the bob of mass m attached at the other end. This system is also known as simple pendulum to distinguish it from spherical pendulum and compound pendula that have more than one degree of freedom.

Potential energy with respect to equilibrium position Θ = 0 is U(Θ) = mgl(1- cosΘ ), where g is acceleration due to gravity, and velocity of bob is v

_{Θ}= lΘ, So that kinetic energy T= 1/2mv_{Θ}^{2}= 1/2ml2Θ2The Lagrangian, L= T - U, is thus

L (Θ, dΘ) = 1/2ml

^{2}dΘ^{2}- mgl(1 - cosΘ)This is also effectively the Lagrangian for particle moving in the sinusoidal spatial potential, so physical pendulum gives a paradigm for problems like motion of the electron in the crystal lattice or of the ion or electron in a plasma wave.

Lagrangian equation of motion is

d/dt(∂L/∂Θ) where ∂L/∂Θ = ml

^{2}Θ and ∂L/∂Θ = -mglsinΘHence Lagrangian equation of motion is ml

^{2}d^{2}Θ/d Θ = -mglsinΘ.ii)

2-D motion with time-varying constraint:Consider the weight rotating about origin on the frictionless horizontal surface and constrained by the thread, initially of length a, that is being pulled steadily downward at speed u through the hole at origin so that radius r = a - ut.

The Lagrangian is L = T= 1/2m [u

^{2}+ (a -ut)^{2}Θ^{2}]So we again have conservation of angular momentum m(a -ut)

^{2}dΘ = l = constantThis can be integrated to provide Θ as the function of t,

Θ = Θ

_{0}+ (l/mu)[1/(a - ut) - 1/a] = Θ_{0}+lt/[ma(a - ut)]iii)

Atwood Machine:Consider two weights of mass m

_{1}and m_{2 }suspended from a frictionless, inertia less pulley of radius a by a rope of fixed length. The height of weight 1 is x with respect to the chosen origin and the holonomic constraint given by rope lets us to express height of weight 2 as -x, so that there is only one degree of freedom for this system. Kinetic and potential energy areT = 1/2(m

_{1}+ m_{2})dx^{2}/dx and U = m_{1}gx - m_{2}gxLagrange's equation of motion now becomes

(m

_{1}+ m_{2})d^{2}x/dx = -(m_{1}- m_{2})gd

^{2}x/dx = [(m_{1}- m_{2})/(m_{1}+ m_{2})]giv) Particle in E.M. field:The fact that Lagrange's equations are the Euler-Lagrange equations for the extraordinarily simple and general Hamilton's principle suggests that Lagrange's equations of motion may have a wider range of validity than simply problems where the force is derivable from a scalar potential.

The equation of motion of a charged particle in an electromagnetic field, under the influence of the Lorentz force,

md

^{2}r/dr = eE(r, t) + edr X B(r,t)Where e is the charge on the particle of mass m. We assume the electric and magnetic fields E and B, respectively, to be given in terms of the scalar potential Φ and vector potential A by the standard relations

E = -∇Φ - ∂

_{t}A, and B = ∇ x ALagrange's equation of motion then becomes

md

^{2}q/dq = e[(-∂Φ/∂q_{i}) - (∂A_{i}/∂t)] + eΣ_{j=1}^{e}dq_{j}[(∂A_{j}/∂q_{i})(∂A_{i}/∂q_{j})]:Transformations of the LagrangianPoint transformations:Given the arbitrary Lagrangian L(q, dq, t) in one generalized coordinate system, q ≡ {q

_{i}|i = 1,n} (e.g. a Cartesian frame), we often want to know the Lagrangian L'(Q, dQ, t ) in another coordinate system, Q ≡ {Q_{i}|i =1,n}, (e.g. polar coordinates). Thus, suppose there exists a set g ≡ {g_{i}|i = 1,n} of twice differentiable functions g_{i}such that q_{i}g_{i}(Q,t ), i = 1,...,n.We need inverse function of g also to be twice differentiable, in which case g: Q → q is said to be C

^{2 }diffeomorphism. Transformation g maps a path in Q-space to the path in q-space. Though, it is physically same path; all we have changed is representation. What we need to discuss how Lagrangian transforms is the coordinate-free formulation of Lagrangian dynamics. Therefore, if we can stateL' so that L'(Q(t), dQ(t), t) = L(q(t), aq(t), t) for any path, then S is automatically invariant under coordinate change and will be stationary for same physical paths, irrespective of what coordinates they are represented in. One can prove that equation provides correct dynamics by computing transformation of Euler-Lagrange equations explicitly and showing that equation signifies

(d/dt)(∂L'/∂Q

_{i}) - (∂L'/∂Q_{i}) = 0Gauge transformations:The Lagrangian gives the function for which Lagrange's equations represent correct dynamical equations of motion. We generally normalize L in the natural way, example by requiring that part linear in mass be equivalent to the kinetic energy, hence this freedom is not encountered much in practice. Though, there is a more significant source of non-uniqueness, called as a gauge transformation of Lagrangian in which L is replaced by L', stated by

L'(q, dq, t) = L(q, dq, t) + ∂M/∂t + dq.∂M/∂q

Now add the gauge term, taking M = 1/2ω

_{0}^{2}x. Then, from equation, new Lagrangian isL' = 1/2(dx

^{2}+2ωxdx - ω_{0}^{2}x^{2})Hamiltonian mechanics:Legendre Transform

The Legendre transform is a method of changing the dependence of a function of one set of variables to another set of variables. In mathematics, it is often desirable to express a functional relationship f(x) as a different function, whose argument is the derivative of f, rather than x. If we let p = df/dx be argument of this new function, then new function is written f *(p) and is known as Legendre transform of original function, named after Adrien-Marie Legendre. Legendre transform f *of the function f is stated as follows: f *( p) = max

_{x }( px - f (x)). Legendre transform is its own inverse. Like the familiar Fourier transform, Legendre transform takes the function f(x) and produces the function of different variable p. Though, while Fourier transform comprises of the integration with the kernel, Legendre transform uses maximization as transformation procedure. Legendre transformation can be generalized to Legendre-Fenchel transformation. It is usually utilized in thermodynamics and in Hamiltonian formulation of classical mechanics.Total differential of Helmholtz free energy is provided by

dA = -dT - PdV.

Where S is entropy, T is temperature, V is volume and p is pressure. From the total differential, it is clear that Helmholtz free energy is stated as the function of temperature and volume; i.e. A = A(T, V ).

:The Classical Hamiltonian and Hamilton's EquationsNow apply the notion of the Legendre transform to the classical Lagrangian. A function of all the generalized coordinates and their respective time derivatives; i.e. L = ({q

_{i}}, {dq_{i}}t). For generality, also include possibility that Lagrangian has explicit time dependence. Such explicit time dependence can happen when external forces applying on the system are time-dependent.Resulting time-dependent potentials can be significant in systems, as for instance, study of interaction of radiation with matter. Light is made up of electric and magnetic fields which oscillate in time, and when light interacts with matter, electrons are subjected to time dependent potentials. For system of particles each having masses m

_{i}described by the set of generalized coordinates q_{i}, the classical Hamiltonian is stated byH = Σ

_{i}P_{i}dq_{i}- L({q_{i}}, {dq_{i}}, t).Hamiltonian is then seen to be the expression for total energy of the conservative system in terms of the generalized coordinates and momenta. With the Hamiltonian expressed in terms of the proper variables, Hamilton's equations of the system as

∂H/∂P

_{x}= P_{x}/m = dx and ∂H/∂x = ∂U/∂x = -Px:Construction of the Hamiltonian in Spherical Polar Coordinates - Central Force MotionFor systems with spherical symmetry (e.g. the rigid rotator), spherical polar coordinates are the most convenient set of generalized coordinates. As depicted above, the spherical polar coordinates are r, Θ and Φ. The coordinate r is the distance from the origin of coordinates to the particle, Θ is the angle the line connecting the origin of the coordinates to the particle position makes with the z-axis and Φ is the angle the projection of the line defining Θ onto the xy-plane makes with the x-axis. The connections between Cartesian coordinates and spherical polar coordinates can be derived readily using trigonometry. The result is

x = rsinΘcosΦ y= rsinΘcosΦ z = rcosΘ

Another significant relation is direct result of Pythagoras theorem

r = √x

^{2}+ y^{2}+ z^{2}Using definition of classical Hamiltonian

= 1/2m [dr

^{2}+ r^{2}dΘ^{2}+ r^{2}sin^{2}ΘΦ^{2}]+ U(r)Finally equation should be transformed to replace velocities with generalized momenta. Lastly get

H = [P

^{2}_{r}/2m + P^{2}_{Θ}/2mr^{2}+ P^{2}_{Φ}/2mr^{2}sin^{2}Θ] +U(r)If needed equations of motion can be attained by applying Hamilton's equations to constructed Hamiltonian.

Tutorsglobe: A way to secure high grade in your curriculum (Online Tutoring)Expand your confidence, grow study skills and improve your grades.

Since 2009, Tutorsglobe has proactively helped millions of students to get better grades in school, college or university and score well in competitive tests with live, one-on-one online tutoring.

Using an advanced developed tutoring system providing little or no wait time, the students are connected on-demand with a tutor at www.tutorsglobe.com. Students work one-on-one, in real-time with a tutor, communicating and studying using a virtual whiteboard technology. Scientific and mathematical notation, symbols, geometric figures, graphing and freehand drawing can be rendered quickly and easily in the advanced whiteboard.

Free to know our price and packages for online physics tutoring. Chat with us or submit request at info@tutorsglobe.com