Reference Trajectory Generation for an Atomic Force Microscope

Size: px
Start display at page:

Download "Reference Trajectory Generation for an Atomic Force Microscope"

Transcription

1 Mechatronics 2010, June 28-30, Swiss Federal Institute of Technology ETH, Zurich Switzerland 1 Reference Trajectory Generation for an Atomic Force Microscope Ivan Furlan and Silvano Balemi 1 Department of innovative technologies, University of Applied Sciences of Southern Switzerland (SUPSI), 6928 Lugano-Manno, Switzerland. ABSTRACT This paper addresses the design of the scanning trajectory for a portable Atomic Force Microscope with large scanning range and high resolution. The proposed reference scanning trajectory presents a low motion bandwidth and is well suited for a discrete-time implementation. The dependence of the trajectory parameters in function of the scanning speed for different scaling policies, has been analyzed. An optimal choice of scaling policy is then proposed. 1. Introduction Atomic force microscopy is typically performed by scanning the considered surface forwards and backwards along many parallel lines while a constant distance between the scanning head and the surface is maintained (see Fig. 1) y ais (m) ais (m) 10 4 Figure. 1: Trajectory for scanning a surface (units in meters). In contact mode, the constant distance between the scanning head and the surface is obtained by controlling the deflection ξ of the cantilever, which remains in contact with the surface. A sketch of the motion in vertical direction (when travelling along a line on the surface) is shown in Fig. 2. The Atomic Force Microscope (AFM) considered in this work moves along the trajectory thanks to a two-degree-offreedom platform based on fleures and actuated by voice-coil motors [1, 2]. In the desired trajectory the scanning head accelerates until the beginning of the first scanning line, then proceeds forwards at a constant scanning speed until the end of the scanning line, thes direction and travels again at the same speed along the same scanning line but backwards until the beginning of the scanning line, where the heads performs simultaneously a turn and a move sidewards to the net scanning line and so on until the last scanning line. This trajectory differs from those considered in other works (e.g. [3]) in that only the transitions between scanning lines can be designed, while the trajectory along the scanning lines is determined by the constant scanning speed. The acceleration/deceleration and the move to the net line must occur rapidly (in order to reduce overall measurement time) but also smoothly (in order to reduce the perturbations of the motion along the scanning line). Moreover, 1 Corresponding author: silvano.balemi@supsi.ch

2 2 12th Mechatronics Forum Biennial International Conference.. z ξ = z w. Cantilever tip. Cantilever w Figure. 2: Atomic force microscope scanning the surface along a line: ξ cantilever deflection, w surface profile, z head elevation. the motion between scanning lines must be easily computable and easily scaled both in size (scanning of surfaces of different sizes) as well as in speed (scanning with different speeds). In section 2 a simple trajectory is proposed which guarantees a smooth transition between scanning lines. The values of the different trajectory parameters are given as a function of the scanning speed v and of a parameter ω. In section 3 the dependence of the trajectory parameters as a function of the scanning speed for different scaling policies is analyzed. An optimal choice of scaling policy is proposed. In section 4 implementation issues are discussed. Conclusions are given in section Trajectory definition and design 2.1. Trajectory definition In order to guarantee smooth transitions between scanning lines the trajectory and its derivatives should be continuous. The ideal choice would include a continuous jerk (fourth derivative of the trajectory). However, the realization of such a trajectory would require a large amount of computations, both off-line (computation of the trajectory parameters) and on-line (eecution of the trajectory). The trajectory chosen in this work only guarantees a continuous acceleration but not a continuous jerk. The transition between two lines is described by the following equation:» (t) y(t)» = R sin (ω t) v side t r sin(2 ω t)» o + y o The trajectory begins at the position ( o, y o) (here for a simpler analysis at time t = 0) and ends at position ( o, y o + v side t turn) where t turn = /ω is the time needed to turn and to move to the net line. In the sequel the transition from one scanning line to the net one is analyzed. Later the other motions before or after a scanning line will be considered Turn and move sidewards When the head changes direction and moves to another scanning line, the speed at the end of the scanning line is v while the initial speed along the net scanning line should be v. Thus the following conditions must be satisfied (1) (0) = o, (t turn) = o ẋ(0) = v, ẋ(t turn) = v y(0) = y o, y(t turn) = y o + ẏ(0) = 0, ẏ(t turn) = 0 (2) where is the distance between the two scanning lines. With the motion equation (1) and conditions (2), the following relations are obtained From equation (3) and with t turn = /ω the mean speed v side is v side t turn = (3) R ω = v (4) v side 2 ω r = 0 (5) v side = = ω t turn (6)

3 Mechatronics 2010, June 28-30, Swiss Federal Institute of Technology ETH, Zurich Switzerland 3 From equations (4), (5) and (6) the radii R and r are R = v ω r = 2 (7) (8) This means that r only depends on the distance between two scanning lines. A sample trajectory during a turn/move motion and the respective speed are shown in Fig. 3. Note that the acceleration at the beginning and at the end of the motion is zero ,y [m] d/dt, dy/dt [m/s] t [s] t [s] Figure. 3: Left: trajectory between two lines: ais (continuous line) and y ais (dashed line). Right: speed between two lines: ais (continuous line) and y ais (dashed line) Turn At the end of a scanning line the head must be decelerated and accelerated backwards. The motion is similar to that of the case Turn and move sidewards, but with the parameters r = 0 and R being negative (no motion along the y ais) Accelerate At the beginning of the first scanning line the head must be accelerated. The motion is similar to that of the case Turn and move sidewards, but with the parameter r = 0 and with half of the motion duration (from t = t turn/2 to t = t turn) Decelerate At the end of the last scanning line the head must be stopped (decelerated). The motion is similar to that of the case Turn and move sidewards, but with the parameter r = 0 and with half of the motion duration (from t = 0 to t = t turn/2). The whole trajectory can thus be generated in a simple manner: the head travels at a constant speed on the scanning lines, the transitions are all governed by the same formula (1) with the necessary parametrizations Trajectory design In the trajectory design the following limitation on the maimal acceleration along the /y aes must be considered ẍ(t) a ma ÿ(t) a ma where a ma is the maimal available acceleration provided by the head actuators. From this limitation the following conditions are obtained ẍ(t) = R ω 2 a ma ÿ(t) = 4 r ω 2 (9) a ma

4 4 12th Mechatronics Forum Biennial International Conference 3. Analysis of the Scaling policies as a function of the scanning speed In this section the evolution of the different trajectory parameters is analyzed for different values of the scanning speed v. Of course the constraints of equation (9) must be satisfied. Note that the maimal acceleration a ma in equation (9) is proportional to the magnetic field generated by the voice-coil motors and thus proportional to the motor currents. However, the currents are limited by the power supply voltage. Therefore, the relation between the power supply voltage and a ma resp. ω has also to be determined Power supply voltage limitations In this section the actuation voltage of the drivers necessary for following the desired trajectory is determined. In order to do this the motion equations for the different parts of the trajectory are determined (accelerate: ACC; forward movement: FW; turn direction: TD; backward movement: BW and move sidewards: MS). Only the first scanning line is considered here. ACC(t) = R cos (ω t) + first FW(t) = R ω `t ω + first TD(t) = R sin `ω `t t ` ω line + FW + t ω line BW(t) = R ω `t 3 t ` ω line + FW + t ω line MS(t) = R sin `ω `t 3 2 t ω line + first y ACC(t) = y FW(t) = y TD(t) = y BW(t) = y first y MS(t) = v side t sin (2 ω t MS) + y first where first, y first = first point of the scanning surface [m] t line = time for covering a scanning line [s] t MS = t 3 ω line If the damping is assumed to be negligible, the relationship between the trajectory in (and similarly in y) and the driver actuation voltage is given by U (s) = X (s) K F2i `s 2 m + k f (s LM + R M) (10) where K F2i = Force to current constant [A/N] 8>< k f = Fleure spring constant [N/m] m = Moving mass [Kg] L M = inductance of the motor [H] >: R M = resistance of the motor [Ω] With equation (10) and the equations for the whole trajectory along the and y aes, the applied voltage at the voice-coil actuators as a function of time can be determined. The maimal supply voltage is needed when the scanning head turns and moves sidewards. The equations for the maimal voltages in this case are u (t) = K F2i R `k f m ω 2 (L M ω cos (ω t) + R M sin(ω t)) + K F2i k f R M u y(t) = K F2i R `4 m ω 2 k f (2 LM ω cos (ω t) + R M sin(ω t)) +K F2i R M k f v side t + y K F2i k f R M + L M K F2i k f v side From these equations, approimated formulas showing the general dependence of the supply voltage from the variable ω can be derived u (t) K F2i L M m R ω 3 u y(t) 8 m KF2i ω Use of maimal acceleration The atomic force microscope must be able to scan large and small surfaces at different scanning speeds (fast speed for smooth surfaces, slow speed for rough surfaces). Then, in order to uniquely determine the trajectory for a given v, one more variable, for instance the parameter ω, must be specified. From (9) and equation (7) the new epression v ω = a

5 Mechatronics 2010, June 28-30, Swiss Federal Institute of Technology ETH, Zurich Switzerland 5 Table 1: Dependencies of quantities in function of scanning speed v and of the maimal acceleration a = k v n. a t turn ω R û û y a y v y k v n k v n+1 k v n 1 1 k v n+2 k v 2n 1 k v 3n 3 2 k v 2n 2 2 k v n 1 indicates that, as an alternative, the maimal acceleration a = ẍ(t) can be specified (which i determines the angular speed ω). A simple choice consists in always eploiting the maimal acceleration/deceleration (or equivalently, the maimal force) in order to reduce the overall measurement time. If the maimal acceleration a = ẍ(t) is assumed to be the same for all scanning speeds v, then ω = ẍ(t) /v increases for decreasing scanning speed v. Also, it means that the side motion speed v side = 2 ẍ(t) /(v ) increases linearly with decreasing v (see equation (6)) but also that the necessary voltage û y = ẍ(t) 3 /v 3 for the side motion increases with a power of three! This policy is obviously useless in practice Policy analysis In order to define uniquely the trajectory for a given value of v, different dependences of the acceleration a in function of the speed v can be analyzed. In particular the parametrization a = k v n in function of k and n is investigated. In Table 1 the epression for various quantities in function of v and of the parameters k and n are given. Different cases can be identified: n = 0: Constant maimal acceleration ẍ(t), n = 1/2: Constant maimal voltage û, n = 1: Constant angular speed ω, constant turn time, constant maimal lateral speed and acceleration for all speeds v, n = 2: Constant radius R. Because the scanning speed may be arbitrarily high, no quantity (voltage) should increase with decreasing speed, that is no quantity in Table 1 may show negative eponent of v. The conditions on ω, û, û y, a y, v y, a are all satisfied with n 1, the conditions on t turn and R are satisfied for n 1, which leads to n = Implementation The generation of the scanning trajectory was implemented in a digital controller at high sampling frequencies (40kHz). Thus an efficient discrete-time implementation had to be developed Generation of the trajectory In order to generate the parts of trajectory describing the transition between two scanning lines or the other motions before or after a scanning line, equation (1) was implemented using the discrete-time state-space equations of the two second order oscillator systems s s = Φ s y s k+1 y k sq = Φ 2 sq s qy s qy and the epression» y k+1» s = C turn s y k k+1» sq + C side s qy k k» 0 + q o k + y o

6 6 12th Mechatronics Forum Biennial International Conference where Φ = C turn = C side = cos(ω ts) sin(ω t s) sin(ω t s) cos(ω t s) R r and where t s is the chosen sampling time. Ideally, after a given number of samples starting from ( o, y o), the correct initial point of the net scanning line should be reached eactly. This implies that t turn for the sideways and for the simple turn motion but also that t turn/2 for the deceleration and the acceleration motion be a multiple of the sampling time t s (see Fig. 4) 2 R = r ( o,y o ) Then t turn = t s with Figure. 4: Eample of resulting trajectory with discretization of the trajectory. tturn = 2 1 = 2 2 t s 2 ω t s Note that also t turn = t 2 2 s is also a multiple of t s since by construction is even. Then, the circular frequency ω becomes ω = t s and the radius R Finally, the rotation matri Φ becomes 4.2. Specific cases Φ = R = v ω 2 4 cos sin = v ts nturn sin cos For the different cases presented above it suffices to correctly initialize the initial conditions of the states [s, s y] T, [s q, s qy] T to generate the trajectory. Turn and move sidewards 3 5 Turn: Accelerate: [s, s y] T = [0, 1] T [s q, s qy] T = [1, 0] T q = [s, s y] T = [0, 1] T [s q, s qy] T = [0, 0] T q = 0 [s, s y] T = [1, 0] T [s q, s qy] T = [0, 0] T q = 0

7 Mechatronics 2010, June 28-30, Swiss Federal Institute of Technology ETH, Zurich Switzerland 7 Decelerate: [s, s y] T = [0, 1] T [s q, s qy] T = [0, 0] T q = 0 Note that C turn and C side do not change (R and r remain positive) Results The method has been etended in a straightforward manner to the case when then trajectory is rotated with respect to the and y orthogonal aes. A real case is presented in Fig. (5) where the reference trajectory is shown together with the actual trajectory obtained with the AFM. The deviation is due to fairly strong couplings between the aes and to unmodeled plant dynamics. The influence of the non-ideal reference trajectory is not relevant y ais (m) ais (m) Figure. 5: Measurement of the scanning motion of the AFM with the new reference trajectory generation method. Units in meters, scanning time per line about 50ms. 5. Conclusions In this paper a reference trajectory generation method and a scaling policy for adaptations to changes in the scanning speed have been proposed. The method can be easily implemented and requires only simple computations, making it suitable also for hardware implementations (e.g. in FPGAs). References [1] S. Henein, C. Aymon, S. Bottinelli, and R. Clavel. Articulated structures with fleible joints dedicated to high precision robotics. In Proc. of International Advanced Robotics Program; Workshop on Micro Robots, Micro Machines and Systems. Moscow, Russia, November, 1999.

8 [2] S. Balemi, J. Moerschell, J-M. Breguet, D. Brändlin, S. Bottinelli, and I. Beltrami. Surface inspection system for industrial applications. In Proc. of Conf. on Robotics and Mechatronics, pages , Aachen, Germany, September [3] A.J. Fleming and A.G. Wills. Optimal input scaling for bandlimited scanning systems. In Proc. of 17th IFAC World Congress, Seoul, Korea, July [4] S. Devasia, E. Eleftheriou, and S.O.R. Moheimani. A survey of control issues in nanopositioning. IEEE Transactions on Control Systems Technology, 15(5):802823, September

Santosh Devasia Mechanical Eng. Dept., UW

Santosh Devasia Mechanical Eng. Dept., UW Nano-positioning Santosh Devasia Mechanical Eng. Dept., UW http://faculty.washington.edu/devasia/ Outline of Talk 1. Why Nano-positioning 2. Sensors for Nano-positioning 3. Actuators for Nano-positioning

More information

Balancing of an Inverted Pendulum with a SCARA Robot

Balancing of an Inverted Pendulum with a SCARA Robot Balancing of an Inverted Pendulum with a SCARA Robot Bernhard Sprenger, Ladislav Kucera, and Safer Mourad Swiss Federal Institute of Technology Zurich (ETHZ Institute of Robotics 89 Zurich, Switzerland

More information

National Quali cations

National Quali cations National Quali cations AH06 X70/77/ Mathematics of Mechanics TUESDAY, 7 MAY :00 PM :00 PM Total marks 00 Attempt ALL questions. You may use a calculator. Full credit will be given only to solutions which

More information

Date: 1 April (1) The only reference material you may use is one 8½x11 crib sheet and a calculator.

Date: 1 April (1) The only reference material you may use is one 8½x11 crib sheet and a calculator. PH1140: Oscillations and Waves Name: Solutions Conference: Date: 1 April 2005 EXAM #1: D2005 INSTRUCTIONS: (1) The only reference material you may use is one 8½x11 crib sheet and a calculator. (2) Show

More information

Date: 31 March (1) The only reference material you may use is one 8½x11 crib sheet and a calculator.

Date: 31 March (1) The only reference material you may use is one 8½x11 crib sheet and a calculator. PH1140: Oscillations and Waves Name: SOLUTIONS AT END Conference: Date: 31 March 2005 EXAM #1: D2006 INSTRUCTIONS: (1) The only reference material you may use is one 8½x11 crib sheet and a calculator.

More information

Trajectory Planning, Setpoint Generation and Feedforward for Motion Systems

Trajectory Planning, Setpoint Generation and Feedforward for Motion Systems 2 Trajectory Planning, Setpoint Generation and Feedforward for Motion Systems Paul Lambrechts Digital Motion Control (4K4), 23 Faculty of Mechanical Engineering, Control Systems Technology Group /42 2

More information

Trajectory planning and feedforward design for electromechanical motion systems version 2

Trajectory planning and feedforward design for electromechanical motion systems version 2 2 Trajectory planning and feedforward design for electromechanical motion systems version 2 Report nr. DCT 2003-8 Paul Lambrechts Email: P.F.Lambrechts@tue.nl April, 2003 Abstract This report considers

More information

Rotational Motion Control Design for Cart-Pendulum System with Lebesgue Sampling

Rotational Motion Control Design for Cart-Pendulum System with Lebesgue Sampling Journal of Mechanical Engineering and Automation, (: 5 DOI:.59/j.jmea.. Rotational Motion Control Design for Cart-Pendulum System with Lebesgue Sampling Hiroshi Ohsaki,*, Masami Iwase, Shoshiro Hatakeyama

More information

Modeling and Analysis of Dynamic Systems

Modeling and Analysis of Dynamic Systems Modeling and Analysis of Dynamic Systems by Dr. Guillaume Ducard Fall 2016 Institute for Dynamic Systems and Control ETH Zurich, Switzerland based on script from: Prof. Dr. Lino Guzzella 1/33 Outline 1

More information

Introduction to Control (034040) lecture no. 2

Introduction to Control (034040) lecture no. 2 Introduction to Control (034040) lecture no. 2 Leonid Mirkin Faculty of Mechanical Engineering Technion IIT Setup: Abstract control problem to begin with y P(s) u where P is a plant u is a control signal

More information

Discrete-time models and control

Discrete-time models and control Discrete-time models and control Silvano Balemi University of Applied Sciences of Southern Switzerland Zürich, 2009-2010 Discrete-time signals 1 Step response of a sampled system Sample and hold 2 Sampling

More information

Fast Seek Control for Flexible Disk Drive Systems

Fast Seek Control for Flexible Disk Drive Systems Fast Seek Control for Flexible Disk Drive Systems with Back EMF and Inductance Chanat La-orpacharapan and Lucy Y. Pao Department of Electrical and Computer Engineering niversity of Colorado, Boulder, CO

More information

Modeling and Experimentation: Mass-Spring-Damper System Dynamics

Modeling and Experimentation: Mass-Spring-Damper System Dynamics Modeling and Experimentation: Mass-Spring-Damper System Dynamics Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin July 20, 2014 Overview 1 This lab is meant to

More information

Position and Velocity Profile Tracking Control for New Generation Servo Track Writing

Position and Velocity Profile Tracking Control for New Generation Servo Track Writing Preprints of the 9th World Congress The International Federation of Automatic Control Cape Town, South Africa. August 24-29, 24 Position and Velocity Profile Tracking Control for New Generation Servo Track

More information

Physics. Practice Questions IEB

Physics. Practice Questions IEB Physics Practice Questions IEB Instructions Individual, eam-style questions The questions contained in this booklet match the style of questions that are typically asked in eams. This booklet is not however,

More information

Springs: Part I Modeling the Action The Mass/Spring System

Springs: Part I Modeling the Action The Mass/Spring System 17 Springs: Part I Second-order differential equations arise in a number of applications We saw one involving a falling object at the beginning of this text (the falling frozen duck example in section

More information

1. Write the relation for the force acting on a charge carrier q moving with velocity through a magnetic field in vector notation. Using this relation, deduce the conditions under which this force will

More information

Downloaded from

Downloaded from Question 4.1: A circular coil of wire consisting of 100 turns, each of radius 8.0 cm carries a current of 0.40 A. What is the magnitude of the magnetic field B at the centre of the coil? Number of turns

More information

APPM 1360 Final Exam Spring 2016

APPM 1360 Final Exam Spring 2016 APPM 36 Final Eam Spring 6. 8 points) State whether each of the following quantities converge or diverge. Eplain your reasoning. a) The sequence a, a, a 3,... where a n ln8n) lnn + ) n!) b) ln d c) arctan

More information

High speed and quick response precise linear stage system using V-shape transducer ultrasonic motors

High speed and quick response precise linear stage system using V-shape transducer ultrasonic motors High speed and quick response precise linear stage system using V-shape transducer ultrasonic motors K. Asumi, T. Fujimura and M. K. Kurosawa Taiheiyo-cement Corporation, Japan Tokyo Institute of Technology,

More information

Active Nonlinear Observers for Mobile Systems

Active Nonlinear Observers for Mobile Systems Active Nonlinear Observers for Mobile Systems Simon Cedervall and Xiaoming Hu Optimization and Systems Theory Royal Institute of Technology SE 00 44 Stockholm, Sweden Abstract For nonlinear systems in

More information

MODELLING OF VIBRATION WITH ABSORBER

MODELLING OF VIBRATION WITH ABSORBER Journal of Machine Engineering, Vol. 10, No.4, 2010 mumerical modelling, Simulink, vibration, machine, absorber Jiří VONDŘICH 1 Evžen THÖNDE 1 Slavomír JIRKŮ 1 MODEING OF VIBRATION WITH ABSORBER Machine

More information

Motion in Two and Three Dimensions

Motion in Two and Three Dimensions chapter 4 Motion in Two and Three Dimensions Projectile motion (Section 4.3) 1. Which target got hit first? Contet of the tetbook: Before Eample 4. 2. Projectile range problem comparable to Eample 7, ecept

More information

PHYS 241 EXAM #2 November 9, 2006

PHYS 241 EXAM #2 November 9, 2006 1. ( 5 points) A resistance R and a 3.9 H inductance are in series across a 60 Hz AC voltage. The voltage across the resistor is 23 V and the voltage across the inductor is 35 V. Assume that all voltages

More information

Example 1: What do you know about the graph of the function

Example 1: What do you know about the graph of the function Section 1.5 Analyzing of Functions In this section, we ll look briefly at four types of functions: polynomial functions, rational functions, eponential functions and logarithmic functions. Eample 1: What

More information

S13 PHY321: Final May 1, NOTE: Show all your work. No credit for unsupported answers. Turn the front page only when advised by the instructor!

S13 PHY321: Final May 1, NOTE: Show all your work. No credit for unsupported answers. Turn the front page only when advised by the instructor! Name: Student ID: S13 PHY321: Final May 1, 2013 NOTE: Show all your work. No credit for unsupported answers. Turn the front page only when advised by the instructor! The exam consists of 6 problems (60

More information

Linear Shaft Motor Sizing Application Note

Linear Shaft Motor Sizing Application Note Linear Shaft Motor Sizing Application Note By Jeramé Chamberlain One of the most straightforward tasks in the design of a linear motion system is to specify a motor and drive combination that can provide

More information

Nonlinear Oscillations and Chaos

Nonlinear Oscillations and Chaos CHAPTER 4 Nonlinear Oscillations and Chaos 4-. l l = l + d s d d l l = l + d m θ m (a) (b) (c) The unetended length of each spring is, as shown in (a). In order to attach the mass m, each spring must be

More information

Lecture 4 Scanning Probe Microscopy (SPM)

Lecture 4 Scanning Probe Microscopy (SPM) Lecture 4 Scanning Probe Microscopy (SPM) General components of SPM; Tip --- the probe; Cantilever --- the indicator of the tip; Tip-sample interaction --- the feedback system; Scanner --- piezoelectric

More information

Chapter 4 - Moving Charges and Magnetism. Magnitude of the magnetic field at the centre of the coil is given by the relation,

Chapter 4 - Moving Charges and Magnetism. Magnitude of the magnetic field at the centre of the coil is given by the relation, Question 4.1: A circular coil of wire consisting of 100 turns, each of radius 8.0 cm carries a current of 0.40 A. What is the magnitude of the magnetic field B at the centre of the coil? Number of turns

More information

13 Path Planning Cubic Path P 2 P 1. θ 2

13 Path Planning Cubic Path P 2 P 1. θ 2 13 Path Planning Path planning includes three tasks: 1 Defining a geometric curve for the end-effector between two points. 2 Defining a rotational motion between two orientations. 3 Defining a time function

More information

Survey of Methods of Combining Velocity Profiles with Position control

Survey of Methods of Combining Velocity Profiles with Position control Survey of Methods of Combining Profiles with control Petter Karlsson Mälardalen University P.O. Box 883 713 Västerås, Sweden pkn91@student.mdh.se ABSTRACT In many applications where some kind of motion

More information

Errors Intensive Computation

Errors Intensive Computation Errors Intensive Computation Annalisa Massini - 2015/2016 OVERVIEW Sources of Approimation Before computation modeling empirical measurements previous computations During computation truncation or discretization

More information

224 g (Total 1 mark) A straight line in a direction opposite to that of the field.

224 g (Total 1 mark) A straight line in a direction opposite to that of the field. Q1. The diagram shows a rigidly-clamped straight horizontal current-carrying wire held mid-way between the poles of a magnet on a top pan balance. The wire is perpendicular to the magnetic field direction.

More information

Mechanics II. Which of the following relations among the forces W, k, N, and F must be true?

Mechanics II. Which of the following relations among the forces W, k, N, and F must be true? Mechanics II 1. By applying a force F on a block, a person pulls a block along a rough surface at constant velocity v (see Figure below; directions, but not necessarily magnitudes, are indicated). Which

More information

spring magnet Fig. 7.1 One end of the magnet hangs inside a coil of wire. The coil is connected in series with a resistor R.

spring magnet Fig. 7.1 One end of the magnet hangs inside a coil of wire. The coil is connected in series with a resistor R. 1 A magnet is suspended vertically from a fixed point by means of a spring, as shown in Fig. 7.1. spring magnet coil R Fig. 7.1 One end of the magnet hangs inside a coil of wire. The coil is connected

More information

1 Fig. 3.1 shows the variation of the magnetic flux linkage with time t for a small generator. magnetic. flux linkage / Wb-turns 1.

1 Fig. 3.1 shows the variation of the magnetic flux linkage with time t for a small generator. magnetic. flux linkage / Wb-turns 1. 1 Fig. 3.1 shows the variation of the magnetic flux linkage with time t for a small generator. 2 magnetic 1 flux linkage / 0 10 2 Wb-turns 1 2 5 10 15 t / 10 3 s Fig. 3.1 The generator has a flat coil

More information

Regulating Web Tension in Tape Systems with Time-varying Radii

Regulating Web Tension in Tape Systems with Time-varying Radii Regulating Web Tension in Tape Systems with Time-varying Radii Hua Zhong and Lucy Y. Pao Abstract A tape system is time-varying as tape winds from one reel to the other. The variations in reel radii consist

More information

The Simplest Tether Control Law in a Small Satellite

The Simplest Tether Control Law in a Small Satellite The Simplest Tether Control Law in a Small Satellite Yosuke Nakamura (Kyushu University) Designing of a small satellite based on the Satellite Design Contest is progressing in Kyushu University. This is

More information

Mathematical Modeling and response analysis of mechanical systems are the subjects of this chapter.

Mathematical Modeling and response analysis of mechanical systems are the subjects of this chapter. Chapter 3 Mechanical Systems A. Bazoune 3.1 INRODUCION Mathematical Modeling and response analysis of mechanical systems are the subjects of this chapter. 3. MECHANICAL ELEMENS Any mechanical system consists

More information

1. A sphere with a radius of 1.7 cm has a volume of: A) m 3 B) m 3 C) m 3 D) 0.11 m 3 E) 21 m 3

1. A sphere with a radius of 1.7 cm has a volume of: A) m 3 B) m 3 C) m 3 D) 0.11 m 3 E) 21 m 3 1. A sphere with a radius of 1.7 cm has a volume of: A) 2.1 10 5 m 3 B) 9.1 10 4 m 3 C) 3.6 10 3 m 3 D) 0.11 m 3 E) 21 m 3 2. A 25-N crate slides down a frictionless incline that is 25 above the horizontal.

More information

The Simple Harmonic Oscillator

The Simple Harmonic Oscillator The Simple Harmonic Oscillator Michael Fowler, University of Virginia Einstein s Solution of the Specific Heat Puzzle The simple harmonic oscillator, a nonrelativistic particle in a potential ½C, is a

More information

PHYSICS 221, FALL 2009 EXAM #1 SOLUTIONS WEDNESDAY, SEPTEMBER 30, 2009

PHYSICS 221, FALL 2009 EXAM #1 SOLUTIONS WEDNESDAY, SEPTEMBER 30, 2009 PHYSICS 221, FALL 2009 EXAM #1 SOLUTIONS WEDNESDAY, SEPTEMBER 30, 2009 Note: The unit vectors in the +x, +y, and +z directions of a right-handed Cartesian coordinate system are î, ĵ, and ˆk, respectively.

More information

Rigid Body Dynamics, SG2150 Solutions to Exam,

Rigid Body Dynamics, SG2150 Solutions to Exam, KTH Mechanics 011 10 Calculational problems Rigid Body Dynamics, SG150 Solutions to Eam, 011 10 Problem 1: A slender homogeneous rod of mass m and length a can rotate in a vertical plane about a fied smooth

More information

ECE 3318 Applied Electricity and Magnetism Spring 2018 Homework #7

ECE 3318 Applied Electricity and Magnetism Spring 2018 Homework #7 EE 3318 Applied Electricity and Magnetism Spring 2018 Homework #7 Date assigned: Tuesday, March 6, 2018 Date due: Tuesday, March 20, 2018 Do Probs. 1, 2, and 7-12. (You are welcome to do the other problems

More information

National Quali cations

National Quali cations National Quali cations AH08 X70/77/ Mathematics of Mechanics TUESDAY, 9 MAY :00 PM :00 PM Total marks 00 Attempt ALL questions. You may use a calculator. Full credit will be given only to solutions which

More information

Positioning Servo Design Example

Positioning Servo Design Example Positioning Servo Design Example 1 Goal. The goal in this design example is to design a control system that will be used in a pick-and-place robot to move the link of a robot between two positions. Usually

More information

Kinematic and Dynamic Models

Kinematic and Dynamic Models Kinematic and Dynamic Models b Advanced Robotics: Navigation and Perception 3/29/2011 1 Reasoning about mobility is important. Image: RedTeamRacing.com 2 ecture Outline Nomenclature Kinematic Constraints

More information

Robotics. Dynamics. Marc Toussaint U Stuttgart

Robotics. Dynamics. Marc Toussaint U Stuttgart Robotics Dynamics 1D point mass, damping & oscillation, PID, dynamics of mechanical systems, Euler-Lagrange equation, Newton-Euler recursion, general robot dynamics, joint space control, reference trajectory

More information

Evaluation of active structural vibration control strategies in milling process

Evaluation of active structural vibration control strategies in milling process Evaluation of active structural vibration control strategies in milling process Monnin, J. (a); Wegener, K. (a) a) Institute of Machine Tools and Manufacturing, Zurich, Switzerland Keywords: Mechatronics,

More information

Exact and Approximate Numbers:

Exact and Approximate Numbers: Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.

More information

6.7 Variation and Problem Solving. OBJECTIVES 1 Solve Problems Involving Direct Variation. 2 Solve Problems Involving Inverse Variation.

6.7 Variation and Problem Solving. OBJECTIVES 1 Solve Problems Involving Direct Variation. 2 Solve Problems Involving Inverse Variation. 390 CHAPTER 6 Rational Epressions 66. A doctor recorded a body-mass inde of 7 on a patient s chart. Later, a nurse notices that the doctor recorded the patient s weight as 0 pounds but neglected to record

More information

HW7: Ch. 26 P 34, 36 Ch.27 Q 2, 4, 8, 18 P 2, 8, 17, 19, 37

HW7: Ch. 26 P 34, 36 Ch.27 Q 2, 4, 8, 18 P 2, 8, 17, 19, 37 Fall 12 PHY 122 Homework Solutions #7 HW7: Ch. 26 P 34, 36 Ch.27 Q 2, 4, 8, 18 P 2, 8, 17, 19, 37 Chapter 26 Problem 34 Determine the magnitudes and directions of the currents in each resistor shown in

More information

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture 6-5 November, 2015

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture 6-5 November, 2015 The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture 6-5 November, 015 Institute of Structural Engineering Method of Finite Elements II 1 Introduction

More information

Lecture 4.2 Finite Difference Approximation

Lecture 4.2 Finite Difference Approximation Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by

More information

Case Study: The Pelican Prototype Robot

Case Study: The Pelican Prototype Robot 5 Case Study: The Pelican Prototype Robot The purpose of this chapter is twofold: first, to present in detail the model of the experimental robot arm of the Robotics lab. from the CICESE Research Center,

More information

0.2 F 0.4 F 0.8 F 4.0 F

0.2 F 0.4 F 0.8 F 4.0 F 1 When a β particle moves at right angles through a uniform magnetic field it experiences a force F. n α particle moves at right angles through a magnetic field of twice the magnetic flux density with

More information

Overview of motors and motion control

Overview of motors and motion control Overview of motors and motion control. Elements of a motion-control system Power upply High-level controller ow-level controller Driver Motor. Types of motors discussed here; Brushed, PM DC Motors Cheap,

More information

CONTROL OF THE NONHOLONOMIC INTEGRATOR

CONTROL OF THE NONHOLONOMIC INTEGRATOR June 6, 25 CONTROL OF THE NONHOLONOMIC INTEGRATOR R. N. Banavar (Work done with V. Sankaranarayanan) Systems & Control Engg. Indian Institute of Technology, Bombay Mumbai -INDIA. banavar@iitb.ac.in Outline

More information

Physics I (Navitas) EXAM #2 Spring 2015

Physics I (Navitas) EXAM #2 Spring 2015 95.141 Physics I (Navitas) EXAM #2 Spring 2015 Name, Last Name First Name Student Identification Number: Write your name at the top of each page in the space provided. Answer all questions, beginning each

More information

4.9 Free Mechanical Vibrations

4.9 Free Mechanical Vibrations 4.9 Free Mechanical Vibrations Spring-Mass Oscillator When the spring is not stretched and the mass m is at rest, the system is at equilibrium. Forces Acting in the System When the mass m is displaced

More information

Phys 111 Exam 1 September 22, 2015

Phys 111 Exam 1 September 22, 2015 Phys 111 Exam 1 September 22, 2015 1. The time T required for one complete oscillation of a mass m on a spring of force constant k is T = 2π m k. Find the dimension of k to be dimensionally correct for

More information

Control for. Maarten Steinbuch Dept. Mechanical Engineering Control Systems Technology Group TU/e

Control for. Maarten Steinbuch Dept. Mechanical Engineering Control Systems Technology Group TU/e Control for Maarten Steinbuch Dept. Mechanical Engineering Control Systems Technology Group TU/e Motion Systems m F Introduction Timedomain tuning Frequency domain & stability Filters Feedforward Servo-oriented

More information

1. An object is dropped from rest. Which of the five following graphs correctly represents its motion? The positive direction is taken to be downward.

1. An object is dropped from rest. Which of the five following graphs correctly represents its motion? The positive direction is taken to be downward. Unless otherwise instructed, use g = 9.8 m/s 2 Rotational Inertia about an axis through com: Hoop about axis(radius=r, mass=m) : MR 2 Hoop about diameter (radius=r, mass=m): 1/2MR 2 Disk/solid cyllinder

More information

Adaptive Robust Control (ARC) for an Altitude Control of a Quadrotor Type UAV Carrying an Unknown Payloads

Adaptive Robust Control (ARC) for an Altitude Control of a Quadrotor Type UAV Carrying an Unknown Payloads 2 th International Conference on Control, Automation and Systems Oct. 26-29, 2 in KINTEX, Gyeonggi-do, Korea Adaptive Robust Control (ARC) for an Altitude Control of a Quadrotor Type UAV Carrying an Unknown

More information

Jerk derivative feedforward control for motion systems

Jerk derivative feedforward control for motion systems Jerk derivative feedforward control for motion systems Matthijs Boerlage Rob Tousain Maarten Steinbuch Abstract This work discusses reference trajectory relevant model based feedforward design. For motion

More information

Wataru Ohnishi a) Student Member, Hiroshi Fujimoto Senior Member Koichi Sakata Member, Kazuhiro Suzuki Non-member Kazuaki Saiki Member.

Wataru Ohnishi a) Student Member, Hiroshi Fujimoto Senior Member Koichi Sakata Member, Kazuhiro Suzuki Non-member Kazuaki Saiki Member. IEEJ International Workshop on Sensing, Actuation, and Motion Control Proposal of Decoupling Control Method for High-Precision Stages using Multiple Actuators considering the Misalignment among the Actuation

More information

The problem of singularity in impacting systems

The problem of singularity in impacting systems The problem of singularity in impacting systems Soumitro Banerjee, Department of Physics Indian Institute of Science Education & Research, Kolkata, India The problem of singularity in impacting systems

More information

Linear and Nonlinear Oscillators (Lecture 2)

Linear and Nonlinear Oscillators (Lecture 2) Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical

More information

Physics 2101, First Exam, Spring 2007

Physics 2101, First Exam, Spring 2007 Physics 2101, First Exam, Spring 2007 January 25, 2007 Name : Section: (Circle one) 1 (Rupnik, MWF 7:40am) 2 (Giammanco, MWF 9:40am) 3 (Rupnik, MWF 11:40am) 4 (Rupnik, MWF 2:40pm) 5 (Giammanco, TTh 10:40am)

More information

INTRODUCTION. Description

INTRODUCTION. Description INTRODUCTION "Acoustic metamaterial" is a label that encompasses for acoustic structures that exhibit acoustic properties not readily available in nature. These properties can be a negative mass density,

More information

Robotics. Dynamics. University of Stuttgart Winter 2018/19

Robotics. Dynamics. University of Stuttgart Winter 2018/19 Robotics Dynamics 1D point mass, damping & oscillation, PID, dynamics of mechanical systems, Euler-Lagrange equation, Newton-Euler, joint space control, reference trajectory following, optimal operational

More information

Question 1. Question 2. Question 3

Question 1. Question 2. Question 3 Question 1 Switch S in in the figure is closed at time t = 0, to begin charging an initially uncharged capacitor of capacitance C = 18.2 μf through a resistor of resistance R = 22.3 Ω. At what time (in

More information

Equal Pitch and Unequal Pitch:

Equal Pitch and Unequal Pitch: Equal Pitch and Unequal Pitch: Equal-Pitch Multiple-Stack Stepper: For each rotor stack, there is a toothed stator segment around it, whose pitch angle is identical to that of the rotor (θs = θr). A stator

More information

25.2. Applications of PDEs. Introduction. Prerequisites. Learning Outcomes

25.2. Applications of PDEs. Introduction. Prerequisites. Learning Outcomes Applications of PDEs 25.2 Introduction In this Section we discuss briefly some of the most important PDEs that arise in various branches of science and engineering. We shall see that some equations can

More information

Fig Use Fig. 3.1 to state the physical properties of this metal. In your answer, you should use appropriate technical terms, spelled correctly.

Fig Use Fig. 3.1 to state the physical properties of this metal. In your answer, you should use appropriate technical terms, spelled correctly. 1 (a) Fig. 3.1 shows the stress against strain graph for a metal X up to its breaking point. stress X 0 0 Fig. 3.1 strain Use Fig. 3.1 to state the physical properties of this metal. In your answer, you

More information

Nonlinear effects on the rotor driven by a motor with limited power

Nonlinear effects on the rotor driven by a motor with limited power Applied and Computational Mechanics 1 (007) 603-61 Nonlinear effects on the rotor driven by a motor with limited power L. Pst Institute of Thermomechanics, Academy of Sciences of CR, Dolejškova 5,18 00

More information

Problem Set 6: Magnetism

Problem Set 6: Magnetism University of Alabama Department of Physics and Astronomy PH 10- / LeClair Spring 008 Problem Set 6: Magnetism 1. 10 points. A wire with a weight per unit length of 0.10 N/m is suspended directly above

More information

Interaction-based Dynamic Measurement of Haptic Characteristics of Control Elements

Interaction-based Dynamic Measurement of Haptic Characteristics of Control Elements Interaction-based Dynamic Measurement of Haptic Characteristics of Control Elements Wenliang Zhou 1, Jörg Reisinger 1, Angelika Peer 2, and Sandra Hirche 3 1 Department HMI Components, Daimler AG, Sindelfingen,

More information

Phys 111 Exam 1 September 19, You cannot use CELL PHONES, ipad, IPOD... Good Luck!!! Name Section University ID

Phys 111 Exam 1 September 19, You cannot use CELL PHONES, ipad, IPOD... Good Luck!!! Name Section University ID Phys 111 Exam 1 September 19, 2017 Name Section University ID Please fill in your computer answer sheet as follows: 1) In the NAME grid, fill in your last name, leave one blank space, then your first name.

More information

Reduction of the effect of floor vibrations in a checkweigher using an electromagnetic force balance system

Reduction of the effect of floor vibrations in a checkweigher using an electromagnetic force balance system ACTA IMEKO ISSN: 1 870X July 017, Volume 6 Number, 65 69 Reduction of the effect of floor vibrations in a checkweigher using an electromagnetic force balance system Yuji Yamakawa 1, Takanori Yamazaki 1

More information

Design of a High Speed and High Precision 3 DOF Linear Direct Drive

Design of a High Speed and High Precision 3 DOF Linear Direct Drive Design of a High Speed and High Precision 3 DOF Linear Direct Drive Bernhard Sprenger Swiss Federal Institute of Technology Zurich (ETHZ) Institute of Robotics 8092 Zurich, Switzerland Email: sprenger@ifr.mavt.ethz.ch

More information

Predictive Cascade Control of DC Motor

Predictive Cascade Control of DC Motor Volume 49, Number, 008 89 Predictive Cascade Control of DC Motor Alexandru MORAR Abstract: The paper deals with the predictive cascade control of an electrical drive intended for positioning applications.

More information

AP Physics C Mechanics Objectives

AP Physics C Mechanics Objectives AP Physics C Mechanics Objectives I. KINEMATICS A. Motion in One Dimension 1. The relationships among position, velocity and acceleration a. Given a graph of position vs. time, identify or sketch a graph

More information

Derivation and Application of a Conserved Orbital Energy for the Inverted Pendulum Bipedal Walking Model

Derivation and Application of a Conserved Orbital Energy for the Inverted Pendulum Bipedal Walking Model Derivation and Application of a Conserved Orbital Energy for the Inverted Pendulum Bipedal Walking Model Jerry E. Pratt and Sergey V. Drakunov Abstract We present an analysis of a point mass, point foot,

More information

Chapter 12. Magnetism and Electromagnetism

Chapter 12. Magnetism and Electromagnetism Chapter 12 Magnetism and Electromagnetism 167 168 AP Physics Multiple Choice Practice Magnetism and Electromagnetism SECTION A Magnetostatics 1. Four infinitely long wires are arranged as shown in the

More information

Chapter Units and Measurement

Chapter Units and Measurement 2 Chapter Units and Measurement 1. Identify the pair whose dimensions are equal [2002] torque and work stress and energy force and stress force and work 2. [2003] [L -1 T] ] [L -2 T 2 ] [L 2 T -2 ] [LT

More information

λ w w. (3) w 2 2 λ The obtained equation (4) may be presented in the stationary approximation as follows (5) with the boundary conditions: for = 0

λ w w. (3) w 2 2 λ The obtained equation (4) may be presented in the stationary approximation as follows (5) with the boundary conditions: for = 0 International Conference on Methods of Aerophysical Research, ICMAR 28 SIMULATION OF POROUS MATERIAL IMPREGNATION A.V. Fedorov, A.A. Zhilin Khristianovich Institute of Theoretical and Applied Mechanics

More information

Automatic Control Motion planning

Automatic Control Motion planning Automatic Control Motion planning (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Motivations 2 Electric motors are used in many different applications,

More information

Intelligent Control of a SPM System Design with Parameter Variations

Intelligent Control of a SPM System Design with Parameter Variations Intelligent Control of a SPM System Design with Parameter Variations Jium-Ming Lin and Po-Kuang Chang Abstract This research is to use fuzzy controller in the outer-loop to reduce the hysteresis effect

More information

State Space and Hidden Markov Models

State Space and Hidden Markov Models State Space and Hidden Markov Models Kunsch H.R. State Space and Hidden Markov Models. ETH- Zurich Zurich; Aliaksandr Hubin Oslo 2014 Contents 1. Introduction 2. Markov Chains 3. Hidden Markov and State

More information

Today s goals So far Today 2.004

Today s goals So far Today 2.004 Today s goals So far Feedback as a means for specifying the dynamic response of a system Root Locus: from the open-loop poles/zeros to the closed-loop poles Moving the closed-loop poles around Today Proportional

More information

UNIVERSITY OF CAMBRIDGE INTERNATIONAL EXAMINATIONS General Certificate of Education Advanced Level

UNIVERSITY OF CAMBRIDGE INTERNATIONAL EXAMINATIONS General Certificate of Education Advanced Level UNIVERSITY OF CAMBRIDGE INTERNATIONAL EXAMINATIONS General Certificate of Education Advanced Level *2082050702* PHYSICS 9702/04 Paper 4 A2 Structured Questions May/June 2007 1 hour 45 minutes Candidates

More information

Physics 23 Exam 3 April 2, 2009

Physics 23 Exam 3 April 2, 2009 1. A string is tied to a doorknob 0.79 m from the hinge as shown in the figure. At the instant shown, the force applied to the string is 5.0 N. What is the torque on the door? A) 3.3 N m B) 2.2 N m C)

More information

a. What is the angular frequency ω of the block in terms of k, l, and m?

a. What is the angular frequency ω of the block in terms of k, l, and m? 1 Problem 1: (4 pts.) Two spherical planets, each of mass M and Radius R, start out at rest with a distance from center to center of 4R. What is the speed of one of the planets at the moment that their

More information

2. a) Explain the equilibrium of i) Concurrent force system, and ii) General force system.

2. a) Explain the equilibrium of i) Concurrent force system, and ii) General force system. Code No: R21031 R10 SET - 1 II B. Tech I Semester Supplementary Examinations Dec 2013 ENGINEERING MECHANICS (Com to ME, AE, AME, MM) Time: 3 hours Max. Marks: 75 Answer any FIVE Questions All Questions

More information

Basic Mathematics and Units

Basic Mathematics and Units Basic Mathematics and Units Rende Steerenberg BE/OP Contents Vectors & Matrices Differential Equations Some Units we use 3 Vectors & Matrices Differential Equations Some Units we use 4 Scalars & Vectors

More information

Theory of Vibrations in Stewart Platforms

Theory of Vibrations in Stewart Platforms Theory of Vibrations in Stewart Platforms J.M. Selig and X. Ding School of Computing, Info. Sys. & Maths. South Bank University London SE1 0AA, U.K. (seligjm@sbu.ac.uk) Abstract This article develops a

More information

Dynamics Kinetics of a particle Section 4: TJW Force-mass-acceleration: Example 1

Dynamics Kinetics of a particle Section 4: TJW Force-mass-acceleration: Example 1 Section 4: TJW Force-mass-acceleration: Example 1 The beam and attached hoisting mechanism have a combined mass of 1200 kg with center of mass at G. If the inertial acceleration a of a point P on the hoisting

More information

Forced Mechanical Vibrations

Forced Mechanical Vibrations Forced Mechanical Vibrations Today we use methods for solving nonhomogeneous second order linear differential equations to study the behavior of mechanical systems.. Forcing: Transient and Steady State

More information