易拉罐形状和尺寸的最优设计模型综述.docx
《易拉罐形状和尺寸的最优设计模型综述.docx》由会员分享，可在线阅读，更多相关《易拉罐形状和尺寸的最优设计模型综述.docx（16页珍藏版）》请在冰点文库上搜索。
易拉罐形状和尺寸的最优设计模型综述
易拉罐形状和尺寸的最优设计模型
查建飞郑娴雅金兰贞
（2006年获全国二等奖）
摘要：
目前，易拉罐饮料在市场上的销量很大，易拉罐的需求也是难以估计的。
而资源是有限的，因此易拉罐的最优设计是非常有必要的。
本文着重从形状和尺寸的角度分析碳酸饮料的铝质易拉罐，在容积确定的条件下以材料最省为目标建立优化模型。
首先对雪碧、可口可乐、蓝带啤酒等易拉罐容器进行测量，获取实测值。
针对易拉罐现有形状和尺寸等数据，进行综合分析，建立了逐渐改进的三个数学模型。
模型I:
把易拉罐近似地看成一个正圆柱体，在易拉罐的容积一定时，以材料最省为目
标，用求极值的方法求得易拉罐高度h与底面半径r之间的关系为力-：
】•・.2r,用实测值进行验证发现比较吻合，但还是有一定误差存在，因此进一步建立模型U进行分析。
模型U:
进一步考虑易拉罐的形状，即罐体上面部分是一个正圆台，下面部分是一个正圆柱体时，利用线性规划方法求得此时易拉罐的最优设计。
通过对模型I中的圆柱型易拉罐的对比，所得模型与实测值更加吻合。
模型川：
以材料最省为主要目的，兼顾易拉罐的舒适度进行设计，建立模型，并给出具体的设计方案。
最后结合本模型的建立过程写对数学建模的认识与数学建模过程的难点。
关键词：
最优设计形状与尺寸合适度
一、问题重述
生活中我们发现饮料量为355毫升的可口可乐、青岛啤酒等销量很大的饮料易拉罐的形状和尺寸几乎都是一样的。
这应该是某种意义下的最优设计。
当然，对于单个的易拉罐来说，这种最优设计可以节省的钱可能是很有限的，但是如果是生产几亿，甚至几十亿个易拉罐的话，可以节约的钱就很可观了。
请通过数学建模来分析上述情况并回答如下问题：
（1）取一个饮料量为355毫升的易拉罐，测量你们认为验证模型所需要的数据，并把数据列表加以说明；如果数据不是你们自己测量得到的，请注明岀处。
（2）设易拉罐是一个正圆柱体。
什么是它的最优设计？
其结果是否可以合理说明你们所测量的易拉罐的形状和尺寸。
（3）设易拉罐的中心纵断面如下图所示，即上面部分是一个正圆台，下面部分是一个正圆柱体。
什么是它的最优设计？
其结果是否可以合理地说明你们所测量的易拉罐的形状和尺寸。
图一
（4）禾I」用你们对所测量的易拉罐的洞察和想象力，做出你们自己的关于易拉罐形状和尺寸的最优设计。
（5）用你们做本题以及以前学习和实践数学建模的亲身体验，写一篇短文，阐述什么是数学建模、它的关键步骤，以及难点。
二、问题分析
要使易拉罐达到最优设计，必须满足以下条件：
（1）保证容量是足够的。
（2）材料要最节省，使生产者在保证质量的情况下，成本能降到最低。
（3）能够保证易拉罐对容器内液体和气体的压力。
（4）易拉罐是大批量生产、运输的，要避免运输过程中瓶罐之间因碰撞造成的损失，必须稳定放置两个易拉罐，才能保证安全运输。
三、模型假设
（1）研究的对象是容量为355ml的碳酸饮料的易拉罐，如可口可乐、雪碧、蓝代啤酒
（2）测量的是无变形、无损坏的易拉罐。
（3）测量的是铝制易拉罐。
b:
易拉罐侧壁的厚度
bi:
易拉罐上底厚度。
四、符号定义
r:
易拉罐柱体部分的内半径
h:
易拉罐的内高度。
H:
易拉罐的总体高度
ai:
上底内咼。
a2:
下底内咼。
hs:
上底内高。
h4:
下底拱高。
-：
:
舒适度。
S:
易拉罐的总面积。
b2:
易拉罐下底厚度。
Z:
所用材料的体积。
:
下底厚与侧壁厚的比值。
>2:
上底厚与侧壁厚的比值。
di：
上底内直径。
d2:
下底内直径。
V:
易拉罐的容量，为一固定值。
上底厚bl
图一
五、模型建立与求解
问题一
1.1测量方法
㈠测厚度：
把易拉罐切开压平，n层的易拉罐壁进行叠加，直至总厚度可达mmm则单个易拉罐
壁厚则为m/nmm同样方法对不同品种的易拉罐进行测量，取平均值.
㈡测外径：
用一条无弹力的绳子水平绕标准易拉罐的柱体部分一圈，再利用直尺测出绳子的长度。
为减少误差，用以上相同方法进行多次测量，取得平均值。
由c=2"：
r,得出工。
㈢其余的数据全由千分尺测得。
近似取值为小数点后两位数。
对易拉罐所测得的数据见下表一：
表一：
自己测量得到的易拉罐所需数据表（单位:
mm）
雪碧
蓝代
可乐
平均数
罐总高度H
122.50
122.50
122.80
122.60
上内咼②
5.20
5.00
4.80
5.30
下内高Q2
5.00
5.00
5.10
5.00
上盖厚b】
0.30
0.30
0.30
0.30
下底厚b2
0.28
0.27
0.30
0.27
侧壁厚b
0.12
0.13
0.13
0.13
罐身直径2r
60.30
60.60
60.50
60.47
上内直径2d,
57.40
57.10
57.20
57.23
下内直径ch
50.40
53.00
50.20
51.20
上外高h3
13.80
13.60
13.50
13.63
下底拱高b
9.90
9.80
10.20
9.90
问题二（模型I）
2.1模型假设:
易拉罐的形状为正圆柱体，如图
-Y・
1
F
图
2.2模型分析：
把易拉罐近似看成一个正圆柱，要求易拉罐内的体积一定时，求能使易拉罐制作所用的材料最省。
2.3模型的建立与求解：
易拉罐侧面所用材料的体积:
-hrb2-Zlhr2=2：
rbhZlhb2,
bG•二hb?
可以忽略•所以侧面材料的体积可以近似看作2-rbh
易拉罐上底面所用材料的体积:
:
f：
r2b
易拉罐下底面所用材料的体积：
一二r2b
Zr,h一二〃厲匕2F•2rh1
V二二r2h
②代入①得：
Z=b2V/［宀’一2r2［
矢口「哼宀
令^=0,解得临界点为r处:
dr
代入②得:
h,「珂2二亠
—b13
2—1\1恥1*2一八
d2
「-学1,r0
r
又因为一§.二4b2二
dr
所以当h二〉
i*2r时，Z取得最优值.
根据测量数据:
-：
＞2;〉22,贝归［心2二4,代入公式得匕二4
r
圆柱型的易拉罐的尺寸虽然与实测数据相对比较吻合，但还是有一定的误差。
通过观察发现实际中易拉罐的上部分有类似圆台的地方，这是为了减少材料还是其他目的呢？
因此建立模型。
问题三（模型n）
3.1
补充假设（如图三）
（1）
1为圆台部分的母线长。
（5）
S2表示圆柱上底面积。
（2）
hi表示圆台部分的高度。
（6）
Vi表示圆台的体积。
（3）
h2表示圆柱上线高度。
（7）
V2表示圆柱上线的体积。
（4）
Si表示圆台上底的面积。
（8）
355-Vi=355-v2即v】
图三
3.2模型分析:
此问题考虑的是把易拉罐的上部分改成圆台使得它跟原先的体积一样，求岀hbh2
使得易拉罐的用材最省。
3.3模型的建立与求解:
圆台的表面积为:
s圆台表二（兀d+兀r）J（rf+h「+ndi2aib
圆柱咼h2的表|侨积为:
其中由模型一可知：
r=2,d=0.13mm,H=122.6mm
化简得：
兀件+r#（r一djf勺（d「一r2）一2rh2（b
max二理一出・rL・&一di$-hi卜二idi-r2_2rh2b
-hichr2dir
Zlr2H-h?
1=355000
E+1/2=i坐
rVhi4
0岂di岂
利用LINGO软件求得di=30.367i6mmR=30.367i6mmhi=0mm
h2=0mm材料为752.8548mm3（见附录一）
从模型所求结果可以看出，圆台的上底半径只有6.ii毫米，这样一来易拉罐罐口
只有ii7.283平方毫米，这与我们所测易拉罐的实际尺寸相差比较大。
因此我们对所假设的模型做如下改进：
在实际中圆台上底半径与易拉罐中间半径相差很小约5；圆台的
高度也是很小约为5。
因此模型修正为：
maxZ-di'力|ir-d"'dj-idj-r2-2rh2b
ch'3h23
—+1/2=:
———
rT帀4
—di2r2二r2122.6-ho1=355000
hi
3hi=5
hi=5.000000mm
dl=5.000000mmr=30.55185mmh2=3.462191mm
料为-450.1384（mm）A3
（见附录二）
由于实际中易拉罐的圆台与罐底，方便运输，而问题三中讨论的易拉罐上部是圆台用料与圆柱的废料问题，因此，浪费点原料来提高人对易拉罐的舒适度。
本文在对问题四的回答中进一步的阐述这个改进。
问题四（模型三）
改进后易拉罐模型的母线为:
改进后易拉罐模型的上底半径:
I二
图四
h/cosv
dj=
htanr
因此，可求得改进后易拉罐模型的表面积:
4二r2hta
h
7
COS”
上下底面积之和为:
hr2十让规定改进后易拉罐模型石容嚥不变，还是"怙$+兀£+hta
H）]h3
（JL+htan「
V,则
nBf+兀r（r+hta
圆柱形的表面积：
rhb
圆柱上下底的面积:
S侧二2二
S底h/r2）2b•二尺」
因此，可得原型易拉罐的面积为：
22
s「二2-rhb+二r：
2br■rb+450
改进后易拉罐模型与原型易拉罐的面积的差值:
rh|br2
S二1/2（4田J+2兀山tanH）|■兀（r+h
nro初一450
cosFl
舒适程度的函数曲线：
2
tan日）o小一2町花b—
「•二f（R——它的图象为可以根据市场调查的散点图形式来拟合曲线，如下图:
综上所述，浪费材料与舒适程度的关系式为：
p二-spi；（k：
舒适程度的增加
钱
的增值；pi：
表示每个单位的价格）
问题五：
浅谈数学建模
在未接触建模时，就己经听说过它了，但不太了解，直到真正接触，才发现原来数学建模就在我们的身边。
早在中学，甚至是小学时就已经用建立数学模型的方法来解决过一些简单或理想化的实际问题。
例如航行问题、速度问题等。
数学建模不只是对实际问题建立数学模型的过程，它还包括了对实际问题进行的解
释、求解、验证并解决的全过程。
数学建模是运用如同MATLABLINGO等数学工具来得
到一个数学结构，用各种数学方程、表格、图形等表现模型的思想。
它是为解决现实中存在的特定的对象，特定的目的，根据特有的内有规律进行必要的简化假设，而设计的模型。
另外，对于同一个客观事物可以有多种数学描述，因此有必要在若干模型中选择一个最简单，最恰当，最易于进行数学处理的模型。
对于模型，可以选择一个最简单，最恰当，最易于进行数学处理的模型。
如下图：
实际1旬题一一A模型假设一》模型建立一•模型求辭
在易拉罐模型中，我们在审题时，结合实际中普遍饮用的易拉罐，大概了解易拉罐的构造，形状，并在某些生产产品的公司网站上查阅了一些有关易拉罐具体尺寸、材料组成等，以及它的制作工艺、过程。
考虑到罐内物质对瓶罐的压力，再查阅了有关于应力等方面的知识。
我们也通过各种方法，测量了易拉罐的直径、高度、厚度等。
其次，建立数学模型。
通过观察可以发现，饮料罐总体上可分为部分：
圆台、圆柱当然，这不是偶然，必定是某种意义上的最优设计。
对模型进行适当的假设与简化，这是建立模型的关键一步，必须对研究对象进行统一，即认定都为饮料量为355毫升的易拉罐。
为减少误差，对于测得的数据用多次测量取平均值。
根据题目要求，明确易拉罐要达到最优，必须满足两个条件：
一、材料最省。
二、能够保证易拉罐的舒服度。
在认为易拉罐是一个长方体绕轴心旋转360度的模型时，求极值问题则可利用求导方法取临界值。
当易拉罐是由圆台和圆柱构成时，可与第一个模型进行比较，从节省材料的角度来考虑，得出最优模型。
在模型中，考虑了舒适度，建立此函数关系时，需要通过顾客心理来描述图形，但是由于没有充裕的时间，整个模型过于简略，不够精确。
六、模型改进
就问题四而言，本文只是针对易拉罐的形状和尺寸进行最优设计，并没有对易拉罐所用材料进行改进。
但是铝质易拉罐的原材料需大量进口，且价格昂贵，加上污染环境等因素，国家己不再批准新建铝质易拉罐生产线。
因此能否找到一种新型材料用来代替铝质易拉罐己成为我国关心的问题。
对此，可以在易拉罐取材上进行分析，将会得到一定的效果的。
七、模型评价
7.1模型优点
（1）模型I运用导数求极值的数学方法，计算方便，简单易懂。
（2）通过利用数学工具和Lingo编程的方法，严格地对模型求解，具有科学性。
（3）建立的模型能与实际紧密联系，结合实际情况对所提出的问题进行求解，使模型更贴近实际，通用性、推广性较强。
7.2模型的缺陷
（1）由于未能查得易拉罐的具体尺寸，只能运用测量工具对易拉罐进行测量，因此所测得的测量值与实际值之间有一定的误差，
（2）模型只考虑了半径与高之比，没有考虑其他，可能有点不全面。
八、参考文献
[1]唐焕文贺明峰，数学模型引论（第三版），高等教育出版社，2005.3
[2]杨启帆方道元，数学建模，浙江大学出版社，1993.3
[3]叶其孝，最优化-导数的应用（极值问题）教学单元，http:
//wwwl.mcm.edu.cn:
8060
[4]上海联合制罐有限公司，http:
//www.suucn/index.htm
附件
model:
max二3.14*（（dl+r）*（（r-dlF2+dlA2）A（l/2）*0.13+0.26*（dlA2-rA2）-0・26*r*h2）;
dl/叶l/2=（3*h2/hl-3/4Fl/2：
1/3*3.14*（dlA2+rA2+dl*r）*hl+3.14*（122.6-h2）*rA2=355000:
h2<=122.6;
dl<=r;
end
Objective
752.8548
value:
Variable
Value
ReducedCost
Dl
30.36716
0.000000
R
30.36716
0.000000
H2
0.000000
18.65100
Hl
0.000000
6.140743
Localoptimalsolutionfoundat
iteration：
406
RowSlackorSurplusDualPrice
1752.8548
1.000000
2-0.5000000E+30
0.000000
3
0
50
.000000
.000000
0.2120718E-02
4
1
22.6000
0.000000
86.77114
附件二
model:
max=3.14*（（dl+r）*（（r-dl）A2+dlA2）A（1/2）*0.13+0.26*（dlA2-rA2）-0.26*r*h2）;
dl/r+l/2=（3*h2/hl-3/4）Al/2;
1/3*3.14*（dlA2+rA2+dl*r）*hl+3.14*（122.6-h2）*rA2=355000;
hl二5;
dl二5;
end
Localoptimalsolutionfoundatiteration：
43636
Objectivevalue:
-450.1384
Variable
Value
ReducedCost
dl
5.000000
0.000000
r
30.55185
0.000000
h2
3.462191
0.000000
hl
5.000000
0.000000
RowSlackorSurplus
DualPrice
-450.1384
1.000000
2-0.1134488E-05
94.65037
3
-0.8106054
-0.1177982E-02
4
0.000000
-25.23358
5
0.000000
-2.453684
升级会员 