Hamilton's principle





























In physics, Hamilton's principle is William Rowan Hamilton's formulation of the principle of stationary action (see that article for historical formulations). It states that the dynamics of a physical system is determined by a variational problem for a functional based on a single function, the Lagrangian, which contains all physical information concerning the system and the forces acting on it. The variational problem is equivalent to and allows for the derivation of the differential equations of motion of the physical system. Although formulated originally for classical mechanics, Hamilton's principle also applies to classical fields such as the electromagnetic and gravitational fields, and plays an important role in quantum mechanics, quantum field theory and criticality theories.




As the system evolves, q traces a path through configuration space (only some are shown). The path taken by the system (red) has a stationary action (δS = 0) under small changes in the configuration of the system (δq).[1]




Contents






  • 1 Mathematical formulation


    • 1.1 Euler–Lagrange equations derived from the action integral


    • 1.2 Canonical momenta and constants of motion


    • 1.3 Example: Free particle in polar coordinates




  • 2 Applied to deformable bodies


  • 3 Comparison with Maupertuis' principle


  • 4 Action principle for fields


    • 4.1 Classical field theory


    • 4.2 Quantum mechanics and quantum field theory




  • 5 See also


  • 6 References





Mathematical formulation


Hamilton's principle states that the true evolution q(t) of a system described by N generalized coordinates q = (q1, q2, ..., qN) between two specified states q1 = q(t1) and q2 = q(t2) at two specified times t1 and t2 is a stationary point (a point where the variation is zero), of the action functional


S[q] =def ∫t1t2L(q(t),q˙(t),t)dt{displaystyle {mathcal {S}}[mathbf {q} ] {stackrel {mathrm {def} }{=}} int _{t_{1}}^{t_{2}}L(mathbf {q} (t),{dot {mathbf {q} }}(t),t),dt}{mathcal  {S}}[{mathbf  {q}}] {stackrel  {{mathrm  {def}}}{=}} int _{{t_{1}}}^{{t_{2}}}L({mathbf  {q}}(t),{dot  {{mathbf  {q}}}}(t),t),dt

where L(q,q˙,t){displaystyle L(mathbf {q} ,{dot {mathbf {q} }},t)}L({mathbf  {q}},{dot  {{mathbf  {q}}}},t) is the Lagrangian function for the system. In other words, any first-order perturbation of the true evolution results in (at most) second-order changes in S{displaystyle {mathcal {S}}}{mathcal {S}}. The action S{displaystyle {mathcal {S}}}{mathcal {S}} is a functional, i.e., something that takes as its input a function and returns a single number, a scalar. In terms of functional analysis, Hamilton's principle states that the true evolution of a physical system is a solution of the functional equation



Hamilton's principle

δq(t)=0{displaystyle {frac {delta {mathcal {S}}}{delta mathbf {q} (t)}}=0}{frac  {delta {mathcal  {S}}}{delta {mathbf  {q}}(t)}}=0




Euler–Lagrange equations derived from the action integral


Requiring that the true trajectory q(t) be a stationary point of the action functional S{displaystyle {mathcal {S}}}{mathcal {S}} is equivalent to a set of differential equations for q(t) (the Euler–Lagrange equations), which may be derived as follows.


Let q(t) represent the true evolution of the system between two specified states q1 = q(t1) and q2 = q(t2) at two specified times t1 and t2, and let ε(t) be a small perturbation that is zero at the endpoints of the trajectory


ε(t1)=ε(t2) =def 0{displaystyle {boldsymbol {varepsilon }}(t_{1})={boldsymbol {varepsilon }}(t_{2}) {stackrel {mathrm {def} }{=}} 0}{boldsymbol  varepsilon }(t_{1})={boldsymbol  varepsilon }(t_{2}) {stackrel  {{mathrm  {def}}}{=}} 0

To first order in the perturbation ε(t), the change in the action functional δS{displaystyle delta {mathcal {S}}}delta {mathcal  {S}} would be


δS=∫t1t2[L(q+ε,q˙˙)−L(q,q˙)]dt=∫t1t2(εL∂q+ε˙L∂)dt{displaystyle delta {mathcal {S}}=int _{t_{1}}^{t_{2}};left[L(mathbf {q} +{boldsymbol {varepsilon }},{dot {mathbf {q} }}+{dot {boldsymbol {varepsilon }}})-L(mathbf {q} ,{dot {mathbf {q} }})right]dt=int _{t_{1}}^{t_{2}};left({boldsymbol {varepsilon }}cdot {frac {partial L}{partial mathbf {q} }}+{dot {boldsymbol {varepsilon }}}cdot {frac {partial L}{partial {dot {mathbf {q} }}}}right),dt}delta {mathcal  {S}}=int _{{t_{1}}}^{{t_{2}}};left[L({mathbf  {q}}+{boldsymbol  varepsilon },{dot  {{mathbf  {q}}}}+{dot  {{boldsymbol  {varepsilon }}}})-L({mathbf  {q}},{dot  {{mathbf  {q}}}})right]dt=int _{{t_{1}}}^{{t_{2}}};left({boldsymbol  varepsilon }cdot {frac  {partial L}{partial {mathbf  {q}}}}+{dot  {{boldsymbol  {varepsilon }}}}cdot {frac  {partial L}{partial {dot  {{mathbf  {q}}}}}}right),dt

where we have expanded the Lagrangian L to first order in the perturbation ε(t).


Applying integration by parts to the last term results in


δS=[εL∂]t1t2+∫t1t2(εL∂q−εddt∂L∂)dt{displaystyle delta {mathcal {S}}=left[{boldsymbol {varepsilon }}cdot {frac {partial L}{partial {dot {mathbf {q} }}}}right]_{t_{1}}^{t_{2}}+int _{t_{1}}^{t_{2}};left({boldsymbol {varepsilon }}cdot {frac {partial L}{partial mathbf {q} }}-{boldsymbol {varepsilon }}cdot {frac {d}{dt}}{frac {partial L}{partial {dot {mathbf {q} }}}}right),dt}delta {mathcal  {S}}=left[{boldsymbol  varepsilon }cdot {frac  {partial L}{partial {dot  {{mathbf  {q}}}}}}right]_{{t_{1}}}^{{t_{2}}}+int _{{t_{1}}}^{{t_{2}}};left({boldsymbol  varepsilon }cdot {frac  {partial L}{partial {mathbf  {q}}}}-{boldsymbol  varepsilon }cdot {frac  {d}{dt}}{frac  {partial L}{partial {dot  {{mathbf  {q}}}}}}right),dt

The boundary conditions ε(t1)=ε(t2) =def 0{displaystyle {boldsymbol {varepsilon }}(t_{1})={boldsymbol {varepsilon }}(t_{2}) {stackrel {mathrm {def} }{=}} 0}{boldsymbol  varepsilon }(t_{1})={boldsymbol  varepsilon }(t_{2}) {stackrel  {{mathrm  {def}}}{=}} 0 causes the first term to vanish


δS=∫t1t2ε(∂L∂q−ddt∂L∂)dt{displaystyle delta {mathcal {S}}=int _{t_{1}}^{t_{2}};{boldsymbol {varepsilon }}cdot left({frac {partial L}{partial mathbf {q} }}-{frac {d}{dt}}{frac {partial L}{partial {dot {mathbf {q} }}}}right),dt}delta {mathcal  {S}}=int _{{t_{1}}}^{{t_{2}}};{boldsymbol  varepsilon }cdot left({frac  {partial L}{partial {mathbf  {q}}}}-{frac  {d}{dt}}{frac  {partial L}{partial {dot  {{mathbf  {q}}}}}}right),dt

Hamilton's principle requires that this first-order change δS{displaystyle delta {mathcal {S}}}delta {mathcal  {S}} is zero for all possible perturbations ε(t), i.e., the true path is a stationary point of the action functional S{displaystyle {mathcal {S}}}{mathcal {S}} (either a minimum, maximum or saddle point). This requirement can be satisfied if and only if



Euler–Lagrange equations

L∂q−ddt∂L∂=0{displaystyle {frac {partial L}{partial mathbf {q} }}-{frac {d}{dt}}{frac {partial L}{partial {dot {mathbf {q} }}}}=0}{frac  {partial L}{partial {mathbf  {q}}}}-{frac  {d}{dt}}{frac  {partial L}{partial {dot  {{mathbf  {q}}}}}}=0



These equations are called the Euler–Lagrange equations for the variational problem.



Canonical momenta and constants of motion


The conjugate momentum pk for a generalized coordinate qk is defined by the equation



pk =def ∂L∂k{displaystyle p_{k} {stackrel {mathrm {def} }{=}} {frac {partial L}{partial {dot {q}}_{k}}}}p_{k} {stackrel  {{mathrm  {def}}}{=}} {frac  {partial L}{partial {dot  {q}}_{k}}}.

An important special case of the Euler–Lagrange equation occurs when L does not contain a generalized coordinate qk explicitly,


L∂qk=0⇒ddt∂L∂k=0⇒dpkdt=0,{displaystyle {frac {partial L}{partial q_{k}}}=0quad Rightarrow quad {frac {d}{dt}}{frac {partial L}{partial {dot {q}}_{k}}}=0quad Rightarrow quad {frac {dp_{k}}{dt}}=0,,}{frac  {partial L}{partial q_{k}}}=0quad Rightarrow quad {frac  {d}{dt}}{frac  {partial L}{partial {dot  {q}}_{k}}}=0quad Rightarrow quad {frac  {dp_{k}}{dt}}=0,,

that is, the conjugate momentum is a constant of the motion.


In such cases, the coordinate qk is called a cyclic coordinate. For example, if we use polar coordinates t, r, θ to describe the planar motion of a particle, and if L does not depend on θ, the conjugate momentum is the conserved angular momentum.



Example: Free particle in polar coordinates


Trivial examples help to appreciate the use of the action principle via the Euler–Lagrange equations. A free particle (mass m and velocity v) in Euclidean space moves in a straight line. Using the Euler–Lagrange equations, this can be shown in polar coordinates as follows. In the absence of a potential, the Lagrangian is simply equal to the kinetic energy


L=12mv2=12m(x˙2+y˙2){displaystyle L={frac {1}{2}}mv^{2}={frac {1}{2}}mleft({dot {x}}^{2}+{dot {y}}^{2}right)}L={frac  {1}{2}}mv^{2}={frac  {1}{2}}mleft({dot  {x}}^{2}+{dot  {y}}^{2}right)

in orthonormal (x,y) coordinates, where the dot represents differentiation with respect to the curve parameter (usually the time, t). Therefore, upon application of the Euler–Lagrange equations,


ddt(∂L∂)−L∂x=0⇒mx¨=0{displaystyle {frac {d}{dt}}left({frac {partial L}{partial {dot {x}}}}right)-{frac {partial L}{partial x}}=0qquad Rightarrow qquad m{ddot {x}}=0}{frac  {d}{dt}}left({frac  {partial L}{partial {dot  {x}}}}right)-{frac  {partial L}{partial x}}=0qquad Rightarrow qquad m{ddot  {x}}=0

And likewise for y. Thus the Euler–Lagrange formulation can be used to derive Newton's laws.


In polar coordinates (r, φ) the kinetic energy and hence the Lagrangian becomes


L=12m(r˙2+r2φ˙2).{displaystyle L={frac {1}{2}}mleft({dot {r}}^{2}+r^{2}{dot {varphi }}^{2}right).}L={frac  {1}{2}}mleft({dot  {r}}^{2}+r^{2}{dot  {varphi }}^{2}right).

The radial r and φ components of the Euler–Lagrange equations become, respectively


ddt(∂L∂)−L∂r=0⇒˙2=0{displaystyle {frac {d}{dt}}left({frac {partial L}{partial {dot {r}}}}right)-{frac {partial L}{partial r}}=0qquad Rightarrow qquad {ddot {r}}-r{dot {varphi }}^{2}=0}{frac  {d}{dt}}left({frac  {partial L}{partial {dot  {r}}}}right)-{frac  {partial L}{partial r}}=0qquad Rightarrow qquad {ddot  {r}}-r{dot  {varphi }}^{2}=0

ddt(∂L∂φ˙)−L∂φ=0⇒φ¨+2rr˙φ˙=0.{displaystyle {frac {d}{dt}}left({frac {partial L}{partial {dot {varphi }}}}right)-{frac {partial L}{partial varphi }}=0qquad Rightarrow qquad {ddot {varphi }}+{frac {2}{r}}{dot {r}}{dot {varphi }}=0.}{frac  {d}{dt}}left({frac  {partial L}{partial {dot  {varphi }}}}right)-{frac  {partial L}{partial varphi }}=0qquad Rightarrow qquad {ddot  {varphi }}+{frac  {2}{r}}{dot  {r}}{dot  {varphi }}=0.

The solution of these two equations is given by



r=(at+b)2+c2{displaystyle r={sqrt {(at+b)^{2}+c^{2}}}}r={sqrt  {(at+b)^{2}+c^{2}}}

φ=tan−1⁡(at+bc)+d{displaystyle varphi =tan ^{-1}left({frac {at+b}{c}}right)+d}varphi =tan ^{{-1}}left({frac  {at+b}{c}}right)+d


for a set of constants a, b, c, d determined by initial conditions.
Thus, indeed, the solution is a straight line given in polar coordinates: a is the velocity, c is the distance of the closest approach to the origin, and d is the angle of motion.



Applied to deformable bodies


Hamilton's principle is an important variational principle in elastodynamics. As opposed to a system composed of rigid bodies, deformable bodies have an infinite number of degrees of freedom and occupy continuous regions of space; consequently, the state of the system is described by using continuous functions of space and time. The extended Hamilton Principle for such bodies is given by


t1t2[δWe+δT−δU]dt=0{displaystyle int _{t_{1}}^{t_{2}}left[delta W_{e}+delta T-delta Uright]dt=0}{displaystyle int _{t_{1}}^{t_{2}}left[delta W_{e}+delta T-delta Uright]dt=0}

where T is the kinetic energy, U is the elastic energy, We is the work done by
external loads on the body, and t1, t2 the initial and final times. If the system is conservative, the work done by external forces may be derived from a scalar potential V. In this case,


δt1t2[T−(U+V)]dt=0.{displaystyle delta int _{t_{1}}^{t_{2}}left[T-(U+V)right]dt=0.}{displaystyle delta int _{t_{1}}^{t_{2}}left[T-(U+V)right]dt=0.}

This is called Hamilton's principle and it is invariant under coordinate transformations.



Comparison with Maupertuis' principle


Hamilton's principle and Maupertuis' principle are occasionally confused and both have been called (incorrectly) the principle of least action. They differ in three important ways:


  • their definition of the action...


Maupertuis' principle uses an integral over the generalized coordinates known as the abbreviated action or reduced action
S0 =def ∫p⋅dq{displaystyle {mathcal {S}}_{0} {stackrel {mathrm {def} }{=}} int mathbf {p} cdot dmathbf {q} }{mathcal  {S}}_{{0}} {stackrel  {{mathrm  {def}}}{=}} int {mathbf  {p}}cdot d{mathbf  {q}}


where p = (p1, p2, ..., pN) are the conjugate momenta defined above. By contrast, Hamilton's principle uses S{displaystyle {mathcal {S}}}{mathcal {S}}, the integral of the Lagrangian over time.


  • the solution that they determine...

Hamilton's principle determines the trajectory q(t) as a function of time, whereas Maupertuis' principle determines only the shape of the trajectory in the generalized coordinates. For example, Maupertuis' principle determines the shape of the ellipse on which a particle moves under the influence of an inverse-square central force such as gravity, but does not describe per se how the particle moves along that trajectory. (However, this time parameterization may be determined from the trajectory itself in subsequent calculations using the conservation of energy). By contrast, Hamilton's principle directly specifies the motion along the ellipse as a function of time.

  • ...and the constraints on the variation.

Maupertuis' principle requires that the two endpoint states q1 and q2 be given and that energy be conserved along every trajectory (same energy for each trajectory). This forces the endpoint times to be varied as well. By contrast, Hamilton's principle does not require the conservation of energy, but does require that the endpoint times t1 and t2 be specified as well as the endpoint states q1 and q2.


Action principle for fields



Classical field theory



The action principle can be extended to obtain the equations of motion for fields, such as the electromagnetic field or gravity.


The Einstein equation utilizes the Einstein–Hilbert action as constrained by a variational principle.


The path of a body in a gravitational field (i.e. free fall in space time, a so-called geodesic) can be found using the action principle.



Quantum mechanics and quantum field theory



In quantum mechanics, the system does not follow a single path whose action is stationary, but the behavior of the system depends on all imaginable paths and the value of their action. The action corresponding to the various paths is used to calculate the path integral, that gives the probability amplitudes of the various outcomes.


Although equivalent in classical mechanics with Newton's laws, the action principle is better suited for generalizations and plays an important role in modern physics. Indeed, this principle is one of the great generalizations in physical science. In particular, it is fully appreciated and best understood within quantum mechanics. Richard Feynman's path integral formulation of quantum mechanics is based on a stationary-action principle, using path integrals. Maxwell's equations can be derived as conditions of stationary action.



See also



  • Analytical mechanics

  • Configuration space

  • Hamilton–Jacobi equation

  • Phase space

  • Geodesics as Hamiltonian flows



References





  1. ^ R. Penrose (2007). The Road to Reality. Vintage books. p. 474. ISBN 0-679-77631-1..mw-parser-output cite.citation{font-style:inherit}.mw-parser-output .citation q{quotes:"""""""'""'"}.mw-parser-output .citation .cs1-lock-free a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/6/65/Lock-green.svg/9px-Lock-green.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .citation .cs1-lock-limited a,.mw-parser-output .citation .cs1-lock-registration a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/d/d6/Lock-gray-alt-2.svg/9px-Lock-gray-alt-2.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .citation .cs1-lock-subscription a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/a/aa/Lock-red-alt-2.svg/9px-Lock-red-alt-2.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .cs1-subscription,.mw-parser-output .cs1-registration{color:#555}.mw-parser-output .cs1-subscription span,.mw-parser-output .cs1-registration span{border-bottom:1px dotted;cursor:help}.mw-parser-output .cs1-ws-icon a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/4/4c/Wikisource-logo.svg/12px-Wikisource-logo.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output code.cs1-code{color:inherit;background:inherit;border:inherit;padding:inherit}.mw-parser-output .cs1-hidden-error{display:none;font-size:100%}.mw-parser-output .cs1-visible-error{font-size:100%}.mw-parser-output .cs1-maint{display:none;color:#33aa33;margin-left:0.3em}.mw-parser-output .cs1-subscription,.mw-parser-output .cs1-registration,.mw-parser-output .cs1-format{font-size:95%}.mw-parser-output .cs1-kern-left,.mw-parser-output .cs1-kern-wl-left{padding-left:0.2em}.mw-parser-output .cs1-kern-right,.mw-parser-output .cs1-kern-wl-right{padding-right:0.2em}




  • W.R. Hamilton, "On a General Method in Dynamics.", Philosophical Transactions of the Royal Society Part II (1834) pp. 247–308; Part I (1835) pp. 95–144. (From the collection Sir William Rowan Hamilton (1805–1865): Mathematical Papers edited by David R. Wilkins, School of Mathematics, Trinity College, Dublin 2, Ireland. (2000); also reviewed as On a General Method in Dynamics)

  • Goldstein H. (1980) Classical Mechanics, 2nd ed., Addison Wesley, pp. 35–69.

  • Landau LD and Lifshitz EM (1976) Mechanics, 3rd. ed., Pergamon Press.
    ISBN 0-08-021022-8 (hardcover) and
    ISBN 0-08-029141-4 (softcover), pp. 2–4.

  • Arnold VI. (1989) Mathematical Methods of Classical Mechanics, 2nd ed., Springer Verlag, pp. 59–61.

  • Cassel, Kevin W.: Variational Methods with Applications in Science and Engineering, Cambridge University Press, 2013.




這個網誌中的熱門文章

Xamarin.form Move up view when keyboard appear

Post-Redirect-Get with Spring WebFlux and Thymeleaf

Anylogic : not able to use stopDelay()