管路订购和运输数学建模

./管路订购和运输[摘要]本文是要研究并确定分别对7个钢管厂订购钢管量,并运输到15个节点以及钢管铺设线路所花费的费用最省。我们建立模型,并通过求解模型能更好的解决次类问题。更能使人一目了然的对每个钢厂订购钢管情况,管路运输路线,以及铺设路线的费用有更深刻的理解,并做出最优化的选择方案。问题一：制定一个主管道钢管的订购和运输计划,使其总费用最省。根据题目要求,以及lingo求解得出对钢厂订购量分别为：S1=800,S2=800,S3=1000,S500=1246,S6=1325,运输计划为：S1→A4=207S2→A2=179S2→A8=111S3→A9=664S5→A10=231S6→A13=333S1→A5=127S2→A3=92S3→A3=129S5→A3=287S5→A11=415S6→A14=621S1→A6=200S2→A4=152S3→A4=98S5→A4=10S6→A10=120S6→A15=165S1→A7=266S2→A5=77S3→A5=109S5→A5=303S6→A12=86最终求得最省总费用为1278632万元。问题二：在问题一的基础下进行分析,并对问题一的模型进一步整理,分别把每个钢管厂的生产上限减少500万单位钢管,得到答案记录并制作图表对比分析,S1钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大。假设每个钢厂钢管的销价分别增加200万元对其运输计划的影响。其它数据不变,并利用LINGO软件求解,对得到的数据会制成表格见表〔3。分析对比表格可以得出S1钢厂钢管的销价的变化对购运计划和总费用最大。关键词：最短路径分部计算法建立模型1问题重述1.1问题背景要铺设一条的输送天然气的主管道,如图一所示<见附录>。可以生产这种主管道钢管的钢厂有,一个钢厂如果承担制造这种钢管,至少需要生产500个单位。钢厂在指定期限能生产该钢管的最大数量为个单位,钢管出厂销价1单位钢管为万元,如下表：1234567800800100020002000200030001601551551601551501601单位钢管的铁路运价如下表：l里程<km>≤300301～350351～400401～450451～500运价<万元>2023262932里程<km>501～600601～700701～800801～900901～1000运价<万元>37445055601000km以上每增加1至100km运价增加5万元。公路运输费用为1单位钢管每公里0.1万元〔不足整公里部分按整公里计算,钢管可由铁路、公路运往铺设地点〔不只是运到点,而是管道全线。1.2问题提出由上述问题的背景,需要我们研究下列问题〔1需要制定一个主管道钢管的订购和运输计划,使得总费用最小。〔2就〔1的模型分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大,哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大。1.3符号设置:钢管厂Si的出场钢管单位价格〔万元:公路上一单位钢管的每公里费用〔d=0.1万元

:为1单位钢管从Si到Aj最小费用

:从Aj到Aj+1的距离〔千米:

钢管运到Aj向左铺设

:钢管运到Aj向右铺设2问题分析因为铁路和公路运输不能直接应用最短路径算法来求解铁路和公路交通路线中任意两点间的最小费用问题。所以我们采用了分部计算法,即运输费为铁路运输费加公路运输费,结合问题图表给出的路程距离可以求出不同路线所需的费用,所有运输线路的单位运输钢管的费用〔单位钢管铁路公路的运费和,在不考虑路程长短只考虑运输费用最省的前提下只保留运输费用最少的路线的费用并记录,记录的数据再乘以每个钢管厂所订购的钢管量,即为运输钢管费用的总和〔P。在对钢管在铺设时在交叉点分别向左向右铺设综上在利用lingo求解。在问题一中,我们要先求出每个钢管厂Si分别所要订购多少钢管和所需订购的费用〔M,以及经铁路与公路运输钢管所需最少费用的总和〔P,还有铺设钢管费用的总和〔T,其中钢管在铺设时在交叉点分别向左向右铺设在铺设费用计算时其放下过的路段费用不重复计算。在问题二中,我们在问题一的基础上对问题及模型进行分析,对模型数据进行调整。分析每个钢管厂的生产上限减少500万单位的钢管,其它数据不变,并且利用问题一的软件进行求解,得到数据并绘制表格〔3得到哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大。在分析每个钢厂钢管的销价增加200万元,其它数据不变,并利用软件求解,对得到的数据会制成表格,分析对比表格可以得出哪个钢厂钢管的销价的变化对购运计划和总费用最大。3问题假设假设一：在运输过程中钢管无损耗无丢失。假设二：在运输过程中只考虑运输和铺设费用不含其它的费用。4问题求解4.1问题一根据题目中已知的数值,结合问题图表给出的数据计算：第一步我们先求出在钢管厂订购钢管最省费为〔订购钢管费用第二步因为路程距离可以求出不同路线所需的费用,即求出钢管厂〔到15个节点〔所有运输线路的单位运输钢管的费用〔单位钢管铁路公路的运费和在不考虑路程长短,只考虑运输费用最省的前提下,只保留运输费用最省路线的费用并记录保存,记录的数据绘制成表格如下表：表1单位钢管从钢管厂运到各个节点的最省费用A1A2A3A4A5A6A7A8A9A10A11A12A13A14A15S1170.7160.3140.298.63820.53.121.264.29296106121.2128142S2215.7205.3190.2171.611195.58671.2114.2142146156171.2178192S3230.7220.3200.2181.6121105.59686.248.2828696111.2118132S4260.7250.3235.2216.6156140.5131116.284.262516176.28397S5255.7245.3225.2206.6146130.5121111.279.257335171.27387S6265.7255.3235.2216.6156140.5131121.284.262514526.21128S7275.7265.3245.2226.6166150.5141131.299.277665638.2262即最省运输钢管费用为第三步是求出钢管铺设的费最省费用总和〔铺设钢管费用的总和最后求出最省总费用<总费用>4.2目标建立根据题目已知条件,建立目标函数：在钢管厂要么订购钢管要么不订购,且一个钢厂如果承担制造这种钢管,至少需要生产500个单位2.钢管在铁路运输中1000km以上每增加1至100km运价增加5万元。3.公路运输费用为1单位钢管每公里0.1万元〔不足整公里部分按整公里计算。4.钢管可由铁路、公路运往铺设地点〔不只是运到点,而是管道全线。4.3模型建立综合上述步骤分析,建立模型如下：4.4模型求解1.根据数学模型,把数值代入模型,利用软件进行求解如下：得到最省费用为1278632〔万元2.运输计划：到的钢管数,根据数学模型,把数值代入模型,利用软件进行求解如下绘制如表3.经计算得出主管道钢管在钢管厂订购主管道的订购计划：=800km=800km=1000km=1246km=1325km4.利用软件进行求解得出每个点向左右各铺的长度经绘制如表：4.5问题二就〔1的模型与软件分析,假设每个钢管厂的生产上限减少500万单位的钢管,其它数据不变,并且利用问题一的软件进行求解,得到数据并绘制表格如下表：从每个点向右铺的长度A1A2A3A4A5A6A7A8A9A10A11A12A13A14A15总和Z0104226468606184.5189.5125505321270751992861653724Y07528209.515.576175159301451113433501447总和0179508468615.5200265.5300664351415863336211655171由此表格分析得到钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大。假设每个钢厂钢管的销价增加200万元对运输计划和运费的影响,其它数据不变,并利用软件求解,对得到的数据会制成表格如下：分析对比表格可以得钢厂钢管的销价的变化对购运计划和总费用最大。4.4结果验证根据软件求解验证S1到A4.5.6.7点运输为207.127.200.266km总和=800kmS2到A2.3.4.5.8点运输为179.92.152.77.300km总和=800kmS3到A3.4.5.9点运输为129.98.109.664km总和=1000kmS5到A3.4.5.10.11点运输为287.10.330.231.415km总和=1246kmS6到A10.12.13.14.15点运输为120.86.333.621.165km总和=1325km附录

问题一

Lingo程序

sets:

aa/1..7/:p,s,t;

bb/1..15/:y,z;

cc<aa,bb>:c,x;

dd/1..14/:d;

endsets

data:

p=160

155

155

160

155

150

160;

s=800

800

1000

2000

2000

2000

3000;

d=104,301,750,606,194,205,201,680,480,300,220,210,420,500;

c=170.7000160.3000140.200098.600038.000020.50003.100021.200064.200092.000096.0000106.0000121.2000128.0000142.0000215.7000205.3000190.2000171.6000111.000095.500086.000071.2000114.2000142.0000146.0000156.0000171.2000178.0000192.0000230.7000220.3000200.2000181.6000121.0000105.500096.000086.200048.200082.000086.000096.0000111.2000118.0000132.0000260.7000250.3000235.2000216.6000156.0000140.5000131.0000116.200084.200062.000051.000061.000076.200083.000097.0000255.7000245.3000225.2000206.6000146.0000130.5000121.0000111.200079.200057.000033.000051.000071.200073.000087.0000265.7000255.3000235.2000216.6000156.0000140.5000131.0000121.200084.200062.000051.000045.000026.200011.000028.0000275.7000265.3000245.2000226.6000166.0000150.5000141.0000131.200099.200077.000066.000056.000038.200026.00002.0000;

enddata

min=sum<cc<i,j>:p<i>*x<i,j>>+sum<cc<i,j>:c<i,j>*x<i,j>>+0.05*sum<bb<j>:z<j>^2+z<j>+y<j>^2+y<j>>;

for<aa<i>:sum<bb<j>:x<i,j>>>=500*t<i>>;

for<aa<i>:sum<bb<j>:x<i,j>><=s<i>*t<i>>;

for<aa<i>:bin<t<i>>>