露天矿生产的车辆安排.docx
《露天矿生产的车辆安排.docx》由会员分享，可在线阅读，更多相关《露天矿生产的车辆安排.docx（29页珍藏版）》请在冰点文库上搜索。
露天矿生产的车辆安排
露天矿生产的车辆安排（CMCM2003B）
摘要
本文要解决的问题主要围绕减少实际中的经济成本而展开的。
在经济运作中，减少成本是作为露天矿生产带来利润最大化的有效手段之一。
而合理安排有效路线和车次，成为了解决减少成本问题的关键。
鉴于铲点到卸点线路的复杂性，我们把问题分成两个层次加以解决。
首先我们采用了整体规划的算法，建立了数学模型以求得最小运输量。
其基本思想是提取重要的约束性条件，对于总产量达最小的目标函数进行约束，运用lingo程序求出其最优解，最后得出最小运输量为85628.62吨，且第5、6、7个铲点没有使用。
对于层次二，通过解决层次一所得出具体流量计算卡车在各个路线上一个班次最多可以运行的次数及各路线上需要的卡车数，从而得出所有路线要出动13辆卡车.
问题重述
此题类似与产地与销地的运输的整数规划问题。
10个矿位，5个卸点，运输矿石和岩石两种产品。
此题的重点在于限制条件的提取，由于题中所给条件较多如：
每个铲位至多能安置一台电铲，电铲的平均装车时间为5分钟。
卸货地有卸矿石的矿石漏、2个铁路倒装场和卸岩石的岩石漏、岩场等。
每个卸点都有各自的产量要求，岩石卸点的品味限制都为29.5%
1%。
电铲和卸点都不能同时为两辆及两辆以上卡车服务。
卡车每次都是满载运输，且在一个班次中不存在卡车等待情况，卡车只在开始时点火一次。
卡车的平均卸车时间为3分钟。
一个班次为8小时等等。
总之，限制条件比较隐晦，需要从题中认真挖掘。
本文需要解决的问题是如何设计一个班次使得总运量（吨公里）最小，同时出动最少的卡车，从而运输成本最小。
一个班次的生产计划应该包含以下内容：
出动几台电铲，分别在哪些铲位上；出动几辆卡车，分别在哪些路线上各运输多少次。
因此问题可转化为：
1、如何在满足所有限制条件的条件下，使得一个班次内所有卡车的总运量最小。
2在总运量最小的情况下，如何设计卡车路线，使得出动的卡车数最少，从而使总成本最低。
模型假设
1.卡车每次都是满载，且不出现堵车情况。
2.发动机点火时需要消耗相当多的电瓶能量，故每个班次每台车只在开时点火一次。
3.不存在卡车等待时间
4.在铲位或卸点处因两条路线（及以上）造成的冲突时，只要平均时间能完成任务即可，不进行排时讨论
符号约定
Xij:
从i号铲位到j号卸点的石料运量，单位：
车·次；
Cij：
从i号铲位到j号卸点的距离，单位：
公里；
Tij：
在i号铲位到j号卸点路线上运行一个周期平均所需时间，单位：
分；
Aij：
从i号铲位到j号卸点最多能同时运行的卡车数，单位：
辆；
Bij：
从i号铲位到j号卸点，一辆车一个班次中最多可以运行次数，单位：
次；
pi：
i号铲位的矿石铁含量
：
j号卸点任务需求
：
i号铲位的铁矿石储量，单位：
万吨；
：
i号铲位的岩石储量，单位：
万吨；
：
描述第i号铲位是否使用的0−1变量，
不使用的铲位为0，使用的为1；
：
从i号铲位到j号卸点所需要的车的最小辆数（实数）数值上等于
/
模型建立与求解
问题一：
总运量最小运输路线的设计
目标函数:
由题意易得目标函数为
（1）道路能力约束：
一个电铲（卸点）不能同时为两辆卡车服务，一条路线上最多能同时运行的卡车数是有限制的。
卡车从i号铲位到j号卸点运行一个周期平均所需时间为
（分钟）。
由于装车时间5分钟大于卸车时间3分钟，所以这条路线上在卡车不等待条件下最多能同时运行的卡车数为：
；其中最后开始发车的一辆卡车一个班次中在这条路线上最多可以运行的次数为（其他卡车可能比此数多1次）
，这里
是开始装车时最后一辆车的延时时间。
一个班次中这条固定路线上最多可能运行的总车次大约为：
，总吨数
。
（2）电铲能力约束：
一台电铲不能同时为两辆卡车服务，所以一台电铲在一个班次中的最大可能产量为8×60/5×154（吨）。
（3）卸点能力约束：
卸点的最大吞吐量为每小时60/3=20车次，于是一个卸点在一个班次中的最大可能产量为8×20×154（吨）。
（4）铲位储量约束：
铲位的矿石和岩石产量都不能超过相应的储藏量。
（5）产量任务约束：
各卸点的产量不小于该卸点的任务要求。
（6）铁含量约束：
各矿石卸点的平均品位要求都在指定的范围内。
（7）电铲数量约束：
电铲数量约束无法用普通不等式表达，可以引入10个0—1变量来标志各个铲位是否有产量。
（8）整数约束：
当把问题作为整数规划模型时，流量xij除以154为非负整数。
（9）卡车数量约束：
不超过20辆。
得到的一种模型为
（0）
（1）
（2）
（3）
（4）
（5）
（6）
为整数（7）
（8）
求解模型
由于限制条件较多，需要用快速解法
利用LINGO程序对以上模型求解
sets:
cai/1..10/:
p,cy,ck,f;
xie/1..5/:
q;
link（cai,xie）:
a,b,c,t,x,che;
endsets
data:
v=28;
p=30282932313332313331;
q=1.21.31.31.91.3;
c=5.26001.90005.89000.64004.4200
5.19000.99005.61001.76003.8600
4.21001.90005.61001.27003.7200
4.00001.13004.56001.83003.1600
2.95001.27003.51002.74002.2500
2.74002.25003.65002.60002.8100
2.46001.48002.46004.21000.7800
1.90002.04002.46003.72001.6200
0.64003.09001.06005.05001.2700
1.27003.51000.57006.10000.5000;
cy=1.251.101.351.051.151.351.051.151.351.25;
ck=0.951.051.001.051.101.251.051.301.351.25;
enddata
@for（link:
t=120*c/v+8;a=@floor（t/5）;b=@floor（（485-5*a）/t））;
min=@sum（link:
x*154*c）;
@for（link:
x<=a*b）;
@for（cai（i）:
@sum（xie（j）:
x（i,j））<=f（i）*96）;
@for（xie（j）:
@sum（cai（i）:
x（i,j））<=160）;
@for（cai（i）:
x（i,1）+x（i,2）+x（i,5）<=ck（i）*10000/154）;
@for（cai（i）:
x（i,3）+x（i,4）<=cy（i）*10000/154）;
@for（xie（j）:
@sum（cai（i）:
x（i,j））>=q（j）*10000/154）;
@sum（cai（i）:
x（i,1）*（p（i）-30.5））<=0;
@sum（cai（i）:
x（i,2）*（p（i）-30.5））<=0;
@sum（cai（i）:
x（i,5）*（p（i）-30.5））<=0;
@sum（cai（i）:
x（i,1）*（p（i）-28.5））>=0;
@sum（cai（i）:
x（i,2）*（p（i）-28.5））>=0;
@sum（link:
x/b）<=20;
@sum（cai:
f）<=7;
@for（link:
@gin（x））;
@for（cai:
@bin（f））;
@for（link:
che=x/b）;
end
运行结果最小运输量为85628.62吨且第5、6、7个铲点没有使用
问题二：
在总运量最小的情况下，如何设计路线，使得卡车的使用辆最小
由问题一模型求解结果的以下数据
各个路线上的最佳运输车次：
铲位1
铲位2
铲位3
铲位4
铲位5
铲位6
铲位7
铲位8
铲位9
铲位10
矿石漏
13
54
11
倒装场Ⅰ
42
43
岩场
70
15
岩石漏
81
43
倒装场Ⅱ
13
2
70
各路线上需要的卡车数（实数）：
铲位1
铲位2
铲位3
铲位4
铲位5
铲位6
铲位7
铲位8
铲位9
铲位10
矿石漏
0.867
1.862
0.314
倒装场Ⅰ
1.077
1.162
岩场
1.892
0.326
岩石漏
1.841
1.229
倒装场Ⅱ
0.684
0.1
1.489
所有路线所需卡车数（实数）的和为12.843。
即要最少出动13辆卡车。
易得：
有7辆卡车在一个班次内在固定路线上运输，它们的分布如下
铲位1
铲位2
铲位3
铲位4
铲位5
铲位6
铲位7
铲位8
铲位9
铲位10
矿石漏
1
倒装场I
1
1
岩场
1
岩石漏
1
1
倒装场II
1
另外，有6量车需要联合运输来满足总运输量最小。
联合运输因尽量以V型路线为主，避免Z型路线，过程比较复杂，需要用到组合数等理论，在此就暂不求解。
模型评价
本文根据题中所给条件，建立了整数规划模型，由于题中的限制条件较多，运用LINGO软件快速求解，减少了运算量，方便快捷。
优点：
1.此模型应运面较广，使用类似的产销，集中产品的生产等线性规划问题都可适用。
2、用LINGO软件求解，方便快捷
缺点：
1.模型建立不太完善。
2、没有求出联合派车的具体路线。
模型改进和应用
此模型在现实生活中比较常见，类似的生产规划问题都能改进此模型，快速求出较优生产计划。
附录
LINGO程序运行结果：
Globaloptimalsolutionfoundatiteration:
604
Objectivevalue:
85628.62
VariableValueReducedCost
V28.000000.000000
P
（1）30.000000.000000
P
（2）28.000000.000000
P（3）29.000000.000000
P（4）32.000000.000000
P（5）31.000000.000000
P（6）33.000000.000000
P（7）32.000000.000000
P（8）31.000000.000000
P（9）33.000000.000000
P（10）31.000000.000000
CY
（1）1.2500000.000000
CY
（2）1.1000000.000000
CY（3）1.3500000.000000
CY（4）1.0500000.000000
CY（5）1.1500000.000000
CY（6）1.3500000.000000
CY（7）1.0500000.000000
CY（8）1.1500000.000000
CY（9）1.3500000.000000
CY（10）1.2500000.000000
CK
（1）0.95000000.000000
CK
（2）1.0500000.000000
CK（3）1.0000000.000000
CK（4）1.0500000.000000
CK（5）1.1000000.000000
CK（6）1.2500000.000000
CK（7）1.0500000.000000
CK（8）1.3000000.000000
CK（9）1.3500000.000000
CK（10）1.2500000.000000
F
（1）1.0000000.000000
F
（2）1.0000000.000000
F（3）1.0000000.000000
F（4）1.0000000.000000
F（5）0.0000000.000000
F（6）0.0000000.000000
F（7）0.0000000.000000
F（8）1.0000000.000000
F（9）1.0000000.000000
F（10）1.0000000.000000
Q
（1）1.2000000.000000
Q
（2）1.3000000.000000
Q（3）1.3000000.000000
Q（4）1.9000000.000000
Q（5）1.3000000.000000
A（1,1）6.0000000.000000
A（1,2）3.0000000.000000
A（1,3）6.0000000.000000
A（1,4）2.0000000.000000
A（1,5）5.0000000.000000
A（2,1）6.0000000.000000
A（2,2）2.0000000.000000
A（2,3）6.0000000.000000
A（2,4）3.0000000.000000
A（2,5）4.0000000.000000
A（3,1）5.0000000.000000
A（3,2）3.0000000.000000
A（3,3）6.0000000.000000
A（3,4）2.0000000.000000
A（3,5）4.0000000.000000
A（4,1）5.0000000.000000
A（4,2）2.0000000.000000
A（4,3）5.0000000.000000
A（4,4）3.0000000.000000
A（4,5）4.0000000.000000
A（5,1）4.0000000.000000
A（5,2）2.0000000.000000
A（5,3）4.0000000.000000
A（5,4）3.0000000.000000
A（5,5）3.0000000.000000
A（6,1）3.0000000.000000
A（6,2）3.0000000.000000
A（6,3）4.0000000.000000
A（6,4）3.0000000.000000
A（6,5）4.0000000.000000
A（7,1）3.0000000.000000
A（7,2）2.0000000.000000
A（7,3）3.0000000.000000
A（7,4）5.0000000.000000
A（7,5）2.0000000.000000
A（8,1）3.0000000.000000
A（8,2）3.0000000.000000
A（8,3）3.0000000.000000
A（8,4）4.0000000.000000
A（8,5）2.0000000.000000
A（9,1）2.0000000.000000
A（9,2）4.0000000.000000
A（9,3）2.0000000.000000
A（9,4）5.0000000.000000
A（9,5）2.0000000.000000
A（10,1）2.0000000.000000
A（10,2）4.0000000.000000
A（10,3）2.0000000.000000
A（10,4）6.0000000.000000
A（10,5）2.0000000.000000
B（1,1）14.000000.000000
B（1,2）29.000000.000000
B（1,3）13.000000.000000
B（1,4）44.000000.000000
B（1,5）17.000000.000000
B（2,1）15.000000.000000
B（2,2）38.000000.000000
B（2,3）14.000000.000000
B（2,4）30.000000.000000
B（2,5）18.000000.000000
B（3,1）17.000000.000000
B（3,2）29.000000.000000
B（3,3）14.000000.000000
B（3,4）35.000000.000000
B（3,5）19.000000.000000
B（4,1）18.000000.000000
B（4,2）36.000000.000000
B（4,3）16.000000.000000
B（4,4）29.000000.000000
B（4,5）21.000000.000000
B（5,1）22.000000.000000
B（5,2）35.000000.000000
B（5,3）20.000000.000000
B（5,4）23.000000.000000
B（5,5）26.000000.000000
B（6,1）23.000000.000000
B（6,2）26.000000.000000
B（6,3）19.000000.000000
B（6,4）24.000000.000000
B（6,5）23.000000.000000
B（7,1）25.000000.000000
B（7,2）33.000000.000000
B（7,3）25.000000.000000
B（7,4）17.000000.000000
B（7,5）41.000000.000000
B（8,1）29.000000.000000
B（8,2）28.000000.000000
B（8,3）25.000000.000000
B（8,4）19.000000.000000
B（8,5）31.000000.000000
B（9,1）44.000000.000000
B（9,2）21.000000.000000
B（9,3）37.000000.000000
B（9,4）15.000000.000000
B（9,5）35.000000.000000
B（10,1）35.000000.000000
B（10,2）20.000000.000000
B（10,3）45.000000.000000
B（10,4）13.000000.000000
B（10,5）46.000000.000000
C（1,1）5.2600000.000000
C（1,2）1.9000000.000000
C（1,3）5.8900000.000000
C（1,4）0.64000000.000000
C（1,5）4.4200000.000000
C（2,1）5.1900000.000000
C（2,2）0.99000000.000000
C（2,3）5.6100000.000000
C（2,4）1.7600000.000000
C（2,5）3.8600000.000000
C（3,1）4.2100000.000000
C（3,2）1.9000000.000000
C（3,3）5.6100000.000000
C（3,4）1.2700000.000000
C（3,5）3.7200000.000000
C（4,1）4.0000000.000000
C（4,2）1.1300000.000000
C（4,3）4.5600000.000000
C（4,4）1.8300000.000000
C（4,5）3.1600000.000000
C（5,1）2.9500000.000000
C（5,2）1.2700000.000000
C（5,3）3.5100000.000000
C（5,4）2.7400000.000000
C（5,5）2.2500000.000000
C（6,1）2.7400000.000000
C（6,2）2.2500000.000000
C（6,3）3.6500000.000000
C（6,4）2.6000000.000000
C（6,5）2.8100000.000000
升级会员 