CN110109350A - A kind of power capture optimization method of wave-power device that catching flame algorithm based on chaos moth - Google Patents
A kind of power capture optimization method of wave-power device that catching flame algorithm based on chaos moth Download PDFInfo
- Publication number
- CN110109350A CN110109350A CN201910252076.6A CN201910252076A CN110109350A CN 110109350 A CN110109350 A CN 110109350A CN 201910252076 A CN201910252076 A CN 201910252076A CN 110109350 A CN110109350 A CN 110109350A
- Authority
- CN
- China
- Prior art keywords
- wave
- power
- flame
- moth
- chaos
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Other Liquid Machine Or Engine Such As Wave Power Use (AREA)
Abstract
The invention discloses a kind of power capture optimization methods of wave-power device that flame algorithm is caught based on chaos moth, based on flame algorithm is caught by moth, the new strategy based on chaos operator is added in innovatory algorithm, utilize chaotic maps, it randomly selects 3 elite individuals and does linear combination, the ergodic and randomness for giving full play to Optimizing operator, effectively improve ability of searching optimum.Meanwhile catching flame optimizing mode in spiral and improving, according to light source fitness weight factor, logarithmic spiral formula optimizing route is adjusted, moth is improved and catches flame algorithm and have the shortcomings that solving precision is poor.It by flame automatic adjusument measure, takes into account two kinds and catches the advantage played in searching process between flame mode, electromagnetic linear motor power control parameter is effectively avoided to fall into local optimum because hydrodynamic model is non-linear, solve the B in global optimum's neighborhoodg、Kg, significantly improved so that chaos moth catches the lower mean power exported of flame algorithm optimization control.
Description
Technical field
The present invention relates to the technical field of wave-activated power generation more particularly to a kind of waves that flame algorithm is caught based on chaos moth
The power capture optimization method of power generator.
Background technique
Ocean wave energy is a kind of cleanliness without any pollution new energy, has a very wide distribution, energy density height, has huge hair
Open up potentiality.And wave-power device is a kind of device for extracting ocean energy, is shown from device research and development, experimental verification to practical application
Model has had certain application conditions.As other new energy such as wind energy, ocean wave energy there are very strong randomness,
How having a rest property and uncertain feature, improve wave energy capture rat, is the key that realize wave-power device commercial applications
Place.Maximal power tracing technology (MPPT) can effectively improve sea wave energy capture efficiency, be the core technology of wave-activated power generation.
For direct-drive wave power generation system, when the intrinsic natural frequency of the mover of linear motor is identical as wave frequencies
When, covibration can occur between the two, the output power of direct-drive wave power generation system is maximum at this time, i.e. the capture rate of wave
It is maximum.But since practical wave has very strong time variation, and wave electric power system is an inertia high order system, the matter fixed
Amount and the more difficult change of frequency, this just needs to coordinate by control strategy the changing rule of both sides.But there is presently no a kind of effective
Control method can coordinate the changing rule of both sides, so that the capture rate of wave is reached maximum.
Summary of the invention
It is an object of the invention to overcome the deficiencies in the prior art, propose that one kind can be improved global search performance and convergence
Velocity optimization method quickly calculates linear motor optimized parameter by provided optimization algorithm, utilizes obtained parameter
The anti-electromagnetic force of control linear motor is gone, so that wave-power device frequency is identical as wave frequencies, resonance is generated, to improve
The efficiency of wave energy capture.
To achieve the above object, technical solution provided by the present invention are as follows:
A kind of power capture optimization method of wave-power device that catching flame algorithm based on chaos moth quickly acquires wave
In the case where the optimal damper parameter and elastic parameter of wave electric generating apparatus, flame algorithm is caught by chaos moth and finds out optimal change
The electromagnetic linear motor power control parameter for becoming wave-power device average output power, so that the motion frequency of wave-power device
It is identical as wave frequencies and generate resonance, thus improve wave energy capture efficiency.
Further, specifically includes the following steps:
S1: the MPPT maximum power point tracking control analysis of wave-power device is carried out, wave-power device is found out and averagely exports
The electromagnetic linear motor power control parameter of the calculation formula of power and changeable wave-power device average output power;
S2: the calculation formula based on the step S1 wave-power device average output power found out is established wave power and is caught
Obtain the objective function of optimization problem;
S3: optimizing objective function using chaos moth algorithm, and it is flat to find out optimal changeable wave-power device
The electromagnetic linear motor power control parameter of equal output power so that the motion frequency of wave-power device it is identical as wave frequencies and
Resonance is generated, to improve the efficiency of wave energy capture.
Further, detailed process is as follows by the step S1:
Assuming that cylindrical body float does heaving in perfect fluid, x-y plane is overlapped with standing level, is with float axis
Z-axis establishes vertical coordinate system two-by-two, the velocity potential that incidence wave is generated from x-axis direction are as follows:
(7) in formula, g is acceleration of gravity, and H is wave height, and z is distance of the float center of gravity to standing level, wave number k (ω)=
ω2/ g, d are the depth of water, and ω is angular wave frequency, and velocity potential can be divided into two parts, φrFor the radiation potential generated by floating motion,
φωThe incident gesture generated for wave free movement;
Using wave motion float boundary condition, the wave-excited force that floating body is subject to is acquired according to Green's second theorem:
(8) in formula, FdFor wave-excited force, ρ is the density of water, and s is floating body immersed surface product,For floating body unit normal
Projection of the vector in z-axis;
In radiation potential φrUnder effect, the radiant force of vertical direction cylindrical float generation are as follows:
(9) in formula, FrFor radiant force,For float acceleration,For float speed, A (ω), Brad(ω) point
Not Wei radiant force act on additional mass and radiation damping that lower floating body generates, numerical solution is carried out to A (ω), Brad (ω),
Expression formula are as follows:
Assuming that cylindrical float is in equilibrium state without under wave sea situation, quiet buoyancy and float displacement are positively correlated, indicate
Are as follows:
Fh=-ρ gS0X (t)=- Ksx(t) (11)
(11) in formula, FhFor quiet buoyancy, S0For cylindrical surface waterline sectional area, equivalent elastic coefficient Ks=ρ gS0;
As formula (8-11) it is found that hydrodynamic force suffered by cylindrical float are as follows:
Fwave=Fd+Fr+Fh (12)
According to Newtonian mechanics, the stress condition of float heaving is analyzed in convolution (7)~(12):
(13) in formula, M is float quality, FgIt (t) is the anti-electromagnetic force of linear motor, FvIt (t) is fluid viscous power, Ff(t)
For frictional force;Since the working medium of research object local environment is seawater, ignore fluid viscous power and frictional force suffered by float;
To improve output power, linear motor is needed to provide anti-electromagnetic force F appropriateg(t), wave-power device is controlled
Movement, makes wave-power device motion frequency match with wave frequencies;
Fg(t) expression formula are as follows:
(14) in formula, x (t) is float displacement, Bg、KgFor electromagnetic linear motor power control parameter, wave-activated power generation is respectively indicated
The active and reactive component that system absorbs;
Ocean incidence wave frequency rate will affect wave electric power system output power, formula (9), (11) be substituted into (13), and right
It carries out Fourier transformation, and by the response condition of frequency-domain analysis float, kinetics equation is as follows:
Formula (15) is similar to the phasor form of Second-Order RLC Filter Circuit circuit equation, dynamic using circuit theory simulated waves electricity generation system
Mechanical equation;
Wherein float speedIt is equivalent to electric current Iv, exciting force Fd(j ω) and anti-electromagnetic force Fg(j ω) is equivalent to respectively
Supply voltage and load voltage;By maximum transmitted theory it is found that when resonance occurs for system, load obtains power maximum, i.e. wave
Unrestrained electricity generation system output power is optimal;
Assuming that linear motor magnetic field linear, can be ignored the electromagnetic consumable inside wave-power device, wave electric power system is inhaled
Receiving instantaneous power and mean power may be expressed as:
Wherein, complex power real part indicates that system absorbs mean power, and it is flat can to obtain wave-power device for joint type (9), (10)
Equal output power:
Float mass M and equivalent elastic coefficient KgBelong to the own physical characterisitic parameter of point sucking type wave-power device, fills
Parameter is then definite value after setting fixed installation.
Further, calculating of the step S2 based on the step S1 wave-power device average output power found out is public
Formula establishes the objective function of wave power Optimal harvest problem;Objective function are as follows:
Wherein, η is objective function, indicates the efficiency of wave energy capture, BgAnd KgRespectively the damping parameter of linear motor and
Elastic parameter, ω are wave frequencies, BradAnd KsThe respectively damping parameter and elastic parameter of wave, ω, BradAnd KsThese three ginsengs
Number is by actually measured.
Further, detailed process is as follows by the step S3:
S3-1: setting basic parameter, including chaos moth algorithm greatest iteration number N are 20, chaos moth algorithm population number m
It is 20, dimension D is 2, is randomly generated initial population x (i :), and the matrix that x (i :) indicates 20 rows 2 column respectively represents
20 groups of linear motor damping parameters and elastic parameter;
S3-2: it executes straight line and catches flame mode;
S3-3: with the number of iterations, flame number automatic adjusument;
S3-4: when flame number be less than current iteration number number, switch to execute spiral catch flame mode;
S3-5: algorithm operation reaches maximum number of iterations, and algorithm terminates;Otherwise it is ranked up, according to flame number and repeatedly
Flame mode is caught in generation number, judgement.
Further, the step S3-2 executes straight line and catches flame mode detailed process is as follows:
Using more elite position grouping strategies, first three excellent individual is chosen, carries out chaos linear combination with each individual:
M in formula (2)jFor the current location of elite, MiFor the position of other moths, rnFor chaos optimization operator, embody
Formula are as follows:
rn=4r3 n-1-3rn-1 (3)
When algorithm starts, every moth is all corresponding with the optimal solution and flame of oneself, as the number of iterations increases, flame
Number automatic adjusument, so that moth gradually converges on global optimum.
Further, in the step S3-3, flame adaptive updating rule are as follows:
N is current iteration number in formula (4);
As the number of iterations increases, algorithm, which catches flame mode and switchs to spiral, catches flame, and in search process, artificial moth Mi is utilized
Phototactic characteristics, perceive from itself preferably flames F exiting i, capture flame, the location updating mechanism of step S3-4 around logatithmic spiral are as follows:
S(Mi,Fi)=Di·ebtcos(2πt)+Fi (5)
In formula (5), S (Mi,Fi) be updated position function, b is logarithmic spiral shape constant, i.e. path coefficient, t be [-
1,1] random number, DiFor i-th moth between optimal flame at a distance from, may be expressed as:
Di=| Mi-Fi| (6)
After the completion of catching flame every time, fitness is calculated, is then sorted;Using the moth sorted as the data of next iteration.
Compared with prior art, this programme principle and advantage is as follows:
This programme passes through chaos in the case where quickly acquiring the optimal damper parameter and elastic parameter of wave-power device
Moth catches the electromagnetic linear motor power control parameter that flame algorithm finds out optimal changeable wave-power device average output power,
So that the motion frequency of wave-power device is identical as wave frequencies and generates resonance, to improve the efficiency of wave energy capture.
The power optimization of wave electric power system controls, and its essence is a nonlinear optimal control problems.Flame is caught with moth
It joined the new strategy based on chaos operator based on algorithm, in innovatory algorithm, using chaotic maps, randomly select 3 elite
Individual does linear combination, has given full play to the ergodic and randomness of Optimizing operator, has effectively improved ability of searching optimum.Meanwhile
Flame optimizing mode is caught in spiral to be improved, and according to light source fitness weight factor, is adjusted logarithmic spiral formula optimizing route, is changed
It has been apt to moth and catches flame algorithm to have the shortcomings that solving precision is poor.By flame automatic adjusument measure, two kinds is taken into account and have caught flame mould
The advantage played in searching process between formula effectively avoids electromagnetic linear motor power control parameter because hydrodynamic model is non-linear
Local optimum is fallen into, the B in global optimum's neighborhood is solvedg、Kg, so that chaos moth catches the lower output of flame algorithm optimization control
Mean power significantly improves.
Detailed description of the invention
Fig. 1 is equivalent circuit diagram;
Average output power curve graph when Fig. 2 is T=5s;
Target evolution curve graph when Fig. 3 is T=5s;
Average output power curve graph when Fig. 4 is T=7s;
Target evolution curve graph when Fig. 5 is T=7s.
Specific embodiment
The present invention is further explained in the light of specific embodiments:
A kind of power capture optimization side of wave-power device that catching flame algorithm based on chaos moth described in the present embodiment
Method, comprising the following steps:
S1: the MPPT maximum power point tracking control analysis of wave-power device is carried out, wave-power device is found out and averagely exports
The electromagnetic linear motor power control parameter of the calculation formula of power and changeable wave-power device average output power;This step
It is rapid that detailed process is as follows:
Assuming that cylindrical body float does heaving in perfect fluid, x-y plane is overlapped with standing level, is with float axis
Z-axis establishes vertical coordinate system two-by-two, the velocity potential that incidence wave is generated from x-axis direction are as follows:
(7) in formula, g is acceleration of gravity, and H is wave height, and z is distance of the float center of gravity to standing level, wave number k (ω)=
ω2/ g, d are the depth of water, and ω is angular wave frequency, and velocity potential can be divided into two parts, φrFor the radiation potential generated by floating motion,
φωThe incident gesture generated for wave free movement;
Using wave motion float boundary condition, the wave-excited force that floating body is subject to is acquired according to Green's second theorem:
(8) in formula, FdFor wave-excited force, ρ is the density of water, and s is floating body immersed surface product,For floating body unit normal
Projection of the vector in z-axis;
In radiation potential φrUnder effect, the radiant force of vertical direction cylindrical float generation are as follows:
(9) in formula, FrFor radiant force,For float acceleration,For float speed, A (ω), Brad(ω) point
Not Wei radiant force act on additional mass and radiation damping that lower floating body generates, numerical solution is carried out to A (ω), Brad (ω),
Expression formula are as follows:
Assuming that cylindrical float is in equilibrium state without under wave sea situation, quiet buoyancy and float displacement are positively correlated, indicate
Are as follows:
Fh=-ρ gS0X (t)=- Ksx(t) (11)
(11) in formula, FhFor quiet buoyancy, S0For cylindrical surface waterline sectional area, equivalent elastic coefficient Ks=ρ gS0;
As formula (8-11) it is found that hydrodynamic force suffered by cylindrical float are as follows:
Fwave=Fd+Fr+Fh
(12)
According to Newtonian mechanics, the stress condition of float heaving is analyzed in convolution (7)~(12):
(13) in formula, M is float quality, FgIt (t) is the anti-electromagnetic force of linear motor, FvIt (t) is fluid viscous power, Ff(t)
For frictional force;Since the working medium of research object local environment is seawater, ignore fluid viscous power and frictional force suffered by float;
To improve output power, linear motor is needed to provide anti-electromagnetic force F appropriateg(t), wave-power device is controlled
Movement, makes wave-power device motion frequency match with wave frequencies;
Fg(t) expression formula are as follows:
(14) in formula, x (t) is float displacement, Bg、KgFor electromagnetic linear motor power control parameter, wave-activated power generation is respectively indicated
The active and reactive component that system absorbs;
Ocean incidence wave frequency rate will affect wave electric power system output power, formula (9), (11) be substituted into (13), and right
It carries out Fourier transformation, and by the response condition of frequency-domain analysis float, kinetics equation is as follows:
Formula (15) is similar to the phasor form of Second-Order RLC Filter Circuit circuit equation, dynamic using circuit theory simulated waves electricity generation system
Mechanical equation, equivalent circuit are as shown in Figure 1;
Wherein float speedIt is equivalent to electric current Iv, exciting force Fd(j ω) and anti-electromagnetic force Fg(j ω) is equivalent to respectively
Supply voltage and load voltage;By maximum transmitted theory it is found that when resonance occurs for system, load obtains power maximum, i.e. wave
Unrestrained electricity generation system output power is optimal;
Assuming that linear motor magnetic field linear, can be ignored the electromagnetic consumable inside wave-power device, wave electric power system is inhaled
Receiving instantaneous power and mean power may be expressed as:
Wherein, complex power real part indicates that system absorbs mean power, and it is flat can to obtain wave-power device for joint type (9), (10)
Equal output power:
Float mass M and equivalent elastic coefficient KgBelong to the own physical characterisitic parameter of point sucking type wave-power device, fills
Parameter is then definite value after setting fixed installation.
S2: the calculation formula based on the step S1 wave-power device average output power found out is established wave power and is caught
Obtain the objective function of optimization problem;
Objective function are as follows:
Wherein, η is objective function, indicates the efficiency of wave energy capture, BgAnd KgRespectively the damping parameter of linear motor and
Elastic parameter, ω are wave frequencies, BradAnd KsThe respectively damping parameter and elastic parameter of wave, ω, BradAnd KsThese three ginsengs
Number is by actually measured.
S3: optimizing objective function using chaos moth algorithm, and it is flat to find out optimal changeable wave-power device
The electromagnetic linear motor power control parameter of equal output power, detailed process is as follows:
S3-1: setting basic parameter, including chaos moth algorithm greatest iteration number N are 20, chaos moth algorithm population number m
It is 20, dimension D is 2, is randomly generated initial population x (i :), and the matrix that x (i :) indicates 20 rows 2 column respectively represents
20 groups of linear motor damping parameters and elastic parameter;
S3-2: it executes straight line and catches flame mode;
S3-3: with the number of iterations, flame number automatic adjusument;
S3-4: when flame number be less than current iteration number number, switch to execute spiral catch flame mode;
S3-5: algorithm operation reaches maximum number of iterations, and algorithm terminates;Otherwise it is ranked up, according to flame number and repeatedly
Flame mode is caught in generation number, judgement.
Among the above, step S3-2 executes straight line and catches flame mode detailed process is as follows:
Using more elite position grouping strategies, first three excellent individual is chosen, carries out chaos linear combination with each individual:
M in formula (2)jFor the current location of elite, MiFor the position of other moths, rnFor chaos optimization operator, embody
Formula are as follows:
rn=4r3 n-1-3rn-1 (3)
When algorithm starts, every moth is all corresponding with the optimal solution and flame of oneself, as the number of iterations increases, flame
Number automatic adjusument, so that moth gradually converges on global optimum.
Step S3-3 flame adaptive updating rule are as follows:
N is current iteration number in formula (4);
As the number of iterations increases, algorithm, which catches flame mode and switchs to spiral, catches flame, and in search process, artificial moth Mi is utilized
Phototactic characteristics, perceive from itself preferably flames F exiting i, capture flame, the location updating mechanism of step S3-4 around logatithmic spiral are as follows:
S(Mi,Fi)=Di·ebtcos(2πt)+Fi (5)
In formula (5), S (Mi,Fi) be updated position function, b is logarithmic spiral shape constant, i.e. path coefficient, t be [-
1,1] random number, DiFor i-th moth between optimal flame at a distance from, may be expressed as:
Di=| Mi-Fi| (6)
After the completion of catching flame every time, fitness is calculated, is then sorted;Using the moth sorted as the data of next iteration.
After the electromagnetic linear motor power control parameter for finding out optimal changeable wave-power device average output power, make
Wave-power device motion frequency it is identical with wave frequencies and generate resonance, thus improve wave energy capture efficiency.
In order to confirm the validity of the above method, below based on MATLAB/Simulink platform, wave-activated power generation system is built
System model, carries out simulation analysis.
Wave-power device parameter is as follows: stator inductance Ls=8.3mH, stator resistance Rs=6.48 Ω, rotor permanent magnet magnetic
Chain ψf=0.417Wb, pole span τ=0.05m, float mass M=1972kg, coefficient of elasticity Ks=6127N/m.Chaos moth catches flame
Algorithm parameter is provided that artificial moth number number N=20, maximum number of iterations m=100, system emulation step-length Ts=
0.01s.Under regular waves input condition, verifying chaos moth catches flame algorithm and moth catches the power optimization effect of flame algorithm
Fruit, and pass through the match condition of analysis wave energy acquisition equipment and wave frequencies.
When being respectively 5s and 7s in the regular waves period, A, B suffered by floatradAnd FdAs shown in table 1.In conjunction with table 1
Data, operation is based on chaos moth catches flame algorithm (CMFO), moth catches flame algorithm (MFO) and particle swarm algorithm in Simulink
(PSO) wave electric power system simulation model.Being averaged under the target evolution curve and optimal control under corresponding wave period
Output power curve is as shown in Figure 2-5.
1 wave period of table is the 5, hydrodynamic parameter under 7s and wave-excited force
In the wave electric power system of regular waves, wave period is respectively 5,7s, compares CMFO, MFO and PSO
The lower mean power exported of algorithm optimization control.By Fig. 2,4 it is found that compared with MFO, PSO algorithm, it is based on the control of CMFO algorithm optimization
Under system, the climbing speed for capturing power is significantly improved, and wave electric power system average output power improves 24.8%, 28.1% and
39.4%, 47%, further illustrate that innovatory algorithm can solve rapidly the electromagnetic linear motor power control parameter under optimal load,
System power tracking performance is good.
Fig. 3,5 show that under the different sea situation of wave period, CMFO algorithm converges to global optimum rapidly in or so 10 generations
In range, and in or so 30 to 50 generation of MFO and PSO algorithm, just converges to stationary value.Innovatory algorithm catches flame by introducing chaos optimization
Strategy, ability of searching optimum and convergence rate have obtained larger improvement.
The examples of implementation of the above are only the preferred embodiments of the invention, and implementation model of the invention is not limited with this
It encloses, therefore all shapes according to the present invention, changes made by principle, should all be included within the scope of protection of the present invention.
Claims (7)
1. a kind of power capture optimization method for the wave-power device for catching flame algorithm based on chaos moth, which is characterized in that
In the case where the optimal damper parameter and elastic parameter that quickly acquire wave-power device, flame algorithm is caught by chaos moth and is found out
The electromagnetic linear motor power control parameter of optimal changeable wave-power device average output power, so that wave-power device
Motion frequency it is identical as wave frequencies and generate resonance, thus improve wave energy capture efficiency.
2. a kind of power capture of wave-power device for being caught flame algorithm based on chaos moth according to claim 1 is optimized
Method, which is characterized in that specifically includes the following steps:
S1: the MPPT maximum power point tracking control analysis of wave-power device is carried out, wave-power device average output power is found out
Calculation formula and changeable wave-power device average output power electromagnetic linear motor power control parameter;
S2: it is excellent to establish wave power capture for the calculation formula based on the step S1 wave-power device average output power found out
The objective function of change problem;
S3: optimizing objective function using chaos moth algorithm, and it is average defeated to find out optimal changeable wave-power device
The electromagnetic linear motor power control parameter of power out, so that the motion frequency of wave-power device is identical as wave frequencies and generates
Resonance, to improve the efficiency of wave energy capture.
3. a kind of power capture of wave-power device for being caught flame algorithm based on chaos moth according to claim 2 is optimized
Method, which is characterized in that detailed process is as follows by the step S1:
Assuming that cylindrical body float does heaving in perfect fluid, x-y plane is overlapped with standing level, using float axis as z-axis,
Establish vertical coordinate system two-by-two, the velocity potential that incidence wave is generated from x-axis direction are as follows:
(7) in formula, g is acceleration of gravity, and H is wave height, and z is distance of the float center of gravity to standing level, wave number k (ω)=ω2/ g, d
For the depth of water, ω is angular wave frequency, and velocity potential can be divided into two parts, φrFor the radiation potential generated by floating motion, φωFor wave
The incident gesture that unrestrained free movement generates;
Using wave motion float boundary condition, the wave-excited force that floating body is subject to is acquired according to Green's second theorem:
(8) in formula, FdFor wave-excited force, ρ is the density of water, and s is floating body immersed surface product,For floating body cooler normal vector
Projection in z-axis;
In radiation potential φrUnder effect, the radiant force of vertical direction cylindrical float generation are as follows:
(9) in formula, FrFor radiant force,For float acceleration,For float speed, A (ω), Brad(ω) is respectively spoke
It penetrates power and acts on additional mass and radiation damping that lower floating body generates, numerical solution, expression formula are carried out to A (ω), Brad (ω)
Are as follows:
Assuming that cylindrical float is in equilibrium state without under wave sea situation, quiet buoyancy and float displacement are positively correlated, indicate are as follows:
Fh=-ρ gS0X (t)=- Ksx(t) (11)
(11) in formula, FhFor quiet buoyancy, S0For cylindrical surface waterline sectional area, equivalent elastic coefficient Ks=ρ gS0;
As formula (8-11) it is found that hydrodynamic force suffered by cylindrical float are as follows:
Fwave=Fd+Fr+Fh (12)
According to Newtonian mechanics, the stress condition of float heaving is analyzed in convolution (7)~(12):
(13) in formula, M is float quality, FgIt (t) is the anti-electromagnetic force of linear motor, FvIt (t) is fluid viscous power, FfIt (t) is friction
Power;Since the working medium of research object local environment is seawater, ignore fluid viscous power and frictional force suffered by float;
To improve output power, linear motor is needed to provide anti-electromagnetic force F appropriateg(t), the movement of wave-power device is controlled,
Wave-power device motion frequency is set to match with wave frequencies;
Fg(t) expression formula are as follows:
(14) in formula, x (t) is float displacement, Bg、KgFor electromagnetic linear motor power control parameter, wave electric power system is respectively indicated
The active and reactive component absorbed;
Ocean incidence wave frequency rate will affect wave electric power system output power, by formula (9), (11) substitute into (13) in, and to its into
Row Fourier transformation, by the response condition of frequency-domain analysis float, kinetics equation is as follows:
Formula (15) is similar to the phasor form of Second-Order RLC Filter Circuit circuit equation, utilizes circuit theory simulated waves electricity generation system dynamics
Equation;
Wherein float speedIt is equivalent to electric current Iv, exciting force Fd(j ω) and anti-electromagnetic force Fg(j ω) is equivalent to power supply respectively
Voltage and load voltage;By maximum transmitted theory it is found that when resonance occurs for system, load obtains power maximum, i.e. wavy hair
Electric system output power is optimal;
Assuming that linear motor magnetic field linear, can be ignored the electromagnetic consumable inside wave-power device, wave electric power system absorbs wink
When power and mean power may be expressed as:
Wherein, complex power real part indicates that system absorbs mean power, and it is average defeated can to obtain wave-power device for joint type (9), (10)
Power out:
Float mass M and equivalent elastic coefficient KgBelong to the own physical characterisitic parameter of point sucking type wave-power device, device is fixed
Parameter is then definite value after installation.
4. a kind of power capture of wave-power device for being caught flame algorithm based on chaos moth according to claim 2 is optimized
Method, which is characterized in that calculation formula of the step S2 based on the step S1 wave-power device average output power found out,
Establish the objective function of wave power Optimal harvest problem;Objective function are as follows:
Wherein, η is objective function, indicates the efficiency of wave energy capture, BgAnd KgThe respectively damping parameter and elasticity of linear motor
Parameter, ω are wave frequencies, BradAnd KsThe respectively damping parameter and elastic parameter of wave, ω, BradAnd KsThese three parameters by
It is actually measured.
5. a kind of power capture of wave-power device for being caught flame algorithm based on chaos moth according to claim 2 is optimized
Method, which is characterized in that detailed process is as follows by the step S3:
S3-1: setting basic parameter, including chaos moth algorithm greatest iteration number N is 20, chaos moth algorithm population number m is
20, dimension D are 2, are randomly generated initial population x (i :), and the matrix that x (i :) indicates 20 rows 2 column respectively represents 20
Group linear motor damping parameter and elastic parameter;
S3-2: it executes straight line and catches flame mode;
S3-3: with the number of iterations, flame number automatic adjusument;
S3-4: when flame number be less than current iteration number number, switch to execute spiral catch flame mode;
S3-5: algorithm operation reaches maximum number of iterations, and algorithm terminates;Otherwise it is ranked up, according to flame number and iteration time
Flame mode is caught in number, judgement.
6. a kind of power capture of wave-power device for being caught flame algorithm based on chaos moth according to claim 5 is optimized
Method, which is characterized in that the step S3-2 execution straight line catches flame mode, and detailed process is as follows:
Using more elite position grouping strategies, first three excellent individual is chosen, carries out chaos linear combination with each individual:
M in formula (2)jFor the current location of elite, MiFor the position of other moths, rnFor chaos optimization operator, expression
Are as follows:
rn=4r3 n-1-3rn-1 (3)
When algorithm starts, every moth is all corresponding with the optimal solution and flame of oneself, as the number of iterations increases, flame number
Automatic adjusument, so that moth gradually converges on global optimum.
7. a kind of power capture of wave-power device for being caught flame algorithm based on chaos moth according to claim 5 is optimized
Method, which is characterized in that in the step S3-3, flame adaptive updating rule are as follows:
N is current iteration number in formula (4);
As the number of iterations increases, algorithm, which catches flame mode and switchs to spiral, catches flame, and in search process, artificial moth Mi is special using the light that becomes
Property, it perceives from itself preferably flames F exiting i, captures flame, the location updating mechanism of step S3-4 around logatithmic spiral are as follows:
S(Mi,Fi)=Di·ebtcos(2πt)+Fi (5)
In formula (5), S (Mi,Fi) it is updated position function, b is logarithmic spiral shape constant, i.e. path coefficient, and t is [- 1,1]
Random number, DiFor i-th moth between optimal flame at a distance from, may be expressed as:
Di=| Mi-Fi| (6)
After the completion of catching flame every time, fitness is calculated, is then sorted;Using the moth sorted as the data of next iteration.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910252076.6A CN110109350A (en) | 2019-03-29 | 2019-03-29 | A kind of power capture optimization method of wave-power device that catching flame algorithm based on chaos moth |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910252076.6A CN110109350A (en) | 2019-03-29 | 2019-03-29 | A kind of power capture optimization method of wave-power device that catching flame algorithm based on chaos moth |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110109350A true CN110109350A (en) | 2019-08-09 |
Family
ID=67484817
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910252076.6A Pending CN110109350A (en) | 2019-03-29 | 2019-03-29 | A kind of power capture optimization method of wave-power device that catching flame algorithm based on chaos moth |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110109350A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115406882A (en) * | 2022-10-31 | 2022-11-29 | 常州安控电器成套设备有限公司 | GBDT and improved MFO-based water quality pollutant detection method |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104090490A (en) * | 2014-07-04 | 2014-10-08 | 北京工业大学 | Input shaper closed-loop control method based on chaotic particle swarm optimization algorithm |
CN106710222A (en) * | 2017-03-22 | 2017-05-24 | 广东工业大学 | Method and device for predicting traffic flow |
CN108062585A (en) * | 2017-12-29 | 2018-05-22 | 航天科工智慧产业发展有限公司 | A kind of method that Function Extreme value is calculated based on a flying moth darts into the fire algorithm |
CN108899936A (en) * | 2018-08-31 | 2018-11-27 | 广东工业大学 | A kind of wave-activated power generation method based on simulated annealing particle swarm algorithm |
CN109118016A (en) * | 2018-08-31 | 2019-01-01 | 广东工业大学 | A kind of oscillating float type wave-activated power generation method based on Multiple-population Genetic Algorithm |
CN109344994A (en) * | 2018-08-28 | 2019-02-15 | 温州大学 | A kind of prediction model method based on improvement moth optimization algorithm |
-
2019
- 2019-03-29 CN CN201910252076.6A patent/CN110109350A/en active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104090490A (en) * | 2014-07-04 | 2014-10-08 | 北京工业大学 | Input shaper closed-loop control method based on chaotic particle swarm optimization algorithm |
CN106710222A (en) * | 2017-03-22 | 2017-05-24 | 广东工业大学 | Method and device for predicting traffic flow |
CN108062585A (en) * | 2017-12-29 | 2018-05-22 | 航天科工智慧产业发展有限公司 | A kind of method that Function Extreme value is calculated based on a flying moth darts into the fire algorithm |
CN109344994A (en) * | 2018-08-28 | 2019-02-15 | 温州大学 | A kind of prediction model method based on improvement moth optimization algorithm |
CN108899936A (en) * | 2018-08-31 | 2018-11-27 | 广东工业大学 | A kind of wave-activated power generation method based on simulated annealing particle swarm algorithm |
CN109118016A (en) * | 2018-08-31 | 2019-01-01 | 广东工业大学 | A kind of oscillating float type wave-activated power generation method based on Multiple-population Genetic Algorithm |
Non-Patent Citations (4)
Title |
---|
MINGJING WANG等: "Toward an optimal kernel extreme learning machine using a chaotic moth-flame optimization strategy with applications in medical diagnoses", 《NEUROCOMPUTING》 * |
吴伟民等: "基于增强个体信息交流的蜻蜓算法", 《计算机工程与应用》 * |
吴伟民等: "改进飞蛾捕焰算法在网络流量预测中的应用", 《计算机工程》 * |
曲泉铀等: "浮摆式波能发电装置浮体系统的数值模拟", 《水力发电学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115406882A (en) * | 2022-10-31 | 2022-11-29 | 常州安控电器成套设备有限公司 | GBDT and improved MFO-based water quality pollutant detection method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hasanien | Whale optimisation algorithm for automatic generation control of interconnected modern power systems including renewable energy sources | |
CN108899936A (en) | A kind of wave-activated power generation method based on simulated annealing particle swarm algorithm | |
CN110071496A (en) | A kind of configuration of direct-current grid power optimization and operation method based on wave-activated power generation | |
CN109118016A (en) | A kind of oscillating float type wave-activated power generation method based on Multiple-population Genetic Algorithm | |
CN110206686A (en) | A kind of adaptive maximum power tracking and controlling method for wind power generating set | |
Mahdy et al. | Modeling and optimal operation of hybrid wave energy and PV system feeding supercharging stations based on golden jackal optimal control strategy | |
CN105631152A (en) | PSO (Particle Swarm Optimization)-based wind energy capturing method of VSVF (Variable-speed Variable-frequency) wind power system | |
CN110109350A (en) | A kind of power capture optimization method of wave-power device that catching flame algorithm based on chaos moth | |
Mahdy et al. | Nonlinear modeling and real-time simulation of a grid-connected AWS wave energy conversion system | |
Jama et al. | Self-tunable fuzzy logic controller for the optimization of heaving wave energy converters | |
Magaña et al. | Maximization of wave energy conversion into electricity using oscillating water columns and nonlinear model predictive control | |
Han et al. | The uncertainty and its influence of wind generated power on power system transient stability under different penetration | |
CN110145432B (en) | Wave power generation device power control method based on Fourier analysis and improved Grey wolf algorithm | |
Martinelli et al. | A wave-2-wire experimental investigation of the new" seabreath† wave energy converter: the hydraulic response | |
Cantarellas et al. | Adaptive vector control of wave energy converters | |
Guha | Non-integer disturbance observer-aided resilient frequency controller applied to hybrid power system | |
CN109950933A (en) | A kind of wind-solar-storage joint peak regulation optimization method based on improvement particle swarm algorithm | |
Vanitha et al. | Solving non-convex economic load dispatch problem by efficient hybrid simulated annealing algorithm | |
CN110136025A (en) | A method of improving sea wave energy capture efficiency | |
Syifa et al. | Frequency stability analysis on optimization of virtual inertia controller settings based on retired electric vehicles battery using firefly algorithm | |
Zhang et al. | Individual pitch control for mitigation of power fluctuation of variable speed wind turbines | |
Ahmed et al. | Optimizing distributed generation sizing and siting using particle swarm optimization: A comparative study | |
Zhang et al. | Mitigation of wind power fluctuation by active current control of variable speed wind turbines | |
Rosselan et al. | Dolphin echolocation algorithm for optimal sizing of grid-connected photovoltaic system | |
Huang et al. | A global maximum power point tracking control strategy based on particle swarm optimization algorithm for point-absorber-type wave energy converters |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190809 |
|
WD01 | Invention patent application deemed withdrawn after publication |