CN108255062A - The energy saving thrust distribution method of dynamic positioning based on improved differential evolution mechanism - Google Patents
The energy saving thrust distribution method of dynamic positioning based on improved differential evolution mechanism Download PDFInfo
- Publication number
- CN108255062A CN108255062A CN201810067356.5A CN201810067356A CN108255062A CN 108255062 A CN108255062 A CN 108255062A CN 201810067356 A CN201810067356 A CN 201810067356A CN 108255062 A CN108255062 A CN 108255062A
- Authority
- CN
- China
- Prior art keywords
- thrust
- dynamic positioning
- propeller
- max
- value
- 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.)
- Granted
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)
- Feedback Control In General (AREA)
Abstract
The present invention provides a kind of energy saving thrust distribution method of dynamic positioning based on improved differential evolution mechanism, includes the following steps:Step S1:Determine dynamic positioning power demand;Step S2:Establish dynamic positioning power distribution mathematical model;Step S3:Initially dissolve for the power distribution problems in dynamic positioning;Step S4:Mutational formats are selected according to mutation constant λ;Step S5:Carry out mutation operation;Step S6:Carry out crossover operation;Step S7:The adaptive value of all solutions is evaluated, and records optimal solution;Step S8:Judge whether to meet end condition, be, perform step S9, otherwise return to step S4;Step S9:Stop cycle, export optimal solution, power distribution is carried out according to optimal solution.By the present invention, thrust needed for dynamic positioning ship is reasonably allocated to each propeller, on the basis of safety of ship is ensured, ship energy consumption is reduced, lowers shipping transport cost, reduce Pollution From Ships.
Description
Technical field
The present invention relates to shipboard automation field more particularly to a kind of dynamic positioning sections based on improved differential evolution mechanism
It can thrust distribution method.
Background technology
The major function of Ship Dynamic Positioning Systems Based (DPS) is to maintain ship in desired position or driving ship along finger
Fixed track navigation.It mainly does according to the departure between the current location status and desired value of ship and external environment
It disturbs, calculates ship required total thrust and torque in real time, then calculate ship further according to thrust distribution logical method
The thrust and its angle of the required generation of each propeller of equipment, and control instruction is transmitted to propeller, propeller according to
The thrust that the instruction generates requirement completes ship's fix.Therefore, thrust distribution be a ring important in DPS, it by controller with
Push system connects to form an entirety.Propeller is no less than 5 under normal conditions, but is not The more the better, also
The factors such as the interaction between two propellers between propeller and hull are must take into account, while also to take into account ship institute
Required maneuverability.Just there are many kinds of disclosure satisfy that the thrust size of positioning requirements and the group in direction to entire propulsion system at this time
It closes, and the major function of thrust distribution logic is exactly to find out a classic combination, and it is real-time to be translated into control instruction
Distribute to each propeller.
Thrust allocation optimization problems are a multi-constraints optimization problems, it requires propeller real-time high-efficiency reasonably to generate ship
The required thrust of oceangoing ship, and meet the requirement for making energy consumption minimum as possible while a degree of maneuverability is kept.Outstanding
Distribution method can not only improve ship's fix precision, but also possess reduction energy expenditure, reduce mechanical wear and noise etc.
Function.
What dynamic positioning ship was generally equipped is full circle swinging thruster, and such propeller and its propulsion system, which are generally gathered around, to be had
Following restrictive condition:
(1) acc power limitation driven, maximum value that there are one the thrusts that each propeller generates, propeller can not generate
Than the thrust of maximum value bigger.
(2) propeller thrust and its rate of change of angle are restricted, can not change the short time very greatly.
(3) all-direction propeller can be deposited around the rotation of 360 ° of its main shaft, but because closing on the interaction between propeller
In thrust forbidden zone, it is therefore desirable to set a variation range to the rotational angle of each propeller.
In addition, the it is also contemplated that appearance of singular structure.For dynamic positioning ship, when its equipment all propellers all
When possessing the same propulsion azimuth, singular structure will occur.It is even lost at this point, Ship Maneuverability will be greatly reduced
Go, and due to the limitation of propeller rotary speed, whole system will longer time restore, this would be possible to relevant
It is produced in engineer operation into serious consequence.Therefore, it in the research and application of dynamic positioning of vessels thrust distribution, should avoid strange
The generation of diverse structure.Therefore in order to keep the balance between energy consumption and operability, a kind of scheme of compromise is proposed, i.e., according to reality
Demand setting adjustment factor, come balance ship to operability and energy consumption demand.
In the prior art, mainly using following several technical solutions:
Patent《Using the thrust force distribution method of power location system of ship of genetic algorithm》, application number
201110253236.2 the assigning process of thrust is carried out using genetic algorithm approach.
Patent《The thrust distribution method of dynamic positioning system》, application number 201310003584.3, using low-pass filtering and
Two kinds of optimization algorithm Distributed Calculation thruster azimuths of Sequential Quadratic Programming method and thrust size.
Patent《Using the dynamic positioning thrust distributor and its distribution method for dynamically forbidding angle》, application number
201210356674.6 the assigning process of thrust is carried out using Novel Algorithm.
The main method of dynamic positioning of vessels thrust distribution at present is to ignore the energy consumption characteristics of thrust distribution, it is difficult to find and push away
The optimal solution of power distribution.
Invention content
The purpose of the present invention is distributing mathematical model according to constructed dynamic positioning thrust, pass through improved differential evolution
Mechanism efficiently finds the optimal solution for distributing each propeller thrust and direction instruction, is ensureing that propeller resultant force satisfaction control will
On the basis of asking, marine propeller energy consumption is reduced.
To achieve the above object, the present invention uses following technical scheme:A kind of power based on improved differential evolution mechanism
Positioning energy-saving thrust distribution method, includes the following steps:Step S1:Determine dynamic positioning power demand;Step S2:It establishes dynamic
Power Positioning power distributes mathematical model;Step S3:Initially dissolve for the power distribution problems in dynamic positioning;Step S4:Root
Mutational formats are selected according to mutation constant λ;Step S5:Carry out mutation operation;Step S6:Carry out crossover operation;Step S7:Evaluation
The adaptive value of all solutions, and record optimal solution;Step S8:Judge whether to meet end condition, be, perform step S9, otherwise return
Return step S4;Step S9:Stop cycle, export optimal solution, power distribution is carried out according to optimal solution.
In an embodiment of the present invention, step S2 includes step in detail below:Step S21:About thrust assignment problem
Nonlinear optimization mathematical model is as follows:
Constraint adjusts as follows:
S=τ-B (α) F
Fmin≤F≤Fmax
△Fmin≤F-F0≤△Fmax
αmin≤α≤αmax
△αmin≤α-α0≤△αmax
Wherein, W is total energy consumption in first item, and P is weight coefficient, FiFor the thrust of i-th of propeller, kiFor calculating parameter;
Section 2 sTQs is penalty term, and s is broad sense thrust error vector;Weight matrix Q is diagonal positive definite matrix, should be taken large values, and is used
To ensure that error s levels off to zero;Section 3 (α-α0)TΩ(α-α0) it is the pace of change that angle is promoted for constraining, wherein α is
The azimuth of the propeller at this moment, α0The azimuth of propeller for previous moment, weight matrix Ω>0 is used for adjusting optimization
Target;Section 4Be for avoiding singular structure, wherein
lxnAnd lynIt is the X-direction coordinate and Y-direction coordinate of n-th of propeller respectively;X-direction is to bow, Y from ship
Direction is to ship starboard from ship;If propulsion system is unusual or levels off to unusual, i.e. det (B (α) B ' (α)) is equal to zero or near
Zero is similar to, then the value of Section 4 can be very big, is equivalent to penalty;ε in formula>0, δ be more than zero, δ be adjustment factor, for balancing
The energy consumption and maneuverability of ship, the bigger maneuverability of δ values is better, and corresponding energy consumption can also increased, and δ values are smaller, then feelings
Condition is exactly the opposite;Step S22:In constraints, τ=(τx,τy,τM)TIt is desired power and torque, F is propulsion system
Thrust matrix, B (α) F is then the practical conjunction thrust and resultant moment generated of propulsion system, and reality is calculated by B (α) calculation formula
Broad sense thrust error vector between expectation;FmaxAnd FminThe maximum value and minimum value of propeller thrust, limitation are represented respectively
The thrust range of propeller;△FmaxWith △ FminThe bound that propeller thrust changes within the unit interval, B are represented respectively
(α) calculation formula defines the range of thrust variation rate;Correspondingly, αmaxAnd αminFor the range of propeller rotation angle, △ αmax
With △ αminIt is the bound of propeller angle change amplitude between two moment.
In an embodiment of the present invention, step S4 includes the following steps:Step S41:In improved differential evolution, according to such as
Artificial bee colony algorithm is added in the mode of differential evolution algorithm as local searching strategy by lower formula;
In formula, xiFor the individual of carry out neighborhood search chosen, xkFor xiAdjacent individual, x 'iFor obtaining after neighborhood search
The particle arrived, i.e. x 'iIt is xiPass through the random individual x adjacent with itskAmendment solution after being compared to each other,It is random in [- 1,1]
The value of generation;K ∈ { 1,2,3 ..., SN } wherein SN is the quantity of population;Step S42:Population is added according to equation below to calculate
The social recognition part of method optimizes differential evolution algorithm:x′i=xr1+F(xr2-xr3)+ψ(xgbest-xr1);Wherein ψ is range
The numerical value generated at random in [0,1], xgbestIt is globally optimal solution;xr1, xr2And xr3Be the variable that randomly selects and r1 ≠ r2 ≠
r3;F be for controlling perturbation total amount in mutation process and improving the zoom factor of convergence rate, value range be [0,
1];New explanation x 'iBy three partially synthetic generations:The target individual x ' that first part is selectedi;Second part is by randomly choosing
Parent individuality difference generate vector;Last part is according in the target individual by being selected and current entire population
The difference of globally optimal solution carries out the solution vector of difference operation generation;Step S43:According to following steps by this step S41 and
Step S42 improved procedures are combined:If rand<λ is then utilizedOtherwise according to x 'i=xr1+F
(xr2-xr3)+ψ(xgbest-xr1) carry out mutation operation;Rand is a random function for generating random number;Wherein λ is that range exists
[0,1] control variable.
Further, according to the following formula optimizing strategy is selected to adjust λ value:
Wherein t is current cycle time, and T is selected maximum cycle.
Step S5 includes the following steps in an embodiment of the present invention:It is utilized respectively formula
ABC improvement strategies or PSO improvement strategies carry out mutation operation;Wherein ABC improvement strategies include the following steps:1. lead bee neighborhood
Search generates new explanation, and calculates its fitness value;2. according to formulaGenerate solution to be evaluated;3. it obtains
Mutant test group;PSO improvement strategies include the following steps:1. the history optimal value Pbest and the overall situation of each individual of update are most
Figure of merit Gbest;2. according to formula x 'i=xr1+F(xr2-xr3)+ψ(xgbest-xr1) generation solution to be evaluated;3. obtain mutant test group
Body.
In an embodiment of the present invention, step S6 includes the following steps:According to the intersection of equation below design dynamic increment
Parameter
CR=CRmin+(CRmax-CRmin)×(t/T)2
Wherein CRmaxAnd CRminIt is the maximum value and minimum value of cross parameter;T is current cycle time, T be it is selected most
Systemic circulation number.
Compared with prior art, artificial bee colony algorithm and particle cluster algorithm are introduced differential evolution algorithm by the method for the present invention
In, increase the diversity of differential evolution algorithm mutational formats;Mutational formats are selected, and pass through and adjust λ using constant λ is mutated
Value carry out the dynamic regulations of mutational formats, search stresses ABC improvement strategies early period, and the search later stage stresses PSO improvement strategies, with
Improve convergence speed of the algorithm and precision;Intersect the dynamic regulation of the value progress crossover operation of constant CR by adjusting, adjustment is prominent
Become the crossover probability between individual and target individual, improve algorithm the convergence speed and enhance the search precision of algorithm.Present invention side
Method according to the optimizing result of dynamic positioning ship thrust distribution model and improved differential evolution algorithm carry out thrust distribution, make its
On the basis of meeting ship thrust requirements guarantee safety of ship, make the distribution of dynamic positioning ship thrust more reasonable.
Description of the drawings
Fig. 1 is the main flow schematic diagram of the present invention.
Specific embodiment
Explanation is further explained to the present invention in the following with reference to the drawings and specific embodiments.
Being briefly described below for differential evolution algorithm is shown:
Mutation:Mutation is a very important link in the cyclostage, it randomly selects three parents in group
Body generates experimental subjects.In mutation operation, experimental subjects variable x 'iIt obtains as follows:
x′i=xr1+F(xr2-xr3) (1)
Wherein, xr1, xr2And xr3It is the variable that randomly selects and r1 ≠ r2 ≠ r3, F are for controlling in mutation process
The total amount that perturbs and the zoom factor for improving convergence rate, its value range are [0,1].
Intersect:Crossover operation is similar to the exchange and recombination that gene carries out two segments according to certain probability, generates filial generation
Individual.The operation is used to increase population diversity.In this operation, offspring individual is by experimental subjects and parent individuality according to following public affairs
It is generated after formula cross-mixing.
Wherein CR is by our customized intersection constants, RjIt is the real number randomly selected in range [0,1].J represents jth
The corresponding array of a component.
Selection:Selection operation is that target individual offspring individual related to it eliminates competition process.Using adaptive value as according to excellent
Win bad eliminate.To choose the minimum standard of adaptive value, selection operation can be represented with equation below:
Wherein f (x) represents the adaptive value of individual.
Being briefly described below for improved differential evolution algorithm is shown:
In improved differential evolution, we are in a manner that artificial bee colony algorithm is added in differential evolution algorithm by equation below
As local searching strategy.
In formula, x 'iIt is xkAmendment solution after being compared to each other by the random value adjacent with its,It is to be given birth at random in [- 1,1]
Into value.K ∈ { 1,2,3 ..., SN } wherein SN is the quantity of population.xiFor the individual of carry out neighborhood search chosen, xkFor xi
Adjacent individual, x 'iFor the obtained particle after neighborhood search.
In the searching process later stage, since the approximate region of optimal solution has been determined, it would be desirable to accelerate the convergence of algorithm
Property in order to optimizing faster.Therefore we add in the social recognition part of particle cluster algorithm to optimize according to equation below
Differential evolution algorithm.
x′i=xr1+F(xr2-xr3)+ψ(xgbest-xr1) (5)
Wherein ψ is the numerical value generated at random in range [0,1], xgbestIt is globally optimal solution.By formula (5) it is found that new explanation
x′iBy three partially synthetic generations:The target individual x ' that first part is selectedi;Second part is by randomly selected parent
The vector that the difference of body generates;Last part is according to global optimum in the target individual by being selected and current entire population
The difference of solution carries out the solution vector of difference operation generation.This three parts constitutes the new paragon for finding globally optimal solution.
Two improvements mode is combined according to following steps.
If(rand<λ)
Mutation operation is carried out using formula (4);
Else
Mutation operation is carried out using formula (5);
End
Wherein λ is control variable of the range in [0,1];
Since the optimizing strategy that formula (4) provides has stronger ability of searching optimum, in entire optimization process
First half, refined Hook Jeeves algorighm mainly use formula (4).When algorithm enters in follow-up searching process, what formula (5) provided seeks
Dominant strategy has higher convergence rate and precision, and refined Hook Jeeves algorighm mainly carries out optimizing using formula (5).Pass through following public affairs
Formula selects optimizing strategy to adjust λ value:
Wherein t is current cycle time, and T is selected maximum cycle.
It is the Important Parameters for determining crossover probability between mutated individual and target individual to intersect constant CR.CR is larger, then tries
It tests that mutated individual accounting in individual is larger, the field search capability of algorithm can be enhanced and improves algorithm the convergence speed;CR is smaller,
It then tests mutated individual in individual and accounts for that smaller and target individual is larger in the individual accounting of experiment, the search essence of algorithm can be improved
Degree.Therefore, according to the cross parameter of equation below design dynamic increment, the search speed of algorithm early period and the search in later stage are improved
Precision.
CR=CRmin+(CRmax-CRmin)×(t/T)2 (7)
Wherein CRmaxAnd CRminIt is the maximum value and minimum value of cross parameter.
As shown in Figure 1, in a specific embodiment of the invention, a kind of dynamic positioning section based on improved differential evolution mechanism
Energy thrust distribution method, specifically includes following steps:
Step 1:Determine dynamic positioning thrust demand;
Step 2:Establish dynamic positioning power distribution mathematical model;
In an embodiment of the present invention, it is as follows about the nonlinear optimization mathematical model of thrust assignment problem:
Constraint:
S=τ-B (α) F (9)
Fmin≤F≤Fmax (10)
△Fmin≤F-F0≤△Fmax (11)
αmin≤α≤αmax (12)
△αmin≤α-α0≤△αmax (13)
In formula, W is total energy consumption in first item, and P is weight coefficient, is mainly used to adjust ratio of the energy consumption in optimization aim
Weight, the calculation of W have provided, wherein FiFor the thrust of i-th of propeller, kiFor calculating parameter, general value is 0.176.
Section 2 sTQs is penalty term, and s is broad sense thrust error vector, is mainly used to ensure the thrust and torque of propeller
It is enough to offset external interference, completes its basic location tasks.Weight matrix Q is diagonal positive definite matrix, should be taken large values, to
Ensure that error s levels off to zero.
Section 3 (α-α0)TΩ(α-α0) it is the pace of change that angle is promoted for constraining.Wherein α is the propulsion at this moment
The azimuth of device, α0The azimuth of propeller for previous moment, weight matrix Ω>0 is used for adjusting optimization aim.
Section 4Be for avoiding singular structure, wherein
lxnAnd lynBe respectively n-th of propeller X-direction coordinate (to bow from ship) and Y-direction coordinate (from ship to
Ship starboard).If propulsion system is unusual or levels off to unusual, i.e. det (B (α) B ' (α)) is equal to zero or is similar to zero, then the 4th
The value of item can be very big, be equivalent to penalty.ε in formula>0, this is that denominator is zero in order to prevent.δ is more than zero, is adjustment factor,
For balancing the energy consumption of ship and maneuverability, the bigger maneuverability of δ values is better, and corresponding energy consumption can also increased, and δ values are got over
Small, then situation is exactly the opposite.
In addition, in constraints, τ=(τx,τy,τM)TIt is desired power and torque, F is the thrust of propulsion system
Matrix, B (α) F are then the practical conjunction thrust and resultant moment generated of propulsion system, and formula (15) calculates practical between expectation
Broad sense thrust error vector.FmaxAnd FminThe maximum value and minimum value of propeller thrust are represented respectively, and formula (15) then limits
The thrust range of propeller.△FmaxWith △ FminThe bound that propeller thrust changes within the unit interval, formula are represented respectively
(15) range of thrust variation rate is defined.Correspondingly, αmaxAnd αminFor the range of propeller rotation angle, △ αmaxWith △ αmin
It is the bound of propeller angle change amplitude between two moment.
Step 3:Initially dissolve for the power distribution problems in dynamic positioning;
Step 4:Mutational formats are selected according to mutation constant λ;
Two improvements mode is combined according to following steps:
If(rand<λ)
Mutation operation is carried out using formula (4);
Else
Mutation operation is carried out using formula (5);
End
And λ value is adjusted by formula (6).
Step 5:Carry out mutation operation;
It is utilized respectively formula (1), ABC improvement strategies or PSO improvement strategies and carries out mutation operation.
ABC improvement strategies include the following steps:
4. bee neighborhood search is led to generate new explanation, and calculate its fitness value;
5. solution to be evaluated is generated according to formula (4);
6. obtain mutant test group.
PSO improvement strategies include the following steps:
1. the history optimal value Pbest and global optimum Gbest of each individual of update;
2. solution to be evaluated is generated according to formula (5);
3. obtain mutant test group.
Step 6:Carry out crossover operation;
Crossover operation is carried out according to formula (2), the wherein value of CR is determined by formula (7).
Step 7:The adaptive value of all solutions is evaluated, and records optimal solution;
Step 8:Judge whether to meet end condition, be, perform step 9, otherwise return to step 4;
End condition is greatest iteration step number.
Step 9:Stop cycle, export optimal solution, thrust distribution is carried out according to optimal solution.
The above are preferred embodiments of the present invention, all any changes made according to the technical solution of the present invention, and generated function is made
During with range without departing from technical solution of the present invention, all belong to the scope of protection of the present invention.
Claims (6)
1. a kind of energy saving thrust distribution method of dynamic positioning based on improved differential evolution mechanism, it is characterised in that:Including following
Step:
Step S1:Determine dynamic positioning power demand;
Step S2:Establish dynamic positioning power distribution mathematical model;
Step S3:Initially dissolve for the power distribution problems in dynamic positioning;
Step S4:Mutational formats are selected according to mutation constant λ;
Step S5:Carry out mutation operation;
Step S6:Carry out crossover operation;
Step S7:The adaptive value of all solutions is evaluated, and records optimal solution;
Step S8:Judge whether to meet end condition, be, perform step S9, otherwise return to step S4;
Step S9:Stop cycle, export optimal solution, power distribution is carried out according to optimal solution.
2. the energy saving thrust distribution method of the dynamic positioning according to claim 1 based on improved differential evolution mechanism, special
Sign is:Step S2 includes step in detail below:
Step S21:Nonlinear optimization mathematical model about thrust assignment problem is as follows:
Constraint adjusts as follows:
S=τ-B (α) F
Fmin≤F≤Fmax
△Fmin≤F-F0≤△Fmax
αmin≤α≤αmax
△αmin≤α-α0≤△αmax
Wherein, W is total energy consumption in first item, and P is weight coefficient, FiFor the thrust of i-th of propeller, kiFor calculating parameter;Second
Item sTQs is penalty term, and s is broad sense thrust error vector;Weight matrix Q is diagonal positive definite matrix, should be taken large values, to protect
Card error s levels off to zero;Section 3 (α-α0)TΩ(α-α0) it is the pace of change that angle is promoted for constraining, wherein α is this when
The azimuth of the propeller at quarter, α0The azimuth of propeller for previous moment, weight matrix Ω>0 is used for adjusting optimization aim;
Section 4Be for avoiding singular structure, wherein
lxnAnd lynIt is the X-direction coordinate and Y-direction coordinate of n-th of propeller respectively;X-direction is is to bow, Y-direction from ship
To ship starboard from ship;If propulsion system is unusual or levels off to unusual, i.e. det (B (α) B ' (α)) is equal to zero or is similar to zero,
Then the value of Section 4 can be very big, is equivalent to penalty;ε in formula>0, δ be more than zero, δ be adjustment factor, for balancing ship
Energy consumption and maneuverability, the bigger maneuverability of δ values is better, and corresponding energy consumption can also increased, and δ values are smaller, then situation is just
On the contrary;
Step S22:In constraints, τ=(τx,τy,τM)TIt is desired power and torque, F is the moment of thrust of propulsion system
Battle array, B (α) F is then the practical conjunction thrust and resultant moment generated of propulsion system, is calculated by B (α) calculation formula practical and expectation
Between broad sense thrust error vector;FmaxAnd FminThe maximum value and minimum value of propeller thrust are represented respectively, limit propeller
Thrust range;△FmaxWith △ FminThe bound that propeller thrust changes within the unit interval is represented respectively, and B (α) calculates public
Formula defines the range of thrust variation rate;Correspondingly, αmaxAnd αminFor the range of propeller rotation angle, △ αmaxWith △ αminIt is
The bound of propeller angle change amplitude between two moment.
3. the energy saving thrust distribution method of the dynamic positioning according to claim 1 based on improved differential evolution mechanism, special
Sign is:Step S4 includes the following steps:
Step S41:In improved differential evolution, in a manner that artificial bee colony algorithm is added in differential evolution algorithm by equation below
As local searching strategy;
In formula, xiFor the individual of carry out neighborhood search chosen, xkFor xiAdjacent individual, xi' for obtaining after neighborhood search
Particle, i.e. x 'iIt is xiPass through the random individual x adjacent with itskAmendment solution after being compared to each other,It is to be generated at random in [- 1,1]
Value;K ∈ { 1,2,3 ..., SN } wherein SN is the quantity of population;
Step S42:The social recognition part of particle cluster algorithm is added according to equation below to optimize differential evolution algorithm,
x′i=xr1+F(xr2-xr3)+ψ(xgbest-xr1)
Wherein ψ is the numerical value generated at random in range [0,1], xgbestIt is globally optimal solution;xr1, xr2And xr3It randomly selects
Variable and r1 ≠ r2 ≠ r3;F is the zoom factor for controlling perturbation total amount and raising convergence rate in mutation process,
Value range is [0,1];New explanation x 'iBy three partially synthetic generations:The target individual x ' that first part is selectedi;Second part
It is by the vector of the difference generation of randomly selected parent individuality;Last part is according to the target individual by being selected and works as
The difference of globally optimal solution carries out the solution vector of difference operation generation in preceding entire population;
Step S43:This step S41 and step S42 improved procedures are combined according to following steps:If rand<λ is then utilizedOtherwise according to x 'i=xr1+F(xr2-xr3)+ψ(xgbest-xr1) carry out mutation operation;Rand is production
One random function of raw random number;Wherein λ is control variable of the range in [0,1].
4. the energy saving thrust distribution method of the dynamic positioning according to claim 3 based on improved differential evolution mechanism, special
Sign is:According to the following formula optimizing strategy is selected to adjust λ value:
Wherein t is current cycle time, and T is selected maximum cycle.
5. the energy saving thrust distribution method of the dynamic positioning according to claim 3 based on improved differential evolution mechanism, special
Sign is:Step S5 includes the following steps:It is utilized respectively formulaABC improvement strategies or PSO are improved
Strategy carries out mutation operation;
Wherein ABC improvement strategies include the following steps:
1. bee neighborhood search is led to generate new explanation, and calculate its fitness value;
2. according to formulaGenerate solution to be evaluated;
3. obtain mutant test group;
PSO improvement strategies include the following steps:
1. the history optimal value Pbest and global optimum Gbest of each individual of update;
2. according to formula x 'i=xr1+F(xr2-xr3)+ψ(xgbest-xr1) generation solution to be evaluated;
3. obtain mutant test group.
6. the energy saving thrust distribution method of the dynamic positioning according to claim 1 based on improved differential evolution mechanism, special
Sign is:
Step S6 includes the following steps:According to the cross parameter of equation below design dynamic increment
CR=CRmin+(CRmax-CRmin)×(t/T)2
Wherein CRmaxAnd CRminIt is the maximum value and minimum value of cross parameter;T is current cycle time, and T is that selected maximum is followed
Ring number.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810067356.5A CN108255062B (en) | 2018-01-22 | 2018-01-22 | Power positioning energy-saving thrust distribution method based on improved differential evolution mechanism |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810067356.5A CN108255062B (en) | 2018-01-22 | 2018-01-22 | Power positioning energy-saving thrust distribution method based on improved differential evolution mechanism |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108255062A true CN108255062A (en) | 2018-07-06 |
CN108255062B CN108255062B (en) | 2021-01-05 |
Family
ID=62742888
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810067356.5A Active CN108255062B (en) | 2018-01-22 | 2018-01-22 | Power positioning energy-saving thrust distribution method based on improved differential evolution mechanism |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108255062B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109709970A (en) * | 2018-12-25 | 2019-05-03 | 哈尔滨工程大学 | A kind of underwater robot six degree of freedom thrust distribution optimization method |
CN109765910A (en) * | 2019-03-08 | 2019-05-17 | 哈尔滨工程大学 | A kind of dynamic positioning of vessels thrust distribution method based on improved differential evolution algorithm |
CN111812976A (en) * | 2020-06-06 | 2020-10-23 | 智慧航海(青岛)智能***工程有限公司 | Ship thrust distribution system and thrust distribution method |
CN112394641A (en) * | 2020-11-25 | 2021-02-23 | 哈尔滨工程大学 | Ship real-time thrust distribution method based on double-thread parallel computation |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508431A (en) * | 2011-11-04 | 2012-06-20 | 江苏科技大学 | Thrust distribution method for power positioning system of offshore drilling platform |
CN103092077A (en) * | 2013-01-06 | 2013-05-08 | 中国海洋石油总公司 | Thrust distributing method of dynamic positioning system |
CN103823372A (en) * | 2014-02-24 | 2014-05-28 | 中国船舶重工集团公司第七○二研究所 | Method for distributing thrust of multiple thrusters of ocean engineering equipment dynamic positioning system |
CN105425610A (en) * | 2015-11-09 | 2016-03-23 | 哈尔滨工程大学 | Thrust control simulation system for work-type underwater robot hydraulic propeller |
EP3034411A1 (en) * | 2014-12-17 | 2016-06-22 | Thales | Method for guiding the stationing of a satellite |
CN105911867A (en) * | 2016-06-16 | 2016-08-31 | 哈尔滨工程大学 | Ship thrust distribution method based on NSGA-II algorithm |
RU2015133049A (en) * | 2015-08-07 | 2017-02-09 | Акционерное общество "Информационные спутниковые системы" имени академика М.Ф. Решетнёва" | METHOD FOR DISCOVERING A SPACE VEHICLE INTO A GEOSTATION ORBIT BY USING ELECTRO-REACTIVE ENGINES |
FR3043432A1 (en) * | 2015-11-05 | 2017-05-12 | Snecma | METHOD AND DEVICE FOR ADJUSTING SET VALUES OF A PLURALITY OF VARIABLE GEOMETRIES DETERMINED OF A TURBOMACHINE |
CN106773722A (en) * | 2017-02-17 | 2017-05-31 | 哈尔滨工程大学 | A kind of thrust force distribution method of power location system of ship based on artificial fish-swarm algorithm |
CN106875063A (en) * | 2017-02-21 | 2017-06-20 | 集美大学 | A kind of dynamic positioning ship energy management optimization method |
CN107221964A (en) * | 2017-06-09 | 2017-09-29 | 集美大学 | A kind of dynamic positioning ocean platform multiple generator group scheduling method |
-
2018
- 2018-01-22 CN CN201810067356.5A patent/CN108255062B/en active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508431A (en) * | 2011-11-04 | 2012-06-20 | 江苏科技大学 | Thrust distribution method for power positioning system of offshore drilling platform |
CN103092077A (en) * | 2013-01-06 | 2013-05-08 | 中国海洋石油总公司 | Thrust distributing method of dynamic positioning system |
CN103823372A (en) * | 2014-02-24 | 2014-05-28 | 中国船舶重工集团公司第七○二研究所 | Method for distributing thrust of multiple thrusters of ocean engineering equipment dynamic positioning system |
EP3034411A1 (en) * | 2014-12-17 | 2016-06-22 | Thales | Method for guiding the stationing of a satellite |
RU2015133049A (en) * | 2015-08-07 | 2017-02-09 | Акционерное общество "Информационные спутниковые системы" имени академика М.Ф. Решетнёва" | METHOD FOR DISCOVERING A SPACE VEHICLE INTO A GEOSTATION ORBIT BY USING ELECTRO-REACTIVE ENGINES |
FR3043432A1 (en) * | 2015-11-05 | 2017-05-12 | Snecma | METHOD AND DEVICE FOR ADJUSTING SET VALUES OF A PLURALITY OF VARIABLE GEOMETRIES DETERMINED OF A TURBOMACHINE |
CN105425610A (en) * | 2015-11-09 | 2016-03-23 | 哈尔滨工程大学 | Thrust control simulation system for work-type underwater robot hydraulic propeller |
CN105911867A (en) * | 2016-06-16 | 2016-08-31 | 哈尔滨工程大学 | Ship thrust distribution method based on NSGA-II algorithm |
CN106773722A (en) * | 2017-02-17 | 2017-05-31 | 哈尔滨工程大学 | A kind of thrust force distribution method of power location system of ship based on artificial fish-swarm algorithm |
CN106875063A (en) * | 2017-02-21 | 2017-06-20 | 集美大学 | A kind of dynamic positioning ship energy management optimization method |
CN107221964A (en) * | 2017-06-09 | 2017-09-29 | 集美大学 | A kind of dynamic positioning ocean platform multiple generator group scheduling method |
Non-Patent Citations (3)
Title |
---|
DEFENG WU: "A dynamic positioning thrust allocation approach based on a hybrid artificial colony bee algorithm with chaotic search", 《INTERNATIONAL JOURNAL OF MODELLING, IDENTIFICATION AND CONTROL》 * |
DEFENG WU: "An energy optimal thrust allocation method for the marine dynamic positioning system based on adaptive hybrid artifi cial bee colony algorithm", 《OCEAN ENGINEERING》 * |
刘学军等: "面向船舶动力定位模拟器的推力分配仿真子***开发", 《船舶工程》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109709970A (en) * | 2018-12-25 | 2019-05-03 | 哈尔滨工程大学 | A kind of underwater robot six degree of freedom thrust distribution optimization method |
CN109709970B (en) * | 2018-12-25 | 2022-01-14 | 哈尔滨工程大学 | Underwater robot six-degree-of-freedom thrust distribution optimization method |
CN109765910A (en) * | 2019-03-08 | 2019-05-17 | 哈尔滨工程大学 | A kind of dynamic positioning of vessels thrust distribution method based on improved differential evolution algorithm |
CN111812976A (en) * | 2020-06-06 | 2020-10-23 | 智慧航海(青岛)智能***工程有限公司 | Ship thrust distribution system and thrust distribution method |
CN112394641A (en) * | 2020-11-25 | 2021-02-23 | 哈尔滨工程大学 | Ship real-time thrust distribution method based on double-thread parallel computation |
CN112394641B (en) * | 2020-11-25 | 2023-09-29 | 哈尔滨工程大学 | Ship real-time thrust distribution method based on double-linear parallel computing |
Also Published As
Publication number | Publication date |
---|---|
CN108255062B (en) | 2021-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108255062A (en) | The energy saving thrust distribution method of dynamic positioning based on improved differential evolution mechanism | |
LU102400B1 (en) | Path planning method and system for unmanned surface vehicle based on improved genetic algorithm | |
CN106527133A (en) | Multi-paddle coordination control distribution method for ship | |
CN108490770A (en) | A kind of thrust force distribution method of power location system of ship based on hybrid algorithm | |
CN110109358B (en) | Feedback-based hybrid multi-agent cooperative control method | |
CN110347181A (en) | Unmanned plane distribution formation control method based on energy consumption | |
CN110165714B (en) | Micro-grid integrated scheduling and control method based on extreme dynamic programming algorithm and computer readable storage medium | |
CN110490422A (en) | A kind of target fighting efficiency method for situation assessment based on game cloud model | |
CN108572550A (en) | A kind of online real-time thrust distribution method based on machine learning | |
CN106845725B (en) | Engineering parameter optimization method and system | |
CN108394429B (en) | Method for generating automatic driving curve for urban rail train group | |
CN111459026A (en) | Multi-target control distribution method for tilt rotor unmanned aerial vehicle based on competition group algorithm | |
CN108594793A (en) | A kind of improved RBF flight control systems fault diagnosis network training method | |
CN116859728A (en) | Ship dynamic positioning system thrust distribution method based on improved multi-target particle swarm algorithm | |
CN108537663A (en) | One B shareB trend forecasting method | |
CN110398961B (en) | Hovercraft track keeping rudder propeller cooperative control and distribution method | |
CN115390452B (en) | LQR transverse controller parameter online self-adaption method and system | |
CN107918393B (en) | Marine Autopilot based on depth confidence network | |
CN113126503B (en) | Intelligent ship thrust distribution method considering propeller faults | |
Liu et al. | Ship collision avoidance path planning strategy based on quantum bacterial foraging algorithm | |
CN117170382B (en) | Double unmanned ship cooperative control method suitable for homodromous real-time point location tracking | |
CN109754108A (en) | Unit Economic load distribution method based on fluctuating acceleration coefficient Chaos-Particle Swarm Optimization | |
CN117111620B (en) | Autonomous decision-making method for task allocation of heterogeneous unmanned system | |
Lin et al. | A self-constructing compensatory neural fuzzy system and its applications | |
CN114089633B (en) | Multi-motor coupling driving control device and method for underwater robot |
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 |