CN110262223A - A kind of comprehensive model modelling approach of hydraulic turbine based on Fractional Order PID speed-regulating system - Google Patents

A kind of comprehensive model modelling approach of hydraulic turbine based on Fractional Order PID speed-regulating system Download PDF

Info

Publication number
CN110262223A
CN110262223A CN201910643191.6A CN201910643191A CN110262223A CN 110262223 A CN110262223 A CN 110262223A CN 201910643191 A CN201910643191 A CN 201910643191A CN 110262223 A CN110262223 A CN 110262223A
Authority
CN
China
Prior art keywords
regulating system
hydraulic turbine
model
drosophila
pid speed
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
CN201910643191.6A
Other languages
Chinese (zh)
Other versions
CN110262223B (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.)
North China University of Water Resources and Electric Power
Original Assignee
North China University of Water Resources and Electric Power
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 North China University of Water Resources and Electric Power filed Critical North China University of Water Resources and Electric Power
Priority to CN201910643191.6A priority Critical patent/CN110262223B/en
Publication of CN110262223A publication Critical patent/CN110262223A/en
Application granted granted Critical
Publication of CN110262223B publication Critical patent/CN110262223B/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
    • G05B11/00Automatic controllers
    • G05B11/01Automatic controllers electric
    • G05B11/36Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • G05B11/42Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P.I., P.I.D.
    • 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
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/20Hydro energy

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

The present invention provides a kind of comprehensive model modelling approach of the hydraulic turbine of adaptive fractional rank PID speed-regulating system, the following steps are included: step 1: establishing the hydrogovernor simulation analysis mathematical model based on Fractional Order PID speed-regulating system, and optimize the objective function based on Fractional Order PID speed-regulating system using hybrid algorithm;Step 2: establishing the machinery hydraulic system and diversion system mathematical model of Turbine Governor System;Step 3: establishing hydraulic turbine generator model;The present invention utilizes MATLAB simulation modeling, classical PID speed-regulating system is improved on the basis of forefathers' research, according to instance analysis, it establishes and is able to reflect Turbine Governor System and the comprehensive model of generator system parameters, in practical applications to the operating condition of practical water turbine set, the variation of hydraulic turbine parameters has good reaction when in face of load fluctuation, provides safeguard for accident forecast with system safety operation.

Description

A kind of comprehensive model modelling approach of hydraulic turbine based on Fractional Order PID speed-regulating system
Technical field:
The present invention relates to power station mechanics, Turbine Governor System principles, and in particular to one kind is adjusted the speed based on Fractional Order PID The comprehensive model modelling approach of the hydraulic turbine of system.
Background technique
Traditional Turbine Governor System control structure generally uses classical Turbine Governor System, i.e. parallel connection PID speed regulation System control, has the advantages such as simple to operation, with the development of electric system, to the stability requirement of Hydropower Unit also not Disconnected to improve, still, the hydraulic turbine system of PID governor system control in parallel has the following problems:
1. stablizing when starting slow;
2. occurring that stability is poor when system loading fluctuation, in face of the emergency event capacity of self-regulation deficiency etc. of electric system Problem.
The shortcomings that in face of classical Turbine Governor System, a large amount of scholars propose advanced intelligentized control method mode, such as mould PID control, BP neural network control, chaotic particle group control etc. are pasted, these researchs advanced optimize Turbine Governor System, There is the speed adjusting performance of the hydraulic turbine to be obviously improved.
In view of the control mode of the domestic hydraulic turbine mostly uses the adjusting based on PID governor system control in parallel control Mode brings smaller impact, the speed regulation system of the hydraulic turbine to electric system as far as possible in order to quick response system requirement The optimization of system just seems most important.
There is remarkable progress to the Nonlinear Mechanism Modeling Research of the hydraulic turbine in terms of simulation modeling, but has not divided deeply Analyse influence of the load fluctuation to Turbine Governor System;Load and generator model are generally equivalent to the transmitting letter of single order at present Number, the practical problem limitation of reflection is bigger, for hydraulic turbine generator spy on the basis of reflecting water turbine governing performance When property, dynamic load model and electric power system fault in the influence of Turbine Governor System up for further analyzing.
Summary of the invention
The hydraulic turbine is each when the present invention is to solve not reflecting dynamic load in the Nonlinear Mechanism modeling of the hydraulic turbine The problem of item important parameter variation, propose a kind of comprehensive model modelling approach of the hydraulic turbine based on Fractional Order PID speed-regulating system.
In order to achieve the above object, the invention adopts the following technical scheme:
A kind of comprehensive model modelling approach of the hydraulic turbine of adaptive fractional rank PID speed-regulating system, comprising the following steps:
Step 1: establishing the hydrogovernor simulation analysis mathematical model based on Fractional Order PID speed-regulating system, and utilize The hybrid algorithm optimization objective function based on Fractional Order PID speed-regulating system;
Step 2: establishing the machinery hydraulic system and diversion system mathematical model of Turbine Governor System;
Step 3: hydraulic turbine generator model is established, by traditional single order hydraulic turbine generator model extension to 5 rank water wheels Machine model, and thus introduce exciter model and dynamic load model;
The step 1 specifically includes the following steps:
Step 1.1: Fractional Order PID speed-regulating system is obtained by classical PID speed-regulating system, method particularly includes:
The differential expressions of classical PID speed-regulating system are as follows:
U (t)=kpe(t)+kiDe(t)+kdDe(t);
Wherein, kp indicates scale parameter, and setting range is 0.5~20;Ki indicates integral parameter, setting range 0.05s- 1~10s-1, kd indicate differential parameter, and setting range is 0~5s;
Tri- adjustment parameters of kp, ki, kd are replaced with into transient state slip coefficient bt, damping time constant Td, acceleration respectively Time constant Tn, specific:
The transmission function of classical PID speed-regulating system are as follows:
Wherein, bt is transient state slip coefficient, takes 0~1.0;Td is damping time constant, and setting range is 2s~20s, Tn For G-time constant, setting range is that 0~2s, pc and P are respectively that the power of the assembling unit is given and the power of the assembling unit, ycDistinguish with y For guide vane opening is given and guide vane degree, fc and f are respectively that frequency is given and unit frequency, Δ f are frequency difference, and Δ f' is through overfrequency Dead zone efFrequency departure afterwards;
The differential expressions of above-mentioned classical PID speed-regulating system are transformed to Fractional Order PID speed-regulating system by Laplace Differential expressions:
U (t)=kpe(t)+kiDαe(t)+kdDμe(t); (4)
Wherein α, μ > 0, α are integral order, and μ is differential order;
The transmission function of Fractional Order PID speed-regulating system are as follows:
Step 1.2: the objective function of Fractional Order PID speed-regulating system is defined, method particularly includes:
It is allocated as using double fault difference-product as the objective function of the Fractional Order PID speed-regulating system described in step 1.1:
Wherein, e (t) indicates the deviation of reality output and desired output, and t is time, ISEIndicate square deviation integral;
Step 1.3: using the objective function based on Fractional Order PID speed-regulating system described in hybrid algorithm Optimization Steps 1.2, Method particularly includes:
By kp、ki、kd, five parameters of α, μ form initial drosophila positions, so that drosophila is carried out optimizing by flavor concentration, utilize BP The deep search characteristic of neural network is to kp、ki、kd, five parameter optimizations of α, μ solve, specifically includes the following steps:
Step 1.3.1: initialization drosophila population scale groupsize (400) and maximum number of iterations maxnum (400), Drosophila population position (X is randomly generatedaxis, Yaxis), iteration step value R is defined as (0.85-1);
Step 1.3.2: randomization drosophila population position and direction, according to drosophila population position and initial point distance DDistDetermine Drosophila population flavor concentration value Si
Step 1.3.3: the flavor concentration of drosophila individual is calculated, by drosophila population flavor concentration value SiFlavor concentration is substituted into sentence Determine function FfunctionIn, and flavor concentration optimum individual is found from each drosophila group;
Retain optimal flavor concentration bbestsmellWith corresponding (Xi, Yi) position, drosophila group SsmellThe coordinate is flown to, bbest indesIndicate optimal drosophila position;
Step 1.3.4: the drosophila position coordinates (X that step 1.3.3 is obtainedi, Yi) input BP neural network hidden layer nnet (2)(k), input and output are defined;
Wherein, i=1,2 ..., Q define n, Q value, Z according to controlled device complexityj (2)For drosophila population after optimization (Xi, Yi) position, wij (2)For hidden layer weighting coefficient, Zi(k)(2)Take activation functions Sigmoid function;
Step 1.3.5: BP neural network output layer n is definednet (3)(k), drosophila individual is updated, chooses E (k) as performance mistake Poor index, while retaining the position coordinates (X of the most dense drosophila individual of tastei, Yi);
Step 1.3.6: the performance indicator calculated by judging previous step exports optimal ginseng when reaching maximum number of iterations Number kp,ki,kd, α, μ, otherwise the process of end is transferred to step 1.3.2, until while reaching maximum number of iterations exports optimized parameter kp,ki,kd, α, μ terminate process.
The step 2 specifically includes the following steps:
Step 2.1: the mathematical model of the machinery hydraulic system of Turbine Governor System is established, following methods are specifically used:
The effect of machinery hydraulic system is transform electric and to zoom into the mechanical displacement letter with certain operating force Number;Due to second level servomotor responsive time constant Ty1Much smaller than main servomotor responsive time constant Ty, usually by it in modeling It is reduced to an inertial element, wherein y is machinery hydraulic system output signal, may also indicate that guide vane opening relative size, according to Water-Turbine Engine Adaption characteristic joined frequency dead band and saturation limiting element, main control valve and servomotor linear segment transmission function Are as follows:
Step 2.2: the mathematical model of the diversion system of Turbine Governor System is established, following methods are specifically used:
According to the difference of diversion penstock length, on simulation modeling, rigidity testing machine model is used in 700m or less, is greater than Elastic water hammer is used when 700m, water flow waterpower change procedure is described by following equation in diversion system:
The equation of motion:
Wherein H is head, and Q is water flow, and S is waterpipe cross-sectional area, and t is the time, and g is acceleration of gravity;
Flow equation:
Wherein a is water flow acceleration;
N=0 and n=1 is taken to obtain rigidity testing machine G through Taylor series expansion water attack equation in pressure diversion pipelineA1With Elastic water hammer GA2Transmission function:
Wherein TWFor rigid water attack time constant, Tr=2L/V is elastic water attack time constant, and L is pressure diversion conduit deferent Road length, V are water flow velocity of wave, H0For head, Q0For flow, V0For water flow base speed.
The step 3 specifically includes the following steps:
Hydraulic turbine generator model is established, by traditional single order hydraulic turbine generator model extension to 5 rank hydraulic turbine models, And thus introduce exciter model and dynamic load model;Method particularly includes:
Step 3.1: the first order modeling that traditional hydraulic turbine generator model uses are as follows:
Wherein, f is unit frequency, and p is unit input power difference, and Tn is unit set inertia time constant, general TnFor 3s ~12s, enFor unit static frequency self-regulation coefficient, enTake 0.5~2.0s;
Above-mentioned first order modeling can not reflect actual Parameters of Hydro-generator and parameters of excitation system, herein by The Synchronous Machine Models that MATLAB/Simulink emulation platform establishes five rank state-space equations are equivalent to hydrogenerator mould Type defines the basic parameter characteristic of the hydraulic turbine, it is contemplated that the dynamic characteristic in stator and rotor magnetic field and Damper Winding, model etc. Imitating circuit indicates in rotor reference system (dq frame).
Wherein, stator voltage equation are as follows:
Wherein UdFor stator d shaft voltage, UqFor stator q shaft voltage, φdAnd φqIndicate dp axis magnetic linkage, idAnd iqFor equivalent dq Electric current under coordinate system, raFor equivalent resistance, dq axis time transient internal voltage is E "dWith E "q, X "dAnd X "qFor dq subtranient reactance;
Rotor f winding, dq winding voltage equation and equation of rotor motion are as follows:
Wherein d axis and q axis time constant (all as unit of s), d axis transient state is opened a way (Tdo') or short-circuit (Td') time Constant, d axis subtransient is opened a way (Tdo ") or short-circuit (Td ") time constant, and q axis transient state is opened a way (Tqo') or short-circuit (Tq') time Constant (only limits round rotor), and q axis subtransient is opened a way (Tqo ") or short-circuit (Tq ") time constant.E'qFor transient internal voltage, EeFor Dq winding electromotive force, X 'dFor d axis transient state reactance, XdAnd XqFor dq axis reactance, W is rotor machinery angular speed, TjFor generating set Inertia time constant, TmFor prime mover machine torque,To act on the uneven torque on armature spindle.
The invention has the benefit that
Hydrogovernor simulation analysis mathematical model of the foundation based on Fractional Order PID speed-regulating system of the present invention, Using MATLAB simulation modeling, classical PID speed-regulating system is improved on the basis of forefathers' research, according to example point Analysis, it is established that Turbine Governor System and the comprehensive model of generator system parameters are able to reflect, for hydraulic turbine machine Pool-size, frequency, rotor windings, stator winding, excitation winding, number of magnetic pole pairs, output voltage current conditions, excitation voltage feelings Condition, load variations situation can be made a concrete analysis of, and temporary stable state when all kinds of power station actual motions can be largely emulated Variation compensates for the defect of traditional hydraulic turbine modeling;In practical applications to the operating condition of practical water turbine set, load is faced The variation of hydraulic turbine parameters has good reaction when fluctuation, provides safeguard for accident forecast with system safety operation.
Detailed description of the invention:
Fig. 1 is flow chart of the method for the present invention;
Fig. 2 is hydrogovernor simulation analysis mathematics of the foundation based on Fractional Order PID speed-regulating system of the present invention Model, and utilize the method flow diagram of the hybrid algorithm optimization objective function based on Fractional Order PID speed-regulating system;
Fig. 3 is the target letter based on Fractional Order PID speed-regulating system described using hybrid algorithm optimization of the present invention Several method flow diagrams;
Classical PID speed-regulating system and Fractional Order PID governor system control Comparative result schematic diagram when Fig. 4 is load disturbance;
FOA-PID speed-regulating system and BP-PID governor system control Comparative result schematic diagram when Fig. 5 is load disturbance;
BP-PID speed-regulating system and BPFOA-FOPID governor system control Comparative result schematic diagram when Fig. 6 is load disturbance;
Hydraulic turbine excitation voltage, active, rotation speed change when Fig. 7 loads switching;
Three-phase current changes when Fig. 8 is load investment;
Fig. 9 is that three-phase current changes when loading section is cut off;
Figure 10 is BPFOA-PID double objective optimal control Comparative result schematic diagram;
Figure 11 is the schematic diagram of the adaptive optimizing result of BPFOA-FOPID.
Specific embodiment:
A kind of comprehensive model of the hydraulic turbine of adaptive fractional rank PID speed-regulating system of the present invention is built as shown in Figure 1: Mould method, which comprises the following steps:
Step 1: establishing the hydrogovernor simulation analysis mathematical model based on Fractional Order PID speed-regulating system, and utilize The hybrid algorithm optimization objective function based on Fractional Order PID speed-regulating system;
Step 2: establishing the machinery hydraulic system and diversion system mathematical model of Turbine Governor System;
Step 3: hydraulic turbine generator model is established, by traditional single order hydraulic turbine generator model extension to 5 rank water wheels Machine model, and thus introduce exciter model and dynamic load model;
It is as shown in Figure 2: the step 1 specifically includes the following steps:
Step 1.1: Fractional Order PID speed-regulating system is obtained by classical PID speed-regulating system, method particularly includes:
The differential expressions of classical PID speed-regulating system are as follows:
U (t)=kpe(t)+kiDe(t)+kdDe(t); (1)
Wherein, kp indicates scale parameter, and setting range is 0.5~20;Ki indicates integral parameter, setting range 0.05s- 1~10s-1;Kd indicates differential parameter, and setting range is 0~5s;
Then by kp、ki、kdThree adjustment parameters replace with transient state slip coefficient b respectivelyt, damping time constant Td, accelerate Spend time constant Tn, it is specific:
The transmission function of classical PID speed-regulating system are as follows:
Wherein, bt is transient state slip coefficient, takes 0~1.0;Td is damping time constant, and setting range is 2s~20s;Tn For G-time constant, setting range is 0~2s;Pc and P is respectively that the power of the assembling unit is given and the power of the assembling unit is (for power tune Section mode);ycIt is respectively that guide vane opening gives and guide vane degree with y;Frequency modulation mode and aperture regulation mode;Fc and f difference It is given and unit (or power grid) frequency for frequency;Δ f is frequency difference;Δ f' is by frequency dead band efFrequency departure afterwards;
It is that Fractional Order PID adjusts the speed after Laplace is converted by the differential expressions (formula (1)) of classical PID speed-regulating system The differential expressions of system are as follows:
U (t)=kpe(t)+kiDαe(t)+kdDμe(t); (4)
Wherein α, μ > 0, α are integral order, and μ is differential order, and when the value of α and μ is all 1, then formula (4) becomes classical PID speed-regulating system (formula (1)), the introducing of fractional calculus make traditional PID adjust more application and flexibility;
The transmission function of Fractional Order PID speed-regulating system are as follows:
Step 1.2: the objective function of Fractional Order PID speed-regulating system is defined, method particularly includes:
In order to optimize the speed adjusting performance of controlled hydraulic turbine system, it is allocated as using double fault difference-product as the score described in step 1.1 The objective function of rank PID speed-regulating system:
Wherein, e (t) indicates the deviation of reality output and desired output, and t is time, ISEIndicate square deviation integral, value It is smaller, illustrate that the response speed of system is faster;ITAEIt is that time and square deviation integrate, value is smaller, and reflection system stability is got over It is good;
It is as shown in Figure 3: step 1.3: using described in hybrid algorithm Optimization Steps 1.2 based on Fractional Order PID speed-regulating system Objective function, method particularly includes:
By kp、ki、kd, five parameters of α, μ form initial drosophila positions, so that drosophila is carried out optimizing by flavor concentration, utilize BP The deep search characteristic of neural network is to kp、ki、kd, five parameter optimizations of α, μ solve, specifically includes the following steps:
Step 1.3.1: initialization drosophila population scale groupsize (400) and maximum number of iterations maxnum (400), Drosophila population position (X is randomly generatedaxis, Yaxis), iteration step value R is defined as (0.85-1);
Step 1.3.2: randomization drosophila population position and direction, according to drosophila population position and initial point distance DDistDetermine Drosophila population flavor concentration value Si
Step 1.3.3: the flavor concentration of drosophila individual is calculated, by drosophila population flavor concentration value SiFlavor concentration is substituted into sentence Determine function FfunctionIn, and flavor concentration optimum individual is found from each drosophila group;
Retain optimal flavor concentration bbestsmellWith corresponding (Xi, Yi) position, drosophila group SsmellThe coordinate is flown to, bbest indesIndicate optimal drosophila position;
Step 1.3.4: the drosophila position coordinates (X that step 1.3.3 is obtainedi, Yi) input BP neural network hidden layer nnet (2)(k), input and output are defined;
Wherein, i=1,2 ..., Q define n, Q value, Z according to controlled device complexityj (2)For drosophila population after optimization (Xi, Yi) position, wij (2)For hidden layer weighting coefficient, Zi(k)(2)Take activation functions Sigmoid function;
Step 1.3.5: BP neural network output layer n is definednet (3)(k), drosophila individual is updated, chooses E (k) as performance mistake Poor index, while retaining the position coordinates (X of the most dense drosophila individual of tastei, Yi);
Step 1.3.6: the performance indicator calculated by judging previous step exports optimal ginseng when reaching maximum number of iterations Number kp,ki,kd, α, μ, otherwise the process of end is transferred to step 1.3.2, until while reaching maximum number of iterations exports optimized parameter kp,ki,kd, α, μ terminate process.
The step 2 specifically includes the following steps:
Step 2.1: the mathematical model of the machinery hydraulic system of Turbine Governor System is established, following methods are specifically used:
The effect of machinery hydraulic system is transform electric and to zoom into the mechanical displacement letter with certain operating force Number;Due to second level servomotor responsive time constant Ty1Much smaller than main servomotor responsive time constant Ty, usually by it in modeling It is reduced to an inertial element, wherein y is machinery hydraulic system output signal, may also indicate that guide vane opening relative size, according to Water-Turbine Engine Adaption characteristic joined frequency dead band and saturation limiting element, main control valve and servomotor linear segment transmission function Are as follows:
Step 2.2: the mathematical model of the diversion system of Turbine Governor System is established, following methods are specifically used:
The mathematical model of the diversion system of Turbine Governor System is established, and according to the difference of diversion penstock length, is imitated Rigidity testing machine model generally is used in 700m or less in true modeling, when greater than 700m for precision the considerations of generally uses elastic water Hammer, water flow waterpower change procedure is described by following equation in diversion system:
The equation of motion:
Wherein H is head, and Q is water flow, and S is waterpipe cross-sectional area, and t is the time, and g is acceleration of gravity;
Flow equation:
Wherein a is water flow acceleration;
N=0 and n=1 is taken to obtain rigidity testing machine G through Taylor series expansion water attack equation in pressure diversion pipelineA1With Elastic water hammer GA2Transmission function:
Wherein TWFor rigid water attack time constant, Tr=2L/V is elastic water attack time constant, and L is pressure diversion conduit deferent Road length, V are water flow velocity of wave, H0For head, Q0For flow, V0For water flow base speed.
The step 3 specifically includes the following steps:
Hydraulic turbine generator model is established, by traditional single order hydraulic turbine generator model extension to 5 rank hydraulic turbine models, And thus introduce exciter model and dynamic load model;
Method particularly includes: traditional hydraulic turbine generator model is usually using first order modeling:
F is unit frequency;P is unit input power difference;Tn is that (dynamic frequency is special for unit (load) inertia time constant Property time constant), general Tn be 3s~12s;En is unit (load) static frequency self-regulation (characteristic) coefficient, with load characteristic Different and different, generally taking en is 0.5~2.0s;
Above-mentioned first order modeling can not reflect actual Parameters of Hydro-generator and parameters of excitation system, herein by The Synchronous Machine Models that MATLAB/Simulink emulation platform establishes five rank state-space equations are equivalent to hydrogenerator mould Type defines the basic parameter characteristic of the hydraulic turbine, it is contemplated that the dynamic characteristic in stator and rotor magnetic field and Damper Winding, model etc. Imitating circuit indicates in rotor reference system (dq frame).
Wherein, stator voltage equation are as follows:
Wherein UdFor stator d shaft voltage, UqFor stator q shaft voltage, φdAnd φqIndicate dp axis magnetic linkage, idAnd iqFor equivalent dq Electric current under coordinate system, raFor equivalent resistance, dq axis time transient internal voltage is E "dWith E "q, X "dWith X "qFor dq subtranient reactance;
Rotor f winding, dq winding voltage equation and equation of rotor motion are as follows:
Wherein d axis and q axis time constant (all as unit of s), d axis transient state is opened a way (Tdo') or short-circuit (Td') time Constant, d axis subtransient is opened a way (Tdo ") or short-circuit (Td ") time constant, and q axis transient state is opened a way (Tqo') or short-circuit (Tq') time Constant (only limits round rotor), and q axis subtransient is opened a way (Tqo ") or short-circuit (Tq ") time constant.E'qFor transient internal voltage, EeFor Dq winding electromotive force, X 'dFor d axis transient state reactance, XdAnd XqFor dq axis reactance, W is rotor machinery angular speed, TjFor generating set Inertia time constant, TmFor prime mover machine torque,To act on the uneven torque on armature spindle;
Excitation system is using IEEE typical DC exciter;When power system load increases suddenly, subtracts, to generator It carries out automatic field forcing, subtract magnetic, to improve the stability of electric system.
The comprehensive model modelling approach of a kind of hydraulic turbine based on Fractional Order PID speed-regulating system proposed by the present invention, according to Instance analysis is tested using MATLAB platform emulation, and Turbine Governor System Transfer Parameters (being shown in Table 1) is taken to be used as simulation example, Due to the rotary speed property of water turbine set, the dead zone of revolving speed ± 0.001 (pu) is provided with to speed-regulating system, wherein hydraulic turbine generator Partial parameters are shown in Table 2.
1 Turbine Governor System parameter of table
Wherein EqyIt is flow deviation relative value to the carry-over factor of servomotor deviation of stroke relative value, EqhFor flow deviation Carry-over factor of the relative value to head deviation relative value, EyIt is opposite to servomotor deviation of stroke for hydraulic turbine torque deviation relative value The carry-over factor of value, EhCarry-over factor of the hydraulic turbine torque deviation relative value to head deviation relative value.
2 hydraulic turbine generator parameter of table
Wherein PnFor rated power, VnFor rated line voltage, FnFor rated frequency, P is number of magnetic pole pairs, and H is inertia coeffeicent, d Axis synchronous reactance xd, transient reactance Xd' and subtransient reactance Xd ", q axis synchronous reactance Xq, transient reactance Xq'(are only when circle turns Son) and subtransient reactance Xq ", it is finally leakage reactance X1, WrefFor rated speed, VrefRated excitation voltage.
Compared from Tables 1 and 2: the present invention carries out Turbine Governor System on the basis of current research state It improves, it is as shown in Figure 4-Figure 6 to compare current mainstream speed-regulating system speed regulation result;Loaded portion according to hydraulic turbine operation characteristic according to With 40% and 75% two kind of load operating condition under investment, excision situation carry out simulation analysis, every simulation result of the hydraulic turbine is such as Shown in Fig. 7-Fig. 9.
Compared to classical Turbine Governor System, a kind of water wheels based on Fractional Order PID speed-regulating system proposed by the present invention The comprehensive model modelling approach of machine, based on Fractional Order PID speed-regulating system increasing due to controllable parameter, in regulating time and overshoot It measures two aspects and is superior to classical Turbine Governor System;Using the Adaptive PID Control of FOA algorithm and BP neural network algorithm Governing time, the overshoot of comparison PID and FOPID control significantly reduce, but above four kinds of control all there is system to shake number More defects;Traditional mode of speed regulation at least shakes 4 system above and just starts to tend towards stability, using mixing drosophila algorithm (BP- FOA) the Bi-objective FOPID optimal controller controlled further reduces regulating time and overshoot while drastically reducing shake Swing number, it is thus only necessary to which just steadily, specific speed regulation optimization is as shown in Figure 10 for 2 subsystems.
It is as shown in figure 11: by FOA algorithm, BP neural network algorithm and herein the BP-FOA algorithm that designs carry out from Optimizing Parameter analysis is adapted to, due to the high iteration step value of setting, the global optimizing ability of drosophila algorithm enhances, but the number of iterations More, optimal time is long;The optimizing ability of BP neural network algorithm is strong, fast convergence rate, but easily falls into local optimum, in face of negative There is wild effect when lotus fluctuation operation;BP-FOA algorithm utilizes BP neural network on the basis of global optimizing ability is strong Fast convergence further reduces adaptive adjustment time, improves 24.58%, shows in face of load fluctuation stronger suitable Ying Xingyu stability.
It is adjusted it can be seen that controlling FOPID adjuster by mixing drosophila algorithm Bi-objective in optimization by table 3, the comparison of table 4 There is obvious advantage in terms of fast time and control overshoot, is reduced on the basis of current mainstream BP neural network control algolithm 46.86% overshoot δ, shortens 36.59% regulating time Ts, wherein TrFor rise time, TpFor time to peak.
Different controller index simulation result comparisons under 3 load disturbance of table
Different controller optimization parameter comparisons under 4 load disturbance of table
It tests hydraulic turbine model initial strip 40% through BPFOA-FOPID biobjective scheduling to load, in 30s investment again 35% Load, 35% load of 60s excision.It is stable when in the switching load of 30s and 60s, the fluctuation of speed is in ± 2%, 3.9s, it mentions significantly High system stabilizing power;Due to the DC exciter of use, generator is used to encourage and subtract by force in load increase and decrease and is encouraged, substantially It is able to reflect exciter operating condition;Generator exports power output and the output basic held stationary of current conditions in load fluctuation, Embody good output performance.
This modeling analysis compared on the basis of establishing comprehensive hydraulic turbine model PID, FOPID, FOA algorithm with And BP neural network controls four kinds of mainstream single goal control methods, proposes hybrid algorithm BPFOA-FOPID two-objective programming control Mode processed tests it when in face of load fluctuation, to the quick-action and stability of Turbine Governor System control.Simulation result Show that hydraulic turbine regulating time and overshoot significantly reduce under BPFOA-FOPID Bi-objective control mode, compare conventional control Mode embodies stronger robustness and adaptability.
This patent utilizes MATLAB simulation modeling, is changed on the basis of forefathers' research to Turbine Governor System Into according to instance analysis, it is established that it is able to reflect Turbine Governor System and the comprehensive model of generator system parameters, For water turbine units capacity, frequency, rotor windings, stator winding, excitation winding, number of magnetic pole pairs, output voltage current conditions, Excitation voltage situation, load variations situation can be made a concrete analysis of, and all kinds of power station actual motions can be largely emulated When the variation of temporary stable state, compensate for the defect of traditional hydraulic turbine modeling.In practical applications to the operation feelings of practical water turbine set Condition, the variation of hydraulic turbine parameters has good reaction when in face of load fluctuation, mentions for accident forecast with system safety operation For ensureing.

Claims (4)

1. a kind of comprehensive model modelling approach of the hydraulic turbine of adaptive fractional rank PID speed-regulating system, which is characterized in that including with Lower step:
Step 1: establishing the hydrogovernor simulation analysis mathematical model based on Fractional Order PID speed-regulating system, and utilize mixing Objective function based on Fractional Order PID speed-regulating system described in algorithm optimization;
Step 2: establishing the machinery hydraulic system and diversion system mathematical model of Turbine Governor System;
Step 3: hydraulic turbine generator model is established, by traditional single order hydraulic turbine generator model extension to 5 rank hydraulic turbine moulds Type, and thus introduce exciter model and dynamic load model.
2. a kind of comprehensive model modeling side of the hydraulic turbine of adaptive fractional rank PID speed-regulating system according to claim 1 Method, which is characterized in that the step 1 specifically includes the following steps:
Step 1.1: Fractional Order PID speed-regulating system is obtained by classical PID speed-regulating system, method particularly includes:
The differential expressions of classical PID speed-regulating system are as follows:
U (t)=kpe(t)+kiDe(t)+kdDe(t);
Wherein, kp indicates scale parameter, and setting range is 0.5~20;Ki indicate integral parameter, setting range be 0.05s-1~ 10s-1, kd indicate differential parameter, and setting range is 0~5s;
By kp、ki、kdThree adjustment parameters replace with transient state slip coefficient b respectivelyt, damping time constant Td, the acceleration time it is normal Number Tn, it is specific:
The transmission function of classical PID speed-regulating system are as follows:
Wherein, bt is transient state slip coefficient, takes 0~1.0;Td is damping time constant, and setting range is 2s~20s, and Tn is to add Velocity Time constant, setting range are that 0~2s, pc and P are respectively that the power of the assembling unit is given and the power of the assembling unit, ycIt is respectively to lead with y Leaf aperture gives and guide vane degree, and fc and f are respectively that frequency is given and unit frequency, Δ f are frequency difference, and Δ f' is by frequency dead band efFrequency departure afterwards;
The differential expressions of above-mentioned classical PID speed-regulating system are transformed to the differential of Fractional Order PID speed-regulating system by Laplace Expression formula:
U (t)=kpe(t)+kiDαe(t)+kdDμe(t); (4)
Wherein α, μ > 0, α are integral order, and μ is differential order;
The transmission function of Fractional Order PID speed-regulating system are as follows:
Step 1.2: the objective function of Fractional Order PID speed-regulating system is defined, method particularly includes:
It is allocated as using double fault difference-product as the objective function of the Fractional Order PID speed-regulating system described in step 1.1:
Wherein, e (t) indicates the deviation of reality output and desired output, and t is time, ISEIndicate square deviation integral;
Step 1.3: using the objective function based on Fractional Order PID speed-regulating system described in hybrid algorithm Optimization Steps 1.2, specifically Method are as follows:
By kp、ki、kd, five parameters of α, μ form initial drosophila positions, so that drosophila is carried out optimizing by flavor concentration, utilize BP nerve The deep search characteristic of network is to kp、ki、kd, five parameter optimizations of α, μ solve, specifically includes the following steps:
Step 1.3.1: initialization drosophila population scale groupsize (400) and maximum number of iterations maxnum (400), at random Generate drosophila population position (Xaxis, Yaxis), iteration step value R is defined as (0.85-1);
Step 1.3.2: randomization drosophila population position and direction, according to drosophila population position and initial point distance DDistDetermine drosophila Population flavor concentration value Si
Step 1.3.3: the flavor concentration of drosophila individual is calculated, by drosophila population flavor concentration value SiIt substitutes into flavor concentration and determines letter Number FfunctionIn, and flavor concentration optimum individual is found from each drosophila group;
Retain optimal flavor concentration bbestsmellWith corresponding (Xi, Yi) position, drosophila group SsmellThe coordinate is flown to, bbestindesIndicate optimal drosophila position;
Step 1.3.4: the drosophila position coordinates (X that step 1.3.3 is obtainedi, Yi) input BP neural network hidden layer nnet (2) (k), input and output are defined;
Wherein, i=1,2 ..., Q define n, Q value, Z according to controlled device complexityj (2)For drosophila population (X after optimizationi, Yi) Position, wij (2)For hidden layer weighting coefficient, Zi(k)(2)Take activation functions Sigmoid function;
Step 1.3.5: BP neural network output layer n is definednet (3)(k), drosophila individual is updated, E (k) is chosen and refers to for performance error Mark, while retaining the position coordinates (X of the most dense drosophila individual of tastei, Yi);
Step 1.3.6: the performance indicator calculated by judging previous step exports optimized parameter k when reaching maximum number of iterationsp, ki,kd, α, μ, otherwise the process of end is transferred to step 1.3.2, until exporting optimized parameter k when reaching maximum number of iterationsp,ki, kd, α, μ terminate process.
3. a kind of comprehensive model modeling side of the hydraulic turbine of adaptive fractional rank PID speed-regulating system according to claim 1 Method, which is characterized in that the step 2 specifically includes the following steps:
Step 2.1: the mathematical model of the machinery hydraulic system of Turbine Governor System is established, following methods are specifically used:
The effect of machinery hydraulic system is transform electric and to zoom into the mechanical displacement signal with certain operating force;By In second level servomotor responsive time constant Ty1Much smaller than main servomotor responsive time constant Ty, usually it is reduced in modeling One inertial element, wherein y is machinery hydraulic system output signal, may also indicate that guide vane opening relative size, according to the hydraulic turbine Control characteristic joined frequency dead band and saturation limiting element, main control valve and servomotor linear segment transmission function are as follows:
Step 2.2: the mathematical model of the diversion system of Turbine Governor System is established, following methods are specifically used:
According to the difference of diversion penstock length, on simulation modeling, rigidity testing machine model is used in 700m or less, is greater than 700m Shi Caiyong elastic water hammer, water flow waterpower change procedure is described by following equation in diversion system:
The equation of motion:
Wherein H is head, and Q is water flow, and S is waterpipe cross-sectional area, and t is the time, and g is acceleration of gravity;
Flow equation:
Wherein a is water flow acceleration;
N=0 and n=1 is taken to obtain rigidity testing machine G through Taylor series expansion water attack equation in pressure diversion pipelineA1And elasticity Water hammer GA2Transmission function:
Wherein TWFor rigid water attack time constant, Tr=2L/V is elastic water attack time constant, and L is pressure diversion conduit deferent road length Degree, V are water flow velocity of wave, H0For head, Q0For flow, V0For water flow base speed.
4. a kind of comprehensive model modeling side of the hydraulic turbine of adaptive fractional rank PID speed-regulating system according to claim 1 Method, which is characterized in that the step 3 specifically includes the following steps:
Hydraulic turbine generator model is established, by traditional single order hydraulic turbine generator model extension to 5 rank hydraulic turbine models, and by This introduces exciter model and dynamic load model;Method particularly includes:
Step 3.1: the first order modeling that traditional hydraulic turbine generator model uses are as follows:
Wherein, f is unit frequency, and p is unit input power difference, and Tn is unit set inertia time constant, general TnFor 3s~12s, enFor unit static frequency self-regulation coefficient, enTake 0.5~2.0s;
Above-mentioned first order modeling can not reflect actual Parameters of Hydro-generator and parameters of excitation system, herein by The Synchronous Machine Models that MATLAB/Simulink emulation platform establishes five rank state-space equations are equivalent to hydrogenerator mould Type defines the basic parameter characteristic of the hydraulic turbine, it is contemplated that the dynamic characteristic in stator and rotor magnetic field and Damper Winding, model etc. Imitating circuit indicates in rotor reference system (dq frame).
Wherein, stator voltage equation are as follows:
Wherein UdFor stator d shaft voltage, UqFor stator q shaft voltage, φdAnd φqIndicate dp axis magnetic linkage, idAnd iqFor equivalent dq coordinate It is lower electric current, raFor equivalent resistance, dq axis time transient internal voltage is E "dWith E "q, X "dAnd X "qFor dq subtranient reactance;
Rotor f winding, dq winding voltage equation and equation of rotor motion are as follows:
Wherein d axis and q axis time constant (all as unit of s), d axis transient state is opened a way (Tdo') or short-circuit (Td') time constant, D axis subtransient is opened a way (Tdo ") or short-circuit (Td ") time constant, and q axis transient state is opened a way (Tqo') or short-circuit (Tq') time constant (only limiting round rotor), q axis subtransient is opened a way (Tqo ") or short-circuit (Tq ") time constant.E'qFor transient internal voltage, EeFor dq around Group electromotive force, X'dFor d axis transient state reactance, XdAnd XqFor dq axis reactance, W is rotor machinery angular speed, TjFor the used of generating set Property time constant, TmFor prime mover machine torque,To act on the uneven torque on armature spindle.
CN201910643191.6A 2019-07-16 2019-07-16 Water turbine comprehensive model modeling method based on fractional PID speed regulation system Active CN110262223B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910643191.6A CN110262223B (en) 2019-07-16 2019-07-16 Water turbine comprehensive model modeling method based on fractional PID speed regulation system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910643191.6A CN110262223B (en) 2019-07-16 2019-07-16 Water turbine comprehensive model modeling method based on fractional PID speed regulation system

Publications (2)

Publication Number Publication Date
CN110262223A true CN110262223A (en) 2019-09-20
CN110262223B CN110262223B (en) 2022-10-18

Family

ID=67926550

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910643191.6A Active CN110262223B (en) 2019-07-16 2019-07-16 Water turbine comprehensive model modeling method based on fractional PID speed regulation system

Country Status (1)

Country Link
CN (1) CN110262223B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110989345A (en) * 2019-11-28 2020-04-10 国网福建省电力有限公司 Water turbine system control parameter optimization method based on biophysical optimization algorithm
CN111259864A (en) * 2020-03-04 2020-06-09 哈尔滨理工大学 Method for identifying running state of water turbine
CN111736466A (en) * 2020-06-08 2020-10-02 武汉理工大学 Optimal control method and system for quick load shedding system of semi-submersible platform
CN112000017A (en) * 2020-09-08 2020-11-27 金陵科技学院 Global stabilization control method of fractional order water turbine adjusting system
CN112733424A (en) * 2020-12-12 2021-04-30 国网新源控股有限公司回龙分公司 Modeling simulation method and system for pumped storage power station generator
CN112966394A (en) * 2021-03-31 2021-06-15 华中科技大学 Simulation method and system for dynamic characteristics of hydroelectric generator group under hydraulic coupling condition
CN117650583A (en) * 2024-01-30 2024-03-05 三峡金沙江云川水电开发有限公司 Hydropower station one-pipe multi-machine grid-connection multi-target coordination optimization control method and system

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010229962A (en) * 2009-03-30 2010-10-14 Mitsubishi Electric Corp Hydraulic turbine and speed governing controller for pump turbine
CN102175445A (en) * 2011-02-24 2011-09-07 华中科技大学 Simulation test device for hydroturbine speed-regulating system
CN102352812A (en) * 2011-07-18 2012-02-15 华北电力大学 Sliding mode-based hydro turbine governing system dead zone nonlinear compensation method
CN103590969A (en) * 2013-11-20 2014-02-19 华中科技大学 PID hydraulic turbine governor parameter optimization method based on multi-working-condition time domain response
CN103942625A (en) * 2014-04-23 2014-07-23 天津大学 Hydropower station plant dam structure vibration response prediction method based on FOA-GRNN
CN106125552A (en) * 2016-08-08 2016-11-16 国家电网公司 Pump-storage generator governing system fuzzy score rank PID control method
CN106849814A (en) * 2017-03-27 2017-06-13 无锡开放大学 Leapfroged Fuzzy Neural PID linear synchronous generator control method based on fruit bat
CN107834610A (en) * 2017-11-29 2018-03-23 西南交通大学 A kind of mains frequency dynamic analysing method for considering hydraulic turbine water hammer effect
CN108564235A (en) * 2018-07-13 2018-09-21 中南民族大学 A kind of improved FOA-BPNN exit times prediction technique
CN109634116A (en) * 2018-09-04 2019-04-16 贵州大学 A kind of acceleration adaptive stabilizing method of fractional order mechanical centrifugal governor system

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010229962A (en) * 2009-03-30 2010-10-14 Mitsubishi Electric Corp Hydraulic turbine and speed governing controller for pump turbine
CN102175445A (en) * 2011-02-24 2011-09-07 华中科技大学 Simulation test device for hydroturbine speed-regulating system
CN102352812A (en) * 2011-07-18 2012-02-15 华北电力大学 Sliding mode-based hydro turbine governing system dead zone nonlinear compensation method
CN103590969A (en) * 2013-11-20 2014-02-19 华中科技大学 PID hydraulic turbine governor parameter optimization method based on multi-working-condition time domain response
CN103942625A (en) * 2014-04-23 2014-07-23 天津大学 Hydropower station plant dam structure vibration response prediction method based on FOA-GRNN
CN106125552A (en) * 2016-08-08 2016-11-16 国家电网公司 Pump-storage generator governing system fuzzy score rank PID control method
CN106849814A (en) * 2017-03-27 2017-06-13 无锡开放大学 Leapfroged Fuzzy Neural PID linear synchronous generator control method based on fruit bat
CN107834610A (en) * 2017-11-29 2018-03-23 西南交通大学 A kind of mains frequency dynamic analysing method for considering hydraulic turbine water hammer effect
CN108564235A (en) * 2018-07-13 2018-09-21 中南民族大学 A kind of improved FOA-BPNN exit times prediction technique
CN109634116A (en) * 2018-09-04 2019-04-16 贵州大学 A kind of acceleration adaptive stabilizing method of fractional order mechanical centrifugal governor system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
蔡超豪等: "具有区域极点偏置的水轮机调速系统的H2/H∞混合控制", 《电力科学与工程》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110989345A (en) * 2019-11-28 2020-04-10 国网福建省电力有限公司 Water turbine system control parameter optimization method based on biophysical optimization algorithm
CN111259864A (en) * 2020-03-04 2020-06-09 哈尔滨理工大学 Method for identifying running state of water turbine
CN111736466A (en) * 2020-06-08 2020-10-02 武汉理工大学 Optimal control method and system for quick load shedding system of semi-submersible platform
CN111736466B (en) * 2020-06-08 2021-09-10 武汉理工大学 Optimal control method and system for quick load shedding system of semi-submersible platform
CN112000017A (en) * 2020-09-08 2020-11-27 金陵科技学院 Global stabilization control method of fractional order water turbine adjusting system
CN112733424A (en) * 2020-12-12 2021-04-30 国网新源控股有限公司回龙分公司 Modeling simulation method and system for pumped storage power station generator
CN112966394A (en) * 2021-03-31 2021-06-15 华中科技大学 Simulation method and system for dynamic characteristics of hydroelectric generator group under hydraulic coupling condition
CN112966394B (en) * 2021-03-31 2024-04-23 华中科技大学 Simulation method and system for dynamic characteristics of hydroelectric generating set under hydraulic coupling condition
CN117650583A (en) * 2024-01-30 2024-03-05 三峡金沙江云川水电开发有限公司 Hydropower station one-pipe multi-machine grid-connection multi-target coordination optimization control method and system
CN117650583B (en) * 2024-01-30 2024-04-26 三峡金沙江云川水电开发有限公司 Hydropower station one-pipe multi-machine grid-connection multi-target coordination optimization control method and system

Also Published As

Publication number Publication date
CN110262223B (en) 2022-10-18

Similar Documents

Publication Publication Date Title
CN110262223A (en) A kind of comprehensive model modelling approach of hydraulic turbine based on Fractional Order PID speed-regulating system
CN109217362B (en) System and method for positioning low-frequency oscillation disturbance source of grid-connected system of double-fed fan
Ghosh et al. Doubly fed induction generator (DFIG)-based wind farm control framework for primary frequency and inertial response application
CN109308390B (en) Combined regulation simulation system and method for wind/light energy and hydroelectric generating set of power grid at transmitting and receiving ends
CN105226723A (en) A kind of double-fed fan motor unit is based on the virtual inertia control method of wind power tracking Automatic adjusument
CN106910142A (en) A kind of power system frequency characteristic computing method containing the active frequency coupling of wind-powered electricity generation
CN105863948B (en) A kind of band, which becomes, rises tailwater tunnel hydrogovernor variable parameter control method
CN111525594A (en) Control method and device for speed regulating system of hydroelectric generating set
CN105591402B (en) A kind of modeling and simulation method and apparatus of direct-drive permanent-magnetism Wind turbines
CN106227950B (en) Wind turbines primary frequency control system modeling method based on pitch control
Yin et al. Integrated pitch control for wind turbine based on a novel pitch control system
CN110195660A (en) Aero-engine control device based on depth Q study
CN107240918A (en) A kind of power system Equivalent Simplification method of wind power integration
CN112542855B (en) Modeling and simulation method for phasor model of double-fed wind power generation system
CN111384730B (en) Method for determining control parameters of virtual inertia of fan
CN111478365B (en) Optimization method and system for control parameters of virtual synchronizer of direct-drive wind turbine generator
CN104201954A (en) Marine electric power system exciting voltage regulation method
CN107743001B (en) Load simulation method, frequency converter, load simulator and static frequency conversion starting system
Lv et al. Characteristic analysis and optimal regulation of primary frequency regulation condition in low water head area based on hydraulic-mechanical-electrical coupling model of pumped storage unit
Neshati et al. Hardware-in-the-loop testing of wind turbine nacelles for electrical certification on a dynamometer test rig
Moon et al. Modified PID load-frequency control with the consideration of valve position limits
CN104196639B (en) Gas turbine control method and device
CN109711013A (en) A kind of Adaptive System of Water-Turbine Engine electrohydraulic actuator model
CN108549207A (en) A kind of method of Adaptive System of Water-Turbine Engine control parameter
Sahin et al. Performance comparison of two turbine blade pitch controller design methods based on equilibrium and frozen wake assumptions

Legal Events

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