0 引言
1 问题描述与模型构建
1.1 问题描述
1.2 模型构建
1.2.1 符号说明
表1 集合、参数和变量定义 |
| 符号 | 定义 | |
|---|---|---|
| 集合 | V | 节点集合,包括仓库节点0和客户节点1~n, |
| K1 | 货车集合 | |
| K2 | 无人机集合 | |
| T | 公交趟次集合 | |
| 参数 | qi1 | 客户的货车直配需求量( |
| qi2 | 客户的无人机直配需求量( | |
| qi3 | 客户的普通配送需求量( | |
| Q1 | 货车的最大装载量 | |
| Q2 | 无人机的最大装载量 | |
| Q3 | 公交的最大装载量 | |
| R | 无人机的最大续航里程 | |
| c1 | 货车的单位距离成本 | |
| c2 | 无人机的单位距离成本 | |
| c3 | 每趟公交的固定成本 | |
| f1 | 货车的固定发车成本 | |
| f2 | 无人机的固定起降成本 | |
| dij | 节点i与节点j之间的距离 | |
| 变量 | | 二元变量,货车k从节点i行驶到节点j时为1, 否则为0( |
| | 二元变量,无人机k从节点i行驶到节点j时为1, 否则为0( | |
| | 二元变量,公交趟次t是否被激活( | |
| | 公交趟次t为客户i配送的普通配送需求量 ( | |
| sik | 货车k为客户i配送的普通配送需求量 ( | |
| | 辅助变量,货车k服务节点i的顺序 | |
| μik | 辅助变量,无人机k服务节点i的顺序 |
1.2.2 模型建立
2 自适应大邻域搜索算法设计
2.1 算法框架
2.2 初始解生成
2.3 破坏算子和修复算子
2.3.1 破坏算子
2.3.2 修复算子
2.4 自适应过程
3 案例分析
3.1 算例仿真
表2 CPLEX求解器结果与ALNS算法求解结果 |
| 算例 节点数 | CPLEX求解结果 | ALNS算法求解结果(10次运行结果) | D1(%) | D2(%) | |||
|---|---|---|---|---|---|---|---|
| 最优解 | 求解时间/s | 最优解 | 平均解 | 平均求解时间/s | |||
| 10 | 832.54 | 361.71 | 832.54 | 832.54 | 6.85 | 0 | 0 |
| 20 | 1 180.52 | 2 782.79 | 1 175.23 | 1 176.08 | 7.14 | -0.45 | 0.07 |
| 30 | — | — | 1 811.46 | 1 832.06 | 7.29 | — | 1.14 |
| 40 | — | — | 2 518.48 | 2 530.28 | 10.16 | — | 0.47 |
| 50 | — | — | 3 098.05 | 3 104.17 | 12.07 | — | 1.98 |
| 60 | — | — | 3 965.34 | 3 974.97 | 13.22 | — | 0.24 |
3.2 案例求解
表3 威坪镇汽车客运站及村级站点的经纬度 |
| 节点索引 | 节点名称 | 经度 | 纬度 |
|---|---|---|---|
| 0 | 威坪汽车客运站 | 118.794 504° | 29.735 831° |
| 1 | 河村村 | 118.850 38° | 29.832 334° |
| 2 | 茶合村 | 118.840 118° | 29.824 674° |
| 3 | 洞源村 | 118.835 394° | 29.837 448° |
| 4 | 妙石村 | 118.832 562° | 29.818 527° |
| 5 | 贤茂村 | 118.827 886° | 29.814 591° |
| 6 | 叶家村 | 118.824 598° | 29.813 383° |
| 7 | 安川村 | 118.809 243° | 29.806 449° |
| 8 | 杨家畈村 | 118.826 213° | 29.804 163° |
| 9 | 岭脚村 | 118.822 269° | 29.791 653° |
| 10 | 洪圻村 | 118.819 461° | 29.781 248° |
| 11 | 邵宅村 | 118.807 635° | 29.779 358° |
| 12 | 三洲村 | 118.812 267° | 29.765 774° |
| 13 | 方宅村 | 118.804 871° | 29.757 763° |
| 14 | 驮岭脚村 | 118.790 176° | 29.741 061° |
表4 不同节点的3种配送需求量 |
| 配送 需求 | 节点 | |||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | |
| q1 | 22 | 30 | 37 | 33 | 0 | 40 | ||||||||
| q2 | 7 | 7 | 10 | 5 | ||||||||||
| q3 | 24 | 13 | 27 | 10 | 21 | 13 | 16 | 10 | 15 | 24 | 26 | 24 | 16 | 20 |
表5 模型参数设置 |
| 参数 | 取值 | 参数 | 取值 | ||
|---|---|---|---|---|---|
| 载 运 工 具 相 关 | 货车最大 载重/kg | 180 | 成 本 相 关 | 货车固定发车 成本/(元/辆) | 100 |
| 无人机最大 载重/kg | 15 | 货车单位距离 成本/(元/km) | 1.5 | ||
| 无人机最大 续航里程/km | 20 | 无人机固定发车 成本/(元/架) | 45 | ||
| 公交最大 载重/kg | 50 | 无人机单位距离 成本/(元/km) | 1.2 | ||
| 公交每天最大 运行趟次/趟 | 2 | 公交趟次固定 成本/(元/趟) | 30 | ||
表6 不同配送模式优化结果对比 |
| 载运工具 | 联合配送 | 货车+无人机 | 货车+公交 | 纯货车 |
|---|---|---|---|---|
| 公交服务 | 趟次1:5-3 趟次2:8-2-1 | 趟次1:5-3 趟次2:8-2-1 | ||
| 货车 | 路径1:0-11-7-6-4-9-0 路径2:0-14-13-10-12-0 | 路径1:0-7-6-5-4-2-1-3-0 路径2:0-12-11-10-9-8-0 路径3:0-14-13-0 | 路径1:0-11-7-6-4-9-0 路径2:0-14-13-10-12-0 路径3:0-4-0 路径4:0-5-0 路径5:0-10-0 路径6:0-11-0 | 路径1:0-7-6-5-4-2-1-3-0 路径2:0-12-11-10-9-8-0 路径3:0-14-13-0 路径4:0-4-0 路径5:0-5-0 路径6:0-10-0 路径7:0-11-0 |
| 无人机 | 路径1:0-11-10-0 路径2:0-5-4-0 | 路径1:0-11-10-0 路径2:0-5-4-0 | ||
| 总成本/元 | 436.22 | 502.12 | 797.77 | 863.67 |