CN105045233A - Optimum design method for PID (Proportion Integration Differentiation) controller based on time dimension in heat-engine plant thermal system - Google Patents

Optimum design method for PID (Proportion Integration Differentiation) controller based on time dimension in heat-engine plant thermal system Download PDF

Info

Publication number
CN105045233A
CN105045233A CN201510404996.7A CN201510404996A CN105045233A CN 105045233 A CN105045233 A CN 105045233A CN 201510404996 A CN201510404996 A CN 201510404996A CN 105045233 A CN105045233 A CN 105045233A
Authority
CN
China
Prior art keywords
pid
time
pid controller
parameter
controlled device
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
CN201510404996.7A
Other languages
Chinese (zh)
Other versions
CN105045233B (en
Inventor
陈宝林
华山
王德华
吴雨浓
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.)
Guodian Science and Technology Research Institute Co Ltd
Original Assignee
Guodian Science and Technology Research Institute Co Ltd
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 Guodian Science and Technology Research Institute Co Ltd filed Critical Guodian Science and Technology Research Institute Co Ltd
Priority to CN201510404996.7A priority Critical patent/CN105045233B/en
Publication of CN105045233A publication Critical patent/CN105045233A/en
Application granted granted Critical
Publication of CN105045233B publication Critical patent/CN105045233B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/418Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM]
    • G05B19/41885Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by modeling, simulation of the manufacturing system
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Feedback Control In General (AREA)
  • General Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)

Abstract

An optimum design method for a PID (Proportion Integration Differentiation) controller based on a time dimension in a heat-engine plant thermal system. Under a condition that the mathematical model of a controlled object is not completely mastered, the transfer function of a PID controller is GPID (s)=Kf(Kp+Ki/s+Kds/(Tds+1)); a first step of setting is to set parameters ranges of Kp, Ti and Kd of the PID controller, and a golden section method is used to determine values of Kp, Ti and Kd; and a second step of setting is to observe system closed-loop output at set parameters of Kp, Ti and Kd of the PID controller in the first step, and performing parameter setting for two or three times using the golden section method to obtain a satisfied control effect. The invention can complete debugging well and fast and can achieve the satisfied control quality.

Description

Based on the Optimization Design of the PID controller of time metric in Power Plant Thermal system
Technical field
The present invention relates to a kind of in Power Plant Thermal system, optimization based on the model-free PID controller of time metric is adjusted method for designing, golden section optimization method is applied in design, with other apply golden section optimization method design PID controller parameter unlike, method for designing of the present invention does not need by the mathematical model of control object, and increases pre-factor K fdesign, this is the principal feature of method for designing of the present invention.It is applied to Power Plant Thermal Systematical control, belongs to automatic control technology field.
Background technology
In industrial process control, PID controller is most widely used controller.Usually, engineers design's PID controller, needs the accurate model knowing controlled device, but in Practical Project, accurate model is often difficult to obtain.When model-free, design and Tuning PID Controller parameter are complicated debug processs, if the engineering technical personnel at scene do not understand the characteristic of controlled device, unsuitable adjustment method may make system disperse, cause industrial accident, this debug process is a job of wasting time and energy, and uncertainly can obtain satisfied Control platform.Present patent application is invented and is devised design and the setting method of measuring the model-free PID controller of indignant thought based on system time, and it better can complete debug process faster, and can obtain satisfied Control platform.Apply golden section optimization method in program, with other apply golden section optimization method design PID controller parameter unlike, present patent application method for designing does not need by the mathematical model of control object, and increases by pre-factor K fdesign, this is the principal feature of present patent application method for designing.
Summary of the invention
Controller in modern control system is all that appliance computer realizes, and each control system has certain control cycle (or sampling period) t s, t srelevant to the predominant frequency of computing machine.
Definition a: analog control system controls (operation) cycle t at it swhen determining, after steering order (being generally step signal) sends, system when not exceeding steering order value, the time T=nt of steering order value that quick and stable completes (reaching) s.Wherein T is system non-overshoot settling time, and n is the time metric of control system.
Design and the setting method of measuring the model-free PID controller of indignant thought based on system time of present patent application invention and design are exactly control cycle (or sampling period) t according to control system sset out.
PID controller parameter based on time metric is adjusted design procedure
Step one: do not comprise unstable limit and pure delay time in controlled device, in the control system built by PID controller, T integral time of PID controller i≈ (1/3 ~ 1/6) nt swherein n is the time metric of control system.If controlled device is G (s), if open cycle system Q (s)=G (s) G after PID compensates pID(s) ≈ 1/T is, then closed-loop system H (s)=1/ (T is+1), this system completes the time T ≈ 6T of steering order (step signal) i.
If the open cycle system after PID compensates (in most cases, the transport function of open cycle system can be of equal value therewith):
Q(s)=G(s)G PID(s)=1/(cs 2+ds)=1/(T is(T is/4+1))
Due to system non-overshoot, damping coefficientζ=1, then c=1/w 2, d=2/w, T i=d.Closed-loop system H (s)=1/ (cs 2+ ds+1), complete the time T ≈ 7/w of steering order, then T i≈ (2/7) T.
In formula, symbol description is as follows: T ≈ 7/w, w are system undamped vibration frequencies.
This is the basis of integral time in Patent design PID controller of the present invention.
The transport function of traditional PID controller is: G pID(s)=K p+ K i/ s+K ds
Step 2: the transport function of the PID controller that the present invention improves is:
G PID(s)=K f(K p+K i/s+K ds/(T ds+1))(1)
Wherein K f, T deffect be very important, select suitable K fthe effect of getting twice the result with half the effort can be obtained to the PID controller design of Large time delay process system.
K f=f (τ/t p)=f (h), h=τ/t p; Wherein, τ is system pure delay time, t pwhen without pure delay time τ, the settling time of system, K fthe monotonic decreasing function about h, 0 < K f≤ 1, (when τ=0, K f=1).The present invention obtains formula below by the method for a large amount of project data and curve:
K f=f(τ/t p)
K f=f(h)=c 0+c 1h+c 2h 2+c 3h 3+c 4h 4+…,(2)
Or K f = f ( h ) = e - ( c 0 + c 1 h + c 2 h 2 + c 3 h 3 + c 4 h 4 + ... ) - - - ( 3 )
Facts have proved, as coefficient h=τ/t pduring increase, in order to the stable of Guarantee control system and quality, K must be reduced f.
Calculating K above fformula non-linear, compensate for just by control object G (s)=k 0e -τ s/ (as 2+ bs+1)=e -τ sg 1non-minimum phase bit position e in (s) -τ sto minimum phase part k/ (as 2+ bs+1) impact.
Differential K in conventional PID controllers ds is improved to K ds/ (T ds+1), effectively can improve the percussive action of differential, pass through T dsetting can increase time of the differential action, reduce the interference that differential brings.
Step 3: when not exclusively knowing controlled device mathematical model, if set the time scale of control system as n, according to above-mentioned proposition, can according to method Tuning PID Controller parameter below.
In practical engineering application, at scattered control system (DCS) control cycle t swhen determining, this PID controller parameter setting method does not need the mathematical model knowing controlled device, only needs 1 basic parameter, i.e. time scale n of roughly certainty annuity.The time T that system completes steering order generally also easily obtains (T=nt at the scene s), which reflects control engineering teacher to the basic understanding of controlled system characteristic and understanding, such as in fuel-burning power plant in the systems such as steam temperature, water level, vapour pressure, load, T is different.
If can not n be estimated, just can only according to the control cycle of system (or sampling period) t sbe optimized design.
(1) if can estimate n, T integral time is set i=n it s, wherein n i=(1/3 ~ 1/6) n.If can not estimate n, T integral time is set i=(60 ~ 200) t s, initial setting pid parameter is as follows:
PID controller pre-factor K f=1, Proportional coefficient K p=(0.3 ~ 1.5); , differential coefficient K d=0, derivative time T d=(1 ~ 5) t s.General, there is large fluctuation in order to anti-locking system exports, need integral time selection larger, get T i=200t s, at Kp, T i, K dgo control system in this group parameter area, can ensure that controlled volume there will not be large fluctuation.Application golden section approach, can organize in parameter area at this and determine a K p, T i, K dvalue.
(2) apply the pid parameter of above-mentioned first step setting, if controlled device exports slow, application Fibonacci method reduces T integral time i, increase Proportional coefficient K pif (can n be estimated, arrange T i=nt s/ (4.5 ~) 6; ) or increase K f, K f=(1.1 ~ 1.7).
Apply the pid parameter of above-mentioned first step setting, if the output of controlled device is too fast and fluctuation is large, illustrate that controlled device comprises pure delay time link, or the damping coefficientζ < 0.76 of controlled device, can apply Fibonacci method increases T integral time i, reduce Proportional coefficient K p, (if can estimate n, T is set i=nt s/ (2.5 ~ 3.5)).If the output of controlled device still has fluctuation, but fluctuation reduces to some extent, can add differential K ds/ (T ds+1) (K is acted on d=T i/ 4, T d=(1 ~ 5) t s); If after adding the differential action, controlled device exports still has fluctuation, can reduce K f, get K f=0.1 ~ 0.7.
(3) if through (1), the design of (2) step still can not obtain satisfied Control platform, needs to get back to the first step and resets K p, T i, K dspan, repeat (1) and (2) step design until obtain be satisfied with Control platform.
According to above-mentioned (1), the step design tuning PID controller parameter of (2) and (3), this pid parameter setting method based on time scale, in (2) step, changes K fbe equivalent to and change other 3 parameter K simultaneously p, T i, K d, this is a kind of simple effective method, namely can obtain satisfied Control platform general 2 ~ 4 times.
Advantage and effect: based on the Optimization Design of the PID controller of time metric in a kind of Power Plant Thermal system of the present invention, its advantage measures design and the setting method of the model-free PID controller of indignant thought, it better can complete debug process faster, and can obtain satisfied Control platform.Apply golden section optimization method in program, with other apply golden section optimization method design PID controller parameter unlike, present patent application method for designing does not need by the mathematical model of control object, and increases by pre-factor K fdesign, this is the principal feature of present patent application method for designing.
Accompanying drawing explanation
Fig. 1 typical PID control system schematic diagram.
The pid parameter control system closed loop output response schematic diagram of Fig. 2 first time design tuning.
The pid parameter control system closed loop output response schematic diagram of Fig. 3 second time design tuning.
The control action schematic diagram of P, I, D in Fig. 4 PID controller.
Fig. 5 a, b are the program interface schematic diagram of design.
Fig. 6 is calculation procedure block scheme.
Embodiment:
The PID controller of adjusting designed by the present invention needs to realize in the upper configuration of the distribution type control system (DCS) of user, then carries out real time execution control.Also can implement on industrial computer.
According to the method for designing in summary of the invention, in conjunction with ratio P in PID controller, the effect (see Figure of description 4) of integration I and differential D, the PID controller parameter program of adjusting of the present invention's design is as shown below:
Apply golden section (0.618) optimization method in program, with other apply golden section optimization method design PID controller parameter unlike, this patent method for designing does not need by the mathematical model of control object, and increases pre-factor K fdesign, this is the principal feature of method for designing of the present invention.
So-called golden section (0.618) optimization method program is:
function[xo,fo]=Opt_Golden(f,a,b,TolX,TolFun,k)
%%%% golden search algorithms asks the optimum solution on interval [a, b]
%f is objective function, and TolX is x threshold value, and TolFun is function threshold, and k is iterations
R=(sqrt (5)-1)/2=0.618; %r is golden section point value,
H=b-a; % interval width
rh=r*h;
%%% gets c, d at 2, and calculates corresponding functional value fc and fd
c=b-rh;
d=a+rh;
fc=feval(f,c);
fd=feval(f,d);
%%% algorithm second step judges whether to stop iteration
%%%% algorithm the 3rd step, carries out new round iteration
Above in program Opt_Golden (f, a, b, TolX, TolFun, k)
Performance index f determines according to the computing method of analysis below
From perturbation analysis, when the disturbance of system is deterministic perturbation, the performance evaluation of its control system cannot directly apply minimum variance index.Usually under deterministic perturbation, the performance of system is divided into dynamic property and steady-state behaviour, is divided into again time domain performance and frequency domain performance from research field.The steady-state behaviour of time domain adopts steady-state error usually.The dynamic performance index of time domain is usually as follows:
Rise time: unit-step response rises to 90% of steady-state value from 10% of steady-state value, the required time.System for system unit step response monotone variation is like this, and to there being the system of concussion, can be defined as the time arriving stable state from initial value to first time.Rise time, reaction be the response speed of system.
Time to peak: for the system having concussion, its unit-step response is passed through final value and is reached first peak value time used.
Regulating time: unit-step response reaches the error band of 5% (2%) of steady-state value and the time of maintenance first.
Overshoot: when the unit step of system has overshoot, the peak-peak of its unit-step response deducts the percentage of the difference of steady-state value and the ratio of steady-state value.
Also a series of index is had to carry out the performance of reactive system in a frequency domain, as reflected the zero-frequency value of systematic steady state precision, the resonance peak of reflection system overshoot situation, the bandwidth of reflection system fast response characteristic and anti-interference filtration characteristic, the phase margin of reflection systematic steady state characteristic and magnitude margin etc.
Error intergal performance index
The various indexs introduced in upper joint are performance index of individual event, and the integrated form some integrated performance indexs being included to error has various forms, and conventional having is following several:
1) error intergal (IE)
I E = &Integral; 0 &infin; e ( t ) d t
2) absolute error index (IAE)
I A E = &Integral; 0 &infin; | e ( t ) | d t
3) integral square error (ISE)
I S E = &Integral; 0 &infin; e 2 ( t ) d t
4) time and Error Absolute Value integration (ITAE)
I T A E = &Integral; 0 &infin; t | e ( t | ) d t
The estimation of above-mentioned various index be situation for whole process but emphasis is had nothing in common with each other.The error gentle (SSE) of more brief and practical in actual performance is evaluated:
S S E = 1 n &Sigma; i = 1 n &epsiv; 2 i
Can the determinacy performance index of define system based on this:
&eta; d = &Sigma; i = 0 d - 1 e 2 ( i ) &Sigma; i = 0 &infin; e 2 ( i )
The Unified Form of performance index
No matter what exist in loop is deterministic perturbation or randomness disturbance, B.Huang is unified H 2performance index under norm meaning:
&eta; = m i n ( | | G c l | | 2 2 ) | | G c l | | 2 2
When system only has randomness disturbance, assuming that the setting value of system is 0, namely as randomness disturbance α t≠ 0, deterministic perturbation time performance evaluation as follows:
Its closed loop H 2norm is defined as follows:
| | G c l | | 2 2 = 1 2 &pi; &Integral; &pi; &pi; t r &lsqb; G c l ( e j w ) G c l T ( e - j w ) &rsqb; d w
Wherein: G clfor control system closed loop transfer function,
Randomness disturbance is write as polynomial form to the transport function of departure:
G c l = &Sigma; i = 0 &infin; f &alpha; ( i ) z - i
Now the stochastic performance index of system is defined as follows:
&eta; s = m i n ( | | G c l | | 2 2 ) | | G c l | | 2 2 = &Sigma; i = 0 d - 1 f &alpha; 2 ( i ) &Sigma; i = 0 &infin; f &alpha; 2 ( i ) = &delta; m v 2 &delta; y
η smore close to 1, think that the performance of system is better.The computing method of its stochastic performance index can be drawn from formula right-hand component.Estimation for minimum variance can have been come by FCOR algorithm or ARIMA modeling.
When system only has setting value disturbance, performance evaluation being carried out to it and it is only there being the evaluation in randomness disturbance situation similar, starting with from departure, defining its tracking performance index.Now, α t=0 formula is easy to get:
Now, its closed loop H 2norm is defined as follows:
| | G | | 2 2 = G S E E = 1 n &Sigma; i = 1 n &epsiv; 2 i
Tracking performance (determinacy performance) index now [44]be defined as follows:
&eta; d = m i n ( | | G c l | | 2 2 ) | | G c l | | 2 2 = &Sigma; i = 0 d - 1 f &zeta; 2 ( i ) &Sigma; i = 0 &infin; f &zeta; 2 ( i ) = &Sigma; i = 0 d - 1 &epsiv; 2 ( i ) &Sigma; i = 0 &infin; &epsiv; 2 ( i )
Easily find out from the definition of this index, step signal on the perfect tracking that system is faster after dead band, then η dmore close to 1, the tracking performance of system is better.
Other performance index
The Performance Evaluating Indexes of control system, except minimum variance, also has other indexs many.
First based on the index of historical data be easily expect and also practicality.Operations staff is by oneself observation and a certain control system of empirical discovery, control effects within certain a period of time is ideal, so based on the service data in this period, calculate value (as variance) now as the benchmark to this system evaluation according to a certain algorithm.In this can be based on the performance index of historical data:
&eta; h i s = J h i s J a c t
But certainly close to 1, this index more thinks that the performance of system performance is more satisfactory.Its shortcoming is that this desired value is limited between [0,1] unlike minimum variance evaluation index, and its scope is indefinite.But its advantage is also clearly, this index overcomes minimum variance performance index and makes control action too fierce in some cases and the shortcoming that cannot realize, because must reach based on the performance number of historical data, and calculating is easily, and easy to understand.
Secondly, the minimum variance (GeneralisedMinimumVariance, GMV) of broad sense is also a conventional index.Calculate MV index time, calculating be system export y tor error e tvariance under control law, and this control law usually causes the control action of system too drastic and cannot realize LMS control in practice.For overcoming this shortcoming, GMV is by controlled quentity controlled variable u twith error e tbe weighted structure new variables, control law is by realizing this new variables least squares optimization.By the selection of the weighting coefficient to new variables, consideration is carried out to the control action of system and avoids control procedure too drastic, thus the control law obtained can be used in practice.Designed program interface is as Fig. 5 a, Fig. 5 b, and Fig. 6 is shown in by calculation procedure block scheme.
The present invention's design (increases pre-factor K artificial experience fwith the estimation to integral time) combine with golden section (0.618) optimization method, devise control system below.
Figure of description 1 is the typical PID control system block diagram designed in Power Plant Thermal process.
Under coal-supplying amount and the indeclinable situation of absorbing quantity, to be the transport function between cigarette vapour oxygen level and air output be controlled device:
G(s)=1.3e -5s/(54s 2+14s+1);
Oxygen amount set-point is 4.3 (%), control cycle t s=0.2 (s);
According to calculation procedure above, process of optimization is as follows:
1. first step parameter designing: get K f=1; t i=80t s=16; k p=0.62; k d=0.0; Td=1;
Figure of description 2 is application methods for designing of the present invention, the pid parameter control system closed loop output response of first step design tuning.
2. as can be seen from Figure 2, the first time pid parameter output of design tuning has overshoot, second time parameter designing: get K f=0.55; Other parameter constant, control system closed loop output response as represented in fig. 3.Second time parameter designing value control effects is fine.
Fig. 4 is the control action of P, I, D and PID in PID controller.
Fig. 5 a, b are the program interface schematic diagram of design, and Fig. 6 is calculation procedure block scheme.

Claims (1)

1. in Power Plant Thermal system based on the Optimization Design of the PID controller of time metric, it is characterized in that: the method concrete steps are as follows:
Step one: do not comprise unstable limit and pure delay time in controlled device, in the control system built by PID controller, T integral time of PID controller i≈ (1/3 ~ 1/6) nt swherein n is the time metric of control system; If controlled device is G (s), if open cycle system Q (s)=G (s) G after PID compensates pID(s) ≈ 1/T is, then closed-loop system H (s)=1/ (T is+1), this system completes the time T ≈ 6T of steering order and step signal i;
If the open cycle system after PID compensates is
Q(s)=G(s)G PID(s)=1/(cs 2+ds)=1/(T is(T is/4+1))
Due to system non-overshoot, damping coefficientζ=1, then c=1/w 2, d=2/w, T i=d; Closed-loop system H (s)=1/ (cs 2+ ds+1), complete the time T ≈ 7/w of steering order, then T i≈ (2/7) T;
In formula, symbol description is as follows: T ≈ 7/w, w are system undamped vibration frequencies;
This is the basis of integral time in design PID controller, and the transport function of traditional PID controller is:
G PID(s)=K p+K i/s+K ds;
Step 2: the transport function of the PID controller of improvement is:
G PID(s)=K f(K p+K i/s+K ds/(T ds+1))(1)
Wherein K f, T dimportant role, select suitable K fvery crucial to the PID controller design of Large time delay process system;
K f=f (τ/t p)=f (h), h=τ/t p; Wherein, τ is system pure delay time, t pwhen without pure delay time τ, the settling time of system, K fthe monotonic decreasing function about h, 0 < K f≤ 1, when τ=0, K f=1; Formula by the method for a large amount of project data and curve obtains below:
K f=f(τ/t p)
K f=f(h)=c 0+c 1h+c 2h 2+c 3h 3+c 4h 4+…,(2)
Or K f = f ( h ) = e - ( c 0 + c 1 h + c 2 h 2 + c 3 h 3 + c 4 h 4 + ... ) - - - ( 3 )
Facts have proved, as coefficient h=τ/t pduring increase, in order to the stable of Guarantee control system and quality, K must be reduced f;
Calculating K above fformula non-linear, compensate for just by control object G (s)=k 0e -τ s/ (as 2+ bs+1)=e -τ sg 1non-minimum phase bit position e in (s) -τ sto minimum phase part k/ (as 2+ bs+1) impact;
Differential K in conventional PID controllers ds is improved to K ds/ (T ds+1), effectively can improve the percussive action of differential, pass through T dsetting can increase time of the differential action, reduce the interference that differential brings;
Step 3: when not exclusively knowing controlled device mathematical model, if set the time scale of control system as n, according to above-mentioned proposition, according to method Tuning PID Controller parameter below;
In practical engineering application, at distributed monitoring control system control cycle t swhen determining, this PID controller parameter setting method does not need the mathematical model knowing controlled device, only needs 1 basic parameter, i.e. time scale n of roughly certainty annuity; The time T that system completes steering order generally also easily obtains T=nt at the scene s, which reflects control engineering teacher to the basic understanding of controlled system characteristic and understanding, in fuel-burning power plant in steam temperature, water level, vapour pressure, load system, T is different;
If can not n be estimated, just can only according to the control cycle of system or sampling period t sbe optimized design;
(1) if can estimate n, T integral time is set i=n it s, wherein n i=(1/3 ~ 1/6) n; If can not estimate n, T integral time is set i=(60 ~ 200) t s, initial setting pid parameter is as follows:
PID controller pre-factor K f=1, Proportional coefficient K p=(0.3 ~ 1.5); , differential coefficient K d=0, derivative time T d=(1 ~ 5) t s; General, there is large fluctuation in order to anti-locking system exports, need integral time selection larger, get T i=200t s, at Kp, T i, K dgo control system in this group parameter area, can ensure that controlled volume there will not be large fluctuation; Application golden section approach, determines a K in this group parameter area p, T i, K dvalue;
(2) apply the pid parameter of above-mentioned first step setting, if controlled device exports slow, application Fibonacci method reduces T integral time i, increase Proportional coefficient K pif can n be estimated, arrange T i=nt s/ (4.5 ~ 6) or increase K f, K f=(1.1 ~ 1.7);
Apply the pid parameter of above-mentioned first step setting, if the output of controlled device is too fast and fluctuation is large, illustrate that controlled device comprises pure delay time link, or the damping coefficientζ < 0.76 of controlled device, application Fibonacci method increases T integral time i, reduce Proportional coefficient K pif can n be estimated, arrange T i=nt s/ (2.5 ~ 3.5)); If the output of controlled device still has fluctuation, but fluctuation reduces to some extent, adds differential K ds/ (T ds+1) K is acted on d=T i/ 4, T d=(1 ~ 5) t s; If after adding the differential action, controlled device exports still has fluctuation, reduces K f, get K f=0.1 ~ 0.7;
(3) if through (1), the design of (2) step still can not obtain satisfied Control platform, needs to get back to the first step and resets K p, T i, K dspan, repeat (1) and (2) step design until obtain be satisfied with Control platform;
According to above-mentioned (1), the design tuning PID controller parameter of (2) and (3) step, this pid parameter setting method based on time scale, in (2) step, changes K fbe equivalent to and change other 3 parameter K simultaneously p, T i, K d, this is a kind of simple effective method, can obtain satisfied Control platform general 2 ~ 4 times.
CN201510404996.7A 2015-07-10 2015-07-10 The Optimization Design of PID controller based on time metric in Power Plant Thermal system Active CN105045233B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510404996.7A CN105045233B (en) 2015-07-10 2015-07-10 The Optimization Design of PID controller based on time metric in Power Plant Thermal system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510404996.7A CN105045233B (en) 2015-07-10 2015-07-10 The Optimization Design of PID controller based on time metric in Power Plant Thermal system

Publications (2)

Publication Number Publication Date
CN105045233A true CN105045233A (en) 2015-11-11
CN105045233B CN105045233B (en) 2017-09-19

Family

ID=54451837

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510404996.7A Active CN105045233B (en) 2015-07-10 2015-07-10 The Optimization Design of PID controller based on time metric in Power Plant Thermal system

Country Status (1)

Country Link
CN (1) CN105045233B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105487375A (en) * 2015-12-31 2016-04-13 中国船舶重工集团公司第七一九研究所 Discrete PID controller parameter setting method
CN105700353A (en) * 2016-01-30 2016-06-22 河南城建学院 A PID controller parameter optimal setting method based on a differential evolution method
CN106054610A (en) * 2016-06-23 2016-10-26 东南大学 Optimized PI (Proportional Integral) controller parameter engineering setting method
CN106094510A (en) * 2016-06-30 2016-11-09 电子科技大学 A kind of pid parameter control method based on interference inverter
CN107361760A (en) * 2017-07-12 2017-11-21 中国科学院上海微系统与信息技术研究所 The magnetocardiograph of magnetocardiograph diagnostic system and the application system
CN110083877A (en) * 2019-03-29 2019-08-02 北京国电龙源环保工程有限公司 A kind of fluctuation heat power engineering system transfer function modeling method
CN113110034A (en) * 2021-05-10 2021-07-13 常州市新港热电有限公司 DCS-based fuzzy PID control system for induced draft fan
CN113139291A (en) * 2021-04-23 2021-07-20 广东电网有限责任公司电力科学研究院 Method and device for obtaining optimal sliding window filtering model of controlled process
CN113325692A (en) * 2021-04-28 2021-08-31 东南大学 PID controller pull-back type setting method based on neighbor equivalence
CN113708952A (en) * 2020-10-22 2021-11-26 天翼智慧家庭科技有限公司 Method and device for guaranteeing web service connection stability

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10302585A1 (en) * 2003-01-22 2004-08-12 Endress + Hauser Wetzer Gmbh + Co Kg Optimization of PID regulation parameters of a PID regulator by setting minimum control values using minimum values in the proportional regions of control values
US20050065621A1 (en) * 2000-06-20 2005-03-24 Danyang Liu Methods of designing optimal linear controllers
CN101069902A (en) * 2007-06-13 2007-11-14 湖南大学 High-strength aluminium, magnesium alloy equi-temperature extrusion method based on numerical simulation
JP2012190364A (en) * 2011-03-11 2012-10-04 Hiroshima Univ Control parameter adjustment device and control parameter adjustment method
CN103439880A (en) * 2013-08-30 2013-12-11 上海电力学院 PID parameter setting method based on MCP standard transfer function
CN104360595A (en) * 2014-11-19 2015-02-18 上海电力学院 Method for MCP-PID parameter setting of delay-free inertia process

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050065621A1 (en) * 2000-06-20 2005-03-24 Danyang Liu Methods of designing optimal linear controllers
DE10302585A1 (en) * 2003-01-22 2004-08-12 Endress + Hauser Wetzer Gmbh + Co Kg Optimization of PID regulation parameters of a PID regulator by setting minimum control values using minimum values in the proportional regions of control values
CN101069902A (en) * 2007-06-13 2007-11-14 湖南大学 High-strength aluminium, magnesium alloy equi-temperature extrusion method based on numerical simulation
JP2012190364A (en) * 2011-03-11 2012-10-04 Hiroshima Univ Control parameter adjustment device and control parameter adjustment method
CN103439880A (en) * 2013-08-30 2013-12-11 上海电力学院 PID parameter setting method based on MCP standard transfer function
CN104360595A (en) * 2014-11-19 2015-02-18 上海电力学院 Method for MCP-PID parameter setting of delay-free inertia process

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
何芝强: "PID控制器参数整定方法及其应用研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 *
刘聪: "锅炉快速控制原型开发及PID参数自整定方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
刘邦利 等: "《电厂化学仪表及自动化》", 28 February 2015, 中国电力出版社 *
朱晓青: "《过程检测控制技术与应用》", 30 June 2002, 冶金工业出版社 *
李国林: "PID控制器参数整定技术研究与优化设计", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105487375B (en) * 2015-12-31 2018-06-22 中国船舶重工集团公司第七一九研究所 A kind of Discrete PID Controller Parameters setting method
CN105487375A (en) * 2015-12-31 2016-04-13 中国船舶重工集团公司第七一九研究所 Discrete PID controller parameter setting method
CN105700353A (en) * 2016-01-30 2016-06-22 河南城建学院 A PID controller parameter optimal setting method based on a differential evolution method
CN106054610A (en) * 2016-06-23 2016-10-26 东南大学 Optimized PI (Proportional Integral) controller parameter engineering setting method
CN106054610B (en) * 2016-06-23 2019-01-25 东南大学 A kind of PI controller parameter practical tuning method of optimization
CN106094510A (en) * 2016-06-30 2016-11-09 电子科技大学 A kind of pid parameter control method based on interference inverter
CN106094510B (en) * 2016-06-30 2019-11-05 电子科技大学 A kind of pid parameter adjusting method based on interference inverter
CN107361760A (en) * 2017-07-12 2017-11-21 中国科学院上海微系统与信息技术研究所 The magnetocardiograph of magnetocardiograph diagnostic system and the application system
CN110083877A (en) * 2019-03-29 2019-08-02 北京国电龙源环保工程有限公司 A kind of fluctuation heat power engineering system transfer function modeling method
CN113708952A (en) * 2020-10-22 2021-11-26 天翼智慧家庭科技有限公司 Method and device for guaranteeing web service connection stability
CN113708952B (en) * 2020-10-22 2024-06-07 天翼数字生活科技有限公司 Method and device for guaranteeing web service connection stability
CN113139291A (en) * 2021-04-23 2021-07-20 广东电网有限责任公司电力科学研究院 Method and device for obtaining optimal sliding window filtering model of controlled process
CN113325692A (en) * 2021-04-28 2021-08-31 东南大学 PID controller pull-back type setting method based on neighbor equivalence
CN113325692B (en) * 2021-04-28 2022-12-27 东南大学 PID controller pull-back type setting method based on neighbor equivalence
CN113110034A (en) * 2021-05-10 2021-07-13 常州市新港热电有限公司 DCS-based fuzzy PID control system for induced draft fan

Also Published As

Publication number Publication date
CN105045233B (en) 2017-09-19

Similar Documents

Publication Publication Date Title
CN105045233A (en) Optimum design method for PID (Proportion Integration Differentiation) controller based on time dimension in heat-engine plant thermal system
JP5005695B2 (en) Adaptive multi-variable MPC controller
CN104238366B (en) The forecast Control Algorithm and device of piezoelectric ceramic actuator based on neuroid
CN113031565B (en) Method and system for predicting primary frequency modulation power response performance of thermal power generating unit
CN100511042C (en) Groove type reactor non-linear predication control method based on multi-kernel support vector machine
CN107450310B (en) Set value excitation closed-loop identification method of continuous process model
CN106530118A (en) Constrained GPC optimal control method for load frequency of multi-zone interconnection power system
Chertovskikh et al. An adaptive PID controller with an online auto-tuning by a pretrained neural network
Anusha et al. Comparison of tuning methods of PID controller
Folorunso et al. Control of a two layered coupled tank: Application of IMC, IMC-PI and pole-placement PI controllers
CN104808486A (en) Predictive control method and predictive control device for piezoelectric ceramic actuator based on fuzzy TS (Takagi-Sugeno) model
Deshpande et al. Design of multi scale PID controller for Temperature process
Yang et al. Application of fuzzy neural network PID algorithm in oil pump control
CN105607486A (en) General control method for series chemical integration time lag process
Škarda et al. Bode-like control loop performance index evaluated for a class of fractional-order processes
Zhang et al. Nonlinear model predictive control of ultra-supercritical once through boiler-turbine unit
Grelewicz et al. Practical Verification of the Advanced Control Algorithms Based on the Virtual Commissioning Methodology-A Case Study
Campo PID control design
Ionesi et al. On-line parameter and state estimation of an air handling unit model: experimental results using the modulating function method
Dulău et al. Conventional control vs. robust control on heat-exchangers
Ruusu et al. Sliding mode SISO control of model parameters for implicit dynamic feedback estimation of industrial tracking simulation systems
Aravind et al. Reliability of pid controller for a real time sensor interfaced temperature process
Qingli et al. Control system of the boiler water level based on normalized PID
Gao et al. Control of system with large parametric uncertainty using multiple robust controllers and switching
Mhankale et al. Robust control analysis using μ-synthesis for nonlinear control system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant