Sliding Mode Flow Rate Observer Design

Similar documents
A Propagating Wave Packet Group Velocity Dispersion

Homotopy perturbation technique

Data Assimilation 1. Alan O Neill National Centre for Earth Observation UK

Full Waveform Inversion Using an Energy-Based Objective Function with Efficient Calculation of the Gradient

Einstein Equations for Tetrad Fields

Higher order derivatives

Cramér-Rao Inequality: Let f(x; θ) be a probability density function with continuous parameter

Studies of Turbulence and Transport in Alcator C-Mod Ohmic Plasmas with Phase Contrast Imaging and Comparisons with GYRO*

Elements of Statistical Thermodynamics

Torque Ripple Minimization in Switched Reluctance Motor Using Passivity-Based Robust Adaptive Control

NEW APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA

Rotor Stationary Control Analysis Based on Coupling KdV Equation Finite Steady Analysis Liu Dalong1,a, Xu Lijuan2,a

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.

Pipe flow friction, small vs. big pipes

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula

MEASURING HEAT FLUX FROM A COMPONENT ON A PCB

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002

Massachusetts Institute of Technology Department of Mechanical Engineering

Exercise 1. Sketch the graph of the following function. (x 2

The Matrix Exponential

Anti-Synchronization of Tigan and Li Systems with Unknown Parameters via Adaptive Control

COMPUTER GENERATED HOLOGRAMS Optical Sciences 627 W.J. Dallas (Monday, April 04, 2005, 8:35 AM) PART I: CHAPTER TWO COMB MATH.

Collisions between electrons and ions

(Upside-Down o Direct Rotation) β - Numbers


Full Order Observer Controller Design for Two Interacting Tank System Based on State Space Approach

INTEGRATION BY PARTS

1 Minimum Cut Problem

4.2 Design of Sections for Flexure

Recall that by Theorems 10.3 and 10.4 together provide us the estimate o(n2 ), S(q) q 9, q=1

The Matrix Exponential

Design and Simulation of PMSM Speed Regulating System Based on SMVSC

2008 AP Calculus BC Multiple Choice Exam

1.2 Faraday s law A changing magnetic field induces an electric field. Their relation is given by:

Two Products Manufacturer s Production Decisions with Carbon Constraint

A Prey-Predator Model with an Alternative Food for the Predator, Harvesting of Both the Species and with A Gestation Period for Interaction

The Relativistic Stern-Gerlach Force C. Tschalär 1. Introduction

Brief Introduction to Statistical Mechanics

Dynamic Modelling of Hoisting Steel Wire Rope. Da-zhi CAO, Wen-zheng DU, Bao-zhu MA *

Iterative learning control with initial rectifying action for nonlinear continuous systems X.-D. Li 1 T.W.S. Chow 2 J.K.L. Ho 3 J.

u 3 = u 3 (x 1, x 2, x 3 )

Where k is either given or determined from the data and c is an arbitrary constant.

Introduction to Condensed Matter Physics

EXST Regression Techniques Page 1

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory

Math 34A. Final Review

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012

Solution of Assignment #2

Exam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator.

ECE602 Exam 1 April 5, You must show ALL of your work for full credit.

Discrete Hilbert Transform. Numeric Algorithms

Construction of asymmetric orthogonal arrays of strength three via a replacement method

EEO 401 Digital Signal Processing Prof. Mark Fowler

Sundials and Linear Algebra

On the irreducibility of some polynomials in two variables

Procdings of IC-IDC0 ( and (, ( ( and (, and (f ( and (, rspctivly. If two input signals ar compltly qual, phas spctra of two signals ar qual. That is

ME469A Numerical Methods for Fluid Mechanics

In this lecture... Subsonic and supersonic nozzles Working of these nozzles Performance parameters for nozzles

SER/BER in a Fading Channel

Recursive Estimation of Dynamic Time-Varying Demand Models

State-space behaviours 2 using eigenvalues

Network Congestion Games

AS 5850 Finite Element Analysis

Contemporary, atomic, nuclear, and particle physics

5.80 Small-Molecule Spectroscopy and Dynamics

Dealing with quantitative data and problem solving life is a story problem! Attacking Quantitative Problems

Numerical considerations regarding the simulation of an aircraft in the approaching phase for landing

Title: Vibrational structure of electronic transition

Evaluating Reliability Systems by Using Weibull & New Weibull Extension Distributions Mushtak A.K. Shiker

Supplementary Materials

That is, we start with a general matrix: And end with a simpler matrix:

Mathematics. Complex Number rectangular form. Quadratic equation. Quadratic equation. Complex number Functions: sinusoids. Differentiation Integration

Neutrino Mass and Forbidden Beta Decays

COUNTING TAMELY RAMIFIED EXTENSIONS OF LOCAL FIELDS UP TO ISOMORPHISM

Estimation of apparent fraction defective: A mathematical approach

SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero.

Lecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields

APPROXIMATION THEORY, II ACADEMIC PRfSS. INC. New York San Francioco London. J. T. aden

The pn junction: 2 Current vs Voltage (IV) characteristics

Inference Methods for Stochastic Volatility Models

Analysis of Chattering in Sliding Mode Control Systems with Continuous Boundary Layer Approximation of Discontinuous Control

Gradient-based step response identification of low-order model for time delay systems Rui Yan, Fengwei Chen, Shijian Dong, Tao Liu*

Another view for a posteriori error estimates for variational inequalities of the second kind

Linear-Phase FIR Transfer Functions. Functions. Functions. Functions. Functions. Functions. Let

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

MCB137: Physical Biology of the Cell Spring 2017 Homework 6: Ligand binding and the MWC model of allostery (Due 3/23/17)

Differential Equations

Propositional Logic. Combinatorial Problem Solving (CPS) Albert Oliveras Enric Rodríguez-Carbonell. May 17, 2018

CO-ORDINATION OF FAST NUMERICAL RELAYS AND CURRENT TRANSFORMERS OVERDIMENSIONING FACTORS AND INFLUENCING PARAMETERS

A Sub-Optimal Log-Domain Decoding Algorithm for Non-Binary LDPC Codes

Computing and Communications -- Network Coding

[ ] 1+ lim G( s) 1+ s + s G s s G s Kacc SYSTEM PERFORMANCE. Since. Lecture 10: Steady-state Errors. Steady-state Errors. Then

Forces. Quantum ElectroDynamics. α = = We have now:

Probability and Stochastic Processes: A Friendly Introduction for Electrical and Computer Engineers Roy D. Yates and David J.

Deift/Zhou Steepest descent, Part I

A Recent Approach to Repetitive Control Strategy for Induced Draft Fan

Calculus II (MAC )

Determination of Vibrational and Electronic Parameters From an Electronic Spectrum of I 2 and a Birge-Sponer Plot

There is an arbitrary overall complex phase that could be added to A, but since this makes no difference we set it to zero and choose A real.

Indoor Autonomous Control of a Two-Wheeled Inverted Pendulum Vehicle Using Ultra Wide Band Technology

Transcription:

Sliding Mod Flow Rat Obsrvr Dsign Song Liu and Bin Yao School of Mchanical Enginring, Purdu Univrsity, Wst Lafaytt, IN797, USA liu(byao)@purdudu Abstract Dynamic flow rat information is ndd in a lot of applications; howvr it is normally not masurabl Calculation of th flow rat from fluid prssur dynamic quations usually rsults in poor stimats du to th vry noisy prssur masurmnts and unavoidabl modlling uncrtaintis This papr proposs a sliding mod dynamic flow rat obsrvr to stimat th mtr-in and mtr-out flow rats of a hydraulic cylindr Thortical convrgnc can b guarantd and simulation and xprimntal rsults show th ffctivnss of th proposd flow rat obsrvr Kywords Dynamic flow rat, sliding mod obsrvr Introduction Dynamic flow rat information is ndd in a lot of applications, such as automatd modlling of valv flow mapping, fault dtction, diagnostics and prognostics Howvr, with a convntional flow mtr which is only capabl of masuring static flow rat, th dynamic flow rat is not masurabl In ordr to obtain th flow rat information, on must sk othr approachs to avoid th dmand of dirct masurmnt In hydraulic systms activatd by hydraulic cylindrs, th flow rat is dirctly rlatd to th cylindr prssur dynamics and th motion of th cylindr rod Th cylindr prssur dynamics can b writtn as [] V P = Ax + Q + β () V P = A x Q + β whr x is th cylindr rod displacmnt, A and A ar th ram aras of th cylindr, P and P rprsnt prssurs at th two chambrs, Q and Q stand for th mtr-in and mtr-out flow rats rspctivly, β is th ffctiv bulk modulus, and V (x) and V (x) ar th total cylindr volums of th had and rod nds including conncting hos volums, and rprsnt th lumpd disturbanc flows that ar not associatd with th control valvs (i, th cylindr intrnal and xtrnal lakag flows, and so on) It is assumd that: a th cylindr rod position x and th vlocity x as wll as th prssurs P and P ar availabl, b th ffctiv bulk modulus β dos not chang or changs slowly, th actual valu of β may not b known but th bounds β min and β max ar known c th lumpd disturbanc flows and ar boundd by som constant known bounds δ and δ, rspctivly From () on may think to calculat flow rats from th prssur dynamic quations, g, V Q = Ax + P () β Though th prssur P is masurabl, it is vry noisy It is not practical to numrically diffrntiat P unlss a low pass filtr is usd, which would introduc phas lag Th accuracy of th calculation basd on () also dpnds on th accuracy of th ffctiv bulk modulus β, which is usually of larg varianc and vry hard to b dtrmind xactly Th apparanc of th disturbanc flow furthr complicats th problm A flow rat obsrvr may b a solution to this problm Howvr, bcaus th flow rat Q i, i= or, is not a stat of th systm dynamics, Lunbrgr obsrvr dos not work in this cas vn whn all paramtrs ar known and th disturbanc flow is ngligibl In ordr to dsign th flow rat obsrvr, on must ovrcom th following difficultis: a Th systm dynamics is not linar du to th chang in total cylindr volums V (x) and V (x) b Som systm paramtrs, such as β may not b xactly known c Th flow rat Q and Q ar inputs to th systm instad of bing systm stats This papr proposs a sliding mod flow rat obsrvr to ovrcom th nonlinarity and modl uncrtaintis

[,3,] Th proposd sliding mod flow rat obsrvr has strong thortical prformanc and robustnss It guarants convrgnc in finit tim whn thr is no modl uncrtaintis and boundd stimation rror in th prsnc of modlling rror Sliding Mod Flow Rat Obsrvr Sinc th obsrvrs for Q and Q ar basically sam, only th on for Q would b dsignd hr, th othr can b workd out in th sam way Sliding mod flow rat obsrvr dsign Lt θ =, th prssur dynamics () can b β rwrittn as V θ P = Ax + Q + (3) Th flow rat obsrvr, basd on th cylindr prssur dynamics (3), is givn by ˆ V θ Pˆ = Ax + Qˆ K( Pˆ P) () whr ˆP, ˆQ and ˆ θ rprsnt th stimats of P, Q and θ, rspctivly, K is a positiv obsrvr gain Subtract (3) from (), on can obtain: ( Pˆ P ) ( ˆ ) V ˆ θ + V θ θ P = Qˆ Q K( Pˆ P) (5) Dfin = Pˆ P as th prssur stimation rror ˆ as th paramtr stimation rror, and θ = θ θ rwrit (5) in trm of, V ˆ θ + K= Qˆ Q V θ P (6) Whn paramtrs lik th bulk modulus β is known and th disturbanc flow is ngligibl (i, assum = ), Eq (6) can b simplifid into: V ˆ θ + K= Qˆ Q (7) Sinc th mtrd flow rat Q must go through th valv orific, Q is always limitd by th valv or pump capability Mathmatically, Q is boundd by a known bound, i, Q Qmax In (7), ˆQ can b considrd as an input to th prssur stimation rror dynamics Thrfor, on may want to choos a propr ˆQ to mak th stimation rror (t) convrg to zro A discontinuous ˆQ is chosn as Qˆ = Q sign() (8) max Not that during th sliding mod whn =, mathmatically, (8) is dfind by th so-calld quivalnt control [,3] and can b approximatd obtaind by applying a low pass filtr to th discontinuous ˆQ : ˆ ˆ Q + ξω Q + ω Qˆ = Qˆ q n q n q whr ξ and ω n ar th damping ratio and natural frquncy of th filtr Thorm: Assuming paramtrs ar known and th disturbanc flow is ngligibl, th flow rat obsrvr () and (8) guarants that th prssur stimation rror (t) rachs zro in finit tim and, subsquntly, th quivalnt flow stimat Q ˆ q convrgs to th actual valv flow Q Proof: Dfin a positiv dfinit function V = V ˆ θ, whr V = V + Ax, V is th cylindr volum including hos and fitting volums whn th cylindr rod is fully rtractd, i, x = Diffrntiating V and noting (7), on can obtain: V = V ˆ θ + V ˆ θ ˆ = Ax θ K + ( Qˆ Q) ˆ = ( K Ax θ) ( Qmaxsign( ) + Q) (9) () Sinc th cylindr rod vlocity is always boundd du to finit flow capacity of any pump, it is asy to choos a K gratr than ˆ Axθ Thrfor th drivativ of V satisfis:

V < ( Q sign( ) + Q ) < η () max for som positiv η Thrfor, (t) rachs zro in finit tim and stays at zro Onc (t) is qual to zro, ˆP P as wll as ˆP P would b zro Hnc Eq (7) indicats th quivalnt flow stimat Q ˆ q would convrg to Q This nds th proof Robust prformanc analysis Whn thr is paramtric uncrtainty θ or th disturbanc flow in Eq (6), th drivativ of V would chang corrspondingly from () to th following: V = ( K Ax ˆ θ) ( Qmaxsign( ) + Q+ V θp + ) as () Sinc both V θ P and ar boundd, as long Q max is chosn larg nough to dominat Q as wll as V P θ and, th drivativ of V still satisfis: V < ( Qmax sign( ) + Q + V θp + ) < µ (3) for som positiv µ It is obvious that th prssur stimation rror still rachs zro in finit tim Aftr that, during th sliding mod whn = and =, from (6), th quivalnt flow stimat is Qˆ = Q + V θ P + () q Thus, thr will b a boundd valv flow stimation rror ˆ ( ) that dpnds on th lvl of of Q q Q = V x θp+ paramtric uncrtaintis and th disturbanc flow Howvr, th convrging rat of th valv flow stimation is similar to th prvious cas whn no modl uncrtaintis xist as sn from (3) whn compard to () 3 Rmarks Th thortically xcllnt robustnss of sliding mod control is prsrvd in th sliding mod obsrvr Th two main drawbacks of sliding mod control, i, larg control authority and control chattring, do not limit th sliding obsrvr s practical application This is du to th fact that th chattring and larg control authority issus in th sliding mod obsrvr dsign ar only linkd to numrical implmntation rathr than hard mchanical limitations [3] Th quivalnt flow stimat Q ˆ q is obtaind by snding ˆQ through a low pass filtr (LPF) to rmov th high frquncy chattring componnts and kp th ffctiv low frquncy componnts Th LPF usd in th sliding mod obsrvr has a diffrnt ffct from th on usd to filtr prssur masurmnts in () Th LPF usd in flow rat calculation is to smooth th masurmnt nois, which is almost whit and covrs th ntir frquncy rang from DC to Nyquist frquncy Th cut off frquncy of that LPF must b chosn vry low to mak th prssur diffrntiation possibl On th othr hand, th LPF usd to obtain quivalnt Q ˆ q is to rmov th chattring in ˆQ, whos frquncy mainly dpnds on th obsrvr gain K, th magnitud of th discontinuous trm Q max and th frquncy of digital implmntation Th largr K, Q max and digital implmntation frquncy ar, th highr th chattring frquncy is Practically, th obsrvr can run at a highr sampling frquncy than th control systm sampling frquncy, which would nabl th chatting at a frquncy vn highr than th control systm sampling frquncy Hnc LPF in th obsrvr dsign can hav a much highr cutoff frquncy and dos not introduc svr phas lag to th flow rat stimation 3 Simulations and Exprimnts Simulations ar don in Simulink A simpl PI controllr is usd to control a singl-rod doubl-acting hydraulic cylindr with a srvo valv, whos dynamics is nglctd Th cylindr rod is controlld to track a point-to-point rfrnc trajctory Whit nois is addd to prssur masurmnts

Th cylindr ram aras A and A ar m and m, rspctivly Th supply and tank prssurs ar st as 69KPa (PSI) and KPa, rspctivly Th ffctiv bulk modulus β is qual to 8 Th cylindr intrnal lakag is addd as th disturbanc flow: = K P P sign( P P ) (5) lak whr K lak is chosn to b -8 Th flow obsrvr paramtrs K and Q max ar st as K=- and Q max =-m 3 /s A scond ordr LPF with cut off frquncy qual to Hz is usd to rmov th chattring in ˆQ convrg to th actual valus vry wll without any phas lag, as prdictd by th thorm Th flow rat stimats also convrg vry wll xcpt for two spiks whr actual flow rats hav discontinuous changs Howvr, this may not b a problm bcaus flow rat would not chang discontinuously in practic du to th valv and flow dynamics P, Pa x 6 3 Simulation prssur with modlling rror Estimatd prssur Actual Prssur Estimation rror - 3 5 6 7 8 9 P, Pa 6 x 6 Simulation prssur w/o modlling rror Estimatd prssur Actual Prssur Estimation rror P, Pa 8 x 6 6 Estimatd prssur Actual Prssur Estimation rror P, Pa Mtr-in flow, m 3 /sc Mtr-out flow, m 3 /sc - 3 5 6 7 8 9 8 x 6 6 Estimatd prssur Actual Prssur Estimation rror - 3 5 6 7 8 9 x - - Tim, scond Fig Simulation prssurs w/o modlling rror Simulation flow rat, w/o modlling rror Estimatd Flow Rat Actual Flow Rat Estimation rror - 3 5 6 7 8 9 x - - - 3 5 6 7 8 9 Tim, scond Estimatd Flow Rat Actual Flow Rat Estimation rror Fig Simulation flow rats w/o modlling rror Figurs and show th prssur and flow rat stimats in th absnc of disturbanc flows It is obvious that th prssur stimats for both chambrs Mtr-in flow, m 3 /sc Mtr-out flow, m 3 /sc - 3 5 6 7 8 9 - Tim, scond Fig 3 Simulation prssurs with modlling rror x - Simulation flow rat, with modlling rror Estimatd Flow Rat Actual Flow Rat Estimation rror - 3 5 6 7 8 9 x - - - 3 5 6 7 8 9 Tim, scond Estimatd Flow Rat Actual Flow Rat Estimation rror Fig : Simulation flow rats with modlling rror Figurs 3 and illustrat th prssur and flow rat stimats in th prsnc of disturbanc flows Th prssur stimats, as prdictd by th thorm, ar not affctd by th disturbanc flows and convrg to th actual valus vry wll Th flow rat stimats, though affctd by th cylindr intrnal lakag, still show vry good convrgnc Exprimnts ar don with th cylindr usd to activat th swing motion of a thr dgr-of-frdom

lctro-hydraulic robot arm A nonlinar controllr is usd to control th swing motion to track an angular trajctory Sinc th cylindr rod motion in th xprimnt is diffrnt from th on in simulation, th xprimnt flow rats look diffrnt from th simulation Sinc thr is no way to know th actual flow rat, w can not compar th stimats with thir actual valu Th flow rat stimats ar shown in Fig 5 Mtr-in flow, m 3 /sc Mtr-out flow, m 3 /sc x - Exprimnt stimatd flow rat - - 3 5 6 7 8 9 x - - - 3 5 6 7 8 9 Tim, scond Fig 5: Exprimnt flow rat stimats Conclusions and Futur Works Th proposd sliding mod dynamic flow rat obsrvr taks advantag of th us of th quivalnt control for th stimation of flow rats with discontinuous trm It succssfully ovrcoms th nonlinar dynamics and uncrtain modlling rror, and thortically guarants that th valv flow rat stimats convrg to thir tru valus in th absnc of th disturbanc flow and th paramtr stimation rror, or a boundd stimation rror in th prsnc of th disturbanc flow and paramtr stimation rror Simulation and xprimntal rsults ar obtaind to illustrat th prformanc of th proposd flow rat obsrvr Th flow rat obsrvr is th first stp toward an automatd modlling mchanism of cartridg valv flow mapping Futur works includ adopting adaptiv robust obsrvr (ARO) tchniqu to dal with th changing paramtrs in th prssur dynamics, such as ffctiv bulk modulus β, and to rduc ffcts of paramtr variation and th disturbanc flow stimation rrors Acknowldgmnt Th work is supportd in part through th National Scinc Foundation undr th grant CMS-79 Rfrncs [] H E Mrritt, Hydraulic control systms, John Wily & Sons, 967 [] V I Utkin, Sliding Mods and Thir Application in Variabl Structur Systms Mir Publishrs, 978 [3] J J E Slotin, J K Hdrick and E A Misawa, On Sliding Obsrvrs for Nonlinar Systms ASME Journal of Dynamics Systms, Masurmnt, and Control, Vol 9, pp 5-5, 987 [] R A McCann, M S Islam, and I Husain, Application of a Sliding-Mod Obsrvr for Position and Spd Estimation in Switchd Rluctanc Motor Drivs IEEE Trans On Industry Applications, Vol 37, No, pp 5-58,