板材数模论文1.docx
- 文档编号:26784401
- 上传时间:2023-06-22
- 格式:DOCX
- 页数:18
- 大小:22.70KB
板材数模论文1.docx
《板材数模论文1.docx》由会员分享,可在线阅读,更多相关《板材数模论文1.docx(18页珍藏版)》请在冰豆网上搜索。
板材数模论文1
数学建模
〔一〕、装箱设计问题
〔二〕、板材玻璃下料问题
组员:
日期:
板材玻璃下料问题
摘要
该问题属于优化问题中的排样问题。
排样下料问题在很多工业领域中都有广泛的应用,解决好排样问题,可以提高材料的利用率。
本文解决的是玻璃板材的最优化下料策略,不同的下料策略形成不同的线性规划模型。
在充分理解题意的基础上,以使用原材料张数最少、材料利用率最高为目标,采用逐级优化的方法,进行下料方案的筛选。
在第一题中,对每块原材料进行两个层次的切割。
首先按照零件需求量选用由大面积到小面积下料的两个方向排料优选的下料策略,成品料的长在原材料的长和宽两个方向上分别排列,求出最优解;其次采用由大面积到小面积下料中成品料的长和宽在原材料的长、宽两个方向套裁排料优选,而对每次切割的余料按同种方法再进行一次切割。
算出所需原材料的块数和利用率,求出最正确下料方案。
按照原材料的利用率,筛选出最正确的下料方案为按照零件需求量,进行几种零件的配套优选下料方案,所选方案是原材料的长对成品的宽,所求需要原材料的块数为579,利用率为95.03%。
第二题的求解以第一题相似,当有两种规格的原材料时,在第一题的基础上,通过控制第一种规格原材料的基础上,来选取两种材料的最正确组合。
求得需要规格为2100cm×1650cm的原材料447块,需要规格为2000cm×1500cm的原材料146块,共计593块,利用率为%。
此模型可以推广到更多板材排样下料领域的应用,通过逐级优化和组合原理,确定各种切割方式,然后再进行优化问题的求解。
关键词:
优化排样板材下料最优化
一·问题重述
在大型建筑工程中,需要大量使用玻璃材料,如门窗等。
在作材料预算时,需要求出原材料的张数。
已知板材玻璃原材料和下料后的成品料均为矩形。
由于玻璃材料特点,切割玻璃时,刀具只能走直线,且中间不能拐弯或停顿,即每切一刀均将玻璃板一分为二。
切割次序和方法的不同、各种规格搭配〔即下料策略〕不同,材料的消耗将不同。
工程实际需要解决如下问题,在给定一组材料规格尺寸后:
〔1〕在原材料只有一种规格的情况下〔例如长为2100cm,宽1650cm〕,给出最优下料策略,时所需要材料张数最少。
〔2〕在原材料为两种规格的情况下〔例如2100cm×1650cm和2000cm×1500cm〕,给出最优下料策略,使所需要材料张数最少,且利用率〔实际使用总面积与原材料总面积之比〕尽量高。
〔3〕下表是一些成品料及所需块数〔长×宽×块数〕,分别以一种原材料2100cm×1650cm及两种原材料规格2100cm×1650cm、2000cm×1500cm为例,分别给出
(1)和〔2〕的算法及数字结果,并给出两种情况下的利用率。
表1:
成品料规格及所需块数
序号
长×宽
块数
序号
长×宽
块数
1
865×857
98
2
857×715
98
3
804×746
196
4
857×675
28
5
857×665
28
6
804×663
224
7
804×661
308
8
804×639
84
9
804×631
56
10
804×563
224
11
804×536
196
12
804×535
392
13
804×551
392
14
865×446
98
15
762×446
196
16
715×446
98
17
680×446
224
18
675×446
28
19
667×446
28
20
655×446
84
21
647×446
56
22
667×426
308
23
580×446
224
24
552×446
196
25
551×446
392
26
527×426
392
二·模型假设
1、假设不考虑在切割板材玻璃的过程中的损耗;
2、假设不考虑人为的损耗;
3、假设不考虑切割工艺的不同;
4、假设不考虑玻璃厚度的影响;
5、假设不考虑刀片的厚度;
6、假设不考虑两种原材料的优先级及成本,只考虑原材料的利用率;
三·符号说明
符号
表示意义
m
当前原材料的长
e
当前原材料的宽
l(i)
成品零件的长,i=1,2,……,26
w(i)
成品零件的宽,i=1,2,……,26
n(i)
所需成品料的块数,i=1,2,……,26
S1
规格为2100cm×1650cm的原材料的面积
S2
规格为2000cm×1500cm的原材料的面积
当前原材料的利用率
两种原材料的利用率
si
第j种切割方式下余料面积
zhang
有两种规格原材料是,所需规格为2100cm×1650cm原材料的块数
zhang1
有两种规格原材料是,所需规格为2000cm×1500cm原材料的块数
N(i)
第i种切割方式使用的次数
U
有两种规格原材料是,所需总原材料的块数
S
所需成品料的总面积
四·问题分析
板材玻璃下料问题属于线性规划中的二维下料〔板材下料〕,玻璃最优化切割中,原材料有两个方向〔长和宽〕,在只有一种原材料的时候,令原材料的长为m,宽为l,零件的长为l(i),宽为w(i),i=1,2,……,26。
从操作方便的角度考虑,一张板材上不宜下过多的零件,但一般来说,参加套裁的零件种类越多,材料的利用率越高,在实际玻璃切割中要兼顾这两方面的情况,既要考虑操作的方便,又要考虑材料的利用率,一般我们讨论零件种数最多为3种的情况。
而且由于玻璃材料特点,切割玻璃时,刀具只能走直线,且中间不能拐弯或停顿,即每切一刀均将玻璃板一分为二.在一张板材上可排少数几种的综合套裁法。
既然原材料有长和宽两个方向,零件也有长和宽两个方向,则每个零件的长可在原材料的长和宽方向上排列,宽也可在原材料的长和宽的方向上排列,这就够成了二维下料方式的多样性,当所需下料的零件种类较多时,下料方式也就相应的比较多,这又为二维下料增加了困难,为了克服这个困难,并考虑到采取组成最优化方案的切割方式的材料利用率都不应太低,因此采取逐级优化的方法,进行优化切割方案的筛选。
对于第一题中有一种原材料的时候,有以下三种切割方案:
(1)第一层次按单一下料成品长对应原材料长,第二层次对余料切割按单一下料在长、宽两个方向套裁排料优选。
(2)第一层次按单一下料成品宽对应原材料长,第二层次对余料切割按单一下料在长、宽两个方向套裁排料优选。
(3)第一层次按单一下料在长、宽两个方向套裁排料优选,进行几种成品料配套优选,选出最正确下料方案。
对于第二题中有两种原材料的时候,在考虑原材料利用率的情况下,我们选用第一题中利用率较高的切割方案,即按照零件需求量,进行几种零件配套优选,选出最正确下料方案。
最后,根据所求的〔1〕题和〔2〕题中原材料的数量,算出每种切割方案的利用率,选出最正确切割方案。
五·模型的建立与求解
下料问题模型主要有密切相关的两部分组成。
第一部分为初始切割方式下的优化选取模型,第二部分为下料方案的优化选取。
问题一模型的建立与求解
第一层次按单一下料成品长对应原材料长,第二层次对余料切割按单一下料在长、宽两个方向套裁排料优选。
对于一张原材料上仅裁一种零件〔即单一下料〕的切割方式,每一种零件可以排出两种单一切割方式。
第一种,成品料的长在原材料长的方向上排列的下料方式;第二种,成品料的长在原材料宽的方向上排列的下料方式。
〔1〕、单一下料长对长排料优选数学模型如下:
n(i)=n(i)-INT((m/l(i))*(e/w(i))),i=1,2,……,26;
=max(Internet(e/l(i))*Internet(2100/w(i))*l(i)*w(i)/210/e,fix(e/w(i))*fix(2100/l(i))*l(i)*w(i)/210/e)>s〔i〕)
Zhang=N(i)
如第一种零件l
(1)在m方向上能排2次,w
(1)在e方向上仅能排1次,因此下第一种零件经这次切割之后所需的总数是n
(1)=98-2=96,如图5.1。
同样对剩余两块余料再进行一次类似的切割,切割的零件种类
〔2〕、由matlab编程求解可得单一下料长对长排料优选结果如下:
通过计算解得,在原材料只有一种规格2100cm×1650cm的情况下,用单一下料两个方向排料优选的方法需要原材料602块,原材料的利用率为数据%。
第一层次按单一下料成品宽对应原材料长,第二层次对余料切割按单一下料在长、宽两个方向套裁排料优选
单一下料宽对长排料优选数学模型如下:
n(i)=n(i)-int((m/w(i))*(e/l(i))),i=1,2,……,26;
=max(int(e/w(i))*int(2100/l(i))*l(i)*w(i)/210/e,int(e/l(i))*int(2100/w(i))*l(i)*w(i)/210/e)>s〔i〕)
Zhang=N(i)
如第一种零件l
(1)在e方向上能排2次,w
(1)在m方向上仅能排1次,因此下第一种零件经这次切割之后所需的总数是n
(1)=98-2=96,如图5.2。
同样对剩余两块余料再进行一次类似的切割,切割的零件种类。
单一下料在长、宽两个方向套裁排料优选结果如下:
通过计算解得,在原材料只有一种规格2100cm×1650cm的情况下,用单一下料在长、宽两个方向套裁排料优选的方法需要原材料579块,原材料的利用率为%。
第一层次按单一下料在长、宽两个方向套裁排料优选,进行几种成品料配套优选,选出最正确下料方案。
按照零件需求量,进行玻璃套裁下料方式的数学模型表示为:
n(i)=n(i)-int((m/w(i))*(e/l(i))),i=1,2,……,26;
=max(int(e/w(i))*int(2100/l(i))*l(i)*w(i)/210/e,int(e/l(i))*int(2100/w(i))*l(i)*w(i)/210/e)>s〔i〕),i=1、2、3……26
Zhang=N(i)
按照零件需求量,进行几种成品料配套优选结果如下:
通过matlab软件,对上述方程进行编程求解,求得在原材料只有一种规格2100cm×1650cm的情况下,用按照零件需求量,进行几种成品料配套优选的方法需要原材料595块,原材料的利用率为%。
所以,通过逐级优化,筛选优化切割方案比较得出,按照零件需求量,进行几种成品料配套优选的下料策略原材料的利用率最高,为%,故应采取这种下料策略,所需要原材料张数最少,为579张。
5.2问题二模型的建立与求解
〔1〕、有两种原材料的最优下料策略的数学模型:
在有两种原材料的情况下,采用与5.1.2中按照零件需求量,进行几种成品料配套优选方法同样的数学模型,只是在玻璃板材的切割方案中,增加一种原材料,单一下料宽对长排料优选数学模型如下:
当第一种原材料用了K(K≤579)张时,各种成品所需还需块数:
n(i)=n(i)-int((m/w(i))*(e/l(i))),i=1,2,……,26;
对个成品还需的块数n(i),用第二种原材料切割。
用同样方法可求得需zhang1张。
=max(S/〔K*S1+zhang1*S2〕),即可求出K和zhang1的最正确组合。
〔K不同时原材料利用率见附录1〕
〔2〕、有两种原材料的最优下料策略的结果如下:
通过用matlab软件求解,求得在有两种规格2100cm×1650cm和2000cm×1500cm原材料的情况下,用按照零件需求量,进行几种成品料配套优选的方法需要规格为2100cm×1650cm的原材料447,需要规格为2000cm×1500cm的原材料146块,总共594,原材料的利用率为%。
虽然总共需要的块数比只有一种原材料大,但规格为2000cm×1500cm的原材料面积比规格为2100cm×1650cm的原材料小,所以原材料的利用率提高了。
六·模型的结果分析
本模型按照逐级优化的方法,进行下料方案的筛选。
当只有一种原材料时,单一下料长对长排料优选的下料策略需要原材料602块,原材料的利用率为%;单一下料在宽对长排料优选的方法需要原材料579块,原材料的利用率为%;按照零件需求量,进行几种成品料配套优选的方法需要原材料595块,原材料的利用率为%。
故最优下料策略是第二种,需要规格为2100cm×1650cm的原材料579块。
在有两种规格2100cm×1650cm和2000cm×1500cm原材料的情况下,需要规格为2100cm×1650cm的原材料447块,需要规格为2000cm×1500cm的原材料146块,总共594块,原材料的利用率为%。
本文采用了两个层次切割,在一定程度提高了利用率。
七·模型的评价与补充
7.1模型的评价
显然,本文主要采用逐级优化的方法,对玻璃板材的最优下料策略进行分析计算的方法是基本成功的。
事实上,对于一般的二维下料〔板材下料〕问题,初始切割方式确实定是下料问题中至关重要的内容,如果确定了最优的切割方式,就可以把下料问题化成一个优化组合问题,编程进行计算求解。
模型的优点
通过逐级优化的方法,在分析各种不同切割方式的情况下,对下料方案进行比较筛选,可以选出最优的下料策略。
模型的缺点
虽然所得结果为最优,但利用率相对还是有点低。
7.2模型的补充
如果考虑到其他Guillotine切割的方法,如用遗传学算法。
遗传学算法有几个优点,建模简单;计算量相对的小,可得到近似最优的可行解;材料的利用率高。
在实际应用时,可根据具体问题的特点,先用基因遗传算法确定较有前途的假设干下料方法,再自动借用常用模型及其算法,确定以这几种方法为基础最优下料方案。
去两种模型之长,避其之短,可能会得到更好的结果。
八·模型的改良和推广
优化下料,提高材料的利用率,是国内外非常活跃的研究课题。
而板材玻璃的下料问题是一类最常见的二维下料问题,就是把一定数量的大矩形板切割成所需要的各种规格大小和数量的小矩形,钢板,木板等的剪切下料也属于此类下料问题。
因此,本模型可以推广到钢板,木板等的下料问题。
参考文献
王秋萍.数学建模讲义[J].西安理工大学,2011.
附录
1、第一种材料为1:
579时原材料利用率
2、问题一中最优策略,原材料的长对成品的宽的程序:
clear;
l=zeros(1,27);w=zeros(1,27);l
(1)=865;l
(2)=857;l(3)=804;
l(4)=857;l(5)=857;l(6)=804;l(7)=804;l(8)=804;l(9)=804;l(10)=804;l(11)=804;
l(12)=804;l(13)=804;l(14)=865;l(15)=762;l(16)=715;l(17)=680;l(18)=675;
l(19)=667;l(20)=655;l(21)=647;l(22)=667;l(23)=580;l(24)=552;l(25)=551;
l(26)=527;
w
(1)=857;w
(2)=715;w(3)=746;w(4)=675;w(5)=665;w(6)=663;w(7)=661;
w(8)=639;w(9)=631;w(10)=563;w(11)=536;w(12)=535;w(13)=551;w(14)=446;
w(15)=446;w(16)=446;w(17)=446;w(18)=446;w(19)=446;w(20)=446;
w(21)=446;w(22)=426;w(23)=446;w(24)=446;w(25)=446;w(26)=426;
zhang=0;n=zeros(26,1);n
(1)=98;n
(2)=98;n(3)=196;n(4)=28;n(5)=28;n(6)=224;
n(7)=308;n(8)=84;n(9)=56;n(10)=224;n(11)=196;n(12)=392;n(13)=392;n(14)=98;
n(15)=196;n(16)=98;n(17)=224;n(18)=28;n(19)=28;n(20)=84;n(21)=56;
n(22)=308;n(23)=224;n(24)=196;n(25)=26;n(26)=392;
zhang1=0;
s1=0;
s2=0;
s3=0;
s4=0;
temp2=0;
kk=0;
forj=1:
26
m=2100-w(j)*fix(2100/w(j));
e=1650-l(j)*fix(1650/l(j));
whilen(j)>0
n(j)=n(j)-fix(2100/w(j))*fix(1650/l(j));
zhang=zhang+1;
kk=0;
temp1=0;
temp2=0;
fori=j:
26
if(e>w(i)&n(i)>fix(e/l(i))*fix(2100/w(i))&n(i)>fix(e/w(i))*fix(2100/l(i))&n(i)>0)
if(max(fix(e/l(i))*fix(2100/w(i))*l(i)*w(i)/210/e,fix(e/w(i))*fix(2100/l(i))*l(i)*w(i)/210/e)>s1)
s1=max(fix(e/l(i))*fix(2100/w(i))*l(i)*w(i)/210/e,fix(e/w(i))*fix(2100/l(i))*l(i)*w(i)/210/e);
temp1=i;
kk=1;
end
end
if(m>w(i)&n(i)>fix(m/l(i))*fix((1650-e)/w(i))&n(i)>fix(m/w(i))*fix((1650-e)/l(i))&n(i)>0)
if(max(fix(m/l(i))*fix((1650-e)/w(i))*l(i)*w(i)/((1650-e)*m),fix(m/w(i))*fix((1650-e)/l(i))*l(i)*w(i)/(1650-e)/m)>s2)
s2=max(fix(m/l(i))*fix((1650-e)/w(i))*l(i)*w(i)/((1650-e)*m),fix(m/w(i))*fix((1650-e)/l(i))*l(i)*w(i)/(1650-e)/m);
temp2=i;
kk=1;
end
end
end
s1=0;
s2=0;
ifkk~=0
if(temp1~=0)
n(temp1)=n(temp1)-max(fix(e/l(temp1))*fix(2100/w(temp1)),fix(e/w(temp1))*fix(2100/l(temp1)));
end
if(temp2~=0)
n(temp2)=n(temp2)-max(fix(m/l(temp2))*fix((1650-e)/w(temp2)),fix(m/w(temp2))*fix((1650-e)/l(temp2)));
end
end
end
end
3.问题二程序:
clear;
l=zeros(1,27);w=zeros(1,27);l
(1)=865;l
(2)=857;l(3)=804;
l(4)=857;l(5)=857;l(6)=804;l(7)=804;l(8)=804;l(9)=804;l(10)=804;l(11)=804;
l(12)=804;l(13)=804;l(14)=865;l(15)=762;l(16)=715;l(17)=680;l(18)=675;
l(19)=667;l(20)=655;l(21)=647;l(22)=667;l(23)=580;l(24)=552;l(25)=551;
l(26)=527;
w
(1)=857;w
(2)=715;w(3)=746;w(4)=675;w(5)=665;w(6)=663;w(7)=661;
w(8)=639;w(9)=631;w(10)=563;w(11)=536;w(12)=535;w(13)=551;w(14)=446;
w(15)=446;w(16)=446;w(17)=446;w(18)=446;w(19)=446;w(20)=446;
w(21)=446;w(22)=426;w(23)=446;w(24)=446;w(25)=446;w(26)=426;
zhang=0;n=zeros(26,1);n
(1)=98;n
(2)=98;n(3)=196;n(4)=28;n(5)=28;n(6)=224;
n(7)=308;n(8)=84;n(9)=56;n(10)=224;n(11)=196;n(12)=392;n(13)=392;n(14)=98;
n(15)=196;n(16)=98;n(17)=224;n(18)=28;n(19)=28;n(20)=84;n(21)=56;
n(22)=308;n(23)=224;n(24)=196;n(25)=26;n(26)=392;
zhang1=0;
s1=0;
s2=0;
s3=0;
s4=0;
temp2=0;
kk=0;
forj=1:
26
m=q-w(j)*fix(q/w(j));
e=p-l(j)*fix(p/l(j));
while(n(j)>0&zhang>0)
n(j)=n(j)-fix(q/w(j))*fix(p/l(j));
zhang=zhang-1;
kk=0;
temp1=0;
temp2=0;
fori=j:
26
if(e>w(i)&n(i)>fix(e/l(i))*fix(q/w(i))&n(i)>fix(e/w(i))*fix(q/l(i))&n(i)>0)
if(max(fix(e/l(i))*fix(q/w(i))*l(i)*w(i)/210/e,fix(e/w(i))*fix(q/l(i))*l(i)*w(i)/210/e)>s1)
s1=max(fix(e/l(i))*fix(q/w(i))*l(i)*w(i)/210/e,fix(e/w(i))*fix(q/l(i))*l(i)*w(i)/210/e);
temp1=i;
kk=1;
end
end
if(m>w(i)&n(i)>fix(m/l(i))*fix((p-e)/w(i))&n(i)>fix(m/w(i))*fix((p-e)/l(i))&n(i)>0)
if(max(fix(m/l(i))*fix((p-e)/w(i))*l(i)*w(i)/((p-e)*m),fix(m/w(i))*fix((p-e)/l(i))*l(i)*w(i)/(p-e)/m)>s2)
s2=max(fix(m/l(i))*fix((p-e)/w(i))*l(i)*w(i)/((p-e)*m),fix(m/w(i))*fix((p-e)/l(i))*l(i)*w(i)/(p-e)/m);
temp2=i;
kk=1;
end
end
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 板材 数模 论文