Tutorial Test 5 2D welding robot

Size: px
Start display at page:

Download "Tutorial Test 5 2D welding robot"

Transcription

1 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. (a) The equations of motion. The kinematic variables are the angular position θ an the raial coorinate r of the weling tip on the arm. These variables have the following time erivatives of interest, t θ = ω, θ = α, t t r = v, t r = a. The time erivatives of r carry a prime, to enote that they are taken wrt the reference frame rotating with the robotic arm. We now set up the equations of motion for the weling robot. To robot consists of two pieces, the weling tip an the arm, which are rigi boies by themselves but o not form a rigi boy jointly. Thus, we have to set up the equations of motion for these two rigi boies separately. Here is a free boy iagram: F F M M ~ (on arm) -F (on tip) (on tip) (on arm) The moment M exerte by motor 1 an the force F exerte by motor are riving the motion. The irection of M is out of the page an the irection of F is parallel to the arm. The moment M causes the force F an the couple moment M on the tip as follows. As a consequence of M, the arm will unergo angular acceleration. But since the weling tip is confine to the arm, it must accelerate with it. That is, the tip will have a component of linear acceleration perpenicular to the arm, an the force F, exerte by the arm on the tip, will generate it. Then, by Newton s thir law, the tip will react with a force F on the arm. This force generates a moment which reuces M an thus reuces the acceleration. There is one more thing to take into account. If the arm accelerates with an angular acceleration α, then the weling tip accelerates with the same α, since it is not a mass point an it is mounte on the arm. A couple moment M is require to generate this angular acceleration α of the weling tip, an the weling tip will react back to the arm with a couple moment M, reucing the original moment M further. Now, the force EoM for the weling tip is F + F = m a, where we use m to enote both the mass of the tip an the bar, since they are equal. The acceleration a of the tip is not irectly one of the kinematic variables, an we use the relative motion formula for rotating axes to relate it to 1

2 the linear acceleration a in the rotating reference frame. Using the fact that the acceleration of the pinpoint A is a A = 0, this gives F + F = m ( α r ω r + ω v + a ). (1) We now ecompose the vector-value Eq. (1) into its scalar components, beginning with the component perpenicular to the arm. To extract it, we take the cross prouct r Eq. (1). We obtain The moment-eom for the arm reas where l is the length of the arm. The moment-eom for the weling tip reas r F = m ( r α + rv ω ) () M r F M = ml α, (3) 3 M = mk G α, (4) where k G is the raius of gyration for the weling tip. In the last three equations, the force F an the couple moment M are unknown. However, by aing up the Eqns. () - (4), we can eliminate them. We obtain l M = m 3 + r + kg α + rv ω. Since α, ω ˆk, the above is a scalar-value equation (posing as a vector-value equation). We multiply with ˆk an obtain M = m l 3 + r + kg α + rv ω. (5) This is our first equation of motion. Note that the factor multiplying α on the rhs is the combine mass moment of inertia of arm an tip. We now return to Eq. (1) an observe that we have not yet extracte the component parallel to the arm. To obtain it, we multiply Eq. (1) with e r, the unit vector along r. We obtain F = m ( ω r + a ), (6) which is our secon equation of motion. Consistency with angular momentum. The moment M measures the rate of change of the angular momentum L, M = t L. The angular momentum L = I total ω for the given situation is l L = m 3 + r + kg ω. Taking the time erivative of this equation, an noting that the raial position r of the tip is timeepenent (/t r = v ), we obtain l M = m 3 + r + kg α + rv ω,

3 which coincies with Eq. (5). Test OK! (b) The situation now simplifies because the weling tip stays at its outer position at all times, r l, v, a = 0. The remaining equation of motion, Eq. (5) therefore simplifies to 4 M = m 3 l + kg α. }{{} I The optimal strategy uner these circumstances is to first accelerate as har as possible, an then ecelerate as har as possible, subject to the two constraints that 1. The arm accumulates a rotation angle of θ = π uring that motion, an. It arrives at the final position with angular velocity zero. Accelerating phase. The epenence of M on ω is M = M 0 (1 ω/ω 0 ), such that we obtain the following linear ifferential equation for ω. t ω = M ) 0 (1 ωω0. I The solution of this LDE, with bounary conition ω(t = 0) = 0, is where the characteristic time scale τ is Integrating Eq. (7) in time, we obtain ω(t) = ω 0 (1 e t/τ ), (7) θ(t) = ω 0 τ τ = Iω 0 M 0. ( ) e t/τ 1 + t/τ. (8) The accelerating phase ens at a time T 1, whose numerical value remains to be etermine. Decelerating phase. The epenence of M on ω now is M = M 0 ( 1 ω/ω 0 ), such that we obtain the following linear ifferential equation for ω. t ω = M ) 0 ( 1 ωω0. I The solution of this LDE, with bounary conition ω(t = T 1 ) = ω(t 1 ) = ω 0 ( 1 e T 1 /τ ), is Integrating Eq. (9) in time, we obtain (( ω(t + T 1 ) = ω ω(t ) ) 1) e t/τ 1. (9) ω 0 θ(t + T 1 ) = θ(t 1 ) + τ(ω 0 + ω(t 1 )) ( 1 e t/τ ) ω 0 t. (10) The ecelerating phase lasts for a time T, whose numerical value remains to be etermine. θ(t 1 ) on the rhs of the last expression is given by θ(t 1 ) = ω 0 τ ( e T 1/τ 1 + T 1 /τ ), cf Eq. (8). 3

4 Numerical solution. We now solve Eqs. (9) an (10) for T 1, T, such that ω(t 1 + T ) = 0 an θ(t 1 + T ) = π. We fin T 1 =.158 s, T = s. The total uration of the motion thus is T 1 + T =.737 s. The angle covere in the accelerating phase is.685 ra = 154, an the ecelerating phase is ra = 6. (c) Optimal strategy in the general case. A goo strategy is to get the weling tip as quickly as possible to the centre A of the arm, such that the mass moment of inertia T total is minimize, allowing for rapi acceleration an eceleration of the arm. Depening on the relative strength of motors 1 an, however, if motor is weak it may happen that the weling tip oes not reach the central position A in time, an must turn aroun sooner. If motor is strong enough for the weling tip to reach position A, the optimal strategy is: Motor 1: 1. Accelerate har for a time T 1.. Decelerate har for a time T Motor : 1. Accelerate har inwars for a time T 3.. Har ecelerate inwar motion for a time T 4 (arrive at A with v = 0). 3. Stay at the centre A for a time T Accelerate har outwar for a time T Har ecelerate outwar motion for a time T 7. The motion is thus escribe by the seven parameters T 1,..., T 7. Given these 7 parameters, the motion follows from (numerically) integrating the couple equations of motion Eqs. (5), (6). We have the following constraints: 1. ω(t 1 + T ) = 0.. θ(t 1 + T ) = π. 3. v (T 3 + T 4 ) = r(t 3 + T 4 ) = v (T 3 + T 4 + T 5 + T 6 + T 7 ) = r(t 3 + T 4 + T 5 + T 6 + T 7 ) = l. 7. T 1 + T = T 3 + T 4 + T 5 + T 6 + T 7. We thus have 7 constraints for 7 parameters. If the values of the given quantities are right (motor sufficiently strong) we can expect a solution. Marking: (a) 5 points, (b) 5 points, (c) points. 4

5 M M 0 F 0 F v 0 ω0 ω v Motor 1 Motor Figure 1: Characteristics of the two motors. Left: Moment M exerte on the arm by Motor 1, as a function of the angular velocity ω of the arm. Right: Force F exerte on the weling tip by Motor, as a function of the velocity v of the weling tip with respect to the arm. Problem 1. We iscuss the simplifie esign of a weling robot in two imensions; see the figure below. The weling robot consists of an arm (= thin ro) capable of rotating about the en point A, an a weling tip, which can freely move along the ro, from the outermost point of the arm all the way to the center point A. The esign invokes two electric motors, one to turn the arm, an one to move the weling tip along the arm. weling tip A robotic arm (initial position) final position In the initial position (Pos. 1), the arm is at rest (ω = 0), an the weling tip is at rest with respect to the arm (v x y = 0). In the final position (Pos. ), the arm is again at rest (ω = 0), an the weling tip is again at rest with respect to the arm (v x y = 0). Furthermore, the arm in its final position is at a 180 egree angle from the initial position. Information about the motors. The robotic arm is riven by two electric motors, one rotating the arm (motor 1), an one moving the weling tip on the arm inwars or outwars (motor ). The moment which motor 1 generates to rotate the robotic arm epens linearly on the angular frequency of the arm. It is M 0 = 30 Nm at ω = 0, an zero at ω 0 = ra/s. The force which motor generates to move the weling tip in an out on the robotic arm epens linearly on the velocity v of the weling tip with respect to the arm. It is F = 0 N at v = 0 an F = 0 at v 0 = m/s. These characteristics are plotte in Fig. 1. We have to consier what happens when the voltage across Motor 1 an/or Motor is reverse. The effect is the same, an so we isplay it below for Motor 1 only. Namely, what happens is that M 0 M 0 an ω 0 ω 0. For Motor 1, the epenency of M on ω for the reverse voltage (= maximum negative voltage) is plotte in Fig., along with the same epenence for the original (= maximum positive) voltage. Analogously, for Motor, F 0 F 0 an v 0 v 0 uner reversal of the voltage. There is one point worth of note about these characteristics: The motor can be use both as an accelerating evice an as a brake. An, it always brakes harer than it accelerates. To see this, consier the vertical ashe line in Fig.. At that point near ω 0, Motor 1 can only eliver a weak accelerating moment, but, with the maximum negative voltage applie, it elivers a large ecelerating moment. Dimensions, masses an constraints. The robotic arm may be consiere as a thin ro. It has a length of l = 1 m, an a mass of 10 kg. The weling tip also has a mass of 10 kg, an a raius of gyration of 10 cm. The weling tip is free to move on the arm for its whole length l, i.e., the centre of mass of the weling tip may move between the pin point A an the outermost point of the arm. The weling tip must be steere such that it is never slamme into its extreme inner an out position on the arm, but rather arrives at these points with zero relative velocity. 5

6 M M 0 maximum positive voltage -ω 0 ω 0 ω -M 0 maximum negative voltage Figure : M vs. ω characteristics (Motor 1) for maximum positive an maximum negative voltage. Questions. (a) [5 points] Derive an state the equations of motion for this weling robot. Hint & consistency check: Consier the situation where Motor 1 elivers zero moment, an the weling tip is moving inwars. What oes conservation of angular momentum say is going to happen? Is this reflecte in your equations of motion? (b) [5 points] Consier the simple case where Motor is not use at all, i.e., the weling tip stays in its outer position at all times. The motion of the arm begins at the angular position θ = 0 with angular velocity ω = 0, an ens at the angular position θ = 180 with angular velocity ω = 0. Devise an optimal (or goo) strategy for moving the arm. Calculate the time neee to move the arm from the initial to the final position. Hint: Analytical arguments go a long way, but in the en you will nee to recur to numerics. (c) [5 points] Now using both motors, evise an optimal (or goo) strategy for moving the arm from its initial position at θ = 0 to its final position at θ = 180. The weling tip starts an ens in the outermost position on the arm, with velocity zero. Explain the working principles behin your strategy. How many free parameters oes your strategy have? How many bounary conitions nee to be matche? Can they be matche by your strategy? ( ) [5 extra points] For the strategy evise in part (c), compute the time neee to move the robotic arm from the initial to the final position. Hint: Rely on numerics. 6

Numerical Integrator. Graphics

Numerical 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 information

Vectors in two dimensions

Vectors 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 information

6. Friction and viscosity in gasses

6. 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 information

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity

1. 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 information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling 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

More information

Statics. There are four fundamental quantities which occur in mechanics:

Statics. There are four fundamental quantities which occur in mechanics: Statics Mechanics isabranchofphysicsinwhichwestuythestate of rest or motion of boies subject to the action of forces. It can be ivie into two logical parts: statics, where we investigate the equilibrium

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 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 information

Quantum Mechanics in Three Dimensions

Quantum 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 information

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012

Free 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 information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 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 information

Schrödinger s equation.

Schrö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 information

Department of Physics University of Maryland College Park, Maryland. Fall 2005 Final Exam Dec. 16, u 2 dt )2, L = m u 2 d θ, ( d θ

Department of Physics University of Maryland College Park, Maryland. Fall 2005 Final Exam Dec. 16, u 2 dt )2, L = m u 2 d θ, ( d θ Department of Physics University of arylan College Park, arylan PHYSICS 4 Prof. S. J. Gates Fall 5 Final Exam Dec. 6, 5 This is a OPEN book examination. Rea the entire examination before you begin to work.

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture 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 information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The 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 information

Cable holds system BUT at t=0 it breaks!! θ=20. Copyright Luis San Andrés (2010) 1

Cable holds system BUT at t=0 it breaks!! θ=20. Copyright Luis San Andrés (2010) 1 EAMPLE # for MEEN 363 SPRING 6 Objectives: a) To erive EOMS of a DOF system b) To unerstan concept of static equilibrium c) To learn the correct usage of physical units (US system) ) To calculate natural

More information

An inductance lookup table application for analysis of reluctance stepper motor model

An inductance lookup table application for analysis of reluctance stepper motor model ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR

More information

Short Intro to Coordinate Transformation

Short 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 information

4. Important theorems in quantum mechanics

4. 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 information

5-4 Electrostatic Boundary Value Problems

5-4 Electrostatic Boundary Value Problems 11/8/4 Section 54 Electrostatic Bounary Value Problems blank 1/ 5-4 Electrostatic Bounary Value Problems Reaing Assignment: pp. 149-157 Q: A: We must solve ifferential equations, an apply bounary conitions

More information

Table of Common Derivatives By David Abraham

Table 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 information

Physics 170 Week 7, Lecture 2

Physics 170 Week 7, Lecture 2 Physics 170 Week 7, Lecture 2 http://www.phas.ubc.ca/ goronws/170 Physics 170 203 Week 7, Lecture 2 1 Textbook Chapter 12:Section 12.2-3 Physics 170 203 Week 7, Lecture 2 2 Learning Goals: Learn about

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

PERMANENT MAGNETS CHAPTER MAGNETIC POLES AND BAR MAGNETS

PERMANENT MAGNETS CHAPTER MAGNETIC POLES AND BAR MAGNETS CHAPTER 6 PERAET AGET 6. AGETIC POLE AD BAR AGET We have seen that a small current-loop carrying a current i, prouces a magnetic fiel B o 4 ji ' at an axial point. Here p ia is the magnetic ipole moment

More information

Euler equations for multiple integrals

Euler 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 information

Chapter 2 Governing Equations

Chapter 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 information

The Principle of Least Action and Designing Fiber Optics

The 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 information

Euler Equations: derivation, basic invariants and formulae

Euler Equations: derivation, basic invariants and formulae Euler Equations: erivation, basic invariants an formulae Mat 529, Lesson 1. 1 Derivation The incompressible Euler equations are couple with t u + u u + p = 0, (1) u = 0. (2) The unknown variable is the

More information

1.2 - Stress Tensor Marine Hydrodynamics Lecture 3

1.2 - Stress Tensor Marine Hydrodynamics Lecture 3 13.021 Marine Hyroynamics, Fall 2004 Lecture 3 Copyright c 2004 MIT - Department of Ocean Engineering, All rights reserve. 1.2 - Stress Tensor 13.021 Marine Hyroynamics Lecture 3 Stress Tensor τ ij:. The

More information

a) Identify the kinematical constraint relating motions Y and X. The cable does NOT slip on the pulley. For items (c) & (e-f-g) use

a) Identify the kinematical constraint relating motions Y and X. The cable does NOT slip on the pulley. For items (c) & (e-f-g) use EAMPLE PROBLEM for MEEN 363 SPRING 6 Objectives: a) To erive EOMS of a DOF system b) To unerstan concept of static equilibrium c) To learn the correct usage of physical units (US system) ) To calculate

More information

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1 Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems

More information

23 Implicit differentiation

23 Implicit differentiation 23 Implicit ifferentiation 23.1 Statement The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x an y. If a value of x is given, then a corresponing value of y is etermine. For

More information

Fluid Mechanics EBS 189a. Winter quarter, 4 units, CRN Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 Veihmeyer Hall.

Fluid Mechanics EBS 189a. Winter quarter, 4 units, CRN Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 Veihmeyer Hall. Flui Mechanics EBS 189a. Winter quarter, 4 units, CRN 52984. Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 eihmeyer Hall. Course Description: xioms of flui mechanics, flui statics,

More information

Calculus of Variations

Calculus 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 information

Experimental Determination of Mechanical Parameters in Sensorless Vector-Controlled Induction Motor Drive

Experimental Determination of Mechanical Parameters in Sensorless Vector-Controlled Induction Motor Drive Experimental Determination of Mechanical Parameters in Sensorless Vector-Controlle Inuction Motor Drive V. S. S. Pavan Kumar Hari, Avanish Tripathi 2 an G.Narayanan 3 Department of Electrical Engineering,

More information

Related Rates. Introduction. We are familiar with a variety of mathematical or quantitative relationships, especially geometric ones.

Related Rates. Introduction. We are familiar with a variety of mathematical or quantitative relationships, especially geometric ones. Relate Rates Introuction We are familiar with a variety of mathematical or quantitative relationships, especially geometric ones For example, for the sies of a right triangle we have a 2 + b 2 = c 2 or

More information

Prep 1. Oregon State University PH 213 Spring Term Suggested finish date: Monday, April 9

Prep 1. Oregon State University PH 213 Spring Term Suggested finish date: Monday, April 9 Oregon State University PH 213 Spring Term 2018 Prep 1 Suggeste finish ate: Monay, April 9 The formats (type, length, scope) of these Prep problems have been purposely create to closely parallel those

More information

Experimental Robustness Study of a Second-Order Sliding Mode Controller

Experimental Robustness Study of a Second-Order Sliding Mode Controller Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans

More information

Calculus in the AP Physics C Course The Derivative

Calculus 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 information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math 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 information

Prof. Dr. Ibraheem Nasser electric_charhe 9/22/2017 ELECTRIC CHARGE

Prof. Dr. Ibraheem Nasser electric_charhe 9/22/2017 ELECTRIC CHARGE ELECTRIC CHARGE Introuction: Orinary matter consists of atoms. Each atom consists of a nucleus, consisting of protons an neutrons, surroune by a number of electrons. In electricity, the electric charge

More information

A Second Time Dimension, Hidden in Plain Sight

A Second Time Dimension, Hidden in Plain Sight A Secon Time Dimension, Hien in Plain Sight Brett A Collins. In this paper I postulate the existence of a secon time imension, making five imensions, three space imensions an two time imensions. I will

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE 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 information

Derivation of angular momentum balance law using the cauchy stress tensor measure. (HW#4, MAE 295. UCI)

Derivation of angular momentum balance law using the cauchy stress tensor measure. (HW#4, MAE 295. UCI) Derivation of angular momentum balance law using the cauchy stress tensor measure. (HW#4, MAE 295. UCI) by Nasser Abbasi February 28, 2006 Problem Derive the angular momentum balance (AMB) equation for

More information

Homework 7 Due 18 November at 6:00 pm

Homework 7 Due 18 November at 6:00 pm Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine

More information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

Analysis of Instantaneous Center of Zero Acceleration of Rigid Body in Planar Motion

Analysis of Instantaneous Center of Zero Acceleration of Rigid Body in Planar Motion Moern pplie Science pril, 009 nalysis of Instantaneous Center of Zero cceleration of Rigi Boy in Planar Motion Haibin Sun (Corresponing author) & Tingting Liu Department of Physics an Electronics Taishan

More information

2.20 Marine Hydrodynamics Lecture 3

2.20 Marine Hydrodynamics Lecture 3 2.20 Marine Hyroynamics, Fall 2018 Lecture 3 Copyright c 2018 MIT - Department of Mechanical Engineering, All rights reserve. 1.7 Stress Tensor 2.20 Marine Hyroynamics Lecture 3 1.7.1 Stress Tensor τ ij

More information

6 Wave equation in spherical polar coordinates

6 Wave equation in spherical polar coordinates 6 Wave equation in spherical polar coorinates We now look at solving problems involving the Laplacian in spherical polar coorinates. The angular epenence of the solutions will be escribe by spherical harmonics.

More information

Interpolated Rigid-Body Motions and Robotics

Interpolated 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 information

Inverting: Representing rotations and translations between coordinate frames of reference. z B. x B x. y B. v = [ x y z ] v = R v B A. y B.

Inverting: Representing rotations and translations between coordinate frames of reference. z B. x B x. y B. v = [ x y z ] v = R v B A. y B. Kinematics Kinematics: Given the joint angles, comute the han osition = Λ( q) Inverse kinematics: Given the han osition, comute the joint angles to attain that osition q = Λ 1 ( ) s usual, inverse roblems

More information

Math 115 Section 018 Course Note

Math 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 information

A new identification method of the supply hole discharge coefficient of gas bearings

A new identification method of the supply hole discharge coefficient of gas bearings Tribology an Design 95 A new ientification metho of the supply hole ischarge coefficient of gas bearings G. Belforte, F. Colombo, T. Raparelli, A. Trivella & V. Viktorov Department of Mechanics, Politecnico

More information

Section 2.7 Derivatives of powers of functions

Section 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 information

6.003 Homework #7 Solutions

6.003 Homework #7 Solutions 6.003 Homework #7 Solutions Problems. Secon-orer systems The impulse response of a secon-orer CT system has the form h(t) = e σt cos(ω t + φ)u(t) where the parameters σ, ω, an φ are relate to the parameters

More information

Introduction to the Vlasov-Poisson system

Introduction 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 information

ACCELERATION, FORCE, MOMENTUM, ENERGY : solutions to higher level questions

ACCELERATION, 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 information

Math 11 Fall 2016 Section 1 Monday, September 19, Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v

Math 11 Fall 2016 Section 1 Monday, September 19, Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v Math Fall 06 Section Monay, September 9, 06 First, some important points from the last class: Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v passing through

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic 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 information

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2 Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z

More information

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory

More information

The Principle of Least Action

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 information

The Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis

The 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 information

Designing Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations

Designing Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations EECS 6B Designing Information Devices an Systems II Fall 07 Note 3 Secon Orer Differential Equations Secon orer ifferential equations appear everywhere in the real worl. In this note, we will walk through

More information

Solutions to MATH 271 Test #3H

Solutions to MATH 271 Test #3H Solutions to MATH 71 Test #3H This is the :4 class s version of the test. See pages 4 7 for the 4:4 class s. (1) (5 points) Let a k = ( 1)k. Is a k increasing? Decreasing? Boune above? Boune k below? Convergant

More information

12.11 Laplace s Equation in Cylindrical and

12.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 information

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.

Lecture 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 information

M E 320 Professor John M. Cimbala Lecture 16

M E 320 Professor John M. Cimbala Lecture 16 M E 320 Professor John M. Cimbala Lecture 16 Toay, we will: Do some more example problems lear CV momentum equation Discuss the control volume equation for angular momentum E. The Lear Momentum Equation

More information

Physics 115C Homework 4

Physics 115C Homework 4 Physics 115C Homework 4 Problem 1 a In the Heisenberg picture, the ynamical equation is the Heisenberg equation of motion: for any operator Q H, we have Q H = 1 t i [Q H,H]+ Q H t where the partial erivative

More information

Static Equilibrium. Theory: The conditions for the mechanical equilibrium of a rigid body are (a) (b)

Static 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 information

MEEN 363. EXAMPLE of ANALYSIS (1 DOF) Luis San Andrés

MEEN 363. EXAMPLE of ANALYSIS (1 DOF) Luis San Andrés MEEN 363. EAMPLE of ANALYSIS (1 DOF) Luis San Anrés Objectives: a) To erive EOM for a 1-DOF (one egree of freeom) system b) To unerstan concept of static equilibrium c) To learn the correct usage of physical

More information

Experiment 2, Physics 2BL

Experiment 2, Physics 2BL Experiment 2, Physics 2BL Deuction of Mass Distributions. Last Upate: 2009-05-03 Preparation Before this experiment, we recommen you review or familiarize yourself with the following: Chapters 4-6 in Taylor

More information

Objective: To introduce the equations of motion and describe the forces that act upon the Atmosphere

Objective: To introduce the equations of motion and describe the forces that act upon the Atmosphere Objective: To introuce the equations of motion an escribe the forces that act upon the Atmosphere Reaing: Rea pp 18 6 in Chapter 1 of Houghton & Hakim Problems: Work 1.1, 1.8, an 1.9 on pp. 6 & 7 at the

More information

Lectures HD#14 Dynamic response of continuous systems. Date: April

Lectures HD#14 Dynamic response of continuous systems. Date: April ectures -3 Date: April 4 7 Toay: Vibrations of continuous systems HD#4 Dynamic response of continuous systems Free vibrations of elastic bars an beams. Properties of normal moe functions. Force response

More information

ECE 422 Power System Operations & Planning 7 Transient Stability

ECE 422 Power System Operations & Planning 7 Transient Stability ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3 Transient Stability The ability of

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization 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 information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

ECE 6310 Spring 2012 Exam 1 Solutions. Balanis The electric fields are given by. E r = ˆxe jβ 0 z

ECE 6310 Spring 2012 Exam 1 Solutions. Balanis The electric fields are given by. E r = ˆxe jβ 0 z ECE 6310 Spring 2012 Exam 1 Solutions Balanis 1.30 The electric fiels are given by E i ˆxe jβ 0 z E r ˆxe jβ 0 z The curl of the electric fiels are the usual cross prouct E i jβ 0 ẑ ˆxe jβ 0 z jβ 0 ŷe

More information

2.5 SOME APPLICATIONS OF THE CHAIN RULE

2.5 SOME APPLICATIONS OF THE CHAIN RULE 2.5 SOME APPLICATIONS OF THE CHAIN RULE The Chain Rule will help us etermine the erivatives of logarithms an exponential functions a x. We will also use it to answer some applie questions an to fin slopes

More information

Applications of the Wronskian to ordinary linear differential equations

Applications of the Wronskian to ordinary linear differential equations Physics 116C Fall 2011 Applications of the Wronskian to orinary linear ifferential equations Consier a of n continuous functions y i (x) [i = 1,2,3,...,n], each of which is ifferentiable at least n times.

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 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 information

Phys 7221, Fall 2006: Homework # 5

Phys 7221, Fall 2006: Homework # 5 Phys 722, Fal006: Homewor # 5 Gabriela González October, 2006 Prob 3-: Collapse of an orbital system Consier two particles falling into each other ue to gravitational forces, starting from rest at a istance

More information

VIBRATIONS OF A CIRCULAR MEMBRANE

VIBRATIONS OF A CIRCULAR MEMBRANE VIBRATIONS OF A CIRCULAR MEMBRANE RAM EKSTROM. Solving the wave equation on the isk The ynamics of vibrations of a two-imensional isk D are given by the wave equation..) c 2 u = u tt, together with the

More information

for which the garbage will fall within the can. If the range is large, the shot is easy; if the range is small, the target is hard to hit.

for which the garbage will fall within the can. If the range is large, the shot is easy; if the range is small, the target is hard to hit. EN40: Dynamics an Vibrations Homework : Soling equations of motion for particles Due Friay Feb 18 011 School of Engineering Brown Uniersity 1. Suppose that you wish to throw a piece of garbage into a garbage

More information

FLUID MECHANICS UNIVERSITY OF LEEDS. May/June Examination for the degree of. BEng/ MEng Civil Engineering. Time allowed: 2 hours

FLUID MECHANICS UNIVERSITY OF LEEDS. May/June Examination for the degree of. BEng/ MEng Civil Engineering. Time allowed: 2 hours This question paper consists of printe pages, each of which is ientifie by the Coe Number CIVE 4 UNIVERSITY OF LEEDS May/June Examination for the egree of BEng/ MEng Civil Engineering FLUID MECANICS Time

More information

Problem Set 2: Solutions

Problem Set 2: Solutions UNIVERSITY OF ALABAMA Department of Physics an Astronomy PH 102 / LeClair Summer II 2010 Problem Set 2: Solutions 1. The en of a charge rubber ro will attract small pellets of Styrofoam that, having mae

More information

TEST 2 (PHY 250) Figure Figure P26.21

TEST 2 (PHY 250) Figure Figure P26.21 TEST 2 (PHY 250) 1. a) Write the efinition (in a full sentence) of electric potential. b) What is a capacitor? c) Relate the electric torque, exerte on a molecule in a uniform electric fiel, with the ipole

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

Further Differentiation and Applications

Further 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 information

1 Lecture 20: Implicit differentiation

1 Lecture 20: Implicit differentiation Lecture 20: Implicit ifferentiation. Outline The technique of implicit ifferentiation Tangent lines to a circle Derivatives of inverse functions by implicit ifferentiation Examples.2 Implicit ifferentiation

More information

Topic 2.3: The Geometry of Derivatives of Vector Functions

Topic 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 information

Lecture 2 - First order linear PDEs and PDEs from physics

Lecture 2 - First order linear PDEs and PDEs from physics 18.15 - Introuction to PEs, Fall 004 Prof. Gigliola Staffilani Lecture - First orer linear PEs an PEs from physics I mentione in the first class some basic PEs of first an secon orer. Toay we illustrate

More information

SELF-ERECTING, ROTARY MOTION INVERTED PENDULUM Quanser Consulting Inc.

SELF-ERECTING, ROTARY MOTION INVERTED PENDULUM Quanser Consulting Inc. SELF-ERECTING, ROTARY MOTION INVERTED PENDULUM Quanser Consulting Inc. 1.0 SYSTEM DESCRIPTION The sel-erecting rotary motion inverte penulum consists of a rotary servo motor system (SRV-02) which rives

More information

CIRCULAR MOTION AND SHM: Solutions to Higher Level Questions

CIRCULAR 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 information

G j dq i + G j. q i. = a jt. and

G 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 information

1.4.3 Elementary solutions to Laplace s equation in the spherical coordinates (Axially symmetric cases) (Griffiths 3.3.2)

1.4.3 Elementary solutions to Laplace s equation in the spherical coordinates (Axially symmetric cases) (Griffiths 3.3.2) 1.4.3 Elementary solutions to Laplace s equation in the spherical coorinates (Axially symmetric cases) (Griffiths 3.3.) In the spherical coorinates (r, θ, φ), the Laplace s equation takes the following

More information

CENTURION UNIVERSITY OF TECHNOLOGY & MANAGEMENT,ODISHA CUEE-2015

CENTURION UNIVERSITY OF TECHNOLOGY & MANAGEMENT,ODISHA CUEE-2015 CENTURION UNIVERSITY OF TECHNOLOGY & MANAGEMENT,ODISHA CUEE-015 PHYSICS 1. The imensional formula of angular momentum is a) ML T - b) MLT - c) MLT -1 ) ML T -1. If A B = B A, then the angle between A an

More information

arxiv: v2 [math.ds] 26 Apr 2017

arxiv: v2 [math.ds] 26 Apr 2017 GAUSSIAN CURVATURE AND GYROSCOPES GRAHAM COX AND MARK LEVI arxiv:1607.03217v2 [math.ds] 26 Apr 2017 Abstract. We relate Gaussian curvature to the gyroscopic force, thus giving a mechanical interpretation

More information

and from it produce the action integral whose variation we set to zero:

and 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 information

Problem 1 (20 points)

Problem 1 (20 points) ME 309 Fall 01 Exam 1 Name: C Problem 1 0 points Short answer questions. Each question is worth 5 points. Don t spen too long writing lengthy answers to these questions. Don t use more space than is given.

More information

Center of Gravity and Center of Mass

Center of Gravity and Center of Mass Center of Gravity an Center of Mass 1 Introuction. Center of mass an center of gravity closely parallel each other: they both work the same way. Center of mass is the more important, but center of gravity

More information