1、 大连海事大学 2016 年 数学建模 竞赛 甲 (乙丙 )组 论文 配送中心选址 姓 名 学 号 学 院 专业 班级 王丽丽 2220061234 信息科学与工程学院 自动化 20061 王丽丽 2220061234 信息科学与工程学院 自动化 20061 王丽丽 2220061234 信息科学与工程学院 自动化 20061 2016 年 11 月 25 日 - 1 - 摘 要 一件产品由若干零件组装而成，标志产品性能的某个参数取决于这些零件的参数。零件参数包括标定值和容差 两部分。进行成批生产时，标定值表示一批零件该参数的平均值，容差则给出了参数偏离标定值的容许范围。若将零件参数视为随机变
2、量，则标定值代表期望值，在生产部门无特殊要求时，容差通常规定为均方差的 3 倍。 进行零件参数设计，就是要确定其标定值和容差。这时要考虑两个方面因素：一是当各零件组装成产品时，如果产品参数偏离预先设定的目标值，就会造成质量损失，偏离越大，损失越大；二是零件容差的大小决定了其制造成本，容差设计得越小，成本越高。 关键词 ： 标定值，容差，非线性规划，正态分布 - 2 - 1. 问题重述 一件产品由若干零件组装 而成，标志产品性能的某个参数取决于这些零件的参数。零件参数包括标定值和容差两部分。进行成批生产时，标定值表示一批零件该参数的平均值，容差则给出了参数偏离标定值的容许范围。若将零件参数视为随
3、机变量，则标定值代表期望值，在生产部门无特殊要求时，容差通常规定为均方差的 3 倍。 进行零件参数设计，就是要确定其标定值和容差。这时要考虑两个方面因素：一是当各零件组装成产品时，如果产品参数偏离预先设定的目标值，就会造成质量损失，偏离越大，损失越大；二是零件容差的大小决定了其制造成本，容差设计得越小，成本越高。 试通过如下的 具体问题给出一般的零件参数设计方法。 粒子分离器某参数 (记作 y)由 7 个零件的参数 (记作 x1, x2, , x7)决定，经验公式为 7616.1242/356.02485.01235136.0162.2142.174 xx xxxxxxxxxy y 的目标值
4、(记作 y0)为 1.50。当 y 偏离 y0 0.1 时，产品为次品，质量损失为 1000 元；当 y 偏离 y0 0.3 时，产品为废品，质量损失为 9000 元。 零件参数的标定值有一定的容许变化范围；容差分为 A、 B、 C 三个等级，用与标定值的相对值表示，A 等为 1%， B 等为 5%， C 等为 10%。 7 个零件参数标定值的容许范围及不 同容差等级零件的成本 (元 )如表 1(符号 / 表示无此等级零件 )所示。 现成批生产，每批产量 1000 个。在原设计中， 7 个零件参数的标定值为： x1 = 0.1， x2 = 0.3， x3 = 0.1，x4 = 0.1， x5
5、= 1.5， x6 = 16， x7 = 0.75；容差均取最便宜的等级。 请你综合考虑 y 偏离 y0 造成的损失和零件成本，重新设计零件参数 (包括标定值和容差 )，并与原设计比较，总费用降低了多少。 表 1 各零件参数标定值的容许范围及不同容差等级零件的成本 标定值容许范围 C 等 B 等 A 等 x1 0.075, 0.125 / 25 / x2 0.225, 0.375 20 50 / x 3 0.075, 0.125 20 50 200 x 4 0.075, 0.125 50 100 500 - 3 - x 5 1.125, 1.875 50 / / x 6 12, 20 10 2
6、5 100 x 7 0.5625, 0.935 / 25 100 2. 基本假设与符号约定 为了简化问题和方便讨论，除问题中给出的假设外，我们进一步做如下的假设和说明： (1) 零件在加工制造过程 中存在多种随机因素。由于每个零件的参数 xi(i = 1, , 7)是由大量相互独立的随机因素综合影响 而 形成的 ，且 每一个个别因素在总的影响中作用都很微小 ，因此，根据中心极限定理可假设 xi 是服从正态分布的随机变量，即设 xi N(xi0, i2 )， xi0 为 xi 的标定值。 (2) 假设组成产品 (粒子分离器 )的 7 个零件在生产过程中互不影响，而且这些零件可以无困难地组装成一件
7、产品。因此， 7 个零件的参数可视作相互独立的正态随机变量。 (3) 由于产品 (粒子分离器 )的参数是由 7 个零件的参数确定的，因此产品的参数也 是随机变量，记为Y，它的取值记为 y。 (4) 假设问题中的经验公式无系统偏差，在给定的零件参数变化范围内均是有效的。 (5) 假设生产过程中没有工艺失误造成的产品损坏，产品的等级 (正品、次品、废品 )仅与产品的参数偏离其目标值的程度有关 。 在此，我们也约定文中所用符号如下： xi0 第 i 个零件参数 xi 的标定值， i = 1, , 7；令 x0 = (x10, x20, , x70)T i 第 i 个零件参数 xi 的均方差， i =
8、 1, , 7 ti 第 i 个零件参数 xi 的容差， i = 1, , 7 ri 第 i 个零件参数 xi 的相对容差 (容差关于标定值 xi0 的相对值， i = 1, , 7) ai， bi 第 i 个零件参数 xi 标定值取值的下界和上界， i = 1, , 7 y0 产品 (粒子分离器 )参数的目标值 y 产品 (粒子分离器 )参数的 数学期望 y 产品 (粒子分离器 )参数的均方差 ty 产品 (粒子分离器 )参数的容差 ry 产品 (粒子分离器 )参数的 相对容差 (容差关于目标值 y0 的相对值 ) Y 描述产品 (粒子分离器 )参数的随机变量，取值为 y - 4 - dij
9、 第 i 个零件参数 xi 采用的 容差等级， i = 1, , 7， j = 1, 2, 3(分别对应 A、B、 C 三个等级 )；如 d12表示 x1 的容差等级为 B cij 第 i 个零件参数 xi 采用第 j 个 容差等级时所需成本， i = 1, , 7， j = 1, 2, 3 Ci(ri) 第 i 个零件的成本 C(d) 每件产品的总成本 L(Y) 质量损失函数 f (x1, x2, , x7) 由 7 个零件参数表示的产品参数的经验公式，也记为 f (x) 3. 问题 的分析与模型的建立 3.1 问题的分析 显然，给定的问题是一个优化问题，目标函数为产品的总费用。分析题意可知
10、 ，产品的总费用由两部分组成，一部分是零件的成本， 另 一部分是 产品的参数 y 偏离其目标值 y0 造成的质量损失 。 这两部分相互制约：零件容差等级越低，成本越低，但产品参数偏离目标值的可能性就越大， 造成的质量损失 也越大；产品参数偏离目标值越小， 造成的质量损失 越小，但要求的零件容差等级就越高，成本也越高。因此，原问题要求在零件的容差等级和产品参数偏离目标值之间找到一种平衡，使得总费用最小。 (1) 由于零件的成本只与容差有关，记第 i 个零件的成本为 Ci(ti)，其中 ti 为 第 i 个零件参数 xi 的 容差(i = 1, , 7)，则 每件产品的总成本为 C1(t1) +
11、+ C7(t7)。 根据题 设 ，容差不 取 连续值，而是分为三个等级 A、 B、 C。如果 我们分别用 1、 2、 3 表示这三个等级 ，并引入布尔变量 dij： 否则 个容差等级采用第参数,0 ,1 jiij xd， i = 1, , 7， j = 1, 2, 3 (1) 则 dij表示 第 i 个零件参数 xi 采用第 j 个 容差等级与否，并且满足如下条件 131 j ijd (i = 1, , 7)及 d11 = d13 = d21 = d51 = d52 = d73 = 0 (2) 再令 cij表示 xi 取到第 j 个 容差等级时所需的成本 (取值 见表 1， i = 1, ,
12、7， j = 1, 2, 3)，于是 每件产品的总成本是零件容差等级的函数，为 71 31)( i j ijij dcdC (3) 其中 d = (dij0)37 是满足 (2)式的 一组取值。 事实上，容差等级共有 1233132 = 108 种组合，故 d 的所有可能取值共有 108 个。 (2) 由假设 (3)，产品的参数也 是随机变量，记为 Y。令产品的质量损失函数为 L(Y)， 则 - 5 - 3.0 |,90003.0 | 1.0,10001.0 |,0)(000yYyYyYYL (4) 其中 y0 = 1.5。由于 Y 是随机变量，因此大批量生产时平均每件产品的损失费用应是损失函
13、数 L(Y)的数学期望 EL(Y)。如果分别用 p1、 p2、 p3 表示产品为正品、次品、废品的概率，即 p1 = P|Y y0 | 0.3 (5) 则 EL(Y) = 1000p2 + 9000p3 (6) 于是，生产一批 1000 件产品的总费用为 Z = 1000C(d) + EL(Y) = 1000C(d) + 1000p2 + 9000p3 (7) 然而 ， EL(Y) = 1000p2 + 9000p3， 其中涉及了产品的次品率和废品率 (随机变量 Y 的概率 )。为此，我们必须讨论 Y 的分布，确定其分布形式和概率密度函数。 附录： 1. 经 验公式的偏导数 题目中给出的经验公
14、式 y = f (x)， x = (x1, x2, , x7)T，主要为乘积和幂的形式，根据微积分知识，先对f (x) 取对数然后再求其偏导。 对 f (x) 取对数有 7616.12456.02412351lnln)()(36.0162.21ln21 )l n ( l n85.0lnln)42.174l n ()(ln23 xxxxxxxxxxxxf 并令， 56.02436.01 xx ， 16.1242/362.21 xx ，于是 121185.01)(ln xxxx xf ， 16.2216.142/36.126.042/112252.13 9 6.085.0)(ln xxxxxxx
15、xf 3385.0)(ln xx xf ， 16.1216.042/36.024.042/1452.13 9 6.0)(ln xxxxx xf 551)(ln xx xf ， 66 21)(ln xx xf ， 77 21)(ln xx xf (F1) 再由 - 6 - ii xxfxfxxf )(ln)()( 即可求得 f (x)/xi， i = 1, , 7。 2. 正态分布随机变量 参 见 参考文献 1。 参考文献 1 刘嘉焜等 ， 应用概率统计 ， 北京 ： 科学 出版社， 2004。 2 沙定国 ， 实用误差 理论与数据处理 ， 北京：北京理工大学出版社， 1993。 3 朱道元等 ， 数学建模案例精选 ， 北京 ： 科学 出版社， 2003。 4 汪国强 ， 数学建模优秀案例选编，广州：华南理工大学出版社， 1998。 5 姜启源 ， “零件的参数设计”模型和评价 ， 数学的实践与认识， 1998， 28(1)： 5457。
