Controller Parameters Dependence on Model Information Through Dimensional Analysis

Transcription

1 Joint 48th IEEE Conference on Decision and Control and 28th Chinese Control Conference Shanghai, PR China, December 16-18, 2009 WeBIn49 Controller Parameters Dependence on Model Information hrough Dimensional Analysis P Balaguer, A Ibeas, C Pedret and S Alcántara Abstract In this article we make use of dimensional analysis in order to investigate the controller parameters dependence on model information ie model parameters he objective is to relate the influence that each model parameter has on each controller parameter In order to accomplish this goal the model transfer function is first analyzed using dimensional analysis and characterized by means of of dimensionless numbers Secondly each controller parameter is related with the model parameters As a result it is derived the underlaying structure that any homogeneous tuning rule must follow he general results are particularized for PID control of first order and second order systems he results can be applied to PID tuning and PID tuning rules comparison I INRODUCION Dimensional analysis [1] is a well known theory applied to physical problems which permits a variety of important results such as i to understand the physical phenomenon, not only the depending quantities but also their relations, ii to represent results in a more compact way by means of dimensionless numbers and iii to generalize experimental results, thus reducing the experimental requirements Control theory can beneficiate from dimensional analysis in a wide variety of ways as can be seen from the bibliography In [2] dimensional analysis is explicitly presented as a theory to be used in control theory related problems In particular dimensional analysis is used to define systems equivalence and to relate system sensitivity to dimensional concepts From a practical point of view dimensionless models of vehicle dynamics are derived and used for controller synthesis In [6] the Buckingham Π theorem is used in order to represent the PID controller parameters of first order plus dead time FOPD model by means of dimensionless numbers, with the aim of performing a numerical optimization of the tuning parameters Practical advances in control design are proposed in [3] In particular it is shown that by means of dimensional analysis it is possible to reduce the scheduling parameters of gain scheduling controllers, thus reducing the analysis and design complexity Dimensional analysis is also used to dump uncertainty representation in dimensionless parameters, what may lead to less conservative results he financial support received from the Spanish CIC programme under grants DPI is greatly recognized P Balaguer is with the Department of Industrial Systems Engineering and Design, Universitat Jaume I de Castelló, Castelló, Spain pbalague@esidujies A Ibeas, C Pedret and S Alcántara are with the Department of elecommunication and Systems Engineering, Autonomous University of Barcelona, Barcelona, Spain Recently, in [5] it is studied the physical dimensions of matrices in state-space models by means of dimensional analysis Moreover, throughout the control bibliography, it is interesting to note that dimensional analysis results are used, although not stated explicitly, in order to tackle some control issues For example in [4] dimensionless parameters of FOPD model are used to compare distinct autotuning methods From the previous review we can see that dimensional analysis permits to tackle distinct control theory issues Although it has been mainly applied to identification and modeling of systems and systems uncertainty, the potential results of dimensional analysis related to controller design and controller comparison seems to have received less attention In this article we apply dimensional analysis on the transfer function framework in order to determine the structural relation among model parameters and controller parameters It is shown that an underlaying structure in fact exists he main result is obtained by firstly representing the transfer function by means of dimensionless numbers, in such a way that the transfer function behavior is characterized using a reduced set of dimensionless parameters he following questions are answered: i Which are the physical dimensions of the transfer function parameters?, ii Is it possible to reduce the number of transfer function parameters by means of dimensionless numbers? If the answer is affirmative, how many parameters can be reduced?, iii What is a proper system of units for the transfer function parameters? Next, on the basis of the preceding results, the dependence of controller parameters on the transfer function dimensionless numbers is established In particular the following questions are answered: i What is the relationship among the controller parameters and the plant model parameters?, ii Do all the model parameters affect in the same way all the controller parameters? and iii Is it possible to state these relationships in a more compact way by means of dimensionless numbers? Finally the general results obtained are particularized for PID control of FOPD and second order models in order to show the practical benefits of the results obtained In particular the results obtained are useful for comparing PID tuning methods as the characterization of controllers and models by means of dimensionless numbers allows the reduction of parameters in the comparison procedure ie the dimensionless number / he contributions of the article are organized as follows /09/$ IEEE 1914

2 WeBIn49 In Section II the basic concepts of dimensional analysis are reviewed Next it is shown in Section III that transfer functions can be represented by means of a reduced set of dimensionless numbers In Section IV it is analyzed the controller parameters dependence on transfer function dimensionless numbers Finally the theoretical results obtained in the preceding Sections are applied to PID control of first and second order systems II DIMENSIONAL ANALSIS Dimensional Analysis [1] is a technique used extensively in physic problems eg fluid dynamics he idea on which dimensional analysis is based is that physical laws are dimensionally homogeneous ie do not depend on the choice of any basic units of measurement, that is the law is valid in any system of dimensions his leads to the fact that functions expressing physical laws have a fundamental property called generalized homogeneity or symmetry his property allows the number of function arguments to be reduced, therefore making the dependence simpler Following we briefly review the fundamental concepts of dimensional analysis and state the fundamental Buckingham Π heorem II-A Physical quantities, Units and Dimensions When we refer to the fact that the length of a cable is 5 meters we are stating the magnitude of a physical quantity However, we could also say that the cable length is 5000 millimeters he only difference being in the units used to describe the physical quantity We then have that: Physical quantity = numerical value Unit 1 hus in both cases we are referring to the same physical quantity ie the cable length is constant, being the only difference the numerical value due to the difference in units eg meters Vs millimeters On the other hand, in both cases we are describing a length he dimension is defined as a qualitative description of a sensory perception of a physical entity It is clear that in both cases we are dealing with the same dimensions We write the dimensions of a quantity φ between square brackets, that is [φ] For example, the quantity lc that describes the length of a cable has length dimensions ie L, that is [lc] = L On the other hand, the velocity dimension is [v] = L 1 where is time dimension he systems of units is the set of fundamental units sufficient for measuring the properties of the class of phenomena under consideration For example in order to study the dynamics of mechanical systems, one system of units could be mass, length and time ie M,L and and a different system of units could be length, force and time ie L,F and However both systems of units are equally valid to study the mechanical system It is necessary to stress that the dimensions are defined by the systems of units employed Finally the concept of independent dimensions is introduced Given a set of quantities, these quantities are said to have independent dimensions if none of these quantities have dimensions which can be represented in terms of a product of powers of the dimensions of the remaining quantities For example, density [ρ]=ml 3, velocity [v]=l 1 and force [f]=ml 2 have independent dimensions In fact there are no x and y that accomplish, for example, the following equation: II-B Buckingham Π heorem [ρ] = [v] x [f] y 2 he Buckingham Π heorem [1] states that if some governed parameter a is a function of n governing parameters, that is: a = fa 1,, a k, a k+1,, a n 3 in such a way that only the set of k parameters given by a 1,, a k have independent dimensions, then there exists n k independent dimensionless parameters Π j j [1, n k] such that the following relation can be established: Π = ΦΠ 1,, Π n k 4 Each one of the dimensionless numbers related with the dependent governing parameters is given by: Π j = a k+j a pj 1 arj k where the exponents p j r j are chosen such that the parameter Π j is dimensionless he dimensionless number related with the governed parameter is: a Π = a p 1 ar k Note that there are n k dimensionless numbers related with each one of the n k dependent governing parameters ie Π j j [1, n k], plus one dimensionless number related with the governed parameter a ie Π As a result the function f equivalently the governed parameter as a = fa 1,,a n can be written in terms of a function Φ of smaller number of variables in the following form: III 5 6 fa 1,, a n = a p 1 ar k ΦΠ 1,, Π n k 7 RANSFER FUNCION DIMENSIONAL ANALSIS In this section we characterize the transfer function by means of dimensional analysis theory he objective is to express a given transfer function by means of dimensionless numbers hus we are able to represent the transfer function information in the most compact way In order to express a transfer function by means of dimensionless numbers, first we find the dimensions of the transfer function parameters Secondly a suitable system of units for transfer function parameters is defined 1915

3 WeBIn49 III-A Analysis of ransfer Function Parameters Dimension he transfer function of a system described by a linear time invariant differential equation is defined as the quotient between the Laplace transform of the output and the Laplace transform of the input under the assumption of zero initial conditions In order to apply the Buckingham Π heorem to a transfer function, it is necessary to know the dimensions of the transfer function parameters his is shown by the following lemma: Lemma 1: Consider a generic transfer function Gs in pole-zero form with c Z and Gs = Ks cg zs g p s e s 8 g z s = g p s = z zj s + 1 j=1 p pi s + 1 i=1 with the following transfer function parameters K, zj, pi and satisfying p z c 0 Note that multiple complex conjugate poles and zeros are also allowed by appropriate complex values of the terms pi and zj hen the sets { pi } and { zj } must be closed under complex conjugation he dimensions of the transfer function parameters are: [K] = U 1 c [ pi ] = [ zj ] = [] = where is the output dimension, U the input dimension and the time dimension Proof: he result of the lemma 1 follows from the following facts First note that [g z s] = 1, that is, dimensionless In fact the dimensions of zj are time dimensions ie they are time constants, thus [ zj ] = On the other hand, the dimensions of the s operator are [s] = 1 It can be seen either in the time domain where the s operator is transformed to operator p defined as p = d/dt or in the frequency domain where s is transformed to jω hus each product zl s yields a dimensionless number he same reasoning applies to g p s and e s Furthermore the dimensions of the transfer function are [Gs] = U 1 As a conclusion of the above analysis it results that [Gs] = [Ks c ], or equivalently [K] = [Gs][s c ] 1 which yields: [K] = U 1 c 9 As a result, the system of units required to study the inputoutput relation given by a transfer function is U ie output dimension, input dimension and time dimension In fact any variable or parameter dimensions can be obtained by a combination of the U dimensions III-B ransfer Function Dimensionless Numbers In this section we make use of the Buckingham Π theorem presented in Section II in order to characterize the transfer function by means of dimensionless numbers As a result we obtain the minimum number of dimensionless numbers required to characterize a transfer function he number of dimensionless numbers is equal to n k where n is the number of governing parameters and k is the number of parameters which have independent dimensions We can consider the following situations regarding the number of physical quantities n considered in the relationship and the number of independent dimensions k: n > k: In this case the number of quantities n appearing in the relating equation is greater than the independent dimensions k As a result we can find n k dimensionless numbers than equivalently relate the above relation n = k: In this case the number of physical quantities equals the number of independent dimensions then n k = 0, so there are no dimensionless numbers n < k: In this case the number of quantities n is less than the number of independent dimensions k his case is not physically possible as the resulting equation is no longer homogenous We now study the number of dimensionless numbers that can be derived from a general transfer function presented in 8 In order to accomplish this goal, it is necessary to calculate the number of governing variables n and the number of variables with independent dimensions k First of all we divide the parameters in three groups: 1 Signals: the output y and the input u 2 ransfer Function parameters: K, zj, pi, and 3 Frequency variable: s If we consider the output as the governed variable, we have: 1916 y = fk,, p1,, pp, z1,, zz, s, u 10 he number of parameters can be calculated as follows: y: he output of a transfer consists of only one parameter, that is #y=1 ie # refers to cardinality u: he input also consists of one parameter then #u=1; K: he system gain is again another single parameter #K=1 : he system delay quantity #=1 if there is time delay In case there is not time delay, #=0 p: he system number of poles with associated parameter is given by p In fact, if c < 0 there are c poles at the origin but with no associated parameter z: he system number of zeros with associated parameter is given by z If c > 0 there are c zeros at the origin but with no associated parameter s: he complex variable is also a parameter governing the output, then #s=1

4 WeBIn49 hen given the relation s = GsUs we find that the number of governing quantities we consider y the governed quantity is given by: n = #u + #K + # + p + z + #s 11 Considering a SISO system we have #u=#k=1 Moreover #s=1 hen it follows that n=3+# +p+z On the other hand as derived in the preceding section, we have that the fundamental system units U he number of governing parameters with independent dimensions is k=3 his is easily seen by considering the input [u] = U, the gain [K] = U 1 c and any other governing parameter, as the rest of them have time dimensions As a result, given any relation described by any transfer function Gs we can expect the following number of dimensionless numbers: n k = 3 + # + m + n 3 = # + p + z 12 hus, we can see that the number of dimensionless numbers is a function of the time parameters, that is the delay and the number of poles and zeros of the transfer function III-C Examples Example 1: Consider a first order system plus time delay: s = Ke s Us 13 s + 1 Firstly we state the governed parameter and the governing parameters hese are: = fu, K,,, s 14 he number of dimensionless numbers is n k = # + p+z In the FOPD case we have that #=1, #p=1 and #z=0, thus n-k=2 As a result we have two dimensionless number: Π 1 =, Π 2 = s 15 Note that the dimensionless numbers are not unique as / and s are also valid here is always an extra dimensionless number related with the governed variable, in this case: Π = K 1 y u 16 hen we have the following relationship Π = ΦΠ 1, Π 2, which can also be written as: = KUΦ, s 17 with the function Φ given by: KU = 1 Π e Π1Π2 18 what can be easily shown by direct substitution In the rest of the article we use the following more meaningful notation to refer to each dimensionless number: yielding finally = Π 1 =, s = Π 2 = s 19 KU = 1 e s 20 s + 1 Example 2: Consider now a second order system without time delay, that is: s = K Us 21 s s + 1 he relations among quantities considered is: = fu, K,, 2, s 22 he number of dimensionless numbers is n k = # +p+ z Now we have that #=0, #p=2 and #z=0, thus the number of dimensionless numbers is again n-k=2 he dimensionless numbers are: Π 1 = 2, Π 2 = s = s 23 he dimensionless number related with the governed variable is the same as the example before, so we have: = KUΦ, s 24 with function Φ given by KU = 1 s + 1Π 1 s Example 3: Consider now a second order system with time delay, that is: s = Ke s Us 26 s s + 1 he relations among quantities considered is: = fu, K,,, 2, s 27 he number of dimensionless numbers is n k = # + p + z In this case we have that #=1, #p=2 and #z=0, thus the number of dimensionless numbers is now n-k=3 he dimensionless number are: Π 1 =, Π 2 = 2, Π 3 = s 28 he dimensionless number related with the governed variable is the same as the example before, so we have: = KUΦ, 2, s 29 with function Φ given by KU = e Π1 s s + 1Π 2 s

5 WeBIn49 IV CONROLLER PARAMEERS DEPENDENCE ON DIMENSIONLESS NUMBERS In this section we apply the dimensional analysis theory in order to analyze the controller parameters dependence on the model transfer function parameters We then consider a generic controller Cs given by: with p Z and Cs = K c s q gc z s g c p s 31 zc gz c s = zj c s + 1 g c ps = j=1 pc pis c + 1 i=1 with the following transfer function parameters K c, zj c and pi c he relative order of the system is pc zc q 0 Note that complex poles and zeros are also allowed by appropriate complex values of the terms pi c and zj c In this case, the sets {pi c } and { zj c } must be closed under complex conjugation hen we are faced with the problem of finding the controller parameters dependence on the model parameters, that is: K c = f 1 K,, p1,, pp, z1,, zz 32 c p1 = f 2 K,, p1,, pp, z1,, zz 33 c ppc = f pc+1 K,, p1,, pp, z1,, zz 34 c z1 = f pc+2 K,, p1,, pp, z1,, zz 35 c zzc = f pc+zc+1 K,, p1,, pp, z1,, zz 36 Remark 1 Note that the governing parameters of each one of the controller parameters is a subset of the governing parameters of the system output as can be seen in equation 10 In fact, the controller parameters do not depend on the input u and on the complex variable s Moreover note that in equation 10 the parameters K and u, due to their dimensions, only affect the governed variable y Remark 2 he governing parameters on equations only have two distinct dimensions, the gain dimension K and the time dimension Moreover there is just one parameter with gain dimension, the transfer function gain K It then follows, due to dimensional homogeneity, that the gain can not appear in the rest of the governing variables all with time dimension in order to make them dimensionless In fact the only parameter that can be made dimensionless by means of the model gain K is the controller gain K c From the above remarks it can be seen that each one of the controller parameters ie governed variables are a function of n k 1 dimensionless numbers hese dimensionless numbers ie Π 1,, Π n k 1 are the ones generated from all transfer function time parameters ie, p1,, pp, z1,, zz As a result we can write KK c = Φ 1 Π 1,, Π n k 1 p1 c = Φ 2 Π 1,, Π n k 1 c ppc z1 c c zzc = Φ pc+1 Π 1,, Π n k 1 = Φ pc+2 Π 1,, Π n k 1 = Φ pc+zc+1 Π 1,, Π n k 1 Note that each one of the controller time constants are divided by his is so because in order to form the dimensionless numbers Π 1,, Π n k 1, the governing parameter was chosen to be the one forming the base of parameters with independent dimension However this selection is arbitrary In the following we particularize this general result for PID controllers of FOPD models and SOPD models V APPLICAION O PID CONROL In what follows we apply the general results stated in Section IV to two distinct model plants to be controlled with PID control V-A First order plus time delay systems In this section we tackle the problem of the PID controller parameter dependence on a FOPD model, with transfer function Gs given by: Gs = Ke s s Recall from example 1 that in this case there is just one dimensionless number related with the transfer function time parameters, which is As a result we have that KK c = Φ 1 i = Φ 2 d = Φ 3 hus, the controller gain depends in an inverse way on the process gain and on an arbitrary function Φ which defines the tuning rule of the dimensionless parameter / he integral time and the derivative time do not depend on the model gain K but are an arbitrary function of the dimensionless parameter / 1918

6 WeBIn49 V-B Second order plus time delay systems In this section we study the PID controller parameters dependence on the second order model plus delay given by: Gs = Ke s s s aking the first time constant as the independent dimension variable we have find the following dimensionless numbers of the transfer function time parameters and hen the controller parameters are obtained as: K c = K 1 Φ 1, 1 i = Φ 2, d = Φ 3, As a result, an equivalent relation as the one obtained in the preceding case is obtained However now the functions Φ i i [1, 2, 3] defining each one of the controller parameters is a function of two dimensionless variables V-C hird order systems In this section we study the PID controller parameters dependence on the third order system given by: Gs = K s s s aking the first time constant as the independent dimension variable we have find the following dimensionless numbers of the transfer function time parameters and 3 hen the controller parameters are obtained as: K c = K 1 Φ 1, 3 i = Φ 2, 3 d = Φ 3, 3 We can see that the controller parameters are again a function of two dimensionless parameters VI APPLICAION EXAMPLE In this section we compare Pade approximations of time delay in FOPD models as the one presented by equation 20 By straightforward calculation the transfer function of the dimensionless approximation error Ē s = G s Ĝ s for a first order Pade approximation is given by Ē s = 1 s + 1 e s 2 s 1 2 s where Ē s is the normalized error transfer function defined s Ŷ s as KU s, with s the FOPD output, Ŷ s the output of the approximate model, U s the model input and K the FOPD gain s is the normalized laplace variable defined as s = s, with the FOPD time constant is the dimensionless time delay defined as = / It can be seen that the normalized error transfer function is a function of just one parameter, the dimensionless time delay As a result the H norm of normalized error transfer function Ē s can be calculated for distinct Pade approximation orders as it is shown in Figure 1 Ē s st Order Pade Approximation 2nd Order Pade Approximation 3rd Order Pade Approximation 4th Order Pade Approximation Ē s / Fig 1 Comparison of Pade Approximations of several Orders VII CONCLUSIONS In this article we have used dimensional analysis to investigate two important aspects of control theory On the one hand it has been established that given a general transfer function, it is possible to characterize it by means of dimensionless numbers he number of dimensionless numbers obtained are equal to the transfer function time parameters On the other hand, dimensional analysis has been used to determine the controller parameters dependence on the model information ie model parameters hese results are general and of interest in order to characterize plant behaviour for analysis purposes as well as to derive PID tuning rules REFERENCES [1] G I Barenblatt Dimensional Analysis Gordon and Breach Science Publishers, 1987 [2] Sean Brennan On Size and Control: he Use of Dimensional Analysis in Controller Design PhD thesis, University of Illinois at Urbana- Champaign, 2002 [3] Haftay Hailu Dimensional ransformation: A Novel Method for Gain Scheduling and Robust Control PhD thesis, he Pennsylvania State University, 2006 [4] A Leva Comparative study of model-based PID autotuning methods In American Control Conference, 2007 [5] H J Palanthandalam-Madapusi, D S Berstein, and R Venugopal Dimensional analysis of matrices IEEE Control Systems Magazine, December: , 2007 [6] S avakoli and M avakoli Optimal tuning of PID controllers for first order plus time delay models using dimensional analysis In International Conference on Control and Automation,