CN110334407B - Double-fed wind motor electromagnetic transient simulation method and system based on FPGA - Google Patents
Double-fed wind motor electromagnetic transient simulation method and system based on FPGA Download PDFInfo
- Publication number
- CN110334407B CN110334407B CN201910503950.9A CN201910503950A CN110334407B CN 110334407 B CN110334407 B CN 110334407B CN 201910503950 A CN201910503950 A CN 201910503950A CN 110334407 B CN110334407 B CN 110334407B
- Authority
- CN
- China
- Prior art keywords
- current
- simulation
- stator
- wind
- wind power
- 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.)
- Expired - Fee Related
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 142
- 230000001052 transient effect Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000004364 calculation method Methods 0.000 claims abstract description 61
- 238000007667 floating Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 21
- 238000003860 storage Methods 0.000 claims description 19
- 230000004907 flux Effects 0.000 claims description 11
- 238000013500 data storage Methods 0.000 claims description 10
- 238000002347 injection Methods 0.000 claims description 8
- 239000007924 injection Substances 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 238000004891 communication Methods 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- 239000005441 aurora Substances 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 238000011426 transformation method Methods 0.000 claims description 2
- 238000011161 development Methods 0.000 abstract description 11
- 238000005516 engineering process Methods 0.000 abstract description 10
- 230000001133 acceleration Effects 0.000 abstract description 2
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- 238000002054 transplantation Methods 0.000 abstract 1
- 239000000243 solution Substances 0.000 description 8
- 238000011160 research Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000003915 air pollution Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 238000010792 warming 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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Human Resources & Organizations (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Public Health (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Control Of Eletrric Generators (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses an electromagnetic transient simulation method and system of a double-fed wind power motor based on an FPGA (field programmable gate array), belonging to the technical field of power system simulation and the field of hardware calculation acceleration. According to the method, the decoupled double-fed wind power motor circuit is simulated on the FPGA in real time according to an equivalent circuit model of the wind power motor; the parameter updating module of the wind motor is quickly realized by utilizing an HLS high-level comprehensive technology, so that the development period of a hardware system is shortened; meanwhile, simulation calculation is carried out by using single-precision floating point numbers in the parameter updating module, so that the utilization rate of FPGA hardware resources is reduced while higher calculation precision is kept; by adopting Park conversion, repeated calculation of a parameter updating area is reduced, and hardware resource consumption is further reduced. The system framework of the invention is clear and definite, the transplantation and the expansion are convenient, the simulation performance is superior, and the real-time simulation can be realized, so the invention is particularly suitable for the electromagnetic transient real-time simulation calculation of the double-fed wind motor on the FPGA parallel framework platform.
Description
Technical Field
The invention relates to the field of fan electromagnetic transient simulation technology and hardware calculation acceleration, in particular to a double-fed wind motor electromagnetic transient simulation method and system based on an FPGA (field programmable gate array).
Background
Energy plays an important role in social progress and productivity improvement, but the global environment is increasingly worsened by the problems of air pollution, climate warming and the like caused by long-term large consumption of fossil energy, and the rapid development of renewable energy has become a consensus all over the world. Wind power is seen by all countries due to the characteristic of low cleanness, and all countries around the world are vigorously developing wind power to improve energy structures and global environments.
With the rapid development of the wind power generation industry, the traditional asynchronous generator can not meet the requirements of a wind power plant due to the defects that the rotating speed can not be changed, the operation efficiency is low and the like. A double-Fed Induction generator (DFIG) is gaining increasing favor due to a series of advantages such as variable speed and constant frequency operation characteristics, parameter capability and reactive power regulation of a system, and is currently becoming a mainstream generator of a wind farm.
With the continuous expansion of the scale of the wind power plant, the problems brought by the wind power network are gradually highlighted. On one hand, wind power has the characteristics of intermittence, randomness and the like, so that the instability of the operation of a main power grid is increased when the wind power is combined into the main power grid, and the problem of large-scale wind power grid connection stability becomes a bottleneck restricting the development of the wind power. On the other hand, although the modeling of the wind turbine and the wind farm has been preliminarily explored and researched in the industry, with the further development of the wind power technology and the improvement of the complexity of the wind farm, the existing modeling method is difficult to comprehensively and effectively analyze the wind power network. In the face of dilemma encountered by current wind power development and wide application prospects of wind power technology, related industrial workers strive to be put into the research of the wind power technology, so that modeling and simulation of wind turbines and wind power plants become hot spots of power grid simulation research.
In the real-time digital simulation of new energy and novel electronic equipment under the background of a large power grid, due to the improvement of the influence degree of a device on the system frequency, the requirement of the electromagnetic transient real-time simulation on the simulation step length reaches the microsecond order, which brings serious challenge to the design of an electromagnetic transient real-time simulation system. With the reduction of the simulation step length, the traditional electromagnetic transient simulation system based on the high-performance server is difficult to meet the real-time requirement, and people begin to search for a new electromagnetic transient simulation scheme. At present, most of domestic and foreign researches on modeling and simulation of wind power plants and fans are established on off-line simulation software based on a PC and subjected to simulation verification.
The reconfigurable and highly parallelized characteristics of an FPGA (Field Programmable Gate Array, FPGA for short) are well in line with the requirements of an electromagnetic transient simulation system, so that the FPGA is increasingly applied to the Field of electromagnetic transient simulation. At present, the simulation of the wind turbine based on the FPGA mainly focuses on the research of the operation control system of the wind turbine. Compared with a high-performance server, the FPGA can be combined with a hardware parallelization framework to enable an electromagnetic transient system to easily reach microsecond simulation step length, so that the FPGA is very suitable for being used as an electromagnetic transient simulation platform of a double-fed wind power motor. However, resources on the FPGA board are still limited, and in order to implement a larger-scale electromagnetic transient simulation system and maintain higher simulation accuracy under the same resources, a proper circuit model needs to be designed according to the characteristics of the FPGA architecture, so that hardware resources on the FPGA board can be efficiently utilized.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention aims to provide a double-fed wind power motor electromagnetic transient simulation method and system based on an FPGA.
1. And (3) design scheme is provided.
In view of the foregoing problems, the present invention designs an electromagnetic transient simulation method for a doubly-fed wind turbine generator based on an FPGA, which performs real-time simulation on the wind turbine generator under the condition of ensuring the single-precision floating point number precision of the IEEE754 standard, and rapidly implements a doubly-fed wind turbine generator parameter update module by using a High Level Synthesis (HLS) technique, thereby significantly shortening the development period of an FPGA hardware system.
According to the method, hardware resources on an FPGA board are efficiently utilized, and electromagnetic transient simulation of the doubly-fed wind power motor is divided into four modules: the wind power motor simulation calculation model has high simulation precision and can meet the simulation experiment requirements of more real application scenes.
The method is characterized in that Verilog hardware description language design development is used in a wind motor data storage module and a wind motor simulation control module, high-level comprehensive technology development is used in a wind motor parameter updating module, and the accuracy and the flexibility of a simulation system are considered.
2. In order to achieve the above object, the present invention adopts the following technical solutions.
An electromagnetic transient simulation method of a double-fed wind power motor based on an FPGA is characterized by comprising the following steps:
step 1-1, establishing a communication link between the FPGA and a server, and performing data interaction through an Aurora protocol;
step 1-2, the FPGA communicates with a server to obtain configuration parameters and initialization data required by simulation of a main circuit containing the doubly-fed wind motor;
wherein the configuration parameters and initialization data include:
-network node initialization voltage data of the simulation network;
-the number of simulated doubly-fed wind turbines;
-number information of each fan in the simulation network and of the network nodes;
-the number of masses per wind motor, the position of the mass on which the rotor is located, the stator resistance, the stator d-axis initial voltage, the stator q-axis initial voltage and the stator q-axis initial current;
step 1-3, acquiring initialization data of a simulation network from a system bus, and storing the number of network matrixes, the number of system subnets, simulation duration, a network initial node voltage array and a network initial node current array in a public storage area;
step 1-4, acquiring initialization data of a doubly-fed wind power motor from a system bus, storing the number of wind power motor elements, the simulation angular frequency, the number of wind power motor mass blocks and the positions of mass blocks of a rotor in a global variable region of a wind power motor data storage module, storing stator d-axis voltage, stator q-axis voltage, stator 0-axis voltage, stator d-axis rotating potential, stator d-axis current, stator q-axis rotating potential, stator q-axis current, rotor d-axis voltage, rotor q-axis voltage and rotor 0-axis voltage in a state variable region of the wind power motor data storage module, and storing a calculation angle of the wind power motor, electromagnetic torque of the wind power motor, mechanical torque of the wind power motor and mechanical power distribution proportion of each mass block of the wind power motor in a non-state variable region of the wind power motor data storage module;
step 3, calculating by a wind power motor parameter updating module according to the stator 0 shaft voltage and the stator 0 shaft impedance to obtain a stator 0 shaft current value of the current simulation step length;
step 4, updating the stator d-axis current and the stator q-axis current value of the current simulation step length by the wind motor parameter updating module;
step 5, the wind motor parameter updating module carries out rotation equation solving, and the concrete process is as follows:
step 5-1, updating and calculating the stator flux linkage and the rotor flux linkage of the current simulation step length;
step 5-2, updating and calculating the mechanical torque and the electromagnetic torque of the motor of the current simulation step length;
step 5-3, updating and calculating the rotor angular speed of the current simulation step length,
step 6, circularly executing the step 5-1 to the step 5-3 until the angular speed of the rotor with the current simulation step size is converged;
step 7, updating the rotor angular velocity predicted value of the current simulation step length by the wind motor parameter updating module;
step 9, updating the equivalent parallel current source value of the wind motor with the current simulation step length by a wind motor parameter updating module;
step 10, merging equivalent parallel current sources of all the wind motors with the current simulation step length by a wind motor node current merging module;
step 11, combining the equivalent branch current of the wind power motor with equivalent current sources of other elements to obtain updated network node current;
step 13, judging whether the current simulation step length reaches a set simulation end step length, if so, ending the simulation; otherwise, returning to the step 2 for iterative calculation.
In the step 1-4, the initialized data of the doubly-fed wind power generator stored on the FPGA uses the data of IEEE754 single-precision floating point standard.
Preferably, in step 2, the wind turbine parameter updating module is implemented by using a high-level integration technology.
Preferably, in step 2, when the stator d-axis voltage, q-axis voltage, and 0-axis voltage are calculated, the calculation is performed by using a hardware-optimized Park matrix transformation method.
Preferably, in step 5, an iterative calculation method is used to calculate a rotation equation of the wind turbine, so that the angular speed of the rotor of the wind turbine converges.
The simulation system comprises a CPU part and an FPGA part, a server of the CPU part comprises a public storage area and a wind power motor storage module, the FPGA part comprises an internal storage module, a control module, a fan parameter updating module and an internal current merging module, and the FPGA part is in communication connection with the server of the CPU part and performs node voltage solving and system current merging.
Preferably, the server of the CPU part receives the state variable initialization value of the wind turbine system bus to perform wind turbine initialization calculation, performs equivalent circuit calculation according to the initialization data to obtain equivalent circuit parameters, performs fan mass calculation according to the quality parameters to obtain mass parameters, performs fan constant coefficient calculation according to the voltage parameters, the rotation parameters, and the matrix parameters to obtain flux linkage parameters, and performs global initialization data calculation to obtain state parameters and non-state parameters; the FPGA part carries out variable updating calculation according to parameters calculated by a server of the CPU part to obtain a state variable and a non-state variable of the wind motor, and carries out injection current calculation according to the state variable and the non-state variable to obtain equivalent current and a new state variable so as to complete injection of the wind motor updating system into a current source.
Preferably, the wind motor storage module comprises a global variable area, a state variable area and a non-state variable area.
Compared with the prior art, the invention has the beneficial effects that:
1. the single-precision floating point number of IEEE754 standard is used in the electromagnetic transient simulation part of the wind motor, so that the simulation precision is ensured, the hardware resource consumption of the double-fed wind motor simulation is reduced, and the Park conversion is adopted in the calculation logic, so that the resource consumption caused by repeated calculation is reduced, and therefore, the invention can realize the real-time simulation of a larger number of double-fed wind motors on one FPGA board.
2. In an electromagnetic transient simulation system based on the FPGA, a high-level comprehensive technology is used for developing a wind motor parameter updating module, so that the design and development time of a hardware circuit is shortened while the simulation step of the wind motor is kept correct.
3. The simulation circuit has clear and definite calculation steps, is convenient to expand and has simple control logic.
Drawings
FIG. 1 is a simulation flow chart of a doubly-fed wind turbine;
FIG. 2 is an electromagnetic transient simulation architecture diagram of a doubly-fed wind turbine;
FIG. 3 is a fast implementation of a double-fed wind turbine hardware model;
FIG. 4 is a hardware computational optimization of Park matrix transformation;
fig. 5 is a flow chart of iterative calculation of a rotational equation.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
The electromagnetic transient simulation method and the electromagnetic transient simulation system of the double-fed wind power motor based on the FPGA, which are provided by the invention, are described by combining a hardware circuit structure and an operation working flow with reference to the figures 1-5:
As shown in FIG. 1, the main calculation flow of the wind turbine can be divided into three steps: the method comprises the steps of firstly, initializing and calculating the wind driven generator, secondly, calculating the state variable and the non-state variable of the wind driven generator, and thirdly, calculating the injection current source of the updating system of the wind driven generator. In the three steps, because the module in the first step is mainly used for calculating parameter information needed by the wind driven generator rotation equation, equivalent circuit solving and the like, and the delay influence factor for hardware simulation is small, the initialization part of the wind driven generator in the first step is calculated on a CPU (central processing unit) of a high-performance server side, on the basis, the calculated initialization data is transmitted to the FPGA through an optical fiber interface after being framed, and is stored in BRAM (block buffer management) on the FPGA in a blocking mode. Because the calculation logics in the second step and the third step are relatively dense, and the logics are all positioned in a key path of the hardware simulation of the wind driven generator, the calculation part of the state variable and the non-state variable of the wind driven generator and the current source injection module of the wind driven generator updating system are quickly integrated into a hardware IP (Internet protocol) by using a high-level integration technology and are integrated in an electromagnetic transient simulation system based on the FPGA (field programmable gate array).
In the simulation process, firstly, configuration parameters and initialization data are sent at a server end, and an FPGA end analyzes and stores the data, and the specific process is as follows:
and 1-1, establishing a communication link between the FPGA and a server, and performing data interaction through an Aurora protocol.
Step 1-2, the FPGA communicates with a server to obtain configuration parameters and initialization data required by simulation of a main circuit containing the doubly-fed wind motor; the configuration parameters and initialization data include: initializing voltage data of a network node of the simulation network; the number of the simulated double-fed wind motors is increased; the number information of each fan in the simulation network and the network node; the number of the mass blocks of each wind motor, the position of the mass block where the rotor is located, the resistance of the stator, the initial voltage of the d axis of the stator, the initial voltage of the q axis of the stator, the initial current of the q axis of the stator and the like.
Step 1-3, acquiring initialization data of the simulation network from a system bus, and storing the number of network matrixes, the number of system subnets, the simulation time, a network initial node voltage array, a network initial node current array and the like in a public storage area.
1-4, acquiring initialization data of a doubly-fed wind power motor from a system bus, storing parameters such as the number of elements of the wind power motor, the simulation angular frequency, the number of mass blocks of the wind power motor, the position of a mass block where a rotor is located in a global variable region of a wind power motor data storage module, storing parameters such as stator d-axis voltage, stator q-axis voltage, stator 0-axis voltage, stator d-axis rotating potential, stator d-axis current, stator q-axis rotating potential, stator q-axis current, rotor d-axis voltage, rotor q-axis voltage and rotor 0-axis voltage in a state variable region of the wind power motor data storage module, and storing parameters such as the calculation angle of the wind power motor, the electromagnetic torque of the wind power motor, the mechanical torque of the wind power motor, and the mechanical power distribution proportion of each mass block of the wind power motor in a non-state variable region of the wind power motor data storage module.
Fig. 3 shows the HLS development process of the parameter update module, and when designing, the C engineering model of the fan is first rewritten into a code style that the HLS tool can synthesize, and the rewritten code is verified for functional correctness; then, performing high-level synthesis of the C project by using an HLS tool, and optimizing the synthesized hardware code and verifying the cooperative consistency of software and hardware; finally, packaging the hardware codes into a hardware IP, and adding the hardware IP into a hardware IP library for a designer to call; and the wind power motor parameter updating module acquires network node voltage according to the node number of the wind power motor, and obtains the stator d-axis voltage, the stator q-axis voltage and the stator 0-axis voltage value of the current simulation step length through Park matrix transformation calculation.
As shown in fig. 4, when the hardware circuit performs Park matrix transformation, the coefficients outside the matrix are equivalent to the inside of the matrix, and constant elements are calculated in advance, and a constant λ is introduced 1 、λ 2 And λ 5. Besides the over-constant elements, two common elements lambda are extracted for the angle transformation calculation part contained in the matrix 3 And λ 4 And repeated calculation is reduced, and hardware resource consumption is reduced.
The Park matrix before optimization is:
introducing a constant λ 1 、λ 2 、λ 5 And variable lambda 3 、λ 4 Then, throughThe hardware-optimized Park matrix is:
wherein the content of the first and second substances,
and 3, calculating by the wind power motor parameter updating module according to parameters such as stator 0 shaft voltage, stator 0 shaft impedance and the like to obtain a stator 0 shaft current value of the current simulation step length.
And 4, updating and calculating the stator d-axis current and the stator q-axis current value of the current simulation step by the wind motor parameter updating module.
Step 5, the wind power motor parameter updating module carries out solution of the rotational equation, as shown in fig. 5, an iterative method is used for carrying out calculation solution, and the specific process is as follows:
step 5-1, updating and calculating the stator flux linkage and the rotor flux linkage of the current simulation step length;
step 5-2, updating and calculating the mechanical torque and the electromagnetic torque of the motor of the current simulation step length;
and 5-3, updating and calculating the rotor angular speed of the current simulation step length.
Specifically, the expression of the rotational equation is as follows:
in the formula, J is rotational inertia, omega is rotor angular velocity, D is damping coefficient considering viscous action and wind friction, and T is m And T e Respectively the mechanical torque and the electromagnetic torque of the motor.
As shown in fig. 5, when solving the rotational equation, firstly, flux linkage calculation and torque calculation are required; then, the calculated mechanical torque and electromagnetic torque are used for solving a rotational equation, and further the angular speed of the rotor is obtained; and then, judging the convergence of the rotating speed, if the convergence is judged, solving the rotating equation, if the convergence is not judged, returning to perform flux linkage calculation, torque calculation and rotating equation solving again in sequence, and judging the convergence of the rotating speed, and repeating the steps until the iteration is completed.
And 6, circularly executing the step 5-1 to the step 5-3 until the angular speed of the rotor with the current simulation step size is converged.
And 7, updating the rotor angular speed predicted value of the current simulation step length by the wind motor parameter updating module.
And 9, updating the equivalent parallel current source value of the wind motor with the current simulation step length by the wind motor parameter updating module.
And step 10, merging equivalent parallel current sources of all the wind motors with the current simulation step length by a wind motor node current merging module.
And step 11, combining the equivalent branch current of the wind power motor with equivalent current sources of other elements to obtain updated network node current.
And step 12, performing matrix vector multiplication calculation on the network admittance matrix and the network node current to obtain an updated network node voltage value.
Step 13, judging whether the current simulation step length reaches a set simulation end step length, if so, ending the simulation; otherwise, returning to the step 2 for iterative calculation.
Referring to fig. 1, the invention further provides a simulation system for implementing the electromagnetic transient simulation method, the simulation system includes a CPU portion and an FPGA portion, a server of the CPU portion includes a common storage area and a wind turbine storage module, the FPGA portion includes an internal storage module, a control module, a fan parameter update module and an internal current merge module, and the FPGA portion is in communication connection with the server of the CPU portion and performs node voltage solving and system current merging.
The wind power motor storage module comprises a global variable area, a state variable area and a non-state variable area.
Further, a server of the CPU part receives a state variable initialization value of a wind motor system bus to perform wind motor initialization calculation, performs equivalent circuit calculation according to initialization data to obtain equivalent circuit parameters, performs fan mass calculation according to quality parameters to obtain mass parameters, performs fan constant coefficient calculation according to voltage parameters, rotation parameters and matrix parameters to obtain flux linkage parameters, and performs global initialization data calculation to obtain state parameters and non-state parameters; the FPGA part carries out variable updating calculation according to parameters calculated by a server of the CPU part to obtain a state variable and a non-state variable of the wind motor, and carries out injection current calculation according to the state variable and the non-state variable to obtain equivalent current and a new state variable so as to complete injection of the wind motor updating system into a current source.
The system carries out electromagnetic transient simulation analysis according to the method provided by the application, can realize a double-fed wind power motor simulation calculation model with higher simulation precision by utilizing limited hardware resources on an FPGA (field programmable gate array) board, and meets the simulation experiment requirements of a large number of real application scenes.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. An electromagnetic transient simulation method of a double-fed wind power motor based on an FPGA is characterized by comprising the following steps:
step 1, a server side sends configuration parameters and initialization data, and an FPGA side analyzes and stores the data, and the specific process is as follows:
step 1-1, establishing a communication link between the FPGA and a server, and performing data interaction through an Aurora protocol;
step 1-2, the FPGA communicates with a server to obtain configuration parameters and initialization data required by simulation of a main circuit containing the doubly-fed wind motor;
wherein the configuration parameters and initialization data include:
-network node initialization voltage data of the simulation network;
-the number of simulated doubly-fed wind turbines;
-number information of each fan in the simulation network and of the network nodes;
-the number of masses per wind motor, the position of the mass on which the rotor is located, the stator resistance, the stator d-axis initial voltage, the stator q-axis initial voltage and the stator q-axis initial current;
step 1-3, acquiring initialization data of a simulation network from a system bus, and storing the number of network matrixes, the number of system subnets, simulation duration, a network initial node voltage array and a network initial node current array in a public storage area;
step 1-4, acquiring initialization data of a doubly-fed wind power motor from a system bus, storing the number of elements of the wind power motor, the simulation angular frequency, the number of mass blocks of the wind power motor and the position of the mass block where a rotor is located in a global variable region of a wind power motor data storage module, storing stator d-axis voltage, stator q-axis voltage, stator 0-axis voltage, stator d-axis rotating potential, stator d-axis current, stator q-axis rotating potential, stator q-axis current, rotor d-axis voltage, rotor q-axis voltage and rotor 0-axis voltage in a state variable region of the wind power motor data storage module, and storing a calculation angle of the wind power motor, electromagnetic torque of the wind power motor, mechanical torque of the wind power motor and mechanical power distribution proportion of each mass block of the wind power motor in a non-state variable region of the wind power motor data storage module;
step 2, the electromagnetic transient simulation part of the doubly-fed wind power motor comprises four modules: the system comprises an internal storage module, a control module, a wind power motor parameter updating module and an internal current merging module; the wind power motor parameter updating module acquires network node voltage according to the node number, and obtains stator d-axis voltage, stator q-axis voltage and stator 0-axis voltage values of the wind power motor in the current simulation step length through Park matrix transformation calculation;
step 3, calculating by the wind motor parameter updating module according to the stator 0 shaft voltage and the stator 0 shaft impedance to obtain a stator 0 shaft current value of the current simulation step length;
step 4, updating the stator d-axis current and the stator q-axis current value of the current simulation step length by the wind motor parameter updating module;
step 5, the wind motor parameter updating module carries out rotation equation solving, and the concrete process is as follows:
step 5-1, updating and calculating the stator flux linkage and the rotor flux linkage of the current simulation step length;
step 5-2, updating and calculating the mechanical torque and the electromagnetic torque of the motor with the current simulation step length;
step 5-3, updating and calculating the rotor angular speed of the current simulation step length,
step 6, circularly executing the step 5-1 to the step 5-3 until the angular speed of the rotor of the current simulation step is converged;
step 7, updating the rotor angular velocity predicted value of the current simulation step length by the wind motor parameter updating module;
step 8, calculating by the wind power motor parameter updating module according to the d-axis current of the stator, the q-axis current of the stator and the 0-axis current of the stator to obtain a stator a-axis current, a stator b-axis current and a stator c-axis current value of the current simulation step length; updating and calculating the rotor three-phase current value of the current simulation step length according to the rotor angle and the rotor current;
step 9, updating the equivalent parallel current source value of the wind motor with the current simulation step length by a wind motor parameter updating module;
step 10, a wind power motor node current merging module merges equivalent parallel current sources of all wind power motors with the current simulation step length;
step 11, combining the equivalent branch current of the wind power motor with equivalent current sources of other elements to obtain updated network node current;
step 12, performing matrix vector multiplication calculation on the network admittance matrix and the network node current to obtain an updated network node voltage value;
step 13, judging whether the current simulation step length reaches a set simulation end step length, if so, ending the simulation; otherwise, returning to the step 2 for iterative computation.
2. The electromagnetic transient simulation method of claim 1, wherein the initialization data of the doubly-fed wind turbine stored on the FPGA in step 1-4 is data of IEEE754 single precision floating point number standard.
3. The electromagnetic transient simulation method of claim 1, wherein in step 2 the wind turbine parameter update module is implemented using high level synthesis techniques.
4. The electromagnetic transient simulation method of claim 3, wherein in step 2, the wind turbine parameter updating module adopts pipeline instruction optimization during high-level integrated design, so as to reduce the delay time calculated by the module.
5. The electromagnetic transient simulation method of claim 2, wherein in step 2, the wind turbine parameter updating module adopts a multi-group parallel strategy, so as to shorten the calculation delay when the number of wind turbines is large.
6. The electromagnetic transient simulation method of claim 1, wherein in step 2, the stator d-axis voltage, q-axis voltage and 0-axis voltage are calculated by using a hardware-optimized Park matrix transformation method.
7. The electromagnetic transient simulation method of claim 1, wherein the rotational equations of the wind turbine are solved in step 5 using iterative calculations.
8. A simulation system for implementing the electromagnetic transient simulation method according to any one of claims 1 to 7, characterized by: the simulation system comprises a CPU part and an FPGA part, a server of the CPU part comprises a public storage area and a wind power generator storage module, the FPGA part comprises an internal storage module, a control module, a fan parameter updating module and an internal current combining module, and the FPGA part is in communication connection with the server of the CPU part and performs node voltage solving and system current combining.
9. The simulation system of claim 8, wherein: the method comprises the following steps that a server of the CPU part receives a state variable initialization value of a wind motor system bus to carry out wind motor initialization calculation, equivalent circuit calculation is carried out according to initialization data to obtain equivalent circuit parameters, fan mass block calculation is carried out according to quality parameters to obtain mass block parameters, fan constant coefficient calculation is carried out according to voltage parameters, rotation parameters and matrix parameters to obtain flux linkage parameters, and global initialization data calculation is carried out to obtain state parameters and non-state parameters; the FPGA part carries out variable updating calculation according to parameters calculated by a server of the CPU part to obtain a state variable and a non-state variable of the wind motor, and carries out injection current calculation according to the state variable and the non-state variable to obtain equivalent current and a new state variable so as to complete injection of the wind motor updating system into a current source.
10. The simulation system of claim 8, wherein: the wind motor storage module comprises a global variable area, a state variable area and a non-state variable area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910503950.9A CN110334407B (en) | 2019-06-12 | 2019-06-12 | Double-fed wind motor electromagnetic transient simulation method and system based on FPGA |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910503950.9A CN110334407B (en) | 2019-06-12 | 2019-06-12 | Double-fed wind motor electromagnetic transient simulation method and system based on FPGA |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110334407A CN110334407A (en) | 2019-10-15 |
| CN110334407B true CN110334407B (en) | 2022-11-04 |
Family
ID=68141115
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910503950.9A Expired - Fee Related CN110334407B (en) | 2019-06-12 | 2019-06-12 | Double-fed wind motor electromagnetic transient simulation method and system based on FPGA |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110334407B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110968973B (en) * | 2019-11-11 | 2023-06-16 | 上海科梁信息科技股份有限公司 | Simulation method, control method, electronic equipment and storage medium of transformer model |
| CN111460648B (en) * | 2020-03-30 | 2022-06-07 | 国网福建省电力有限公司电力科学研究院 | Electromagnetic transient real-time simulation data communication method for power system |
| CN112636648B (en) * | 2020-12-11 | 2022-07-08 | 上海科梁信息科技股份有限公司 | Modeling method and system of double-fed fan based on FPGA |
| CN113987821A (en) * | 2021-11-04 | 2022-01-28 | 上海远宽能源科技有限公司 | Multi-type motor real-time simulation method and system based on FPGA |
| CN115034087A (en) * | 2022-07-05 | 2022-09-09 | 南方电网电力科技股份有限公司 | Modeling method and device for wind turbine generator, electronic equipment and storage medium |
| CN115598994A (en) * | 2022-09-07 | 2023-01-13 | 贵州电网有限责任公司（Cn） | Multi-rate real-time simulation method of permanent magnet direct drive type wind power generation system |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104123172A (en) * | 2014-07-17 | 2014-10-29 | 南车株洲电力机车研究所有限公司 | Double-fed wind power simulation system and circuit simulation module thereof |
| CN104317577A (en) * | 2014-10-08 | 2015-01-28 | 陕西高新实业有限公司 | Design method for realizing openCV by using ViadoHLS in FPGA (Field Programmable Gate Array) development |
| CN106250614A (en) * | 2016-07-29 | 2016-12-21 | 上海交通大学 | It is applicable to the data processing method of FPGA platform electro-magnetic transient real-time simulation |
| CN106844900A (en) * | 2016-12-31 | 2017-06-13 | 上海交通大学 | The erection method of electromagnetic transient simulation system |
| US9690894B1 (en) * | 2015-11-02 | 2017-06-27 | Altera Corporation | Safety features for high level design |
| CN107395074A (en) * | 2017-07-24 | 2017-11-24 | 吉林大学 | The permanent magnet DC motor multiresolution controller realized based on FPGA |
| CN108549758A (en) * | 2018-04-03 | 2018-09-18 | 上海交通大学 | A kind of modularization multi-level converter real-time emulation method suitable for FPGA |
| CN108768220A (en) * | 2018-07-23 | 2018-11-06 | 南京工程学院 | A kind of integrative coordinated control device and method of articulated robot |
| CN109800455A (en) * | 2018-12-13 | 2019-05-24 | 中国电力科学研究院有限公司 | A kind of idle characteristic Simulation method and apparatus of double-fed fan motor unit transient state |
-
2019
- 2019-06-12 CN CN201910503950.9A patent/CN110334407B/en not_active Expired - Fee Related
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104123172A (en) * | 2014-07-17 | 2014-10-29 | 南车株洲电力机车研究所有限公司 | Double-fed wind power simulation system and circuit simulation module thereof |
| CN104317577A (en) * | 2014-10-08 | 2015-01-28 | 陕西高新实业有限公司 | Design method for realizing openCV by using ViadoHLS in FPGA (Field Programmable Gate Array) development |
| US9690894B1 (en) * | 2015-11-02 | 2017-06-27 | Altera Corporation | Safety features for high level design |
| CN106250614A (en) * | 2016-07-29 | 2016-12-21 | 上海交通大学 | It is applicable to the data processing method of FPGA platform electro-magnetic transient real-time simulation |
| CN106844900A (en) * | 2016-12-31 | 2017-06-13 | 上海交通大学 | The erection method of electromagnetic transient simulation system |
| CN107395074A (en) * | 2017-07-24 | 2017-11-24 | 吉林大学 | The permanent magnet DC motor multiresolution controller realized based on FPGA |
| CN108549758A (en) * | 2018-04-03 | 2018-09-18 | 上海交通大学 | A kind of modularization multi-level converter real-time emulation method suitable for FPGA |
| CN108768220A (en) * | 2018-07-23 | 2018-11-06 | 南京工程学院 | A kind of integrative coordinated control device and method of articulated robot |
| CN109800455A (en) * | 2018-12-13 | 2019-05-24 | 中国电力科学研究院有限公司 | A kind of idle characteristic Simulation method and apparatus of double-fed fan motor unit transient state |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110334407A (en) | 2019-10-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110334407B (en) | Double-fed wind motor electromagnetic transient simulation method and system based on FPGA | |
| Ciri et al. | Model-free control of wind farms: A comparative study between individual and coordinated extremum seeking | |
| CN104600756B (en) | A kind of group of planes equivalent modeling method of middle-size and small-size hydro-generating Unit | |
| CN103810646B (en) | Improved projection integral algorithm based active power distribution system dynamic simulation method | |
| CN103700036B (en) | A kind of transient stability projecting integral method being suitable to power system Multiple Time Scales | |
| Ni et al. | Model order reduction based dynamic equivalence of a wind farm | |
| CN105591402B (en) | A kind of modeling and simulation method and apparatus of direct-drive permanent-magnetism Wind turbines | |
| Khaouch et al. | Mechatronic modeling of a 750 kW fixed-speed wind energy conversion system using the Bond Graph Approach | |
| CN103678900B (en) | A kind of network decoupling calculation method for regional distribution network real-time simulation | |
| CN103219725A (en) | Wind power plant equivalent modeling method based on real-time operation data | |
| CN107947228B (en) | Stochastic stability analysis method for power system containing wind power based on Markov theory | |
| CN204314716U (en) | A kind of for wind power station control system hardware-in―the-loop test platform | |
| CN106250614A (en) | It is applicable to the data processing method of FPGA platform electro-magnetic transient real-time simulation | |
| CN106777499B (en) | Complete machine dynamic modeling method for double-fed asynchronous wind generating set | |
| KR20220100975A (en) | Cable route planning methods, systems, media and electronic devices in wind farms | |
| CN108365620B (en) | Method and device for simulating energy storage system based on equivalent model | |
| CN110765594B (en) | Method and device for controlling internal interaction interface of wind power plant real-time simulator | |
| Moussa et al. | kW wind turbine emulator implementation on FPGA using matlab/simulink | |
| CN113673114A (en) | Collaborative coupling solving method for power distribution network information physical system | |
| Mohammadi et al. | Platform for design, simulation, and experimental evaluation of small wind turbines | |
| CN116467877A (en) | Floating wind turbine generator system platform dynamic response determination method and device and electronic equipment | |
| Tang et al. | Creating multi-timescale digital twin models for regional multiple energy systems on CloudPSS | |
| Zhang et al. | Review of modeling methods in electromagnetic and thermal design of permanent magnet generators for wind turbines | |
| CN110597084A (en) | Wind power plant simulation test method and device and wind power plant simulation model | |
| Solcia et al. | Efficient multirate simulation of complex multibody systems based on free software |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20221104 |