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 PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000013461 design Methods 0.000 title claims abstract description 41
- 230000010354 integration Effects 0.000 title abstract description 4
- 230000004069 differentiation Effects 0.000 title abstract 2
- 238000013178 mathematical model Methods 0.000 claims abstract description 9
- 238000005457 optimization Methods 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 9
- 230000009471 action Effects 0.000 claims description 7
- 238000013016 damping Methods 0.000 claims description 4
- 230000009699 differential effect Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 2
- 230000008859 change Effects 0.000 claims description 2
- 230000003247 decreasing effect Effects 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 2
- 238000012544 monitoring process Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 6
- 238000012546 transfer Methods 0.000 abstract description 2
- 230000004044 response Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 230000009514 concussion Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 206010000372 Accident at work Diseases 0.000 description 1
- 101100502522 Mus musculus Fcor gene Proteins 0.000 description 1
- 235000019504 cigarettes Nutrition 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- YHXISWVBGDMDLQ-UHFFFAOYSA-N moclobemide Chemical compound C1=CC(Cl)=CC=C1C(=O)NCCN1CCOCC1 YHXISWVBGDMDLQ-UHFFFAOYSA-N 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/418—Total 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/41885—Total 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
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Engineering & Computer Science (AREA)
- Feedback Control In General (AREA)
- Manufacturing & Machinery (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
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
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)
2) absolute error index (IAE)
3) integral square error (ISE)
4) time and Error Absolute Value integration (ITAE)
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:
Can the determinacy performance index of define system based on this:
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:
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:
Wherein: G
clfor control system closed loop transfer function,
Randomness disturbance is write as polynomial form to the transport function of departure:
Now the stochastic performance index of system is defined as follows:
η
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:
Tracking performance (determinacy performance) index now
[44]be defined as follows:
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:
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
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.
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)
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)
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 |
-
2015
- 2015-07-10 CN CN201510404996.7A patent/CN105045233B/en active Active
Patent Citations (6)
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)
Title |
---|
何芝强: "PID控制器参数整定方法及其应用研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 * |
刘聪: "锅炉快速控制原型开发及PID参数自整定方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
刘邦利 等: "《电厂化学仪表及自动化》", 28 February 2015, 中国电力出版社 * |
朱晓青: "《过程检测控制技术与应用》", 30 June 2002, 冶金工业出版社 * |
李国林: "PID控制器参数整定技术研究与优化设计", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (15)
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 | |
CN103558755B (en) | Fractional order integration PID controller setting and self-setting method | |
CN103676651B (en) | Steam Temperature for Boiler forecast Control Algorithm based on state observation model | |
CN103984234A (en) | Electro hydraulic servo system self-correction fuzzy PID control method | |
CN103406909B (en) | Tracking control device and method of mechanical arm system | |
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 | |
CN105867138A (en) | Stable platform control method and device based on PID controller | |
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 | |
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 |