CN112257244B - Modeling method of micro gas turbine power generation system - Google Patents
Modeling method of micro gas turbine power generation system Download PDFInfo
- Publication number
- CN112257244B CN112257244B CN202011105307.XA CN202011105307A CN112257244B CN 112257244 B CN112257244 B CN 112257244B CN 202011105307 A CN202011105307 A CN 202011105307A CN 112257244 B CN112257244 B CN 112257244B
- Authority
- CN
- China
- Prior art keywords
- calculating
- follows
- calculation formula
- formula
- specific calculation
- 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.)
- Active
Links
- 238000010248 power generation Methods 0.000 title claims abstract description 14
- 238000000034 method Methods 0.000 title claims abstract description 12
- 239000000446 fuel Substances 0.000 claims abstract description 30
- 230000001133 acceleration Effects 0.000 claims abstract description 13
- 238000004088 simulation Methods 0.000 claims abstract description 4
- 230000001419 dependent effect Effects 0.000 claims description 17
- 230000004907 flux Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims 23
- 239000007789 gas Substances 0.000 abstract description 16
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract description 14
- 239000003345 natural gas Substances 0.000 abstract description 7
- 238000011161 development Methods 0.000 abstract description 5
- 238000006243 chemical reaction Methods 0.000 abstract description 4
- 230000007613 environmental effect Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 230000002860 competitive effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000013468 resource allocation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Business, Economics & Management (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Public Health (AREA)
- Geometry (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- General Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- Human Resources & Organizations (AREA)
- Power Engineering (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Algebra (AREA)
- Evolutionary Computation (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Control Of Electric Motors In General (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
A modeling method for a micro gas turbine power generation system, comprising the steps of: s1, setting model parameters; s2, calculating a speed-related fuel input signal; s3, calculating an acceleration-related fuel input signal; s4, calculating a temperature-related fuel input signal; s5, calculating a fuel control signal; s6, calculating the number of forward points; s7, calculating and controlling delay output; s8, calculating a fuel flow signal; s9, calculating forward points; s10, calculating flow delay output; s11, calculating the temperature; s12, calculating forward points; s13, calculating flow delay output; s14, calculating torque input; s15, calculating torque; s16, calculating the electromagnetic torque of the motor; s17, calculating the acceleration of the motor and the like; and finally, judging whether the simulation is finished or not, if so, exiting, and otherwise, enabling k = k +1 to enter S2 for continuous iteration. The invention can realize the conversion of natural gas energy to electric energy, promote the development and utilization of clean energy, optimize energy distribution and scheduling and improve the energy utilization rate.
Description
Technical Field
The invention relates to the technical field of comprehensive energy system planning, in particular to a modeling method of a micro gas turbine power generation system.
Background
With the development of technologies such as power electronics, new energy, automatic control, energy storage and the like, a comprehensive energy system is not popular, has the advantages of cleanness, environmental protection, low carbon, flexibility and the like, greatly optimizes resource allocation, and becomes an important solution for dealing with energy crisis and environmental problems.
The comprehensive energy system integrates various access modes such as electric power, natural gas, wind power, photovoltaic and the like, and provides resources such as electricity, heat, cold and the like for users. Various forms of coupling and conversion exist among different energy sources, energy distribution and scheduling can be optimized, the energy utilization rate is improved, related experts and scholars put forward concepts of an energy concentrator and an energy router, and wide attention is drawn. Compared with different types of energy, the natural gas has the advantages of obvious advantages, high efficiency, cleanness and environmental protection. With the continuous development of natural gas systems, research for combining the natural gas systems with power systems is more and more common. In the energy coupling unit, the gas turbine set has short construction time, high power generation efficiency and high regulation speed, and is widely applied. The micro gas turbine, one of the most commercially competitive distributed power sources at present, has the advantages of abundant available fuels, small volume, light weight, little pollution, simple operation and maintenance, easy realization of combined cooling, heating and power, and the like, and thus has attracted increasing attention. Therefore, in a large environment that encourages high penetration of distributed power, it is of great significance to study modeling of micro gas turbine power generation systems.
In the literature, a chopper is adopted to simulate the electrical characteristics of a micro gas turbine power generation system, so that the electrical characteristics of the system under different working conditions and operation modes can be reproduced within a certain precision range, and thermal power units such as air inlet, heat return, combustion, exhaust and the like of the micro gas turbine power generation system can be omitted; two methods for researching the characteristics of the power generation system of the micro gas turbine are given in the literature, but the methods do not model the micro gas turbine; in some documents, a modeling control method for each part of a micro gas turbine is studied, and a thermal unit of the micro gas turbine is analyzed in detail, but the above method is not suitable for analyzing the electrical characteristics of a power generation system of the micro gas turbine.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provide a modeling method of a micro gas turbine power generation system, so that the micro gas turbine power generation system can realize the conversion of natural gas energy to electric power energy, promote the development and utilization of clean energy, optimize energy distribution and scheduling and improve the energy utilization rate.
The invention adopts the following technical scheme for solving the technical problems:
a modeling method of a micro gas turbine power generation system, characterized by comprising the following S:
s1, setting model parameters, including: velocity correlation coefficient { a } ω0 、b ω0 、b ω1 }, reference rotation speed omega r Calculating the period T S Temperature dependent coefficient { b } T0 、b T1 Control coefficient kc and control delay T d1 Traffic delay { T } d2 、T d3 H, fuel related coefficient a f0 、a f1 、b f0 Torque correlation coefficient { a } A0 、b A0 J moment of inertia, N number of pole pairs p Stator resistance R, inductance L, rotor flux linkage psi f Etc.;
s2, calculating a speed-dependent fuel input signal f ω,k As shown in formula
f ω,k ＝a ω0 ·f ω,k-1 +b ω0 ·ω k-1 +b ω1 ·ω k-2 -(b ω0 +b ω1 )·ω r (1)
In the formula, ω k-1 、ω k-2 The rotating speed at the previous moment and the previous two moments; subscripts k, k-1, k-2 respectively represent the current time, the previous two times, and so on;
s3, calculating an acceleration-dependent fuel input signal f α,k As shown in the following formula
f α,k ＝f α,k-1 +T s ·(α k-1 -0.01) (2)
In the formula, alpha k-1 Acceleration at the previous moment;
s4, calculating a temperature-dependent fuel input signal f T,k As shown in the following formula
f T,k ＝f T,k-1 +b T0 ·T k-1 +b T1 ·T k-2 -950·(b T0 +b T1 ) (3)
In the formula, T k-1 、T k-2 The temperatures at the previous moment and the previous two moments respectively;
s5, calculating a fuel control signal f c,k As shown in the following formula
f c,k ＝k c ·ω k-1 ·min(f ω,k ,f α,k ,f T,k ) (4)
In the formula, min represents the minimum value;
s6, calculating the number m of forward points 1 As shown in the following formula
m 1 ＝INT(T d1 /T s ) (5)
In the formula, the symbol INT represents an integer part;
s7, calculating and controlling the delay output f d,k As shown in the following formula
In which the subscript k-m 1 、k-1-m 1 Respectively represents m 1 、m 1 Before +1 time;
s8, calculating a fuel flow signal W f,k As shown in the following formula
W f,k ＝a f0 ·W f,k-1 +a f1 ·W f,k-2 +b f0 ·(f d,k +0.23) (7)
S9, calculating the number m of forward push points 2 As shown in the following formula
m 2 ＝INT(T d2 /T s ) (8)
S10, calculating flow delay output W d1,k As shown in the following formula
S11, calculating the temperature T k As shown in the following formula
T k ＝950-700(1-W f1,k )+550(1-ω k-1 )
(10)
S12, calculating the number m of forward push points 3 As shown in the following formula
m 3 ＝INT(T d3 /T s )
(11)
S13, calculating flow delay output W d2,k As shown in the following formula
S14, calculating a torque input A k As shown in the following formula
A k ＝a A0 ·A k-1 +b A0 ·W d2,k
(13)
S15, calculating the torque T m,k As shown in the following formula
T m,k ＝1.3·(A k -0.23)+0.5·(1-ω k-1 )
(14)
S16, calculating the electromagnetic torque T of the motor e,k As shown in the following formula
T e,k ＝P L,k /ω k-1
(15)
In the formula, P L,k The load active power at the current moment;
s17, calculating the acceleration alpha of the motor k As shown in the following formula
α k ＝(T m,k -T e,k )/J
(16)
S18, calculating the rotating speed omega of the motor k As shown in the following formula
ω k ＝ω k-1 +T s ·α k
(17)
S19, calculating the electrical angle theta of the motor k As shown in the following formula
θ k ＝θ k-1 +N p ·ω k
(18)
S20, calculating the quadrature-direct axis voltage u q,k 、u d,k As shown in the following formula
u d,k ＝R·i d,k-1 +L·(i d,k-1 -i d,k-2 )/T s -N p ·ω k-1 ·L·i q,k-1
(19)
u q,k ＝R·i q,k-1 +L·(i q,k-1 -i q,k-2 )/T s +N p ·ω k-1 ·ψ f
(20)
S21, calculating quadrature axis current i q,k As shown in the following formula
i q,k ＝T e,k /(N p ·ψ f )
(21)
S22, calculating the direct axis current i d,k As shown in the following formula
i d,k ＝(P L,k -u q,k ·i q,k )/u d,k
(22)
S23, calculating the three-phase voltage { u } of the motor a,k 、u b,k 、u c,k }, three-phase currents { i a,k 、i b,k 、i c,k }, as shown in the following formula
And S24, judging whether the simulation is finished, if so, exiting, and otherwise, enabling k = k +1 to enter S2 for continuous iteration.
The invention has the beneficial effects that: by adopting the method, the conversion of natural gas energy to electric power energy can be realized, the development and utilization of clean energy are promoted, the energy distribution and dispatching are optimized, and the energy utilization rate is improved.
Drawings
FIG. 1 is a flow chart of the model algorithm of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the specification, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, belong to the scope of the present invention.
As shown in fig. 1, the specific implementation S of the present invention is:
s1, setting model parameters, including: velocity correlation coefficient { a } ω0 、b ω0 、b ω1 }, reference rotation speed omega r Calculating the period T S Temperature dependent coefficient { b } T0 、b T1 }, fuel control coefficient kc and control delay T d1 Traffic delay { T } d2 、T d3 Fuel related coefficient { a } f0 、a f1 、b f0 Torque correlation coefficient { a } A0 、b A0 J moment of inertia, N number of pole pairs p Stator resistance R, inductance L, rotor flux linkage psi f Etc.;
s2, calculating a speed-dependent fuel input signal f ω,k As shown in the formula
f ω,k ＝a ω0 ·f ω,k-1 +b ω0 ·ω k-1 +b ω1 ·ω k-2 -(b ω0 +b ω1 )·ω r (1)
In the formula, ω k-1 、ω k-2 The rotating speed at the previous moment and the previous two moments; subscripts k, k-1, k-2 respectively represent the current time, the previous two times, and so on;
s3, calculating an acceleration-dependent fuel input signal f α,k As shown in the following formula
f α,k ＝f α,k-1 +T s ·(α k-1 -0.01) (2)
In the formula, alpha k-1 Acceleration at the previous moment;
s4, calculating a temperature-dependent fuel input signal f T,k As shown in the following formula
f T,k ＝f T,k-1 +b T0 ·T k-1 +b T1 ·T k-2 -950·(b T0 +b T1 ) (3)
In the formula, T k-1 、T k-2 The temperatures at the previous moment and the previous two moments respectively;
s5, calculating a fuel control signal f c,k As shown in the following formula
f c,k ＝k c ·ω k-1 ·min(f ω,k ,f α,k ,f T,k ) (4)
In the formula, min represents the minimum value;
s6, calculating the number m of forward points 1 As shown in the following formula
m 1 ＝INT(T d1 /T s ) (5)
In the formula, the symbol INT represents an integer part;
s7, calculating and controlling a delay output f d,k As shown in the following formula
In which the subscript k-m 1 、k-1-m 1 Respectively represents m 1 、m 1 Before +1 time;
s8, calculating a fuel flow signal W f,k As shown in the following formula
W f,k ＝a f0 ·W f,k-1 +a f1 ·W f,k-2 +b f0 ·(f d,k +0.23) (7)
S9, calculating the number m of forward push points 2 As shown in the following formula
m 2 ＝INT(T d2 /T s ) (8)
S10, calculating flow delay output W d1,k As shown in the following formula
S11, calculating the temperature T k As shown in the following formula
T k ＝950-700(1-W f1,k )+550(1-ω k-1 )
(10)
S12, calculating the number m of forward points 3 As shown in the following formula
m 3 ＝INT(T d3 /T s )
(11)
S13, calculating flow delay output W d2,k As shown in the following formula
S14, calculating a torque input A k As shown in the following formula
A k ＝a A0 ·A k-1 +b A0 ·W d2,k
(13)
S15, calculating the torque T m,k As shown in the following formula
T m,k ＝1.3·(A k -0.23)+0.5·(1-ω k-1 )
(14)
S16, calculating the electromagnetic torque T of the motor e,k As shown in the following formula
T e,k ＝P L,k /ω k-1
(15)
In the formula, P L,k The load active power at the current moment;
s17, calculating the acceleration alpha of the motor k As shown in the following formula
α k ＝(T m,k -T e,k )/J
(16)
S18, calculating the rotating speed omega of the motor k As shown in the following formula
ω k ＝ω k-1 +T s ·α k
(17)
S19, calculating the electrical angle theta of the motor k As shown in the following formula
θ k ＝θ k-1 +N p ·ω k
(18)
S20, calculating the voltage u of the quadrature-direct axis q,k 、u d,k As shown in the following formula
u d,k ＝R·i d,k-1 +L·(i d,k-1 -i d,k-2 )/T s -N p ·ω k-1 ·L·i q,k-1
(19)
u q,k ＝R·i q,k-1 +L·(i q,k-1 -i q,k-2 )/T s +N p ·ω k-1 ·ψ f
(20)
S21, calculating quadrature axis current i q,k As shown in the following formula
i q,k ＝T e,k /(N p ·ψ f )
(21)
S22, calculating the direct axis current i d,k As shown in the following formula
i d,k ＝(P L,k -u q,k ·i q,k )/u d,k
(22)
S23, calculating the three-phase voltage { u } of the motor a,k 、u b,k 、u c,k Triple phase current i a,k 、i b,k 、i c,k As shown in the following formula
And S24, judging whether the simulation is finished, if so, exiting, and otherwise, enabling k = k +1 to enter S2 for continuous iteration.
Although the present invention has been described in detail with reference to the foregoing embodiments, it will be apparent to those skilled in the art that various changes and modifications can be made therein without departing from the spirit and scope of the invention.
Claims (1)
1. A method of modeling a micro gas turbine power generation system, comprising:
s1, setting model parameters;
s2, calculating a speed-dependent fuel input signal f ω,k ；
S3, calculating an acceleration-dependent fuel input signal f α,k ；
S4, calculating a temperature-dependent fuel input signal f T,k ；
S5, calculating a fuel control signal f c,k ；
S6, calculating the number m of forward points 1 ；
S7, calculating and controlling the delay output f d,k ；
S8, calculating a fuel flow signal W f,k ；
S9, calculating the number m of forward points 2 ；
S10, calculating flow delay output W d1,k ；
S11, calculating the temperature T k ；
S12, calculating the number m of forward points 3 ；
S13, calculating flow delay output W d2,k ；
S14, calculating a torque input A k ；
S15, calculating the torque T m,k ；
S16, calculating the electromagnetic torque T of the motor e,k ；
S17, calculating the acceleration alpha of the motor k ；
S18, calculating the rotating speed omega of the motor k ；
S19, calculating the electric angle theta of the motor k ；
S20, calculating the voltage u of the quadrature-direct axis q,k 、u d,k ；
S21, calculating quadrature axis current i q,k ；
S22, calculating the direct axis current i d,k ；
S23, calculating the three-phase voltage { u } of the motor a,k 、u b,k 、u c,k Triple phase current i a,k 、i b,k 、i c,k }；
S24, judging whether the simulation is finished or not, if so, exiting, and otherwise, enabling k = k +1 to enter S2 for continuous iteration;
speed dependent fuel input signal f of S2 ω,k The specific calculation formula is as follows,
f ω,k ＝a ω0 ·f ω,k-1 +b ω0 ·ω k-1 +b ω1 ·ω k-2 -(b ω0 +b ω1 )·ω r
in the formula, a ω0 、b ω0 、b ω1 As a velocity-related coefficient, ω r Is the reference rotational speed, omega k-1 、ω k-2 The rotating speed at the previous moment and the previous two moments; subscripts k, k-1, k-2 respectively represent the current time, the previous time, and the previous two times;
acceleration-dependent fuel input signal f of S3 α,k The specific calculation formula is as follows,
f α,k ＝f α,k-1 +T s ·(α k-1 -0.01)
in the formula, alpha k-1 Acceleration at the previous moment, T S Is a calculation cycle;
s4 temperature dependent Fuel input Signal f T,k The specific calculation formula is as follows,
f T,k ＝f T,k-1 +b T0 ·T k-1 +b T1 ·T k-2 -950·(b T0 +b T1 )
in the formula, b T0 、b T1 Is a temperature dependent coefficient, T k-1 、T k-2 The temperatures at the previous moment and the previous two moments respectively;
fuel control signal f of S5 c,k The specific calculation formula is as follows,
f c,k ＝k c ·ω k-1 ·min(f ω,k ,f α,k ,f T,k )
in the formula, min represents the minimum value;
number m of forward points of S6 1 The specific calculation formula is as follows,
m 1 ＝INT(T d1 /T s )
in the formula, T d1 For controlling the delay, the symbol INT represents the integer part;
s7 control delay output f d,k The specific calculation formula is as follows,
in which the subscript k-m 1 、k-1-m 1 Respectively represents m 1 、m 1 Before +1 time;
s8 Fuel flow Signal W f,k The specific calculation formula is as follows,
W f,k ＝a f0 ·W f,k-1 +a f1 ·W f,k-2 +b f0 ·(f d,k +0.23)
in the formula, a f0 、a f1 、b f0 Is a fuel correlation coefficient;
number m of pushforward points of S9 2 The specific calculation formula is as follows,
m 2 ＝INT(T d2 /T s )
in the formula, T d2 A flow delay that is temperature dependent;
s10 flow delay output W d1,k The specific calculation formula is as follows,
temperature T of S11 k The specific calculation formula is as follows,
T k ＝950-700(1-W d1,k )+550(1-ω k-1 )；
number m of forward points of S12 3 The specific calculation formula is as follows,
m 3 ＝INT(T d3 /T s )
in the formula, T d3 A torque dependent flow delay;
s13 flow delay output W d2,k In particular toThe calculation formula is as follows,
torque input A of S14 k The specific calculation formula is as follows,
A k ＝a A0 ·A k-1 +b A0 ·W d2,k
in the formula, a A0 、b A0 Is a torque correlation coefficient;
torque T of S15 m,k The specific calculation formula is as follows,
T m,k ＝1.3·(A k -0.23)+0.5/(1-ω k-1 )；
electromagnetic torque T of S16 e,k The specific calculation formula is as follows,
T e,k ＝P L,k /ω k-1
in the formula, P L,k The load active power at the current moment;
s17 motor acceleration α k The specific calculation formula is as follows,
α k ＝(T m,k -T e,k )/J
in the formula, T m,k Is torque, T e,k Is the electromagnetic torque, and J is the moment of inertia;
s18 motor rotation speed ω k The specific calculation formula is as follows,
ω k ＝ω k-1 +T s ·α k
s19 electric angle theta of motor k The specific calculation formula is as follows,
θ k ＝θ k-1 +N p ·ω k
in the formula, N p Is the number of pole pairs;
s20 quadrature-direct axis voltage u q,k 、u d,k The specific calculation formula is as follows,
u d,k ＝R·i d,k-1 +L·(i d,k-1 -i d,k-2 )/T s -N p ·ω k-1 ·L·i q,k-1
u q,k ＝R·i q,k-1 +L·(i q,k-1 -i q,k-2 )/T s +N p ·ω k-1 ·ψ f
wherein R is stator resistance, L is quadrature-direct axis inductance, psi f Is a rotor flux linkage;
quadrature axis current i of S21 q,k The specific calculation formula is as follows,
i q,k ＝T e,k /(N p ·ψ f )
direct axis current i of S22 d,k The specific calculation formula is as follows,
i d,k ＝(P L,k -u q,k ·i q,k )/u d,k
s23 motor three-phase voltage { u } a,k 、u b,k 、u c,k }, three-phase currents { i a,k 、i b,k 、i c,k The specific calculation formula is as follows,
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011105307.XA CN112257244B (en) | 2020-10-15 | 2020-10-15 | Modeling method of micro gas turbine power generation system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011105307.XA CN112257244B (en) | 2020-10-15 | 2020-10-15 | Modeling method of micro gas turbine power generation system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112257244A CN112257244A (en) | 2021-01-22 |
| CN112257244B true CN112257244B (en) | 2022-10-18 |
Family
ID=74243563
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011105307.XA Active CN112257244B (en) | 2020-10-15 | 2020-10-15 | Modeling method of micro gas turbine power generation system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112257244B (en) |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103490398B (en) * | 2013-10-17 | 2015-02-18 | 南车株洲电力机车研究所有限公司 | Direct axis current protection method and device for PMSM (permanent magnet synchronous motor) transmission system |
| CN108768232B (en) * | 2018-06-13 | 2019-12-20 | 青岛科技大学 | Control method and system for propulsion motor of semi-submersible vessel dynamic positioning system |
| CN109194229B (en) * | 2018-09-27 | 2020-07-17 | 北京理工大学 | Torque closed-loop-based MTPA control system and method for permanent magnet synchronous motor |
| CN110112781A (en) * | 2019-05-17 | 2019-08-09 | 河海大学 | A kind of comprehensive energy electric system distribution interconnection modeling and analysis methods |
| CN110502859B (en) * | 2019-08-28 | 2023-04-07 | 南方电网科学研究院有限责任公司 | Multi-rate dynamic simulation method for electric coupling park comprehensive energy system |
-
2020
- 2020-10-15 CN CN202011105307.XA patent/CN112257244B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN112257244A (en) | 2021-01-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102355006B (en) | Joint wind-solar-storage joint grid-connected power generation coordination control method | |
| CN107358345B (en) | Distributed combined cooling heating and power system optimization operation method considering demand side management | |
| Salameh et al. | A novel numerical simulation model for the PVT water system in the GCC region | |
| Kadri et al. | MPPT control methods in wind energy conversion system using DFIG | |
| Guan et al. | Optimal configuration and operation of multi-energy complementary distributed energy systems | |
| CN113255198A (en) | Multi-objective optimization method for combined cooling, heating and power supply micro-grid with virtual energy storage | |
| Hossain | Application of wind energy into the transportation sector | |
| CN103352802A (en) | Power generation device with vertical axis wind power and solar energy hot airflow complemented | |
| CN112257244B (en) | Modeling method of micro gas turbine power generation system | |
| Zhang et al. | Scheduling optimization for rural micro energy grid multi-energy flow based on improved crossbreeding particle swarm algorithm | |
| Bansal et al. | Blade pitch angle and tip speed ratio control schemes for constant power generation of WECS | |
| Zhou et al. | Digital Twin Modeling and Simulation of Distributed and Centralized Integrated Energy System on CloudPSS-IESLab | |
| Liu et al. | Effect of impeller solidity on the generating performance for solar power generation | |
| Chen et al. | A multi-energy microgrid modelling and optimization method based on exergy theory | |
| Rivero et al. | Wind park electric power estimation based on the Jensen wake model | |
| Kadri et al. | New grid connected hybrid generation system using wind turbine and solar Dish/Stirling engine | |
| CN216162639U (en) | Industrial wastewater and waste gas heat energy power device | |
| Jie et al. | The Multi-micronet Integrated Energy System Based on Multi-target Algorithm optimizes the Scheduling Strategy | |
| Kassem et al. | Power generation with different types of the vertical wind turbine for domestic use in northern cyprus: a case study | |
| Moorthy et al. | Wind energy recovery from a cooling tower with the help of a wind turbine generator | |
| Zheng et al. | Design and Performance Analysis of Distributed Equal Angle Spiral Vertical Axis Wind Turbine | |
| Kim et al. | Increasing the Efficiency of Thin-Film Silicon Solar Panels | |
| Shchur et al. | Stabilization of the coaxial counter-rotating vertical axis wind turbine via torque balancing by special double rotor PMSG | |
| Jia et al. | Research on integrated energy systems | |
| Hossain et al. | Zero-Emission Vehicles |
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 |