钢管订购和运输论文

承诺书我们仔细阅读了全国大学生数学建模的竞赛规则（）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与本队以外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们愿意承担由此引起的一切后果。我们的参赛（报名）队号为：32参赛组别（研究生或本科）：本科参赛队员：兰潇根、柳达强、汪锡平钢管订购和运输摘要：本文拟建立一个最合理的钢管运输与铺设方案模型。利用离散数学和数据结构中图论相关知识，应用最短路径的floyd算法和灵敏度分析法建立一个以总费用为目标函数的非线性规划模型，对于钢管订购和运输的总费用，分为三部分：购买钢管费用，由钢厂运送到站点的费用以及由站点开始铺设的费用，对于由钢厂运送到站点的费用，用Floyd算法，求出铁路网和公路网的最短路径，然后转化为最少运输费用，之后利用Lingo软件编程，求解分析，解决问题。关键词：Floyd算法，非线性规划，Lingo1一问题重述要铺设一条1521AAA的输送天然气的主管道,如题图一所示。经筛选后可以生产这种主管道钢管的钢厂有721,,SSS。图中粗线表示铁路，单细线表示公路，双细线表示要铺设的管道(假设沿管道或者原来有公路，或者建有施工公路)，圆圈表示火车站，每段铁路、公路和管道旁的阿拉伯数字表示里程(单位km)。为方便计，1km主管道钢管称为1单位钢管。一个钢厂如果承担制造这种钢管，至少需要生产500个单位。钢厂iS在指定期限内能生产该钢管的最大数量为is个单位，钢管出厂销价1单位钢管为ip万元，如下表：i1234567is80080010002000200020003000ip1601551551601551501601单位钢管的铁路运价如下表：里程(km)≤300301～350351～400401～450451～500运价(万元)2023262932里程(km)501～600601～700701～800801～900901～1000运价(万元)37445055601000km以上每增加1至100km运价增加5万元。公路运输费用为1单位钢管每公里0.1万元（不足整公里部分按整公里计算）。钢管可由铁路、公路运往铺设地点（不只是运到点1521,,,AAA，而是管道全线）。（1）请制定一个主管道钢管的订购和运输计划，使总费用最小（给出总费用)。（2）请就（1）的模型分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果。（3）如果要铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，请就这种更一般的情形给出一种解决办法，并对题图二按（1）的要求给出模型和结果。2二问题分析1.问题一所有的钢管必须通过铁路运送到铺设线路上的站点1521AAA，之后再通过公路运输向左或右铺设。因此，总的费用由三部分组成：一部分为购买所有主管道钢管的总费用，一部分为由钢管厂运送到各个站点时的铁路运费和公路运费的总和，最后一部分为由站点向左右两边铺设时的运输费用。对于从钢管厂到各个站点的最小运费，由于在铁路和公路上的运费计算方法不同，所以，可以先用Floyd算法，求出钢管厂到铁路上任意结点的最小距离和路线，得到相应的单位钢管运费，同理再求出各个站点到公路上任意结点的最小距离和路线，得到相应的单位钢管运费，再将两运费求和求出最小值，于是就得到从某钢厂到某铺设地点运输单位钢管的最少运输费用。2.问题二题目中“哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大”可以理解为，当该模型达到最优解时，钢管销价或者产量上限变化一个单位时，对购运计划和总费用的影响的大小问题。可以利用Lingo编程运行得到结果。3.问题三要铺设的管道是一个树形图，是题图一的一种延拓，通过观察可知，只有9、11、17站点的铺设方向有三个，其它站点的铺设方向只有左右，因此，可以沿用问题一里的思路，在问题一的基础上再增加一个变量middle（j），用于表示向第三方向铺设的钢管数量。三模型的假设与符号说明1、模型的假设⑴.沿管道或者原来有公路或者建有施工公路。⑵钢管全部由这7个钢厂生产，一个钢厂如果承担制造这种钢管，至少需要生产500个单位。⑶公路运输费用为1单位钢管每公里0.1万元（不足整公里部分按整公里计算）。⑷由于公路运输费中，不足整公里部分按整公里计算，因此，从站点向左右两边运送钢管时，不应该是边运送边卸下钢管，这样也不符合实际，应当是走一个单位的公路，卸下一个单位的钢管。32、符号说明符号说明is钢厂iS在指定期限内能生产该钢管的最大数量ip钢厂iS的钢管出厂单位销价（单位：万元）cost(i,j)单位钢管从钢厂iS运到jA的最小费用（单位：万元）l(j)从jA到1jA之间的距离（单位：千米）n(i,j)从钢厂iS运到jA的钢管数量left(j)运到jA站点向左铺设的钢管数量right(j)运到jA地的钢管向右铺设的钢管数量middle（j）运到jA站点的钢管向第三方向铺设的钢管数量c(i)=0钢厂iS不提供钢管c(i)=1钢厂iS提供钢管四模型的建立与求解（一）、问题一的模型：采用Floyd算法，用matlab编程可以求出单位钢管从iS运输到jA的最小运输费用，数据如下表：4A1A2A3A4A5A6A7A8A9A10A11A12A13A14A15S1170.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.276665638.2262目标函数为总费用W，包括三个部分，购买所有主管道钢管的费用1W，将钢管从钢厂运到各个站点的费用2W，将钢管从站点运到铺设地点的费用3WW=1W+2W+3W其中),(*711511jinpWiji),(*),(cos711512jinjitWij2)()(12)()(1*1.01513jleftjleftjrightjrightWj则目标函数：minW=),(*)(71151jinipij+),(*),(cos71151jinjitij+2)()(12)()(1*1.0151jleftjleftjrightjrightj约束条件：1.钢厂的钢管产量：)(*),()(*500151icsjinicij52.运到各个站点的钢管刚好用完：)()(),(71jrightjleftjini（j=1,2…15）3.jA与1jA之间的钢管：)()()1(jljrightjleft,（j=1,2,…,14）4.钢管数量的非负性：n(i,j)≥0,left(j)≥0,right(j)≥0(i=1,2,…,7,j=1,2,…,15)5°钢管数量的整数性：n(i,j)∈N运用数学软件Lingo编程求解问题一的结果最优最小费用1278632W(万元)（二）、问题二的模型：用Lingo对问题一求解后，即可根据Lingo的结果对问题二进行解答。各钢厂销价的变化：p(1)p(2)p(3)p(4)p(5)p(6)p(7)对偶价格-800-800-10000-1320-1250.990对偶价格表示，在最优解的情况下，各钢厂钢管销价减少一个单位时，对总费用的影响。根据表中的数据，S(5)钢厂钢管的销价对购运计划和总费用影响最大。产量上限：s(1)s(2)s(3)s(4)s(5)s(6)s(7)对偶价格10335253.330016对偶价格表示，在最优解的情况下，各钢厂钢管生产上限每增加一个单位时，对总费用的影响。根据表中的数据，得S(1)钢厂钢管的产量上限的变化对购运计划和总费用的影响最大。（三）、问题三的模型：6采用Floyd算法，用matlab编程求出单位钢管从iS运输到jA的最小运输费用，数据如下表：A1A2A3A4A5A6A7A8A9A10S1170.7160.3140.298.63820.53.121.264.292S2215.7205.3190.2171.611195.58671.2114.2142S3230.7220.3200.2181.6121105.59686.248.282S4260.7250.3235.2216.6156140.5131116.284.262S5255.7245.3225.2206.6146130.5121111.279.257S6260.7250.3235.2216.6156140.5128.1116.284.261S7275.7265.3245.2226.6166150.5141131.299.276（续表）A11A12A13A14A15A16A17A18A19A20A21S196106121.21281426095100105115125S2146156171.2178192110145150155165175S38696111.211813244859095105115S4516176.28397805055607080S5335171.27387753245506575S6473716.2112880463336100S7645638.2262956350553226由于树形图的出现，发现在站点9、11、17处出现了3条支路的情况。则模型一中模型的变量left(j)，right(j)不再适用，此时可考虑增加一个支路变量middle（j），相应的增加约束条件，在目标函数中增加相应的从站点运到铺设地点的费用。目标函数：),(*),(cos),(*)(min7121171211jinjitjinipWijij)20,19,17,17,11,9(21)()(2)(1)m(2)()(12)()(1*1.0212141nmnrightnrightmmiddlemiddlejleftjleftjrightjrightjj）（）（约束条件：1.钢厂的钢管产量：)(*i),()(*500211icsjinicj）（)7,...,1(i72.运到各个站点的钢管刚好用完：)()(),(71jrightjleftjini(j=1,…,21且j≠9,11,17))()()(),(71jmiddlejrightjleftjini（j=9,11,17）3.jA与1jA之间的钢管：)()()1(jljrightjleft（j=1,2,…,14）middle(9)+left(16)=42middle(11)+middle(17)=10left(17)+left(18)=130right(17)+left(19)=190right(19)+left(20)=260right(20)+left(21)=1004