易拉罐下料问题: 用数学建模优化生产管理问题

最新推荐文章于 2023-01-19 00:30:00 发布

wamg潇潇

最新推荐文章于 2023-01-19 00:30:00 发布

阅读量6.7k

点赞数 3

分类专栏： matlab数学建模

本文链接：https://blog.csdn.net/qq_29831163/article/details/89392435

版权

matlab数学建模专栏收录该内容

109 篇文章 1253 订阅

订阅专栏

问题描述

某公司采用一套冲压设备生产一种罐装饮料的易拉罐，这种易拉罐是用镀锡板冲压制成的（参见图 2）。易拉罐为圆柱形，包括罐身、上盖和下底，罐身高 10 cm ，上盖和下底的直径均为 5cm。该公司使用两种不同规格的镀锡板原料，规格 1 的镀锡板为正方形，边长 24cm；规格 2 的镀锡板为长方形，长、宽分别为 32cm 和 28cm。由于生产设备和生产工艺的限制，对于规格 1 的镀锡板原料，只可以按照图 3 中的模式 1、 2 或 3 进行冲压；对于规格 2 的镀锡板原料只能按照模式 4 进行冲压。使用模式 1、2、 3、4 进行每次冲压所需要的时间分别为 1.5s、2s、1s、3s。

该工厂每周工作 40h，每周可供使用的规格 1、2 的镀锡板原料分别为 5 万张和 2 万张。目前每只易拉罐的利润为 0.10 元，原料余料损失为 0.001 元/ $cm^{2}$ (如果周末有罐身、上盖或下底不能配套组装成易拉罐出售，也看做是原料余料损失）。工厂应如何安排每周的生产？

1 问题分析

与钢管下料问题不同的是，这里的切割模式已经确定，只需计算各种模式下的余料损失。已知上盖和下底的直径 d=5 cm，可得其面积为 $S_{1}=\pi d^{2}/4\approx 19.6\: cm^{2}$ 。周长为 $L=\pi d\approx 15.7\: cm$ ；已知罐身高 h =10 cm，可得其面积为 $S_{2}=hL\approx 157.1\: cm^{2}$ 。于是模式 1 下的余料损失为 $24^{2}-9S_{1}-S_{2}\approx 242.2$ 。同理计算其它模式下的余料损失，
并可将 4 种冲压模式的特征归纳如表 4

问题的目标显然应是易拉罐的利润扣除原料余料损失后的净利润最大，约束条件除每周工作时间和原料数量外，还要考虑罐身和底、盖的配套组装。

2 模型建立

决策变量：

用 $x_{i}$ 表示按照第i种模式的冲压次数（ i=1,2,3,4 ）， $y_{1}$ 表示一周生产的易拉罐个数。为计算不能配套组装的罐身和底、盖造成的原料损失，用 $y_{2}$ 表示不配套的罐身个数， $y_{3}$ 表示不配套的底、盖个数。虽然实际上 $x_{i}$ 和 $y_{1},\: y_{2},\: y_{3}$ 应该是整数，但是由于生产量相当大，可以把它们看成是实数，从而用线性规划模型处理。我们的计量单位是万。

决策目标：

假设每周生产的易拉罐能够全部售出，公司每周的销售利润是 0.1 $y_{1}$ 。原料余料损失包括两部分：4 种冲压模式下的余料损失，和不配套的罐身和底、盖造成的原料损失。按照前面的计算及表 4 的结果，总损失为

$0.001\left ( 242.2x_{1} +202.9x_{2}+340.4x_{3}+189.1x_{4}+157.1y_{2}+19.6y_{3}\right )$

于是，决策目标为

$\textup{max}\left \ { \ \: \: \: 0.1y_{1}- 0.001\left ( 242.2x_{1} +202.9x_{2}+340.4x_{3}+189.1x_{4}+157.1y_{2}+19.6y_{3}\right )$ ( 1 )

约束条件

1）时间约束。每周工作时间不超过 40h=144000s=14.4（万秒），由表 4 最后一列得

$1.5x_{1}+2x_{2}+x_{3}+3x_{4}\leqslant 14.4$ ( 2 )

2）原料约束。每周可供使用的规格 1、2 的镀锡板原料分别为 5 万张和 2 万张，即

$x_{1} +x_{2}+x_{3}\leqslant 5 ;\: \: x_{4}\leqslant 2$ ( 3 )

3）配套约束。由表 4 知，一周生产的罐身个数为 $x_{1}+2x_{2}+4x_{4}$ ，一周生产的底、
盖个数为 $9x_{1}+3x_{2}+12x_{3}+4x_{4}$ ，因为应尽可能将它们配套组装成易拉罐销售。所以 $y_{1}$
满足 $y_{1}=\textup{min } \: \: \left \{ \left ( x_{1}+2x_{2}+4x_{4} \right ) ,\: \left ( 9x_{1}+3x_{2}+12x_{3}+4x_{4} \right ) /2 \right \}$ ( 4 )

这时不配套的罐身个数 $y_{2}$ 和不配套的底、盖个数 $y_{3}$ 应为

$y_{2}=x_{1}+2x_{2}+4x_{4}-y_{1}\: ;\: y_{3}=9x_{1}+3x_{2}+12x_{3}+4x_{4}-2y_{1}$ ( 5 )

式（1）～（5）就是我们得到的模型，其中式（4）式一个非线性关系，不易直接处理，但是它可以等价为以下两个线性不等式： $y_{1} \leqslant x_{1}+2x_{2}+4x_{4} \right )\: ,\: y_{1}\leqslant \left ( 9x_{1}+3x_{2}+12x_{3}+4x_{4} \right ) /2$ ( 6 )

3 模型求解

在 LINGO 程序中，我们没有必要把式（20）线性化，LINGO 能够自动线性化的。编写 LINGO 程序如下：

max=0.1*y1-0.2422*x1-0.2029*x2-0.3404*x3-0.1891*x4-0.1571*y2-0.01 96*y3; 1.5*x1+2*x2+x3+3*x4<14.4; 
x1+x2+x3<5; 
x4<4; 
y1=@smin(x1+2*x2+4*x4,(9*x1+3*x2+12*x3+4*x4)/2); 
y2=x1+2*x2+4*x4-y1; 
y3=9*x1+3*x2+12*x3+4*x4-2*y1;

计算结果略。

4 评注

下料问题的建模主要有两部分组成，一是确定下料模式，二是构造优化模型。确定下料模式尚无通用的方法，对于钢管下料这样的一维问题，当需要下料的规格不太多时，可以枚举出下料模式，建立整数线性规划模型；否则就要构造整数非线性规划模型, 而这种模型求解比较困难。本节介绍的增加约束条件的方法是将原来的可行域“割去” 一部分，但要保证剩下的可行域中仍存在原问题的最优解。而像易拉罐这样的二维问题，就要复杂多了。