数学建模钢管订购和运输.doc
- 文档编号:153258
- 上传时间:2022-10-04
- 格式:DOC
- 页数:26
- 大小:482KB
数学建模钢管订购和运输.doc
《数学建模钢管订购和运输.doc》由会员分享,可在线阅读,更多相关《数学建模钢管订购和运输.doc(26页珍藏版)》请在冰豆网上搜索。
钢管的订购和运输优化模型
摘要
本文建立的多元非线性优化模型。
问题一在保证天然气管道铺设可以顺利实施的情况下,给出了钢管的订购与运输总费用最小的方案。
在求钢管由钢厂运输到站点的费用和铺设钢管时产生的运输费,根据图一,我们通过深度优先遍历的方法对整个图一进行路径搜索,然后根据每条搜索到的路径上的铁路和公路上的不同权重,找到了各个钢厂到各个天然气管道上的站点的最佳路径。
对于整个优化过程我们给出了相关的算法,并用matlab软件编程,经过一系列计算之后,得出了最优的订购与运输方案。
对于问题1,我们求得的最优解为(具体方案见表五):
总费用
800
800
1000
0
1190
1181
0
对于问题2我们经过计算比较得出:
钢管销价的变化对购运计划和总费用影响最大。
的生产上限的变化购运计划和总费用影响最大。
对于问题3,当天然气管道呈现的是一个树状图的时候,我们得到的最优解为(具体方案见表六):
总费用
800
800
1000
0
1450
1853
0
关键字:
非线性优化深度优先遍历最佳路径
一、问题重述
要铺设一条的输送天然气的主管道,如图一所示(见下页)。
经筛选后可以生产这种主管道钢管的钢厂有。
图中粗线表示铁路,单细线表示公路,双细线表示要铺设的管道(假设沿管道或者原来有公路,或者建有施工公路),圆圈表示火车站,每段铁路、公路和管道旁的阿拉伯数字表示里程(单位km)。
为方便计,1km主管道钢管称为1单位钢管。
一个钢厂如果承担制造这种钢管,至少需要生产500个单位。
钢厂在指定期限内能生产该钢管的最大数量为个单位,钢管出厂销价1单位钢管为万元,如下表:
1
2
3
4
5
6
7
800
800
1000
2000
2000
2000
3000
160
155
155
160
155
150
160
1单位钢管的铁路运价如下表:
里程(km)
≤300
301~350
351~400
401~450
451~500
运价(万元)
20
23
26
29
32
里程(km)
501~600
601~700
701~800
801~900
901~1000
运价(万元)
37
44
50
55
60
1000km以上每增加1至100km运价增加5万元。
公路运输费用为1单位钢管每公里0.1万元(不足整公里部分按整公里计算)。
钢管可由铁路、公路运往铺设地点(不只是运到点,而是管道全线)。
(1)请制定一个主管道钢管的订购和运输计划,使总费用最小(给出总费用)。
(2)请就
(1)的模型分析:
哪个钢厂钢管的销价的变化对购运计划和总费用影响最大,哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大,并给出相应的数字结果。
(3)如果要铺设的管道不是一条线,而是一个树形图,铁路、公路和管道构成网络,请就这种更一般的情形给出一种解决办法,并对图二按
(1)的要求给出模型和结果。
A1
3
2
5
80
10
10
31
20
12
42
70
10
88
10
70
62
70
30
20
20
30
450
104
301
750
606
194
205
201
680
480
300
220
210
420
500
600
3060
195
202
720
690
520
170
690
462
160
320
160
110
290
1150
1100
1200
A2
A3
A4
A5
A6A11
A711A11
A8A11
A911A11
A10
A11
A12
A13
A14
A15
S1
S2
S3
S4
S5
S6
S7
图一
A1
3
2
5
80
10
10
31
20
12
42
70
10
88
10
70
62
70
30
20
20
30
450
104
301
750
606
194
205
201
680
480
300
220
210
420
500
600
3060
195
202
720
690
520
170
690
462
160
320
160
110
290
1150
1100
1200
A19
130
190
260
100
A2
A3
A4
A5
A6
A7
A8A11
A9
A10
A11
A12
A13
A14
A15
S1
S2
S3
S4
S5
S6
S7
A16
A17
A18
A20
(A21)
图二
二、模型假设
1、假设沿管道或者原来有公路,或者建有施工公路;
2、运费只按铁路、公路里程收取,即不考虑火车、汽车由于停靠站等其他一切外因带来的费用;
3、钢管在铺设过程中以1km为单位进行铺设;
4、钢管可由铁路、公路运往铺设路线任一地点;
5、所有钢管在指定期限内都能按时生产并运送指定地点;
6、钢管铺设过程中由站点向左右两边进行铺设。
三、符号说明
:
第个厂;
:
第个站点;
:
向运送的钢管量单位(km);
:
在指定期限内的最大生产量单位(km);
:
向右铺设的钢管量单位(km);
:
向左铺设的钢管量单位(km);
:
到间的距离单位(km);
:
管道全线总长单位(km);
:
钢管出厂销价单位(万元/单位);
:
向运送一单位钢管所需的铁路费单位(万元/单位);
:
向运送一单位钢管所需的公路费单位(万元/单位);
:
购买钢管所花的总费用;
:
由厂到站点所需运输总费;
:
由站点到铺设地点所需运输总费;
:
订购和运输钢管所需总费用单位(万元)。
四、问题分析
问题一是在一定约束条件下的非线性优化问题,由题意知,拟建立以总费用为目标函数来寻求最优解。
总费用由钢管的购买费、厂到站点的运输费以及站点到铺设地点的运输费三部分组成。
一、钢管的购买费可由在每个厂的购买量与每个厂的出厂销价的线性运算得到。
在每个厂购买的钢管量必须大于500km,否则则不在该厂购买。
可以构造一个的矩阵,那么当为0时,表示不在第个钢厂购买,否则则在第个钢厂购买大于500km的钢量。
二、要求得每个钢厂到站点的运输费需先知道每个厂到各个站点的钢管输送量,以及所选择的路线即铁路总长和公路总长,所以需要首先计算出各个钢厂到每个站点的最佳运输路径,使得平均单位公里的运输费用最小。
但是由于铁路每公里的运输费用不是线性变化,而是变化不均匀的分段函数。
在这里,我们利用深度优先遍历,找到某个厂到达各个站点的所有路径,然后根据每条路径的铁路和公路里程数计算出平均每公里运输费用最小的一条。
以此类推,计算出所有钢厂到所有站点的最佳路径。
三、在站点到铺设地点的运输费问题上,如果我们认为车边向前走边进行铺设,即边走边将钢管放下,那么就需要通过积分来计算。
但是,尽管用积分算下来结果会很精确,但在实际中不可能这样实施。
另外,这也与题目中不足整公里的按整公里计算相矛盾。
所以,我们假设以1km为单位进行铺设,即铺设中车每向前开1km便将1km的钢管放下。
由于铺设管道是线型的,除了两个端点外,每个站点需要往两边进行铺设管道。
所以,假设第个站点往左、右边铺设管道为和公里,则由站点到铺设地点的运输费就可以通过等差数列求和得到。
问题二即为对问题一中模型的灵敏度分析,在讨论各厂的钢管销价和生产上限对购运计划和总费用的影响时,只让其中一个量变化,其他一切条件皆不变,即逐个变量单独分析。
问题三即为问题一中模型的推广,在问题一的基础上将站点向左右两边铺设变为向三个方向铺设,按问题一处理即可。
五、模型建立(问题一)
总费用由钢管的购买费、厂到站点的运输费以及站点到铺设地点的运输费三部分组成,则
在第个厂的购买费应为15个站点在第个厂的购买总量与该厂销价的乘积总和,即,则总购买费
第个厂向第个站点的运输费为运送量与运送1单位所需铁路费和公路费的和的乘积,第个厂向各个站点运送钢管的总运费即为,则各厂到站点的运输费
要算出钢管由站点运送到铺设地点的费用需知道钢管按何种方式进行铺设的。
在问题分析里一讨论边走边铺与实际不符,且有违题目条件,所以我们假设钢管在铺设过程中以1km为单位进行铺设,且由站点向两边进行铺设,则可由等差数列求和公式得到,即
由于一个钢厂如果承担制造这种钢管,至少需要生产500个单位,且各厂在指定期限内有生产上线,则在第个厂的购买总量需满足
或
钢管由站点向左右两边进行铺设,则第个站点向右铺设部分与第个站点向左铺设部分之和应为两站点之间的管道长度,且第一个站点向左铺设部分与最后一站点向右铺设部分都为0,即
第站点向左铺设部分与向右铺设部分之和应为七个厂向第站点输送钢管总量,即
综合考虑钢管的购买费、厂到站点的运输费以及站点到铺设地点的运输费,钢管的订购和运输优化模型建立如下:
目标函数
min++()
或
s.t
六、模型求解
由于铁路、公路相互交错,无法直接选出钢厂到站点的费用最小路线,所以此处我们采用深度优先遍历方法。
首先建立一个39维数组,将图一中39个交点两两之间有铁路、公路连接的用具体路线长写入数组,且铁路用负数表示,公路用正数表示,而没有路线连接的用无穷大代替,最后换算成到各站点的铁路、公路总费。
全过程通过matlab编程完成(程序见附录),。
表一到的最小费用(单位:
万元/单位)
170.7
215.7
230.7
260.7
255.7
265.7
275.7
160.3
205.3
220.3
250.3
245.3
255.3
265.3
140.2
190.2
200.2
235.2
225.2
235.2
245.2
98.6
171.6
181.6
216.6
206.6
216.6
226.6
38
111
121
156
146
156
166
20.5
64.6
105.5
139.6
130.5
140.5
150.5
3.1
86
96
131
121
131
141
21.2
71.2
86.2
116.2
111.2
121.2
131.2
64.2
114.2
48.2
84.2
79.2
84.2
99.2
92
142
82
62
57
62
76
96
146
86
51
33
51
66
106
156
96
61
51
45
56
121.2
171.2
111.2
76.2
71.2
26.2
38.2
128
178
118
83
73
11
26
142
192
132
97
87
28
2
因为matlab无法直接对约束条件或进行处理,所以我们先将此条件改为,则原模型变为
min++()
s.t
通过matlab编程(程序见附录)计算结果见表二
表二各厂的生产量及总费用(生产量可小于500)(单位:
单位、万元)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学 建模 钢管 订购 运输