2018年全国大学生数学建模竞赛A题高温服装设计（含word论文和源代码资源）

零、相关内容

2018年全国大学生数学建模竞赛的A题“高温服装设计”要求参赛者在给定的高温工作环境中，通过数学建模确定适合的专用服装设计，以确保工作时的人身安全。以下是对该问题的详细分析与解答：

1. 问题理解与模型建立

   • 参数分析：需要根据附件1提供的专用服装材料参数值，对不同层次的材料特性进行深入理解，包括导热性、厚度等对热传导的影响[¹]。
   • 热传导方程：建立数学模型时，要考虑到热通过多层材料传导的物理规律，使用相应的热传导方程来计算温度分布，并验证其准确性[¹]。
   • 实验数据处理：利用附件2中提供的实验数据，校准模型，确保模型能够准确预测假人皮肤外侧的温度变化[¹]。

2. 具体问题解决

   • 计算温度分布：针对环境温度75ºC、II层厚度6mm、IV层厚度5mm的条件，计算并生成温度分布的Excel文件，以可视化假人皮肤外侧的温度变化[¹]。
   • 最优厚度确定：当环境温度降至65ºC、IV层厚度调至5.5mm时，通过模型优化II层的厚度，确保工作时间60分钟内，假人皮肤外侧温度控制在安全范围内[¹]。
   • 双层厚度优化：将环境温度提高至80ºC，同时优化II层和IV层的厚度，以保持工作30分钟内假人皮肤外侧温度不超过47ºC，并尽可能缩短超过44ºC的时间[¹]。

3. 模型校验与优化

   • 敏感性分析：对建立的模型进行敏感性分析，判断哪些参数对温度影响较大，从而进行更精细的调整和优化[¹]。
   • 优化算法应用：可能会涉及到复杂的优化算法，如遗传算法、模拟退火算法等，以寻找最优的服装设计方案[¹]。

4. 实际应用与评估

   • 实际工作环境测试：在实验室条件下测试得到的最优解需在实际应用环境中进一步测试，确保其在真实工作场景中的有效性[¹]。
   • 成本与效率考量：在设计过程中不仅要考虑安全性，还需关注成本控制和工作效率，以满足实际生产需求[¹]。

5. 模型的扩展与改进

   • 动态环境适应性：随着工作环境温度的不断变化，模型应当能够动态调整参数，适应不同的工作条件[¹]。
   • 长期稳定性研究：除了短期内的温度控制外，还需评估长时间工作的稳定性和服装材料的耐久性[¹]。

总结来说，通过对2018年高教社杯全国大学生数学建模竞赛A题“高温服装设计”的分析，可以看出这是一个涉及热传导理论、材料科学及优化算法的复杂问题。在实际操作中，除了上述的数学建模和优化过程，还需考虑实验数据的可靠性、模型的现实适用性以及成本效益分析。

一、部分题目

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

A 题 高温作业专用服装设计

在高温环境下工作时，人们需要穿着专用服装以避免灼伤。专用服装通常由 三层织物材料构成，记为I、II、III 层，其中I 层与外界环境接触，III 层与皮肤 之间还存在空隙，将此空隙记为IV 层。 为设计专用服装，将体内温度控制在37ºC 的假人放置在实验室的高温环境 中，测量假人皮肤外侧的温度。为了降低研发成本、缩短研发周期，请你们利用 数学模型来确定假人皮肤外侧的温度变化情况，并解决以下问题：

请建立数学模型解决以下问题

问题1 专用服装材料的某些参数值由附件1 给出，对环境温度为75ºC、II 层厚 度为6 mm、IV 层厚度为5 mm、工作时间为90 分钟的情形开展实验，测量得到 假人皮肤外侧的温度（见附件2）。建立数学模型，计算温度分布，并生成温度 分布的Excel 文件（文件名为problem1.xlsx）。

问题2 当环境温度为65ºC、IV 层的厚度为5.5 mm 时，确定II 层的最优厚度， 确保工作60 分钟时，假人皮肤外侧温度不超过47ºC，且超过44ºC 的时间不超 过5 分钟。

问题3 当环境温度为80ºC 时，确定II 层和IV 层的最优厚度，确保工作30 分 钟时，假人皮肤外侧温度不超过47ºC，且超过44ºC 的时间不超过5 分钟。

二、部分论文

三、部分Matlab源代码

问题1

1 求解h1h2

完整代码 https://github.com/yan-fanyu/CUMCM-Paper-And-SourceCode
%搜索h1、h2的具体值
%实际操作发现h1敏感度不大，而h2敏感度非常大
%先h1保持不变，搜索出h2=8.374，再保持h2不变，找出h1的具体数值
clc;clear;close all;
%  [FJ2,FJ2_text]=xlsread('A.xlsx','附件2');
% T0=FJ2(:,2); %初始各点的温度，T0(1)是0时的温度
%  save FJ_2
 load FJ_2.mat
%有限差分法，隐式
arg_min=1e6;
L1=0.6;L2=6;L3=3.6;L4=5;
T_en=75;tt=5400;
h2=8.374;
for h1=116:0.01:118
[T]=qiujie(L1,L2,L3,L4,h1,h2,T_en,tt);
TT=[FJ2(:,1),T(end,:)'];%模拟人的温度变化
arg=norm(TT(:,2)-FJ2(:,2))^2;
if arg<arg_min
    arg_min=arg;
    h_1=h1;
    h_2=h2;
    TT_2=TT;
end
end
%绘图
figure('Color','w')
plot(1:5401,TT_2(:,2),'b','LineWidth',1.5);
hold on
plot(1:5401,FJ2(:,2),'r:','LineWidth',1.5);
legend('模拟值','实验值')
xlabel('时间(秒)');ylabel('温度(摄氏度)')

2 已知h1h2求解温度分布

完整代码 https://github.com/yan-fanyu/CUMCM-Paper-And-SourceCode
%已知h1、h2，求解温度分布
clc;clear;close all;
 load FJ_2.mat
%有限差分法，隐式
L1=0.6;L2=6;L3=3.6;L4=5;
T_en=75;tt=5400;
h1=116.79;h2=8.374;
[T]=qiujie(L1,L2,L3,L4,h1,h2,T_en,tt);
%绘最终温度分布图
figure('Color','w')
x=[0:length(T(:,1))-1]/10;
plot(x,T(:,end)')
hold on
text(L1,T(L1*10+1,end),'|','color','r')
text(L1+L2,T((L1+L2)*10+1,end),'|','color','r')
text(L1+L2+L3,T((L1+L2+L3)*10+1,end),'|','color','r')
text(L1+L2+L3+L4,T((L1+L2+L3+L4)*10+1,end),'|','color','r')
grid on
xlabel('L(mm)');ylabel('T(°)')
%拟合优度计算
SST=norm(FJ2(:,2)-mean(FJ2(:,2)))^2;
SSE=arg_min;
R=1-SSE/SST;
%绘三维温度分布图
t=0:5400;
figure('Color','w')
surf(x,t,T')
shading interp
xlabel('L(mm)');ylabel('t(s)');zlabel('T(°)')
axis([0 15.2 0 5400 35 80 ])
%储存温度分布
xlswrite('problem1.xlsx',T')

问题2

1 求解第二层最佳厚度

完整代码 https://github.com/yan-fanyu/CUMCM-Paper-And-SourceCode
%求解二层的最优厚度
clc;clear;close all;
L1=0.6;L3=3.6;L4=5.5;
h1=116.790;h2=8.374;
T_en=65;tt=3600;
%遍历找到L2
min_L2=25;
for L2=0.6:0.1:25
[T]=qiujie(L1,L2,L3,L4,h1,h2,T_en,tt);
TT=T(end,:);%提取人体表面温度
count=find(TT>44);%找出温度超44°的点
if TT(end,end)<=47&&length(count)<=5*60&&L2<min_L2
    min_L2=L2;
    T_re=TT;
    T_max=TT(end);
    pass=length(count);%超出时间
end
end
%绘图,皮肤随时间的温度变化
figure('Color','w')
x=0:tt;
plot(x,T_re,'b-','LineWidth',1.7)
hold on
plot(x,ones(1,length(x))*T_max,'k:','LineWidth',1.5);
xlabel('t(秒)');ylabel('T(摄氏度)')
legend('L2=17.5')
axis([0 tt 35 45])

四、完整word论文和源代码