Adaptive Control - Autotuning Structure of presentation: Relay feedback autotuning outline Relay feedback autotuning details How close is the estimate of the ultimate gain and period to the actual ultimate gain and period? The effect of noise Determination of the gain and phase margins using relay autotuning Static load disturbance Commercial autotuning controller Question and Answer. Relay feedback autotuning - outline. Place the controller in proportional mode only (i.e. set T i to a maximum and to a minimum). T d. Increase K c until the closed loop system output goes marginally stable ; record K (calling it K, the ultimate gain), and the ultimate period, T. PI controller settings: K c 0.45K u c T 0.83T Ideal PID controller settings: K c 0.6K u T d 0.5T u Reminder of ultimate cycle tuning i T 0.5T i u u u Reference: Ziegler, J. and Nichols, N. (94). Optimum settings for automatic controllers, Transactions of the ASME, 64, 759-768. u Modified ultimate cycle method relay auto-tuning Åström and Hägglund (984) have developed an attractive alternative to the ultimate cycle method. In the relay auto-tuning method, a simple experimental test is used to determine K u and T u. For this test, the feedback controller is temporarily replaced by an on-off controller (or relay). After the control loop is closed, the controlled variable exhibits a sustained oscillation that is characteristic of on-off control. The operation of the relay auto-tuner includes a dead band; the dead band is used to avoid frequent switching caused by measurement noise. Reference: Åström, K. and Hägglund, T. (984). Automatic tuning of simple regulators with specification on the gain and phase margins, Automatica, 0, 645. 3 The relay auto-tuning method has several important advantages over the ultimate cycle method:. Only a single experiment test is required instead of a trial-and-error procedure.. The amplitude of the process output a can be restricted by adjusting relay amplitude d. 3. The process is not forced to a stability limit. 4. The experimental test is easily automated using commercial products. 4
. Relay feedback autotuning - details Relay feedback autotuning - details This oscillation is almost sinusoidal, depending on the filtering properties of G p (s). A comparison can be drawn to the oscillation obtained from the ultimate cycle experiment: When a relay is switched in, a sustained oscillation at c is observed: 5 6 Relay feedback autotuning - details Knowledge of K u and T u allows the controller to be retuned using simple tuning methods (such as those of Ziegler and Nichols). The relay autotuner will supply an approximate estimate of T u (labelled ). An equivalent approximate K u (labelled ) may be deduced from the relay autotuner using describing function analysis: Relay feedback autotuning - details The limit cycle is defined when G p (jω) intersects. N(A) 4d From the tables of describing functions, N(A), A πa half-peak amplitude of limit cycle output. 7 The relay autotuner can be considered equivalent to: 8
Relay feedback autotuning - details 3. How close is the estimate of the ultimate gain and period to the actual ultimate gain and period? An oscillatory output at c exists when 9 0 Simulation of relay autotuner
Estimate of ultimate gain and ultimate period from simulation Example Peak to peak amplitude ( A).34 i.e. A 0.67; d (in the simulation); Therefore,.90. From the simulated output, the estimate of the ultimate period 4 seconds. Overall: Such large % errors are related to the sinusoidal nature of the limit cycle output. 3 4 Simulation of relay autotuner 5 6
Estimate of ultimate gain and ultimate period from simulation Peak to peak amplitude ( A).3 i.e. A 0.57; d (in the simulation); 4. The effect of noise The presence of noise can cause difficulties in measuring the amplitude and period of the limit cycle output. In a simulation (noise amplitude: max. 0.; min. -0.): Therefore,.5. From the simulated output, the estimate of the ultimate period 4.3 seconds. Overall: More sinusoidal limit cycle output -> less percentage error. 7 8 One possibility: Use hysteresis on the relay Generally, the hysteresis width is made larger than the maximum noise level; if the maximum noise amplitude is ± 0., set up hysteresis on the relay to be ± 0.5, say. However, introducing hysteresis changes the amplitude and frequency of the controlled variable 9 0
Measuring amplitude and period With or without hysteresis, the challenge is to measure accurately the amplitude and period of the controlled variable. Without hysteresis: Estimated peak-peak amplitude 0.96, estimated period 4.6 s. With hysteresis: Estimated peak-peak amplitude.53, estimated period 6.5 s. Estimating ultimate gain If hysteresis is absent, the ultimate gain may be estimated as.65. The following table summarises the results: Ultimate gain Ultimate period s e G p(s) ( +.5s)( + 3s)( + 4s) Actual.3 3. Relay, no noise.5 (-3%) 4.3 (+9%) Relay, noise present.65 (+4%) 4.6 (+%) If hysteresis is present, N(A) [from describing function analysis] changes see table of describing functions previously. Thus, the ultimate gain is estimated by a different formula. This is left as an exercise. 5. Determination of the gain and phase margins using relay autotuning Reminder: Gain and phase margin Gain margin relay autotuning A relay autotuner may be used to estimate the gain margin of compensated systems. The method is set up as follows: 3 4
Gain margin relay autotuning Gain margin relay autotuning is approximated by: 5 6 Gain margin relay autotuning Example: A process, given by s ( +.5s)e Gp(s) + 8.5s +.5s + 8s.87s.95e is modelled as: Gm(s) + 6.69s 3 PI controller: G c (s).36 + 5.9s The PI controller was designed, using a standard method, to achieve the specifications of ± % settling time of 5 seconds 0 and a phase margin of 45. Gain margin relay autotuning We will use MATLAB to determine the gain margin of this compensated system, and compare the result to the estimate of the gain margin obtained from the relay test. 7 8
Gain margin relay autotuning Gain margin relay autotuning d (in the simulation); A.34 i.e. A.7; 4() + π(.7) Therefore, gain margin.09 (estimate) π(.7) Gain margin (from MATLAB simulation).035 Error in gain margin estimate +3%. 9 30 Phase margin relay autotuning A relay autotuner may be used to estimate the gain margin of compensated systems. The method is set up as follows: Phase margin relay autotuning The transfer function of the system in series with the relay is: This is equivalent to the following: 3 3
Phase margin relay autotuning Phase margin relay autotuning 33 34 Phase margin relay autotuning Phase margin relay autotuning Example: A process, given by s ( +.5s)e Gp(s) + 8.5s +.5s + 8s.87s.95e is modelled as: Gm(s) + 6.69s 3 PI controller: G + c (s).36 5.9s The PI controller was designed, using a standard method, to achieve the specifications of ± % settling time of 5 seconds 0 and a phase margin of 45. 35 We will use MATLAB to determine the phase margin of this compensated system, and compare the result to the estimate of the phase margin obtained from the relay test. Phase margin (see slide 8). 0 3 36
Phase margin relay autotuning Phase margin relay autotuning From plot, phase lag 6. seconds Period of waveform 4 seconds Phase lag (in degrees) (6./4.5)360 54 degrees Estimated phase margin 80 54 6 degrees. Phase margin (from MATLAB simulation) 30.9 degrees Error in phase margin estimate -6%. 37 38 6. Static load disturbance Static load disturbances during the relay tuning experiment introduce errors in the estimates of the ultimate gain and ultimate period. This section shows how an automatic bias can be introduced to overcome the problem. Static load disturbance A typical static load disturbance occurs in a heating and ventilation system if environmental conditions change: We consider the effects of a static load disturbance when the autotuner is a relay without hysteresis: Reference: Hang, C.C., Astrom, K.J. and Ho, W.K. (993). Relay auto-tuning in the presence of static load disturbance, Automatica, Vol. 9, No., pp. 563-564; also Hang, C.C., Lee, T.H. and Ho, W.K. (993), 39 Adaptive Control, Chapter 4. Thus, if outside temperature increases (for example), the heater does not need to be on as long to maintain the desired temperature. The ultimate period changes (as does the ultimate gain), which has a knock-on effect on the subsequent controller tuning. 40
Static load disturbance Static load disturbance The load disturbance may be determined, as follows: The DC component of the manipulated variable (process input), m dc, is The DC component of the controlled variable (process output), y dc is t + t ydc r + (y r) dt (t + t ) 0 DC component with no load disturbance DC component of the relay waveform 4 DC component without static load Average value of the oscillation 4 Static load disturbance Now, i.e. i.e. K y K p dc r K p p m dc t t + l + d r + t + t (t + t ) t t Kpl + Kpd t + t ( t + t ) t + t t + t 0 ( y r ) 0 (y r)dt dt 43 i.e. Static load disturbance t t t t l d + (y r) t + t K dt p + ( t + t ) To cancel the effect of the static load during the autotuning test, a bias, u b, equal to the negative of the estimated load should be added to the relay output: t t ub t + t d ^ K p 0 t t + ( ) t + t 0 (y r) dt Note the estimated value of process gain used. This bias term may be automatically incorporated into an existing relay autotuner. 44
s Static load e A process is given by: Gp(s) disturbance ( + s) When an autotuning relay is incorporated into the loop, the following response is determined (static load introduced after approximately 0 seconds, corrective bias term introduced after seconds): s Static load e A process is given by: Gp(s) disturbance ( + s) Autotuning relay is incorporated into the loop; oscillations with symmetrical positive and negative half-cycles occur as shown in the first 0 seconds of the data: 45 46 s Static load e A process is given by: Gp(s) disturbance ( + s) A small static load disturbance (of 0.08) is introduced at approximately t 0; the oscillations become asymmetrical for the next seconds. The error that would result in the estimated ultimate gain and ultimate period is +4% and +%, respectively. s Static load e A process is given by: Gp(s) disturbance ( + s) At t, a bias is applied based on an initial estimate of the process gain of 0.5 (noting that the actual process gain ). Exact symmetrical oscillations are not achieved; however, the asymmetry is slight, and the error that would result in the estimated ultimate gain and ultimate period is +% and -%, respectively. 47 48
7. Commercial autotuning controller Example: Fisher DPR 900 Controller Commercial autotuning controller Some features:. If no process knowledge exists, autotuning is preformed as follows: The process is brought to a desired operating point, either by the operator in manual mode or by the controller in automatic mode When the loop is stationary, the operator presses a tuning button. The PID controller is temporarily disconnected and the noise level is measured. Reference: Hang, C.C., Lee, T.H. and Ho, W.K. (993), Adaptive Control, Chapter 4. 49 50 Commercial autotuning controller After a short period, a relay with hysteresis is introduced into the loop. The hysteresis width of the relay is determined automatically from the noise level. During the oscillation, the relay amplitude is adjusted, so that a desired level of the oscillation amplitude is obtained. When an oscillation with constant amplitude and period is obtained, the relay experiment is interrupted and estimates of the ultimate gain and ultimate period will be made. Commercial autotuning controller. The PID controller parameters are then calculated from the ultimate gain and the ultimate period. Fast, medium or slow responses are available; for example, the medium tuning formula is K 0.35, T.3T and T 0.0T. Typical performance; G p ( + 5s)( + s) c K u i u d u 5 5
8. Question and Answer Answer 53 54 Answer Answer 55 56
Answer Answer 57 58 Answer 59