CN109139368B - Wind power magnetic suspension cabin suspension eddy current damping optimization method - Google Patents

Wind power magnetic suspension cabin suspension eddy current damping optimization method Download PDF

Info

Publication number
CN109139368B
CN109139368B CN201811022687.3A CN201811022687A CN109139368B CN 109139368 B CN109139368 B CN 109139368B CN 201811022687 A CN201811022687 A CN 201811022687A CN 109139368 B CN109139368 B CN 109139368B
Authority
CN
China
Prior art keywords
suspension
engine room
suspended
air gap
aluminum plate
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811022687.3A
Other languages
Chinese (zh)
Other versions
CN109139368A (en
Inventor
褚晓广
蔡彬
孔英
衣学涛
王强
董同乐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qufu Normal University
Original Assignee
Qufu Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Qufu Normal University filed Critical Qufu Normal University
Priority to CN201811022687.3A priority Critical patent/CN109139368B/en
Publication of CN109139368A publication Critical patent/CN109139368A/en
Application granted granted Critical
Publication of CN109139368B publication Critical patent/CN109139368B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03D—WIND MOTORS
    • F03D7/00—Controlling wind motors 
    • F03D7/02—Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor
    • F03D7/0204—Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor for orientation in relation to wind direction
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/72—Wind turbines with rotation axis in wind direction

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Sustainable Development (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • Sustainable Energy (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • Vibration Prevention Devices (AREA)
  • Magnetic Bearings And Hydrostatic Bearings (AREA)

Abstract

The wind power magnetic suspension cabin disclosed by the invention can yaw and face wind under suspension, so that the yaw power consumption is greatly reduced, but the problems of large fluctuation of a suspended air gap of the cabin, high power consumption of a suspended converter, high failure rate of suspended equipment and the like are easily caused due to random fluctuation of high-frequency wind speed, and the eddy current damping system comprising the suspended winding and the annular aluminum plate is provided, so that the suspended power consumption and the converter burden are effectively reduced, and the suspended oscillation is passively stabilized. The method comprises the following steps of designing a vortex damping aluminum plate, firstly giving a vortex damping force and an engine room suspension force, and constructing an engine room suspension dynamic model; constructing a closed loop Lyapunov function based on the mechanical energy of the suspension of the engine room, and designing the reference setting of the suspension current of the engine room; and taking the pressure under the engine room generated by the rated wind speed as a working condition, comprehensively considering the performance indexes of the air gap fluctuation and the suspension power consumption of the engine room suspension, and optimally designing the thickness of the vortex aluminum plate. The implementation of the invention can effectively improve the performance of the suspended air gap of the engine room and reduce the control burden and the suspended power consumption of the suspended converter of the engine room.

Description

Wind power magnetic suspension cabin suspension eddy current damping optimization method
Technical Field
The invention discloses a method for optimizing the suspended eddy damping of a wind power magnetic suspension engine room, which is applied to a system with large suspended matter mass, more external interference and extremely easy oscillation of a suspended air gap, provides an effective optimization method of an eddy damping system, and belongs to the field of electrical engineering control.
Background
The fan yaw system is a core component of a large and medium-sized wind turbine generator, so that the fan blades face the wind in the front direction, and the maximum wind energy is obtained. At present, a fan yaw system adopts a multi-motor and multi-gear driving mode to realize yaw windward of a cabin, heavier weight of the cabin and a multi-gear driving mechanism, so that the yaw power consumption of the cabin is high and the failure rate is high. The new energy research institute of the university of Qufu David introduces magnetic suspension and disc type motor technology into a fan yaw system, and provides a wind power magnetic suspension yaw system which comprises blades, a suspension winding, a yaw stator and a yaw rotating body integrating a cabin and the suspension winding, wherein when the wind direction changes, the suspension winding is electrified to generate electromagnetic suction to suspend the cabin rotating body, and then under the action of three-phase alternating current of the yaw stator, electromagnetic torque is generated to drive the cabin rotating body to yaw without friction to wind.
The magnetic suspension technology is widely applied to the fields of magnetic suspension trains, magnetic suspension bearings, centrifugal compressors, flywheel energy storage and the like because of the advantages of no contact, no friction, low maintenance cost and the like; however, the problems of nonlinearity, weak damping, unstable open loop and the like of the magnetic suspension essence make the suspension control challenging, and currently, the non-linear methods such as Back steering, self-adaption, sliding mode control and the like are mostly adopted to realize good tracking of the suspension air gap reference, but the problem of suspension force regulation lag caused by the inductance of the suspension winding is weakened; in addition, the problem of weak damping of levitation is also often realized by introducing air gap speed or acceleration in control, but interference noise is easily introduced, especially the problem of hysteresis of levitation current is easily introduced, so that there are many limitations on improving levitation performance. The eddy current damping is introduced to passively stabilize the suspension oscillation of the engine room, and the method is an effective measure for improving the suspension stability and reducing the burden of the suspension converter. The eddy current damping research dates back to 1993 at the earliest, and MacLatchy researches the speed reduction problem of a permanent magnet which vertically falls in a copper pipe; KWak and the like adopt a permanent magnet at the upper end of a cantilever and a cantilever aluminum plate to form an eddy damping system to quickly stabilize the vibration of the cantilever; the Elbuken et al place the aluminum plate under the suspension of the magnet, so as to realize suspension oscillation suppression; the design and the like research the damping arrangement of the copper disc between the high-temperature superconductor and the permanent magnet, and compare the oscillation convergence characteristics from the aspects of damping ratio, joule heat, the oscillation frequency of the permanent magnet and the like; the room is built by introducing a damper consisting of an aluminum plate and a Halbach magnet into a high-speed compressor, so that the passive stabilization of the radial displacement fluctuation of the shaft of the compressor is realized. It should be noted that the eddy current damping is composed of copper or aluminum conductor plate and permanent magnet, and the eddy current damping force mainly comes from the relative displacement change of the permanent magnet. However, the introduction of permanent magnets into the nacelle levitation system undoubtedly increases the levitation weight and the levitation power consumption.
Disclosure of Invention
The main purposes of the invention are as follows: aiming at the defects and blanks in the prior art, the invention provides a wind power magnetic suspension cabin suspension eddy current damping optimization method which is characterized in that an optimized objective function of the cabin aluminum plate thickness is determined by combining an optimized weight with an engine cabin suspension air gap fluctuation index and a suspension power consumption index, and the optimized disc-shaped eddy current aluminum plate thickness is further obtained; the eddy current damping system is used for inhibiting the suspension fluctuation of the wind power magnetic suspension engine room, and a disc-shaped eddy current aluminum plate is arranged on the lower side of the yaw stator and cooperates with the suspension winding on the engine room rotating body to form the eddy current damping system; the core of the eddy current damping system is eddy current damping force which is generated by the axial movement of the suspended winding of the rotating body of the engine room, generates eddy current on the disc-shaped eddy current aluminum plate and is generated under the action of a suspended magnetic field, and is closely related to the suspended current of the engine room, the fluctuation speed of a suspended air gap and the thickness of the aluminum plate. The fluctuation of the suspension air gap of the engine room comes from the pressure change under the engine room caused by high-frequency wind speed, and is calculated by the rigidity of the suspension closed loop of the engine room and the pressure under the engine room under the rated wind speed; the suspension power consumption mainly comes from copper consumption generated by a suspension winding, and the copper consumption calculation working condition is obtained by the product of the square of the suspension current and the internal resistance of the suspension winding under the rated wind speed and the effective suspension air gap.
The design steps of the wind power magnetic suspension engine room suspension eddy current damping optimization method are as follows:
step 1, constructing a wind power magnetic suspension cabin suspension dynamic model
a) Method for calculating eddy damping force F containing disc-shaped eddy aluminum plate by adopting current infinitesimal methodzAnd axial suspension force F
Figure GDA0002309828920000031
Figure GDA0002309828920000032
Wherein: saH is the effective air gap of the engine room suspension, a is the radius of the suspension winding, I is the suspension current, N is the number of turns of the suspension winding, c is the thickness of the disc-shaped eddy current damping plate, sigma is the conductivity of the aluminum plate, mu0And v is the suspension speed of the cabin.
b) The engine room suspension motion equation is
Figure GDA0002309828920000033
Wherein: m is the total weight of the nacelle rotor, fdFor axial downward pressure disturbance。
c) The nacelle suspension converter model is
Wherein: saIs the total area of the winding; rsIs the internal resistance of the winding; u shapedcThe bus voltage u is the duty cycle;
Figure GDA0002309828920000035
d) based on the equilibrium point (h)0,I0) And F (h)0,I0) Linearizing formula (3) and formula (4) as mg
Figure GDA0002309828920000036
Wherein:
Figure GDA0002309828920000037
Figure GDA0002309828920000041
step 2, designing a cabin suspension controller
a) Constructing a Lyapunov function as
Wherein the content of the first and second substances,
Figure GDA0002309828920000043
e＝href-h；ke、kdand kpRespectively positive control coefficients; h ismaxIs the maximum value of the suspended air gap of the engine room.
b) By differentiating the equation (6), the rate of change of the energy function can be obtained as
Figure GDA0002309828920000044
Is provided with
Figure GDA0002309828920000045
Can realize that
Figure GDA0002309828920000046
The nacelle levitation force can be expressed as:
Figure GDA0002309828920000047
the nacelle levitation current is referenced as
Figure GDA0002309828920000048
Wherein:
Figure GDA0002309828920000049
step 3, designing a cabin suspension current tracking controller
a) The formula (9) i is realized by adopting a control strategy combining counter electromotive force compensation and current PI tracking for the formula (5)refReference current tracking, introducing dummy variables
Figure GDA00023098289200000410
And substituting the control equation into (5) to obtain the control equation of the suspension current as
b) The suspended current PI tracking controller is set asWherein k isPITo proportional gain, c0For integral coefficient, the closed loop transfer function of the tracking control of the nacelle suspension current is
Figure GDA0002309828920000051
c) Is provided with
Figure GDA0002309828920000052
kPI＝2ξωnL0-RsWherein damping ξ is 0.3, ωn＝500wad/s。
The output of the nacelle suspension converter is
Figure GDA0002309828920000053
Wherein: e.g. of the typei＝iref-i。
Step 3, thickness optimization method for disc-shaped vortex aluminum plate
a) The fluctuation index of the suspension air gap of the engine room is calculated by the suspension rigidity and the pressure under the engine room at the rated wind speed, and the fluctuation of the suspension air gap of the engine room at the rated wind speed is
Figure GDA0002309828920000054
Wherein: a is the rated wind speed amplitude, and w is the rated wind speed angle radian.
Adopting a frequency domain method, substituting s-j omega into formula (13), and the amplitude of the fluctuation of the suspended air gap of the engine room is
Figure GDA0002309828920000055
Where ω is the angular radian.
b) In view of the relatively small iron loss, the performance index levitation power consumption is set as the copper loss under the effective levitation air gap
Figure GDA0002309828920000056
c) Comprehensively considering the fluctuation of the suspension air gap and the suspension power consumption under the rated wind speed of the engine room, simultaneously considering the difference between the fluctuation of the suspension air gap and the magnitude order of the suspension power consumption, and combining the optimized weight rho1(0.6) and ρ2(0.4) an objective function for optimizing the thickness of the aluminum plate of the nacelle is
Figure GDA0002309828920000061
Wherein: deltaB＝hmax；0＜c≤cmax，0＜I≤Imax；
Figure GDA0002309828920000062
d) And (3) obtaining an extreme value of the formula (16) by adopting a Lagrange extreme value method to obtain the optimized thickness of the disc-shaped vortex aluminum plate.
The invention has the following beneficial effects:
the invention provides an eddy current damping system based on integration of a suspension winding and a yawing stator aluminum plate, wherein the eddy current damping force is caused by air gap fluctuation and suspension weight change. The real-time performance of the nacelle suspension control is improved by the designed nacelle suspension controller based on the Lyapunov function constructed by the nacelle suspension mechanical energy; the introduced eddy current damping system passively stabilizes the suspension oscillation of the engine room, can greatly improve the suspension performance of the engine room, and greatly reduces the burden of the suspension converter.
Drawings
Fig. 1 is a schematic diagram of a wind power magnetic suspension system.
FIG. 2 is a wind nacelle levitation control block diagram.
FIG. 3 is a graph of an optimization objective function of thickness of an aluminum plate of a nacelle as a function of thickness of the aluminum plate.
FIG. 4 is a graph of nacelle suspension performance as a function of aluminum plate thickness.
In the figure: 1 fan paddle, 2 fan cabins, 3 stator cores that driftage, 4 cabin rotators, 5 suspension windings, 6 pylon, 7 suspension converters, 8 discoid vortex aluminum plate.
Detailed Description
The invention is further illustrated by the following figures and examples.
The invention provides a wind power magnetic suspension engine room suspension eddy current damping optimization method which is characterized in that eddy current damping is used for inhibiting suspension fluctuation of a wind power magnetic suspension engine room 2, a disc-shaped eddy current aluminum plate 8 is arranged on the lower side of a yaw stator iron core 3, and is cooperated with a suspension winding 5 on an engine room rotator 4 to form an eddy current damping system; the core of the eddy current damping system is eddy current damping force which is generated by the axial movement of the suspended winding of the rotating body of the engine room, generates eddy current on the disc-shaped eddy current aluminum plate and is generated under the action of a suspended magnetic field, and is closely related to the suspended current of the engine room, the fluctuation speed of a suspended air gap and the thickness of the aluminum plate.
The design steps of the wind power magnetic suspension engine room suspension eddy damping optimization method are as follows:
step 1, constructing a wind power magnetic suspension cabin suspension dynamic model
a) Method for calculating eddy current damping force F of aluminum plate containing eddy current by adopting current infinitesimal methodzAnd axial suspension force F
Figure GDA0002309828920000071
Figure GDA0002309828920000072
Wherein: saH is the effective air gap of the engine room suspension, a is the radius of the suspension winding, I is the suspension current, N is the number of turns of the suspension winding, c is the thickness of the disc-shaped eddy current damping plate, sigma is the conductivity of the aluminum plate, mu0And v is the suspension speed of the cabin.
b) The engine room suspension motion equation is
Figure GDA0002309828920000073
Wherein: m is the total weight of the nacelle rotor, fdIs an axial downward pressure disturbance.
c) The nacelle suspension converter model is
Figure GDA0002309828920000074
Wherein: saIs the total area of the winding; rsIs the internal resistance of the winding; u shapedcThe bus voltage u is the duty cycle;
Figure GDA0002309828920000075
d) based on the equilibrium point (h)0,I0) And F (h)0,I0) Linearizing formula (19) and formula (20) as mg
Figure GDA0002309828920000081
Wherein:
Figure GDA0002309828920000082
Figure GDA0002309828920000083
step 2, designing a cabin suspension controller
a) Constructing a Lyapunov function as
Wherein the content of the first and second substances,
Figure GDA0002309828920000085
e＝href-h；ke、kdand kpRespectively positive control coefficients; h ismaxIs the maximum value of the suspended air gap of the engine room.
b) By differentiating the equation (22), the rate of change of the energy function is obtained
Figure GDA0002309828920000086
Is provided with
Figure GDA0002309828920000087
Can realize that
Figure GDA0002309828920000088
The nacelle levitation force can be expressed as:
Figure GDA0002309828920000089
the nacelle levitation current is referenced as
Wherein:
Figure GDA00023098289200000811
step 3, designing a cabin suspension current tracking controller
a) The formula (25) i is realized by adopting a control strategy combining counter electromotive force compensation and current PI tracking for the formula (21)refReference current tracking, introducing dummy variables
Figure GDA00023098289200000812
And substituting the control equation into (5) to obtain the control equation of the suspension current as
Figure GDA0002309828920000091
b) The suspended current PI tracking controller is set as
Figure GDA0002309828920000092
Wherein k isPITo proportional gain, c0For integral coefficient, the closed loop transfer function of the tracking control of the nacelle suspension current is
Figure GDA0002309828920000093
c) Is provided withkPI＝2ξωnL0-RsWherein damping ξ is 0.3, ωn＝500wad/s。
The output of the nacelle suspension converter is
Figure GDA0002309828920000095
Wherein: e.g. of the typei＝iref-i。
Step 3, thickness optimization method for disc-shaped vortex aluminum plate
a) The fluctuation index of the suspension air gap of the engine room is calculated by the suspension rigidity and the pressure under the engine room at the rated wind speed, and the fluctuation of the suspension air gap of the engine room at the rated wind speed is
Figure GDA0002309828920000096
Wherein: a is the rated wind speed amplitude, and w is the rated wind speed angle radian.
Adopting a frequency domain method, substituting s-j omega into formula (29), and the amplitude of the fluctuation of the suspended air gap of the engine room is
Figure GDA0002309828920000097
Where ω is the angular radian.
b) In view of the relatively small iron loss, the performance index levitation power consumption is set as the copper loss under the effective levitation air gap
Figure GDA0002309828920000098
c) Comprehensively considering the fluctuation of the suspension air gap and the suspension power consumption under the rated wind speed of the engine room, simultaneously considering the difference between the fluctuation of the suspension air gap and the magnitude order of the suspension power consumption, and combining the optimized weight rho1(0.6) and ρ2(0.4) an objective function for optimizing the thickness of the aluminum plate of the nacelle is
Figure GDA0002309828920000101
Wherein: deltaB＝hmax；0＜c≤cmax，0＜I≤Imax；
Figure GDA0002309828920000102
d) And (3) obtaining an extreme value of the formula (32) by adopting a Lagrange extreme value method to obtain the optimized thickness of the disc-shaped vortex aluminum plate.
The invention will be further described below with reference to a preferred embodiment.
Based on a wind power magnetic suspension yaw prototype shown in figure 1 and a 3kW cabin suspension converter, a cabin suspension test platform is built, the cabin suspension weight is 484kg, the optimized range of the aluminum plate thickness is 0-20 mm, and the internal resistance of the converter is 8.4 ohms. The nacelle suspension control is based on a nacelle suspension controller and a suspension current tracking controller shown in fig. 2, and stable suspension of the nacelle is completed.
FIG. 3 shows the trajectory of the variation of the objective function for optimizing the thickness of the aluminum plate, when the thickness of the aluminum plate is 4.95mm, the minimum value of the objective function J is 0.69, and the objective function is further increased as the thickness continues to increase, because the effective levitation air gap is increased and the levitation power consumption is increased due to the increase of the thickness of the aluminum plate; when the thickness of the aluminum plate is less than 4.95mm, eddy damping is relatively small, and the fluctuation of the suspended air gap under the action of the pressure interference under the cabin is relatively large.
FIG. 4 shows the cabin down force f caused by high frequency windd＝2000sin(70t+β0) Under the action, the suspension air gap fluctuates, the suspension power consumption and the damping force change curve. The engine room only takes 1s to reach the suspended target air gap of 10mm, and the effective air gap h is realized at the moment0And a suspension height delta0The engine room is 10mm, the starting without overshoot of the engine room is realized, and the steady state error is only 0.01 mm; in order to verify the influence of the aluminum plate on the suspension performance of the engine room, 2s begins to add pressure interference under the engine room, the thickness of the aluminum plate is linearly increased from 0mm to 10mm, and the suspension fluctuation amount is gradually reduced from 0.5mm to 0.03mm, which is caused by the increase of the thickness of the aluminum plate and the increase of the suspension damping, the eddy current damping force is gradually increased from 100N (0mm) to 700N (4.95 mm in thickness) until 1400N is reached when the thickness is 10mm, and the convergence of air gap fluctuation is accelerated by the eddy current damping force; and the levitation power consumption is gradually increased along with the increase of the thickness of the aluminum plate, the levitation power consumption is 350W at 2s and is 10s (thickness)Degree 10mm) the suspension power consumption is increased to 750W, but the fluctuation amount of the suspension power consumption is opposite to the power consumption change, the fluctuation amount is greatly reduced due to the increase of the thickness of the aluminum plate, the fluctuation amount of the power consumption reaches 600W at 2s, the fluctuation amount is only 80W at 10mm, and the reduction of the fluctuation amount of the suspension power consumption effectively represents the reduction of the control burden of the suspension converter. The comprehensive suspension power consumption and air gap fluctuation change can be known as follows: the optimal thickness of the aluminum plate is 4.95mm, so that the suspension power consumption, air gap fluctuation and the load of the suspension converter are integrally optimal.

Claims (2)

1. A wind power magnetic suspension engine room suspension eddy current damping optimization method is characterized in that eddy current damping is used for inhibiting wind power magnetic suspension engine room suspension fluctuation, a disc-shaped eddy current aluminum plate is arranged on the lower side of a yaw stator, and is cooperated with a suspension winding on an engine room rotating body to form an eddy current damping system; the core of the eddy current damping system is eddy current damping force which is generated by the axial movement of the suspended winding of the engine room rotator, generates eddy current on the disc-shaped eddy current aluminum plate and is generated under the action of a suspended magnetic field, and is closely related to the engine room suspended current, the suspended air gap fluctuation speed and the aluminum plate thickness; determining an optimized objective function of the thickness of the aluminum plate of the engine room by combining the fluctuation index of the suspended air gap of the engine room and the suspended power consumption index with the optimized weight, and further obtaining the optimized thickness of the disc-shaped vortex aluminum plate; the fluctuation of the suspension air gap of the engine room comes from the pressure change under the engine room caused by high-frequency wind speed, and is calculated by the rigidity of the suspension closed loop of the engine room and the pressure under the engine room under the rated wind speed; the suspension power consumption mainly comes from copper consumption generated by a suspension winding, and the copper consumption calculation working condition is obtained by the product of the square of the suspension current and the internal resistance of the suspension winding under the rated wind speed and the effective suspension air gap.
2. The wind power magnetic suspension engine room suspension eddy current damping optimization method according to claim 1, characterized by comprising the following design steps:
step 1, constructing a wind power magnetic suspension cabin suspension dynamic model
a) Method for calculating eddy damping force F containing disc-shaped eddy aluminum plate by adopting current infinitesimal methodzAnd axial suspension force F
Figure FDA0002309828910000011
Wherein: saH is the effective air gap of the engine room suspension, a is the radius of the suspension winding, I is the suspension current, N is the number of turns of the suspension winding, c is the thickness of the disc-shaped eddy current damping plate, sigma is the conductivity of the aluminum plate, mu0Is the vacuum magnetic permeability, v is the cabin suspension velocity;
b) the engine room suspension motion equation is
Figure FDA0002309828910000013
Wherein: m is the total weight of the nacelle rotor, fdAxial downforce interference;
c) the nacelle suspension converter model is
Figure FDA0002309828910000021
Wherein: saIs the total area of the winding, RsIs the internal resistance of the winding, UdcThe bus voltage, u is the duty cycle,
d) based on the equilibrium point (h)0,I0) And F (h)0,I0) Linearizing formula (3) and formula (4) as mg
Figure FDA0002309828910000023
Wherein:
Figure FDA0002309828910000024
step 2, designing a cabin suspension controller
a) Constructing a Lyapunov function as
Figure FDA0002309828910000026
Wherein the content of the first and second substances,
Figure FDA0002309828910000027
e＝href-h，ke、kdand kpRespectively, positive control coefficient, hmaxThe maximum value of the suspended air gap of the engine room;
b) by differentiating the equation (6), the rate of change of the energy function can be obtained as
Figure FDA0002309828910000028
Is provided withCan realize thatThe nacelle levitation force can be expressed as:
Figure FDA0002309828910000031
the nacelle levitation current is referenced as
Figure FDA0002309828910000032
Wherein:
Figure FDA0002309828910000033
step 3, designing a cabin suspension current tracking controller
a) The formula (9) i is realized by adopting a control strategy combining counter electromotive force compensation and current PI tracking for the formula (5)refReference current tracking, introducing dummy variables
Figure FDA0002309828910000034
And substituting the control equation into (5) to obtain the control equation of the suspension current as
Figure FDA0002309828910000035
b) The suspended current PI tracking controller is set as
Figure FDA0002309828910000036
Wherein k isPITo proportional gain, c0For integral coefficient, the closed loop transfer function of the tracking control of the nacelle suspension current is
Figure FDA0002309828910000037
c) Is provided with
Figure FDA0002309828910000038
kPI＝2ξωnL0-RsWherein damping ξ is 0.3, ωn500wad/s, the output of the nacelle suspension converter is
Figure FDA0002309828910000039
Wherein: e.g. of the typei＝iref-i；
Step 3, thickness optimization method for disc-shaped vortex aluminum plate
a) The fluctuation index of the suspension air gap of the engine room is calculated by the suspension rigidity and the pressure under the engine room at the rated wind speed, and the fluctuation of the suspension air gap of the engine room at the rated wind speed is
Figure FDA00023098289100000310
Wherein: a is a rated wind speed amplitude, and w is a rated wind speed angle radian;
adopting a frequency domain method, substituting s-j omega into formula (13), and the amplitude of the fluctuation of the suspended air gap of the engine room is
Figure FDA0002309828910000041
Wherein ω is the angular radian;
b) in view of the relatively small iron loss, the performance index levitation power consumption is set as the copper loss under the effective levitation air gap
c) Comprehensively considering the fluctuation of the suspension air gap and the suspension power consumption under the rated wind speed of the engine room, simultaneously considering the difference between the fluctuation of the suspension air gap and the magnitude order of the suspension power consumption, and combining the optimized weight rho1(0.6) and ρ2(0.4) an objective function for optimizing the thickness of the aluminum plate of the nacelle is
Wherein: deltaB＝hmax；0＜c≤cmax，0＜I≤Imax，
Figure FDA0002309828910000044
d) And (3) obtaining an extreme value of the formula (16) by adopting a Lagrange extreme value method to obtain the optimized thickness of the disc-shaped vortex aluminum plate.
CN201811022687.3A 2018-09-04 2018-09-04 Wind power magnetic suspension cabin suspension eddy current damping optimization method Active CN109139368B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811022687.3A CN109139368B (en) 2018-09-04 2018-09-04 Wind power magnetic suspension cabin suspension eddy current damping optimization method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811022687.3A CN109139368B (en) 2018-09-04 2018-09-04 Wind power magnetic suspension cabin suspension eddy current damping optimization method

Publications (2)

Publication Number Publication Date
CN109139368A CN109139368A (en) 2019-01-04
CN109139368B true CN109139368B (en) 2020-02-07

Family

ID=64826439

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811022687.3A Active CN109139368B (en) 2018-09-04 2018-09-04 Wind power magnetic suspension cabin suspension eddy current damping optimization method

Country Status (1)

Country Link
CN (1) CN109139368B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112798982B (en) * 2020-12-29 2021-12-14 武汉大学 Model-based open-circuit fault diagnosis method and system for three-phase converter power tube

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201486756U (en) * 2009-07-21 2010-05-26 曲阜师范大学日照校区 Magnetic suspension direction regulating device for wind turbine
CN102661245B (en) * 2012-05-15 2013-12-25 中国科学院宁波材料技术与工程研究所 Wind power generator set and yawing-braking magnetic thixotropic flexible braking device of wind power generator set
KR101390277B1 (en) * 2012-11-16 2014-04-29 삼성중공업 주식회사 Wind power generator
CN108019316B (en) * 2018-01-22 2023-04-25 曲阜师范大学 Magnetic suspension wind power yaw system with active and passive cooperative regulation

Also Published As

Publication number Publication date
CN109139368A (en) 2019-01-04

Similar Documents

Publication Publication Date Title
CN108999745B (en) Active and passive suspension control method of wind power magnetic suspension engine room
CN110195686B (en) Two-point suspension type ultra-low speed yaw control method for horizontal axis wind power engine room
CN108183650B (en) Wind power magnetic suspension yaw motor control method based on model prediction control
CN107070342A (en) A kind of control system for permanent-magnet synchronous motor of bringing onto load state observer
CN109236572B (en) Low-wind-speed high-power magnetic suspension vertical axis wind turbine generator and control method thereof
Zhang et al. Direct control of radial displacement for bearingless permanent-magnet-type synchronous motors
Shahgholian et al. State space modeling and eigenvalue analysis of the permanent magnet DC motor drive system
CN109113929B (en) Suspension control method for wind power magnetic suspension engine room
CN104300861A (en) Method for controlling three-phase permanent magnet synchronous motor
Zhang et al. Research on the Vector Control System Based on the Difference Frequency of Wind Turbine Generator
CN107165779B (en) Vertical axis suspension permanent magnet wind driven generator and control method thereof
Nasiri Full digital current control of permanent magnet synchronous motors for vehicular applications
Zhu et al. Rotor vibration control of a bearingless induction motor based on unbalanced force feed-forward compensation and current compensation
CN109139368B (en) Wind power magnetic suspension cabin suspension eddy current damping optimization method
CN202043069U (en) Decoupling controller of five-degree-freedom bearingless synchronous reluctance motor
CN102136822A (en) Five-DOF (freedom of degree) bearingless synchronous reluctance motor decoupling controller and construction method thereof
Kascak et al. Motoring performance of a conical pole-pair separated bearingless electric machine
Tang et al. Active vibration control of the flexible rotor to pass the first bending critical speed in high energy density magnetically suspended motor
CN111622899B (en) Synchronous suspension control method for horizontal axis wind power engine room
CN106330041B (en) Composite excitation formula generator without bearing direct predictive control device based on flux observation
Fang et al. Simulation of speed-control system for PMSM based on sliding mode control
CN103427754A (en) Direct controller of radial displacement of bearing-less asynchronous motor rotor
CN113775474B (en) Vertical axis wind turbine generator suspension control method based on adaptive neural network finite time control
Kante et al. A review paper on modeling and simulation of permanent magnet synchronous generator based on wind energy conversion system
Zhang et al. Fast response sliding mode linear speed control of a magnetic screw motor based on feedforward compensation strategy

Legal Events

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