# Theoretical Mechanics - PHY6200 - PowerPoint PPT Presentation

Theoretical Mechanics - PHY6200

1 / 37
Theoretical Mechanics - PHY6200

## Theoretical Mechanics - PHY6200

- - - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - - -
##### Presentation Transcript

1. Theoretical Mechanics - PHY6200 Chapter 6 Introduction to the calculus of variations Prof. Claude A Pruneau, Physics and Astronomy Department Wayne State University

2. Introduction • Many mechanics problem more easily analyzed/solved by means of calculus of variation, Lagrange Eqs., etc. • We will: • Consider some general principles • Omit formal existence proofs • As an example, consider Fermat’s principle.

3. air water Fermat’s Principle • Light travels between two points along the path that takes the least amount of time. L

4. Snell’s Law By construction (geometry): Total travel time: Extremum for : Simplification/Substitution : Snell’s Law :

5. 2 1 Statement of the problem • In classical mechanics, the problem amounts to finding a trajectory between an initial and a final point (boundary conditions). • Express the problem to solve as a minimization problem. • Find a functional that once minimized will yield an optimal trajectory which corresponds to the physical trajectory (according to Newton’s laws of motion).

6. 2 1 Basic Formulation (1 dimension) • Determine a function, y(x), such that the following integral is an extremum. • f(y,y’;x) is a function considered as “given” • Limits of integration are fixed. • y(x), the trajectory is to be varied until an extremum is found for J. • y’(x)=dy(x)/dx • x : independent variable

7. Basic Formulation (cont’d) • If the optimal trajectory y(x) produces a minimal value for J, then any neighboring trajectory ya(x) corresponds to a large value of J. • Hence by minimizing J, one finds the optimal trajectory. 2 1

8. Parametric Representation • The function y(x) can be considered with a parametric representation • ya(x) = y(a,x) • Where the variable a denotes a parameter that describes functions close to y(x) but that differ by an “amount” proportional to a. • By definition, a=0 corresponds to the optimal trajectory y(x). • y0(x) = y(0,x)=y(x)

9. Parametric Representation (cont’d) • Consider the specific representation • y(a,x)=y(0,x)+ (x) • where  (x) is some “arbitrary” function of x, with continuous 1st derivative, and vanishing at the boundaries • (x1) =  (x2) = 0

10. Parametric Representation (cont’d) • The functional J is written • This integral has an extremum if • This must be true for all functions h(x) • Note: This is a necessary but not sufficient condition.

11. Example 1 • Consider the function f(y’;x)=(dy/dx)2. • Where y(x)=x. • Consider h(x)=sin(x) • Find J() between x=0, and x=2. • Show that the stationary value (extremum) of J occurs for  =0.

12. Example 1 - Solution • Let: • y(a,x) = x +  sin(x) • Note that by choice/construction, we have • (0) = 0 • (2) = 0 • Proceed to calculate…

13. Example 1 - Solution (cont’d)

14. Derivation of Euler’s Equation • Calculate the derivative of J • For fixed integration limits, one can change the order of operations

15. By construction: • Thus:

16. 0 • Integrate 2nd term by parts • The function h(x) being totally arbitrary. • We get • Provided the integrand itself vanishes. • So…

17. Euler’s Equation This equation yields a solution for y(x) and y’(x) which produces an extremum for J.

18. Example 2 - Brachistochrone • Consider a particle moving in a constant force field, starting at rest from some point (x1,y1) and ending at (x2,y2). • Find the path (trajectory) that allows the particle to accomplish the motion/transit in the least amount of time.

19. Example 2 - Solution (1) • Constant force, neglect friction, implies a conservative system. • T+U=constant. • Let U(x=0)=0, and T(0)+U(0)=0. • T=0.5 mv2 • U=-Fx=-mgx y (x1,y1) (x2,y2) x

20. Example 2 - Solution (2) • Time required for transit:

21. Example 2 - Solution (3) • Use Euler’s Equation • Thus…

22. y (0,2pa) (0,pa) (0,0) q (2a,0) cycloid x Example 2 - Brachistochrone - Answer Solution passes through (0,0)

23. 2nd Form of Euler’s Equation • For a function “f” that does not explicitly depend on “x”. 1) 2)

24. 2nd Form of Euler’s Equation (cont’d) • We find: 2nd Form of Euler’s Equation is thus If f does not depend on “x”, there is a conserved quantity

25. Function with Several Dependent Variables • Consider a more general case where f is a functional of several dependent variables. • E.g. Motion of one particle in 3 dimensions, or many particles…

26. Several Dependent Variables… • Repeat our previous reasoning… yi(a,x)=yi (0,x)+ ahi(x)

27. Several Dependent Variables… • So we have • Functions hi are all independent. • Each term of the sum must be null for the above derivative to vanish at a=0. Euler’s Equations for many dependent variables

28. Euler Equations with Auxiliary Conditions • In many problems, additional constrains come into play. • Consider e.g. motion constrained to a spherical shell, or some other type of curved surface. • Then obviously the path (motion) must be on the surface, and thereby satisfy the equation of the surface, which can be generally be written: g{yi;x}=0 • One therefore introduce constraint equations.

29. and are thus no longer independent. The variations Euler Equations with Auxiliary Conditions Consider a case: Also include constraints:

30. No terms in x appear since We write: The constraints equation becomes: or:

31. The left hand side involves only derivatives of f and g with respect to y and y’, whereas the right hand side involves only derivatives of f and g with respect to z and z’. Because y and z are both functions of x, the two sides may be set equal to some function of x, which we note: -l(x). The complete solution requires one finds three functions: y(x), z(x), and l(x). Fortunately, three relations may be used: the above two equations and the constraints Eq.

32. The function l(x) is known as Lagrange underdetermined multiplier. In general, if there are many external constraints, one has: m equations m+n unknowns n equations Note that the constraints can also be written: Often more “useful” in problem solutions….

33. The d Notation • To simplify writing in calculus of variations, one uses a shorthand notation. Consider:

34. The d Notation • Condition of extremum • Swap d and integral: • Note: • So…

35. The  Notation • Do keep in mind that the d notation is only a shorthand expression for a more detailed and precise quantity. • You can visualize the varied path y as a virtual displacement from the actual path consistent with all forces and constraints. It is to be distinguished from dy (differential displacement) by the condition that dt=0 (fixed time)

36. Geodesic on a sphere