《易拉罐形状和尺寸的最优设计》由会员分享，可在线阅读，更多相关《易拉罐形状和尺寸的最优设计（17页珍藏版）》请在金锄头文库上搜索。
1、07级数学实验探索实验报告易拉罐的形状和尺寸的最优设计问题高袁 屠凤华 姚鹏 成治尧 统计 071 分工：高袁：编程屠凤华：数据测量和文章处理姚鹏 成治尧：分析结论摘要饮料灌装是饮料生产中十分重要的一环，饮料灌装容器的设计不仅直接关系到生 产企业的制造成本，同是也决定着饮料产品的品质和价值。理想的饮料灌装容器应能起 到以下作用：保护内在质量、免受物理损坏、使用方便、便于运输、和促进销售。在日 常生活中,我们总会买些易拉罐装的饮料和食品,殊不知,易拉罐的设计便包含了一定的 物理、数学知识。对易拉罐的设计，生产者总会考虑让它成本最低，并且功能最强。如： 设计一个体积固定为V的圆柱形易拉罐，什么样的
2、设计方案最优？首先我们根据测的一组数据得直径和高的比值接近黄金分割点。本文基于用铝材料 做成一个容积一定的圆柱形的容器用料最省问题，我们分析说明表面积最小是正圆柱体 的最优设计。再从实际情况出发，注意到罐的顶盖比其他部分都要厚，我们引入了厚度 因子a,并结合模型一的结论r:h=l:4，考虑用材料的体积SV，建立模型二，得出 a=3.再以此为基础，建立模型三:Min S=2兀 x R x H + 兀 x R2 +3兀 x r2 + 兀 x (R + r) x Jh2 + (0.3h)2 x bS.t.x (R 3 一 r 3)R=r+0.3h设定从顶盖到胖体部分的斜率为a.并代入工程生产中普遍认
3、定的斜率0.3，运用 Mathematica软件求解，得出h=4r的结论，这与我们在第一问中用游标卡尺所测得的数 据吻合对此时的SV进行求偏导数，得出极值点为h=5.36221,r=1.49597, R=3.1046,H=10.8017 问题四我们用曲面积分思想建立了模型四Min SV =冗 x R2 x H +兀 x R2 + 3兀r2 + 2m x R x (JR2 一 r2 - R2 一R2) x b1 00 1 0m11S.t V=m x R2 x H + m x R2 x h 一 x ( R2 R2 + h)3 一(十R 一 R2 )303 v 0 0得出我们设计的易拉罐H=6.54
4、 h=2.54 R=3.82直径：高度=2R： (H+h)一、问题的提出：我们只要稍加留意就会发现销量很大的饮料 （例如饮料量为355毫升的可口可乐、 青岛啤酒等） 的饮料罐（即易拉罐）的形状和尺寸几乎都是一样的。看来，这并非偶然， 这应该是某种意义下的最优设计。当然，对于单个的易拉罐来说，这种最优设计可以节 省的钱可能是很有限的，但是如果是生产几亿，甚至几十亿个易拉罐的话，可以节约的 钱就很可观了。现在就请你们小组来研究易拉罐的形状和尺寸的最优设计问题。具体说，请你们完 成以下的任务：1取一个饮料量为355毫升的易拉罐，例如355毫升的可口可乐饮料罐，测量你们认 为验证模型所需要的数据，例如
5、易拉罐各部分的直径、高度，厚度等，并把数据列 表加以说明；如果数据不是你们自己测量得到的，那么你们必须注明出处。2设易拉罐是一个正圆柱体。什么是它的最优设计？其结果是否可以合理地说明你们 所测量的易拉罐的形状和尺寸，例如说，半径和高之比，等等。3设易拉罐的中心纵断面如下图所示，即上面部分是一个正圆台，下面部分是一个正 圆柱体。什么是它的最优设计？其结果是否可以合理地说明你们所测量的易拉罐的形状和 尺寸。4利用你们对所测量的易拉罐的洞察和想象力，做出你们自己的关于易拉罐形状和尺 寸的最优设计。5用你们做本题以及以前学习和实践数学建模的亲身体验，写一篇短文，阐述什么是 数学建模、它的关键步骤，以及
6、难点。二、模型的基本假设1、材料厚度远远小于罐的半径;2、除易拉罐的上底盖外，罐通体的厚度相同;3、硬度体现在同样材料的厚度上;4、易拉罐的材料为同样的铝材料;5、易拉罐的上端卷口材料的影响很小,可不计;6 、易拉罐底端的曲面简化成为一个平面 ;7、易拉罐拉环不加考虑。三、符号说明h:易拉罐的圆台高度(cm)r:易拉罐的上顶盖半径(cm)V :易拉罐的内部体积(cm3 )b :材料的厚度(cm)SV:所用易拉罐材料的体积( cm3)a :上顶盖的厚度因子CV(r): 罐侧面所用材料的体积为( cm 3)R:圆柱体半径R:0H:上部所接球冠半径下部圆柱高四、模型的建立与求解对问题一易拉罐所需数据
7、的测量及结果:hlDIh2nil吠355毫升可口可乐易拉罐测量数据单位:毫米圆台 上表 面直 径(D1)胖体 圆柱 直径 (D2)胖体部 分 的高度 (h2)圆台 咼度 (h1)总的高 度 (h1+h2)侧面 厚度 (m1)上顶 盖厚 度(m2)卷口 咼度 (h3)下底盖 厚度(m3)160.0266.00102.0219.96122.020.1000.2940.0250.281260.0465.98101.9620.00122.040.1030.2950.0260.283360.0266.04102.0020.02121.980.1010.2950.0260.282460.0066.9810
8、2.0220.00122.000.0980.2960.0240.282559.9866.00102.0019.98122.980.0990.2940.0270.280659.9665.96102.9820.04122.000.1000.2960.0260.283759.8866.02102.0419.98121.960.0990.2940.0280.283平均值59.9866.00102.0020.00122.000.1000.2950.0260.282从上述测量数据中，胖体圆柱直径：胖体部分高度=D2:h20.647 和黄金分割点比值0.618虽然很接近。对问题二的模型建立与求解：分析：设易
9、拉罐的高为h,底面半径为r。要求饮料罐内体积V定的前提下，根 据制造过程中消耗材料的多少来判别优劣，即最优易拉罐应该具有的最小表面S,从而 求出能使易拉罐制作所用的材料最省的上顶盖的直径和从上底盖到底部的高之比。F= AbsoluteThickness 1, Line -3.1, 0, 3.1,0,3.1,12.4,-3.1,0,-3,0.1,3,0.1,3,12.1,-3,12.1,-3,0.1mygrapg = ShowGraphicF,AxesLabel-x, y,AspectRatio-Automatic, PlotRange- 0, 12.5模型由圆柱的体积公式：V=V (r,h)二
10、兀 x r2 x h，得h = V兀x r 2又易拉罐的表面积：S - S(r) 一 2兀 x r 2 + 2兀 x r x h(1)将 h=V/pir2代入(1)式得aS(r) - 2兀 x r 2 + 2v (2)r对S(r)求一阶导数得出它的极值点：dSdr得出r二哙3h=上 X 3；（竺）2兀VI4k 2 X V 3也就是说，设计成等边圆柱时，表面积最小.但是在实际中,用手摸一上顶盖就能感觉到它的硬度要比其他的材料要硬 (厚, 因 为要使劲拉),所以以上的模型,理论上是可以的,但在实际中却存在误差.假设除易拉罐的顶盖外，罐的厚度相同，记作b,硬度体现在同样材料的厚度上，并 设上顶盖的厚
11、度为ab, a这时必须考虑所用材料的体积SV对此，我们建立以下模型模型二设饮料罐的半径为r,罐的高为h.罐内体积为V. b为除顶盖外的材料的厚度. 所用材料的体积为SV, b和V是固定参数,a是待定参数. 罐侧面所用材料的体积为：CV(r,h)二兀 x (r + B)2 兀 x r2 x H + (1 + a) x B=2k x r x H x B + 2 兀 x r x (1 + a) x B 2 + H x 兀 x B 2 + 兀 x (1 + a) x B 3饮料罐顶盖所用材料的体积为 axBxk x r2饮料罐底部所用材料的体积为 Bxk x r2所以，SV和V分别为，SV(r,H)=
12、CV(r,H) + axBxk xr2=2兀 x r x h x b + 兀(1 + a) x b 3 +h x 兀 x b 2 + 2兀 x r x (1 + a) x b 2 + (1 + a) x 兀 x b x r 2因为b0， h 0其中S(r,h)是目标函数，f (r,h) -0是约束条件，V是已知的(即罐内体积一定)，即 要在体积一定的条件下，求罐的体积最小的r, h和a使得r, h和测量结果吻合.这 是一个求条件极值的问题.采用化简消元法对上述模型的求解:从f (r, h)-n x r2 x h 一 v - 0解出h - V ,n x r2代入S,则使原问题化为求b、h,使S最
13、小，也即求r，使:2vS(r,h(r) -b x +n (1 + a)2r最小.令其导数为零，求临界点:得dS = 2b x (1 + a) xn x 8r - drr 2解得临界点为r-33 (1 + a) xn因此I v=2 (1+a) x(3(ixr)h=2,r1+a=4a=33 倍.由模型一可得:=(1+a) x r即顶盖的厚度是其他材料厚度的为验证这个r确实使S达到极小。计算S的二阶导数也=4b x 2n x (1 + a) + 兰0dr 2r3因为r0，所以这个r确实使 极小.S达到局部极小，因为临界点只有一个，因此也是全局对题目中问题二进行求解:饮料罐的体积：V = h xn x r2作饮料罐用材的总面积：S(r,h)=3 兀 x r2 x b + 兀 x r2 x b + 2兀 x r x h x b=(4 兀 x r 2 + 2兀 x r x h) x b对S(r)求一阶导数，得dS =2 xnx b x drv4r - = 0nxr得到：r= 3： V xn3 4| 4兀 h=泣4-)h=4r v=355代入数据可得：h=4x二128r=4=3.05直径：高度=2r:h=0.5这个比值和问题一的解答中计算出的实际直径和高度的
