物流网络选址与路径优化问题的模型与启发式解法
物流网络选址与路径优化问题的模型与启发式解法120交通运输工程学报2006年(i∈T)(10表1小规模问题的计算结果Tab. 1 Computation result of smalF-scaled problem(xh+xbk)-z≤1(i∈T,∈Ck∈K)(11)最优解∈TUC(j∈C问题规模目标值运行时间/s(12)NG= 3, NT=3, Nc= 83885k≤s1(≤c1121=∑体k∈(13NG=3,Nr=3,Nc=10421.221U-U+Nx≤N-1(i,∈S,k∈K)(14)NG=3,Nr=3,Nc=126850580∈TNG=3,NT=3,Nc=148741162000注:Nc为供应商数量,Nr为配送中心数量,Nc为客户数量ba≤B(g∈G)(16用启发式算法求得初期解xi∈/Q1(i,j∈T∪C,k∈K)(17)l0;∈/0,1∈T(18通过交换配送中心间的路径进行第1次解的改善∈/O,1∈T,j∈C)(19)yk∈/0,1(i∈C,k∈K)(20)通过交换同一路径中客户的位置(2-OPT法)进行第2次解的改善式中:G为供应商的集合;T为配送中心的集合;C为客户的集合;K为车辆的集合;Qk为车辆k的最通过交换不同路径中的客户进行第3次解的改善大载质量:S为C的部分集合;C为从点到点j的距离;B为供应商g的最大供应量;V为配送中心SA模块的最大货物通过量;D,为客户j的需求量;H2为配送中心i的固定费用;L为从供应商g到配送中是否满足终止条件心i的单位运输费用;Fk为车辆k的固定费用;a为解的输出与通过量有关的系数;xk为对于车辆k,如果点i以后的访问点是点j,即为1,否则为0;y为如果点j图1基于SA的混合启发式算法的货物由车辆k配送,即为1,否则为0;v;为如果FiMixed heuristic algorit hm based on sa使用配送中心i,即为1,否则为0;z为如果客户j传统启发式算法与智能启发式算法相结合的混合算由配送中心提供服务,即为1,否则为0;pa为从供法,以期在短时间内求得全局最优解应商g到配送中心的供应量。其中x、3计算分析和pg为决策变量。2问题的求解为了对SA的参数进行设定,进行了预备实验并确定参数如下:初始温度To为200冷却率α为为了验证上述数学模型的正确性,用数理规划07,与温度相关的循环次数调整参数β为.1,最商用软件 LINGO8.0对小规模问题进行了数值计大循环次数将按照问题规模的大小做适当的设定,算,结果见表1。可以看出,随着问题的规模扩大,数据采用人工生成数据,在200km×200km的区计算时间急剧增加;当客户达到14个时,计算时间域内随机生成指定个数的供应商、配送中心和客户,长达45h,显然无法满足解决现实问题的需要。为并生成距离矩阵和客户需求量;采用C语言编程,了满足解决现实问题的需要,有必要开发岀一种能计算结果见表2。从表2中的结果比较可以看出,够在合理的时间内求得准最优解的近似算法。表2最优解与近似解比较传统启发式算法能够在短时间内求得局部最优Tab 2 Comparison of optimal solution and approxi mate solution解,但往往容易陷于局部最优,而无法求得全局最优最优解近似解问题规模解。智能启发式算法能够求得全局最优解,但计算时目标值运行时间/s目标值运行时间/s%间相对较长。如果能够将两者结合起来,既可以防止N=3N=3,Nc=838533965求解过程陷于局部搜索无法跳出,保证全局解的搜3N=3NG=10+2122112100索,又可以缩短搜索时间达到在短时间内求得全局N-3-3Nc12|60s0620110最优解的目的。基于上述考虑本文提出了图1的将ublishrigfoustAingnsestrvcu,trttp/www.urrhi.rctNG=3,Nr=3,Nc=1487411620008895第3期陈松岩,等:物流网络选址与路径优化问题的模型与启发式解法121本算法求得的近似解与 LINGO80计算的最优解心,再从配送中心到客户这一典型的物流过程,涉及之间的误差很小,但计算时间却大大缩短了。依据运输与配送2个阶段和供应商、配送中心和客户3结果虽然无法判定所提岀的混合启发式算法对于大个层次,提出了多供应商、多物流中心情况下的物流规模问题的有效性,但可以看岀,对于求解小规模问路径与配送路径优化问题,给岀了问题的数学模型,题是有效和良好的。对于大规模问题,将利用实例利用传统启发式算法与模拟退火法开发了混合近似进行计算验证。解法,通过人工生成数据和实例计算验证,可以看出4应用实例所提出的数学模型可以准确地描述此类问题,具有良好的适应性,所提出的混合近似解法能够在短时在应用实例中,将港口作为供应商来考虑,以进间内求解问题并得到接近于最优解的近似解,具有口货物从港口经配送中心配送到客户过程中发生的较高的实用价值。但本模型没有考虑库存问题与供费用最小化为目的,以港口的数量和位置、配送中心应商的成本问题,无法达到物流网络中各个环节的数量和位置作为对象进行优化。实例的区域选择山整体优化,有待于今后进一步研究。东省,候选港口为天津港、烟台港、威海港、青岛港参考文献日照港和连云港等6处,候选配送中心设置于山东省除港口城市以外的13个地级市,设定客户分别位References于90个县(包括县级市)。为了分析候选港口和配1 I anen g Flpo C. Spatial de composit ion for a multI送中心的数量及位置与目标值之间的关系,在计算of Production Economics, 2000, 64(1/2/3): 177186过程中,候选港口和配送中心的数量分别从1开始(21 Melkote s, Das kin m s. a n in te grated model of facil ity loca tion增加(港口的位置为随机选择),计算结果见图2、3and transportat ion netw ork design[ J. Transport ation Research part A,2001,35(6):5155386[3] Goets chalckx M, Vidal C J Dogan K Modeling and design ofglobal logistics s yst ems: a review of integrated strat e gic an dt act ical models and design algorithm s[ J European Journal of005◇◇0◇◇◇◇◇◇◇Operat ional Research, 2002, 143 (1):118135791113[4] H wang H S Design of suppl y chain logist ics sy stem con sider-港口数量配送中心数量ing service level[ J. Computers and Industrial Engineerin g,图2港口数量与图3配送中心数量与2002,43(1/2):283297目标函数值关系目标函数值关系[5] WuT H, Low C, Bai J W. Heurist ic s ol ution s to mu ltt depotFig2 Relation of ob ject value Fig 3 Relation of ob ject valuelocatioN rou tin g pr ob lems[ J]. Computers and Operations Re-and ports numberand changing depots num besearch,2002,29(10):13931415可以看出,随着候选港口数量的增加.目标值呈[6 Syam ss. A model and met hodologies for the locat ion p rob lem下降趋势,说明可供选择的港口越多,求得最优解的with logist ical components[ J]. Computers and Operations Re机会越大,但本例中,当候选港口的数量增加到5个se arch,200)2,29(9):l173-1193时,目标值达到最小(实际被选中的港口为3个);候[7 Amiri A. Designing a distribution network in a suppl y chainsystem[ J]. European Jou rnal of O perat ional Research, 2004选配送中心数量的变化也呈相同的趋势,当候选配171(2):567576送中心数量达到9个时,目标值达到最小(实际被选8 G ena m, Syarif a. Hybrid genetic algorit hm for mult+ time pe-中的配送中心为8个);本实例的计算时间都在8sriod production/ distribution planning[ J1. Computers and Ir以内,虽然无法判断所求解为最优解,但从计算结果dustrial Engineering, 2005, 48(4): 799809来看,基本接近最优,因此可以认为本算法对于求解9工丰元,潘福全张丽霞等基于交通限制的路网最优路径算大规模问题也是有效和良好的法J.交通运输工程学报,2005,5(1):9295Wang Feng yuan, Pan Fuquan, Zhang Li xia, et al. Opti mal5结语path algorithm of road netw ork with traffic rest riction[ JIJourn al of Traffic and Trans port ation Engin eering, 2005, 5(1)本文将研究范围界定在商品从供应商到配送中9295.(in Ch ineseo1994-2012ChinaacAdemicJournalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
- 2020-12-11下载
- 积分：1
FLUENT经典算例合集7个
案例分析的电子书《数值计算与工程仿真》专刊一 FLUENT HELP算例精选中文版(→)目录前言..2目录FLUENT经典算例翻译之一算例1介绍如何使用 Fluent算例4非定常可压缩流动模型.52算例5辐射与自然对流模拟·:······99FLUENT经典算例翻译之二算例13使用非预混燃烧模型151算例15蒸发性液体喷雾建模.∴214算例18使用混合物多相模型和欧拉多相模型..∴252算例21使用欧拉多相粒子传热模型垂··D垂垂垂。看垂垂垂看D看垂·D4。垂278ww.myCFDcn《数值计算与工程仿真》专刊一 FLUENT HELP算例精选中文版(→)算例1介绍如何使用 Fluent引言此向导通过图例说明了一个发生在混合弯管处的两维湍流流动和传热问题的求解方法和过程。这种混合弯管的结构常见于发电厂和化工厂的管道系统中。正确计算出弯管内流体交汇处附近的流场和温度场分布,对于设计合适的入口管道位置具有重要意义通过此向导,用户可学会以下内容●在 Fluent中输入网格文件使用混合单位制去定义儿何体和流体的属性设置强制对流的湍流流动的流体物性和边界条件迭代计算并使用残差监视器监测计算过程及其收敛性●使用隔离求解器进行求解使用等势图检察流场和温度场●运用二阶离散化方法重新计算以获得更佳的温度分布对网格进行温度梯度自适应,进一步求解更佳的温度场分布前提条件在学习此向导之前,假设用户还没有使用 Fluent的给验,不过,已经学习过用户指南第一章中的简单算例,并且熟悉 Fluent的界面及其指南中的规约可题描述问题如图1-1所小。一股温度为26℃的冷流体流入大管道,在弯管处与另股温度为40℃热流体混合。管道的长度单位为英寸,而流体的属性和边界条件则使用国际单位。入口管道的雷诺数为2.03×105,因此,选择湍流流动模型。4ww.myCFDcn《数值计算与工程仿真》专刊一 FLUENT HELP算例精选中文版(→)P =100 kg/viscosity8x14FConductivity: k=0.537 Y4 miK2SpC Ic Heat 9=4216 Jkg-k6酽.2mT1121m图图1-1问题说明准备工作1.从 Fluent的文件光盘中拷贝文件 elbow/ elbow, msh到电脑的 Fluent作日录中对于Unⅸx系统,当把文件光盘放入电脑光驱后,可以在以下目录找到这个文件:/ edrom/uent61/help/ tuttles述 cdrom为电脑的光驱目录对于 windows系统,当把文件光盘放入电脑光驱后,可以在以下冂录找到个文件cdrom: fluent 6.1 help tutfiles上述 cdrom为电脑的光驱目录2.启动 Fluent,选择2D求解器。ww.myCFDcn《数值计算与工程仿真》专刊一 FLUENT HELP算例精选中文版(→)第1步:与网格相关的操作读取网格文件 elbow, mshFilc→Read+}CascSelect fileFiterh,于d,SH,GRdirectoriesFileshome user tutorial/elbow. msh/home userautorialaCase Filehome /user/tutorialoKFiterCancela)在 Files项中点击选中 elbow. msh,然后点击OK完成操作。注意当 Fluent读取网格文件的同时,信息会不断显示在反馈窗口内,报告网格转化的过程。当读取网格文件完毕, Fluent的反馈窗凵会显示共读取了918个三角形的流体单元,以及许多带着不同分区标识符的边界面。2.网格检查。Grid→} Check6ww.myCFDcn《数值计算与工程仿真》专刊一 FLUENT HELP算例精选中文版(→)F1uent的信息反馈窗口会显示如下信息:Grid checkDomain extentscoordinate:min(m)=0.0000009+00,max(m)=6.400001e+01y- coordinate:min(m)=-4.538534e+00,max(m)=6.400000e+01Volume slatisticsminimum volume (m3): 2.782193c-01maximum volume (m3):3.926232e+00total volume (m3):1.682930e+03Face area statisticsminimum face area (m2):8.015718e-01maximum face area(m2):4. 118252e+00Checking number of nodes per cellChecking number of faces per cellChecking thread pointers.Checking number of cells per faceChecking face cellsChecking bridge facesChecking right-handed cellsChecking face handednessChecking element type consistencyChecking boundary typesChecking face pairs.Checking periodic boundaries.Checking node countChecking nosolve cell countChecking nosolve face countChecking face childrenChecking cell childrenChecking storageww.myCFDcn《数值计算与工程仿真》专刊一 FLUENT HELP算例精选中文版(→)Done注意网格检查结東后,信息反馈窗∏会以默认的SI单位制给出网格在ⅹ轴和Y轴上的最大和最小值,并将报告出网格的賦它特性。网格检查还会报告出有关网格的任何错误。需要特别注意的是,确保最小体积不能是负值,否则 Fluent无法进行计算。在SI单位制中,默认单位是m,若想改变单位制,使用 in ches:可以打开 Scale grid对话框。3.平滑(或者交换)网格Grid→} Smooth/swap…Smooth/Swap GridsmoothSwap InfoMethodNumber SwappedskewnessMinimum skewnesscumber visitedNumher of lerationsSmoothCloseFluent读取网格文件后,平滑三角形或四边形网格是一个良好的习惯,那样能确保使用质量铰好的网格进行计算。a)点击按钮 Smooth,再点击按钮Swap,重复上述操作,直到 Fluent报告没有需要交换的面为止。若 Fluent再无法通过交换改善网格质量,则没有平面可被交换了。b)点击 Close关闭对话框。ww.myCFDcn《数值计算与工程仿真》专刊一 FLUENT HELP算例精选中文版(→)4.更改网格的长度单位Grid→ scalea)在 Units conversion(单位转换)项的 Grid was created ln(网格长度单位)的右侧下拉列表中选择In(代表选择了英寸b)点击 Scale按钮,更改长度单位。在 Domain extents栏中采用了默认的SI单位制,长度单位为mc)点击按钮 Change length Units,设定 inches(英寸)为此次计算采用的长度单位确保Xmax(in)和Ymax(in)中数值为64英尺。(如图1.1)Scale gridScale factorsUnit a conversion00254rid Was Created In inY0+0254Change Length UnitsDomain extentsXmax [in]E4400001Ymin (in534Ymax(in) 64scaleUnscaleCloseHeld)计算采川的长度单位已被吏改为 inches(英寸),此时便能正确反映网格的几何尺寸注意此算例的求解过程中,除了长度外,其它单位均采用SI制。一般来说,没有必要对其它单位进行改动。按照上述的操作,长度单位已被确定为 inches若用户想采用别的单位制作为长度单位,如mm,可以在 Define的下拉菜单中打开 Set units对话框,进行更改。ww.myCFDcn《数值计算与工程仿真》专刊一 FLUENT HELP算例精选中文版(→)5.显示网格。(图1.2)Display→Grid,,Grid DisplOptionsEdge TypeSurfaces三彐p Mores今Atermal-3pressure-outlet-7H Edgesv feature velocity-inlet-5p facesy outlineuelocily-inlet-rall-4Partitionswall-iiShrink factor到终安Surface Types且彐Surface eame patternclip-surfHatchfanOutline InteriorDIsplayColorsCloseHelpa)确保在 surfaces项屮的所有表面都被选屮,然后点击 Display。ww.myCFDcn
- 2020-11-30下载
- 积分：1
