Chapter 2 Lagrangian Modeling
|
|
- Jordan Cannon
- 5 years ago
- Views:
Transcription
1 Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie key concepts to moel-relate physical phenomenon. The Lagrangian formulation of moeling erives from the basic work energy principle an Newton s laws of motion. The basic law states that the force acting on a boy is irectly proportional to its acceleration equate with constant mass. F ¼ ma ¼ m v ¼ m 2 x 2 ð2:1þ In this equation, force is irectly proportional to the erivative of velocity v or ouble erivative of isplacement x with respect to time. If a boy moves from point A to point B then the work one by a boy is a ot prouct of force an isplace path integrate from point A to B. W ¼ ð B A W ¼ ð B A F x m 2 x x ¼ 2 ð B A m x x ð2:2þ ð2:3þ We now erive the equation as xx ¼ _x x x ¼ _x ¼ _x_x ð2:4þ Springer International Publishing Switzerlan 2016 A.M. Mughal, Real Time Moeling, Simulation an Control of Dynamical Systems, DOI / _2 19
2 20 2 Lagrangian Moeling So Eq. (2.3) formulates as ð B W ¼ m _x_x ¼ m _x 2 B ¼ m A 2 A 2 _x 2 B _x 2 A W ¼ K B K A ¼ ΔK ð2:5þ ð2:6þ It is important to know that here we are iscussing motion in a conservative fiel where work one is inepenent of path an epens upon the ifference between the initial an final value of kinetic energy K. We note that work one is represente as change in K between two points. By the law of conversation of energy an work energy principle we know that a change in kinetic energy shoul also change the potential energy U, as the total energy of the system remains constant. So a change in kinetic an change in potential energies shoul have a net effect of zero as ΔK þ ΔU ¼ 0 ð2:7þ So ΔU ¼W ¼ ð B A F x ð2:8þ Using the funamental theorem of calculus for antierivatives, we express Eq. (2.8) as F ¼ U x ð2:9þ This shows forces require to change the potential energy, whereas the force require to change the kinetic energy is etermine from the efinition of K. K ¼ 1 2 mv2 ¼ 1 2 m _x 2 ð2:10þ So K _x ¼ m _x ð2:11þ Differentiating Eq. (2.11) gives a force of kinetic energy F ¼ K _x ¼ ð m _x Þ ¼ ma ð2:12þ
3 2.1 Moeling in Three Axes 21 As sum of forces equal to zero so K _x U x ¼ 0 ð2:13þ Now we efine a Lagrangian function L as the ifference of kinetic an potential energies L ¼ K U ð2:14þ We observe that kinetic energy K is a function of velocity _x, an potential energy U is a function of isplacement x. Accoringly the Lagrangian function L has two terms as a function of velocity x _ : an isplacement x inepenently. As we know that ¼ K an _x _x x ¼ U ð2:15þ x Lagrangian function L now represents Eq. (2.13) as _x x ¼ 0 ð2:16þ This is Lagrangian formulation for a system with motion only in single axis. 2.1 Moeling in Three Axes Now consier a motion in three axes an every equation now to be represente in three imensional vector notations. F * ¼ m 2 x * 2 ð2:17þ The force F * has three components, F x, F y, an F z in x, y an z axes respectively. The potential energy U is a scalar fiel in three imensions an the negative graient of this fiel provies a force vector in three axes. Now Eq. (2.9) represents as F * ¼ Ux; ð y; zþ ¼ U T U U ð2:18þ x y z
4 22 2 Lagrangian Moeling The kinetic energy K in three axes is given as K ¼ 1 2 m _x* 2 ¼ 1 T 2 m _x* * _x ¼ 1 2 m _x 2 þ _y 2 þ _z 2 ð2:19þ This provies the three components of a force through kinetic energy Eq. (2.19) as F * ¼ T K K K ð2:20þ _x _y _z By equating forces we can write Eq. (2.13) for three axes as T K K K U T U U ¼ 0 ð2:21þ _x _y _z x y z This gives us Lagrangian formulation of rigi boy (particle) for each axis in Cartesian coorinates as _x _y _z x ¼ 0 y ¼ 0 z ¼ 0 ð2:22þ 2.2 Moeling in Cartesian an Spherical/Polar Coorinates In the above section, we evelope a moeling scheme in only a Cartesian coorinate system but often the moel is represente in spherical coorinates in 3D space or with polar coorinates in 2D space. The translation an rotational motion also become part of a moel in ynamical equations. Now, we iscuss the metho in 3D space with spherical coorinates, which can also be use in polar coorinates by eliminating an axis. We know that x ¼ fðr; θ; φ; tþ y ¼ fðr; θ; φ; tþ z ¼ fðr; θ; φ; tþ ð2:23þ We efine our system in Cartesian coorinates by using the transformation in Eq. (2.23). Now let us consier that we want to generalize the concept of spherical,
5 2.2 Moeling in Cartesian an Spherical/Polar Coorinates 23 polar, or cylinrical coorinates transforme into Cartesian coorinates by using a generalize set of variables instea of r, θ, φ. We now efine a generalize vector * q ¼ ½ q1 q 2 q 3 T an the position in Cartesian coorinates with a vector * p ¼ ½x y z T ¼ ½p 1 p 2 p 3 T. 2 3 p 1 ¼ fðq 1 ; q 2 ; q 3 ; tþ p 2 ¼ fðq 1 ; q 2 ; q 3 ; tþ5 ¼ * p ¼ f * q ; t ð2:24þ p 3 ¼ fðq 1 ; q 2 ; q 3 ; tþ The functions in Eq. (2.24) follow the chain rule for respective erivatives because x, y, z are ifferentiable with respect to time as well as with respect to each generalize coorinates. Now we follow the convention that _x ¼ x ; a ot on top of the variables represents its ifferentiation with respect to time for both Cartesian an generalize coorinates. The change with respect to time in Cartesian coorinates is efine as _p 1 ¼ p 1 q 1 _q 1 þ p 1 q 2 _q 2 þ p 1 q 3 _q 3 _p 2 ¼ p 2 q 1 _q 1 þ p 2 q 2 _q 2 þ p 2 q 3 _q 3 ð2:25þ _p 3 ¼ p 3 q 1 _q 1 þ p 3 q 2 _q 2 þ p 3 q 3 _q 3 Now consier a generalize notation for any element p i (i.e. x, y, z) of a position vector * p with respect to time erivates an generalize coorinates _p i ¼ X3 j¼1 _q j ð2:26þ The change with respect to generalize coorinates is given as δp i ¼ X3 j¼1 δq j ð2:27þ We observe that _p i ¼ p i þ _p i ð2:28þ
6 24 2 Lagrangian Moeling Solving for the secon erivative of position coorinate, we get p i ¼ _p i _p i ð2:29þ Differentiating Eq. (2.26) with respect to q j where j ¼ f1; 2; 3g _p i _q j ¼ ð2:30þ As p i is a function of q j then the partial erivatives with respect to time are given as ¼ X3 k¼1 q k _q q k j ð2:31þ We know that in complex or real fiels, secon-orer cross partial erivatives are equal for any position variable with respect to other generalize coorinates q k ¼ q k ð2:32þ Now taking partial erivative of Eq. (2.26) with respect to q j we get _p i ¼ X3 k¼1 _q q k ¼ X3 k By comparing Eqs. (2.31) an (2.33) we get k¼1 ¼ _p i q k Now substituting Eqs. (2.30) an (2.34) in Eq. (2.29) p i ¼ _p i _p i _q j _p i _p i _q q k j ð2:33þ ð2:34þ ð2:35þ Now we write Eq. (2.35) as p i ¼ 2 _p i 2 _p i _q j 2 2 ð2:36þ
7 2.3 Work an Energy Formulation Work an Energy Formulation We know that change in work one is the ot prouct of force applie an corresponing change in position given as δw ¼ F * δp * ð2:37þ Equivalently, force vector F * has three components F P1, F P2 ; an F P3, such as δw ¼ X3 F Pi δ Pi ð2:38þ Each component of force is expresse with Newton s law as Equating δw in Eq. (2.38) we have F Pi ¼ m p i ð2:39þ δw ¼ X3 m p i δ Pi ¼ X3 F Pi δ Pi ð2:40þ Substituting Eq. (2.27) in Eq. (2.40) we get δw ¼ X3 j¼1 X 3 m p i Pi δ qj ¼ X3 j¼1 X 3 F Pi Pi qj δ qj ð2:41þ Let F qj ¼ X3 F Pi Pi qj ð2:42þ Equation (2.41) becomes δw ¼ X3 j¼1 X 3 m p i Pi δ qj ¼ X3 j¼1 F qj δ qj ð2:43þ
8 26 2 Lagrangian Moeling Now by substituting Eq. (2.36) in Eq. (2.43) we get δw ¼ X3 j¼1 m X3 2 _p i 2 _p i δ qj ¼ X3 F qj δ qj _q j 2 2 j¼1 ð2:44þ By efinition of kinetic energy of a particle of mass m is given as K ¼ m 2 X 3 _p 2 i ð2:45þ There are three Lagrange equations which generate from Eq. (2.44) for each axis ðj ¼ 1, 2, 3Þ as m X3 2 _p i 2 _p i ¼ F qj ð2:46þ _q j 2 2 Each value of j represents a Lagrange equation in its axis. The purpose of keeping δw in Eq. (2.44) is that we nee to equate change in energy through work one by other physical interpretations. The Lagrange equation eals with the kinetic energy of an object in motion in Eq. (2.44). Work one is also represente by change in potential energy as given in Eq. (2.8). There may be non-conservative forces that are acting on the boy or internal energy of the system, but ue to the law of conservation of energy these forces must balance each other. The negative graient of potential energy provies forces acting upon the boy which may be causing it to move or stop. So if we efine potential energy U in Cartesian coorinates, forces in Cartesian coorinates are represente as F Pi ¼ U ð2:47þ Equation (2.42) represent as F qj ¼ X3 U Pi qj ¼ U Now we write Eq. (2.46) as a Lagrange equation in q j coorinates m X3 2 _p i 2 _p i ¼ U _q j 2 2 ð2:48þ ð2:49þ Again using the efinition of kinetic energy from Eq. (2.45) to represent Eq. (2.49) in energy variables
9 2.3 Work an Energy Formulation 27 K K ¼ U _q j ð2:50þ Kinetic energy is a function of both position an velocity of a particle in generalize coorinate q j an the potential energy is only a function of generalize position coorinates, ¼ _q j K _q j an ¼ U ð2:51þ Now using a Lagrange variable L from Eq. (2.14) in three axes system an using Eq. (2.51) in each coorinates we get the following equation for a conservative system in each q j coorinate _q j ¼ 0 ð2:52þ In a conservative fiel Eq. (2.7) hols for change in both potential an kinetic energy, but in a non-conservative fiel the change is also ue to non-conservative, external, or internal forces in each coorinate system. E nc is a non-conservative energy ue to all of these forces acting upon the boy, so the energy equation becomes ΔK þ ΔU ¼ E nc ð2:53þ In this case the Lagrange equation (2.52) oes not satisfy the ue non-conservative force acting upon a boy. Let F ncqj is the sum of all non-conservative forces acting in q j irection. The generalize Lagrange equation in generalize coorinates is _q j ¼ F ncqj ð2:54þ Depening upon the space in which coorinates are efine, the value of j can be change. For a motion with two egrees of freeom there will be only two Lagrange equations but in a motion with three egrees of freeom there will be three Lagrange equations to escribe the system. In certain cases, there is motion in both Cartesian coorinates an generalize coorinates; then accoringly a irection of motion in Cartesian is also treate as a variable q j.
10 28 2 Lagrangian Moeling 2.4 State Space Representation The Lagrange equation itself is a nonlinear or linear state space representation. All non-conservative forces are treate as input to the system whether these are controllable or exogenous. The generalize coorinates an their first-orer erivatives constitute the state vector. Example 2.1: A Simple Penulum The simple penulum is a boy with mass m attache at length l of massless string from the origin (or groun) to move between two en points as shown in Fig If no external force is applie then the movement will epen upon the starting point an eventually comes to en at mass exactly below the origin. In Cartesian coorinates x an y mass move in a semicircle, which can also be relate through polar coorinates l an θ, which are fixe length of a string an angle from the origin respectively. There is only one egree of freeom, so there is only one generalize coorinate q 1 ¼ θ. Conventionally these simple examples are represente with θ instea of q 1. Now we efine the position an velocity of Cartesian coorinate as a function of a generalize coorinate as xðθþ ¼ l sin θ yðθþ ¼ l cos θ _x θ; θ _ ¼ l θ _ cos θ _y θ; _θ ¼ l _θ sin θ ð2:55þ ð2:56þ The kinetic energy of the mass m in Cartesian coorinate is given as K ¼ 1 2 mv2 ¼ 1 2 m _x 2 þ _y 2 1 ¼ 2 ml2 _θ 2 ð2:57þ The potential energy U is ue to gravitational pull an it is lowest when θ ¼ 0 where kinetic energy is highest an vice versa when θ ¼ π 2 or π 2.So or U ¼ mglð1 cos θþ ð2:58þ Fig. 2.1 A simple penulum
11 2.4 State Space Representation 29 The Lagrangian L ¼ K U is given as L ¼ 1 2 ml2 θ _ 2 mglð1 cos θþ ð2:59þ The Lagrange equation of the system erives as follows _θ θ ¼ 0 ð2:60þ ml 2 _θ mgl sin θ ¼ 0 ð2:61aþ ml 2 θ þ mgl sin θ ¼ 0 ð2:61bþ The final equation with one egree of freeom now appears as θ þ g l sin θ ¼ 0 ð2:62þ State space formulation of the system represent Eq. (2.62) with efinition of state variables given as x 1 ¼ θ an x 2 ¼ _θ. We have state vector * x ¼ θ _θ an the nonlinear state space representation as _x 1 ¼ x 2 _x 2 ¼ g l sin x 1 ð2:63þ The system can be linearize at relaxe equilibrium point " # θ e T _θ e ¼ ½0 0 with only one matrix A ¼ 0 1 g 0 l Example 2.2: Inverte Penulum on a Moving Cart A cart of mass M is moving by an applie force F acting along x-axis. An inverte penulum of mass m is attache on the center of cart with a ro of length l (Fig. 2.2). The motion of the cart is efine from origin an the angle is measure from vertical y-axis. The system can either be represente with horizontal an vertical variables of cart an bob or by using horizontal isplacement of cart an angular position of bob. There are only two egrees of freeom; an angle can be measure from the horizontal but in orer to be conversant with remaining literature we chose the angle θ measure from y-axis. So for 2-DoF system, there are two generalize coorinates x 1 ¼ x, the cart s isplacement an θ the angle of a bob from vertical. The position of bob (x 2, y 2 ) can be given in terms of x 1 an θ as
12 30 2 Lagrangian Moeling Fig. 2.2 Inverte penulum on a moving cart 2 m 2 1 ( ) M x 2 ¼ x þ l sin θ y 2 ¼ l cos θ ð2:64þ The velocity components of bob are given as _x 2 ¼ _x þ l _θ cos θ _y 2 ¼l _θ sin θ ð2:65þ The kinetic energies of both masses are given as K 1 ¼ 1 2 M _x 2 K 2 ¼ 1 2 m _x þ _y 2 K 2 ¼ 1 h 2 m _x þ l 2 2 i _θ cos θ þ l _θ sin θ K 2 ¼ 1 2 m _x 2 þ 2l _x θ _ cos θ þ l 2 θ _ 2 ð2:66þ ð2:67aþ ð2:67bþ ð2:67cþ So the total kinetic energy of the system from Eqs. (2.66) an (2.67c) is given as K ¼ K 1 þ K 2 ¼ 1 2 ðm þ mþ_x 2 þ ml _x _θ cos θ þ 1 2 ml2 _θ 2 ð2:68þ The potential energy in a bob is given as U ¼ mgy 2 ¼ mgl cos θ ð2:69þ
13 2.4 State Space Representation 31 The Lagrangian L is given as L ¼ K U ¼ 1 2 ðm þ mþ_x 2 þ ml _x θ _ cos θ þ 1 2 ml2 θ _ 2 mgl cos θ ð2:70þ we represent two equations of motion sing the variable L an two generalize coorinates (x, θ) as _x x ¼ F ð2:71þ an _θ θ ¼ 0 ð2:72þ From Eqs. (2.70), (2.71), an (2.72) are solve as So equations of motion are _x ¼ þ mþ_x θ _ cos θ ð2:73aþ _θ ¼ ml _x cos θ þ ml2 _θ ð2:73bþ x ¼ 0 ¼ mgl sin θ θ ð2:74þ ð2:73cþ ðm þ mþ x þ ml θ cos θ ml _θ 2 sin θ ¼ F ð2:75þ ml x cos θ ml _x _θ sin θ þ ml 2 θ mgl sin θ ¼ 0 ð2:76þ These are nonlinear equations for the system of a moving cart with an inverte penulum. In orer to make a state space system, first we nee to solve variables with ouble erivatives simultaneously. Each generalize coorinate of the system is secon-orer so a total orer of the system is 4. We have both x an θ, so we nee to solve these equations simultaneously. Eqs. (2.75) an (2.76) are represente as ml 2 l θ ¼ F ðm þ mþ x þ ml θ cos θ _ 2 sin θ ð2:77þ x ¼ _x _θ tan θ l θ cos θ þ g tan θ ð2:78þ
14 32 2 Lagrangian Moeling Inserting value of x from Eq. (2.78) ineq.(2.75) an then inserting the value of ml 2 θ from Eq. (2.77) in Eq. (2.76) simultaneously solve the equations for state space T representation of a state vector x θ _x _θ as x ¼ m cos θ M þ m 1 m _x _θ sin θ þ mg sin θ þ ml _θ 2 tan θ þ F cos θ cos θ θ ¼ 1 l M cos θ Problems m cos θ þ m cos θ ð2:79þ 1 ðm þ mþ _x θ _ tan θ þ g tan θ þ ml θ _ 2 sin θ þ F P2.1 Consier spherical coorinates of mass m in motion x ¼ r cos θ sin φ y ¼ r sin θ sin φ z ¼ r cos φ ð2:80þ Formulate Lagrangian equation of motion in three axis for three generalize coorinates an represent as state space moel. P2.2 A penulum given in Example 2.1 is rotating with angular velocity ω in a circular path of raius r. The motion is uner the influence of gravity with no external forces acting upon it. The Cartesian coorinates are given as xt ðþ¼rcos ðωtþþl cos θ yt ðþ¼rsin ðωtþþ l sin θ Fin the state space representation of a system through Lagrangian formulation. P2.3 Define a generalize coorinate system for the problem given in Example 2.2 an represent its nonlinear state space formulation. Use the symbolic math toolbox to obtain Jacobian matrices A an B given in Eq. (1.19) an Eq. (1.20) by linearizing at relaxe position. P2.4 Nonlinear systems are also linearize by approximating nonlinear terms within a specifie range on equilibrium point. The inverte penulum on a cart can be linearize by approximating a small angle an negligible velocities at equilibrium point for trigonometric functions i.e., sin θ θ, cos θ 1. Fin a linear state space equation by these assumptions an obtain a state space representation of the inverte penulum on a cart system.
15 References 33 P2.5 The position of bob of the inverte penulum on a moving cart can be controlle by applying external torque τ at the hinge changing the Eq. (2.78) as ml x cos θ ml _x _θ sin θ þ ml 2 θ mgl sin θ ¼ τ Linearize the complete system by the assumption given in P2.4 an obtain state space representation to monitor all state variables at output inepenently. References 1. Frielan, Bernar Control System Design An Introuction to State Space Methos. Mineola: Dover Publications. 2. Zak, Stanislaw H Systems an Control. Oxfor: Oxfor University Press. 3. Brizar, Alain J An Introuction to Lagrangian Mechanics. Singapore: Worl Scientific Publishing. 4. Wellstea, P.E Introuction to Physical System Moelling. Lonon: Acaemic Press. 5. Amirouche, Fari Funementals of Multiboy Dynamics Theory an Applications. Basel: Birkhäuser.
16
The Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationG j dq i + G j. q i. = a jt. and
Lagrange Multipliers Wenesay, 8 September 011 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationand from it produce the action integral whose variation we set to zero:
Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationShort Intro to Coordinate Transformation
Short Intro to Coorinate Transformation 1 A Vector A vector can basically be seen as an arrow in space pointing in a specific irection with a specific length. The following problem arises: How o we represent
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationVectors in two dimensions
Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication
More informationGyroscopic matrices of the right beams and the discs
Titre : Matrice gyroscopique es poutres roites et es i[...] Date : 15/07/2014 Page : 1/16 Gyroscopic matrices of the right beams an the iscs Summary: This ocument presents the formulation of the matrices
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationLecture 27: Generalized Coordinates and Lagrange s Equations of Motion
Lecture 27: Generalize Coorinates an Lagrange s Equations of Motion Calculating T an V in terms of generalize coorinates. Example: Penulum attache to a movable support 6 Cartesian Coorinates: (X, Y, Z)
More informationChapter 2 Governing Equations
Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement
More informationTutorial Test 5 2D welding robot
Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.
More informationProblem set 2: Solutions Math 207B, Winter 2016
Problem set : Solutions Math 07B, Winter 016 1. A particle of mass m with position x(t) at time t has potential energy V ( x) an kinetic energy T = 1 m x t. The action of the particle over times t t 1
More informationGoldstein Chapter 1 Exercises
Golstein Chapter 1 Exercises Michael Goo July 17, 2004 1 Exercises 11. Consier a uniform thin isk that rolls without slipping on a horizontal plane. A horizontal force is applie to the center of the isk
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationMath Chapter 2 Essentials of Calculus by James Stewart Prepared by Jason Gaddis
Math 231 - Chapter 2 Essentials of Calculus by James Stewart Prepare by Jason Gais Chapter 2 - Derivatives 21 - Derivatives an Rates of Change Definition A tangent to a curve is a line that intersects
More informationAN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A
AN INTRODUCTION TO AIRCRAFT WIN FLUTTER Revision A By Tom Irvine Email: tomirvine@aol.com January 8, 000 Introuction Certain aircraft wings have experience violent oscillations uring high spee flight.
More informationNumerical Integrator. Graphics
1 Introuction CS229 Dynamics Hanout The question of the week is how owe write a ynamic simulator for particles, rigi boies, or an articulate character such as a human figure?" In their SIGGRPH course notes,
More informationFree rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012
Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate
More informationTopic 2.3: The Geometry of Derivatives of Vector Functions
BSU Math 275 Notes Topic 2.3: The Geometry of Derivatives of Vector Functions Textbook Sections: 13.2 From the Toolbox (what you nee from previous classes): Be able to compute erivatives scalar-value functions
More information1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity
AP Physics Multiple Choice Practice Electrostatics 1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity. A soli conucting sphere is given a positive charge Q.
More informationLinear and quadratic approximation
Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function
More information6. Friction and viscosity in gasses
IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner
More information12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes
Differentiation of vectors 12.5 Introuction The area known as vector calculus is use to moel mathematically a vast range of engineering phenomena incluing electrostatics, electromagnetic fiels, air flow
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationStatic Equilibrium. Theory: The conditions for the mechanical equilibrium of a rigid body are (a) (b)
LPC Physics A 00 Las Positas College, Physics Department Staff Purpose: To etermine that, for a boy in equilibrium, the following are true: The sum of the torques about any point is zero The sum of forces
More informationu t v t v t c a u t b a v t u t v t b a
Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationSummary: Differentiation
Techniques of Differentiation. Inverse Trigonometric functions The basic formulas (available in MF5 are: Summary: Differentiation ( sin ( cos The basic formula can be generalize as follows: Note: ( sin
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationLecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.
b. Partial erivatives Lecture b Differential operators an orthogonal coorinates Recall from our calculus courses that the erivative of a function can be efine as f ()=lim 0 or using the central ifference
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationAngles-Only Orbit Determination Copyright 2006 Michel Santos Page 1
Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More informationPartial Differential Equations
Chapter Partial Differential Equations. Introuction Have solve orinary ifferential equations, i.e. ones where there is one inepenent an one epenent variable. Only orinary ifferentiation is therefore involve.
More informationInterpolated Rigid-Body Motions and Robotics
Interpolate Rigi-Boy Motions an Robotics J.M. Selig Faculty of Business, Computing an Info. Management. Lonon South Bank University, Lonon SE AA, U.K. seligjm@lsbu.ac.uk Yaunquing Wu Dept. Mechanical Engineering.
More informationThe Ehrenfest Theorems
The Ehrenfest Theorems Robert Gilmore Classical Preliminaries A classical system with n egrees of freeom is escribe by n secon orer orinary ifferential equations on the configuration space (n inepenent
More informationMath 115 Section 018 Course Note
Course Note 1 General Functions Definition 1.1. A function is a rule that takes certain numbers as inputs an assigns to each a efinite output number. The set of all input numbers is calle the omain of
More informationLagrangian and Hamiltonian Dynamics
Lagrangian an Hamiltonian Dynamics Volker Perlick (Lancaster University) Lecture 1 The Passage from Newtonian to Lagrangian Dynamics (Cockcroft Institute, 22 February 2010) Subjects covere Lecture 2: Discussion
More informationCalculus of Variations
Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,
More informationG4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const.
G4003 Avance Mechanics 1 The Noether theorem We alreay saw that if q is a cyclic variable, the associate conjugate momentum is conserve, q = 0 p q = const. (1) This is the simplest incarnation of Noether
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationIMPLICIT DIFFERENTIATION
IMPLICIT DIFFERENTIATION CALCULUS 3 INU0115/515 (MATHS 2) Dr Arian Jannetta MIMA CMath FRAS Implicit Differentiation 1/ 11 Arian Jannetta Explicit an implicit functions Explicit functions An explicit function
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationII. First variation of functionals
II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent
More informationThe Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis
The Kepler Problem For the Newtonian 1/r force law, a miracle occurs all of the solutions are perioic instea of just quasiperioic. To put it another way, the two-imensional tori are further ecompose into
More informationACCELERATION, FORCE, MOMENTUM, ENERGY : solutions to higher level questions
ACCELERATION, FORCE, MOMENTUM, ENERGY : solutions to higher level questions 015 Question 1 (a) (i) State Newton s secon law of motion. Force is proportional to rate of change of momentum (ii) What is the
More informationMath 1271 Solutions for Fall 2005 Final Exam
Math 7 Solutions for Fall 5 Final Eam ) Since the equation + y = e y cannot be rearrange algebraically in orer to write y as an eplicit function of, we must instea ifferentiate this relation implicitly
More informationSection 7.1: Integration by Parts
Section 7.1: Integration by Parts 1. Introuction to Integration Techniques Unlike ifferentiation where there are a large number of rules which allow you (in principle) to ifferentiate any function, the
More information4. Important theorems in quantum mechanics
TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this
More informationTHE DISPLACEMENT GRADIENT AND THE LAGRANGIAN STRAIN TENSOR Revision B
HE DISPLACEMEN GRADIEN AND HE LAGRANGIAN SRAIN ENSOR Revision B By om Irvine Email: tom@irvinemail.org Febrary, 05 Displacement Graient Sppose a boy having a particlar configration at some reference time
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationTHERE ARE BASICALLY TWO approaches to conservative mechanics, due principally
CHAPTER 2 LAGRANGE S EQUATIONS THERE ARE BASICALLY TWO approaches to conservative mechanics, ue principally to Lagrange an Hamilton respectively Each has avantages over the other, although the Lagrangian
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More information1 Heisenberg Representation
1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.
More informationIntroduction to variational calculus: Lecture notes 1
October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational
More informationThe Principle of Least Action and Designing Fiber Optics
University of Southampton Department of Physics & Astronomy Year 2 Theory Labs The Principle of Least Action an Designing Fiber Optics 1 Purpose of this Moule We will be intereste in esigning fiber optic
More informationOptimization of Geometries by Energy Minimization
Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.
More informationDOUBLE PENDULUM VIBRATION MOTION IN FLUID FLOW
ENGINEERING FOR RURA DEEOPMENT Jelgava,.-.5.5. DOUBE PENDUUM IBRATION MOTION IN FUID FOW Janis iba, Maris Eiuks, Martins Irbe Riga Technical University, atvia janis.viba@rtu.lv, maris.eiuks@rtu.lv, martins.irbe@rtu.lv
More informationIntegration Review. May 11, 2013
Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationExam 2 Review Solutions
Exam Review Solutions 1. True or False, an explain: (a) There exists a function f with continuous secon partial erivatives such that f x (x, y) = x + y f y = x y False. If the function has continuous secon
More informationChapter Primer on Differentiation
Capter 0.01 Primer on Differentiation After reaing tis capter, you soul be able to: 1. unerstan te basics of ifferentiation,. relate te slopes of te secant line an tangent line to te erivative of a function,.
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationFurther Differentiation and Applications
Avance Higher Notes (Unit ) Prerequisites: Inverse function property; prouct, quotient an chain rules; inflexion points. Maths Applications: Concavity; ifferentiability. Real-Worl Applications: Particle
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationProblem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs
Problem Sheet 2: Eigenvalues an eigenvectors an their use in solving linear ODEs If you fin any typos/errors in this problem sheet please email jk28@icacuk The material in this problem sheet is not examinable
More informationAverage value of position for the anharmonic oscillator: Classical versus quantum results
verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive
More informationClass Notes for Advanced Dynamics (MEAM535)
Class Notes f Avance Dynamics MEAM535 Michael A. Carchii December 9, 9 Equilibrium Points, Small Perturbations Linear Stability Analysis The following notes were initially base aroun the text entitle:
More informationTotal Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables*
51st IEEE Conference on Decision an Control December 1-13 212. Maui Hawaii USA Total Energy Shaping of a Class of Uneractuate Port-Hamiltonian Systems using a New Set of Close-Loop Potential Shape Variables*
More informationOptimization of a point-mass walking model using direct collocation and sequential quadratic programming
Optimization of a point-mass walking moel using irect collocation an sequential quaratic programming Chris Dembia June 5, 5 Telescoping actuator y Stance leg Point-mass boy m (x,y) Swing leg x Leg uring
More informationGravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann
Gravitation as the result of the reintegration of migrate electrons an positrons to their atomic nuclei. Osvalo Domann oomann@yahoo.com (This paper is an extract of [6] liste in section Bibliography.)
More informationDifferentiation Rules Derivatives of Polynomials and Exponential Functions
Derivatives of Polynomials an Exponential Functions Differentiation Rules Derivatives of Polynomials an Exponential Functions Let s start with the simplest of all functions, the constant function f(x)
More informationDot trajectories in the superposition of random screens: analysis and synthesis
1472 J. Opt. Soc. Am. A/ Vol. 21, No. 8/ August 2004 Isaac Amiror Dot trajectories in the superposition of ranom screens: analysis an synthesis Isaac Amiror Laboratoire e Systèmes Périphériques, Ecole
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationCalculus and optimization
Calculus an optimization These notes essentially correspon to mathematical appenix 2 in the text. 1 Functions of a single variable Now that we have e ne functions we turn our attention to calculus. A function
More informationSolutions to Math 41 Second Exam November 4, 2010
Solutions to Math 41 Secon Exam November 4, 2010 1. (13 points) Differentiate, using the metho of your choice. (a) p(t) = ln(sec t + tan t) + log 2 (2 + t) (4 points) Using the rule for the erivative of
More informationStudents need encouragement. So if a student gets an answer right, tell them it was a lucky guess. That way, they develop a good, lucky feeling.
Chapter 8 Analytic Functions Stuents nee encouragement. So if a stuent gets an answer right, tell them it was a lucky guess. That way, they evelop a goo, lucky feeling. 1 8.1 Complex Derivatives -Jack
More information18 EVEN MORE CALCULUS
8 EVEN MORE CALCULUS Chapter 8 Even More Calculus Objectives After stuing this chapter you shoul be able to ifferentiate an integrate basic trigonometric functions; unerstan how to calculate rates of change;
More informationSection 2.7 Derivatives of powers of functions
Section 2.7 Derivatives of powers of functions (3/19/08) Overview: In this section we iscuss the Chain Rule formula for the erivatives of composite functions that are forme by taking powers of other functions.
More informationAnalytic Scaling Formulas for Crossed Laser Acceleration in Vacuum
October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationMathematics 1210 PRACTICE EXAM II Fall 2018 ANSWER KEY
Mathematics 1210 PRACTICE EXAM II Fall 2018 ANSWER KEY 1. Calculate the following: a. 2 x, x(t) = A sin(ωt φ) t2 Solution: Using the chain rule, we have x (t) = A cos(ωt φ)ω = ωa cos(ωt φ) x (t) = ω 2
More informationMany problems in physics, engineering, and chemistry fall in a general class of equations of the form. d dx. d dx
Math 53 Notes on turm-liouville equations Many problems in physics, engineering, an chemistry fall in a general class of equations of the form w(x)p(x) u ] + (q(x) λ) u = w(x) on an interval a, b], plus
More informationCIRCULAR MOTION AND SHM: Solutions to Higher Level Questions
CIRCULAR MOTION AND SHM: Solutions to Higher Level Questions ****ALL QUESTIONS-ANSWERS ARE HIGHER LEVEL*** Solutions 015 Question 6 (i) Explain what is meant by centripetal force. The force - acting in
More informationMath 1272 Solutions for Spring 2005 Final Exam. asked to find the limit of the sequence. This is equivalent to evaluating lim. lim.
Math 7 Solutions for Spring 5 Final Exam ) We are gien an infinite sequence for which the general term is a n 3 + 5n n + n an are 3 + 5n aske to fin the limit of the sequence. This is equialent to ealuating
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More informationCalculus BC Section II PART A A GRAPHING CALCULATOR IS REQUIRED FOR SOME PROBLEMS OR PARTS OF PROBLEMS
Calculus BC Section II PART A A GRAPHING CALCULATOR IS REQUIRED FOR SOME PROBLEMS OR PARTS OF PROBLEMS. An isosceles triangle, whose base is the interval from (0, 0) to (c, 0), has its verte on the graph
More informationSTUDENT S COMPANIONS IN BASIC MATH: THE FOURTH. Trigonometric Functions
STUDENT S COMPANIONS IN BASIC MATH: THE FOURTH Trigonometric Functions Let me quote a few sentences at the beginning of the preface to a book by Davi Kammler entitle A First Course in Fourier Analysis
More informationCalculus in the AP Physics C Course The Derivative
Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.
More informationCHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold
CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the
More informationAP Calculus AB One Last Mega Review Packet of Stuff. Take the derivative of the following. 1.) 3.) 5.) 7.) Determine the limit of the following.
AP Calculus AB One Last Mega Review Packet of Stuff Name: Date: Block: Take the erivative of the following. 1.) x (sin (5x)).) x (etan(x) ) 3.) x (sin 1 ( x3 )) 4.) x (x3 5x) 4 5.) x ( ex sin(x) ) 6.)
More informationExamining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing
Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June
More informationNoether s theorem applied to classical electrodynamics
Noether s theorem applie to classical electroynamics Thomas B. Mieling Faculty of Physics, University of ienna Boltzmanngasse 5, 090 ienna, Austria (Date: November 8, 207) The consequences of gauge invariance
More information