CN111191813A - Vehicle distribution path optimization method based on cargo load and soft time window limitation - Google Patents

Vehicle distribution path optimization method based on cargo load and soft time window limitation Download PDF

Info

Publication number
CN111191813A
CN111191813A CN201910974101.1A CN201910974101A CN111191813A CN 111191813 A CN111191813 A CN 111191813A CN 201910974101 A CN201910974101 A CN 201910974101A CN 111191813 A CN111191813 A CN 111191813A
Authority
CN
China
Prior art keywords
cost
node
vehicle
electric vehicle
distribution
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
Application number
CN201910974101.1A
Other languages
Chinese (zh)
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.)
Chongqing University of Post and Telecommunications
Original Assignee
Chongqing University of Post and Telecommunications
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 Chongqing University of Post and Telecommunications filed Critical Chongqing University of Post and Telecommunications
Priority to CN201910974101.1A priority Critical patent/CN111191813A/en
Publication of CN111191813A publication Critical patent/CN111191813A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G06Q10/047Optimisation of routes or paths, e.g. travelling salesman problem
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/12Computing arrangements based on biological models using genetic models
    • G06N3/126Evolutionary algorithms, e.g. genetic algorithms or genetic programming
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/08Logistics, e.g. warehousing, loading or distribution; Inventory or stock management
    • G06Q10/083Shipping
    • G06Q10/0835Relationships between shipper or supplier and carriers
    • G06Q10/08355Routing methods
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Biophysics (AREA)
  • Strategic Management (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • General Business, Economics & Management (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Tourism & Hospitality (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Development Economics (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Genetics & Genomics (AREA)
  • Game Theory and Decision Science (AREA)
  • Biomedical Technology (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Physiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Electric Propulsion And Braking For Vehicles (AREA)

Abstract

The invention requests to protect a vehicle distribution path optimization method based on cargo load and soft time window limitation, on the basis of the existing electric vehicle path optimization, the influence of the vehicle load on the energy consumption of a battery is mainly considered, soft time window constraint is also considered, an electric vehicle path optimization model with the minimized distribution cost (departure cost, driver cost, charging cost and punishment time window cost) as a target is established, and then a genetic algorithm is adopted for solving to obtain a distribution path optimization scheme. The invention establishes a charging cost function of the electric vehicle based on energy consumption; establishing a soft time window penalty cost function of the electric vehicle; and can formulate the delivery route according to customer's goods weight, demand of charging etc. in the actual conditions to reach the delivery cost minimum.

Description

Vehicle distribution path optimization method based on cargo load and soft time window limitation
Technical Field
The invention belongs to the field of logistics distribution, and also belongs to the field of computer application and the technical field of complex function optimization, and particularly relates to an electric vehicle distribution path optimization method based on cargo load and soft time window limitation.
Background
In recent years, with the rapid development of economy, the consumption level and the purchasing demand of users are increasingly improved, and the logistics activity is increased sharply. However, a series of environmental pollution problems caused by energy consumption, carbon emission and the like of logistics vehicles also appear, and how to apply new technology, new process and new material to reduce the negative effect of logistics activities becomes an important problem in the development of modern logistics industry. The electric automobile has the characteristics of low consumption, low emission and low pollution, and is one of effective ways for solving the problem, so that a series of policy and regulations are successively issued by the country to promote the development of the electric automobile industry and supporting service facilities thereof, and popular leaders in the industries of Jingdong, Alibaba and the like also vigorously implement the electromotization of the automobile. Compared with the traditional fuel vehicle, the electric vehicle has shorter driving distance and needs to go to the charging station for charging irregularly in the distribution process, and because the charging facilities are relatively less at the present stage, the situation that the electric vehicle goes around to the charging station often occurs, so that the normal logistics distribution process is influenced, the time in transit of the distribution vehicle is prolonged, and the logistics service quality is reduced. Therefore, how to reasonably optimize the distribution route of the electric vehicle becomes an important problem to be solved urgently.
The electric vehicle routing problem is developed from the vehicle routing problem. In recent years, most of research on the distribution problem of electric vehicles focuses on the charging technology and the location problem of a charging station, and the research results of the path optimization problem based on electric logistics vehicles are relatively few. The existing electric vehicle path planning research mainly considers constraint conditions such as battery electric quantity and optimal charging path, and constructs a model by considering different influence factors (low carbon, driving range, time window and the like). At present, constraint conditions such as low carbon, driving range, charging strategies and the like are mainly considered in an electric logistics distribution route optimization model, only the influence of a transport distance on electric energy consumption of a vehicle is considered in the aspect of electric energy consumption, although the influence of other factors such as the position of the electric vehicle on the electric energy consumption is rarely considered in the near future, the influence of cargo load on the electric energy consumption of the electric vehicle is not considered, so that the actual electric energy consumption of the vehicle during running deviates from the theory, the condition of electric energy consumption of the vehicle midway occurs, and the scientificity of the electric vehicle distribution route is further influenced. Electric vehicles are also adopted in logistics distribution to provide distribution service for each customer point, and each customer point has a distribution time limit. A certain penalty will be accepted if the delivery fails to reach the client within the time window. The total cost formed by the charging cost, the punishment cost, the departure cost and the driver cost in logistics distribution is minimum, and scientific and reasonable planning needs to be carried out on the vehicle running path in distribution. The invention carries out deep analysis on the researched problems and constructs a mathematical model corresponding to the problems on the premise of certain basic assumption.
Disclosure of Invention
The present invention is directed to solving the above problems of the prior art. The method is capable of making a distribution route according to the weight of the client goods, the charging requirement and the like in actual conditions, and therefore the distribution cost is minimum. The technical scheme of the invention is as follows:
a vehicle distribution path optimization method based on cargo load and soft time window limitation comprises the following steps:
step1, acquiring actual conditions of travel constraint, cargo load and time constraint of an electric vehicle, and constructing an improved electric vehicle path optimization model with minimized distribution cost, wherein the distribution cost comprises vehicle exit cost, driver cost, charging cost and penalty time window cost, an introduced energy consumption formula is introduced into a target function of the improved electric vehicle path optimization model, the influence of a soft time window and vehicle load on battery energy consumption is added into energy consumption, and constraint conditions are set for the target function of the improved electric vehicle path optimization model;
and 2, designing a genetic algorithm coding scheme by adopting a natural number coding mode, fusing the genetic algorithm coding scheme and a target function into a fitness function of the algorithm because the coding scheme of the genetic algorithm cannot completely reflect the constraint condition of the model, and adopting a new improved cross operator to retain excellent parent sub-path information as much as possible in the cross operation and retain and inherit the excellent sub-paths in the parent chromosomes into the child chromosomes in order to improve the optimizing capability of the algorithm in the later evolution stage. And (3) solving an objective function of the improved electric vehicle path optimization model in the step1 to obtain an optimal distribution path.
Further, the objective function of the improved delivery path optimization objective model of the electric vehicle in step1 is
Figure BDA0002233052790000031
Figure BDA0002233052790000032
Shows the vehicle cost C1
Figure BDA0002233052790000033
Indicating driver cost C2
Figure BDA0002233052790000034
Represents a charging cost C3,PC(i)Represents a penalty time window cost C4(ii) a Z represents the total cost of vehicle distribution.
Wherein, V is the set of all nodes, V ═ O ∪ C ∪ F ∪ O ', which is composed of the customer set C, the charging station set F, the distribution center O and the virtual distribution center O' together, K ═ K {1, 2.., K } represents the set of electric vehicles1、C2、C3、C4Respectively representing the use cost, the driver cost, the charging cost and the punishment cost of the vehicle; c0、Cl、CeThe unit vehicle use cost, the unit running cost and the unit electricity price are represented; sikService time, S, of electric vehicle k at customer node iikService time of the electric vehicle k at the charging node i; lijDistance between node i and node j, SikService time of the electric vehicle k at the charging node i;
Tikthe moment when the electric vehicle k reaches the node i; t isi early、Ti delayRepresents the standard earliest, latest service time of client i; m isijRepresenting the energy consumed by the electric vehicle on node i to node j; x is the number ofijk: variable 0-1, xijkA value of 1 indicates that the electric vehicle k directly drives to the point j when reaching the node i, and otherwise, the value is 0; y isik: variable 0-1, yikA value of 1 indicates that the electric vehicle k arrives at the node i to be charged, and otherwise is 0.
The charging cost in the objective function is equal to the cost for compensating the energy consumption of the battery of the electric logistics vehicle, namely, the load of the goods is introduced, and the influence of the load of the vehicle on the energy consumption of the battery is mainly considered;
energy consumption formula: m isij=(αij(ω+fij)lij+βvij 2lij)/efSubstitution into
Figure BDA0002233052790000035
Wherein, αij=a+gsinθij+gCrcosθijIs a specific constant, β ═ 0.5CdA ρ represents a vehicle specific constant, efRepresenting engine efficiency, fijRepresenting the load of the electric vehicle k on the node i and the node j; v. ofijRepresenting the distribution speed of the electric vehicle at the node i and the node j; ω represents vehicle weight at full (tons), a represents acceleration (m/s2), g represents gravitational constant (m/s2), θ represents road angle, A represents frontal area of vehicle (m2), ρ represents air density (kg/m3), and CrDenotes the coefficient of rolling resistance, CdThe rolling resistance coefficient is shown.
The expression for the charging cost function is found as follows:
Figure BDA0002233052790000041
the objective function also takes into account soft time window constraints, assuming that each electric vehicle has a delivery time window Ti early,Ti delay]When the vehicle is delivered in the time window, no penalty cost is paid, if the vehicle isWhen the user arrives early or late, a certain punishment time cost is generated, the punishment time cost is in a linear relation with the time of arriving the client point early or late, and a punishment cost function is obtained as follows:
Figure BDA0002233052790000042
can be expressed as
Figure BDA0002233052790000043
Wherein, CaRepresents the cost per unit time, C, of the early arrivalbRepresents the cost per unit time, T, of late arrivali early、Ti delayIndicating the standard earliest, latest service time for client i.
Further, the objective function setting constraint conditions of the improved electric vehicle path optimization model are as follows:
Figure BDA0002233052790000044
Figure BDA0002233052790000045
Figure BDA0002233052790000046
Figure BDA0002233052790000047
Figure BDA0002233052790000048
Figure BDA0002233052790000049
Figure BDA00022330527900000410
Figure BDA00022330527900000411
Figure BDA00022330527900000412
Figure BDA0002233052790000051
Figure BDA0002233052790000052
Figure BDA0002233052790000053
Figure BDA0002233052790000054
Figure BDA0002233052790000055
in the above constraints:
equations (1) (2) represent the assurance that each subscribing client is accessible;
the formula (3) represents a flow conservation criterion, and the electric vehicle can leave after reaching a certain node;
the formulas (4) and (5) show that the electric vehicle starts from the distribution center and returns to the distribution center after completing the task;
formula (6) represents that the number of electric vehicles for the distribution task is less than or equal to the total number of electric vehicles in the distribution center;
equation (7) represents the load capacity when the electric vehicle leaves the distribution center;
equation (8) represents that the load capacity of the electric vehicle when leaving the distribution center is less than the rated load capacity;
the expressions (9) and (10) are expressed from the access node i to the nodeThe consumption of the remaining battery at j should be reduced by mij
Equation (11) represents a time calculation equation of the next access node j;
equation (12) represents that the remaining capacity of the electric vehicle at each node is greater than 0;
the formulae (13) and (14) represent yik、xijkThe 0-1 variable of (1);
variables and parameter symbols in each formula define:
in the formula, Wmax、WokRespectively representing the rated load capacity of the electric vehicle and the load capacity when leaving the distribution center; diDemand of node i; q: battery capacity of the electric vehicle; y isijk: the remaining capacity of the electric vehicle k from the node i to the node j.
Further, the step2 is to solve the objective function of the electric vehicle path optimization model by using a basic genetic algorithm, and comprises the following steps:
(1) chromosomal coding
The method adopts a natural number coding form, wherein a chromosome consists of N nodes, and then the coding length is N, and the genes are randomly arranged by integers from 1 to N;
(2) population initialization
Randomly generating a random permutation of length N, let Wii+1Represents the load capacity of the vehicle leaving the ith customer point to the (i + 1) th customer point in the chromosome, and judges Wii+1Whether or not it is greater than Q, when Vi-1,i≤Wmax、Wi,i+1>WmaxWhen the gene is inserted, 0 is inserted in front of the i position of the chromosome; restart calculation of W from the position after insertion of 0ii+1Values, repeated until the sequence ends; meanwhile, calculating consumed electric quantity from a first node, judging that when the node i is reached, if the residual electric quantity cannot enable the transport vehicle to reach the next node or cannot reach the nearest charging station from the next node, inserting the number of the charging station nearest to the node behind the node i, wherein the number of the node is N +1, N +2, … and N + F, and circulating to the last node; then, inserting 0 into the head and the tail of the chromosome respectively to finally form an initial chromosome; finally, the above process is repeatedGenerating a plurality of chromosomes to form an initial population;
(3) a fitness function, which is set as a target function CiThe inverse of (d), i.e. the fitness function, is: fi=1/Ci
(4) Selection operation
Selecting individuals with high fitness of chromosomes by adopting a roulette mode, and using the individuals as parents to be inherited into the next generation;
(5) crossover operation
There are many methods for interleaving natural number codes, such as sequential interleaving, circular interleaving, etc. These intersection methods are applicable to the TSP problem and are not well suited to the vehicle path optimization problem for multiple delivery vehicles and multiple sub-routes. In order to improve the optimizing capability of the algorithm in the later period of evolution, a new improved crossover operator is adopted, and the excellent sub-paths of the parent generation are reserved to the greatest extent.
(6) Mutation operation
And carrying out mutation operation on the gene sequence according to a set mutation probability Pm by adopting two-point reciprocity.
(7): and (5) evolution reversion operation. Firstly, generating random natural numbers r1 and r 2; and then all nodes between the two points are numbered in reverse. This operation is only valid if fitness increases, and then proceeds to the next generation until the maximum number of iterations is reached.
Further, the specific steps of the interleaving operation in step (5) include:
step1, selecting two parent chromosomes, and randomly selecting a path on the parent chromosomes;
step2, leading the selected sub-road section;
step3, taking the sub-path A of the parent chromosome 1 as a part of the child chromosome 1, adding the codes which the sub-path A does not have in the parent chromosome 2 to the back of the sub-path A according to the sequence in the parent chromosome 2, and adding the code 0 at the end of the chromosome;
step4, for the offspring chromosome 1, adding 1 code 0 at any one of 7 positions after the path A, and then generating 7 chromosomes, and calculating the fitness of the 7 chromosomes, wherein the maximum fitness value is the offspring chromosome 1; other daughter chromosomes are obtained in the same way.
Further, the mutation operation of step (6): and (3) performing mutation operation on the gene sequence by adopting two-point reciprocity according to a set mutation probability Pm.
The invention has the following advantages and beneficial effects:
the existing research mainly considers the constraint conditions of low carbon, driving range, charging strategy and the like for an electric logistics distribution route optimization model, only considers the influence of the transport distance on the electric energy consumption of a vehicle on the electric energy consumption, although the influence of other factors such as the position of the electric vehicle on the electric energy consumption is rarely considered in the recent research, the influence of cargo load on the electric energy consumption of the electric vehicle is not considered, so that the actual electric energy consumption of the vehicle during running deviates from the theory, the condition of electric energy consumption of the vehicle midway occurs, and the scientificity of the electric vehicle distribution route is further influenced. Meanwhile, most studies rarely consider path optimization under the constraint of a soft time window when studying the electric vehicle path problem, which is not consistent with the actual delivery situation. Therefore, on the basis of the existing electric vehicle path optimization, the influence of the vehicle load on the battery energy consumption is mainly considered, the soft time window constraint is also considered, an electric vehicle path optimization model with the aim of minimizing the distribution cost (including the vehicle-out cost, the driver cost, the charging cost and the punishment time window cost) is established, and then a genetic algorithm is adopted for solving to obtain a distribution path optimization scheme. In the process of solving the electric vehicle path optimization model, the distribution path meeting the actual requirement can be quickly obtained.
Compared with the prior art, the method can realize the path optimization of the electric vehicle considering the cargo load and the soft time window, and has important guiding significance for improving the distribution and transportation efficiency, improving the economic benefit of enterprises and reducing the distribution cost.
Drawings
FIG. 1 is a preferred embodiment electric vehicle distribution route provided by the present invention;
FIG. 2 is a flow chart of the present invention for solving a model using a genetic algorithm;
FIG. 3 is a diagram of the operation of the genetic algorithm crossover strategy of the present invention;
fig. 4 is a flowchart of a vehicle distribution route optimization method based on cargo load and soft time window constraints according to the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be described in detail and clearly with reference to the accompanying drawings. The described embodiments are only some of the embodiments of the present invention.
The technical scheme for solving the technical problems is as follows:
as shown in fig. 4, the method for optimizing the distribution route of the electric vehicle, wherein the step of establishing the optimization model of the route of the electric vehicle with the aim of minimizing the cost comprises the following steps:
(1) in order to minimize distribution cost, various factors such as vehicle capacity, transportation time and transportation cost need to be comprehensively considered.
Specifically, according to the actual conditions of the electric vehicle travel constraint, the cargo load and the time constraint, an electric vehicle path optimization model with minimized distribution cost is constructed, and an objective function is as follows:
Figure BDA0002233052790000081
Figure BDA0002233052790000082
shows the vehicle cost C1
Figure BDA0002233052790000083
Indicating driver cost C2
Figure BDA0002233052790000084
Represents a charging cost C3,PC(i)Represents a penalty time window cost C4(ii) a Z represents the total cost of vehicle distribution.
Wherein V is the set of all nodes, V is O ∪ C ∪ F ∪ O', itThe system consists of a customer set C, a charging station set F, a distribution center O and a virtual distribution center O'. K ═ 1, 2.., K } represents a set of electric vehicles. C1、C2、C3、C4Respectively representing the use cost, the driver cost, the charging cost and the punishment cost of the vehicle; c0、Cl、CeThe unit vehicle use cost, the unit running cost and the unit electricity price are represented; sikService time, S, of electric vehicle k at customer node iikService time of the electric vehicle k at the charging node i; lijDistance between node i and node j, SikService time of the electric vehicle k at the charging node i;
Tikthe moment when the electric vehicle k reaches the node i; t isi early、Ti delayRepresents the standard earliest, latest service time of client i; m isijRepresenting the energy consumed by the electric vehicle on node i to node j; x is the number ofijk: variable 0-1, xijkA value of 1 indicates that the electric vehicle k directly drives to the point j when reaching the node i, and otherwise, the value is 0; y isik: variable 0-1, yikA value of 1 indicates that the electric vehicle k arrives at the node i to be charged, and otherwise is 0.
The charging cost function specifies:
the charging cost is equal to the cost generated by compensating the energy consumed by the battery of the electric logistics vehicle, namely, the cargo load is introduced, and the influence of the vehicle load on the energy consumption of the battery is mainly considered.
Energy consumption formula: m isij=(αij(ω+fij)lij+βvij 2lij)/efSubstitution into
Figure BDA0002233052790000091
Wherein, αij=a+gsinθij+gCrcosθijIs a specific constant, β ═ 0.5CdA ρ represents a vehicle specific constant, efRepresenting engine efficiency, fijRepresenting the load of the electric vehicle k on the node i and the node j; v. ofijIndicating electric vehiclesDelivery rates at node i and node j; ω represents vehicle weight at full (tons), a represents acceleration (m/s2), g represents gravitational constant (m/s2), θ represents road angle, A represents frontal area of vehicle (m2), ρ represents air density (kg/m3), and CrDenotes the coefficient of rolling resistance, CdThe rolling resistance coefficient is shown.
The expression for the charging cost function is found as follows:
Figure BDA0002233052790000092
the penalty function specifies:
suppose that each electric vehicle has a delivery time window Ti early,Ti delay]When the vehicle is delivered in the time window, no penalty cost is paid, if the vehicle arrives early or late, a certain penalty time cost is generated, and the magnitude of the penalty time cost is in a linear relation with the time of arriving at the customer point early or late. The penalty cost function is obtained as:
Figure BDA0002233052790000101
can be expressed as
Figure BDA0002233052790000102
Wherein, CaRepresents the cost per unit time, C, of the early arrivalbRepresents the cost per unit time, T, of late arrivali early、Ti delayIndicating the standard earliest, latest service time for client i.
(2) Constraint conditions
s.t
Figure BDA0002233052790000103
Figure BDA0002233052790000104
Figure BDA0002233052790000105
Figure BDA0002233052790000106
Figure BDA0002233052790000107
Figure BDA0002233052790000108
Figure BDA0002233052790000109
Figure BDA00022330527900001010
Figure BDA00022330527900001011
Figure BDA00022330527900001012
Figure BDA00022330527900001013
Figure BDA00022330527900001014
Figure BDA0002233052790000111
Figure BDA0002233052790000112
The above constraints are illustrated as follows:
equations (1) (2) represent the assurance that each subscribing client is accessible;
the formula (3) represents a flow conservation criterion, and the electric vehicle can leave after reaching a certain node;
the formulas (4) and (5) show that the electric vehicle starts from the distribution center and returns to the distribution center after completing the task;
formula (6) represents that the number of electric vehicles for the distribution task is less than or equal to the total number of electric vehicles in the distribution center;
equation (7) represents the load capacity when the electric vehicle leaves the distribution center;
equation (8) represents that the load capacity of the electric vehicle when leaving the distribution center is less than the rated load capacity;
the expressions (9) and (10) indicate that the consumption of the remaining battery should be reduced by m when accessing the nodes i to jij
Equation (11) represents a time calculation equation of the next access node j;
equation (12) represents that the remaining capacity of the electric vehicle at each node is greater than 0;
the formulae (13) and (14) represent yik、xijkThe 0-1 variable of (a).
In the formula, Wmax、WokRespectively representing the rated load capacity of the electric vehicle and the load capacity when leaving the distribution center; diDemand of node i; q: battery capacity of the electric vehicle; y isijk: the remaining capacity of the electric vehicle k from the node i to the node j.
2. Solving is carried out by using a genetic algorithm, and referring to fig. 2, the basic design idea is as follows:
(1) chromosomal coding
In the form of natural number coding, where a chromosome is composed of N nodes, the coding length is N, and the genes are randomly arranged in integers from 1 to N, e.g., 8 demand points, (23156487) is a chromosome.
(2) Population initialization
Randomly generating a random permutation of length N, let Wii+1Indicating that the vehicle leaves the ith customer site in the chromosome to go to the (i + 1) th customerWeight of load at the time of point, Wii+1Whether or not it is greater than Q, when W isi-1,i≤Wmax、Wi,i+1>WmaxWhen the gene is inserted, 0 is inserted in front of the i position of the chromosome; restart calculation of W from the position after insertion of 0ii+1Values, repeated until the sequence ends; meanwhile, calculating consumed electric quantity from a first node, judging that when the node i is reached, if the residual electric quantity cannot enable the transport vehicle to reach the next node or cannot reach the nearest charging station from the next node, inserting the number of the charging station nearest to the node behind the node i, wherein the number of the node is N +1, N +2, … and N + F, and circulating to the last node; then, a 0 is inserted into the head and the tail of the chromosome respectively, and finally an initial chromosome is formed. For example, if there are 8 demand points, 3 electric vehicles, visit 2 charging stations during distribution, then (02358940171060) is a legal chromosome; finally, the above process is repeated to produce several chromosomes constituting the initial population.
(3) A fitness function. The objective function in this document is to solve the minimum problem, so the fitness function in this document can be set as the objective function CiThe inverse of (d), i.e. the fitness function, is: fi=1/Ci
(4) Selection operation
Individuals with high fitness chromosomes are selected by roulette and inherited as parents into the next generation.
(5) Crossover operation
In order to reserve the excellent sub-paths of the parent, the crossover operator is improved, and the specific process is as shown in fig. 3:
step1, selecting two parent chromosomes, and randomly selecting a path on the parent chromosomes;
step2, leading the selected sub-road section;
step3, taking the child path A of the parent chromosome 1 as a part of the child chromosome 1, adding the codes which the child path A does not have in the parent chromosome 2 to the back of the child path A according to the sequence in the parent chromosome 2, and adding the code 0 at the end of the chromosome;
step4, for the offspring chromosome 1, adding 1 code 0 at any one of 7 positions after the path A, and then generating 7 chromosomes, and calculating the fitness of the 7 chromosomes, wherein the maximum fitness value is the offspring chromosome 1; other daughter chromosomes are obtained in the same way.
(6) Mutation operation
And carrying out mutation operation on the gene sequence according to a set mutation probability Pm by adopting two-point reciprocity.
(7) Operation of evolution reversion
Firstly, generating random natural numbers r1 and r 2; and then all node numbers between the two points are reversed. This operation is only valid if fitness increases, and then proceeds to the next generation until the maximum number of iterations is reached.
In summary, the invention constructs a path optimization model of the electric vehicle based on the general path optimization of the electric vehicle, considering the soft time window constraint and the influence of the vehicle load on the battery energy consumption. In the running process of the vehicle, the load of the goods can influence the running mileage of the electric vehicle and further influence the whole distribution path, so that the model introduces an energy consumption formula into an objective function, the influence of a soft time window and the load of the vehicle on the energy consumption of the battery is mainly considered, and finally, a reasonable running route of the distribution vehicle is designed, so that the whole distribution cost is minimum.
The above examples are to be construed as merely illustrative and not limitative of the remainder of the disclosure. After reading the description of the invention, the skilled person can make various changes or modifications to the invention, and these equivalent changes and modifications also fall into the scope of the invention defined by the claims.

Claims (6)

1. A vehicle distribution path optimization method based on cargo load and soft time window limitation is characterized by comprising the following steps:
step1, acquiring actual conditions of travel constraint, cargo load and time constraint of an electric vehicle, and constructing an improved electric vehicle path optimization model with minimized distribution cost, wherein the distribution cost comprises departure cost, driver cost, charging cost and penalty time window cost, an introduced energy consumption formula is introduced into an objective function of the improved electric vehicle path optimization model, the influence of a soft time window and vehicle load on battery energy consumption is added into energy consumption, and constraint conditions are set for the objective function of the improved electric vehicle path optimization model;
and 2, designing a genetic algorithm coding scheme by adopting a natural number coding mode, fusing the genetic algorithm coding scheme and an objective function into a fitness function of the algorithm because the genetic algorithm coding scheme cannot completely reflect the constraint conditions of the model, and adopting a new improved crossover operator to retain excellent parent sub-path information as much as possible in crossover operation and retain and inherit the excellent sub-paths in parent chromosomes into the child chromosomes in order to improve the optimizing capability of the algorithm in the later evolution stage. And (3) solving an objective function of the improved electric vehicle path optimization model in the step1 by using an improved genetic algorithm to obtain an optimal distribution path.
2. The method as claimed in claim 1, wherein the objective function of the step1 is an objective function of the improved electric vehicle distribution route optimization objective model
Figure FDA0002233052780000011
Figure FDA0002233052780000012
Shows the vehicle cost C1
Figure FDA0002233052780000013
Indicating driver cost C2
Figure FDA0002233052780000014
Represents a charging cost C3,PC(i)Represents a penalty time window cost C4(ii) a Z denotes vehicle deliveryThe total cost.
Wherein, V is the set of all nodes, V ═ O ∪ C ∪ F ∪ O ', which is composed of the customer set C, the charging station set F, the distribution center O and the virtual distribution center O' together, K ═ {1, 2.., K } represents the set of electric vehicles1、C2、C3、C4Respectively representing the use cost, the driver cost, the charging cost and the punishment cost of the vehicle; c0、Cl、CeThe unit vehicle use cost, the unit running cost and the unit electricity price are represented; sikService time, S, of electric vehicle k at customer node iikService time of the electric vehicle k at the charging node i; lijDistance between node i and node j, SikService time of the electric vehicle k at the charging node i;
Tikthe moment when the electric vehicle k reaches the node i; t isi early、Ti delayRepresents the standard earliest, latest service time of client i; m isijRepresenting the energy consumed by the electric vehicle on node i to node j; x is the number ofijk: variable 0-1, xijkA value of 1 indicates that the electric vehicle k directly drives to the point j when reaching the node i, and otherwise, the value is 0; y isik: variable 0-1, yikA value of 1 indicates that the electric vehicle k arrives at the node i to be charged, and otherwise is 0.
The charging cost in the objective function is equal to the cost for compensating the energy consumption of the battery of the electric logistics vehicle, namely, the load of the goods is introduced, and the influence of the load of the vehicle on the energy consumption of the battery is mainly considered;
energy consumption formula: m isij=(αij(ω+fij)lij+βvij 2lij)/efSubstitution into
Figure FDA0002233052780000021
Wherein, αij=a+gsinθij+gCrcosθijIs a specific constant, β ═ 0.5CdA ρ represents a vehicle specific constant, efRepresenting engine efficiency, fijIndicating that the electric vehicle k is at the node i and the nodeThe load at point j; v. ofijRepresenting the distribution speed of the electric vehicle at the node i and the node j; ω represents vehicle weight at full (tons), a represents acceleration (m/s2), g represents gravitational constant (m/s2), θ represents road angle, A represents frontal area of vehicle (m2), ρ represents air density (kg/m3), and CrDenotes the coefficient of rolling resistance, CdRepresents a rolling resistance coefficient;
the expression for the charging cost function is found as follows:
Figure FDA0002233052780000022
the objective function also takes into account soft time window constraints, assuming that each electric vehicle has a delivery time window Ti early,Ti delay]When the vehicle delivers goods in the time window, no penalty cost is paid, if the vehicle arrives early or late, a certain penalty time cost is generated, the magnitude of the penalty time cost is in a linear relation with the time of arriving at the customer point early or late, and the penalty cost function is obtained as follows:
Figure FDA0002233052780000023
can be expressed as
Figure FDA0002233052780000024
Wherein, CaRepresents the cost per unit time, C, of the early arrivalbRepresents the cost per unit time, T, of late arrivali early、Ti delayIndicating the standard earliest, latest service time for client i.
3. The method of claim 2, wherein the objective function of the improved electric vehicle path optimization model sets constraints as follows:
Figure FDA0002233052780000031
Figure FDA0002233052780000032
Figure FDA0002233052780000033
Figure FDA0002233052780000034
Figure FDA0002233052780000035
Figure FDA0002233052780000036
Figure FDA0002233052780000037
Figure FDA0002233052780000038
Figure FDA0002233052780000039
Figure FDA00022330527800000310
Figure FDA00022330527800000311
Figure FDA00022330527800000312
Figure FDA00022330527800000313
Figure FDA00022330527800000314
in the above constraints:
equations (1) (2) represent the assurance that each subscribing client is accessible;
the formula (3) represents a flow conservation criterion, and the electric vehicle can leave after reaching a certain node;
the formulas (4) and (5) show that the electric vehicle starts from the distribution center and returns to the distribution center after completing the task;
formula (6) represents that the number of electric vehicles for the distribution task is less than or equal to the total number of electric vehicles in the distribution center;
equation (7) represents the load capacity when the electric vehicle leaves the distribution center;
equation (8) represents that the load capacity of the electric vehicle when leaving the distribution center is less than the rated load capacity;
the expressions (9) and (10) indicate that the consumption of the remaining battery should be reduced by m when accessing the nodes i to jij
Equation (11) represents a time calculation equation of the next access node j;
equation (12) represents that the remaining capacity of the electric vehicle at each node is greater than 0;
the formulae (13) and (14) represent yik、xijkThe 0-1 variable of (1);
in the formula, Wmax、WokRespectively representing the rated load capacity of the electric vehicle and the load capacity when leaving the distribution center; diDemand of node i; q: battery capacity of the electric vehicle; y isijk: the remaining capacity of the electric vehicle k from the node i to the node j.
4. The method for optimizing vehicle distribution path based on cargo load and soft time window limit as claimed in claim 3, wherein said step2 is to solve the objective function of said electric vehicle path optimization model by using basic genetic algorithm, comprising the following steps:
(1) chromosomal coding
The method adopts a natural number coding form, wherein a chromosome consists of N nodes, the coding length is N, and the genes are randomly arranged by integers from 1 to N;
(2) population initialization
Randomly generating a random permutation of length N, let Wii+1Represents the load capacity of the vehicle leaving the ith customer point to the (i + 1) th customer point in the chromosome, and judges Wii+1Whether or not it is greater than Q, when W isi-1,i≤Wmax、Wi,i+1>WmaxWhen the gene is inserted, 0 is inserted in front of the i position of the chromosome; restart calculation of W from the position after insertion of 0ii+1Values, repeated until the sequence ends; meanwhile, calculating consumed electric quantity from a first node, judging that when the node i is reached, if the residual electric quantity cannot enable the transport vehicle to reach the next node or cannot reach the nearest charging station from the next node, inserting the number of the charging station nearest to the node behind the node i, wherein the number of the node is N +1, N +2, … and N + F, and circulating to the last node; then, inserting 0 into the head and the tail of the chromosome respectively to finally form an initial chromosome; finally, repeating the above process to generate a plurality of chromosomes to form an initial population;
(3) a fitness function, which is set as a target function CiThe inverse of (d), i.e. the fitness function, is: fi=1/Ci
(4) Selection operation
Selecting individuals with high fitness of chromosomes by adopting a roulette mode, and using the individuals as parents to be inherited into the next generation;
(5) crossover operation
There are many methods for interleaving natural number codes, such as sequential interleaving, circular interleaving, etc. These intersection methods are applicable to the TSP problem and are not well suited to the vehicle path optimization problem for multiple delivery vehicles, multiple sub-routes. In order to improve the optimizing capability of the algorithm in the later period of evolution, a new improved crossover operator is adopted, and the excellent sub-paths of the parent generation are reserved to the greatest extent.
(6) Mutation operation
And carrying out mutation operation on the gene sequence according to a set mutation probability Pm by adopting two-point reciprocity.
(7): and (5) evolution reversion operation. Firstly, generating random natural numbers r1 and r 2; and then all node numbers between the two points are reversed. This operation is only valid if fitness increases, and then proceeds to the next generation until the maximum number of iterations is reached.
5. The method for optimizing a vehicle distribution path based on cargo load and soft time window limit as claimed in claim 4, wherein the step (5) of interleaving comprises the specific steps of:
step1, selecting two parent chromosomes, and randomly selecting a path on the parent chromosomes;
step2, leading the selected sub-road section;
step3, taking the sub-path A of the parent chromosome 1 as a part of the child chromosome 1, adding the codes which the sub-path A does not have in the parent chromosome 2 to the back of the sub-path A according to the sequence in the parent chromosome 2, and adding the code 0 at the end of the chromosome;
step4, for the offspring chromosome 1, adding 1 code 0 at any one of 7 positions after the path A, and then generating 7 chromosomes, and calculating the fitness of the 7 chromosomes, wherein the maximum fitness value is the offspring chromosome 1; other daughter chromosomes are obtained in the same way.
6. The cargo load and soft time window limit-based vehicle distribution path optimization method of claim 4, wherein the (6) mutation operation: and (3) performing mutation operation on the gene sequence by adopting two-point reciprocity according to a set mutation probability Pm.
CN201910974101.1A 2019-10-14 2019-10-14 Vehicle distribution path optimization method based on cargo load and soft time window limitation Pending CN111191813A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910974101.1A CN111191813A (en) 2019-10-14 2019-10-14 Vehicle distribution path optimization method based on cargo load and soft time window limitation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910974101.1A CN111191813A (en) 2019-10-14 2019-10-14 Vehicle distribution path optimization method based on cargo load and soft time window limitation

Publications (1)

Publication Number Publication Date
CN111191813A true CN111191813A (en) 2020-05-22

Family

ID=70707201

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910974101.1A Pending CN111191813A (en) 2019-10-14 2019-10-14 Vehicle distribution path optimization method based on cargo load and soft time window limitation

Country Status (1)

Country Link
CN (1) CN111191813A (en)

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111966111A (en) * 2020-10-23 2020-11-20 北京国新智电新能源科技有限责任公司 Automatic power distribution based mobile charging equipment formation control method, system and device
CN111985676A (en) * 2020-06-28 2020-11-24 济南浪潮高新科技投资发展有限公司 Method and equipment for planning transportation line of electric truck
CN112200367A (en) * 2020-10-09 2021-01-08 河北工业大学 Electric vehicle distribution path optimization method supporting charge-discharge strategy
CN112270500A (en) * 2020-11-17 2021-01-26 深圳市兆航物流有限公司 Intelligent supply chain logistics scheduling method and system
CN112488358A (en) * 2020-10-31 2021-03-12 海南电网有限责任公司 Electric vehicle charging path planning method and storage medium
CN112508478A (en) * 2020-12-02 2021-03-16 北京理工大学 Flexible logistics distribution task allocation method based on self-organizing automated guided vehicle
CN112836892A (en) * 2021-02-26 2021-05-25 山东科技大学 Multi-target vehicle distribution path determining method and system based on improved genetic algorithm
CN113029144A (en) * 2021-03-01 2021-06-25 汇链通供应链科技(上海)有限公司 Sub-heuristic algorithm path planning method for collaborative transportation
CN113077106A (en) * 2021-04-16 2021-07-06 北京京东振世信息技术有限公司 Time window-based article transportation method and device
CN113222275A (en) * 2021-05-26 2021-08-06 大连海事大学 Vehicle path optimization method considering space-time distance under time-varying road network
CN113222272A (en) * 2021-05-26 2021-08-06 合肥工业大学 Emergency material transportation and loading cooperative optimization method based on double-layer genetic coding
CN113313285A (en) * 2021-04-21 2021-08-27 山东师范大学 Multi-constraint vehicle path optimization method, system, storage medium and equipment
CN113379115A (en) * 2021-06-04 2021-09-10 大连海事大学 Time-dependent green vehicle path optimization method with time window under fuzzy requirement
CN113780676A (en) * 2021-09-23 2021-12-10 河南科技大学 Method for optimizing distribution path of bottled liquefied gas vehicle
CN114118621A (en) * 2021-12-07 2022-03-01 东华大学 Multi-objective low-carbon logistics scheduling optimization method based on improved Knea
CN114186758A (en) * 2022-02-15 2022-03-15 杭州杰牌传动科技有限公司 Cost-optimal-oriented in-plant logistics distribution method
CN114254822A (en) * 2021-12-19 2022-03-29 浙江工业大学 Unmanned aerial vehicle distribution network optimization model based on Internet of things technology and solution algorithm thereof
CN115660762A (en) * 2022-12-29 2023-01-31 深圳市普拉托科技有限公司 Tray renting method and device and electronic equipment
CN115829170A (en) * 2023-02-17 2023-03-21 鱼快创领智能科技(南京)有限公司 Driving scheme optimization method, system and storage medium
CN116703291A (en) * 2023-06-15 2023-09-05 北京化工大学 Mixed energy vehicle team delivery path optimization method

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102884401A (en) * 2010-05-06 2013-01-16 莱卡地球系统公开股份有限公司 Method and guidance-unit for guiding battery-operated transportation means to reconditioning stations
CN104809549A (en) * 2015-04-02 2015-07-29 常州奥迈信息技术有限公司 Scheduling method of goods vehicle planned driving lines
CN106403968A (en) * 2016-06-06 2017-02-15 四川大学 Planning method for charging of wireless rechargeable sensor networks (WRSNs) with heterogeneous mobile charging vehicles
CN106503836A (en) * 2016-10-09 2017-03-15 电子科技大学 A kind of pure electric automobile logistics distribution Optimization Scheduling of multiple-objection optimization
CN107145971A (en) * 2017-04-18 2017-09-08 苏州工业职业技术学院 A kind of express delivery dispatching optimization method of dynamic adjustment
CN109211257A (en) * 2018-08-31 2019-01-15 北京图森未来科技有限公司 A kind of automatic driving car travel route method and system for planning
CN109583770A (en) * 2018-11-28 2019-04-05 清华四川能源互联网研究院 Vehicle dispatching method and device
CN110059934A (en) * 2019-03-27 2019-07-26 浙江工商大学 The method of fuel vehicle and the scheduling of new energy vehicle coperating distribution

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102884401A (en) * 2010-05-06 2013-01-16 莱卡地球系统公开股份有限公司 Method and guidance-unit for guiding battery-operated transportation means to reconditioning stations
CN104809549A (en) * 2015-04-02 2015-07-29 常州奥迈信息技术有限公司 Scheduling method of goods vehicle planned driving lines
CN106403968A (en) * 2016-06-06 2017-02-15 四川大学 Planning method for charging of wireless rechargeable sensor networks (WRSNs) with heterogeneous mobile charging vehicles
CN106503836A (en) * 2016-10-09 2017-03-15 电子科技大学 A kind of pure electric automobile logistics distribution Optimization Scheduling of multiple-objection optimization
CN107145971A (en) * 2017-04-18 2017-09-08 苏州工业职业技术学院 A kind of express delivery dispatching optimization method of dynamic adjustment
CN109211257A (en) * 2018-08-31 2019-01-15 北京图森未来科技有限公司 A kind of automatic driving car travel route method and system for planning
CN109583770A (en) * 2018-11-28 2019-04-05 清华四川能源互联网研究院 Vehicle dispatching method and device
CN110059934A (en) * 2019-03-27 2019-07-26 浙江工商大学 The method of fuel vehicle and the scheduling of new energy vehicle coperating distribution

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李敏: "基于差异化充电的城市电动物流车配送优化研究" *
沈续昌: "考虑货物重量的电动物流车路径规划研究" *

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111985676A (en) * 2020-06-28 2020-11-24 济南浪潮高新科技投资发展有限公司 Method and equipment for planning transportation line of electric truck
CN112200367A (en) * 2020-10-09 2021-01-08 河北工业大学 Electric vehicle distribution path optimization method supporting charge-discharge strategy
CN111966111A (en) * 2020-10-23 2020-11-20 北京国新智电新能源科技有限责任公司 Automatic power distribution based mobile charging equipment formation control method, system and device
CN112488358A (en) * 2020-10-31 2021-03-12 海南电网有限责任公司 Electric vehicle charging path planning method and storage medium
CN112270500A (en) * 2020-11-17 2021-01-26 深圳市兆航物流有限公司 Intelligent supply chain logistics scheduling method and system
CN112508478A (en) * 2020-12-02 2021-03-16 北京理工大学 Flexible logistics distribution task allocation method based on self-organizing automated guided vehicle
CN112836892A (en) * 2021-02-26 2021-05-25 山东科技大学 Multi-target vehicle distribution path determining method and system based on improved genetic algorithm
CN113029144A (en) * 2021-03-01 2021-06-25 汇链通供应链科技(上海)有限公司 Sub-heuristic algorithm path planning method for collaborative transportation
CN113077106A (en) * 2021-04-16 2021-07-06 北京京东振世信息技术有限公司 Time window-based article transportation method and device
CN113313285A (en) * 2021-04-21 2021-08-27 山东师范大学 Multi-constraint vehicle path optimization method, system, storage medium and equipment
CN113222275A (en) * 2021-05-26 2021-08-06 大连海事大学 Vehicle path optimization method considering space-time distance under time-varying road network
CN113222272A (en) * 2021-05-26 2021-08-06 合肥工业大学 Emergency material transportation and loading cooperative optimization method based on double-layer genetic coding
CN113222272B (en) * 2021-05-26 2022-09-20 合肥工业大学 Emergency material transportation and loading cooperative optimization method based on double-layer genetic coding
CN113379115A (en) * 2021-06-04 2021-09-10 大连海事大学 Time-dependent green vehicle path optimization method with time window under fuzzy requirement
CN113780676A (en) * 2021-09-23 2021-12-10 河南科技大学 Method for optimizing distribution path of bottled liquefied gas vehicle
CN113780676B (en) * 2021-09-23 2023-06-23 河南科技大学 Method for optimizing distribution path of bottled liquefied gas vehicle
CN114118621A (en) * 2021-12-07 2022-03-01 东华大学 Multi-objective low-carbon logistics scheduling optimization method based on improved Knea
CN114118621B (en) * 2021-12-07 2024-04-23 东华大学 Optimization method for multi-target low-carbon logistics scheduling based on improvement Knea
CN114254822A (en) * 2021-12-19 2022-03-29 浙江工业大学 Unmanned aerial vehicle distribution network optimization model based on Internet of things technology and solution algorithm thereof
CN114254822B (en) * 2021-12-19 2024-05-03 浙江工业大学 Unmanned aerial vehicle distribution network optimization model based on Internet of things technology and solving algorithm thereof
CN114186758A (en) * 2022-02-15 2022-03-15 杭州杰牌传动科技有限公司 Cost-optimal-oriented in-plant logistics distribution method
CN115660762A (en) * 2022-12-29 2023-01-31 深圳市普拉托科技有限公司 Tray renting method and device and electronic equipment
CN115829170A (en) * 2023-02-17 2023-03-21 鱼快创领智能科技(南京)有限公司 Driving scheme optimization method, system and storage medium
CN116703291A (en) * 2023-06-15 2023-09-05 北京化工大学 Mixed energy vehicle team delivery path optimization method
CN116703291B (en) * 2023-06-15 2024-01-05 北京化工大学 Mixed energy vehicle team delivery path optimization method

Similar Documents

Publication Publication Date Title
CN111191813A (en) Vehicle distribution path optimization method based on cargo load and soft time window limitation
Shao et al. Electric vehicle‐routing problem with charging demands and energy consumption
Liao et al. Multi-objective green meal delivery routing problem based on a two-stage solution strategy
Keskin et al. Partial recharge strategies for the electric vehicle routing problem with time windows
CN112200367B (en) Electric vehicle distribution path optimization method supporting charge-discharge strategy
Li et al. Public charging station localization and route planning of electric vehicles considering the operational strategy: A bi-level optimizing approach
Zhenfeng et al. The electric vehicle routing problem with time windows using genetic algorithm
CN110472792B (en) Logistics distribution vehicle route optimization method based on discrete bat algorithm
CN112836892A (en) Multi-target vehicle distribution path determining method and system based on improved genetic algorithm
CN111222705B (en) Nonlinear charging vehicle path optimization method
CN115953104B (en) Hybrid fleet scheduling method based on dung beetle optimization algorithm
CN113919557A (en) Logistics route optimization method and system based on self-adaptive NSGAII
CN113822461A (en) Track traffic cross-line operation optimization method, system, equipment and storage medium
CN114444843A (en) Agricultural product green logistics distribution vehicle scheduling method and system based on large-scale variable neighborhood search strategy
CN113780961B (en) Low-carbon vaccine cold chain optimization distribution method of multi-target firework algorithm
Zhao et al. Research on emergency distribution optimization of mobile power for electric vehicle in photovoltaic-energy storage-charging supply chain under the energy blockchain
Zhao et al. Bi-objective optimization for vehicle routing problems with a mixed fleet of conventional and electric vehicles and soft time windows
CN110956317A (en) Modeling and solving method for unmanned vehicle and fuel vehicle mixed distribution mode
Wu et al. Multi-objective reinforcement learning-based energy management for fuel cell vehicles considering lifecycle costs
CN113887782A (en) Genetic-firework mixing method and system for maintenance resource distribution scheduling
CN116358593B (en) Electric vehicle path planning method, device and equipment considering nonlinear energy consumption
Wu et al. A neighborhood comprehensive learning particle swarm optimization for the vehicle routing problem with time windows
Boukhater et al. An intelligent and fair GA carpooling scheduler as a social solution for greener transportation
Park et al. Multiobjective approach to the transit network design problem with variable demand considering transit equity
CN111639822B (en) Express distribution method based on 0-1 knapsack problem analysis

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20200522

RJ01 Rejection of invention patent application after publication