第四届“华数杯”全国大学生数学建模竞赛 A题优秀论文

没有更新完之前，专栏价格为59，更新完毕之后恢复到99. 专栏内包含2024年所有数学建模比赛思路和代码，有些重要比赛着重更新（华数杯、国赛、美赛），小比赛可能会有chatgpt4更新，只需订阅一次。有些文章没有完整代码，请到专栏内查找最新代码和思路。如果比赛结束后没有更新代码（可能会有事情来不及更新）赛后我会统一退款。

2024年第五届“华数杯”专栏地址：https://blog.csdn.net/m0_52343631/category_12482955.html?spm=1001.2014.3001.5482

目录

摘要

一、 问题重述

二、 问题分析

三、 模型假设

四、 符号说明

五、 模型准备

六、 问题一模型建立以及结果分析

七、 问题二模型建立以及结果分析

八、 问题三模型的建立以及结果分析

九、 模型评价

十、 参考文献

十一、 附件

---

2023年华数杯A题优秀论文

摘要

本文研究了平纹织物整体的热导率和单根 A 纤维的热导率之间模型，对题 目所给的表面温度分布进行插值，对于纤维材料建立了整体分析和微观最小单位 分析两种模型，通过建立热传导方程 和 边界条件 并且利用两种模型的参数进行温 度分布的计算，使得两种模型下的温度分布最接近则建立了整体热导率和单根 A 纤维热导率之间的模型。

对于问题一，根据热力学定律在厚度方向上建立 一维热传导方程 ，利用 有限 差分法 对时间和厚度进行差分，建立 向后隐式差分 格式，利用 数值积分 的方法求 出温度分布，求得整体分析的温度分布图如图 3 所示 。 根据平纹织物的几何形状 建立微观最小单位几何模型 ，根据几何形状求得单根 A 纤维和该单位中空气的 体积分布，根据比热和质量比例的关系 求得单根 A 纤维的密度和比热的乘积。 建立初始条件和边界条件，用有限差分法对方程和条件进行差分，数值积分求出

温度分布。计算两个模型下各点温度的残差平方并建立 单目标最优化模型 ，利用 黄金分割法 对单根 A 纤维的热导率进行搜索，求得最优情况下的温度分布如图 4 所示。求得最优的残差平方和为 2.2796 ，最后可以求得单根 A 纤维的热导率为 0.0564 。

对于问题二，根据微观最小单位模型建立整体织物密度比热的乘积和纤维弯 曲倾角、纤维直径、纤维尺寸密度之间的关系，带入空气的比热、密度和问题一 计算得到的单根 A 纤维的比热和密度的乘积建立具体的关系式。建立两种整体 分析模型和微观最小单位模型的热传导方程和边界条件，将方程和条件进行差分， 将残差平方和和整体织物的热导率作为目标函数 ，建立 双目标优化模型 求解最 小的整体导热率，建立约束条件，利用分层序列 的思想先把残差平方和作为目标， 再将整体织物的热导率作为目标，对纤维弯曲倾角、纤维直径、纤维尺寸密度进 行遍历搜索。求得最优的整体织物热导率为 0.0317 。

对于问题三，对两问建立的热传导方程更改厚度为 0 时的边界条件，利用 牛 顿冷却定律 建立织物表层和空气接触面的边界条件，利用问题一和问题二的模型 建立部分参数的计算方法，利用有限差分法将边界条件离散化，建立向后隐式差 分格式，利用数值积分的方法对热传导方程进行求解，新问题一求解得到的最优 残差平方和为 47.6978 ，单根纤维 A 的热导率为 0.8335 。新问题二求解得 到的最优整体织物热导率为 0.0306 。

关键词： 热传导方程、有限差分法、最优化模型、黄金分割法、分层序列法

一、 问题重述

1.1 、 问题背景

随着科学技术的快速发展，航天、军工、石化、建筑等多个行业都急需良好

的隔热材料。隔热材料通过其具有绝缘性能、对热流可起到屏蔽作用的性质，从

而具有导热系数小的特点，可以用隔热材料来防止热工设备及管道的热量散失也

能对于用于家居打造冬冷夏凉的家居生活环境。

现有一种新型隔热材料 A ，用其编织形成一种平纹织物，对于平纹织物进行

热传导考虑，要求在考虑单根织物的以及对于整个织物的热传导进行考虑，在已

知纤维弯曲角度、织物厚度、经密等不同时，会对于织物的导热性能进行影响，

由此解决下述三个问题。

1.2 、 提出问题

问题一：在已知纤维弯曲角度、织物厚度等方面将会对整个织物的热传导造

成影响的前提下，附件一已有热源侧织物的表面温度，现要求建立平纹织物整体

热传导率与单根纤维热传导率之间关系的数学模型，并且用附件 2 的实验样本参

数测得平纹织物的整体热传导率为

，要求使用建立的平纹织物整

体热传导率与单根纤维热传导率之间关系的数学模型来得出单根 A 纤维的热导

率。

问题二：假设：制成织物的任意单根 A 纤维的直径在

并且

织物位于热源一侧表面温度随时间的变化的数据依旧参考附件 1 。由于温度和织

物结构造成的织物整体密度和比热的变化可以忽略。现在要求使用合适的方法选

用单根 A 纤维的直径及调整织物的经密、纬密、弯曲角度，使得织物的整体热

导率最低。

问题三：现对于附件 1 进行修改，现认为附件 1 的温度实际是热源侧织物表

面空气的温度，此时需要考虑该侧发生对流换热，假设织物表面的对流换热系数

为

，要求重新解答问题一和问题二。

二、 问题分析

本文是对于隔热材料进行考虑，即要求使用热力学相关知识来解决材料问题，

在已知附件一的温度下考虑织物的热导率以及单根纤维的热导率之间关系。

2.1 、 问题一

对于问题一，在已知热源侧织物的表面温度要求建立平纹织物整体热导率与

单根纤维热导率之间关系的数学模型 [1] 。首先，根据题意热传导发生在空气之间、

空气与织物之间以及织物内部单根纤维之间，所以要想求得整体织物热导率与单

根纤维热导率之间关系的数学模型，我们需要考虑纤维传热和空隙间气体传热。

在有上述分析后，由于静态空气热导率，经向热导率以及纬向热导率大小不一样，

所以需要进行考虑，但是无论哪种材料，对于本文的作用效果都是为了传热，所

以我们可以对于本文织物空间中的温度进行考虑。我们可以分别对于整体织物进

行热传递考虑，以及对于局部纤维之间进行热传递考虑，再分别观察两种方程下

温度分布，将问题转化为优化类问题，即要对单根纤维热导率进行遍历，使得局

部考虑下对比整体考虑下，织物空间温度分布大小差值最小。温度差值最小情况

下，得出对应织物整体热扩散率与单根纤维热扩散率之间的关系。在有热扩散率

之间的关系后，由于空间既有空气分布又有纤维材料分布，导致密度分布不均匀，

所以要想求得热导率之间的关系，还需要得出密度关系，再根据比热的定义式，

计算出比热，最终得出本题所要求的热导率之间的关系模型。得到关系模型后，

再结合附件 2 的样品参数以及整体热导率，得出最后要求的单根纤维的热导率。

2.2 、 问题二

对于问题二，仍然将问题视作最优化问题，即对于优化问题最后的目标是整

体热导率最低，局部考虑热传递作用和整体考虑热传递作用结合第一问的整体热

导率与单根纤维热导率之间的数学模型。问题中纤维的直径以及织物的经密、纬

密、弯曲角度，这些因素会对密度以及比热造成影响 [2] ，从而热扩散率换算成热

导率会受到影响。结合热传导方程以及文章的约束条件对密度和比热进行考虑。

2.3 、 问题三

对于问题三进行分析。问题三对于前两问，附件 1 的温度变成热源侧织物表

面空气的温度，即此时对于热传导方程中边界条件以及初始条件发生了变化，所

以需要对于前两问的热传递模型进行改变，其余分析步骤几乎一致。

三、 模型假设

1 、为了方便计算，将圆柱体的边界看成圆弧和直线的简易几何组合。

2 、由于对于编织物来说，深度方向的热传导应远大于横向的热传导，故忽

略横向热传导。

3 、由于空气密度，比热在温度变化较小的范围内变化较小，故直接将空气

的密度和比热看成常数

四、 符号说明

五、 模型准备

5.1 、 热传导方程的推导

热传学以热力学第一定律和第二定律为基础，热量是始终从高温热源传递导

低温热源，对于本文来说，不考率外界能量的加入，从而认为本文中热能始终是

守恒的。通常情况下，对于热传递过程，主要考虑热传导以及热对流，对于热辐

射，其作用效果的幅度远小于前两者，从而在问题中可以不对其进行考虑。

热传递是一个时空的过程，即既需要对于时间进行考虑又要对于空间进行考

虑，从而要想得出热传导方程需要进行偏微分方程处理。要想得出热传导方程，

需要对于下述三个条件进行考虑：初始条件、边界条件、边界面条件。其中初始

条件和边界条件主要通过附件一的温度进行考虑，对于边界条件，由于本文中既

有流体空气又有固体纤维，所以认为本文所用的边界面条件为第三类边界条件。

热传导方程 [3] 具体推导如下所示：

传递过程中吸收的热量

将物体微元化处理，对于一个物体的微元，在空间上

范围内，时

间上

范围内，吸收的热量可以表示为：

其中 Q1为热传递过程中吸收的热量， 为比热， 为物体密度。根据微分的定义：

当自变量变化很小时，对应的因变量变化值可以用导数进行表示从而可以得出如

下表达式：

5.2 、 热扩散率和热导率

热扩散率是热导率与比热容和密度的乘积之比，根据质量与体积之比得出对

应的密度，即本文需要度弯曲角度以及厚度进行考虑得出本文中的密度。对于比

热进行求解，可以根据整体的比热容以及空气的比热容进行求解。

六、 问题一模型建立以及结果分析

6.1 、 模型的建立

根据热力学知识，可知热传递主要是通过热传导、热辐射和热对流三种情况

进行的，其中热辐射是由物体内部微观粒子在运动状态改变时所激发出来的，但

是对于本文来说热传导和热对流的作用效果远大于热辐射。热传导方程在模型准

备中已经给出。根据问题分析，我们可以知道，对于某一时刻该系统空间上的温

度分布，可以使用两种热传导方程进行考虑，可以整体对于整个织物进行热传导

考虑，可以得出对应空间上的温度分布，再对于局部材料进行热传导考虑后，也

可以得出另外的温度分布情况，而对于热传导方程，即需要热导率，或热扩散率。

对于本题由热扩散率得出热导率还需进行比热容以及密度考虑所以先建立热扩

散率之间的关系。已知可以使用两种方式得出最后的温度分布情况，即可以根据

温度分布，建立处优化模型，目标函数即是两种不同的热传导方程下温度分布的

残差平方和的最小值。

对整体织物进行热传导方程求解空间中温度分布情况：

对于整体织物在截面方向建立坐标系，如下图所示：

七、 问题二模型建立以及结果分析

7.1 、 模型的建立

对于第二问，可以使用第一问得到的信息并建立优化模型来进行求解。在第

一问中可知单根纤维与整体织物之间存在数学关系，则该关系成为第二问优化模

型的约束条件。因为热导率主要而物质的本身性质决定，而对于本文由于材料没

有改变并且热源没有改变所以可以认为第一问得出的单根纤维的热导率没有发

生变化。但是在第二问中，题目要求就单根纤维的直径及整体织物的经密、纬密、

弯曲角度进行考虑，这几个因素将会改变热导率、热扩散率关系转换中的密度的

大小，并且密度改变也会影响比热从而影响热传导方程。综上，本题主要需要对

于密度进行考虑，可以建立优化模型，然后通过改变上面四个物理量来对密度以

及比热进行遍历，得出最后的最小热导率。

由问题（ 1 ）所建立的密度和比热的参数计算模型来对第二问进行求解分析，

经度方向的截面图如上面问腿一模型建立中图 2 所示。

取相邻经纱正中间和相邻纬纱正中间的部分为一最小单位建立下述模型：

7.2 、 问题二的求解

对于问题二，问题一中已经得出单根纤维材料的热导率，材料没有发生变化，

则认为对于问题二中的单根纤维的热导率可以使用第一问得出的结果进行运算。

根据热导率与热扩散率在热传导方程中的关系表达式可以知道对于经密、纬密、

经向弯角、纬向弯角以及直径的大小都会对密度造成影响，所以需要对这五个变

量分别设置步长再进行全局搜索，在问题一所建立的整体织物热导率与单根纤维

热导率之间的数学模型约束下搜索最优解。

使用 Matlab 进行编程运算后得出最后的几组数据如下所示：

八、 问题三模型的建立以及结果分析

九、 模型评价

十、 参考文献

十一、 附件

1. clear

2. clc

3. load( 'u1.mat' )

4. load( 'ua.mat' )

5. cs=0.05*10^6;

6. pa=1.1691;

7. ca=1012;

8. ka=0.0296;

9. B=100000000;

10. ptct=1.42e+05;

11. for d=3*10^-4:6*10^-4:6*10^-4

12. for aw=100:800:900

13. for aj=100:800:900

14. for ow=0.175:0.5:0.464

15. for oj=0.175:0.5:0.464

16. for ks=0.0296:0.02:0.0564

17. vw=((1/aj-3*d*ow)/cos(ow)+3*d*ow)*pi*d^2/4;

18. vj=((1/aw-3*d*oj)/cos(oj)+3*d*oj)*pi*d^2/4;

19. vs=2*d/(aw*aj);

20. ma=(vs-vw-vj)*pa;

21. aa=ka/(ca*pa);

22. pscs=aw*aj*(ptct*(vj+vw)+ca*pa*(aw*aj/2*d-vj

vw));

23. as=ks/pscs;

24. load( 'datatu.mat' )

25. fai=datatu(:,2);

26. uc=zeros(11,641);

27. u=zeros(11,641);

28. U=zeros(11,641);

29. deltaz=1.2/1000/10;

30. deltat=0.00015625;

31. t1=0:0.02:0.1;

32. t2=0:0.01:0.1;

33. t3=0:0.005:0.1;

34. t4=0:0.0025:0.1;

35. t5=0:0.00125:0.1;

36. t6=0:0.000625:0.1;

37. t7=0:0.0003125:0.1;

38. t=0:0.00015625:0.1;

39. z=0:deltaz:(1.2/1000);

40. interfai=interp1(t1,fai,[0.01,0.03,0.05,0.07,0.09], 'spline'

);

41. fai1=[];

42. for i=1:5

43. fai1=[fai1,fai(i),interfai(i),];

44. end

45. fai1=[fai1,fai(6)];

46. interfai1=interp1(t2,fai1,[0.005 0.015,0.025,0.035,0.04

5,0.055 0.065 0.075 0.085 0.095], 'spline' );

47. fai2=[];

48. for i=1:10

49. fai2=[fai2,fai1(i),interfai1(i),];

50. end

51. fai2=[fai2,fai1(11)];

52. interfai2=interp1(t3,fai2,[0.0025 0.0075,0.0125,0.0175,

0.0225,0.0275 0.0325 0.0375 0.0425 0.0475 0.0525 0.0575 0.0625 0.0675 0.07

25 0.0775 0.0825 0.0875 0.0925 0.0975], 'spline' );

53. fai3=[];

54. for i=1:20

55. fai3=[fai3,fai2(i),interfai2(i),];

56. end

57. fai3=[fai3,fai2(21)];

58. interfai3=interp1(t4,fai3,0.00125:0.0025:0.1, 'spline' );

59. fai4=[];

60. for i=1:40

61. fai4=[fai4,fai3(i),interfai3(i),];

62. end

63. fai4=[fai4,fai3(41)];

64

65. interfai4=interp1(t5,fai4,0.000625:0.00125:0.1, 'spline' );

66. fai5=[];

67. for i=1:80

68. fai5=[fai5,fai4(i),interfai4(i),];

69. end

70. fai5=[fai5,fai4(81)];

71. interfai5=interp1(t6,fai5,0.0003125:0.000625:0.1, 'splin

e' );

72. fai6=[];

73. for i=1:160

74. fai6=[fai6,fai5(i),interfai5(i),];

75. end

76. fai6=[fai6,fai5(161)];

77. interfai6=interp1(t7,fai6,0.00015625:0.0003125:0.1, 'spl

ine' );

78. fai7=[];

79. for i=1:320

80. fai7=[fai7,fai6(i),interfai6(i),];

81. end

82. fai7=[fai7,fai6(321)];

83. uc(:,1)=25+273;

84. uc(1,:)=fai7'+273*ones(641,1);

85. v=zeros(11,641);

86. kt=0.0564;

87. v=zeros(11,641);

88. rs=as*deltat/(deltaz^2);

89. for i=1:640

90. uc(2:10,i+1)=(1-

2*rs)*uc(2:10,i)+rs*((uc(1:9,i)+uc(3:11,i))) ；

91. plot(uc(:,i));

92. uc(11,i+1)=uc(10,i+1);

93. end

94. v=v+uc;

95. u=v/6;

96. % u=0.84*uc+0.16*ua;

97. U=(u-u1).^2;

98. b=sum(sum(U));

99. if B>b

100. B=b;

101. A=at;

102. K=ks;

103. g=u;

104. D=d;

105. AW=aw;

106. AJ=aj;

107. OW=ow;

108. OJ=oj;

109. end

110. end

111. end

112. end

113. end

114. end

115. end

116. disp([ ' 寻找到的最佳位置为： ' ,num2str(K)])

117. disp([ ' 寻找到的最佳位置为 :' ,num2str(A)])

118. g=g-273*ones(11,641);

119. mesh(t,z,g)

120. title( ' 微观分析的温度分布 ' )

121. xlabel( ' 时间 ' )

122. ylabel( ' 厚度 ' )

123. zlabel( ' 温度 ' )