Error Analysis and Correction of a High Precision Co-ordinate Table

Size: px
Start display at page:

Download "Error Analysis and Correction of a High Precision Co-ordinate Table"

Transcription

1 Error Analysis an Correction of a High Precision Co-orinate Table Gyula Hermann BrainWare Lt. Völgy utca 13/A, H-1021 Buapest Hungary hermgy@iif.hu Abstract: A small scale, high precision co-orinate table was evelope as part of a project to solve measurement problems originating from microtechnic. Construction with a single reference surface was selecte an vacuum preloae air bearing applie in orer to enable smooth motion an eliminate stick slip. Laser interferometers were applie to measure both the isplacement an the angular error. The paper escribes the mathematical evaluation of the measuring system s error an how that is taken into account in the evaluation software. Keywors: Co-orinate table, measuring system, evaluation software, error computation 1 Introuction Co-orinate metrology is toay a firmly establishe technique in inustry. The universal applicability an high egree of automation accounts for it s succes. Two other factors are equally important: tolerances are becomming closer as a result of interchangeable manufacturing an the increase of subcontracting is obliging manufacturers to assume responsibility for the imensional integrity of their parts. With the appearance of micromachine mechanical elements increase eman can be observe for higher resolution an more accurate inspection of these parts. 2 The Table Construction For many reasons ultra precise X-Y tables are guie by a flat reference surface. These tables usually employ vacuum or magnetically preloae air bearings. Their avantage is that only one flat guie surface is necessary, whereas the oppose

2 bearing preloaing requires two flat guie surfaces that are parallel. Other avantages of the air bearings are the zero static friction, which makes infinite resolution an the very high repeatability possible. Moreover porous air bearings average the errors of the guie surface finish an irregularities. In our case, the table has four legs an rests on four vacuum preloae flat porous air bearings. The porous carbon in the bearing prouces a uniform air pressure across the face of the bearing while automatically restricting an amping the air flow at the same time. As a consequence the possibility of collapse is reuce an a higher pitch moment stiffness is reache. As a reference surface a 000 class granite[1] control surface plate, having an overall flatness of less then 2 µm, is use an the axes are riven by linear motors. The winow in the moving part of the table enables the installation of a back lighting LED panel. Guieway Linear motors Winow Reference surface Preloae airbearings Carriage Figure 1 The mechanical construction

3 3 The Measuring System The table position is etermine by two plane mirror laser interferometers, a single- an a ualbeam, relative to the measuring frame fixe to the granite base. The laser sources are couple to the sensor heas via fiberoptic cables. The single beam interferometer enables a isplacement measurement along one axis with a resolution of 1 nm (uner ieal environmental conitions) while the ualbeam interferometer allows the simultaneous isplacement an angular measurement. The angular resolution is approximately 0.05 arcsec. The measuring range is 100 mm an ±2 arcmin respectively. singlebeam laser interferometer ualbeam laser interferometer electronic unit X off Y off fiberoptic link the moving part of the table Figure 2 The outline of the measuring system

4 4 The Geometric Moel of the Measuring System Y ref Y table T size T size L x1 Y off L x2 X c,y c + (p x,p y ).θ X table L y 0,0 X off X ref Figure 3 The geometric moel Where: X ref an Y ref are X the Y axis of the measuring frame co-orinate system X table an Y table are X the Y axis of the table co-orinate system X c an Y c are the co-orinates of the table coorinate system s origo X off an Y off are the positions of the interferometers along the measuring frame axis is the separation between the two beams in the ualbeam interferometer L x1, L x2 an L y are the isplacements measure by the interferometer T size is the table size Let us in the first instance assume that the roll an pitch error of the table is negligable. From Fig. 3 the following two equation for the origo of the table coorinate system axis can be seen: X c -L x2 Y off -Y c = tgθ (1)

5 L y - Y c X off - X c = tgθ (2) Where tgθ = L x1-l x2 (3) Hereout follows for the co-orinates of the origo: X c = ( X off - L x2 )tg 2 θ 1+tg 2 θ Y c = Y off.tg 2 θ 1+tg 2 θ + tgθ(l x2 - X off ) Ly 1+tg 2 θ + 1+tg 2 θ + (Y off -L y )tgθ 1+tg 2 θ + L x2 (4) The co-orinates of an arbitrary point on the table can be expresse by using the following equation: P x P y = cos -sin sin cos p x p + y X c Y ===> c cos sin -sin cos P x -X c P y -Y = c p x p (6) y (5) 5 The Error Moel In the above outline measuring system we have to take into account the following error sources: isplacement measurement error as a result of wavelength variation angle measurement error perpenicularity error of the laserbeams an the reflecting mirrors flatness error of the mirrors reference surface flatness error resulting in role an pitch error on the table Let us consier these error components one by one: The isplacement of the reflecting mirror is given by the following equation: D x = Where: c 2f (N + φ 2π ) (7) c is the spee of light f is the laser frequency

6 N is the number 2πs an Φ is the phase angle The laser frequency stability, after warmup is better than 3 x The spee of light is influence by the temperature, the humiity, the turbulence an the contamination of the air. The frequency of the laser can be stabilize an the effect of temperature an humiity can be compensate igitaly by measuring both air an object temperature an the humiity. The error equation for the angle is: Δtgθ Δ(L x2-l x1 ) + Δ tgθ (8) Hereout follows that the angular error is etermine by the relative error of the beamistance of the ualbeam interferometer. The eviation from the prepenicularity of the reflecting mirrors an that of the irections of the interferometers contribute to the overall error of the system. Also the perpenicularity error of the reference system (X ref,y ref ) has a similar effect. The angle between the mirrors is φ, whereas the angles between the reference axis an the axis of the interferometers are given by α an β respectively. From Fig. 4 it can be immeiatly seen that: Δ = 90 - φ = θ - ψ (9) Where Δ = 90º- φ is the perpenicularity error of the table (the reflecting mirrors)

7 YBrefB TBsizeB Y table T size L x1 L x2.φ X table.θ Y off.α X c,y c.ψ ß 0,0 X off L y X ref Figure 4 Perpenicularity error moel X c - L x2 cosα Y off - L x2 sinα - Y c = tg ψ (10) L y cosβ - Y c X off - L y sinβ - X c = tg θ (11) (L x2 - L x1 ) cosα + (L x2 - L x1 ) sinα = tg ψ (12) From hereout follows by substituting sinα = α, cosα = 1, sin β = β, cos β = 1, an neglecting the higher orer small terms we get: tg θ meas = L x2-l x1 an tg θ = tg(90º-φ) + tgψ 1 -tg(90º-φ)tgψ (L x2 -L x1 )cosα θ = 90 - ( φ + ψ) = Δ - arctg + (L x2 -L x1 )sinα = Δ - arctg tgθ meascosα 1+tgθ meas sinα (13) L x2 X c = 1+αtgθ + tgθ meas (Δ+tgθ meas )(X off -L x2 ) meas 1+(3α+Δ)tgθ meas +tg 2 θ meas + - (1+2αtgθ meas)tgθ meas L y 1+(3α+Δ)tgθ meas +tg 2 θ meas + αtg 3 θ meas X off 1+(3α+Δ)tgθ meas +tg 2 θ meas - (1+(2α+Δ)tgθ meas-δtgθ meas )tgθ meas Y off 1+(3α+Δ)tgθ meas +tg 2 θ meas (14)

8 (1+2αtgθ meas ) L y Y c = ( 1+(2α+Δ)tgθ meas +tg 2 θ meas ) + (Δ+tgθ meas )(L x2 - X off ) (1+(2α+Δ)tgθ meas +tg 2 θ meas ) + αtg 2 θ meas X off + ( 1+(2α+Δ)tgθ meas +tg 2 θ meas ) + tgθ meas (Δ+tgθ meas ) Y off (1+(2α+Δ)tgθ meas +tg 2 θ meas ) Using these results the co-orinates of a point in the table co-orinate system can be exprese using equation [6]. The geometric constant values Δ, α an β use in the computation can be measure by conventional optical techniques. If we consier the laserbeam as ray with nearly zero iameter then the flatness error of the mirrors results in an error of the isplacement measurement. However the finite iameter of the beam has an integration property, so the local flatness error can be compensate to a certain exten. Waveness errors larger then the beam iameter have to be taken into account. A map of the mirror surface is the necessary an sufficient information for the compensation of this type of errors. The map can be capture by interferometric techniques. The applie technologies in manufacturing optical surfaces ensure that the functions escribing the surface is continoues into the secon erivatives. Therefore a spline moel of the surface is space an time saving. (15) Y ref Y table Best fit straightline L x1 Angle between the best fit straightlines L x2 X table Y off X c,y c Best fit straightline L y 0,0 X off X ref Figure 5 Moel of the straightness (flatness) error

9 6 The Software Implementation The interferometer elivers through USB ports two current co-orinate an a current angle value. These signals are connecte to a personal computer responsible for the calculations. As high presion measurements are connecte with relavtively slow isplacements the computational time is not a critical factor an as we have seen the amount of ata involve is also small. The compensation for the spee of light eviation is performe by the electronic unit of the interferometer. There is no nee for any aitional harware. The program performing the above mentione calculations is written in C++ an runs on the Winow XP operating system. Base on an analysis of the formulas at some point multilenght arithmetic is use in orer to ensure the require accuracy. The measurement result wil be publishe in separate paper. Conclusions The research a part of which presente here prove the usefulness of software error compensation in co-orinate measurement. Accoring to our expectation combine with new construction the accuracy can be improve with an orer of magnitue. The initial results are promising an the technique can extene into the thir imension. Acknowlegment This research was supporte by the Hungarian Scientific Research Foun through the contract T References [1] G. Belforte, B. Bona, E. Canuto, F. Donati, F. Ferraris, I. Gorini, S. Morei, M. Peisino, S. Sartori Coorinate Measuring Machines an Machine Tools Selfcalibration an Error Correction Annals of the CIRP Vol. 36/1/1987, pp [2] B. P. Demiovich, I. A. Maron Computational Mathematics MIR Publishers Moscow 1973 [3] K. F. Eman, B. T. Wu, A Generalize Geometric Error Moel for Multi- Axis Machines Annals of the CIRP Vol. 36/1/1987, pp [4] P. M. Ferreire, C. R. Liu An Analitical Quaratic Moel for the Geometric Error of a Machine Tool Journal of Manufacturing Systems Vol. 5, No. 1, pp [5] B. Meulenbel, Analyse Uitgeverij Stam N. V. Haarlem 1962 [6] F. C. e Nijs Moelling of a Coorinate Measuring Machine for Analysis of Its Dynamic Behaviour Annals of the CIRP Vol. 37/2/1988, pp

10 [7] J. W. M. C. Teeuiwsen, J. A. Soons, P. H. J Schellekens A General Methoe for Error Description of CMMs Using Polinomial Fitting Proceures Annals of the CIRP Vol. 38/1/1989, pp [8] M. M. Weck, Geometric an kinematic errors. Technology of Machine Tools, 1980, 5, 9-12 [9] M. M. Weck, Werkzeugmaschinen Ban 2 VDI-Verlag Düsselorf 1981 [10] G. Zhang, R. Veale, B. Borchart an R. Hocken Error Compensation of Coorinate Measuring Machines Annals of the CIRP Vol. 34/1/1985, pp

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

sensors ISSN by MDPI

sensors ISSN by MDPI Sensors 006, 6, 546-556 Special Issue Intelligent Sensors Eite by Wilmar Hernanez Full Research Paper sensors ISSN 144-80 006 by MDPI http://www.mpi.org/sensors Differential Laser Doppler base Non-Contact

More information

International Conference KNOWLEDGE-BASED ORGANIZATION Vol. XXIII No

International Conference KNOWLEDGE-BASED ORGANIZATION Vol. XXIII No International Conference KNOWLEDGE-BAED ORGANIZATION Vol. XXIII No 3 2017 METHOD FOR DETERMINATION OF THE PARAMETER OF A MACHINE GUN UPENION MOUNTED ON AN ARMOURED VEHICLE vilen PIRDONOV, vilen TEFANOV

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

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

Solutions to the Exercises of Chapter 9

Solutions to the Exercises of Chapter 9 9A. Vectors an Forces Solutions to the Exercises of Chapter 9. F = 5 sin 5.9 an F = 5 cos 5 4.8.. a. By the Pythagorean theorem, the length of the vector from to (, ) is + = 5. So the magnitue of the force

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

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

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

State-Space Model for a Multi-Machine System

State-Space Model for a Multi-Machine System State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal

More information

Vibration Analysis of Railway Tracks Forced by Distributed Moving Loads

Vibration Analysis of Railway Tracks Forced by Distributed Moving Loads IJR International Journal of Railway Vol. 6, No. 4 / December 13, pp. 155-159 The Korean Society for Railway Vibration Analysis of Railway Tracks Force by Distribute Moving Loas Sinyeob Lee*, Dongkyu Kim*,

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

General Data. Types of bearings 6. Standardization and interchangeability 12. Dimensions and part numbers 14. Bearing manufacturing precision 18

General Data. Types of bearings 6. Standardization and interchangeability 12. Dimensions and part numbers 14. Bearing manufacturing precision 18 General ata Types of bearings 6 efinitions 6 Vocabulary 8 Capabilities 9 Stanarization an interchangeability 12 The Stanars 12 Interchangeability 12 imensions an part numbers 14 General esignations 14

More information

EXPERIMENTAL INVESTIGATION ON PNEUMATIC COMPONENTS

EXPERIMENTAL INVESTIGATION ON PNEUMATIC COMPONENTS Conference on Moelling Flui Flow (CMFF 03) The 12 th International Conference on Flui Flow Technologies Buapest, Hungary, September 3-6, 2003 EXPERIMENTAL INVESTIGATION ON PNEUMATIC COMPONENTS Zoltán MÓZER,

More information

ELECTRON DIFFRACTION

ELECTRON DIFFRACTION ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the

More information

Tutorial Test 5 2D welding robot

Tutorial Test 5 2D welding robot Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.

More 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

Adaptive Back-Stepping Control of Automotive Electronic Control Throttle

Adaptive Back-Stepping Control of Automotive Electronic Control Throttle Journal of Software Engineering an Applications, 07, 0, 4-55 http://www.scirp.org/journal/jsea ISSN Online: 945-34 ISSN Print: 945-36 Aaptive Back-Stepping Control of Automotive Electronic Control Throttle

More information

A New Approach in Analytical Analysis of Eddy Currents in Laminated Core

A New Approach in Analytical Analysis of Eddy Currents in Laminated Core J. Basic. Appl. Sci. Res., (7)741-745, 1 1, TextRoa Publication ISSN 9-434 Journal of Basic an Applie Scientific Research www.textroa.com A New Approach in Analtical Analsis of E Currents in Laminate Core

More information

Sensors & Transducers 2015 by IFSA Publishing, S. L.

Sensors & Transducers 2015 by IFSA Publishing, S. L. Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity

More information

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A AN INTRODUCTION TO AIRCRAFT WIN FLUTTER Revision A By Tom Irvine Email: tomirvine@aol.com January 8, 000 Introuction Certain aircraft wings have experience violent oscillations uring high spee flight.

More information

Pure Further Mathematics 1. Revision Notes

Pure Further Mathematics 1. Revision Notes Pure Further Mathematics Revision Notes June 20 2 FP JUNE 20 SDB Further Pure Complex Numbers... 3 Definitions an arithmetical operations... 3 Complex conjugate... 3 Properties... 3 Complex number plane,

More information

Machine Positioning Uncertainty with Laser Interferometer Feedback

Machine Positioning Uncertainty with Laser Interferometer Feedback Machine Positioning Uncertainty with Laser Interferometer Feedback The purpose of this discussion is to explain the major contributors to machine positioning uncertainty in systems with laser interferometer

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

An Approach for Design of Multi-element USBL Systems

An Approach for Design of Multi-element USBL Systems An Approach for Design of Multi-element USBL Systems MIKHAIL ARKHIPOV Department of Postgrauate Stuies Technological University of the Mixteca Carretera a Acatlima Km. 2.5 Huajuapan e Leon Oaxaca 69000

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

Evaporating droplets tracking by holographic high speed video in turbulent flow

Evaporating droplets tracking by holographic high speed video in turbulent flow Evaporating roplets tracking by holographic high spee vieo in turbulent flow Loïc Méès 1*, Thibaut Tronchin 1, Nathalie Grosjean 1, Jean-Louis Marié 1 an Corinne Fournier 1: Laboratoire e Mécanique es

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

Introduction Basic principles Finite element formulation Nonlinear algorithms Validation & examples Oofelie::MEMS, driven by SAMCEF Field Perspectives

Introduction Basic principles Finite element formulation Nonlinear algorithms Validation & examples Oofelie::MEMS, driven by SAMCEF Field Perspectives Non linear behavior of electrostatically actuate micro-structures Dr. Ir. Stéphane Paquay, Open Engineering SA Dr. Ir. Véronique Rochus, ULg (LTAS-VIS) Dr. Ir. Stefanie Gutschmit, ULg (LTAS-VIS) Outline

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

Crossed Roller Bearings

Crossed Roller Bearings Special Selection & http://www.ikont.eu Official Distributer in Unite Kingom rosse Bearings Unit rackley Way Peartree Lane Duley, West Milans DY UW Tel : 9 / E-mail : uley@goiva-bearings.co.uk for AT-EG

More information

Modeling and analysis of parallel connected permanent magnet synchronous generators in a small hydropower plant

Modeling and analysis of parallel connected permanent magnet synchronous generators in a small hydropower plant Proceeings of the 2006 IASME/WSEAS International Conference on Energy & Environmental Systems, Chalkia, Greece, May 8-10, 2006 (pp83-88) Moeling an analysis of parallel connecte permanent magnet synchronous

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

Both the ASME B and the draft VDI/VDE 2617 have strengths and

Both the ASME B and the draft VDI/VDE 2617 have strengths and Choosing Test Positions for Laser Tracker Evaluation an Future Stanars Development ala Muralikrishnan 1, Daniel Sawyer 1, Christopher lackburn 1, Steven Phillips 1, Craig Shakarji 1, E Morse 2, an Robert

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

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann Gravitation as the result of the reintegration of migrate electrons an positrons to their atomic nuclei. Osvalo Domann oomann@yahoo.com (This paper is an extract of [6] liste in section Bibliography.)

More information

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7. Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency

More information

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer Key Engineering Materials Online: 4-8-5 I: 66-9795, Vols. 7-73, pp 38-33 oi:.48/www.scientific.net/kem.7-73.38 4 rans ech ublications, witzerlan Citation & Copyright (to be inserte by the publisher imulation

More information

Chapter 9 Method of Weighted Residuals

Chapter 9 Method of Weighted Residuals Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe

More information

A Design Method for USBL Systems with Skew Three-element Arrays

A Design Method for USBL Systems with Skew Three-element Arrays A Design Metho for USBL Sstems with Skew Three-element Arras MIKHAIL ARKHIPOV Department of Postgrauate Stuies Technological Universit of the Mixteca Carretera a Acatlima Km. 2.5, Huajuapan e Leon, Oaxaca,

More information

IPMSM Inductances Calculation Using FEA

IPMSM Inductances Calculation Using FEA X International Symposium on Inustrial Electronics INDEL 24, Banja Luka, November 68, 24 IPMSM Inuctances Calculation Using FEA Dejan G. Jerkan, Marko A. Gecić an Darko P. Marčetić Department for Power,

More information

Single Arm, Centrifugal, Water Turbine for Low Head and Low Flow Application: Part 1- Theory and Design

Single Arm, Centrifugal, Water Turbine for Low Head and Low Flow Application: Part 1- Theory and Design Energy an Power 2018, 8(2): 51-55 DOI: 10.5923/j.ep.20180802.03 Single Arm, Centrifugal, Water Turbine for Low ea an Low Flow Application: Part 1- Theory an Design Kiplangat C. Kononen 1, Augustine B.

More information

DOUBLE PENDULUM VIBRATION MOTION IN FLUID FLOW

DOUBLE PENDULUM VIBRATION MOTION IN FLUID FLOW ENGINEERING FOR RURA DEEOPMENT Jelgava,.-.5.5. DOUBE PENDUUM IBRATION MOTION IN FUID FOW Janis iba, Maris Eiuks, Martins Irbe Riga Technical University, atvia janis.viba@rtu.lv, maris.eiuks@rtu.lv, martins.irbe@rtu.lv

More information

Simple Electromagnetic Motor Model for Torsional Analysis of Variable Speed Drives with an Induction Motor

Simple Electromagnetic Motor Model for Torsional Analysis of Variable Speed Drives with an Induction Motor DOI: 10.24352/UB.OVGU-2017-110 TECHNISCHE MECHANIK, 37, 2-5, (2017), 347-357 submitte: June 15, 2017 Simple Electromagnetic Motor Moel for Torsional Analysis of Variable Spee Drives with an Inuction Motor

More information

Optimization of Geometries by Energy Minimization

Optimization of Geometries by Energy Minimization Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.

More information

Mathematical model for determination of strand twist angle and diameter in stranded-wire helical springs

Mathematical model for determination of strand twist angle and diameter in stranded-wire helical springs Journal of Mechanical cience an Technology 4 (6) (1) 13~11 www.springerlink.com/content/1738-494x DOI 1.17/s16-1-333-4 Mathematical moel for etermination of stran twist angle an iameter in strane-wire

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

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables*

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables* 51st IEEE Conference on Decision an Control December 1-13 212. Maui Hawaii USA Total Energy Shaping of a Class of Uneractuate Port-Hamiltonian Systems using a New Set of Close-Loop Potential Shape Variables*

More information

Existence of equilibria in articulated bearings in presence of cavity

Existence of equilibria in articulated bearings in presence of cavity J. Math. Anal. Appl. 335 2007) 841 859 www.elsevier.com/locate/jmaa Existence of equilibria in articulate bearings in presence of cavity G. Buscaglia a, I. Ciuperca b, I. Hafii c,m.jai c, a Centro Atómico

More information

Polynomial Inclusion Functions

Polynomial Inclusion Functions Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl

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

Active control based on finite element analysis. Li Hongying

Active control based on finite element analysis. Li Hongying 5th International Conference on ucation, Management, Information an Meicine (MIM 5) Active control base on finite element analysis Li Hongying Automotive ngineering institute, Jiangxi University of Technology,

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

Multivariable Calculus: Chapter 13: Topic Guide and Formulas (pgs ) * line segment notation above a variable indicates vector

Multivariable Calculus: Chapter 13: Topic Guide and Formulas (pgs ) * line segment notation above a variable indicates vector Multivariable Calculus: Chapter 13: Topic Guie an Formulas (pgs 800 851) * line segment notation above a variable inicates vector The 3D Coorinate System: Distance Formula: (x 2 x ) 2 1 + ( y ) ) 2 y 2

More information

Lecture 6: Control of Three-Phase Inverters

Lecture 6: Control of Three-Phase Inverters Yoash Levron The Anrew an Erna Viterbi Faculty of Electrical Engineering, Technion Israel Institute of Technology, Haifa 323, Israel yoashl@ee.technion.ac.il Juri Belikov Department of Computer Systems,

More information

Situation awareness of power system based on static voltage security region

Situation awareness of power system based on static voltage security region The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran

More information

Attitude Control System Design of UAV Guo Li1, a, Xiaoliang Lv2, b, Yongqing Zeng3, c

Attitude Control System Design of UAV Guo Li1, a, Xiaoliang Lv2, b, Yongqing Zeng3, c 4th National Conference on Electrical, Electronics an Computer Engineering (NCEECE 205) Attitue Control ystem Design of UAV Guo Li, a, Xiaoliang Lv2, b, Yongqing eng3, c Automation chool, University of

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

Dynamics of the synchronous machine

Dynamics of the synchronous machine ELEC0047 - Power system ynamics, control an stability Dynamics of the synchronous machine Thierry Van Cutsem t.vancutsem@ulg.ac.be www.montefiore.ulg.ac.be/~vct These slies follow those presente in course

More information

Formulation of the Equilibrium Equations of Transversely Loaded Elements Taking Beam-Column Effect into Consideration

Formulation of the Equilibrium Equations of Transversely Loaded Elements Taking Beam-Column Effect into Consideration Journal of Emerging Trens in Engineering an Applie Sciences (JETEAS) 3 (): 0-07 Scholarlink Research Institute Journals, 0 (ISSN: 4-70) jeteas.scholarlinkresearch.org Journal of Emerging Trens in Engineering

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

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain Nonlinear Aaptive Ship Course Tracking Control Base on Backstepping an Nussbaum Gain Jialu Du, Chen Guo Abstract A nonlinear aaptive controller combining aaptive Backstepping algorithm with Nussbaum gain

More information

Design of an Industrial Distributed Controller Near Spatial Domain Boundaries

Design of an Industrial Distributed Controller Near Spatial Domain Boundaries Design of an Inustrial Distribute Controller Near Spatial Domain Bounaries Stevo Mijanovic, Greg E. Stewart, Guy A. Dumont, an Michael S. Davies ABSTRACT This paper proposes moifications to an inustrial

More information

Role of parameters in the stochastic dynamics of a stick-slip oscillator

Role of parameters in the stochastic dynamics of a stick-slip oscillator Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society

More information

CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS. Jože Kutin *, Ivan Bajsić

CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS. Jože Kutin *, Ivan Bajsić Sensors an Actuators A: Physical 168 (211) 149-154 oi: 1.116/j.sna.211..7 211 Elsevier B.V. CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS Jože Kutin *, Ivan Bajsić Laboratory of

More information

The Efficiency Optimization of Permanent Magnet Synchronous Machine DTC for Electric Vehicles Applications Based on Loss Model

The Efficiency Optimization of Permanent Magnet Synchronous Machine DTC for Electric Vehicles Applications Based on Loss Model International Power, Electronics an Materials Engineering Conference (IPEMEC 015) The Efficiency Optimization of Permanent Magnet Synchronous Machine DTC for Electric Vehicles Applications Base on Loss

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

Gyroscopic matrices of the right beams and the discs

Gyroscopic matrices of the right beams and the discs Titre : Matrice gyroscopique es poutres roites et es i[...] Date : 15/07/2014 Page : 1/16 Gyroscopic matrices of the right beams an the iscs Summary: This ocument presents the formulation of the matrices

More information

Electromagnet Gripping in Iron Foundry Automation Part II: Simulation

Electromagnet Gripping in Iron Foundry Automation Part II: Simulation www.ijcsi.org 238 Electromagnet Gripping in Iron Founry Automation Part II: Simulation Rhythm-Suren Wahwa Department of Prouction an Quality Engineering, NTNU Tronheim, 7051, Norway Abstract This paper

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

SYNTHESIS ON THE ASSESSMENT OF CHIPS CONTRACTION COEFFICIENT C d

SYNTHESIS ON THE ASSESSMENT OF CHIPS CONTRACTION COEFFICIENT C d SYNTHESIS ON THE ASSESSMENT OF HIPS ONTRATION OEFFIIENT Marius MILEA Technical University Gheorghe Asachi of Iași, Machine Manufacturing an Inustrial Management Faculty, B-ul D Mangeron, nr 59, A, 700050,

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

Damage identification based on incomplete modal data and constrained nonlinear multivariable function

Damage identification based on incomplete modal data and constrained nonlinear multivariable function Journal of Physics: Conference Series PAPER OPEN ACCESS Damage ientification base on incomplete moal ata an constraine nonlinear multivariable function To cite this article: S S Kourehli 215 J. Phys.:

More information

STABILITY CONTROL FOR SIX-WHEEL DRIVE ARTICULATED VEHICLE BASED ON DIRECT YAW MOMENT CONTROL METHOD

STABILITY CONTROL FOR SIX-WHEEL DRIVE ARTICULATED VEHICLE BASED ON DIRECT YAW MOMENT CONTROL METHOD STABILITY CONTROL FOR SIX-WHEEL DRIVE ARTICULATED VEHICLE BASED ON DIRECT YAW MOMENT CONTROL METHOD 1 YITING KANG, 2 WENMING ZHANG 1 Doctoral Caniate, School of Mechanical Engineering, University of Science

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

Lecture Notes: March C.D. Lin Attosecond X-ray pulses issues:

Lecture Notes: March C.D. Lin Attosecond X-ray pulses issues: Lecture Notes: March 2003-- C.D. Lin Attosecon X-ray pulses issues: 1. Generation: Nee short pulses (less than 7 fs) to generate HHG HHG in the frequency omain HHG in the time omain Issues of attosecon

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

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More 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

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

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

Journal of Engineering Science and Technology Review 7 (1) (2014) Research Article

Journal of Engineering Science and Technology Review 7 (1) (2014) Research Article Jestr Journal of Engineering Science an Technology Review 7 () (04) 8 Research Article JOURNA OF Engineering Science an Technology Review www.jestr.org Research on the Measurement Error of MWIR Average

More information

Math 1272 Solutions for Spring 2005 Final Exam. asked to find the limit of the sequence. This is equivalent to evaluating lim. lim.

Math 1272 Solutions for Spring 2005 Final Exam. asked to find the limit of the sequence. This is equivalent to evaluating lim. lim. Math 7 Solutions for Spring 5 Final Exam ) We are gien an infinite sequence for which the general term is a n 3 + 5n n + n an are 3 + 5n aske to fin the limit of the sequence. This is equialent to ealuating

More 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

Developing a Method for Increasing Accuracy and Precision in Measurement System Analysis: A Fuzzy Approach

Developing a Method for Increasing Accuracy and Precision in Measurement System Analysis: A Fuzzy Approach Journal of Inustrial Engineering 6(00 5-3 Developing a Metho for Increasing Accuracy an Precision in Measurement System Analysis: A Fuzzy Approach Abolfazl Kazemi a, Hassan Haleh a, Vahi Hajipour b, Seye

More information

Modelling of an electroactive polymer actuator

Modelling of an electroactive polymer actuator vailable online at www.scienceirect.com Proceia Engineering 48 (1 ) 1 9 MMaM 1 Moelling of an electroactive polymer actuator Ákos czél a * a zéchenyi István University, Egyetem tér 1,Gy r, H-96, Hungary

More information

Modelling dowel action of discrete reinforcing bars in cracked concrete structures

Modelling dowel action of discrete reinforcing bars in cracked concrete structures Title Moelling owel action of iscrete reinforcing bars in cracke concrete structures Author(s) Kwan, AKH; Ng, PL; Lam, JYK Citation The 2n International Symposium on Computational Mechanics an the 12th

More information

ector ition 1. Scalar: physical quantity having only magnitue but no irection is calle a scalar. eg: Time, mass, istance, spee, electric charge, etc.. ector: physical quantity having both magnitue an irection

More information

DESIGN, CONSTRUCTION AND TEST OF A TWO-STAGE PLANETARY TRACTION SPEED REDUCER

DESIGN, CONSTRUCTION AND TEST OF A TWO-STAGE PLANETARY TRACTION SPEED REDUCER DEIGN, CONTRUCTION AND TET OF A TWO-TAGE LANETARY TRACTION EED REDUCER Euaro Lobo Lustosa Cabral Escola olitécnica a Universiae e ão aulo Departamento e Engenharia Mecatrônica e istemas Mecânicos Av. Mello

More information

The new concepts of measurement error s regularities and effect characteristics

The new concepts of measurement error s regularities and effect characteristics The new concepts of measurement error s regularities an effect characteristics Ye Xiaoming[1,] Liu Haibo [3,,] Ling Mo[3] Xiao Xuebin [5] [1] School of Geoesy an Geomatics, Wuhan University, Wuhan, Hubei,

More information

State Space Analysis of Power System Stability Enhancement with Used the STATCOM

State Space Analysis of Power System Stability Enhancement with Used the STATCOM tate pace Analysis of Power ystem tability Enhancement with Use the ACOM M. Mahavian () - G. hahgholian () () Department of Electrical Engineering, Islamic Aza University, Naein Branch, Esfahan, Iran ()

More information

VIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK

VIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK AIAA Guiance, Navigation, an Control Conference an Exhibit 5-8 August, Monterey, California AIAA -9 VIRTUAL STRUCTURE BASED SPACECRAT ORMATION CONTROL WITH ORMATION EEDBACK Wei Ren Ranal W. Bear Department

More information

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann Gravitation as the result of the reintegration of migrate electrons an positrons to their atomic nuclei. Osvalo Domann oomann@yahoo.com (This paper is an extract of [6] liste in section Bibliography.)

More information

In the usual geometric derivation of Bragg s Law one assumes that crystalline

In the usual geometric derivation of Bragg s Law one assumes that crystalline Diffraction Principles In the usual geometric erivation of ragg s Law one assumes that crystalline arrays of atoms iffract X-rays just as the regularly etche lines of a grating iffract light. While this

More information

IN most inverter fed electrical drives, asynchronous

IN most inverter fed electrical drives, asynchronous Comparison of Different Designs of ynchronous Reluctance Machines with High-Anisotropy Rotors Erich chmit, Member, IEEE Abstract In terms of torue capability, power factor an efficiency, synchronous reluctance

More information

Open Access An Exponential Reaching Law Sliding Mode Observer for PMSM in Rotating Frame

Open Access An Exponential Reaching Law Sliding Mode Observer for PMSM in Rotating Frame Sen Orers for Reprints to reprints@benthamscience.ae The Open Automation an Control Systems Journal, 25, 7, 599-66 599 Open Access An Exponential Reaching Law Sliing Moe Observer for PMSM in Rotating Frame

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

DYNAMIC PERFORMANCE OF RELUCTANCE SYNCHRONOUS MACHINES

DYNAMIC PERFORMANCE OF RELUCTANCE SYNCHRONOUS MACHINES Annals of the University of Craiova, Electrical Engineering series, No 33, 9; ISSN 184-485 7 TH INTERNATIONAL CONFERENCE ON ELECTROMECHANICAL AN POWER SYSTEMS October 8-9, 9 - Iaşi, Romania YNAMIC PERFORMANCE

More information