2011年全国大学生数学建模竞赛D题天然肠衣搭配（含word论文和源代码资源）

一、部分题目

2011高教社杯全国大学生数学建模竞赛题目

D题 天然肠衣搭配问题

天然肠衣（以下简称肠衣）制作加工是我国的一个传统产业，出口量占世界首位。肠衣经过清洗整理后被分割成长度不等的小段（原料），进入组装工序。传统的生产方式依靠人工，边丈量原料长度边心算，将原材料按指定根数和总长度组装出成品（捆）。

原料按长度分档，通常以0.5米为一档，如：3-3.4米按3米计算，3.5米-3.9米按3.5米计算，其余的依此类推。表1是几种常见成品的规格，长度单位为米，∞表示没有上限，但实际长度小于26米。

根据以上成品和原料描述，设计一个原料搭配方案，工人根据这个方案“照方抓药”进行生产。 公司对搭配方案有以下具体要求：

问题1 对于给定的一批原料，装出的成品捆数越多越好；

问题2 对于成品捆数相同的方案，最短长度最长的成品越多，方案越好；

问题3 为提高原料使用率，总长度允许有± 0.5米的误差，总根数允许比标准少1根；

问题4 某种规格对应原料如果出现剩余，可以降级使用。如长度为14米的原料可以和长度介于7-13.5米的进行捆扎，成品属于7-13.5米的规格；

问题5 为了食品保鲜,要求在30分钟内产生方案。

请建立上述问题的数学模型，给出求解方法，并对表1、表2给出的实际数据进行求解，给出搭配方案。

二、部分论文

三、部分源代码

第一原料可最多生产捆数

model:
sets:
cailiao/1..8/:L,N;
cixu/1..14/;
links(cailiao,cixu):x;
endsets
data:
L=3 3.5 4 4.5 5 5.5 6 6.5;
N=43 59 39 41 27 28 34 21;
enddata
@for(cixu(j):@sum(links(i,j):x(i,j))<=20);
@for(cixu(j):@sum(links(i,j):x(i,j))>=19);
@for(cixu(j):@sum(links(i,j):x(i,j)*L(i))<=89.5);
@for(cixu(j):@sum(links(i,j):x(i,j)*L(i))>=88.5);
@for(cailiao(i):@sum(cixu(j):x(i,j))<=N(i));
@for(links:@gin(x));
end

第二原料可最多生产捆数

完整代码 https://github.com/yan-fanyu/CUMCM-Paper-And-SourceCode
model:
sets:
cailiao/1..14/:L,N;
cixu/1..37/;
links(cailiao,cixu):x;
endsets
data:
L=7 7.5 8 8.5 9 9.5 10 10.5 11 11.5 12 12.5 13 3.5;
N=24 24 20 25 21 23  21 18 31 23 22 59 18 25；
enddata
@for(cixu(j):@sum(links(i,j):x(i,j))<=8);
@for(cixu(j):@sum(links(i,j):x(i,j))>=7);
@for(cixu(j):@sum(links(i,j):x(i,j)*L(i))<=89.5);
@for(cixu(j):@sum(links(i,j):x(i,j)*L(i))>=88.5);
@for(cailiao(i):@sum(cixu(j):x(i,j))<=N(i));
@for(links:@gin(x));
End

第三原料可最多生产捆数

完整代码 https://github.com/yan-fanyu/CUMCM-Paper-And-SourceCode
model:
sets:
cailiao/1..24/:L,N;
cixu/1..135/;
links(cailiao,cixu):x;
endsets
data:
L=14 14.5 15 15.5	16 16.5 17	 17.5 18 18.5 19 19.5 20	 20.5 21 21.5	22	22.5 23	23.5 24 24.5 25 25.5;
N=35 29 30	42	28	42	45	49	50	64	52	63	49	35	27	16	12	2	0	6	0	0	0 1;
enddata
@for(cixu(j):@sum(links(i,j):x(i,j))<=5);
@for(cixu(j):@sum(links(i,j):x(i,j))>=4);
@for(cixu(j):@sum(links(i,j):x(i,j)*L(i))<=89.5);
@for(cixu(j):@sum(links(i,j):x(i,j)*L(i))>=88.5);
@for(cailiao(i):@sum(cixu(j):x(i,j))<=N(i));
@for(links:@gin(x));
End

原材料三剩余根数降级为规格二使用时最多捆数

完整代码 https://github.com/yan-fanyu/CUMCM-Paper-And-SourceCode
model:
sets:
cailiao/1..14/:L,N;
cailiao2/1..14/:z;
number/1..14/:w;
cixu/1..3/;
links(cailiao2,cixu):x;
endsets
data:
  w = 24 24 6 5 0 0 0 0 0 0 0 0 0 0 ;
L=7 7.5 8 8.5 9 9.5 10 10.5 11 11.5 12 12.5 13 13.5;
enddata
@sum(cailiao(i):N(i))=15;
@for(cailiao2(i):z1=w(1)+N(1));
@for(cailiao2(i):z2=w(2)+N(2));
@for(cailiao2(i):z3=w(3)+N(3));
@for(cailiao2(i):z4=w(4)+N(4));
@for(cailiao2(i):z5=w(5)+N(5));
@for(cailiao2(i):z6=w(6)+N(6));
@for(cailiao2(i):z7=w(7)+N(7));
@for(cailiao2(i):z8=w(8)+N(8));
@for(cailiao2(i):z9=w(9)+N(9));
@for(cailiao2(i):z10=w(10)+N(10));
@for(cailiao2(i):z11=w(11)+N(11));
@for(cailiao2(i):z12=w(12)+N(12));
@for(cailiao2(i):z13=w(13)+N(13));
@for(cailiao2(i):z14=w(14)+N(14));
@for(cixu(j):@sum(links(i,j):x(i,j))<=8);
@for(cixu(j):@sum(links(i,j):x(i,j))>=7);

@for(cixu(j):@sum(links(i,j):x(i,j)*L(i))<=89.5);
@for(cixu(j):@sum(links(i,j):x(i,j)*L(i))>=88.5);
@for(cailiao2(i):@sum(links(i,j):x(1,j))<=z1);
@for(cailiao2(i):@sum(links(i,j):x(2,j))<=z2);
@for(cailiao2(i):@sum(links(i,j):x(3,j))<=z3);
@for(cailiao2(i):@sum(links(i,j):x(4,j))<=z4);
@for(cailiao2(i):@sum(links(i,j):x(5,j))<=z5);
@for(cailiao2(i):@sum(links(i,j):x(6,j))<=z6);
@for(cailiao2(i):@sum(links(i,j):x(7,j))<=z7);
@for(cailiao2(i):@sum(links(i,j):x(8,j))<=z8);
@for(cailiao2(i):@sum(links(i,j):x(9,j))<=z9);
@for(cailiao2(i):@sum(links(i,j):x(10,j))<=z10);
@for(cailiao2(i):@sum(links(i,j):x(11,j))<=z11);
@for(cailiao2(i):@sum(links(i,j):x(12,j))<=z12);
@for(cailiao2(i):@sum(links(i,j):x(13,j))<=z13);
@for(cailiao2(i):@sum(links(i,j):x(14,j))<=z14);
@for(links:@gin(x));
End

原材料二剩余根数降级为规格一使用时最多捆数

完整代码 https://github.com/yan-fanyu/CUMCM-Paper-And-SourceCode
model:
sets:
cailiao/1..8/:L,N;
cailiao2/1..8/:z;
number/1..8/:w;
cixu/1..3/;
links(cailiao2,cixu):x;
endsets
data:
  w =  16 3 0 1 0 0 0 0 ;
L=3 3.5 4 4.5 5 5.5 6 6.5;
enddata
@sum(cailiao(i):N(i))=50;
@for(cailiao2(i):z1=w(1)+N(1));
@for(cailiao2(i):z2=w(2)+N(2));
@for(cailiao2(i):z3=w(3)+N(3));
@for(cailiao2(i):z4=w(4)+N(4));
@for(cailiao2(i):z5=w(5)+N(5));
@for(cailiao2(i):z6=w(6)+N(6));
@for(cailiao2(i):z7=w(7)+N(7));
@for(cailiao2(i):z8=w(8)+N(8));
@for(cixu(j):@sum(links(i,j):x(i,j))<=20);
@for(cixu(j):@sum(links(i,j):x(i,j))>=19);
@for(cixu(j):@sum(links(i,j):x(i,j)*L(i))<=89.5);
@for(cixu(j):@sum(links(i,j):x(i,j)*L(i))>=88.5);
@for(cailiao2(i):@sum(links(i,j):x(1,j))<=z1);
@for(cailiao2(i):@sum(links(i,j):x(2,j))<=z2);
@for(cailiao2(i):@sum(links(i,j):x(3,j))<=z3);
@for(cailiao2(i):@sum(links(i,j):x(4,j))<=z4);
@for(cailiao2(i):@sum(links(i,j):x(5,j))<=z5);
@for(cailiao2(i):@sum(links(i,j):x(6,j))<=z6);
@for(cailiao2(i):@sum(links(i,j):x(7,j))<=z7);
@for(cailiao2(i):@sum(links(i,j):x(8,j))<=z8);
@for(links:@gin(x));
end

四、2010-2023年相关题目

2015D众筹筑屋
2015B太阳能小屋设计
2011D 天然肠衣搭配
2011C 养老金问题
2011B 服务平台设计
2022B 无人机定位
2015B 出租车资源配置