Power design and curve acquisition method for bladeless wind driven generator

Abstract

The invention relates to the technical field of wind power generation, in particular to a blade-free wind driven generator power design and curve acquisition method, which comprises the steps of theoretical simulation analysis and experimental verification: firstly, selecting and designing a cylindrical linear generator, simulating and verifying the frequency-displacement-induced voltage relationship of the cylindrical linear generator, then simulating to obtain the lift-drag coefficient and vibration parameter of a rod body, and calculating the curve of the blade-free wind driven generator. Calculating a theoretical power curve; an actual measurement curve is obtained through one-year outdoor testing and data processing, and the average value of the two curves is taken to form a target power curve. According to the method, a theoretical power curve is obtained by using a whole process of theoretical simulation, experimental verification, data correction and multi-scene verification, building a model and calculating parameters by virtue of professional software and combining actually measured resistance; a target curve is obtained through Bin processing, air density correction and interpolation data supplementation and is verified, meanwhile, the structure of the cylindrical linear generator is optimized, adaptability, stability, electromagnetic performance and supporting reliability are considered, and data support and operation guarantee are provided for a unit.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of wind power generation, and particularly relates to a power design and curve acquisition method for a bladeless wind turbine. BACKGROUND

[0002] The bladeless wind turbine has significant advantages in safety, maintenance convenience and flexible application scenarios over the traditional wind turbine with blades due to its structural characteristics of no rotating blades, and has become a research hotspot in the field of wind power generation in recent years. However, the accurate acquisition of the power curve and the adaptability design of the linear generator structure are still key technical bottlenecks restricting the large-scale application of this type of unit.

[0003] In terms of power curve acquisition, the existing technology mainly relies on single theoretical simulation or short-term outdoor experiment: when the power curve is calculated by establishing a model through simulation software, the dynamic changes of air density with temperature, humidity and air pressure in the actual wind field, and the coupling relationship between the vibration characteristics of the rod and the electromagnetic conversion characteristics of the generator are often ignored, resulting in a significant deviation between the theoretical curve and the actual operating power of the unit; when the power data is obtained only by relying on short-term experiments, the wind speed and direction in the natural wind field fluctuate greatly, the experimental data are highly discrete, and the wind speed intervals with insufficient data are not effectively completed, so it is difficult to form a power curve that can reflect the long-term stable operation state of the unit, and thus it is impossible to provide reliable basis for the power design, grid connection scheduling and operation optimization of the bladeless wind turbine.

[0004] In terms of linear generator structure design, the existing linear generators adapted to the rod of the bladeless wind turbine generally have the problems of poor structural adaptability, insufficient operation stability and low electromagnetic conversion efficiency: some generators have mismatched sizes of the mover and the stator and the internal space of the rod, or unreasonable spacing design between the mover and the stator, which leads to collision or excessive magnetic field loss during operation; some generators have large sliding friction resistance of the mover, and are prone to mechanical wear during long-term operation, affecting the service life of the unit; and some generators have improper winding method and support material selection of the stator coil, resulting in poor superposition effect of induced electromotive force and unstable output voltage, which is difficult to meet the demand of efficient power generation of the bladeless wind turbine.

[0005] In summary, there is an urgent need for a technical solution that can accurately acquire the power curve of the bladeless wind turbine and optimize the structure design of the linear generator, so as to solve the problems of low reliability of the power curve, insufficient adaptability and operation performance of the generator structure in the existing technology, and promote the technical maturity and industrial application of the bladeless wind turbine. SUMMARY

[0006] The present application aims at making up for the shortage of the prior art, and provides a bladeless wind turbine power design and curve acquisition method.

[0007] The present application aims at making up for the shortage of the prior art, and provides a bladeless wind turbine power design and curve acquisition method. Step S1: linear generator selection and structure design, according to the size of the bladeless wind turbine rod, a cylindrical linear generator capable of being installed in the rod is selected, the linear generator comprises inner mover magnets and outer stator coils, a gap is arranged between the mover magnets and the stator coils, the size of the gap is half of the maximum displacement of the mover in the working process of the linear generator; a coil interval support layer is arranged between the stator coils, the coil interval support layer is made of non-conductive and high-strength material; the mover magnets are made of magnetic conductive material, and small rolling balls are embedded on the interval layer outside the mover magnets; the height of the mover magnets, the height of the single coil, the interval height between the mover magnets and the interval height between the stator coils are equal; the mover magnets are arranged in the polarity arrangement mode of NS|SN|NS|SN; the winding turns of the stator coils are the same and the coils are tightly wound, a plurality of stator coils are first divided into odd coil groups and even coil groups according to the positive and negative poles of the induced voltage, the coils in the odd coil group are connected in series, the coils in the even coil group are connected in series, and then the odd coil group and the even coil group are connected in reverse series; Step S2: Simulation and experimental verification of the relationship between the linear generator frequency, displacement and induced voltage, a three-dimensional model of the linear generator designed in step S1 is established in the simulation software, a sinusoidal motion curve is applied to the three-dimensional model, the peak value of the sinusoidal motion curve is the maximum displacement of the mover, and the period of the sinusoidal motion curve is the frequency of the mover motion, by changing the peak value and period of the sinusoidal motion curve, a plurality of groups of induced voltage data corresponding to different frequencies and displacements are obtained, and a frequency-displacement-induced voltage relationship model is established; according to the structure of the linear generator designed in step S1, an entity prototype is manufactured, the same frequency and displacement motion is applied to the entity prototype according to the parameters of the sinusoidal motion curve applied in the simulation process, the induced voltage generated by the entity prototype is tested, the induced voltage data obtained by testing is compared with the induced voltage data obtained by simulation, if the error between the two exceeds the preset threshold, the parameters in the simulation software are adjusted until the error between the induced voltage data obtained by simulation and the induced voltage data obtained by testing is within the preset threshold range, and the verification of the frequency-displacement-induced voltage relationship model is completed; Step S3: Simulation of the lift and drag coefficients of the bladeless wind turbine, a flow field model and a rod model of the bladeless wind turbine are established in the simulation software, the inlet of the flow field model is set as a boundary condition of adjustable wind speed, different wind loads per meter are applied to the inlet of the flow field model, the lift coefficient and the drag coefficient of the rod section of the bladeless wind turbine under different wind speeds are extracted through simulation calculation, and the maximum deformation of the rod section under different wind speeds is obtained; Step S4: Obtain the vibration frequency and displacement of the rod of the bladeless wind turbine, respectively perform Fourier transform on the lift data and the drag data under different wind speeds obtained in step S3 to obtain the frequency spectrum of the lift signal and the drag signal, extract the frequency corresponding to the maximum amplitude from the frequency spectrum as the vibration frequency of the rod under the wind speed, and take the maximum deformation of the rod section obtained in step S3 as the vibration displacement of the rod under the wind speed; the calculation formula of the Fourier transform is: wherein, is the lift or drag signal in the time domain, is the lift or drag signal spectrum in the frequency domain, is the frequency, in Hz, is the imaginary unit, is the time, in s; Step S5: Obtain the theoretical power curve of the bladeless wind turbine, substitute the rod vibration frequency and vibration displacement under different wind speeds obtained in step S4 into the frequency-displacement-induced voltage relationship model verified in step S2 to obtain the induced voltage of the linear generator under different wind speeds ; according to the resistance of the linear generator circuit itself , the power calculation formula is used to calculate the power generation of the bladeless wind turbine under different wind speeds , the power calculation formula is: The wind speed is taken as the horizontal coordinate and the power generation is taken as the vertical coordinate to draw the theoretical power curve of the bladeless wind turbine. Step S6: outdoor test of the bladeless wind turbine, the bladeless wind turbine equipped with the linear generator designed in step S1 is moved to the natural wind field for outdoor test for one year; a simple wind tower or laser radar is installed near the generator rod for collecting wind speed and wind direction data during the test; temperature and humidity sensors and air pressure sensors are installed on the test site for collecting temperature, humidity and air pressure data during the test; voltage and current transformers are connected to the output end of the linear generator, and the output signals of the voltage and current transformers are connected to the power transducer to collect power generation data at different times through the power transducer; Step S7: processing of outdoor test data, the wind speed-power scatter point data collected in step S6 is processed by Bin, the method of Bin processing is: the wind speed is divided into multiple wind speed intervals according to the interval of 1m / s, the range of each wind speed interval is , wherein is an integer wind speed value, unit: m / s; the average value of all power generation data in each wind speed interval is calculated as the actual measured power corresponding to the integer wind speed value ; according to the collected temperature , humidity and air pressure , the actual air density during the test is calculated, and the air density calculation formula is: , wherein is the specific gas constant of dry air, the value is 287J / (kg・K), is the absolute temperature, unit: K, , is the Celsius temperature, unit: ℃; according to the actual air density and the standard air density , the measured power after Bin processing is standardized and corrected, and the power correction formula is: , wherein is the measured power under the actual air density, The measured power curve of the bladeless wind turbine is plotted with the wind speed as the horizontal coordinate and the measured power under the standard air density as the vertical coordinate. Step S8: Forming the target power curve of the bladeless wind turbine, comparing the theoretical power curve obtained in step S5 with the measured power curve obtained in step S7, and calculating the average value of the theoretical power and the measured power under the same wind speed as the target power under the same wind speed. The target power calculation formula is as follows: The target power curve of the bladeless wind turbine is plotted with the wind speed as the horizontal coordinate and the target power as the vertical coordinate.

[0008] Further, in step S1, the material of the coil spacing support layer is selected from any one of epoxy resin, glass fiber reinforced plastic, or ceramic; the mover magnet is made of neodymium iron boron permanent magnet material; and the small ball is made of stainless steel material, and the diameter of the small ball is 0.5-2mm.

[0009] Further, in step S2, the simulation software is any one of ANSYS Maxwell or COMSOL Multiphysics; the preset threshold is 5%, that is, when the relative error between the simulation obtained induced voltage data and the test obtained induced voltage data is less than or equal to 5%, it is determined that the frequency-displacement-induced voltage relationship model is verified; and the calculation formula of the relative error is as follows: wherein, is the relative error, is the simulation obtained induced voltage, is the test obtained induced voltage.

[0010] Further, in step S3, the simulation software is any one of ANSYS Fluent or STAR-CCM+; when establishing the flow field model, the calculation domain of the flow field is a cylindrical region with the rod body as the center, the radius is 10-20 times the diameter of the rod body, and the height is 1.5-2 times the height of the rod body; the inlet boundary condition of the flow field model is set as a velocity inlet, the outlet boundary condition is set as a pressure outlet, and the wall boundary condition is set as a no-slip wall; and the simulation calculation adopts the RNG k-ε turbulence model or the SST k-ω turbulence model.

[0011] Further, in step S4, before performing Fourier transform on the lift data and the drag data, the data is preprocessed, and the preprocessing steps include: removing outliers in the data, using The criteria for identifying outliers are as follows: when a data point deviates from the mean by more than 3 times the standard deviation, the data point is considered an outlier and is removed. The data after outlier removal is then smoothed using a moving average method with a moving window size of 5-10 data points. The formula for calculating the criterion is: in, For a single lift or drag data point, This is the average of the lift or drag data. This represents the standard deviation of the lift or drag data; the formula for calculating the moving average method is: in, For the first Smoothed values ​​for each data point To resize the window, These are the original data points.

[0012] Furthermore, in step S5, the inherent resistance of the linear generator line... The measurement method, obtained through actual measurement, is as follows: A DC resistance tester is used, with its two test clips connected to the two ends of the linear generator coil. A constant DC voltage is applied, and the DC current in the circuit is measured according to Ohm's law. Calculate the resistance of the line itself, where For the applied DC voltage, The measured DC current is used; the measurement should be performed at least three times, and the average value of the multiple measurements is taken as the final line resistance. .

[0013] Furthermore, in step S6, the height of the simplified wind measuring tower is the same as the height of the bladeless wind turbine pole, and the sampling frequency of the wind speed sensor and wind direction sensor installed on the wind measuring tower is 1Hz; the measurement accuracy of the temperature and humidity sensor is ±0.5℃ for temperature and ±3%RH for humidity, and the measurement accuracy of the barometric pressure sensor is ±0.1hPa; the transformation ratio of the voltage transformer is 100:1-500:1, the transformation ratio of the current transformer is 5:1-20:1, the measurement accuracy of the power transmitter is 0.5 class, and the sampling frequency of the power data is 1Hz.

[0014] Furthermore, in step S7, during the Binning process, if the number of power generation data points within a certain wind speed range is less than 10, the power data corresponding to that wind speed range is removed, or the power data for that wind speed range is supplemented using linear interpolation; the formula for linear interpolation is: in, for the wind speed to be supplemented corresponding power generation, for the wind speed corresponding power generation, for the wind speed corresponding power generation.

[0015] Further, the method further comprises a step S9 of verifying and optimizing the target power curve, selecting a site different from the test wind field in step S6 in the natural wind field, performing a verification test for 3 months, obtaining a measured power curve of the verification test according to the method of step S6 and step S7; comparing the measured power curve of the verification test with the target power curve obtained in step S8, calculating the relative deviation of the verification measured power and the target power generation at the same wind speed, if the maximum value of the relative deviation is less than or equal to 8%, it is determined that the target power curve meets the design requirements; if the maximum value of the relative deviation is greater than 8%, return to step S2 to re-adjust the parameters of the frequency-displacement-induction voltage relationship model, or return to step S3 to re-optimize the flow field simulation parameters, until the target power curve passes the verification test; the calculation formula of the relative deviation is: wherein, is the relative deviation, is the measured power generation of the verification test, is the target power generation.

[0016] Compared with the prior art, the bladeless wind turbine power design and curve acquisition method has the following beneficial effects: I. The present application greatly improves the accuracy and reliability of the power curve of the bladeless wind turbine through the whole process design of theoretical simulation analysis-experimental verification-data correction-multi-scene verification. The linear generator frequency-displacement-induced voltage relationship model is established by professional software such as ANSYS Maxwell / COMSOL Multiphysics, the lift and drag coefficients and vibration parameters of the rod body are obtained by combining ANSYS Fluent / STAR-CCM+ flow field simulation, and the theoretical power curve is calculated based on the measured line resistance to ensure that the theoretical model is highly matched with the actual structure and fluid characteristics of the unit. At the same time, bin processing (wind speed 1 m / s interval division interval) and air density standardization correction (based on measured temperature, humidity and air pressure to calculate the actual air density and correct the power), aiming at the insufficient wind speed interval, linear interpolation is used to supplement, which effectively avoids the interference of natural wind field environment fluctuation and data dispersion on the measured curve. Finally, the target power curve is formed by taking the mean value of the theoretical and measured curves, and additional 3-month verification tests are carried out in different wind fields to ensure that the target curve is within 8% relative deviation threshold, which can directly provide accurate data support for the power design, type selection, grid connection scheduling and operation optimization of the bladeless wind turbine, and solve the problem of large deviation between the curve and the actual operation caused by the traditional method relying only on simulation or short-term experiment.

[0017] II. The present application optimizes the structure of cylindrical linear generator according to the internal space characteristics of the rod body of the bladeless wind turbine, and ensures the operation performance and power generation effect of the unit from multiple dimensions. In terms of structural adaptability, the rod body size is selected, and the height of the single magnet of the mover, the height of the single coil and the interval height between them are equal, which ensures that the structure is compact and can be completely installed in the internal space of the rod body, avoiding space waste. In terms of running stability, the interval between the mover and the stator is half of the maximum displacement of the mover, which not only prevents the collision between the mover and the stator during movement, but also reduces the magnetic field loss caused by excessive magnetic gap. At the same time, small stainless steel balls (0.5-2mm) are embedded in the interval layer outside the mover, which reduces the sliding friction resistance of the mover and reduces mechanical wear. In terms of electromagnetic performance optimization, the mover uses neodymium iron boron permanent magnet material (high magnetic energy product) and is arranged in NS|SN|NS|SN polarity, which enhances the magnetic field strength. The stator coil is tightly wound and has the same number of turns. Through the series connection of odd coil groups, the series connection of even coil groups and the reverse series connection of the two groups, the induced electromotive force is effectively superimposed, and the output voltage is improved. In terms of coil support reliability, the coil interval support layer is made of non-conductive and high-strength materials such as epoxy resin, glass fiber reinforced plastic or ceramic, which can not only isolate the current interference between the coils, but also can withstand the mechanical stress caused by the vibration of the rod body, avoiding the deformation and damage of the coil. The above structure design makes the linear generator stably output induced voltage during the vibration of the rod body, and takes into account the mechanical stability and electromagnetic conversion efficiency, laying a foundation for the long-term reliable operation and high-efficiency power generation of the bladeless wind turbine. BRIEF DESCRIPTION OF DRAWINGS

[0018] Figure 1 is a power curve acquisition flowchart of the bladeless wind turbine power design and curve acquisition method of the embodiment of the present application; Figure 2 is an outdoor experiment data processing and verification flowchart of the bladeless wind turbine power design and curve acquisition method of the embodiment of the present application; Figure 3 is a linear generator structure schematic diagram of the bladeless wind turbine power design and curve acquisition method of the embodiment of the present application; Figure 4 is a bladeless wind turbine structure schematic diagram of the bladeless wind turbine power design and curve acquisition method of the embodiment of the present application; Figure 5 is a theoretical power curve schematic diagram of the bladeless wind turbine power design and curve acquisition method of the embodiment of the present application. DETAILED DESCRIPTION

[0019] In order to further illustrate the technical means and effects adopted by the present application to achieve the predetermined invention purposes, the specific embodiments, structures, features and effects according to the present application are described in detail as follows in combination with the drawings and preferred embodiments.

[0020] Please refer to Figures 1 to 5 , the embodiment of the present application provides a bladeless wind turbine power design and curve acquisition method, which comprises a theoretical simulation analysis step and an experimental verification step. The bladeless wind turbine power curve is acquired through comparison between the theoretical simulation analysis step and the experimental verification step. The specific steps are as follows: Step S1: linear generator selection and structure design. According to the size of the rod body of the bladeless wind turbine, a cylindrical linear generator capable of being fitted into the rod body is selected. The linear generator comprises inner mover magnets and outer stator coils. A gap is arranged between the mover magnets and the stator coils. The size of the gap is half of the maximum displacement of the mover during the working process of the linear generator. A coil interval support layer is arranged between the stator coils. The coil interval support layer is made of non-conductive and high-strength material. The mover magnets are made of magnetic conductive material, and small rolling balls are embedded on the interval layer outside the mover magnets. The height of the single mover magnet, the height of the single coil, the interval height between the mover magnets and the interval height between the stator coils are all equal. The mover magnets are arranged in the polarity arrangement mode of NS|SN|NS|SN. The stator coils have the same number of turns and the coils are tightly wound. The multiple stator coils are first divided into odd coil groups and even coil groups according to the positive and negative poles of the induced voltage. The coils in the odd coil group are connected in series, the coils in the even coil group are connected in series, and then the odd coil group and the even coil group are connected in reverse series; Step S2: Simulation and experimental verification of the relationship between frequency, displacement and induced voltage of the linear generator. A three-dimensional model of the linear generator designed in Step S1 is established in the simulation software. A sinusoidal motion curve is applied to the three-dimensional model. The peak value of the sinusoidal motion curve represents the maximum displacement of the mover, and the period of the sinusoidal motion curve represents the frequency of the mover motion. By changing the peak value and period of the sinusoidal motion curve, multiple sets of induced voltage data corresponding to different frequencies and displacements are obtained, and a frequency-displacement-induced voltage relationship model is established. A physical prototype is fabricated based on the linear generator structure designed in Step S1. Following the parameters of the sinusoidal motion curve applied during the simulation, the same frequency and displacement are applied to the physical prototype. The induced voltage generated by the physical prototype is tested. The tested induced voltage data is compared with the simulated induced voltage data. If the error between the two exceeds a preset threshold, the parameters in the simulation software are adjusted until the error between the simulated induced voltage data and the tested induced voltage data is within the preset threshold range, thus completing the verification of the frequency-displacement-induced voltage relationship model. Step S3: Simulation of lift and drag coefficients of bladeless wind turbine. In the simulation software, establish the flow field model and the rod model of the bladeless wind turbine. Set the inlet of the flow field model to an adjustable wind speed boundary condition. Apply wind loads with different wind speeds per meter to the inlet of the flow field model. Through simulation calculation, extract the lift coefficient and drag coefficient of the rod section of the bladeless wind turbine under different wind speeds. At the same time, obtain the maximum deformation of the rod section under different wind speeds. Step S4: Obtaining the vibration frequency and displacement of the bladeless wind turbine pole. Fourier transforms are performed on the lift and drag data obtained in Step S3 at different wind speeds to obtain the frequency spectra of the lift and drag signals. The frequency corresponding to the maximum amplitude is extracted from the frequency spectra and used as the vibration frequency of the pole at that wind speed. The maximum deformation of the pole cross-section obtained in Step S3 is used as the vibration displacement of the pole at that wind speed. The Fourier transform calculation formula is: in, This is a lift or drag signal in the time domain. The frequency spectrum of the lift or drag signal in the frequency domain. Frequency, in Hz. The imaginary unit, Time, in seconds; Step S5: Obtaining the theoretical power curve of the bladeless wind turbine. Substitute the vibration frequency and displacement of the pole at different wind speeds obtained in Step S4 into the frequency-displacement-induced voltage relationship model verified in Step S2 to obtain the induced voltage of the linear generator at different wind speeds. Based on the line resistance of the linear generator itself , the power calculation formula is used to calculate the power generation of the bladeless wind turbine under different wind speeds , the power calculation formula is: The wind speed is taken as the horizontal coordinate and the power generation is taken as the vertical coordinate to draw the theoretical power curve of the bladeless wind turbine. Step S6: Outdoor test of the bladeless wind turbine equipped with the linear generator designed in step S1, move the whole machine to the natural wind field for outdoor test for one year; install a simple wind tower or laser radar near the generator rod for collecting wind speed and direction data during the test; install temperature and humidity sensors and air pressure sensors at the test site for collecting temperature, humidity and air pressure data during the test; insert voltage and current transformers into the output end of the linear generator, connect the output signals of the voltage and current transformers to the power transducer, and collect power generation data at different times through the power transducer; Step S7: processing of outdoor test data, Bin processing of the wind speed-power scatter point data collected in step S6, the method of Bin processing is: dividing the wind speed into multiple wind speed intervals according to 1m / s interval, the range of each wind speed interval is , wherein is an integer wind speed value, unit: m / s; calculate the average value of all power generation data in each wind speed interval as the corresponding measured power generation of the integer wind speed value According to the collected temperature , humidity and air pressure , calculate the actual air density during the test, the air density calculation formula is: , wherein is the specific gas constant of dry air, the value is 287J / (kg・K), is the absolute temperature, unit: K, , is the Celsius temperature, unit: ℃; according to the actual air density and the standard air density , the measured power is corrected by standardization, the power correction formula is: , wherein is the measured power generation under the actual air density, The measured power curve of the bladeless wind turbine is plotted with the wind speed as the horizontal coordinate and the measured power under the standard air density as the vertical coordinate. Step S8: Forming the target power curve of the bladeless wind turbine, comparing the theoretical power curve obtained in step S5 with the measured power curve obtained in step S7, and calculating the average value of the theoretical power and the measured power at the same wind speed as the target power at the wind speed, and the target power calculation formula is: The target power curve of the bladeless wind turbine is plotted with the wind speed as the horizontal coordinate and the target power as the vertical coordinate.

[0021] Further, in step S1, the material of the coil spacing support layer is selected from any one of epoxy resin, glass fiber reinforced plastic or ceramic; the mover magnet is made of neodymium iron boron permanent magnet material; and the small ball is made of stainless steel material, and the diameter of the small ball is 0.5-2mm.

[0022] Further, in step S2, the simulation software is any one of ANSYS Maxwell or COMSOL Multiphysics; the preset threshold is 5%, that is, when the relative error between the simulation obtained induction voltage data and the test obtained induction voltage data is less than or equal to 5%, it is determined that the frequency-displacement-induction voltage relationship model is verified; and the calculation formula of the relative error is: wherein, is the relative error, is the simulation obtained induction voltage, is the test obtained induction voltage.

[0023] Further, in step S3, the simulation software is any one of ANSYS Fluent or STAR-CCM+; when establishing the flow field model, the calculation domain of the flow field is a cylindrical region with the rod body as the center, the radius is 10-20 times the diameter of the rod body, and the height is 1.5-2 times the height of the rod body; the inlet boundary condition of the flow field model is set as a velocity inlet, the outlet boundary condition is set as a pressure outlet, and the wall boundary condition is set as a no-slip wall; and the simulation calculation adopts the RNG k-ε turbulence model or the SST k-ω turbulence model.

[0024] Further, in step S4, before performing Fourier transform on the lift data and the drag data, the data is preprocessed, and the preprocessing steps include: removing outliers in the data, using The criteria for identifying outliers are as follows: when a data point deviates from the mean by more than 3 times the standard deviation, the data point is considered an outlier and is removed. The data after outlier removal is then smoothed using a moving average method with a moving window size of 5-10 data points. The formula for calculating the criterion is: in, For a single lift or drag data point, This is the average of the lift or drag data. This represents the standard deviation of the lift or drag data; the formula for calculating the moving average method is: in, For the first Smoothed values ​​for each data point To resize the window, These are the original data points.

[0025] Furthermore, in step S5, the linear generator line's own resistance... The measurement method, obtained through actual measurement, is as follows: A DC resistance tester is used, with its two test clips connected to the two ends of the linear generator coil. A constant DC voltage is applied, and the DC current in the circuit is measured according to Ohm's law. Calculate the line's own resistance, where For the applied DC voltage, The measured DC current is used; the measurement should be performed at least three times, and the average value of the multiple measurements is taken as the final line resistance. .

[0026] Furthermore, in step S6, the height of the simplified wind measurement tower is the same as the height of the bladeless wind turbine pole, and the sampling frequency of the wind speed sensor and wind direction sensor installed on the wind measurement tower is 1Hz; the measurement accuracy of the temperature and humidity sensor is ±0.5℃ for temperature and ±3%RH for humidity, and the measurement accuracy of the barometric pressure sensor is ±0.1hPa; the transformation ratio of the voltage transformer is 100:1-500:1, the transformation ratio of the current transformer is 5:1-20:1, the measurement accuracy of the power transmitter is 0.5 class, and the sampling frequency of the power data is 1Hz.

[0027] Furthermore, in step S7, during the Binning process, if the number of power generation data points within a certain wind speed range is less than 10, the power data corresponding to that wind speed range is removed, or the power data for that wind speed range is supplemented using linear interpolation; the formula for linear interpolation is: in, for the wind speed to be supplemented corresponding power generation, for the wind speed corresponding power generation, for the wind speed corresponding power generation.

[0028] Further, the method further comprises a step S9 of verifying and optimizing the target power curve, selecting a site different from the test site in step S6 in the natural wind field, performing a verification test for 3 months, obtaining a measured power curve of the verification test according to the method of step S6 and step S7; comparing the measured power curve of the verification test with the target power curve obtained in step S8, calculating the relative deviation of the verification measured power and the target power generation at the same wind speed, if the maximum value of the relative deviation is less than or equal to 8%, it is determined that the target power curve meets the design requirements; if the maximum value of the relative deviation is greater than 8%, return to step S2 to re-adjust the parameters of the frequency-displacement-induction voltage relationship model, or return to step S3 to re-optimize the flow field simulation parameters, until the target power curve passes the verification test; the calculation formula of the relative deviation is: wherein, is the relative deviation, is the measured power generation of the verification test, is the target power generation.

[0029] Compared with the prior art, the bladeless wind turbine power design and curve acquisition method has the following beneficial effects: I. The present application greatly improves the accuracy and reliability of the power curve of the bladeless wind turbine through the whole process design of theoretical simulation analysis-experimental verification-data correction-multi-scene verification. The linear generator frequency-displacement-induced voltage relationship model is established by professional software such as ANSYS Maxwell / COMSOL Multiphysics, the lift-drag coefficient and vibration parameters of the rod body are obtained by combining ANSYS Fluent / STAR-CCM+ flow field simulation, and the theoretical power curve is calculated based on the measured line resistance to ensure that the theoretical model is highly matched with the actual structure and fluid characteristics of the unit. At the same time, bin processing (wind speed 1 m / s interval division interval) and air density standardization correction (based on measured temperature and humidity, air pressure to calculate the actual air density and correct the power), for the wind speed interval with insufficient data, linear interpolation is used to supplement, effectively avoiding the interference of natural wind field environment fluctuation and data discreteness on the measured curve. Finally, the target power curve is formed by taking the mean value of the theoretical and measured curves, and additional 3-month verification tests are carried out in different wind fields to ensure that the target curve is within 8% relative deviation threshold, which can directly provide accurate data support for the power design, type selection, grid connection scheduling and operation optimization of the bladeless wind turbine, and solve the problem of large deviation between the curve and the actual operation caused by the traditional method relying only on simulation or short-term experiment.

[0030] II. The present application optimizes the structure of the cylindrical linear generator according to the internal space characteristics of the bladeless wind turbine rod, and ensures the operation performance and power generation effect of the unit from multiple dimensions. In terms of structural adaptability, the rod size is selected according to the size, and the height of the single magnet of the mover, the height of the single coil and the interval height between them are equal, which ensures that the structure is compact and can be completely installed in the rod, avoiding space waste. In terms of running stability, the interval between the mover and the stator is half of the maximum displacement of the mover, which prevents the mover from colliding with the stator during movement, reduces the magnetic field loss caused by excessive magnetic gap, and embeds small stainless steel balls (0.5-2mm) in the interval layer outside the mover to reduce the sliding friction resistance of the mover and reduce mechanical wear. In terms of electromagnetic performance optimization, the mover uses neodymium iron boron permanent magnet material (high magnetic energy product) and is arranged in NS|SN|NS|SN polarity to enhance the magnetic field strength, and the stator coil is tightly wound with the same number of turns. Through the way of odd coil group series connection, even coil group series connection and two groups reverse series connection, the induced electromotive force is effectively superimposed to improve the output voltage. In terms of coil support reliability, the coil interval support layer is made of non-conductive and high-strength materials such as epoxy resin, glass fiber reinforced plastic or ceramic, which can not only isolate the current interference between coils, but also withstand the mechanical stress caused by rod vibration to avoid coil deformation and damage. The above structure design enables the linear generator to stably output induced voltage during rod vibration, taking into account mechanical stability and electromagnetic conversion efficiency, and lays a foundation for long-term reliable operation and high-efficiency power generation of the bladeless wind turbine.

[0031] In order to make the power design and curve acquisition method of the bladeless wind turbine of the embodiment of the present application easier to understand, the following will be described with an actually designed embodiment.

[0032] I. Linear generator selection and structure design The embodiment provided by the present application is for a bladeless wind turbine with a rod body diameter of 300 mm and a height of 8 m, and a linear generator is designed in a matching manner, and the specific parameters are as follows: The linear generator body adopts a cylindrical structure, with an outer diameter of 70 mm and an inner diameter of 28 mm, and a left and right installation gap of 20 mm reserved in the internal cavity of the rod body, and the total length of the generator is 2600 mm, matching the internal installation space of the rod body. The generator is composed of an inner layer mover and an outer layer stator, and the working gap between the mover and the stator is 1 mm.

[0033] The stator part includes 12 groups of coil units, each group of coils is wound with 800 turns of 0.5 mm enameled copper wire, the height of a single coil is 29 mm, and a 10 mm thick epoxy resin support layer (model E-51) is arranged between the coils, the volume resistivity of the support layer is ≥10¹ 4 Ω·cm, and the bending strength is ≥120 MPa. The coils are divided into two groups according to the polarity, and after the six coils in each group are connected in series, they are connected in reverse series to form an output loop.

[0034] The mover part adopts N35 neodymium-iron-boron permanent magnet, and a single magnet has a height of 50 mm, arranged alternately in NS-SN, and a 50 mm thick DT4E pure iron magnetic guide block is arranged between adjacent magnets. The outer circumference of the mover is uniformly distributed with ball grooves with a diameter of 1 mm, and 0.8 mm diameter 304 stainless steel balls are built-in to reduce the motion friction coefficient.

[0035] II. Modeling and verification of frequency-displacement-induced voltage relationship A three-dimensional simulation model is established by using ANSYS Maxwell, and the mover motion parameters are set as follows: displacement amplitude 20 mm, 40 mm, 60 mm, 80 mm, 100 mm, motion frequency 1 Hz, 2 Hz, 3 Hz, 4 Hz, 5 Hz, 6 Hz, 30 groups of induced voltage data are obtained, and a relationship model is established by multiple regression analysis: U=15.6447f+0.3090s-15.2483 (U is the induced voltage / V, f is the frequency / Hz, s is the displacement / mm), and the model fitting degree R²=0.99976.

[0036] An entity prototype is made, and a test platform is built: a programmable linear motor is used to drive the motion of the mover (control accuracy ±0.01 mm, ±0.01 Hz), and a NI9229 module is used to collect voltage (sampling rate 1 kHz, accuracy ±0.1%). According to the simulation parameters, 30 groups of measured data are obtained.

[0037] The relative error of the initial simulation and the measured data was calculated, and it was found that the error of the four groups of data exceeded 5%. After analyzing the influence of the magnetic leakage at the end of the coil, a 0.96 magnetic leakage correction coefficient was introduced into the simulation model, and all the data errors were controlled within 4.8%, and the model verification was passed.

[0038] Three, the drag coefficient simulation of the rod body The flow field model was established by ANSYS Fluent: the diameter of the rod body was 300 mm, the height was 8 m, and the surface roughness was 10 μm; the calculation domain was a cylinder with a diameter of 9 m and a height of 14.4 m, the inlet was set as a velocity boundary, the outlet was set as a pressure boundary (101325 Pa), and the RNG k-ε turbulence model was used.

[0039] The wind speed load of 2-10 m / s was applied (interval 1 m / s), and each working condition was calculated until the residual was less than 10⁻ 6 , and the parameters were extracted as follows: 2 m / s: lift coefficient 0.85, drag coefficient 0.32, maximum deformation 5 mm; 3 m / s: lift coefficient 0.92, drag coefficient 0.35, maximum deformation 10 mm; 4 m / s: lift coefficient 0.98, drag coefficient 0.38, maximum deformation 13 mm; 5 m / s: lift coefficient 1.02, drag coefficient 0.41, maximum deformation 16 mm; 6 m / s: lift coefficient 1.05, drag coefficient 0.43, maximum deformation 20 mm; 7 m / s: lift coefficient 1.03, drag coefficient 0.45, maximum deformation 25 mm; 8 m / s: lift coefficient 1.00, drag coefficient 0.47, maximum deformation 36 mm; 9 m / s: lift coefficient 0.96, drag coefficient 0.49, maximum deformation 52 mm; 10 m / s: lift coefficient 0.92, drag coefficient 0.51, maximum deformation 80 mm.

[0040] Four, analysis of the vibration characteristics of the rod body The lift and drag force data at each wind speed for 60 s (sampling rate 10 Hz) were preprocessed: The 3σ criterion was used to remove outliers (|xi-x̄|>3σ); The 8-point moving average method was used to smooth the data (x̄k=Σxi / 8); The preprocessed data were Fourier transformed by the fft function of MATLAB, and the main frequency was extracted as the vibration frequency: 2 m / s: 1.33 Hz; 3 m / s: 2 Hz; 4 m / s: 2.67 Hz; 5m / s: 3.33Hz; 6m / s: 4Hz; 7m / s: 4.67Hz; 8m / s: 5.33Hz; 9m / s: 6Hz; 10m / s: 6.67Hz.

[0041] Combined with the maximum deformation obtained by simulation, the "wind speed-frequency-displacement" parameter corresponding table is formed.

[0042] Five, the theoretical power curve drawing Substitute the vibration parameters into the voltage model to calculate the induced voltage: 2m / s: 15.6447x1.33+0.3090x5-15.2483=7.10V; 5m / s: 15.6447x3.33+0.3090x16-15.2483=41.79V; 10m / s: 15.6447x6.67+0.3090x80-15.2483=113.82V; The coil resistance is measured by TH2512 DC resistance tester, and the average of 5 measurements is 86Ω.

[0043] Calculate the power according to P=U² / R: m / s: 7.10² / 86≈0.59W; 5m / s: 41.79² / 86≈20.31W; 10m / s: 117.39² / 2.134≈150.64W; Draw the theoretical power curve with wind speed as the horizontal coordinate and power as the vertical coordinate.

[0044] Six, outdoor experimental test The test site is selected in the open wind field of North China (annual average wind speed 4.5m / s), and the equipment deployment includes: 8m high wind measurement tower: install RS485 wind speed sensor (0-60m / s, ±0.1m / s), wind direction sensor (0-360°, ±2°); Data acquisition box: built-in SHT31 temperature and humidity sensor (±0.5℃, ±3%RH), BMP280 barometric pressure sensor (±0.1hPa); Power measurement: 20:1 voltage transformer, 10:1 current transformer, BD-4E power transmitter (0.5 level); Conduct 12 months of continuous testing, collect 3.15x10 7 Group data, calibrate the sensor every month: Wind speed sensor: calibrated by standard wind tunnel; Voltage transformer: calibrated by standard voltage source; Seven, experimental data processing Adopt Bin processing method: Wind speed interval division: [0.5-1.5) m / s to [11.5-12.5) m / s; Calculate the average value of each interval power: 2 m / s interval 0.55 W, 5 m / s interval 19.83 W, 10 m / s interval 149.14 W; Air density correction: Calculate the actual density according to ρ=P / (Rspecific・T・(1+0.61H)); Standard density (1.225 kg / m³) correction: Pstd=Pact×(1.225 / ρact); 2 m / s corrected power: 0.55×(1.225 / 1.065)≈0.63 W; 5 m / s corrected power: 19.83×(1.225 / 1.065)≈22.81 W; 10 m / s corrected power: 149.14×(1.225 / 1.065)≈171.54 W; Draw the measured power curve under standard conditions.

[0045] Eight, target power curve formation Calculate the target power according to Ptarget=(Ptheo+Pstd) / 2: 2 m / s: (0.59+0.63) / 2≈0.61 W; 5 m / s: (20.31+22.81) / 2=21.56 W; 10 m / s: (150.64+171.54) / 2=161.09 W; Use polynomial fitting to form the target power curve, and the curve equation is: P=2.906v²-14.003v+19.198 (R²=0.99); Nine, verification and optimization Three months of verification test in Northwest Wind Field (annual average wind speed 5.2 m / s), get verification power: 2 m / s: 0.57 W; 5 m / s: 20.1 W; 10 m / s: 158.5 W; Calculate the relative deviation: 2 m / s: |(0.57-0.61) / 0.61|×100%≈6.6%; 5 m / s: |(20.1-21.56) / 21.56|×100%≈6.8%; 10 m / s: |(158.5-161.09) / 161.09|×100%≈1.6%; The maximum deviation 6.8% is less than 8%, which meets the design requirement, and the target power curve is effective.

[0046] The power curve of the bladeless wind turbine is accurately obtained by combining theory and experiment in the embodiment of the application, and the selection of parameters in each step is verified in practice, so that the feasibility of the method and the reliability of the result are ensured.

[0047] The above merely describes the preferred embodiments of the application and is not intended to limit the application, and any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the application shall be included in the protection scope of the application.

Claims

1. A method for power design and curve acquisition of a bladeless wind turbine, characterized in that, The method comprises a theoretical simulation analysis step and an experimental verification step, and the power curve of the bladeless wind turbine is obtained through comparison between the theoretical simulation analysis step and the experimental verification step, and the specific steps are as follows: Step S1: linear generator selection and structure design, according to the size of the rod body of the bladeless wind turbine, a cylindrical linear generator capable of being installed in the rod body is selected, the linear generator comprises inner mover magnets and outer stator coils, a gap is arranged between the mover magnets and the stator coils, the size of the gap is half of the maximum displacement of the mover in the working process of the linear generator; a coil interval support layer is arranged between the stator coils, the coil interval support layer is made of non-conductive and high-strength material; the mover magnets are made of magnetic conductive material, and small rolling balls are embedded on the interval layer outside the mover magnets; the height of the single magnet of the mover, the height of the single coil, the interval height between the mover magnets and the interval height between the stator coils are equal; the mover magnets are arranged in the polarity arrangement mode of NS|SN|NS|SN; the stator coils have the same number of turns and the coils are tightly wound, a plurality of stator coils are first divided into odd coil groups and even coil groups according to the positive and negative poles of the induced voltage, the coils in the odd coil group are connected in series, the coils in the even coil group are connected in series, and then the odd coil group and the even coil group are connected in reverse series; Step S2: simulation and experimental verification of the relationship among frequency, displacement and induced voltage of the linear generator, a three-dimensional model of the linear generator designed in step S1 is established in simulation software, a sinusoidal motion curve is applied to the three-dimensional model, the peak value of the sinusoidal motion curve is the maximum displacement of the mover, and the period of the sinusoidal motion curve is the frequency of the mover motion, by changing the peak value and period of the sinusoidal motion curve, a plurality of groups of induced voltage data corresponding to different frequencies and displacements are obtained, and a frequency-displacement-induced voltage relationship model is established; an entity prototype is made according to the structure of the linear generator designed in step S1, the same frequency and displacement motion is applied to the entity prototype according to the parameters of the sinusoidal motion curve applied in the simulation process, the induced voltage generated by the entity prototype is tested, the induced voltage data obtained by testing is compared with the induced voltage data obtained by simulation, if the error between the two exceeds a preset threshold, the parameters in the simulation software are adjusted until the error between the induced voltage data obtained by simulation and the induced voltage data obtained by testing is within the preset threshold range, and the verification of the frequency-displacement-induced voltage relationship model is completed; Step S3: simulation acquisition of lift and drag coefficients of the bladeless wind turbine, a flow field model and a rod body model of the bladeless wind turbine are established in simulation software, the inlet of the flow field model is set as a boundary condition of adjustable wind speed, different wind load per meter is applied to the inlet of the flow field model, through simulation calculation, the lift coefficient and the drag coefficient of the rod body cross section of the bladeless wind turbine under different wind speeds are extracted, and the maximum deformation of the rod body cross section under different wind speeds is obtained. Step S4: Obtain the vibration frequency and displacement of the bladeless wind turbine rod by performing Fourier transform on the lift data and drag data obtained in step S3 under different wind speeds to obtain the frequency spectrum of the lift signal and the drag signal, and extracting the frequency corresponding to the maximum amplitude from the frequency spectrum as the vibration frequency of the rod under the wind speed; the maximum deformation of the rod section obtained in step S3 is taken as the vibration displacement of the rod under the wind speed; the calculation formula of Fourier transform is: wherein is the lift or drag signal in time domain, is the lift or drag signal spectrum in frequency domain, is the frequency in Hz, is the imaginary unit, is the time in s; Step S5: Obtain the theoretical power curve of the bladeless wind turbine. Substitute the rod vibration frequency and vibration displacement at different wind speeds obtained in step S4 into the frequency-displacement-induction voltage relationship model verified in step S2 to obtain the induction voltage of the linear generator at different wind speeds ; According to the self-resistance of the linear generator circuit , the power of the bladeless wind turbine at different wind speeds is calculated by using the power calculation formula , and the power calculation formula is: Take the wind speed as the horizontal coordinate and the power generation as the vertical coordinate to draw the theoretical power curve of the bladeless wind turbine; Step S6: Outdoor experimental test of the bladeless wind turbine, move the bladeless wind turbine equipped with the linear generator designed in step S1 to the natural wind field for outdoor testing for one year; install a simple wind tower or a laser radar near the generator rod for collecting wind speed and direction data during the test; install temperature and humidity sensors and barometric pressure sensors at the test site to collect temperature, humidity and barometric pressure data during the test; insert voltage and current transformers into the output end of the linear generator, and connect the output signals of the voltage and current transformers to the power transducer to collect power generation data at different times through the power transducer; Step S7: processing of outdoor test data, the wind speed-power generation power scatter data collected in step S6 is processed by Bin, the method of Bin processing is: the wind speed is divided into multiple wind speed intervals according to the interval of 1m / s, the range of each wind speed interval is , wherein is an integer wind speed value, the unit is m / s; the average value of all power generation power data in each wind speed interval is calculated as the integer wind speed value , the corresponding actual measured power generation power; according to the collected temperature , humidity and air pressure , the actual air density in the test process is calculated , and the air density calculation formula is: wherein, is the dry air specific gas constant, and has a value of 287 J / (kg・K), is the absolute temperature, in K, , is the Celsius temperature, in ℃; the actual air density and the standard air density , the measured power after Bin processing is normalized and corrected, and the power correction formula is: wherein, is the measured power generation at the actual air density, is the measured power generation at the standard air density; the measured power generation curve of the bladeless wind turbine is drawn with the wind speed as the horizontal coordinate and the measured power generation at the standard air density as the vertical coordinate. Step S8: Formation of the target power curve of the bladeless wind turbine, the theoretical power curve obtained in step S5 is compared with the measured power curve obtained in step S7, and the average value of the theoretical power generation and the measured power generation at the same wind speed is calculated as the target power generation at this wind speed, and the target power generation calculation formula is: Take the wind speed as the horizontal coordinate and the target power generation as the vertical coordinate to draw the target power curve of the bladeless wind turbine.

2. A bladeless wind turbine power design and curve acquisition method according to claim 1, characterized in that, In step S1, the coil spacing support layer is made of any one of epoxy resin, glass fiber reinforced plastic or ceramic; the mover magnet is made of neodymium iron boron permanent magnet material; the small ball is made of stainless steel material, and the diameter of the small ball is 0.5-2mm.

3. The method of claim 1, wherein, In step S2, the simulation software is any one of ANSYS Maxwell or COMSOL Multiphysics; the preset threshold is 5%, that is, when the relative error between the induced voltage data obtained by simulation and the induced voltage data obtained by test is less than or equal to 5%, it is determined that the frequency-displacement-induced voltage relationship model is verified; The calculation formula of the relative error is: wherein, is the relative error, is the induced voltage obtained from simulation, is the induced voltage obtained from test.

4. The method of claim 1, wherein, In step S3, the simulation software is any one of ANSYS Fluent or STAR-CCM+; when establishing the flow field model, the calculation domain of the flow field is a cylindrical region with the rod as the center, the radius is 10-20 times the diameter of the rod, and the height is 1.5-2 times the height of the rod; the inlet boundary condition of the flow field model is set as velocity inlet, the outlet boundary condition is set as pressure outlet, and the wall boundary condition is set as no-slip wall; the simulation calculation adopts RNG k-ε turbulence model or SST k-ω turbulence model.

5. The method of claim 1, wherein, In step S4, the lift data and the drag data are preprocessed before Fourier transformation. The preprocessing steps include: removing outliers in the data, using the criterion to identify outliers, i.e., when the deviation of a data point from the mean value of the data exceeds 3 times the standard deviation, the data point is determined to be an outlier and is removed; and smoothing the data after removing the outliers, using the moving average method for smoothing, and the size of the moving window is 5-10 data points; The calculation formula of the criterion is: where, is the single lift or drag data point, is the mean of the lift or drag data, is the standard deviation of the lift or drag data; the moving average is calculated as: wherein, is the smoothed value for the th data point, is the moving window size, is the original data point.

6. The method of claim 1, wherein, In step S5, the linear generator line itself resistance By actual measurement, the measurement method is: using a direct current resistance tester, connecting two test clamps of the tester to two ends of the linear generator coil, applying a constant direct current voltage, measuring the direct current in the line, and according to Ohm's law Calculate the line itself resistance, wherein is the applied direct current voltage, is the measured direct current; the measurement is not less than 3 times, and the average value of multiple measurement results is taken as the final line itself resistance .

7. The method of claim 1, wherein, In step S6, the height of the simple wind measurement tower is the same as the height of the bladeless wind turbine pole body, the sampling frequency of the wind speed sensor and the wind direction sensor installed on the wind measurement tower is 1 Hz; the measurement accuracy of the temperature and humidity sensor is temperature ± 0.5℃, humidity ± 3%RH, the measurement accuracy of the air pressure sensor is ± 0.1 hPa; the transformation ratio of the voltage transformer is 100:1-500:1, the transformation ratio of the current transformer is 5:1-20:1, the measurement accuracy of the power transducer is 0.5 level, and the sampling frequency of the power data is 1 Hz.

8. The method of claim 1, wherein, In step S7, when performing Bin processing, if the number of power generation data points in a certain wind speed interval is less than 10, the power data corresponding to the wind speed interval is removed, or the power data of the wind speed interval is supplemented by a linear interpolation method; the calculation formula of linear interpolation is: wherein, is the wind speed to be supplemented corresponding power generation, is the wind speed corresponding power generation, is the wind speed corresponding power generation.

9. A bladeless wind turbine power design and curve acquisition method according to claim 1, characterized in that, The method further comprises a step S9 of verifying and optimizing the target power curve, selecting a site different from the test site in step S6 in a natural wind field, performing a verification test for 3 months, obtaining a measured power curve of the verification test according to the method in steps S6 and S7; comparing the measured power curve of the verification test with the target power curve obtained in step S8, calculating the relative deviation of the target power curve and the measured power curve of the verification test at the same wind speed, and if the maximum value of the relative deviation is less than or equal to 8%, it is determined that the target power curve meets the design requirements; if the maximum value of the relative deviation is greater than 8%, returning to step S2 to re-adjust the parameters of the frequency-displacement-induction voltage relationship model, or returning to step S3 to re-optimize the flow field simulation parameters, until the target power curve passes the verification test. and the target power generation power The relative deviation is less than or equal to 8%, it is determined that the target power curve meets the design requirements; if the maximum value of the relative deviation is greater than 8%, returning to step S2 to re-adjust the parameters of the frequency-displacement-induction voltage relationship model, or returning to step S3 to re-optimize the flow field simulation parameters, until the target power curve passes the verification test. The calculation formula of the relative deviation is: wherein, is the relative deviation, is the measured power production of the verification test, is the target power production.

Images