Simulation calculation method and system for double-wind-wheel wind turbine generator

Abstract

The invention discloses a double-wind-wheel wind turbine generator simulation calculation method and system, and the method comprises the steps: calculating a blade target pitch angle and a target generator torque of a front wind wheel and a rear wind wheel based on basic parameters of the front wind wheel and the rear wind wheel, and generating a control instruction to control a double-wind-wheel wind turbine generator simulation model; according to the target blade pitch angles and the target generator torques of the front and rear wind wheels, calculating blade aerodynamic forces and structural loads of the front and rear wind wheels in the double-wind-wheel wind turbine generator simulation model; obtaining blade aerodynamic moments of the front and rear wind wheels based on the blade aerodynamic forces of the front and rear wind wheels, and calculating to obtain current rotating speeds of the front and rear wind wheels according to the blade aerodynamic moments of the front and rear wind wheels, the target generator torque, the rotational inertia and the rotational inertia of the transmission system for updating basic parameters of the front and rear wind wheels; and repeating the above operation to iteratively calculate the blade aerodynamic force and the structural load of the front and rear wind wheels until the simulation duration reaches a preset value, and completing the time domain simulation of the full-coupling dynamics of the double-wind-wheel wind turbine generator.

Description

Technical Field

[0001] The present application relates to the technical field of wind turbines, and in particular to a simulation calculation method and system for a dual-rotor wind turbine. Background Art

[0002] The tandem twin-rotor wind turbine is a new type of wind energy utilization equipment. By installing a wind rotor in front and behind the nacelle, the rear wind rotor can further utilize the wake wind energy of the front wind rotor, thereby improving the wind energy utilization efficiency of the entire wind turbine and reducing the overall wind energy utilization cost. It has great application prospects.

[0003] There is a significant interaction between the front and rear wind rotors of a tandem twin-rotor wind turbine. The rotation of the front wind rotor will produce an uneven wake, which affects the aerodynamic performance of the rear wind rotor. The aerodynamic force generated by the rear wind rotor will cause greater deformation and vibration of the tower, thus affecting the aerodynamic load of the front wind rotor. Therefore, the interaction between the front and rear twin wind rotors needs to be considered when designing the structure of a tandem twin-rotor wind turbine.

[0004] In the prior art, the calculation of aerodynamic loads for tandem twin-rotor wind turbines generally adopts computational fluid dynamics (CFD) methods. By establishing the flow field of the tandem twin-rotor wind turbine and combining the turbulence model, the NS equations of the flow field are solved to obtain the pressure distribution on the surfaces of the front and rear twin-rotor blades, thereby calculating the aerodynamic forces of the front and rear twin-rotors.

[0005] The CFD method can consider the impact of the front wind rotor wake on the aerodynamic performance of the rear wind rotor, but the CFD method cannot consider the blade deformation and the coordinated control of the front and rear wind rotors, and ignores the blade elasticity, the rotor and the nacelle inertia force, etc., and cannot accurately obtain the dynamic response of the dual-rotor wind turbine under turbulent wind conditions. In addition, the CFD method requires huge computing resources, and the computing efficiency is difficult to meet engineering applications. Summary of the invention

[0006] The purpose of the present application is to provide a simulation calculation method and system for a dual-rotor wind turbine set, taking into account the deformation of blades and the coordinated control of front and rear wind wheels, and realizing the simulation calculation of the aerodynamic loads of the front and rear wind wheels.

[0007] To achieve the above objectives, this application provides the following technical solutions:

[0008] In a first aspect, an embodiment of the present application provides a dual-rotor wind turbine simulation calculation based on the basic parameters of the front and rear wind turbines, calculates the target blade pitch angles of the front and rear wind turbines and the target generator torques of the front and rear wind turbines, and generates control instructions to control the dual-rotor wind turbine simulation model;

[0009] According to the target blade pitch angles of the front and rear wind rotors and the target generator torques of the front and rear wind rotors, the blade aerodynamic forces and structural loads of the front and rear wind rotors in the simulation model of the double-rotor wind turbine are calculated;

[0010] The blade aerodynamic moments of the front and rear wind wheels are obtained based on the blade aerodynamic forces of the front and rear wind wheels, and the current rotation speeds of the front and rear wind wheels in the basic parameters of the front and rear wind wheels are calculated and obtained according to the blade aerodynamic moments of the front and rear wind wheels, the target generator torque, the moment of inertia and the moment of inertia of the transmission system, so as to update the basic parameters of the front and rear wind wheels;

[0011] Repeat the above operations to iteratively calculate the blade aerodynamic force and structural load of the front and rear wind rotors until the simulation time reaches the preset value, completing the time domain simulation of the fully coupled dynamics of the dual-rotor wind turbine set;

[0012] Among them, the dual-rotor wind turbine simulation model includes a front wind wheel, a rear wind wheel, a nacelle and a tower. The front wind wheel and the rear wind wheel are respectively connected to the two ends of the nacelle, and the nacelle is placed on the top of the tower; the diameter of the front wind wheel is larger than that of the rear wind wheel.

[0013] In a second aspect, an embodiment of the present application provides a dual-wind-rotor wind turbine simulation computing system, comprising:

[0014] A calculation control unit, used to calculate the target pitch angles and target generator torques corresponding to the front and rear wind rotors based on the basic parameters of the front and rear wind rotors, and generate control instructions to control the simulation model of the dual-wind-rotor wind turbine set;

[0015] A calculation unit, used to calculate the blade aerodynamic force and structural load of the front and rear wind rotors in the simulation model of the double-wind-rotor wind turbine according to the target pitch angles and target generator torques corresponding to the front and rear wind rotors;

[0016] A parameter updating unit is used to obtain the blade aerodynamic torque of the front and rear wind wheels based on the blade aerodynamic force of the front and rear wind wheels, and calculate the current rotation speed of the front and rear wind wheels in the basic parameters of the front and rear wind wheels according to the blade aerodynamic torque of the front and rear wind wheels, the target generator torque, the moment of inertia and the moment of inertia of the transmission system, so as to update the basic parameters of the front and rear wind wheels;

[0017] The repeated simulation unit is used to repeatedly calculate the aerodynamic force and structural load of the front and rear wind rotors until the simulation time reaches the set value, thus completing the time domain simulation of the fully coupled dynamics of the dual-rotor wind turbine set;

[0018] Among them, the dual-rotor wind turbine simulation model also includes a front wind wheel, a rear wind wheel, a nacelle and a tower. The front wind wheel and the rear wind wheel are respectively connected to the two ends of the nacelle, and the nacelle is placed on the top of the tower; the diameter of the front wind wheel is larger than that of the rear wind wheel.

[0019] In a third aspect, an embodiment of the present application further provides an electronic device, comprising a memory, a processor and a computer program stored in the memory, wherein the processor executes the computer program or instructions to implement the steps of the aforementioned dual-rotor wind turbine simulation calculation method.

[0020] In a fourth aspect, an embodiment of the present application further provides a computer storage medium, in which a computer program or instruction is stored. When the computer program or instruction is executed by a processor, the steps of the aforementioned dual-rotor wind turbine simulation calculation method are implemented.

[0021] In a fifth aspect, an embodiment of the present application further provides a computer program product, including a computer program or instructions, which, when executed by a processor, implement the steps of the aforementioned dual-rotor wind turbine simulation calculation method.

[0022] Technical effects and advantages of the present application: The simulation calculation method of the present application can not only consider the influence of the wake of the front wind wheel on the aerodynamic force of the rear wind wheel, but also calculate the elastic deformation of the blades and tower under the premise of considering the control effect. The present application overcomes the problem that the traditional CFD method cannot consider the elasticity of the blades, and at the same time greatly improves the aerodynamic performance evaluation of the dual-rotor wind turbine and the calculation efficiency of the dynamic response of the whole machine.

[0023] Other features and advantages of the present application will be described in the following description, and partly become apparent from the description, or be understood by practicing the present application. The purpose and other advantages of the present application can be realized and obtained by the structures indicated in the description, claims and drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0024] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0025] Figure 1 This is a flow chart of a simulation calculation method for a dual-wind-rotor wind turbine set according to an embodiment of the present application;

[0026] Figure 2 This is a schematic diagram of a model of a double-wind-rotor wind turbine generator system according to an embodiment of the present application;

[0027] Figure 3 This is a schematic diagram of the structure of a simulation calculation system for a dual-wind-rotor wind turbine set according to an embodiment of the present application;

[0028] Figure 4 This is a schematic diagram of the structure of an electronic device according to an embodiment of the present application;

[0029] Figure 5 A schematic diagram showing a comparison of wind wheel powers calculated by the method of the present application and the CFD method in an embodiment of the present application;

[0030] Figure 6 A schematic diagram comparing the wind wheel thrust calculated by the method of the present application and the CFD method in the embodiment of the present application;

[0031] Figure 7 This is a schematic diagram comparing the calculation time consumption of the method of the present application and the CFD method in the embodiments of the present application;

[0032] Figure 8 Schematic diagram of the comparison of the wind wheel thrust of the front wind wheel and the rear wind wheel under turbulent wind conditions in an embodiment of the present application;

[0033] Fig. 9 This is a schematic diagram comparing the power generation of the front wind wheel and the rear wind wheel under turbulent wind conditions in an embodiment of the present application;

[0034] Fig.10 This is a schematic diagram of wind speed comparison at the hub of the front wind wheel and the rear wind wheel under turbulent wind conditions in an embodiment of the present application;

[0035] Fig.11 This is a schematic diagram comparing the horizontal thrust at the top of the tower of a double wind rotor and a single wind rotor under the same working condition in an embodiment of the present application;

[0036] Fig.12 This is a schematic diagram comparing the vertical thrust at the top of the tower of a double wind rotor and a single wind rotor under the same working condition in an embodiment of the present application. DETAILED DESCRIPTION

[0037] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0038] In order to solve the deficiencies of the prior art, the present application discloses a simulation calculation method for a double-wind-rotor wind turbine. Figure 1 As shown, the following steps are included:

[0039] S1: Based on the basic parameters of the front and rear wind rotors, the target blade pitch angles of the front and rear wind rotors and the target generator torques of the front and rear wind rotors are calculated, and control instructions are generated to control the simulation model of the dual-rotor wind turbine set;

[0040] S2: Calculate the aerodynamic force and structural load of the front and rear wind rotors in the simulation model of the dual-rotor wind turbine according to the target blade pitch angles of the front and rear wind rotors and the target generator torques of the front and rear wind rotors;

[0041] S3: obtaining the blade aerodynamic torque of the front and rear wind wheels based on the blade aerodynamic forces of the front and rear wind wheels, and calculating the current rotation speeds of the front and rear wind wheels in the basic parameters of the front and rear wind wheels according to the blade aerodynamic torque of the front and rear wind wheels, the target generator torque, the moment of inertia and the moment of inertia of the transmission system for updating the basic parameters of the front and rear wind wheels;

[0042] Step S4: repeating the operations of step S1 to step S3 to iteratively calculate the blade aerodynamic force and structural load of the front and rear wind rotors until the simulation time reaches the set value, completing the time domain simulation of the fully coupled dynamics of the dual-wind-rotor wind turbine set;

[0043] Among them, Figure 2 As shown, the simulation model of the dual-rotor wind turbine set includes a front wind rotor, a rear wind rotor, a nacelle and a tower, the front wind rotor and the rear wind rotor are respectively connected to the two ends of the nacelle, and the nacelle is placed on the top of the tower; the diameter of the front wind rotor is larger than that of the rear wind rotor; the simulation model of the dual-rotor wind turbine set also includes a transmission system, which includes an engine, a gearbox and a wind rotor rotating shaft.

[0044] In some specific embodiments, before calculating the target blade pitch angles of the front and rear wind rotors and the target generator torques of the front and rear wind rotors for the first time, setting the design parameters of the dual-wind-rotor wind turbine simulation model;

[0045] The design parameters of the dual-rotor wind turbine simulation model include: aerodynamic calculation parameters and structural dynamic calculation parameters of the wind rotor blades;

[0046] The wind rotor blade aerodynamic calculation parameters include: chord length, twist angle and airfoil lift and drag coefficient at different radii of the front and rear wind rotor blades;

[0047] The structural dynamics calculation parameters include: modal vibration shapes of each order mode of the blade and the corresponding damping ratio, cross-sectional mass at different radii of the blade, flapping stiffness at different radii of the blade, swing stiffness at different radii of the blade, hub position coordinates and hub mass of the front and rear wind wheels, moment of inertia of the transmission system, cabin mass, cabin moment of inertia, tower mass, distribution information of tower stiffness along the tower height, modal vibration shape of the tower, and damping ratio of the tower.

[0048] In some specific embodiments, step S1: based on the basic parameters of the front and rear wind rotors, calculating the target blade pitch angles of the front and rear wind rotors and the target generator torques of the front and rear wind rotors, and generating control instructions to control the dual-wind rotor wind turbine simulation model, includes the following steps:

[0049] The basic parameters of the front and rear wind wheels include: the current blade root loads of the front and rear wind wheels, the current blade deformations of the front and rear wind wheels, the current rotation speeds of the front and rear wind wheels, and the current generator torques of the front and rear wind wheels;

[0050] The basic parameters of the front and rear wind rotors are input into the simulation model control system (i.e., the control system of the dual-rotor wind turbine simulation model), and the target pitch angles and target generator torques corresponding to the front and rear wind rotors are obtained through algorithm calculation based on the state space of the entire dual-rotor wind turbine;

[0051] The simulation model control system generates a control instruction according to the target pitch angle and the target generator torque to control the dual-rotor wind turbine simulation model;

[0052] Among them, the control closed loop of the front wind wheel needs to consider the influence factors of the rear wind wheel and the tower, and the control closed loop of the rear wind wheel needs to consider the influence factors of the front wind wheel and the tower;

[0053] The influencing factors of the rear wind wheel include the rear wind wheel speed, the rear wind wheel power generation and the rear wind wheel thrust; the influencing factors of the front wind wheel include the front wind wheel speed, the front wind wheel power generation and the front wind wheel thrust; the influencing factors of the tower include the tower motion state and the tower load.

[0054] In some specific embodiments, step S2: calculating the blade aerodynamic force and structural load of the front and rear wind rotors in the dual-wind-rotor wind turbine simulation model according to the target pitch angles and target generator torques corresponding to the front and rear wind rotors, comprises the following steps:

[0055] Step S21: According to the defined front wind rotor wind field environment data, respectively calculate the blade aerodynamic force of the front wind rotor and the near-field wake coordinates of the front wind rotor, specifically including:

[0056] The front wind wheel wind field environmental data can cover the time variation data of the three directional components of the wind speed in the plane of the front wind wheel and the tower; the type of the front wind wheel wind field is turbulent wind or steady wind or other wind that satisfies the wind shear law.

[0057] (1) Based on the defined front wind rotor wind field environment data, the aerodynamic force of the front wind rotor blades is calculated using the free vortex wake method and blade element momentum theory, including:

[0058] According to the front wind rotor wind field environment data, a real-time inflow velocity at each blade cross section of the front wind rotor is obtained, and based on the real-time inflow velocity and the current rotation speed of the front wind rotor, an induced velocity at each blade cross section control point of the front wind rotor is obtained by solving the free vortex wake method;

[0059] Determine the velocity triangle at each blade section of the front wind wheel according to the initial rotation speed of the front wind wheel and the induced speed, and calculate the actual angle of attack of each blade section of the front wind wheel based on the velocity triangle;

[0060] According to the lift and drag coefficient of the front wind rotor in the blade aerodynamic calculation parameters, the lift coefficient and drag coefficient at the actual angle of attack of each blade section are retrieved, and the aerodynamic force of the corresponding blade section in the front wind rotor is calculated using the blade element momentum theory;

[0061] The aerodynamic force of each blade cross section of the front wind wheel is superimposed to obtain the blade aerodynamic force of the entire front wind wheel.

[0062] (2) Calculate the near-field wake coordinates of the front wind rotor according to the defined front wind rotor wind field environment data, including:

[0063] According to the front wind wheel wind field environment data, the real-time inflow velocity at each blade section of the front wind wheel is obtained, and the induced velocity at any wake behind the front wind wheel plane is obtained by solving the free vortex wake method;

[0064] The induced velocity field at any wake behind the front wind wheel plane is superimposed with the environmental space field of the real-time inflow velocity to obtain the three-dimensional velocity space coordinates of the near-field wake of the front wind wheel, that is, the coordinates of the near-field wake of the front wind wheel are obtained as follows:

[0065]

[0066] In the formula, represents the three-dimensional velocity space coordinates of the near-field wake of the front wind rotor, It represents the induced speed at any wake behind the front wind rotor plane. Usually this value is negative, indicating the amount of wind speed drop. Indicates the real-time inflow velocity at each blade section of the front wind wheel.

[0067] Step S22: Considering the influence of the rear wind rotor, the dynamic response of the front wind rotor blades and the tower structure is calculated to obtain the displacement, velocity and acceleration of the tower top; specifically including:

[0068] The Kane method is used to establish a front wind rotor multi-body system dynamics model; wherein the front wind rotor multi-body system dynamics model is a multi-body system dynamics model of flexible bodies such as front wind rotor blades and towers, and rigid bodies such as front wind rotor hubs and nacelles;

[0069] Considering the gravity and motion effects of the remaining structures except the wind rotor blades, as well as the load effect of the rear wind rotor, the dynamic response of the front wind rotor blades and the tower structure is calculated using the front wind rotor multi-body system dynamics model to obtain the displacement, velocity and acceleration of the tower top.

[0070] Among them, except for the wind rotor blades, the remaining structures include: transmission system, nacelle, wind rotor hub and tower structure;

[0071] The load action of the rear wind wheel is the combined action of the rear wind wheel aerodynamic force, the rear wind wheel inertia force and the rear wind wheel elastic force, and does not include the action of the cabin gravity and the cabin inertia force;

[0072] When calculating the dynamic response of the front wind rotor blades and tower structure, the modal method is used for the structural deformation of the wind rotor blades and tower, thereby accelerating the efficiency of structural dynamics solution.

[0073] Step S23: Based on the displacement, velocity and acceleration of the tower top and the near-field wake coordinates of the front wind rotor, the real-time inflow velocity of each blade section of the rear wind rotor under the influence of the tower movement is obtained; specifically, it includes:

[0074] According to the displacement, velocity and acceleration of the tower top, the positional relationship between the rear wind wheel hub and the tower top is obtained;

[0075] According to the positional relationship between the hub of the rear wind wheel and the top of the tower, the spatial coordinates, movement speed and acceleration of each blade section of the rear wind wheel under the influence of the tower movement are calculated;

[0076] Based on the spatial coordinates, motion velocity and acceleration of each blade section of the rear wind rotor, as well as the near-field wake coordinates of the front wind rotor, the actual inflow velocity of each blade section of the rear wind rotor is obtained by interpolation;

[0077] Among them, the spatial coordinates, motion speed and acceleration of each blade section of the rear wind wheel are calculated by the following formula:

[0078] Z 后 [x,y,z]=T mat Z 初始 [X,Y,X],

[0079]

[0080] In the formula, Z 后 [x, y, z] represents the spatial coordinates of any section of the rear wind turbine blade in the x, y, and z directions after considering the deformation of the tower top. 初始 [X,Y,X] represents the spatial coordinates of any section of the rear wind turbine blade before considering the deformation of the tower top; V 后 [x,y,z] represents the translational velocity of the rear wind wheel blade section in the x, y, and z directions, V 塔 [x,y,z] represents the translational velocity of the top of the tower in the x, y, and z directions, V 塔 [Rx, Ry, Rz] represents the rotation speed of the top of the tower in the x, y, and z directions; A represents the position vector of the rear wind rotor blade section and the top of the tower;后 [x, y, z] represents the translational acceleration of the rear wind wheel blade section in the x, y, and z directions, A 塔 [x, y, z] represents the translational acceleration of the top of the tower in the x, y, and z directions, A 塔 [Rx, Ry, Rz] represents the rotational acceleration of the top of the tower in the x, y, and z directions; T mat represents the motion conversion matrix, θ1, θ2 and θ3 represent the rotation angles of the top of the tower in the x, y and z directions respectively, and s is an intermediate variable with no real meaning.

[0081] Step S24: Calculate the aerodynamic force of the rear wind rotor blades and solve the structural dynamic response of the wind rotor blades according to the actual inflow velocity of each blade section of the rear wind rotor, the target pitch angle of the blades of the rear wind rotor, and the target generator torque of the rear wind rotor, and obtain the load vector at the top position of the rear wind rotor tower; specifically including:

[0082] Based on the actual inflow velocity of each blade section of the rear wind rotor and the target pitch angle of the rear wind rotor blade output by the simulation model control system, the aerodynamic force of the rear wind rotor blade is calculated using the generalized dynamic inflow theory or the blade element momentum theory;

[0083] Based on the aerodynamic force of the rear wind rotor blade, the modal method is used to calculate the structural mechanical response of the rear wind rotor blade.

[0084] The Kane method is used to establish a multi-body system dynamics model of the rear wind wheel; wherein the multi-body system dynamics model of the rear wind wheel is a multi-body system dynamics model of the rear wind wheel hub, the rear wind wheel blades and the nacelle;

[0085] Based on the target generator torque of the rear wind wheel output by the simulation model control system, the rear wind wheel load vector acting on the top of the tower is calculated using the rear wind wheel multi-body system dynamics model, that is, the load vector at the top of the rear wind wheel tower is obtained;

[0086] Wherein, the multi-body system dynamics model of the rear wind wheel is a multi-body system dynamics model of the rear wind wheel hub, the rear wind wheel blades and the nacelle;

[0087] It should be noted that before the rear wind rotor load calculated in step S24 is used to solve the dynamic response of the front wind rotor blades and tower structure in step S22 in the next round of repeated calculation, the rear wind rotor load needs to be subtracted from the nacelle inertia force and nacelle gravity.

[0088] In some specific embodiments, in the process of repeatedly iteratively calculating the blade aerodynamic forces and structural loads of the front and rear wind wheels to realize time domain simulation, the embodiment of the present application adopts the fourth-order Adams-Bashforth-Moulton prediction-correction time integration method, and the structural dynamics solution of the front and rear wind wheels is performed once in the prediction step and the correction step to ensure the stability of the entire simulation time advancement.

[0089] The present application also provides a dual-wind-rotor wind turbine simulation computing system. Figure 3 As shown, including:

[0090] A calculation control unit, used to calculate the target pitch angles and target generator torques corresponding to the front and rear wind rotors based on the basic parameters of the front and rear wind rotors, and generate control instructions to control the simulation model of the dual-wind-rotor wind turbine set;

[0091] A calculation unit, used to calculate the blade aerodynamic force and structural load of the front and rear wind rotors in the simulation model of the double-wind-rotor wind turbine according to the target pitch angles and target generator torques corresponding to the front and rear wind rotors;

[0092] A parameter updating unit is used to obtain the blade aerodynamic torque of the front and rear wind wheels based on the blade aerodynamic force of the front and rear wind wheels, and calculate the current rotation speed of the front and rear wind wheels in the basic parameters of the front and rear wind wheels according to the blade aerodynamic torque of the front and rear wind wheels, the target generator torque, the moment of inertia and the moment of inertia of the transmission system, so as to update the basic parameters of the front and rear wind wheels;

[0093] The repeated simulation unit is used to repeatedly calculate the aerodynamic force and structural load of the front and rear wind rotors until the simulation time reaches the set value, thus completing the time domain simulation of the fully coupled dynamics of the dual-rotor wind turbine set;

[0094] Among them, the dual-rotor wind turbine simulation model also includes a front wind wheel, a rear wind wheel, a nacelle and a tower. The front wind wheel and the rear wind wheel are respectively connected to the two ends of the nacelle, and the nacelle is placed on the top of the tower; the diameter of the front wind wheel is larger than that of the rear wind wheel.

[0095] Regarding the system in the above embodiment, the specific manner in which each unit module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0096] Based on the same inventive concept, the present application also provides an electronic device, whose structure is as follows: Figure 4 As shown, it includes a memory, a processor and a computer program stored in the memory, and the processor executes the computer program or instructions to implement the aforementioned dual-wind-rotor wind turbine simulation calculation method.

[0097] Based on the same inventive concept, an embodiment of the present application also provides a computer storage medium, in which a computer program or instruction is stored. When the computer program or instruction is executed by a processor, the steps of the aforementioned dual-rotor wind turbine simulation calculation method are implemented.

[0098] Based on the same inventive concept, an embodiment of the present application also provides a computer program product, including a computer program or instructions, which, when executed by a processor, implements the steps of the aforementioned dual-rotor wind turbine simulation calculation method.

[0099] In order to verify that the calculation method proposed in this application can take into account the full coupling effect of the tandem twin-rotor wind turbine, the NREL 5MW wind turbine is installed on one side of the IEA15MW nacelle to form a 20MW tandem twin-rotor wind turbine. Figure 5 and Figure 6 The thrust and power of the front and rear wind wheels at different wind speeds calculated by the technical solution of the present application (i.e., the present method) are compared with those calculated by the traditional CFD method, and as shown in FIG. Figure 7 The figure shows a comparison of the calculation time of the two methods. It can be found that the results obtained by the two calculation methods are very consistent, which shows that the technical solution of the present application can effectively consider the wake interference effect of the front and rear wind wheels, and the method has higher calculation efficiency and the time required is only 1 / 169 of the CFD method.

[0100] Figure 8 , Fig. 9 and Fig.10 The calculation results of the rotor thrust, power generation and hub wind speed of the front and rear rotors of the double-rotor wind turbine composed of IEA15MW and NREL5MW rotors under turbulent wind conditions calculated based on this method are compared; from the results, it can be seen that there are significant differences between the results of the front and rear rotors, and the rotor thrust of the rear rotor oscillates greatly at the initial moment, indicating that the aerodynamic performance of the rear rotor is affected by the front rotor. The wind speed at the hub of the rear rotor is greater than that of the front rotor, which is consistent with the actual measured results of the wind farm. In addition, although the wind speed at the hub is higher, the wind energy utilization efficiency is lower, because the wind speed in the rest of the rotor plane is smaller, which is consistent with the actual situation. The calculation results are very reasonable.

[0101] Fig.11 and Fig.12 The results of the dynamic response of the double-rotor wind turbine composed of IEA15MW and NREL5MW wind turbines and the IEA15MW single-rotor wind turbine under turbulent wind conditions are compared (i.e., the tower top load comparison). It can be seen that the horizontal thrust of the double-rotor wind turbine is larger, and the increase is close to the thrust of the NREL 5MW rear wind turbine. The vertical load on the top of the tower is larger, and the increase is close to the gravity of the NREL5MW wind turbine. The calculation results are very logical.

[0102] Finally, it should be noted that the above is only a preferred embodiment of the present application and is not intended to limit the present application. Although the present application has been described in detail with reference to the aforementioned embodiments, it is still possible for those skilled in the art to modify the technical solutions described in the aforementioned embodiments or to make equivalent substitutions for some of the technical features therein. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present application shall be included in the protection scope of the present application.

Claims

1. A simulation calculation method for a double-wind-rotor wind turbine, characterized in that: include: S1: Based on the basic parameters of the front and rear wind rotors, the target blade pitch angles of the front and rear wind rotors and the target generator torques of the front and rear wind rotors are calculated, and control instructions are generated to control the simulation model of the dual-rotor wind turbine set; S2: Calculate the blade aerodynamic force and structural load of the front and rear wind rotors in the simulation model of the dual-rotor wind turbine according to the target blade pitch angles of the front and rear wind rotors and the target generator torques of the front and rear wind rotors; S3: obtaining the blade aerodynamic torque of the front and rear wind wheels based on the blade aerodynamic forces of the front and rear wind wheels, and calculating the current rotation speeds of the front and rear wind wheels in the basic parameters of the front and rear wind wheels according to the blade aerodynamic torque of the front and rear wind wheels, the target generator torque, the moment of inertia and the moment of inertia of the transmission system for updating the basic parameters of the front and rear wind wheels; Repeat the operations of step S1 to step S3 to iteratively calculate the blade aerodynamic force and structural load of the front and rear wind rotors until the simulation time reaches a preset value, completing the time domain simulation of the fully coupled dynamics of the dual-wind-rotor wind turbine set; Among them, the dual-rotor wind turbine simulation model includes a front wind wheel, a rear wind wheel, a nacelle and a tower. The front wind wheel and the rear wind wheel are respectively connected to the two ends of the nacelle, and the nacelle is placed on the top of the tower; the diameter of the front wind wheel is larger than that of the rear wind wheel.

2. A simulation calculation method for a dual-wind-rotor wind turbine according to claim 1, characterized in that: Before calculating the target blade pitch angles of the front and rear wind rotors and the target generator torques of the front and rear wind rotors for the first time, setting the design parameters of the dual-wind-rotor wind turbine simulation model; The design parameters of the dual-rotor wind turbine simulation model include: aerodynamic calculation parameters and structural dynamic calculation parameters of the wind rotor blades; The wind rotor blade aerodynamic calculation parameters include: chord length, twist angle and airfoil lift and drag coefficient at different radii of the front and rear wind rotor blades; The structural dynamics calculation parameters include: modal vibration shapes of each order mode of the blade and the corresponding damping ratio, cross-sectional mass at different radii of the blade, flapping stiffness at different radii of the blade, swing stiffness at different radii of the blade, hub position coordinates and hub mass of the front and rear wind wheels, moment of inertia of the transmission system, cabin mass, cabin moment of inertia, tower mass, distribution information of tower stiffness along the tower height, modal vibration shape of the tower, and damping ratio of the tower.

3. The simulation calculation method of a dual-wind-rotor wind turbine set according to claim 1, characterized in that: Based on the basic parameters of the front and rear wind rotors, the target blade pitch angles of the front and rear wind rotors and the target generator torques of the front and rear wind rotors are calculated, and control instructions are generated to control the simulation model of the dual-rotor wind turbine, including: The basic parameters of the front and rear wind wheels include: the current blade root loads of the front and rear wind wheels, the current blade deformations of the front and rear wind wheels, the current rotation speeds of the front and rear wind wheels, and the current generator torques of the front and rear wind wheels; The basic parameters of the front and rear wind rotors are input into the simulation model control system, and based on the state space of the entire dual-rotor wind turbine, the target pitch angle and target generator torque corresponding to the front and rear wind rotors are calculated by the algorithm; The simulation model control system generates a control instruction according to the target pitch angle and the target generator torque to control the dual-rotor wind turbine simulation model; Among them, the control closed loop of the front wind wheel needs to consider the influence factors of the rear wind wheel and the tower, and the control closed loop of the rear wind wheel needs to consider the influence factors of the front wind wheel and the tower; The influencing factors of the rear wind wheel include the rear wind wheel speed, the rear wind wheel power generation and the rear wind wheel thrust; the influencing factors of the front wind wheel include the front wind wheel speed, the front wind wheel power generation and the front wind wheel thrust; the influencing factors of the tower include the tower motion state and the tower load.

4. The simulation calculation method for a dual-wind-rotor wind turbine according to claim 1, characterized in that: According to the target pitch angles and target generator torques corresponding to the front and rear wind rotors, the blade aerodynamic forces and structural loads of the front and rear wind rotors in the simulation model of the dual-rotor wind turbine are calculated, including: According to the defined front wind rotor wind field environment data, the blade aerodynamic force of the front wind rotor and the near-field wake coordinates of the front wind rotor are calculated respectively; Considering the influence of the rear wind rotor, the dynamic response of the front wind rotor blades and tower structure is calculated to obtain the displacement, velocity and acceleration of the tower top; Based on the displacement and velocity of the tower top and the near-field wake coordinates of the front wind rotor, the real-time inflow velocity of each blade section of the rear wind rotor under the influence of the tower movement is obtained; According to the actual inflow velocity of each blade section of the rear wind rotor, the target pitch angle of the blades of the rear wind rotor and the target generator torque of the rear wind rotor, the blade aerodynamic force of the rear wind rotor is calculated and the structural dynamic response of the rear wind rotor blades is solved to obtain the load vector at the top position of the rear wind rotor tower.

5. A simulation calculation method for a dual-wind-rotor wind turbine according to claim 3 or 4, characterized in that: According to the defined front wind rotor wind field environment data, the front wind rotor blade aerodynamic force and the front wind rotor near-field wake coordinates are calculated respectively, including: The front wind rotor wind field environment data is the data of changes in the three directional components of wind speed in the plane of the front wind rotor and the tower over time; The type of the front wind wheel wind field is turbulent wind or steady wind or other wind that satisfies the wind shear law; Among them, according to the defined front wind rotor wind field environment data, the free vortex wake method and blade element momentum theory are used to calculate the aerodynamic force of the front wind rotor blade, including: According to the front wind rotor wind field environment data, a real-time inflow velocity at each blade cross section of the front wind rotor is obtained, and based on the real-time inflow velocity and the current rotation speed of the front wind rotor, an induced velocity at each blade cross section control point of the front wind rotor is obtained by solving the free vortex wake method; Determine the velocity triangle at each blade section of the front wind wheel according to the initial rotation speed of the front wind wheel and the induced speed, and calculate the actual angle of attack of each blade section of the front wind wheel based on the velocity triangle; According to the lift and drag coefficient of the airfoil of the front wind rotor, the lift coefficient and drag coefficient at the actual angle of attack of each blade section are retrieved, and the aerodynamic force of the corresponding blade section in the front wind rotor is calculated using the blade element momentum theory; The aerodynamic force of each blade cross section of the front wind wheel is superimposed to obtain the blade aerodynamic force of the entire front wind wheel.

6. A simulation calculation method for a dual-wind-rotor wind turbine according to claim 5, characterized in that: According to the defined front wind rotor wind field environment data, calculate the front wind rotor near-field wake coordinates, including: According to the front wind wheel wind field environment data, the real-time inflow velocity at each blade section of the front wind wheel is obtained, and the induced velocity at any wake behind the front wind wheel plane is obtained by solving the free vortex wake method; The induced velocity field at any wake behind the front wind wheel plane is superimposed with the environmental space field of the real-time inflow velocity to obtain the three-dimensional velocity space coordinates of the near-field wake of the front wind wheel, that is, to obtain the near-field wake coordinates of the front wind wheel; Wherein, the near-field wake coordinates of the front wind rotor are the three-dimensional velocity space coordinates of the near-field wake of the front wind rotor; The three-dimensional velocity space coordinates of the front wind wheel near-field wake are: In the formula, represents the three-dimensional velocity space coordinates of the near-field wake of the front wind rotor, represents the induced velocity at any wake behind the front wind rotor plane; Indicates the real-time inflow velocity at each blade section of the front wind wheel.

7. The simulation calculation method for a dual-wind-rotor wind turbine according to claim 3, characterized in that: Considering the influence of the rear wind rotor, the dynamic response of the front wind rotor blades and tower structure is calculated to obtain the displacement, velocity and acceleration of the tower top, including: The Kane method is used to establish a front wind rotor multi-body system dynamics model; wherein the front wind rotor multi-body system dynamics model is a multi-body system dynamics model of the front wind rotor blades, the front wind rotor hub, the tower and the nacelle; Considering the gravity and motion effects of the remaining structures except the wind rotor blades, as well as the load effect of the rear wind rotor, the dynamic response of the front wind rotor blades and the tower structure is calculated using the front wind rotor multi-body system dynamics model to obtain the displacement, velocity and acceleration of the tower top. Among them, except for the wind rotor blades, the remaining structures include: transmission system, nacelle, wind rotor hub and tower; the load action of the rear wind rotor is the comprehensive action of the rear wind rotor aerodynamic force, the rear wind rotor inertia force and the rear wind rotor elastic force, and does not include the action of the nacelle gravity and the nacelle inertia force; When calculating the dynamic response of the front wind rotor blades and tower structure, the modal method is used for the structural deformation of the wind rotor blades and tower.

8. A simulation calculation method for a dual-wind-rotor wind turbine according to claim 3 or 7, characterized in that: Based on the displacement, velocity and acceleration of the tower top and the near-field wake coordinates of the front wind rotor, the actual inflow velocity of each blade section of the rear wind rotor under the influence of the tower motion is obtained, including: According to the displacement, velocity and acceleration of the tower top, the positional relationship between the rear wind wheel hub and the tower top is obtained; According to the positional relationship between the hub of the rear wind wheel and the top of the tower, the spatial coordinates, movement speed and acceleration of each blade section of the rear wind wheel under the influence of the tower movement are calculated; Based on the spatial coordinates, motion velocity and acceleration of each blade section of the rear wind rotor, as well as the near-field wake coordinates of the front wind rotor, the actual inflow velocity of each blade section of the rear wind rotor is obtained by interpolation; Among them, the spatial coordinates, motion speed and acceleration of each blade section of the rear wind wheel are calculated by the following formula: In the formula, Z 后 [x, y, z] represents the spatial coordinates of any section of the rear wind turbine blade in the x, y, and z directions after considering the deformation of the tower top. 初始 [X,Y,X] represents the spatial coordinates of any section of the rear wind turbine blade before considering the deformation of the tower top; V 后 [x,y,z] represents the translational velocity of the rear wind wheel blade section in the x, y, and z directions, V 塔 [x,y,z] represents the translational velocity of the top of the tower in the x, y, and z directions, V 塔 [Rx, Ry, Rz] represents the rotation speed of the tower top in the x, y, and z directions; P represents the position vector of the rear wind wheel blade section and the tower top; A 后 [x, y, z] represents the translational acceleration of the rear wind wheel blade section in the x, y, and z directions, A 塔 [x, y, z] represents the translational acceleration of the top of the tower in the x, y, and z directions, A 塔 [Rx, Ry, Rz] represents the rotational acceleration of the top of the tower in the x, y, and z directions; T mat represents the motion conversion matrix, θ1, θ2 and θ3 represent the rotation angles of the top of the tower in the x, y and z directions respectively, and s is an intermediate variable with no real meaning.

9. The simulation calculation method for a dual-wind-rotor wind turbine according to claim 3, characterized in that: According to the actual inflow velocity of each blade section of the rear wind rotor, the target blade pitch angle of the rear wind rotor and the target generator torque of the rear wind rotor, the aerodynamic force of the rear wind rotor blade is calculated and the structural dynamic response of the wind rotor blade is solved to obtain the load vector at the top of the rear wind rotor tower, including: Based on the actual inflow velocity of each blade section of the rear wind rotor and the target pitch angle of the blades of the rear wind rotor, the aerodynamic force of the rear wind rotor blades is calculated using the generalized dynamic inflow theory or the blade element momentum theory; Based on the aerodynamic force of the rear wind rotor blade, the modal method is used to calculate the structural mechanical response of the rear wind rotor blade. The Kane method is used to establish a multi-body system dynamics model of the rear wind wheel; wherein the multi-body system dynamics model of the rear wind wheel is a multi-body system dynamics model of the rear wind wheel hub, the rear wind wheel blades and the nacelle; Based on the target generator torque of the rear wind wheel and the mechanical response of the rear wind wheel blade structure, the rear wind wheel load vector acting on the top of the tower is calculated using the rear wind wheel multi-body system dynamics model, that is, the load vector at the top of the rear wind wheel tower is obtained.

10. A dual-wind-rotor wind turbine simulation calculation system, characterized in that: include: A calculation control unit, used to calculate the target pitch angles and target generator torques corresponding to the front and rear wind rotors based on the basic parameters of the front and rear wind rotors, and generate control instructions to control the simulation model of the dual-wind-rotor wind turbine set; A calculation unit, used to calculate the blade aerodynamic force and structural load of the front and rear wind rotors in the simulation model of the double-wind-rotor wind turbine according to the target pitch angles and target generator torques corresponding to the front and rear wind rotors; A parameter updating unit is used to obtain the blade aerodynamic torque of the front and rear wind wheels based on the blade aerodynamic force of the front and rear wind wheels, and calculate the current rotation speed of the front and rear wind wheels in the basic parameters of the front and rear wind wheels according to the blade aerodynamic torque of the front and rear wind wheels, the target generator torque, the moment of inertia and the moment of inertia of the transmission system, so as to update the basic parameters of the front and rear wind wheels; The repeated simulation unit is used to repeatedly calculate the aerodynamic force and structural load of the front and rear wind rotors until the simulation time reaches the set value, thus completing the time domain simulation of the fully coupled dynamics of the dual-rotor wind turbine set; Among them, the dual-rotor wind turbine simulation model also includes a front wind wheel, a rear wind wheel, a nacelle and a tower. The front wind wheel and the rear wind wheel are respectively connected to the two ends of the nacelle, and the nacelle is placed on the top of the tower; the diameter of the front wind wheel is larger than that of the rear wind wheel.

Images