Position predictive measurement method for time grating CNC rotary table

Similar documents
FORECASTS GENERATING FOR ARCH-GARCH PROCESSES USING THE MATLAB PROCEDURES

Detection of Tire Lateral Force Based on a Resolver Mechanism

- The whole joint distribution is independent of the date at which it is measured and depends only on the lag.

Material Resistance and Friction in Cold Rolling

Interpretation of special relativity as applied to earth-centered locally inertial

Augmented Reality II - Kalman Filters - Gudrun Klinker May 25, 2004

Stopping Brownian Motion without Anticipation as Close as Possible to its Ultimate Maximum

Haar Wavelet Operational Matrix Method for Solving Fractional Partial Differential Equations

WATER LEVEL TRACKING WITH CONDENSATION ALGORITHM

Study and Evaluation of Performances of the Digital Multimeter

STATE-SPACE MODELLING. A mass balance across the tank gives:

Hall effect. Formulae :- 1) Hall coefficient RH = cm / Coulumb. 2) Magnetic induction BY 2

4. Electric field lines with respect to equipotential surfaces are

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

14 Autoregressive Moving Average Models

Description of the MS-Regress R package (Rmetrics)

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems

Chapter 5 Digital PID control algorithm. Hesheng Wang Department of Automation,SJTU 2016,03

Determination of the Sampling Period Required for a Fast Dynamic Response of DC-Motors

Institute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Notes on Kalman Filtering

The Paradox of Twins Described in a Three-dimensional Space-time Frame

Sub Module 2.6. Measurement of transient temperature

G. =, etc.

Non-uniform circular motion *

Air Traffic Forecast Empirical Research Based on the MCMC Method

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

Probabilistic Robotics SLAM

A Specification Test for Linear Dynamic Stochastic General Equilibrium Models

Fractional Method of Characteristics for Fractional Partial Differential Equations

Proposal of atomic clock in motion: Time in moving clock

1. Kinematics I: Position and Velocity

Chickens vs. Eggs: Replicating Thurman and Fisher (1988) by Arianto A. Patunru Department of Economics, University of Indonesia 2004

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

Lecture 1: Contents of the course. Advanced Digital Control. IT tools CCSDEMO

Direct teaching on industrial robot without force sensor Wei Gu1,a, Lei Shi1,b and Xincheng Tian1,c

An Analysis of the Impact of Valve Closure Time on the Course of Water Hammer

The equation to any straight line can be expressed in the form:

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

DERIVATION OF MASS, STIFFNESS AND DAMPING MATRICES OF A MODEL WING FROM EXPERIMENTAL MODAL DATA

Our main purpose in this section is to undertake an examination of the stock

Estimation of Kinetic Friction Coefficient for Sliding Rigid Block Nonstructural Components

Computer Vision. Motion Extraction

Tracking. Announcements

On asymptotic behavior of composite integers n = pq Yasufumi Hashimoto

MODIFIED COMPOUND SMOOTHER IN MEDIAN ALGORITHM OF SPAN SIZE 42

A New Perturbative Approach in Nonlinear Singularity Analysis

Chapter 7 Response of First-order RL and RC Circuits

A Video Vehicle Detection Algorithm Based on Improved Adaboost Algorithm Weiguang Liu and Qian Zhang*

( ) T = A New Approach to Identify Key Quality Characteristics in Multistage Manufacturing Process with Variation Transmission Model.

Doppler Method for Identification of Noise Sources on Underwater Moving Target

Time Domain Transfer Function of the Induction Motor

NEWTON S SECOND LAW OF MOTION

OBJECTIVES OF TIME SERIES ANALYSIS

Computing with diode model

Mean-square Stability Control for Networked Systems with Stochastic Time Delay

Embedded Systems and Software. A Simple Introduction to Embedded Control Systems (PID Control)

FAULT PROGNOSTICS AND RELIABILITY ESTIMATION OF DC MOTOR USING TIME SERIES ANALYSIS BASED ON DEGRADATION DATA

Sensors, Signals and Noise

Anti-Disturbance Control for Multiple Disturbances

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

Stationary Time Series

WHY YULE-WALKER SHOULD NOT BE USED FOR AUTOREGRESSIVE MODELLING

Lecture 2: Telegrapher Equations For Transmission Lines. Power Flow.

4.6 One Dimensional Kinematics and Integration

BOX-JENKINS MODEL NOTATION. The Box-Jenkins ARMA(p,q) model is denoted by the equation. pwhile the moving average (MA) part of the model is θ1at

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

RC, RL and RLC circuits

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

Position, Velocity, and Acceleration

Applications in Industry (Extended) Kalman Filter. Week Date Lecture Title

Licenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A

Tracking or Rejecting Rotational-Angle Dependent Signals Using Time Varying Repetitive Control

5. Stochastic processes (1)

JOY: The Journal of Yoga Summer 2008, Volume 7, Number 2

4.5 Constant Acceleration

Most Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.

An Invariance for (2+1)-Extension of Burgers Equation and Formulae to Obtain Solutions of KP Equation

Probabilistic Robotics SLAM

Dynamic Analysis of Damped Driven Pendulum using Laplace Transform Method

1. VELOCITY AND ACCELERATION

SELF-TUNING FUZZY-SLIDING CONTROLLER DESIGN OF VARIABLE DISPLACEMENT HYDRAULIC MOTOR SYSTEM WITHOUT POSITION FEEDBACK OF REGULATING CYLINDER

Block Diagram of a DCS in 411

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.

A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS

Class Meeting # 10: Introduction to the Wave Equation

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

PARAMETER ESTIMATION OF NEARLY NON- STATIONARY AUTOREGRESSIVE PROCESSES

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN

NUMERICAL CALIBRATION FOR 3-AXIS ACCELEROMETERS AND MAGNETOMETERS

2017 3rd International Conference on E-commerce and Contemporary Economic Development (ECED 2017) ISBN:

Chapter 3 Common Families of Distributions

International Industrial Informatics and Computer Engineering Conference (IIICEC 2015)

Cumulative Damage Evaluation based on Energy Balance Equation

ON DETERMINATION OF SOME CHARACTERISTICS OF SEMI-MARKOV PROCESS FOR DIFFERENT DISTRIBUTIONS OF TRANSIENT PROBABILITIES ABSTRACT

IE1206 Embedded Electronics

Silicon Controlled Rectifiers UNIT-1

Transcription:

Posiion redicive measuremen mehod for ime graing CC roary able Liu Xiaokang a, Peng Donglin a, Yang Wei a and Fei Yeai b a Engineering Research Cener of Mechanical Tesing Technology and Equimen, Minisry of Educaion, Chongqing Insiue of Technology, Chongqing 45, China b School of Insrumen Science and Oo-Elecronic Engineering, Hefei Universiy of Technology, Hefei, 239 China Tel.: +86-23-6256355 Fax: +86-23-6256353 E-mail: lxk@cqi.edu.cu Absrac Time graing sensor ransforms sace domain informaion o ime domain informaion and measures saial dislacemen wih ime. To develo high recision ime graing CC roary able and reduce he dynamic osiion feedback error of he able, circular osiion redicive measuremen mehod is roosed for ransforming ime domain informaion back o sace domain informaion based on ime-sace ransformaion echnology. Prediced values are obained by modeling he measured values wih ime series heory, and he las redicion error is correced in real ime using he nex measured values. Modeling mehod and arameer esimaion algorihm are resened. To conform he validiy of he osiion redicion mehod, an exerimenal sysem is designed. The angle dislacemen redicion error of he roary able is 2, and recise redicing is achieved. Keywords: Predicive measuremen, ime-sace ransformaion, ime graing, CC roary able. Inroducion Time graing is a novel dislacemen sensor develoed by he conce of measuring sace wih ime, and angular dislacemen is measured very accuraely wihou he use of a convenional graduaed mechanical scale, which ermis he grealy lowered cos of manufacure []. Time graing ransforms sace domain informaion o ime domain informaion based on ime-sace ransformaion echnology and measures saial osiion every uniform ime. In full closed-loo conrol mode CC sysem requires osiion feeding every uniform sace. To use ime graing as feedback comonen for CC osiioning servo conrol and develo high recision ime graing CC roary able, he new heories and mehods should o be develoed for ransforming ime domain informaion back o sace domain informaion based on ime-sace ransformaion echnology. Predicion is currenly used for indusrial conrol successfully, redicive conrol has become a yical rocess conrol mehod, and redicive conrol heory and algorihm are widely sudied [2]. Predicion conceion was firsly inroduced o measuremen field [3], and his aer will develo circular osiion redicive measuremen mehod o reduce he dynamic osiion feedback error of ime graing CC roary able. 2. Princile of redicive measuremen According o he measuremen rincile of ime graing [, 4], roaing magneic field is used o consruc a moving coordinae sysem S wih consan velociy V, he saial angular dislacemen beween a moving robe P a a random seed v and a fixed robe P b can be measured by deecing he ime difference T ha he roaing magneic field scans he wo robes (Fig. ). I.e. =VT=(36º/scanning eriod) ime difference. In he course of dynamic measuring, he ime inerval of S scans P b is fixed, and he signal of P b riggers daa samling, 2-278

hus ime graing ges a measured angle value every scanning eriod T i.e. measuring eriod. v V P a Pb P a P b T ime difference An absolue angle value of CC roary able is measured by ime graing every ime inerval T. The measured angle values θ i-+, θ i-+2,, θ i from ime oin T i-+ o T i can be reaed as a ime series. Time series model is based on he correlaion of he daa, which describes he dynamic characerisic and variaion rule of he series. Thus he fuure values of he ime series can be rediced by using ime series model [5]. The rincile of redicive measuremen mehod is shown in Fig. 2. A curren ime oin T i, he angle dislacemen of O θ i+ θ i θ i- v θ i-+2 θ i-+ T scanning eriod Fig.. Time graing measuremen model. able osiion redicing absolue angle T i measured value T i-+ T i-+2 T i- T i T i+ i rediced value incremenal ulses Fig. 2. Princile of redicive measuremen mehod. he roary able θ i during he nex measuring eriod T (from ime oin T i o T i+ ) can be rediced by modeling he absolue measured angle values θ i-+, θ i-+2,, θ i, hen incremenal ulses are ouued reresening he angle dislacemen θ i during he nex measuring eriod. In his way, discree absolue measured angle values can be ransformed o coninuous incremenal ulses ha are fed o CC sysem. The angle dislacemen of he roary able θ k during he kh measuring eriod (from ime oin T k- o T k ) is:. k k k () Using ime series heory he rediced angle dislacemen value of he roary able ˆi during he nex measuring eriod (from ime oin T i o T i+ ) can be obained as: ˆ L (,,, ). (2) i i i i i The number of ulse ouued in PWM (Pulse-Widh Modulaion) mode during he nex measuring eriod (from ime oin T i o T i+ ) is: 2-279

P ( ˆ e) Q, (3) i i i where Q is he ulse equivalen and e i is he redicion error of he las measuring eriod (from ime oin T i- o T i ). Time graing will ge a new measured angle value θ i+ when ime oin T i+ come. θ i+ can be used o calculae he redicion error beween he rediced value and he measured value. Eq. (3) shows ha he redicion error of he las measuring eriod is correced when curren redicing is conduced. And he redicion error will no be cumulaed, which assures he aforemenioned mehod of high recision. 3. Algorihm for osiion redicion Predicion model is based on saionary daa series. Firsly i is necessary o rerocess he dynamic daa measured by ime graing, remove he deerminisic rends and ransform non-saionary daa series o saionariy [5, 6]. And we aly saionary daa series o modeling, analysis and redicion for he residual daa hen finish modeling and analysis for curren and as angle measured values of ime graing ogeher wih he deerminisic rends. Finally redicing angle dislacemen of he CC roary able is achieved. The exression of h order auoregressive model AR() for ime series {X } can be described as: X a X,, (4) j j j where { } is whie noise W(, 2 ), real number a=(a, a 2,, a ) T is he auoregressive coefficien and a. X n+ can be successively rediced wih daa X n, X n-,, X n-+. The exression for oimal linear redicion is: n n n n n j n j j Xˆ L( X X, X,, X ) a X. (5) Parameer esimaion algorihm is as follows for auoregressive coefficien a. Firsly we need o rerocess he observed daa X, X 2,, X o obain a new ime series {Y } wih zero mean value: Y X X,,2,, X j X j. (6) Then we model he {Y } wih AR(). The esimaion for auocovariance funcion of he rerocessed daa is: k y y, k,,,. (7) k j jk j The square esimaion for auoregressive coefficien and whie noise T ( a 2, a2, a ), are deermined by Yule-Walker equaion a 2 2 a2 2 a (8) 2-28

angular acceleraion(/ms 2 ) angular acceleraion(/ms 2 ) angle() angular velociy(/ms) and 2 ( a a2 2 a ). (9) 4. Exerimen seu and resuls Exerimen seu is shown in Fig. 3. An AC servomoor drives a roary able insalled wih a.8 ime graing sensor acing as osiion feedback comonen. Siemens 82D digial CC servo sysem works wih SMC3 encoder inerface module for full closed-loo osiion conrol. The CC sysem receives he feedback signals from ime graing and hen Fig. 3. Exerimen seu. recisely conrols he osiion of he roary able. HEIDEHAI angular encoder ROD88 wih 36 line coun and accuracy of is mouned for esing he dynamic redicion error and he osiioning accuracy of he roary able. Inerolaion and digiizing elecronics IBV66B is alied for convering he ROD88 sinusoidal signals o TTL square-wave ulse rains as well as 4-fold subdivision. Consequenly he measuring se of he esing sysem 角度 7 6 5 4 3 2 2 3 4 5 6 7 8 2 3 4 5 6 7 8 9 ime(s) Fig. 4. Tesed osiion curve of roary able. 8 6 4 2 8 8 6 6 2 3 4 4 4 2 2 2 3 4 5 6 7 8 2 3 4 5 6 7 8 9 ime(s) Fig. 5. Velociy curve of roary able..35.3.25.2.5..5 -.5 -. 2 3 4 5 6 7 8 9 ime(s) Fig. 6. Acceleraion curve of roary able.. -. -.2 -.3 -.4 -.5 measured value rediced value redicion error -.6 -.7 5. 5.5 6. 6.5 7. 7.5 ime(s) Fig. 7. Predicion resuls of negaive acceleraion moion. 2-28

is.9. A couner card is develoed o rocess he ouued signals from IBV66B and he incremenal ulses from ime graing and can simulaneously lock he wo couning resuls for synchronous dislacemen comarison. The card is conneced o a comuer via a RS-232C or, and he comuer can auomaically acquire he wo locked daa. The comuer can also ge he measured angle value of ime graing direcly wih anoher RS-232C or. Figure 4 shows a esed osiion curve of roary able. The CC sysem drove he able roaing from 4.985 o 6.9872 wihin 7 seconds. The comuer samled he osiion of he roary able en imes er second. There were four differen moion saus including res, osiive acceleraion moion, aroximae uniform moion and negaive acceleraion moion (Fig. 5). The angular acceleraion curve of roary able is shown in Fig. 6. Figure 7 illusraes he redicion resuls of negaive acceleraion moion (he fourh moion course denoed in Fig. 5) using AR(3) model. The acceleraion redicion error varies form -.23/ms 2 o.9/ms 2. And he redicion error goes bigger along wih he larger acceleraion variaion. According o he relaionshi beween dislacemen and acceleraion, he angle dislacemen redicion error of he roary able is 2. 5. Conclusion This aer briefly inroduced he rincile of ime graing. Circular osiion redicive measuremen mehod was roosed o reduce he dynamic osiion feedback error of ime graing CC roary able. Prediced values were calculaed by modeling he measured values wih auoregressive model, and he redicion coefficien comuing algorihm was resened. The las redicion error was correced in real ime using he nex measured values. Discree absolue measured angle values were ransformed o coninuous incremenal ulses ha were fed o CC sysem. Exerimenal resuls have shown ha recise osiion redicing for ime graing roary able was achieved. Furher research will focus on coninuous redicing for comosie moion. 6. Acknowledgemen This rojec is funded by aional aural Science Foundaion of China (o. 5855). References. X.K. Liu, D.L Peng, X.H. Zhang and X.H. Chen. Solid Sae Phenomena: Mecharonic Sysems and Maerials. 26, vol. 3,. 435-44. 2. S. Joe Qin and Thomas A. Badgwell. Conrol Engineering Pracice. 23, vol.,. 733-764. 3. X.K. Liu, Y. Fei, D. Peng and X. Wang. Measuremen Technology and Inelligen Insrumens VIII: Key Engineering Maerials. 28, vols. 38-382,.43-46. 4. Liu Xiaokang, Peng Donglin, Zhu Ge, e al. Chinese Journal of Mechanical Engineering. 28, vol. 2,. 2-5. 5. Shuyuan He. Alicaions for Time Series Analysis. Beijing Universiy Press. China 27. 6. S. Franisek. Predicions in Time Series Using Regression Models. Sringer-Verlag, ew York. 22. 2-282