# Patent

## Title
Model-free adaptive control of quality variables

## Abstract

A Model-Free Adaptive Quality Variable control system is disclosed for effectively controlling quality variables on-line in closed-loop fashion. It is able to automatically control quality variables under the conditions where there are significant varying time delays and disturbances in the process. Because of its unique capability, the control system is useful for building flexible and adaptive production systems, achieving Six Sigma quality control goals, and fulfilling the on demand manufacturing needs in the new e-commerce environment.

## Background
This is a divisional of Application No. 09/832,688, filed Apr. 11, 2001, now U.S. Pat. No. 6,684,115, which claims the benefit of Provisional Application No. 60/196,206, filed Apr. 11, 2000.


FIELD OF THE INVENTION
The invention relates to automatic control of industrial processes, equipment, and facilities, and more particularly to a method and apparatus for adaptively controlling product quality related process variables so that “Six Sigma” or “zero-defects” quality goals can be reached.
BACKGROUND OF THE INVENTION
“Six Sigma” has become a new standard for quality control in the industrial world. General Electric alone has saved 5 billion dollars within 3 years since they introduced the program. Sigma (σ) is a Greek alphabet symbol used to represent the standard deviation in statistics. It is a good measurement for the process variability in statistical quality control.
Compared to the previous “Three Sigma” industrial standard that requires a process to produce 99.933% good parts, the “Six Sigma” standard requires the process to produce 99.9999966% good parts. This means that for every one million products, a process with Six Sigma capability will produce almost zero defective parts.
Industry surveys show that a process will typically shift 1.5 Sigma from its center. When this happens, a process with Three Sigma capability will start to produce thousands of bad parts per million while a Six Sigma process will only produce a few bad parts. The following table lists Sigma Numbers and their Defects per Million (DPM) with a 1.5σ shift.














Sigma (σ)
Defects per Million (DPM)



















2.0
308,300



3.0
67,000



4.0
6,220



5.0
233



6.0
3.4










These numbers are astonishing since many industrial processes are probably running at a Two Sigma level since their quality related process variables are under manual control. That means, these processes are producing 15% to 30% waste in their normal operations. If quality variables are under manual control, quality and efficiency will depend very much on the individual operator's skill, experience, and work attitude. The operator has to be able to tweak all the magic knobs, wait, and hope the final products produced are within the specifications.
1. Quality Control History and Status
Quality control has gone through 4 stages that can be summarized in the following table:












Stage
Period
Description and Tools







1
1900 to 1940
Product Inspection = Quality Control


2
1940 to 1960
Statistical Quality Control (SQC)


3
1960 to 1986
Total Quality Control (TQC)


4
1987 to Present
(i) Online SQC using computer systems




(ii) CSPC = Conventional Process Control +




Statistical Process Control









In the beginning of the 20th century, American engineer T. W. Taylor proposed the idea of product inspection in his scientific management theory. Products are tested against their specifications to pass the good products and reject the bad ones. Although a very large percentage of products today are still made in this way, the technique has major shortcomings including: (1) It does not address the concept of quality; (2) Defective parts are rejected not prevented.
In the early 1940s, the United States required a large volume of high-quality goods for World War II efforts. The government hired a group of scientists to implement a series of quality control standards, thereby forcing the industry to Statistical Quality Control (SQC) methods proposed by Dr. Shewhart and others. SQC is a quality control method based on statistics, which can distinguish the common causes and special causes of quality variations so that quality inspections can be simplified and quality problems can be prevented. [1]
In the 1960s, Japan rebuilt its industry with the help of SQC introduced by Dr. Deming, an America quality control legend. Japan not only applied SQC to its production, but also improved it with a new name, Total Quality Control. TQC embodies a number of new concepts including: (1) QC is everyone's job and (2) Apply QC to all areas possible. [2]
During the 1980s, the huge success of the Japanese quality-based business model forced American industries to learn and apply TQC, and the other quality control and management methods including the well-known Deming's fourteen point management philosophy. [3] In the late 1980s, personal computers made online SQC applications popular and helped industries to improve product quality tremendously, especially in the discrete manufacturing industry. [4]
Unfortunately, however, Statistical Quality Control methods are not sufficient for “true quality control.” Today, a large percentage of processes are still running in 2 to 3 Sigma conditions even though online SQC systems are used. This is because although SQC/SPC tools can identify a process's unwanted variations, it cannot tell how to eliminate or reduce them. In practice, there can be multiple reasons for unwanted variations including (1) human errors, (2) lack of process capability, and (3) poor control of process variables.
In the early 1990s, the author of this patent application observed this reality and proposed the idea of CSPC, which states that quality control relies on the combination of Conventional Process Control and Statistical Quality Control. [5][6][7] The idea is very simple. SQC can find the abnormal causes of variations but it is up to the automatic control system to reduce the variations and keep them under control. However, at the time, PID based control systems could not provide adequate control for many quality variables. This led the author to spend years of work seeking more effective control methods and eventually developed the Model-Free Adaptive (MFA) control technology. The first MFA control patent was filed on Oct. 6, 1997 with Ser. No. 08/944,450 and allowed in January 2000.
2. Control Quality Variables Automatically
To conclude, The most effective way to solve quality problems permanently is to control the quality variables through an automatic process control system. Therefore, it is desirable to develop an automatic control system that can control quality variables automatically and force the quality variables to track their setpoints while the process has to go through various production changeovers in batch, recipe, and product size, etc.
The Model-Free Adaptive (MFA) control methodology described in U.S. patent Ser. No. 08/944,450, and patent application Ser. Nos. 09/143,165 and 09/174,156 are able to deal with various complex systems and effectively control quality related process variables. However, we still face the following difficulties when controlling quality variables:

    
    
        a) The process has significant time-delays and the delay time varies due to production speed changes;
        b) The process also has large disturbances due to “wild” product inflow changes, load changes, etc.; and
        c) The quality variables cannot be measured online and the measurement for the quality is based on off-line lab test data or other off-line measurement methods; and
        d) There is a big change in the system dynamics so that a regular MFA controller is unable to provide prompt and adequate control action to meet the control performance criteria.
    
    


To describe the application in more detail, a zinc galvanizing process is studied in the following. A continuous galvanizing process applies a thin surface coating of zinc to steel products and is a critical operating unit in steel sheet production. Steel sheets are popular products used for cans, refrigerators, and automobiles, etc. The thickness of the zinc layer is an important quality variable. Too thin a layer can cause corrosion and damage the product, and too thick a layer will waste too much zinc, a precious metal.
From a control point of view, a continuous galvanizing process has the following behavior: (1) nonlinear, (2) large and varying time delays, (3) multivariable, and (4) frequent production changeovers. According to an industry survey, the thickness of most continuous galvanizing processes is still under manual control resulting in lower efficiency, wasted manpower and materials, and inconsistent product quality.
In this patent application, we introduce a Model-Free Adaptive control system to automatically control the quality variables and quality related process variables.
SUMMARY OF THE INVENTION
The present invention overcomes the above-identified limitations of the prior art by using a specially designed Model-Free Adaptive control system, which includes an Anti-delay MFA controller, a set of Feedforward MFA controllers, and a signal emulator. The MFA quality variable control system is able to automatically control quality variables under the conditions where there are significant varying time delays and disturbances in the process. Because of its unique capability, the control system is useful for building flexible and adaptive production systems, achieving Six Sigma quality control goals, and fulfilling the on demand manufacturing needs in the new e-commerce environment.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a block diagram illustrating a single-variable model-free adaptive quality variable control system according to this invention.
FIG. 2 is a block diagram illustrating a multivariable model-free adaptive quality variable control system according to this invention.
FIG. 3 is a block diagram illustrating a single-variable model-free adaptive quality variable control system with a signal emulator according to this invention.
FIG. 4 is a block diagram illustrating a multivariable model-free adaptive quality variable control system with a set of signal emulators according to this invention.



DESCRIPTION OF THE PREFERRED EMBODIMENTS
A. Single-Variable Model-Free Adaptive Quality Control
As illustrated in FIG. 1, a single variable Model-Free Adaptive quality variable control system is introduced. The system includes a standard MFA controller 10, a single-input-single-output (SISO) process 12, a delay predictor 14, signal adders 16 and 18. In addition, the system consists of a feedforward process 20, a feedforward MFA controller 22, a delay process 24, and a signal adder 26. The signals shown in FIG. 1 are as follows:












r(t)-Setpoint.


y(t)-Process Variable or the Measured Variable.


x(t)-Output of Process Gp.


yf(t)-Output of the Feedforward Process Gf.


yx(t)-Input to the Delay Process Gd, yx(t) = x(t) + yt(t).


yc(t)-Output of the Delay Predictor.


u(t)-Feedback MFA controller output.


ul(t)-Feedforward MFA controller output.


uc(t)-Combined control output, uc(t) = u(t) + ul(t).


d(t)-Disturbance.


e(t)-Error between the Setpoint and Predictor Output, e(t) = r(t) − yc(t).









The functions of elements 10, 12, 14, and 16 have been described in the patent application Ser. No. 08/944,450 and Ser. No. 09/143,165, which are herein incorporated by reference.
As disclosed in more detail in the above patent applications, the MFA controller 10 consists of a learning neural network with a multi-node input layer to which time- delayed functions of an error signal e(t) are applied, a multi-neuron hidden layer which sums the outputs of the input layer nodes individually weighted by weighting factors wij(n), and a single-neuron output layer which sums functions of the outputs of the hidden layer neurons individually weighted by weighting factors hj(n). The output u(t) of the neural network is a function of the output of the output layer neuron.
The special delay predictor 14 is described in the Anti-delay MFA controller section of in the patent application Ser. Nos. 08/944,450 and 09/143,165. The predictor is designed to produce a dynamic signal yc(t) to replace the measured variable y(t) as the feedback signal. Then, the input to controller 10 is calculated through adder 16 as

e(t)=r(t)−yc(t)  (1)

The idea of this design is to produce an e(t) signal for the controller and let it “feel” its control effect without much delay so that it will keep producing proper control signals.
Since the MFA controller in the system has powerful adaptive capability, the delay predictor is designed in a simple first-order-lag-plus-delay (FOLPD) form represented by the following Laplace transfer function: 



  
    
      
        
          
            
              
                
                  Y
                  c
                
                ⁡
                
                  (
                  S
                  )
                
              
              =
              
                
                  Y
                  ⁡
                  
                    (
                    S
                    )
                  
                
                +
                
                  
                    Y
                    p
                  
                  ⁡
                  
                    (
                    S
                    )
                  
                
              
            
          
        
        
          
            
              
                =
                
                  
                    Y
                    ⁡
                    
                      (
                      S
                      )
                    
                  
                  +
                  
                    
                      
                        K
                        ⁡
                        
                          (
                          
                            1
                            -
                            
                              e
                              
                                
                                  -
                                  τ
                                
                                ⁢
                                
                                  
                                
                                ⁢
                                S
                              
                            
                          
                          )
                        
                      
                      
                        TS
                        +
                        1
                      
                    
                    ⁢
                    
                      U
                      ⁡
                      
                        (
                        S
                        )
                      
                    
                  
                
              
              ,
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




where Y(S), Yp(S), U(S), and Yc(S) are the Laplace transform of signals y(t), yp(t), u(t) and yc(t), respectively; y(t) is the process variable; yp(t) is the predictive signal; and yc(t) is the output of the predictor; K, T, τ are the parameters of the predictor.

The process static gain can be set as 



  
    
      
        
          K
          =
          
            1
            
              K
              c
            
          
        
        ,
      
    
    
      
        (
        3
        )
      
    
  




where Kc is the MFA controller gain, which is entered by the user.

The predictor time constant can be selected as

T=Tc  (4)

where Tc is the estimated process time constant.

The process delay time τ is set based on a rough estimation of process delay time provided by the user.
1. Feedback MFA Controller Parameters
Although the above referenced Anti-Delay MFA controller is simple and works for many processes with large time delays, it lacks the flexibility of adjusting control performance since the controller gain Kc and Time Constant Tc are also tied to the Delay Predictor. For instance, if we increase the Controller Gain Kc for more active control action, since the estimated process static gain K is calculated based on Kc according to Equation (3), the change will affect the output of the Delay Predictor yc(t). Since e(t) is calculated with yc(t) by Equation (1), it will also change and then affect the controller output. This circular relationship forbids us to freely adjust MFA controller parameters.
One way to achieve flexibility is to allow free entry for controller parameters Kc and Tc, as well as estimated process constants K, T, and τ. That means, the user needs to enter 5 parameters when configuring an Anti-Delay MFA controller.
2. Performance Index
In order to simplify the controller configuration procedure, we can design a Performance Index as follows:


0.01≦Ip≦100,  (5)



  
    
      
        
          
            K
            c
          
          =
          
            Ip
            K
          
        
        ,
      
    
    
      
        (
        6
        )
      
    
  






  
    
      
        
          
            T
            c
          
          =
          
            T
            Ip
          
        
        ,
      
    
    
      
        (
        7
        )
      
    
  




where Ip is a real number as the Performance Index with default value of 1.0, Kc is the Controller Gain and Tc is the Time Constant, K and T are estimated process parameters entered by the user. In this way, the user can simply enter K, T, τ to estimate the process dynamics and use Ip to adjust the control performance. Notice that the estimated delay time τ is not related to the Performance Index.

3. Feedforward MFA Controller Parameters

The patent application Ser. No. 09/143,165 disclosed two Feedforward MFA controller designs based on a first-order dynamic block. In order to compensate for the large time delays included in the feedforward process Gf, we introduce a new design for the Feedforward MFA controller in the following Laplace transfer function: 



  
    
      
        
          
            
              G
              fc
            
            ⁡
            
              (
              S
              )
            
          
          =
          
            
              
                K
                fc
              
              ⁢
              
                e
                
                  
                    -
                    
                      τ
                      fc
                    
                  
                  ⁢
                  S
                
              
            
            
              
                
                  T
                  fc
                
                ⁢
                S
              
              +
              1
            
          
        
        ,
      
    
    
      
        (
        8
        )
      
    
  




where its parameters are defined in the following:

    
    
        Kfc—Feedforward Controller Gain,
        Tfc—Feedforward Time Constant.
        Tfc—Feedforward Delay Time.

The configuration procedures for these parameters are described in the following section.

4. Feedforward MFA Control for Process with Large Time Delays

    
    


A process with large time delays can also have significant yet measurable disturbances. In this case, we can add a feedforward MFA controller to reduce the effect of the disturbance to the loop before the feedback loop takes corrective action, thereby improving the control system performance quite economically.
However, applying a feedforward controller to a process with large time delays can be quite complicated. As shown in FIG. 1, the potential time delays can be in the processes Gp, Gf, and a separate process Gd with pure delay time τd.
Typically, for a quality variable that cannot be measured online, delay time τd is the time in between each manual sample or lab test of the quality measurement. In practical applications, this delay tine can be as big as 4 to 8 hours.
In order to describe the process with time delays in more detail, let us define the following symbols for the different delay times:

    
    
        τp—Delay Time in Process Gp 
        τf—Delay Time in Feedforward Gf 
        Σd—Delay Time in Delay-time Process Gd 

If the process has a large delay time, the MFA control system should be configured according to the following cases:

    
    
    Case 1. τd is much bigger than τp and τf.

Solution 1.

    1. For the Anti-delay MFA controller, configure K, T, and Ip as described in the previous section. Estimate τd and let

τ=τd.  (9)

    2. For the Feedforward MFA controller, ignore τp and τf by letting

τfc=0.  (10)

Configure feed forward MFA controller gain based on the following Equation: 



  
    
      
        
          
            K
            fc
          
          =
          
            -
            
              
                K
                f
              
              
                K
                p
              
            
          
        
        ,
      
    
    
      
        (
        11
        )
      
    
  




where Kfc is the MFA Feedforward controller gain, Kf is the estimated static gain of process Gf, and Kp is the estimated static gain of process Gp. Then estimate Tf and let

Tfc=Tf,  (12)

where Tf is the feedforward time constant of feedforward process Gf.

    Case 2. τd is approximately zero, τp is bigger than τf.

Solution 2: Since τd is 0, the disturbance signal from d(t) will show up quickly at y(t). Then the feedback controller will get this signal with almost no delay. Therefore, there is no need to include feedforward control in this case. The feedback Anti-delay controller can be used as usual.

    Case 3. τd is approximately zero, τf is bigger than or equal to τp.

Solution 3. The feedforward MFA controller is required to compensate for the large delay time in Gf based on the following formula:

τfc=τf−τp,  (13)

where τfc is the delay time of the feedforward MFA controller. That is, the feedforward control action of uf needs to be delayed by τfc in order to cancel out the feedforward disturbance of yf to the system.

B. MIMO Model-Free Adaptive Quality Variable Control



Without losing generality, FIG. 2 illustrates a multi-input-multi-output (MIMO) Model-Free Adaptive Quality Variable Control System, which includes a 2-input-2-output MFA Controller 30, a 2-input-2-output (2×2) process 32, two Delay Predictors 34, and two pair of signal adders 36 and 38. Similar to the single variable case, the system also consists of two Feedforward processes 40, two Feedforward MFA controllers 42, two Delay processes 44, and one pair of signal adders 46.
The 2×2 process consists of two main processes G11 and G22, and two sub- processes G21 and G12. The measured Process Variables y1 and y2 are used as the feedback signals of the main control loops. They are fed into the Delay Predictors 34 along with the controller outputs u1 and u2 to generate the synthesized signals yc1(t) and yc2(t). These two signals are compared with the setpoints r1 and r2 at adders 36 to produce errors e1 and e2, which are the inputs to the 2×2 MFA Controller.
The output of each sub-process is cross-added by adders 48 to produce process outputs x1 and x2. Notice that in real applications the outputs from the sub-processes are not measurable and only their combined signals x1 and x2 can be measured. Thus, by the nature of the 2×2 process, the inputs uc1 and uc2 to the process are interconnected with its outputs x1 and x2. One input change wilt cause both outputs to change.
The signals shown in FIG. 2 are as follows:

    
    
        r1(t), r2(t)—Setpoint of controllers C1 and C2, respectively.
        y1(t), y2(t)—Measured Process Variables of the 2×2 process.
        d1(t), d2(t)—Disturbance to y1 and y2, respectively.
        x1(t) and x2(t)—Output of the 2×2 Process.
        yf1(t) and y12(t)—Output of the Feedforward Processes G11 and G12.
        yx1(t) and yx2(t)—Sum of x1(t) and yf1(t), and sum of x2(t) and y12(t), respectively.
        yc1(t) and yc2(t)—Output of the Delay Predictors D1 and D2.
        u1(t) and u2(t)—Output of the feedback MFA controllers C1 and C2.
        uf1(t) and u12(t)—Output of the feedforward MFA controllers Gfc1 and GfC2.
        uc1(t) and uc2(t)—Combined control outputs.
        e1(t) and e2(t)—Error between the Setpoint and Predictor Output for C1 and C2.
    
    


The functions of elements 30, 32, 34, and 36 have been described in the patent application Ser. No. 08/944,450 and Ser. No. 09/143,165. However, a multivariable control system for controlling a multivariable process with large time delays and disturbances is much more complex than its single variable case.
The first important thing in configuring such a control system is the variable pairing between process inputs and outputs. Process inputs are the so-called manipulated variables since they are manipulated by the controller outputs. Process outputs are the controlled Process Variables. The variable pairing rules are described as follows:

    1. Each process in the main loops has to be controllable, open-loop stable, and either reverse or direct acting.
    2. A process with a large static gain should be included in the main loop as the main process; and a process with a small static gain should be treated as the sub-process.
    3. A faster process should be paired as the main process; and a slower process and processes with time delays should be treated as the sub-processes.


A MIMO Anti-Delay MFA controller is required if at least one main process in a multivariable system has significant time delays. According to the pairing rules, we usually do not want to include a process with large time delays in a main loop. However, since the Pairing Rules 2 and 3 may be in conflict, a tradeoff may be the only option.
For the MIMO MFA Quality Variable control system, the design aspects of the feedback MFA controllers and feedforward MFA controllers are the same as described in the single variable case. The key is to focus on each individual main loop after the pairing and interaction problems are handled.
C. SISO MFA Quality Variable Control System with a Signal Emulator
In real applications, we often face the reality that the quality variables are not measured online and the measurement for the quality is based on off-line lab test data or other off-line measurement methods. It is a major challenge to control the quality variables in a closed-loop fashion.
FIG. 3 illustrates an MFA Quality Variable Control System with a Signal Emulator. It consists of elements 10 to 26 as described in Section A, and a Signal Emulator 28, which is designed to provide a smooth and adequate measurement signal of the Process Variable (PV) to the controller under the following conditions: (1) When no online measurement for the Process Variable is available; and (2) The measurement signal is noisy.
When a Signal Emulator is used, the Delay Predictor's input signal becomes yd(t) and its Laplace transfer function is revised into the following form: 



  
    
      
        
          
            
              
                
                  Y
                  c
                
                ⁡
                
                  (
                  S
                  )
                
              
              =
              
                
                  
                    Y
                    d
                  
                  ⁡
                  
                    (
                    S
                    )
                  
                
                +
                
                  
                    Y
                    p
                  
                  ⁡
                  
                    (
                    S
                    )
                  
                
              
            
          
        
        
          
            
              
                =
                
                  
                    
                      Y
                      d
                    
                    ⁡
                    
                      (
                      S
                      )
                    
                  
                  +
                  
                    
                      
                        K
                        ⁡
                        
                          (
                          
                            1
                            -
                            
                              e
                              
                                
                                  -
                                  τ
                                
                                ⁢
                                
                                  
                                
                                ⁢
                                S
                              
                            
                          
                          )
                        
                      
                      
                        TS
                        +
                        1
                      
                    
                    ⁢
                    
                      U
                      ⁡
                      
                        (
                        S
                        )
                      
                    
                  
                
              
              ,
            
          
        
      
    
    
      
        (
        14
        )
      
    
  




where Yd(S), Yp(S), U(S), and Yc(S) are the Laplace transform of signals yd(t), yp(t), u(t) and yc(t), respectively; yd(t) is the output of the Signal Emulator, yp(t) is the predictive signal; yc(t) is the output of the predictor; K, T, τ are the parameters of the predictor. Notice that here the process variable y(t) is replaced by the output of the Signal Emulator yd(t) assuming the signal y(t) is filtered or emulated by the Signal Emulator 28.

There are different ways to build a Signal Emulator depending on the application and required engineering effort. Some of the techniques are described including the following.
1. Soft-Sensor
Developing a Soft-Sensor is probably the most sophisticated approach for building the Signal Emulator. A Soft-Sensor is typically a computer program based process model that can produce an estimated signal to approximate the Process Variable which cannot be measured online. There is no general-purpose Soft-Sensor available because each Process Variable can be so different. There are different techniques for building Soft-Sensors including the first principle model, identification based model, and rule based model, etc.
2. Noise Filter
If the online or offline measurement data include noise, the Signal Emulator can be designed as a noise filter. Depending on the characteristics of the noise, high-pass, low-pass, and other types of filters can be easily implemented.
3. Data Qualifier
IF the measurement is based on off-line tab test data and needs to be entered manually, it is wise to add a mechanism to qualify the data entered in order to reduce or eliminate human error. This is because any bad data entered into the control system can cause major disturbances and must be screened out. A data qualifier can be as simple as a set of high and low limits that can reject any data falling outside the limits. One can also design a more sophisticated algorithm or rule-based mechanism to screen the data based on data patterns, historical data, and a set of rules.
4. Setting of Sample Interval
The simplest way to control a process with only offline measurement data is to run the controllers slowly based on each fresh sampled data. That is, there is no need to calculate the control outputs during the time period when there is no update feedback signal. In this case, the control systems disclosed in Sections A and B are still valid. The only important matter is to select a proper Sample Interval, which is part of the control algorithm in its digital form.
D. MIMO MFA Quality Variable Control System with a Signal Emulators
FIG. 4 illustrates an MIMO MFA Quality Variable Control System with Signal Emulators. It consists of elements 30 to 46 as described in Section B, and Signal Emulators 50. The 2×2 process is represented in a more generic format with main processes Gp1 and Gp2. The design aspects of the Signal Emulators are the same as described for the single variable case.

## Claims

1. A method of adjusting the performance of a model-free adaptive controller having a delay predictor, comprising the steps of:
a) establishing a selectable performance index Ip ranging from substantially 0.01 to substantially 100;
b) setting the controller gain as Kc=Ip/K, where Kc is the controller gain, K is the estimated process static gain, and Ip is the selected performance index; and
c) setting the controller time constant as Tc=T/Ip, where Tc is the controller time constant, T is the estimated time constant of said delay predictor, and Ip is the selected performance index.