US20030139825A1 - Method of ascertaining control parameters for a control system - Google Patents

Method of ascertaining control parameters for a control system Download PDF

Info

Publication number
US20030139825A1
US20030139825A1 US10/323,707 US32370702A US2003139825A1 US 20030139825 A1 US20030139825 A1 US 20030139825A1 US 32370702 A US32370702 A US 32370702A US 2003139825 A1 US2003139825 A1 US 2003139825A1
Authority
US
United States
Prior art keywords
controller
model
plant
performance index
loop
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
US10/323,707
Other versions
US6917840B2 (en
Inventor
Richard Lund
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
MTS Systems Corp
Original Assignee
MTS Systems Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by MTS Systems Corp filed Critical MTS Systems Corp
Priority to US10/323,707 priority Critical patent/US6917840B2/en
Assigned to MTS SYSTEMS CORPORATION reassignment MTS SYSTEMS CORPORATION ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: LUND, RICHARD A.
Publication of US20030139825A1 publication Critical patent/US20030139825A1/en
Application granted granted Critical
Publication of US6917840B2 publication Critical patent/US6917840B2/en
Assigned to JPMORGAN CHASE BANK, N.A., AS ADMINISTRATIVE AGENT reassignment JPMORGAN CHASE BANK, N.A., AS ADMINISTRATIVE AGENT SUPPLEMENTAL CONFIRMATORY GRANT OF SECURITY INTEREST IN UNITED STATES PATENTS Assignors: MTS SYSTEMS CORPORATION
Assigned to MTS SYSTEMS CORPORATION reassignment MTS SYSTEMS CORPORATION RELEASE BY SECURED PARTY (SEE DOCUMENT FOR DETAILS). Assignors: JPMORGAN CHASE BANK, N.A., AS ADMINISTRATIVE AGENT
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B5/00Anti-hunting arrangements
    • G05B5/01Anti-hunting arrangements electric
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Definitions

  • the present invention relates to control of a system machine or process. More particularly, the present invention relates to ascertaining control parameters for a controller to tune a system and achieve acceptable closed-loop control.
  • Control systems are used in a wide variety of fields, such as industrial processes, environmental systems, electronic systems, mechanical systems, and any other system where system output variables representing measurements and user specified desired outputs are processed to generate signals that control devices which change the system.
  • vibration systems which are well known, apply loads and/or motions to test specimens.
  • Vibration systems are widely used for performance evaluation, durability tests, and various other purposes, as they are highly effective in the development of products. For instance, it is quite common in the development of automobiles, motorcycles, or the like, to subject the vehicle or a sub-structure thereof to a laboratory environment that simulates operating conditions such as a road or test track.
  • Physical simulation in the laboratory involves a well-known method of data acquisition and analysis in order to develop drive signals that can be applied to the vibration system to reproduce the operating environment.
  • This method includes instrumenting the vehicle with transducers “remote” to the physical inputs of the operating environment.
  • Common remote transducers include, but are not limited to, strain gauges, accelerometers, and displacement sensors, which implicitly define the operating environment of interest.
  • the vehicle is then driven in the operating environment, while remote transducer responses (internal loads and/or motions) are recorded.
  • actuators of the vibration system are driven so as to reproduce the recorded remote transducer responses on the vehicle in the laboratory.
  • Tuning is the process of setting the controller's internal logic, circuitry and/or variables so that the controller's output signals cause the desired effect given the controller's feedbacks.
  • the overall system can be organized or defined as having two principle components, i.e., the “controller” and the “plant”.
  • the controller receives inputs or otherwise has stored information pertaining to desired operation of the plant.
  • the controller provides actuation signals to the plant and receives therefrom feedback signals, which can include the remote sensors as discussed above pertaining to the desired response of the plant, and/or intermediate feedback signals deemed necessary to achieve the desired response of the plant.
  • the plant comprises all components that are not part of the controller.
  • the plant or physical system would include the servo valves, which receive the actuation signals from the controller that in turn are used to operate actuators to impart forces or motions upon the test specimen.
  • the plant would also include any necessary linkages between the actuators and the test specimen as well as test specimen itself and sensors (intermediate and/or remote) used to provide feedback signals back to the controller.
  • the commissioning of control systems is commonly accomplished by constructing a controller of suitable architecture, connecting it to the plant that is to be controlled, exciting the system while adjusting the control parameters (tune the system) until acceptable closed-loop control is achieved.
  • the tuning process typically requires a significant skill level and is often quite tedious in the case of Single Input Single Output (SISO) systems, and can be overwhelming in the case of Multiple Input Multiple Output (MIMO) systems with a high degree of cross-coupling.
  • a method of ascertaining control parameters for a control system having a controller and a plant includes establishing a model of the plant and a model of the controller and calculating a performance index for a closed-loop system as a function of controller parameters, accounting for selected stability margins.
  • FIG. 1 is a block diagram of an exemplary environment for practicing the present invention.
  • FIG. 2 is a computer for implementing the present invention.
  • FIG. 3 is a block diagram illustrating components of a closed-loop control system for the environment of FIG. 1.
  • FIG. 4 is a block diagram illustrating a tuning assembly and other aspects of the present invention.
  • FIG. 5 is a z-plane contour.
  • FIG. 6 is a s-plane contour.
  • FIG. 7 is a flow diagram of a method for creating a model of the plant.
  • FIG. 8 is a plot of an impulse response function and a tapering function.
  • FIG. 9 is a two-dimensional Nyquist plot of a single channel.
  • FIGS. 10 and 11 are three-dimensional Nyquist plots.
  • FIG. 12 is a plot of magnitude v. frequency for a plurality of channels.
  • FIG. 13 is a flow diagram for a tuning procedure.
  • FIG. 14 is a flow diagram for calculating a performance index.
  • FIG. 15 is a block diagram illustrating the tuning assembly with multiple plant configurations.
  • FIG. 16 is a block diagram illustrating components of a second form of a closed-loop control system for the environment of FIG. 1.
  • FIG. 17 is a block diagram illustrating a cascade control-loop.
  • FIG. 1 illustrates a plant or physical system 10 .
  • the physical system 10 generally includes a vibration system 13 comprising a servo controller 14 and an actuator 15 .
  • the actuator 15 represents one or more actuators that are coupled through a suitable mechanical interface 16 to a test specimen 18 .
  • the servo controller 14 provides an actuator command signal 19 to a servo valve 25 to operate the actuator 15 , which in turn, excites the test specimen 18 .
  • Suitable feedback 15 A can be provided from the actuator 15 to the servo controller 14 or from other sensors.
  • One or more remote transducers 20 on the test specimen 18 such as displacement sensors, strain gauges, accelerometers, or the like, provide a measured or actual response 21 .
  • a physical system controller 23 receives the actual response 21 as feedback in a response to a drive 17 as input to the servo controller 14 .
  • signal 17 is a reference signal
  • signal 19 is a manipulated variable (command to actuated device)
  • signal 15 A is a feedback variable. This relationship is also illustrated in FIG. 3.
  • FIG. 1 for the single channel case, multiple channel embodiments with signal 15 A comprising N feedback components and the signal 19 comprising M manipulated variable components are typical and considered another embodiment of the present invention.
  • the plant includes the manufacturing machines (e.g. presses, molding apparatus, forming machines, etc.) and the manipulated variable signal 19 provides command signals to said machines, and the feedback signal 15 A comprises manual or automatic measured parameters of the plant.
  • the feedback signal 15 A comprises manual or automatic measured parameters of the plant.
  • Another example includes an oil refinery where the plant is the process plant and the feedback signal 15 A comprises intermediate or final parameters related to its operation.
  • FIG. 2 and the related discussion provide a brief, general description of a suitable computing environment in which the invention may be implemented.
  • the servo controller 14 will be described, at least in part, in the general context of computer-executable instructions, such as program modules, being executed by a computer 30 .
  • program modules include routine programs, objects, components, data structures, etc., which perform particular tasks or implement particular abstract data types.
  • the program modules are illustrated below using block diagrams and flowcharts. Those skilled in the art can implement the block diagrams and flowcharts to computer-executable instructions storable on a computer readable medium.
  • the invention may be practiced with other computer system configurations, including multi-processor systems, networked personal computers, mini computers, main frame computers, and the like.
  • the invention may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network.
  • program modules may be located in both local and remote memory storage devices.
  • the computer 30 illustrated in FIG. 2 comprises a conventional personal or desktop computer having a central processing unit (CPU) 32 , memory 34 and a system bus 36 , which couples various system components, including the memory 34 to the CPU 32 .
  • the system bus 36 may be any of several types of bus structures including a memory bus or a memory controller, a peripheral bus, and a local bus using any of a variety of bus architectures.
  • the memory 34 includes read only memory (ROM) and random access memory (RAM).
  • ROM read only memory
  • RAM random access memory
  • Storage devices 38 such as a hard disk, a floppy disk drive, an optical disk drive, etc., are coupled to the system bus 36 and are used for storage of programs and data. It should be appreciated by those skilled in the art that other types of computer readable media that are accessible by a computer, such as magnetic cassettes, flash memory cards, digital video disks, random access memories, read only memories, and the like, may also be used as storage devices. Commonly, programs are loaded into memory 34 from at least one of the storage devices 38 with or without accompanying data.
  • An input device 40 such as a keyboard, pointing device (mouse), or the like, allows the user to provide commands to the computer 30 .
  • a monitor 42 or other type of output device is further connected to the system bus 36 via a suitable interface and provides feedback to the user.
  • the reference signal 17 can be provided as an input to the computer 30 through a communications link, such as a modem, or through the removable media of the storage devices 38 .
  • the manipulated variable signals 19 are provided to the plant 10 of FIG. 1 based on program modules executed by the computer 30 and through a suitable interface 44 coupling the computer 30 to the vibration system 13 .
  • the interface 44 also receives the feedback signals 15 A.
  • the servo controller 14 can also comprise an analog controller with or without digital supervision as is well known.
  • Controller Drive 17 may be continuously provided or stored in servo controller 14 , depending on the capabilities of the servo controller 14 .
  • Functions of controller 23 and controller 14 can be combined into one computer system.
  • controller 14 is a single board computer operable on a network bus of another computer, which could be controller 23 or another supervisory computer.
  • the schematic diagram of FIG. 2 is intended to generally represent these and other suitable computing environments.
  • test vehicle Although described below with respect to a test vehicle, it should be understood that the present invention discussed below is not confined to testing only vehicles, but can be used on other processes, types of test specimens and substructures or components thereof.
  • FIG. 4 illustrates a servo controller 14 and the plant 10 along with a tuning or optimization assembly 60 .
  • the tuning assembly 60 ascertains or determines control parameters or variables of the servo controller 14 based on user constraints 101 .
  • the control parameters can include one or more values that are used generally in a closed-loop control architecture to generate signal 19 in part based on feedback 15 A during actual operation of the plant 10 .
  • the tuning assembly 60 can ascertain the control parameters offline without repeated operation of the plant 10 , thereby, saving time and minimizing wear or damage to the plant 10 .
  • the tuning assembly 60 uses a model 64 of the system controller 14 and a model 70 of the plant 10 , in one embodiment, both in the frequency domain, in order to estimate a performance index as a function of controller parameter settings and user constraints. Optimization of the performance index using an optimization algorithm such as Nelder-Mead yields the best or optimum set of controller parameters for the servo controller 14 .
  • “best” or “optimum” refers to those values of the controller parameters that have been ascertained by the tuning assembly 60 given user constraints. As appreciated by those skilled in the art, further manual adjustments to the ascertained controller parameters may be possibly made by the user to yield some further improvements in system performance. Nevertheless, the ascertained controller parameters provided by the tuning assembly 60 would achieve desired system performance based on the user constraints, and thus, can be considered “best” or optimum”.
  • the tuning assembly 60 includes an optimization module 62 and a performance calculator 63 .
  • the performance calculator 63 uses the model 64 of the system controller 14 and the model 70 of the plant 10 . Before further describing operation of the performance calculator 63 and the optimization module 64 , it may be helpful to describe construction of the servo controller model 64 and plant model 70 .
  • the servo controller 14 can be embodied in a digital computer such as computer 30 illustrated in FIG. 2.
  • the control architecture is embodied as modules or instructions implementing control functions such as summation, subtraction, integration, differentiation and filtering, just to name a few.
  • Signals 19 are computed using the modules or instructions and converted, if necessary, to analog values through digital-to-analog converters as is well known.
  • the reference signal 17 and the feedback signals 15 A if provided in analog form, can be converted to digital values through suitable analog-to-digital converters. Digital values of the reference signal 17 and the feedback signals 15 A are then used by the modules defining the control architecture.
  • control architecture of servo controller 14 will thus be well defined. Constructing a mathematical representation or model of the servo controller 14 architecture in the frequency domain is well known. A simple exemplary controller structure for illustrative purposes can be
  • Gc Km *( Kp+Kd *( z ⁇ 1)/ Tz+Kdd *(( z ⁇ 1)/ Tz ) ⁇ circumflex over ( ) ⁇ 2)
  • Km is a scalar and Kp, Kd, Kdd are vectors of length equal to the number of channels being controlled.
  • This form is for the case of an integrating plant such as a servo-hydraulic system.
  • construction includes substituting discrete values for z according to the plot of FIG. 5 for each of the controller components.
  • linear spacing along the quarter circles 81 and log spacing along the unit circle 83 realizes an accurate model with a minimum number of data points.
  • this form of the model can be referred to as a transfer function (TF).
  • TF transfer function
  • a vector of complex values becomes the frequency domain representation of the corresponding controller element.
  • This model can be arranged in a matrix where such vectors commonly, but not exclusively, are disposed along a matrix diagonal.
  • servo controller 14 is constructed of analog components
  • the servo controller model 64 is formed by substituting discrete values for s according to the plot of FIG. 6 for each of the controller components.
  • This form is for the case of an integrating plant such as a servo-hydraulic system.
  • linear spacing along the quarter circles 85 and log spacing along the imaginary axis 87 in the s-plane can be used.
  • the models based on the s or z contours can be much more efficient for subsequent calculations (10 times fewer data values), and equally significant is that stability analysis from this data is traceable to first principals of complex analysis. Simplification of this approach to the direct use of the FRF representation can be beneficial in some situations. Actual controller structures that need to be accommodated are typically more complex, but the method of substitution of z and s contours is the same.
  • the plant model 70 can be obtained by model creation module 74 from empirical data using method 80 of FIG. 7.
  • servo controller 14 of FIG. 1 is operated to generate manipulated variable signal 19 as a function of a drive 17 having wide spectrum content.
  • drive 17 can include transients, and/or repeating periodic signals having random phase.
  • a suitable method for formation of drive 17 is described in U.S. Pat. No. 6,385,564, incorporated herein by reference.
  • the manipulated variable signal 19 and the corresponding feedback signal 15 during operation of the plant 10 is recorded, which is schematically illustrated at 84 in FIG. 4.
  • the FRF frequency response function
  • the inverse transform of the FRF is computed to obtain the impulse response function at step 88 .
  • a fundamental attribute of the data collected in storage 84 is to provide one estimate of the plant 10 dynamics with sufficient fidelity and certainty. Further, storage 84 represents collection of such data, whether actually stored or used in real time. Thus, depending on the capabilities of the computer executing the tuning assembly 60 and/or controller 14 , it is possible to ascertain control parameters “online” (when the plant 10 is operating) or “offline” without further operation of the plant 10 .
  • the inverse response functions can be tapered at step 89 to increase model fidelity.
  • FIG. 8 illustrates a representative impulse response function at 90 .
  • a tapering function 92 can be applied to the impulse function 90 to remove content known to be noise.
  • Characteristics of the tapering function 92 i.e. duration and shape of the taper
  • step 96 perform forward transformation of the impulse response functions (possibly, tapered) at discrete values along the complex z domain contour (FIG. 5), or the complex s domain contour (FIG. 6), thereby producing a frequency domain model 70 consistent with the servo controller model 64 , described above.
  • This process minimizes the number of data points and produces a model that is consistent with the servo controller model 64 .
  • elements of the frequency domain model 70 can be normalized in order to ease interpretation of values used during optimization discussed below. For instance, an element representing displacement feedback can be normalized to correspond to 1/s at low frequencies.
  • the plant model 70 and the normalization factors are. stored.
  • the tuning assembly 60 can be implemented on a suitable computer or computers, for instance, the computer 30 discussed above. Generally, as discussed above, the tuning assembly 60 uses a model of the system controller 14 and a model of the plant 10 , both in the frequency domain, in order to estimate a performance index as a function of controller parameter settings and user constraints. Optimization of the performance index using an optimization algorithm such as Nelder-Mead yields the best set of controller parameters for the servo controller 14 .
  • an optimization algorithm such as Nelder-Mead yields the best set of controller parameters for the servo controller 14 .
  • the closed-loop (cl) frequency response can be calculated as:
  • Gcl [I+Km*H*Gp*Gc] ⁇ 1 *Km*Gp*Gc (1)
  • stability can be determined by analysis of the determinant, represented in this case as a scalar valued function of complex frequency.
  • a Nyquist plot of the determinant can thus be used to ascertain stability (note encirclement of the origin in this case). However, this does not produce a measure of stability margin.
  • a method for establishing a stability margin is to decompose the term H*Gp*Gc (the open-loop TF with Km factored out) using a Schur decomposition to obtain:
  • GHt is a triangular matrix of functions of frequency
  • U is a unitary matrix of functions
  • denotes conjugate transposition
  • GHd contains the diagonal elements of GHt (eigen-functions).
  • the tuning assembly 60 optimizes the controller parameters based in part on user constraints, which in the block diagram of FIG. 4 is illustrated at 101 .
  • user constraints can include constraints on independent variables, i.e., the parameters of the controller and constraints on system performance (e.g. stability margin and peaking of the closed-loop frequency response).
  • FIG. 9 illustrates a convenient representation for analyzing stability margin.
  • the circle 102 around “ ⁇ 1,0” represents the stability margin boundary user constraint.
  • the shape is user adjustable. Although a constant radius is shown, other shapes of the stability margin boundary constraint can be used.
  • a single channel Eigen-value trajectory is plotted.
  • Nyquist stability can be plotted in three dimensions with frequency as the third dimension and the stability margin is illustrated as a tube or column 104 , herein having a constant radius centered on “ ⁇ 1,0”.
  • FIG. 10 is an example of such an illustration.
  • FIG. 10 illustrates a multi-channel example. With this form of illustration, a user can easily determine the important frequencies affecting stability as viewed by plotting the Eigen-value trajectories.
  • the stability margin boundary constraint can be allowed to vary as a function of frequency.
  • FIG. 11 illustrates a stability margin boundary (herein a tube 106 having portions of different radii centered at varying positions on the real-imaginary plane) allowing increased stability margin in regions of greater uncertainty such as at higher frequencies.
  • viewing of the three-dimensional renderings of FIGS. 10 and 11 can be changed as desired to allow the trajectories to be readily seen with respect to the stability margin.
  • the shape of the stability margin boundary 106 can be adjusted, predetermined shapes can be stored and selected by the user.
  • FIG. 12 graphically illustrates an absolute upper limit 108 throughout the entire frequency spectrum.
  • FIG. 12 also illustrates a “tuning bandwidth” 110 that is user selectable. The tuning bandwidth defines the frequency region of primary interest for closed-loop performance.
  • a tuning procedure is generally illustrated in FIG. 13 at 130 .
  • the controller model 64 and the plant model 70 are obtained as discussed above.
  • the user can select all or a subset of channels (manipulated variables and feedback responses) to tune or optimize.
  • channels manipulated variables and feedback responses
  • cross-coupling between channels can be more prevalent for some channels than for others.
  • loads and motions can be applied to each of the vehicle spindles from separate loading assemblies coupled to each of the vehicle spindles.
  • Each loading assembly comprises separate actuators and linkages that are operated to apply forces and motions with respect to an orthogonal set of axes, generally with one axis corresponding to the spindle axis.
  • Each loading assembly includes a set of manipulated variable signals for controlling the servo valves of the actuators and a set of corresponding feedback signals.
  • each vehicle spindle includes a loading assembly
  • four sets of channels are provided to the servo controller.
  • cross-coupling may be minimal between some sets of loading assemblies, or even some channels in any one loading assembly. For instance, cross-coupling is generally low between the front and rear vehicle spindle loading assemblies.
  • Step 133 allows the user to select which channels (a subset or all) to tune or optimize. In this manner, if desired, the user can individually select single input and single output channels to optimize, or select a subset of channels to optimize. Generally, it is faster to optimize a plurality of subsets of channels or portions of the plant, than to optimize the complete plant at once. By processing the plant in subsets, some cross-coupling is ignored, but as stated above, this may be acceptable given the configuration of the plant.
  • plant diagnostics can be provided. For instance, it may be apparent from the plant model or from the optimization procedure that channels are totally dysfunctional or poorly performing. In this manner, the optimization or tuning procedure can be used to diagnose problems in the plant. This is particularly true when using single channel tuning repeated over all the channels since this determines the adequacy of performance of each channel operating without interaction with other channels.
  • user constraint values include stability margins, if desired, as a function of frequency, a frequency domain peaking limit and the tuning bandwidth.
  • the gains in the controller model 64 are normalized by frequency by a nominal frequency of operation.
  • the user constraint values can also include an expected system bandwidth 110 (FIG. 12).
  • step 136 initial values of the tuning parameters are provided. If the tuning parameters were normalized, each of the initial values can be set to their nominal value of one.
  • the optimization routine such as Nelder-Mead, is then invoked at step 138 .
  • the optimization routine repeatedly provides a set of controller parameters 140 to the performance calculator 63 , and a performance index 142 is returned from the performance calculator 63 until an optimum value is achieved (i.e. within a selected tolerance range).
  • the optimization routines are readily available and well known, such as in the routines available from Mathworks, Inc. of Natick, Mass., U.S.A.
  • a set of tuning parameters (independent variables) is received by the performance calculator 63 from the optimizer module 62 .
  • the tuning parameters comprise the derivative gains, the second derivative gains (if present), various feedback and/or forward loop filter settings.
  • a function of the proportional gains can be included.
  • the proportional gains can be incorporated into the performance index, i.e., maximizing the proportional gain can be a desired parameter for a single channel system.
  • the open-loop TF model (H*Gp*Gc) is computed as a function of the tuning parameters accounting for any constraints on the independent variables including filter settings. It should be noted that such constraints can also be conveniently handled by a optimization routine, which handles constraints directly.
  • the open-loop TF model is decomposed frequency-by-frequency using the Schur decomposition that produces the Eigen-functions. If desired, Eigen value decomposition can be used.
  • the master gain Km is adjusted (typically an iterative algorithm) at step 158 such that the Eigen-function trajectories contact the stability margin region (nominally a region, for example, a circle, ellipse, etc. that includes ( ⁇ 1,0) of the complex plane).
  • the stability margin region typically a region, for example, a circle, ellipse, etc. that includes ( ⁇ 1,0) of the complex plane.
  • the closed-loop plant response is computed as [I+Km*H*Gp*Gc] ⁇ 1 *Km*Gp*Gc. If peaking of the closed-loop response exceeds a specified maximum value, lower Km until the specified maximum value is achieved (typically an iterative algorithm). A property of Km is that the stability margin is increased and so the constraint on stability is still met.
  • the plant can comprise multiple configurations indicated at 10 , 10 A and 10 B.
  • the different configurations 10 , 10 A and 10 B can be for different test specimens or no specimen, plant operating levels, or other variations of plant operations.
  • a corresponding plant model 70 , 70 A and 70 B would be obtained.
  • the tuning assembly 60 can ascertain the best set of controller parameters that would be useable over all of the various plant configurations. In particular, one set of tuning parameters would be provided at step 152 , but each of the steps 154 , 156 , 158 and 160 would be performed for each of the plant configurations.
  • step 158 This entails at step 158 to adjust the master gain Km until each plant configuration satisfies the stability margin for all system model configurations.
  • the master gain Km is reduced further, if necessary, so that the peaking limit is met across all channels for all plant configurations.
  • the peaking limit can be a function of each channel and/or a function of frequency.
  • visual renderings can be provided to the user to show stability and system performance using the optimized parameters.
  • the renderings comprising another aspect of the present invention, can include any of the three-dimensional Nyquist plots described above. By plotting the Eigen-value trajectories, the user can readily confirm that system stability has been been met because the trajectories will not visually intersect with the selected stability margins. Likewise frequency peaking can be visually rendered such as illustrated in FIG. 12.
  • step 172 visual renderings similar to step 170 for system stability and system performance can be provided to the user for the complete system. If any of the Eigen-value trajectories visually intersect with the selected stability margins, system stability may not be acceptable. Likewise frequency peaking can be checked by visually rendering system performance such as illustrated in FIG. 12. It should be noted visual rendering is not required, but is a convenient form for interpreting system stability and/or system performance. Of course, other forms of indications (e.g. numerical tables, audible alarms, etc.) can be provided to the user to convey violations of system model stability or system model performance.
  • indications e.g. numerical tables, audible alarms, etc.
  • control parameters for the controller 14 are calculated as a function of the optimum set of controller parameters determined by the optimization routine.
  • the “performance index” is a measure of “goodness” of the servo-system. This generally is a measure of control bandwidth. There is a convenient relationship between controller parameters and control bandwidth when using a normalized plant (scaled to 1/s), that is:
  • Kp can be treated as having units of rad/sec, the loop bandwidth.
  • Kp is a vector of length equal to the number of channels (Km is always a scalar) . If Kp is a vector of equal constants (if normalized, i.e. the proportional gains being of a selected relationship among each other), then all channels will have equal loop bandwidth with the loop bandwidth value determined by Km*Kp. Thus maximizing Km over the tuning parameters (derivative gains, filter settings, etc.) makes sense. All channels operating with the same loop bandwidth is highly desirable in systems such as multi-axis tracking systems. In this case, system performance is limited by the weakest channel by necessity.
  • the weakest channel in a system is determined by individual loop dynamics as well as cross-coupling to other channels.
  • cross-coupling When cross-coupling is present, channels are sharing stability in some sense. Thus increasing loop-bandwidth of one channel can require that loop-bandwidth be lowered on one or more other channels.
  • Distributing loop-bandwidth amongst channels can be accomplished be setting the elements of the Kp vector to different values before optimizing Km.
  • the channel least able to make its relative loop-bandwidth will produce the value of Km at the optimum.
  • Performance ⁇ ⁇ Index ⁇ i ⁇ ⁇ ( K m ⁇ K p , ⁇ i ⁇ ref ) x ( 7 )
  • Another useful variation of the performance index is to include a factor that puts more focus on a particular frequency range.
  • An example is the minimum (graphically illustrated at 109 in FIG. 12) of the closed-loop frequency response over a user specified frequency range. This can be applied to each individual channel (diagonal elements) or taken as the minimum over all diagonal elements.
  • Another variation of the performance index is to apply Eq. 7 to preselected groups of channels rather than to each channel individually.
  • Each term in the sum would be a measure of the performance of the group to be balanced against the group. This can be especially useful in the case of symmetry, side-to-side symmetry of channels of the plant, for example.
  • the corresponding left-side and the right-side channels can be treated as a pair and form a group in the overall process.
  • FIG. 16 A wide variety of multivariable servo system models can be represented using the diagram illustrated in FIG. 3. However in other multivariable systems an alternative embodiment of FIG. 16 can be beneficial.
  • the plant 10 is represented as a combination of a plurality (e.g. four) controller blocks.
  • each of the models G1, G2, G3 and G4 represents a model of each of the controlled variables as a function the manipulated variable all in a multivariable system sense.
  • the outputs of G1-G4 in a vibration system can correspond to the controlled variable, displacement feedback, velocity feedback and acceleration feedback, respectively.
  • Kp, Kd and Kdd are adjustable vector of gains for G2, G3 and G4, respectively.
  • Km is the master gain similar to that discussed above.
  • FIG. 16 This configuration of FIG. 16 can be tuned using the tuning assembly 60 in an optimal manner similar to FIG. 3, where the closed loop frequency response can be represented as:
  • Gcl G 1 *[I+Km* ( KpG 2 +KdG 3 +KddG 4)] ⁇ 1 KmKp (9)
  • the present invention can also be applied to servo control configurations involving a loop within a loop.
  • One example is the well-known Cascade configuration having an inner loop 180 and an outer loop 182 , the outer loop 182 including the inner loop 180 as illustrated in FIG. 16.
  • the tuning procedure can be first applied to the inner loop 180 . With the control parameters of inner loop 180 ascertained, the inner loop 180 becomes part of the plant for the outer loop 182 , at which point the tuning procedure can then be applied to the outer loop 182 to ascertain the control parameters thereof.
  • the tuning or optimization procedure of the outer loop 182 can include as a portion of its performance index calculation, the tuning procedure of the inner loop 180 .

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Evolutionary Computation (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Feedback Control In General (AREA)
  • Control Of Eletrric Generators (AREA)
  • Control Of Steam Boilers And Waste-Gas Boilers (AREA)

Abstract

A method of ascertaining control parameters for a control system having a controller and a plant includes establishing a model of the plant and a model of the controller and calculating a performance index for a closed-loop system as a function of controller parameters, accounting for selected stability margins.

Description

    CROSS-REFERENCE TO RELATED APPLICATION
  • This application claims the benefit of U.S. [0001] Provisional Patent Application 60/341,887, filed Dec. 18, 2001, the content of which is hereby incorporated by reference in its entirety.
  • BACKGROUND OF THE INVENTION
  • The present invention relates to control of a system machine or process. More particularly, the present invention relates to ascertaining control parameters for a controller to tune a system and achieve acceptable closed-loop control. [0002]
  • Control systems are used in a wide variety of fields, such as industrial processes, environmental systems, electronic systems, mechanical systems, and any other system where system output variables representing measurements and user specified desired outputs are processed to generate signals that control devices which change the system. For example, vibration systems, which are well known, apply loads and/or motions to test specimens. Vibration systems are widely used for performance evaluation, durability tests, and various other purposes, as they are highly effective in the development of products. For instance, it is quite common in the development of automobiles, motorcycles, or the like, to subject the vehicle or a sub-structure thereof to a laboratory environment that simulates operating conditions such as a road or test track. [0003]
  • Physical simulation in the laboratory involves a well-known method of data acquisition and analysis in order to develop drive signals that can be applied to the vibration system to reproduce the operating environment. This method includes instrumenting the vehicle with transducers “remote” to the physical inputs of the operating environment. Common remote transducers include, but are not limited to, strain gauges, accelerometers, and displacement sensors, which implicitly define the operating environment of interest. The vehicle is then driven in the operating environment, while remote transducer responses (internal loads and/or motions) are recorded. During simulation with the vehicle mounted to the vibration system, actuators of the vibration system are driven so as to reproduce the recorded remote transducer responses on the vehicle in the laboratory. [0004]
  • However, before simulation testing can occur, the control system must be tuned. Tuning is the process of setting the controller's internal logic, circuitry and/or variables so that the controller's output signals cause the desired effect given the controller's feedbacks. [0005]
  • Generally, the overall system can be organized or defined as having two principle components, i.e., the “controller” and the “plant”. Operating in closed-loop control, the controller receives inputs or otherwise has stored information pertaining to desired operation of the plant. The controller provides actuation signals to the plant and receives therefrom feedback signals, which can include the remote sensors as discussed above pertaining to the desired response of the plant, and/or intermediate feedback signals deemed necessary to achieve the desired response of the plant. The plant comprises all components that are not part of the controller. Using the vibration system discussed above by way of example, the plant or physical system would include the servo valves, which receive the actuation signals from the controller that in turn are used to operate actuators to impart forces or motions upon the test specimen. The plant would also include any necessary linkages between the actuators and the test specimen as well as test specimen itself and sensors (intermediate and/or remote) used to provide feedback signals back to the controller. [0006]
  • The commissioning of control systems is commonly accomplished by constructing a controller of suitable architecture, connecting it to the plant that is to be controlled, exciting the system while adjusting the control parameters (tune the system) until acceptable closed-loop control is achieved. The tuning process typically requires a significant skill level and is often quite tedious in the case of Single Input Single Output (SISO) systems, and can be overwhelming in the case of Multiple Input Multiple Output (MIMO) systems with a high degree of cross-coupling. [0007]
  • Accordingly, there is a continuing need for an improved method of tuning control systems. A method or system that is easy to use and requires minimal actual plant operation would be very beneficial. [0008]
  • SUMMARY OF THE INVENTION
  • A method of ascertaining control parameters for a control system having a controller and a plant includes establishing a model of the plant and a model of the controller and calculating a performance index for a closed-loop system as a function of controller parameters, accounting for selected stability margins.[0009]
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • FIG. 1 is a block diagram of an exemplary environment for practicing the present invention. [0010]
  • FIG. 2 is a computer for implementing the present invention. [0011]
  • FIG. 3 is a block diagram illustrating components of a closed-loop control system for the environment of FIG. 1. [0012]
  • FIG. 4 is a block diagram illustrating a tuning assembly and other aspects of the present invention. [0013]
  • FIG. 5 is a z-plane contour. [0014]
  • FIG. 6 is a s-plane contour. [0015]
  • FIG. 7 is a flow diagram of a method for creating a model of the plant. [0016]
  • FIG. 8 is a plot of an impulse response function and a tapering function. [0017]
  • FIG. 9 is a two-dimensional Nyquist plot of a single channel. [0018]
  • FIGS. 10 and 11 are three-dimensional Nyquist plots. [0019]
  • FIG. 12 is a plot of magnitude v. frequency for a plurality of channels. [0020]
  • FIG. 13 is a flow diagram for a tuning procedure. [0021]
  • FIG. 14 is a flow diagram for calculating a performance index. [0022]
  • FIG. 15 is a block diagram illustrating the tuning assembly with multiple plant configurations. [0023]
  • FIG. 16 is a block diagram illustrating components of a second form of a closed-loop control system for the environment of FIG. 1. [0024]
  • FIG. 17 is a block diagram illustrating a cascade control-loop.[0025]
  • DETAILED DESCRIPTION OF THE ILLUSTRATIVE EMBODIMENT
  • FIG. 1 illustrates a plant or [0026] physical system 10. In the exemplary embodiment, the physical system 10 generally includes a vibration system 13 comprising a servo controller 14 and an actuator 15. In the schematic illustration of FIG. 1, the actuator 15 represents one or more actuators that are coupled through a suitable mechanical interface 16 to a test specimen 18. The servo controller 14 provides an actuator command signal 19 to a servo valve 25 to operate the actuator 15, which in turn, excites the test specimen 18. Suitable feedback 15A can be provided from the actuator 15 to the servo controller 14 or from other sensors. One or more remote transducers 20 on the test specimen 18, such as displacement sensors, strain gauges, accelerometers, or the like, provide a measured or actual response 21. A physical system controller 23 receives the actual response 21 as feedback in a response to a drive 17 as input to the servo controller 14. In the illustration of FIG. 1, signal 17 is a reference signal, signal 19 is a manipulated variable (command to actuated device) and signal 15A is a feedback variable. This relationship is also illustrated in FIG. 3. Although illustrated in FIG. 1 for the single channel case, multiple channel embodiments with signal 15A comprising N feedback components and the signal 19 comprising M manipulated variable components are typical and considered another embodiment of the present invention.
  • Although described herein where the environment comprises the [0027] vibration system 13, aspects of the present invention described below can be applied to other plants. For instance, in a manufacturing process, the plant includes the manufacturing machines (e.g. presses, molding apparatus, forming machines, etc.) and the manipulated variable signal 19 provides command signals to said machines, and the feedback signal 15A comprises manual or automatic measured parameters of the plant. Another example includes an oil refinery where the plant is the process plant and the feedback signal 15A comprises intermediate or final parameters related to its operation.
  • FIG. 2 and the related discussion provide a brief, general description of a suitable computing environment in which the invention may be implemented. Although not required, the [0028] servo controller 14 will be described, at least in part, in the general context of computer-executable instructions, such as program modules, being executed by a computer 30. Generally, program modules include routine programs, objects, components, data structures, etc., which perform particular tasks or implement particular abstract data types. The program modules are illustrated below using block diagrams and flowcharts. Those skilled in the art can implement the block diagrams and flowcharts to computer-executable instructions storable on a computer readable medium. Moreover, those skilled in the art will appreciate that the invention may be practiced with other computer system configurations, including multi-processor systems, networked personal computers, mini computers, main frame computers, and the like. The invention may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computer environment, program modules may be located in both local and remote memory storage devices.
  • The [0029] computer 30 illustrated in FIG. 2 comprises a conventional personal or desktop computer having a central processing unit (CPU) 32, memory 34 and a system bus 36, which couples various system components, including the memory 34 to the CPU 32. The system bus 36 may be any of several types of bus structures including a memory bus or a memory controller, a peripheral bus, and a local bus using any of a variety of bus architectures. The memory 34 includes read only memory (ROM) and random access memory (RAM). A basic input/output (BIOS) containing the basic routine that helps to transfer information between elements within the computer 30, such as during start-up, is stored in ROM. Storage devices 38, such as a hard disk, a floppy disk drive, an optical disk drive, etc., are coupled to the system bus 36 and are used for storage of programs and data. It should be appreciated by those skilled in the art that other types of computer readable media that are accessible by a computer, such as magnetic cassettes, flash memory cards, digital video disks, random access memories, read only memories, and the like, may also be used as storage devices. Commonly, programs are loaded into memory 34 from at least one of the storage devices 38 with or without accompanying data.
  • An [0030] input device 40 such as a keyboard, pointing device (mouse), or the like, allows the user to provide commands to the computer 30. A monitor 42 or other type of output device is further connected to the system bus 36 via a suitable interface and provides feedback to the user. The reference signal 17 can be provided as an input to the computer 30 through a communications link, such as a modem, or through the removable media of the storage devices 38. The manipulated variable signals 19 are provided to the plant 10 of FIG. 1 based on program modules executed by the computer 30 and through a suitable interface 44 coupling the computer 30 to the vibration system 13. The interface 44 also receives the feedback signals 15A. Nevertheless, the servo controller 14 can also comprise an analog controller with or without digital supervision as is well known. Drive 17 may be continuously provided or stored in servo controller 14, depending on the capabilities of the servo controller 14. Functions of controller 23 and controller 14 can be combined into one computer system. In another computing environment, controller 14 is a single board computer operable on a network bus of another computer, which could be controller 23 or another supervisory computer. The schematic diagram of FIG. 2 is intended to generally represent these and other suitable computing environments.
  • Although described below with respect to a test vehicle, it should be understood that the present invention discussed below is not confined to testing only vehicles, but can be used on other processes, types of test specimens and substructures or components thereof. [0031]
  • FIG. 4 illustrates a [0032] servo controller 14 and the plant 10 along with a tuning or optimization assembly 60. In general, the tuning assembly 60 ascertains or determines control parameters or variables of the servo controller 14 based on user constraints 101. For example, the control parameters can include one or more values that are used generally in a closed-loop control architecture to generate signal 19 in part based on feedback 15 A during actual operation of the plant 10. However, the tuning assembly 60 can ascertain the control parameters offline without repeated operation of the plant 10, thereby, saving time and minimizing wear or damage to the plant 10.
  • In general, the tuning [0033] assembly 60 uses a model 64 of the system controller 14 and a model 70 of the plant 10, in one embodiment, both in the frequency domain, in order to estimate a performance index as a function of controller parameter settings and user constraints. Optimization of the performance index using an optimization algorithm such as Nelder-Mead yields the best or optimum set of controller parameters for the servo controller 14. As used herein “best” or “optimum” refers to those values of the controller parameters that have been ascertained by the tuning assembly 60 given user constraints. As appreciated by those skilled in the art, further manual adjustments to the ascertained controller parameters may be possibly made by the user to yield some further improvements in system performance. Nevertheless, the ascertained controller parameters provided by the tuning assembly 60 would achieve desired system performance based on the user constraints, and thus, can be considered “best” or optimum”.
  • The tuning [0034] assembly 60 includes an optimization module 62 and a performance calculator 63. The performance calculator 63 uses the model 64 of the system controller 14 and the model 70 of the plant 10. Before further describing operation of the performance calculator 63 and the optimization module 64, it may be helpful to describe construction of the servo controller model 64 and plant model 70.
  • As discussed above, the [0035] servo controller 14 can be embodied in a digital computer such as computer 30 illustrated in FIG. 2. In many cases, the control architecture is embodied as modules or instructions implementing control functions such as summation, subtraction, integration, differentiation and filtering, just to name a few. Signals 19 are computed using the modules or instructions and converted, if necessary, to analog values through digital-to-analog converters as is well known. Likewise, the reference signal 17 and the feedback signals 15A, if provided in analog form, can be converted to digital values through suitable analog-to-digital converters. Digital values of the reference signal 17 and the feedback signals 15 A are then used by the modules defining the control architecture.
  • The control architecture of [0036] servo controller 14 will thus be well defined. Constructing a mathematical representation or model of the servo controller 14 architecture in the frequency domain is well known. A simple exemplary controller structure for illustrative purposes can be
  • Gc=Km*(Kp+Kd*(z−1)/Tz+Kdd*((z−1)/Tz){circumflex over ( )}2),
  • where Km is a scalar and Kp, Kd, Kdd are vectors of length equal to the number of channels being controlled. This form is for the case of an integrating plant such as a servo-hydraulic system. Generally, construction includes substituting discrete values for z according to the plot of FIG. 5 for each of the controller components. In one embodiment, linear spacing along the quarter circles [0037] 81 and log spacing along the unit circle 83 realizes an accurate model with a minimum number of data points. It should also be noted that this form of the model can be referred to as a transfer function (TF). A vector of complex values becomes the frequency domain representation of the corresponding controller element. This model can be arranged in a matrix where such vectors commonly, but not exclusively, are disposed along a matrix diagonal.
  • It should be noted that if [0038] servo controller 14 is constructed of analog components, the servo controller model 64 is formed by substituting discrete values for s according to the plot of FIG. 6 for each of the controller components. A simple illustrative controller structure for this purpose is Gc=Km*(KP+Kd*s+Kdd*s{circumflex over ( )}2), where Km, Kp, Kd and Kdd are as stated above. This form is for the case of an integrating plant such as a servo-hydraulic system. To minimize the number of data points, linear spacing along the quarter circles 85 and log spacing along the imaginary axis 87 in the s-plane can be used. The models based on the s or z contours can be much more efficient for subsequent calculations (10 times fewer data values), and equally significant is that stability analysis from this data is traceable to first principals of complex analysis. Simplification of this approach to the direct use of the FRF representation can be beneficial in some situations. Actual controller structures that need to be accommodated are typically more complex, but the method of substitution of z and s contours is the same.
  • In a first embodiment, the [0039] plant model 70 can be obtained by model creation module 74 from empirical data using method 80 of FIG. 7. Beginning at step 82, servo controller 14 of FIG. 1 is operated to generate manipulated variable signal 19 as a function of a drive 17 having wide spectrum content. For instance, drive 17 can include transients, and/or repeating periodic signals having random phase. A suitable method for formation of drive 17 is described in U.S. Pat. No. 6,385,564, incorporated herein by reference. The manipulated variable signal 19 and the corresponding feedback signal 15 during operation of the plant 10 is recorded, which is schematically illustrated at 84 in FIG. 4.
  • At [0040] step 86, the FRF (frequency response function) is computed from the recorded manipulated variable signal 19 and recorded feedback signal 15A. The inverse transform of the FRF is computed to obtain the impulse response function at step 88. It should be noted that a fundamental attribute of the data collected in storage 84 is to provide one estimate of the plant 10 dynamics with sufficient fidelity and certainty. Further, storage 84 represents collection of such data, whether actually stored or used in real time. Thus, depending on the capabilities of the computer executing the tuning assembly 60 and/or controller 14, it is possible to ascertain control parameters “online” (when the plant 10 is operating) or “offline” without further operation of the plant 10.
  • To reduce the noise, which may be present in the signals, the inverse response functions can be tapered at [0041] step 89 to increase model fidelity. FIG. 8 illustrates a representative impulse response function at 90. A tapering function 92 can be applied to the impulse function 90 to remove content known to be noise. Characteristics of the tapering function 92 (i.e. duration and shape of the taper) can be adjusted depending on noise, damping, etc. characteristics of the system, but represent Fourier truncated physically realizable elements.
  • At [0042] step 96, perform forward transformation of the impulse response functions (possibly, tapered) at discrete values along the complex z domain contour (FIG. 5), or the complex s domain contour (FIG. 6), thereby producing a frequency domain model 70 consistent with the servo controller model 64, described above. This process minimizes the number of data points and produces a model that is consistent with the servo controller model 64.
  • At [0043] step 98, elements of the frequency domain model 70 can be normalized in order to ease interpretation of values used during optimization discussed below. For instance, an element representing displacement feedback can be normalized to correspond to 1/s at low frequencies. At step 100, the plant model 70 and the normalization factors are. stored.
  • The tuning [0044] assembly 60 can be implemented on a suitable computer or computers, for instance, the computer 30 discussed above. Generally, as discussed above, the tuning assembly 60 uses a model of the system controller 14 and a model of the plant 10, both in the frequency domain, in order to estimate a performance index as a function of controller parameter settings and user constraints. Optimization of the performance index using an optimization algorithm such as Nelder-Mead yields the best set of controller parameters for the servo controller 14.
  • Before discussing in detail operations performed by the tuning [0045] assembly 60, a brief review of multi-variable control theory may be helpful.
  • With all the terms illustrated in FIG. 3 determined, the closed-loop (cl) frequency response can be calculated as: [0046]
  • Gcl=[I+Km*H*Gp*Gc] −1 *Km*Gp*Gc   (1)
  • From this relationship it can be observed that, if the elements were stable (have no “rhp” (right-half-plane) poles), then any rhp poles in the closed-loop solution would result from the inverse operation. If the inverse is represented as the adjoint/determinant, it can be seen that any rhp poles would result from rhp roots of the equation: [0047]
  • determinant([I+Km*H*Gp*Gc])=0   (2)
  • Thus stability can be determined by analysis of the determinant, represented in this case as a scalar valued function of complex frequency. [0048]
  • A Nyquist plot of the determinant can thus be used to ascertain stability (note encirclement of the origin in this case). However, this does not produce a measure of stability margin. A method for establishing a stability margin is to decompose the term H*Gp*Gc (the open-loop TF with Km factored out) using a Schur decomposition to obtain: [0049]
  • determinant(U*[I+Km*GHt]*U′)=0   (3)
  • where GHt is a triangular matrix of functions of frequency, U is a unitary matrix of functions, and “′” denotes conjugate transposition. [0050]
  • Since the determinant of the product is the product of the determinants and the determinant of a diagonal matrix is the product of the diagonals, this results in [0051]
  • determinant([I+Km*GHd])=0   (4)
  • where GHd contains the diagonal elements of GHt (eigen-functions). [0052]
  • With this written as [0053]
  • (1+Km*GHd(1))*(1+Km*GHd(2))* . . . =0   (5)
  • it can be seen that the term of the product closest to zero (Km*GHd(i) closest to minus one) will determine the stability margin. It should also be observed that a plot of all Km*GHd(i) terms grows linearly with Km. [0054]
  • As indicated above, the tuning [0055] assembly 60 optimizes the controller parameters based in part on user constraints, which in the block diagram of FIG. 4 is illustrated at 101. Generally, user constraints can include constraints on independent variables, i.e., the parameters of the controller and constraints on system performance (e.g. stability margin and peaking of the closed-loop frequency response).
  • FIG. 9 illustrates a convenient representation for analyzing stability margin. In this example, the [0056] circle 102 around “−1,0” represents the stability margin boundary user constraint. The shape is user adjustable. Although a constant radius is shown, other shapes of the stability margin boundary constraint can be used. In this example, a single channel Eigen-value trajectory is plotted.
  • In yet a further embodiment, Nyquist stability can be plotted in three dimensions with frequency as the third dimension and the stability margin is illustrated as a tube or [0057] column 104, herein having a constant radius centered on “−1,0”. FIG. 10 is an example of such an illustration.
  • It should be noted that FIG. 10 illustrates a multi-channel example. With this form of illustration, a user can easily determine the important frequencies affecting stability as viewed by plotting the Eigen-value trajectories. [0058]
  • In yet a further embodiment, the stability margin boundary constraint can be allowed to vary as a function of frequency. FIG. 11 illustrates a stability margin boundary (herein a [0059] tube 106 having portions of different radii centered at varying positions on the real-imaginary plane) allowing increased stability margin in regions of greater uncertainty such as at higher frequencies. As appreciated by those skilled in the art, viewing of the three-dimensional renderings of FIGS. 10 and 11 can be changed as desired to allow the trajectories to be readily seen with respect to the stability margin. Although the shape of the stability margin boundary 106 can be adjusted, predetermined shapes can be stored and selected by the user.
  • As discussed above, besides accounting for stability, other possible system performance user constraints can include frequency domain peaking. FIG. 12 graphically illustrates an absolute [0060] upper limit 108 throughout the entire frequency spectrum. FIG. 12 also illustrates a “tuning bandwidth” 110 that is user selectable. The tuning bandwidth defines the frequency region of primary interest for closed-loop performance.
  • Referring back to the tuning [0061] assembly 60, a tuning procedure is generally illustrated in FIG. 13 at 130. At step 132, the controller model 64 and the plant model 70 are obtained as discussed above.
  • At [0062] step 133, the user can select all or a subset of channels (manipulated variables and feedback responses) to tune or optimize. In many system environments, cross-coupling between channels can be more prevalent for some channels than for others. For instance, in a road simulator discussed generally in the background section, loads and motions can be applied to each of the vehicle spindles from separate loading assemblies coupled to each of the vehicle spindles. Each loading assembly comprises separate actuators and linkages that are operated to apply forces and motions with respect to an orthogonal set of axes, generally with one axis corresponding to the spindle axis. Each loading assembly includes a set of manipulated variable signals for controlling the servo valves of the actuators and a set of corresponding feedback signals. If each vehicle spindle includes a loading assembly, then four sets of channels (manipulated variable signals and feedback signals) are provided to the servo controller. However, cross-coupling may be minimal between some sets of loading assemblies, or even some channels in any one loading assembly. For instance, cross-coupling is generally low between the front and rear vehicle spindle loading assemblies.
  • [0063] Step 133 allows the user to select which channels (a subset or all) to tune or optimize. In this manner, if desired, the user can individually select single input and single output channels to optimize, or select a subset of channels to optimize. Generally, it is faster to optimize a plurality of subsets of channels or portions of the plant, than to optimize the complete plant at once. By processing the plant in subsets, some cross-coupling is ignored, but as stated above, this may be acceptable given the configuration of the plant.
  • It should also be noted that by allowing the user to select individual or subset of channels to optimize, plant diagnostics can be provided. For instance, it may be apparent from the plant model or from the optimization procedure that channels are totally dysfunctional or poorly performing. In this manner, the optimization or tuning procedure can be used to diagnose problems in the plant. This is particularly true when using single channel tuning repeated over all the channels since this determines the adequacy of performance of each channel operating without interaction with other channels. [0064]
  • At [0065] step 134, the user provides user constraint values. As described above, user constraint values include stability margins, if desired, as a function of frequency, a frequency domain peaking limit and the tuning bandwidth.
  • In one embodiment, the gains in the [0066] controller model 64 are normalized by frequency by a nominal frequency of operation. At step 134, the user constraint values can also include an expected system bandwidth 110 (FIG. 12).
  • At [0067] step 136, initial values of the tuning parameters are provided. If the tuning parameters were normalized, each of the initial values can be set to their nominal value of one.
  • The optimization routine, such as Nelder-Mead, is then invoked at [0068] step 138. The optimization routine repeatedly provides a set of controller parameters 140 to the performance calculator 63, and a performance index 142 is returned from the performance calculator 63 until an optimum value is achieved (i.e. within a selected tolerance range). The optimization routines are readily available and well known, such as in the routines available from Mathworks, Inc. of Natick, Mass., U.S.A.
  • Calculation of the performance index by the [0069] performance calculator 63 is illustrated in FIG. 9 as method 150. At step 152, a set of tuning parameters (independent variables) is received by the performance calculator 63 from the optimizer module 62. Generally, the tuning parameters comprise the derivative gains, the second derivative gains (if present), various feedback and/or forward loop filter settings. Optionally, a function of the proportional gains can be included. In many instances, the proportional gains can be incorporated into the performance index, i.e., maximizing the proportional gain can be a desired parameter for a single channel system.
  • At [0070] step 154, the open-loop TF model (H*Gp*Gc) is computed as a function of the tuning parameters accounting for any constraints on the independent variables including filter settings. It should be noted that such constraints can also be conveniently handled by a optimization routine, which handles constraints directly.
  • At [0071] step 156, the open-loop TF model is decomposed frequency-by-frequency using the Schur decomposition that produces the Eigen-functions. If desired, Eigen value decomposition can be used.
  • The master gain Km is adjusted (typically an iterative algorithm) at [0072] step 158 such that the Eigen-function trajectories contact the stability margin region (nominally a region, for example, a circle, ellipse, etc. that includes (−1,0) of the complex plane).
  • At [0073] step 160, the closed-loop plant response is computed as [I+Km*H*Gp*Gc]−1 *Km*Gp*Gc. If peaking of the closed-loop response exceeds a specified maximum value, lower Km until the specified maximum value is achieved (typically an iterative algorithm). A property of Km is that the stability margin is increased and so the constraint on stability is still met.
  • The value obtained for Km in [0074] step 160 plus optionally an additive weighted performance value is returned at step 162 to the optimizer module 62.
  • It should be noted at this time in a further embodiment illustrated in FIG. 15 that the plant can comprise multiple configurations indicated at [0075] 10, 10A and 10B. For example, the different configurations 10, 10A and 10B can be for different test specimens or no specimen, plant operating levels, or other variations of plant operations. For each of the plant configurations 10, 10A and 10B, a corresponding plant model 70, 70A and 70B would be obtained. The tuning assembly 60 can ascertain the best set of controller parameters that would be useable over all of the various plant configurations. In particular, one set of tuning parameters would be provided at step 152, but each of the steps 154, 156, 158 and 160 would be performed for each of the plant configurations. This entails at step 158 to adjust the master gain Km until each plant configuration satisfies the stability margin for all system model configurations. Once the maximum master gain Km has been determined that produces stability in all models, at step 160 the master gain Km is reduced further, if necessary, so that the peaking limit is met across all channels for all plant configurations. In further embodiments, the peaking limit can be a function of each channel and/or a function of frequency.
  • Referring back to FIG. 13, at [0076] step 170, visual renderings can be provided to the user to show stability and system performance using the optimized parameters. The renderings, comprising another aspect of the present invention, can include any of the three-dimensional Nyquist plots described above. By plotting the Eigen-value trajectories, the user can readily confirm that system stability has been been met because the trajectories will not visually intersect with the selected stability margins. Likewise frequency peaking can be visually rendered such as illustrated in FIG. 12.
  • If the plant was optimized based on processing a plurality of subsets of channels, it may be beneficial to check the stability and performance of the system model using the controller parameters ascertained based on optimizing the parameters in subset calculations. At [0077] step 172, visual renderings similar to step 170 for system stability and system performance can be provided to the user for the complete system. If any of the Eigen-value trajectories visually intersect with the selected stability margins, system stability may not be acceptable. Likewise frequency peaking can be checked by visually rendering system performance such as illustrated in FIG. 12. It should be noted visual rendering is not required, but is a convenient form for interpreting system stability and/or system performance. Of course, other forms of indications (e.g. numerical tables, audible alarms, etc.) can be provided to the user to convey violations of system model stability or system model performance.
  • At [0078] step 174, using the stored normalization values of the plant model and the controller model, the control parameters for the controller 14 are calculated as a function of the optimum set of controller parameters determined by the optimization routine.
  • The “performance index” is a measure of “goodness” of the servo-system. This generally is a measure of control bandwidth. There is a convenient relationship between controller parameters and control bandwidth when using a normalized plant (scaled to 1/s), that is: [0079]
  • F BW =Km*Kp*Hz   (6)
  • Note that with the normalized plant, Kp can be treated as having units of rad/sec, the loop bandwidth. [0080]
  • In a multi-channel system Kp is a vector of length equal to the number of channels (Km is always a scalar) . If Kp is a vector of equal constants (if normalized, i.e. the proportional gains being of a selected relationship among each other), then all channels will have equal loop bandwidth with the loop bandwidth value determined by Km*Kp. Thus maximizing Km over the tuning parameters (derivative gains, filter settings, etc.) makes sense. All channels operating with the same loop bandwidth is highly desirable in systems such as multi-axis tracking systems. In this case, system performance is limited by the weakest channel by necessity. [0081]
  • The weakest channel in a system is determined by individual loop dynamics as well as cross-coupling to other channels. When cross-coupling is present, channels are sharing stability in some sense. Thus increasing loop-bandwidth of one channel can require that loop-bandwidth be lowered on one or more other channels. [0082]
  • Distributing loop-bandwidth amongst channels can be accomplished be setting the elements of the Kp vector to different values before optimizing Km. The channel least able to make its relative loop-bandwidth will produce the value of Km at the optimum. [0083]
  • In cases where it is desirable to have each channel operating at as high a loop-bandwidth as possible together with all other channels, the optimization process can take on a more complex form. A useful example is: [0084] Performance Index = i ( K m K p , i ω ref ) x ( 7 )
    Figure US20030139825A1-20030724-M00001
  • with ω[0085] ref set to the nominal loop-bandwidth and x having a value between 0 and 1 with a nominal value of 0.5. This performance index puts the greatest emphasis on the weakest channels, while encouraging channels that are not compromised by cross-coupling to have as high a loop-bandwidth as possible. With this performance index, optimization must be performed over the Kp's as well as the tuning parameters previously mentioned. Values for x greater than 0.5 puts even greater emphasis on the weaker channels, while values less than 0.5 puts less emphasis on the weaker channels.
  • Another useful variation of the performance index is to include a factor that puts more focus on a particular frequency range. An example is the minimum (graphically illustrated at [0086] 109 in FIG. 12) of the closed-loop frequency response over a user specified frequency range. This can be applied to each individual channel (diagonal elements) or taken as the minimum over all diagonal elements.
  • When applied to individual channels (diagonal elements) this effect can logically be included in the loop balancing performance index of Eq. 7. [0087]
  • Another variation of the performance index is to apply Eq. 7 to preselected groups of channels rather than to each channel individually. Each term in the sum would be a measure of the performance of the group to be balanced against the group. This can be especially useful in the case of symmetry, side-to-side symmetry of channels of the plant, for example. In this case, the corresponding left-side and the right-side channels can be treated as a pair and form a group in the overall process. [0088]
  • A wide variety of multivariable servo system models can be represented using the diagram illustrated in FIG. 3. However in other multivariable systems an alternative embodiment of FIG. 16 can be beneficial. In this configuration, the [0089] plant 10 is represented as a combination of a plurality (e.g. four) controller blocks. In particular, there is a separate block or model for each feedback variable type. In the illustrated example, each of the models G1, G2, G3 and G4 represents a model of each of the controlled variables as a function the manipulated variable all in a multivariable system sense. For instance, the outputs of G1-G4 in a vibration system, can correspond to the controlled variable, displacement feedback, velocity feedback and acceleration feedback, respectively. Kp, Kd and Kdd are adjustable vector of gains for G2, G3 and G4, respectively. Km is the master gain similar to that discussed above.
  • This configuration of FIG. 16 can be tuned using the [0090] tuning assembly 60 in an optimal manner similar to FIG. 3, where the closed loop frequency response can be represented as:
  • Gcl=G1*[I+Km*(KpG2+KdG3+KddG4)]−1 KmKp   (9)
  • and where the characteristic equation is: [0091]
  • Determinant([I+Km*(KpG2+KdG3+KddG4)]=0   (10)
  • It should also be noted that the present invention can also be applied to servo control configurations involving a loop within a loop. One example is the well-known Cascade configuration having an [0092] inner loop 180 and an outer loop 182, the outer loop 182 including the inner loop 180 as illustrated in FIG. 16. The tuning procedure can be first applied to the inner loop 180. With the control parameters of inner loop 180 ascertained, the inner loop 180 becomes part of the plant for the outer loop 182, at which point the tuning procedure can then be applied to the outer loop 182 to ascertain the control parameters thereof. In yet a further embodiment, the tuning or optimization procedure of the outer loop 182 can include as a portion of its performance index calculation, the tuning procedure of the inner loop 180.
  • While the description thus far covers closed-loop tuning/optimization it is straightforward to extend the system/model to include open-loop series compensation since stability is not affected by open-loop series compensation. [0093]
  • Although the present invention has been described with reference to preferred embodiments, workers skilled in the art will recognize that changes may be made in form and detail without departing from the spirit and scope of the invention. For instance, while frequency domain calculations are often optimal, in some situations, it may be beneficial to use other model forms such as finite element models or kinematic models. [0094]

Claims (10)

What is claimed is:
1. A method of ascertaining control parameters for a control system having a controller and a plant, the method comprises:
establishing a model of the plant and a model of the controller; and
calculating a performance index for a closed-loop system as a function of controller parameters, accounting for selected stability margins.
2. The method of claim 1 and further including: providing an optimization routine, the optimization routine repeatedly performing the step of calculating the performance index.
3. The method of claim 2 wherein the optimization routine repeatedly performs the step of calculating the performance index until an optimum value is achieved.
4. The method of claim 3 wherein the plant model and the controller model each comprise multiple channels, and wherein the controller parameters comprise a vector Kp having elements pertaining to proportional gains for a each channel of the controller model, the proportional gains being of a selected relationship among each other, and wherein the performance index is a function of a master proportional gain Km that is applied to each channel of the controller model, and wherein the performance index is repeatedly calculated to maximize Km.
5. The method of claim 3 wherein the plant model and the controller model each comprise multiple channels, and wherein the controller parameters comprise a vector Kp having elements pertaining to proportional gains for a each channel of the controller model, and wherein the performance index is a function of a master proportional gain Km that is applied to each channel of the controller model, and wherein the performance index is calculated as a function of:
i ( K m K p , i ω ref ) x
Figure US20030139825A1-20030724-M00002
where with ωref is set to a nominal loop-bandwidth, x comprises a value between 0 and 1 and i is a running index over all channels of the controller model.
6. The method of claim 1 wherein the stability margin is a function of frequency and represented in three dimensions.
7. The method of claim 1 wherein calculating the performance index includes accounting for a limit on frequency domain peaking.
8. The method of claim 7 wherein the limit comprises a constant.
9. The method of claim 1 wherein calculating the performance index includes accounting for a minimum value of frequency response over a specified range of frequencies.
10. The method of claim 1 wherein establishing a controller model and plant model includes establishing the controller model and the plant model in the frequency domain.
US10/323,707 2001-12-18 2002-12-18 Method of ascertaining control parameters for a control system Expired - Lifetime US6917840B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US10/323,707 US6917840B2 (en) 2001-12-18 2002-12-18 Method of ascertaining control parameters for a control system

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US34188701P 2001-12-18 2001-12-18
US10/323,707 US6917840B2 (en) 2001-12-18 2002-12-18 Method of ascertaining control parameters for a control system

Publications (2)

Publication Number Publication Date
US20030139825A1 true US20030139825A1 (en) 2003-07-24
US6917840B2 US6917840B2 (en) 2005-07-12

Family

ID=23339433

Family Applications (1)

Application Number Title Priority Date Filing Date
US10/323,707 Expired - Lifetime US6917840B2 (en) 2001-12-18 2002-12-18 Method of ascertaining control parameters for a control system

Country Status (9)

Country Link
US (1) US6917840B2 (en)
EP (1) EP1466219B1 (en)
JP (2) JP2005513601A (en)
KR (1) KR100989395B1 (en)
CN (1) CN100380254C (en)
AT (1) ATE413633T1 (en)
AU (1) AU2002364184A1 (en)
DE (1) DE60229782D1 (en)
WO (1) WO2003052529A1 (en)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030199997A1 (en) * 2002-04-18 2003-10-23 Zhiqiang Gao Scaling and parameterizing a controller
US20070078529A1 (en) * 2005-09-30 2007-04-05 Fisher-Rosemount Systems, Inc. On-line adaptive model predictive control in a process control system
US20070118237A1 (en) * 2005-11-18 2007-05-24 Chun-Chieh Wang Autocontrol simulating system and method
US7317953B2 (en) * 2003-12-03 2008-01-08 Fisher-Rosemount Systems, Inc. Adaptive multivariable process controller using model switching and attribute interpolation
US7346402B1 (en) * 2004-10-22 2008-03-18 Expertcontrol, Llc Technique for an integrated and automated model generation and controller tuning process
US20090005886A1 (en) * 2002-04-18 2009-01-01 Cleveland State University Extended Active Disturbance Rejection Controller
US7853431B2 (en) 2006-09-29 2010-12-14 Fisher-Rosemount Systems, Inc. On-line monitoring and diagnostics of a process using multivariate statistical analysis
US20110288846A1 (en) * 2010-05-21 2011-11-24 Honeywell International Inc. Technique and tool for efficient testing of controllers in development (h-act project)
US8710777B2 (en) 2012-04-20 2014-04-29 Linestream Technologies Method for automatically estimating inertia in a mechanical system
US20150015631A1 (en) * 2013-07-11 2015-01-15 Dan C. Lyman Adaptive control of continuous inkjet parameters
US9041337B2 (en) 2012-05-18 2015-05-26 Linestream Technologies Motion profile generator
US9158291B1 (en) * 2010-06-04 2015-10-13 The Mathworks, Inc. Automatic H-infinity optimization for multivariable control systems
WO2016118701A1 (en) * 2015-01-21 2016-07-28 Linestream Technologies Cascaded active disturbance rejection controllers
US10061275B2 (en) 2014-07-29 2018-08-28 Linestream Technologies Optimized parameterization of active disturbance rejection control
US10126202B2 (en) 2015-09-11 2018-11-13 Linestream Technologies Method for automatically estimating inertia, coulomb friction, and viscous friction in a mechanical system
CN112415897A (en) * 2020-11-30 2021-02-26 北京罗克维尔斯科技有限公司 Control system optimization method, device, system, vehicle, medium and equipment

Families Citing this family (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030035376A1 (en) * 2001-08-20 2003-02-20 Xiaofen Chen Derivation of composite step-function response
JP4223894B2 (en) * 2003-08-21 2009-02-12 株式会社山武 PID parameter adjustment device
US7363094B2 (en) * 2006-01-09 2008-04-22 General Electric Company Multivariable controller design method for multiple input/outputs systems with multiple input/output constraints
US7987005B2 (en) * 2006-12-19 2011-07-26 Chevron U.S.A. Inc. System, method and program for dynamic control and optimization of a process having manipulated and controlled variables
KR100925874B1 (en) * 2007-12-10 2009-11-09 에스케이텔레시스 주식회사 Computer simulator of RF interference cancellation repeaters and implementing method therefor
JP2010039966A (en) 2008-08-08 2010-02-18 Hitachi Ltd Data management method
US7990818B2 (en) * 2008-12-23 2011-08-02 General Electric Company Cascaded control of a pick-up head for multi-layer optical data storage
US8631046B2 (en) 2009-01-07 2014-01-14 Oracle International Corporation Generic ontology based semantic business policy engine
US9672478B2 (en) 2009-02-26 2017-06-06 Oracle International Corporation Techniques for semantic business policy composition
US8949236B2 (en) 2010-02-26 2015-02-03 Oracle International Corporation Techniques for analyzing data from multiple sources
CN102279563A (en) * 2010-06-13 2011-12-14 天津理工大学 Method for establishing parameterized kinematics model library for nonstandard digital equipment operator
US9400958B2 (en) 2010-06-30 2016-07-26 Oracle International Corporation Techniques for display of information related to policies
US8954309B2 (en) * 2011-05-31 2015-02-10 Oracle International Corporation Techniques for application tuning
US8898096B2 (en) 2011-05-31 2014-11-25 Oracle International Corporation Application configuration generation
US20130150984A1 (en) * 2011-12-13 2013-06-13 Richard M. Nigro Test system with configurable closed loop
US9574505B2 (en) * 2012-08-29 2017-02-21 Toyota Jidosha Kabushiki Kaisha Plant control device
FR3039130B1 (en) * 2015-07-21 2019-05-31 Airbus Operations METHOD AND DEVICE FOR AUTOMATICALLY MANAGING A ACTUATOR CONTROLLED BY A SERVO-VALVE
TWI681274B (en) * 2018-11-14 2020-01-01 財團法人工業技術研究院 Method for adjusting frequency response parameter of machine tool and system using the same
US11204848B1 (en) 2020-12-15 2021-12-21 International Business Machines Corporation System testing infrastructure with hidden variable, hidden attribute, and hidden value detection
US11188453B1 (en) 2020-12-15 2021-11-30 International Business Machines Corporation Verification of software test quality using hidden variables
US11379352B1 (en) 2020-12-15 2022-07-05 International Business Machines Corporation System testing infrastructure with hidden variable, hidden attribute, and hidden value detection
US11132273B1 (en) 2020-12-15 2021-09-28 International Business Machines Corporation System testing infrastructure with hidden variable, hidden attribute, and hidden value detection
US11113167B1 (en) 2020-12-15 2021-09-07 International Business Machines Corporation System testing infrastructure with hidden variable, hidden attribute, and hidden value detection
CN113325699B (en) * 2021-05-25 2022-08-12 上海机电工程研究所 Parameter adjusting method and system suitable for composite stability control system
CN115079577B (en) * 2022-07-22 2022-11-11 浙江中控技术股份有限公司 Closed loop step test method and test device based on real-time control performance evaluation

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4195337A (en) * 1977-05-26 1980-03-25 Agence Nationale De Valorisation De La Recherche (Anvar) Control method utilizing a model control scheme
US4214300A (en) * 1977-05-17 1980-07-22 Kenneth Robert Jones Three term (PID) controllers
US4234142A (en) * 1978-06-08 1980-11-18 The United States Of America As Represented By The Secretary Of The Navy High angle-of-attack missile control system for aerodynamically controlled missiles
US4882526A (en) * 1986-08-12 1989-11-21 Kabushiki Kaisha Toshiba Adaptive process control system
US5329443A (en) * 1992-06-16 1994-07-12 Praxair Technology, Inc. Two-phase method for real time process control
US5400247A (en) * 1992-06-22 1995-03-21 Measurex Corporation, Inc. Adaptive cross-directional decoupling control systems
US5623402A (en) * 1994-02-10 1997-04-22 Schenck Pegasus Corporation Multi-channel inverse control using adaptive finite impulse response filters
US5689427A (en) * 1995-09-01 1997-11-18 Industrial Technology Research Institute Method and apparatus for optimal feedrate control on wire-cutting electric discharging machine (WEDM)
US6128541A (en) * 1997-10-15 2000-10-03 Fisher Controls International, Inc. Optimal auto-tuner for use in a process control network
US6185470B1 (en) * 1997-11-07 2001-02-06 Mcdonnell Douglas Corporation Neural network predictive control method and system
US6195596B1 (en) * 1997-10-24 2001-02-27 Daewoo Electronics Co., Ltd. Motion effector controller for use in a virtual reality system
USRE37374E1 (en) * 1995-10-26 2001-09-18 Cybernet Haptic Systems Corporation Gyro-stabilized platforms for force-feedback applications
US20020111758A1 (en) * 2000-10-18 2002-08-15 Qing-Guo Wang Robust process identification and auto-tuning control
US6640638B1 (en) * 1999-09-03 2003-11-04 Mts Systems Corporation Loading assembly for a vehicle spindle test fixture
US6701195B2 (en) * 1998-12-17 2004-03-02 Siemens Aktiengesellschaft Sensor prediction system utilizing case based reasoning

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4537077A (en) 1984-02-08 1985-08-27 Mts Systems Corporation Load dynamics compensation circuit for servohydraulic control systems
JPH02308305A (en) * 1989-05-24 1990-12-21 Toshiba Corp System modeling device
JP2954660B2 (en) * 1990-05-30 1999-09-27 株式会社東芝 Model predictive controller
JP3502208B2 (en) * 1995-10-25 2004-03-02 太平洋セメント株式会社 Process identification method and apparatus
EP0944865A1 (en) * 1996-12-11 1999-09-29 Siemens Aktiengesellschaft Automatic setting method for a proportional-integral-differential controller (pid) regulating a linear servo system
US6253113B1 (en) 1998-08-20 2001-06-26 Honeywell International Inc Controllers that determine optimal tuning parameters for use in process control systems and methods of operating the same
JP3469136B2 (en) * 1999-08-20 2003-11-25 三菱重工業株式会社 Waveform control device

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4214300A (en) * 1977-05-17 1980-07-22 Kenneth Robert Jones Three term (PID) controllers
US4195337A (en) * 1977-05-26 1980-03-25 Agence Nationale De Valorisation De La Recherche (Anvar) Control method utilizing a model control scheme
US4234142A (en) * 1978-06-08 1980-11-18 The United States Of America As Represented By The Secretary Of The Navy High angle-of-attack missile control system for aerodynamically controlled missiles
US4882526A (en) * 1986-08-12 1989-11-21 Kabushiki Kaisha Toshiba Adaptive process control system
US5329443A (en) * 1992-06-16 1994-07-12 Praxair Technology, Inc. Two-phase method for real time process control
US5400247A (en) * 1992-06-22 1995-03-21 Measurex Corporation, Inc. Adaptive cross-directional decoupling control systems
US5623402A (en) * 1994-02-10 1997-04-22 Schenck Pegasus Corporation Multi-channel inverse control using adaptive finite impulse response filters
US5689427A (en) * 1995-09-01 1997-11-18 Industrial Technology Research Institute Method and apparatus for optimal feedrate control on wire-cutting electric discharging machine (WEDM)
USRE37374E1 (en) * 1995-10-26 2001-09-18 Cybernet Haptic Systems Corporation Gyro-stabilized platforms for force-feedback applications
US6128541A (en) * 1997-10-15 2000-10-03 Fisher Controls International, Inc. Optimal auto-tuner for use in a process control network
US6195596B1 (en) * 1997-10-24 2001-02-27 Daewoo Electronics Co., Ltd. Motion effector controller for use in a virtual reality system
US6185470B1 (en) * 1997-11-07 2001-02-06 Mcdonnell Douglas Corporation Neural network predictive control method and system
US6701195B2 (en) * 1998-12-17 2004-03-02 Siemens Aktiengesellschaft Sensor prediction system utilizing case based reasoning
US6640638B1 (en) * 1999-09-03 2003-11-04 Mts Systems Corporation Loading assembly for a vehicle spindle test fixture
US20020111758A1 (en) * 2000-10-18 2002-08-15 Qing-Guo Wang Robust process identification and auto-tuning control

Cited By (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140195013A1 (en) * 2002-04-18 2014-07-10 Cleveland State University Extended active disturbance rejection controller
US8644963B2 (en) 2002-04-18 2014-02-04 Cleveland State University Extended active disturbance rejection controller
US20030199997A1 (en) * 2002-04-18 2003-10-23 Zhiqiang Gao Scaling and parameterizing a controller
US8180464B2 (en) 2002-04-18 2012-05-15 Cleveland State University Extended active disturbance rejection controller
US8571691B2 (en) 2002-04-18 2013-10-29 Cleveland State University Scaling and parameterizing a controller
US20090005886A1 (en) * 2002-04-18 2009-01-01 Cleveland State University Extended Active Disturbance Rejection Controller
US8041436B2 (en) * 2002-04-18 2011-10-18 Cleveland State University Scaling and parameterizing a controller
US9268315B2 (en) * 2002-04-18 2016-02-23 Cleveland State University Extended active disturbance rejection controller
US8406905B2 (en) 2002-04-18 2013-03-26 Cleveland State University Scaling and parameterizing a controller
US9052705B2 (en) 2002-04-18 2015-06-09 Cleveland State University Scaling and parameterizing a controller
US9720387B2 (en) 2002-04-18 2017-08-01 Cleveland State University Scaling and parameterizing a controller
US7317953B2 (en) * 2003-12-03 2008-01-08 Fisher-Rosemount Systems, Inc. Adaptive multivariable process controller using model switching and attribute interpolation
US7346402B1 (en) * 2004-10-22 2008-03-18 Expertcontrol, Llc Technique for an integrated and automated model generation and controller tuning process
US20070078529A1 (en) * 2005-09-30 2007-04-05 Fisher-Rosemount Systems, Inc. On-line adaptive model predictive control in a process control system
US20090143872A1 (en) * 2005-09-30 2009-06-04 Fisher-Rosemount Systems, Inc. On-Line Adaptive Model Predictive Control in a Process Control System
US7856281B2 (en) 2005-09-30 2010-12-21 Fisher-Rosemount Systems, Inc. On-line adaptive model predictive control in a process control system
US20070118237A1 (en) * 2005-11-18 2007-05-24 Chun-Chieh Wang Autocontrol simulating system and method
US7853339B2 (en) * 2006-09-29 2010-12-14 Fisher-Rosemount Systems, Inc. Statistical signatures used with multivariate analysis for steady-state detection in a process
US7853431B2 (en) 2006-09-29 2010-12-14 Fisher-Rosemount Systems, Inc. On-line monitoring and diagnostics of a process using multivariate statistical analysis
US8489360B2 (en) 2006-09-29 2013-07-16 Fisher-Rosemount Systems, Inc. Multivariate monitoring and diagnostics of process variable data
US8014880B2 (en) 2006-09-29 2011-09-06 Fisher-Rosemount Systems, Inc. On-line multivariate analysis in a distributed process control system
US7917240B2 (en) 2006-09-29 2011-03-29 Fisher-Rosemount Systems, Inc. Univariate method for monitoring and analysis of multivariate data
US7966149B2 (en) 2006-09-29 2011-06-21 Fisher-Rosemount Systems, Inc. Multivariate detection of transient regions in a process control system
US7937164B2 (en) 2006-09-29 2011-05-03 Fisher-Rosemount Systems, Inc. Multivariate detection of abnormal conditions in a process plant
US20110288846A1 (en) * 2010-05-21 2011-11-24 Honeywell International Inc. Technique and tool for efficient testing of controllers in development (h-act project)
US9760073B2 (en) * 2010-05-21 2017-09-12 Honeywell International Inc. Technique and tool for efficient testing of controllers in development
US9158291B1 (en) * 2010-06-04 2015-10-13 The Mathworks, Inc. Automatic H-infinity optimization for multivariable control systems
US9268316B2 (en) 2012-04-20 2016-02-23 Linestream Technologies Method for automatically estimating a friction coefficient in a mechanical system
US8710777B2 (en) 2012-04-20 2014-04-29 Linestream Technologies Method for automatically estimating inertia in a mechanical system
US9041337B2 (en) 2012-05-18 2015-05-26 Linestream Technologies Motion profile generator
US20150015631A1 (en) * 2013-07-11 2015-01-15 Dan C. Lyman Adaptive control of continuous inkjet parameters
US10061275B2 (en) 2014-07-29 2018-08-28 Linestream Technologies Optimized parameterization of active disturbance rejection control
WO2016118701A1 (en) * 2015-01-21 2016-07-28 Linestream Technologies Cascaded active disturbance rejection controllers
US10126202B2 (en) 2015-09-11 2018-11-13 Linestream Technologies Method for automatically estimating inertia, coulomb friction, and viscous friction in a mechanical system
CN112415897A (en) * 2020-11-30 2021-02-26 北京罗克维尔斯科技有限公司 Control system optimization method, device, system, vehicle, medium and equipment

Also Published As

Publication number Publication date
EP1466219A1 (en) 2004-10-13
JP5048693B2 (en) 2012-10-17
JP2005513601A (en) 2005-05-12
CN100380254C (en) 2008-04-09
US6917840B2 (en) 2005-07-12
JP2009104659A (en) 2009-05-14
WO2003052529A1 (en) 2003-06-26
AU2002364184A1 (en) 2003-06-30
KR20040068305A (en) 2004-07-30
KR100989395B1 (en) 2010-10-25
DE60229782D1 (en) 2008-12-18
EP1466219B1 (en) 2008-11-05
CN1618044A (en) 2005-05-18
ATE413633T1 (en) 2008-11-15

Similar Documents

Publication Publication Date Title
US6917840B2 (en) Method of ascertaining control parameters for a control system
US8055358B2 (en) Multi-objective predictive process optimization with concurrent process simulation
KR100518292B1 (en) Process control system
KR100371728B1 (en) Feedback method for controlling non-linear processes
DE69321952T2 (en) SYSTEM AND METHOD FOR CONTROLLING A PLANT
EP2423766B1 (en) Model predictive controller solution analysis process
US6564194B1 (en) Method and apparatus for automatic synthesis controllers
CA2305334C (en) Optimal auto-tuner for use in a process control network
DE69701878T2 (en) AUTOMATIC DEVICE BASED ON NON-LINEAR APPROXIMATION METHOD
US20030199997A1 (en) Scaling and parameterizing a controller
AU2001240799A1 (en) Process control system
JP4908433B2 (en) Control parameter adjustment method and control parameter adjustment program
Narayanan et al. Stochastic optimal active control of a 2-dof quarter car model with non-linear passive suspension elements
Spiewak et al. A predictive monitoring and diagnosis system for manufacturing
EP1861752B1 (en) Control of processes
O'Reilly Multivariable control for industrial applications
Soderling et al. Servo controller compensation methods selection of the correct technique for test applications
McAvoy Model predictive statistical process control of chemical plants
Koza et al. Evolution of a controller with a free variable using genetic programming
Müller et al. Optimization of inverse model identification for multi-axial test rig control
US6847851B1 (en) Apparatus for improved general-purpose PID and non-PID controllers
Dursun et al. Tracking control solution for road simulators: model-based iterative learning control approach improved by time-domain modelling
Aguiar et al. Comparative study of tuning techniques for fractional PID controllers for FOPDT systems
Strand et al. Control of the ball and plate system: A comparative study of controllers
Koza et al. Automatic Synthesis of Controllers

Legal Events

Date Code Title Description
AS Assignment

Owner name: MTS SYSTEMS CORPORATION, MINNESOTA

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:LUND, RICHARD A.;REEL/FRAME:013932/0621

Effective date: 20030319

STCF Information on status: patent grant

Free format text: PATENTED CASE

FPAY Fee payment

Year of fee payment: 4

FPAY Fee payment

Year of fee payment: 8

AS Assignment

Owner name: JPMORGAN CHASE BANK, N.A., AS ADMINISTRATIVE AGENT, ILLINOIS

Free format text: SUPPLEMENTAL CONFIRMATORY GRANT OF SECURITY INTEREST IN UNITED STATES PATENTS;ASSIGNOR:MTS SYSTEMS CORPORATION;REEL/FRAME:039258/0587

Effective date: 20160705

Owner name: JPMORGAN CHASE BANK, N.A., AS ADMINISTRATIVE AGENT

Free format text: SUPPLEMENTAL CONFIRMATORY GRANT OF SECURITY INTEREST IN UNITED STATES PATENTS;ASSIGNOR:MTS SYSTEMS CORPORATION;REEL/FRAME:039258/0587

Effective date: 20160705

FPAY Fee payment

Year of fee payment: 12

AS Assignment

Owner name: MTS SYSTEMS CORPORATION, MINNESOTA

Free format text: RELEASE BY SECURED PARTY;ASSIGNOR:JPMORGAN CHASE BANK, N.A., AS ADMINISTRATIVE AGENT;REEL/FRAME:055911/0315

Effective date: 20210407