高温作业专用服装设计

1、精选优质文档-倾情为你奉上精选优质文档-倾情为你奉上专心-专注-专业专心-专注-专业精选优质文档-倾情为你奉上专心-专注-专业高温作业专用服装设计一、问题重述1.1问题背景高温作业专用服装属于个体劳动保护专用服装类，广泛适用于冶金、电力、矿山、消防、国防等高温作业场所个人防暑降温，可保障安全和健康。由于环境温度过高导致人体表面温度升高是损害工作人员身体健康的主要原因。近年来由于生产制造业的飞速发展，新工艺、新设备的大量采用，劳动群体也日渐庞大。因此，研制一种适合高温作业人员使用专用服装是我们的首要任务。1.2问题的提出问题一：已知专用服装材料的某些参数值并且当环境温度为75C、II层厚度为6

2、mm、IV层厚度为5 mm、工作时间为90分钟时，假人皮肤外侧的温度数据表已知。通过建立数学模型，计算假人皮肤外侧温度分布，并生成温度分布的Excel文件。问题二：当环境温度为65C、IV层的厚度为5.5 mm，工作60分钟时，假人皮肤外侧温度不超过47C，且超过44C的时间不超过5分钟。求解II层的最优厚度。问题三：当环境温度为80时，工作30分钟时，假人皮肤外侧温度不超过47C，且超过44C的时间不超过5分钟，求解II层和IV层的最优厚度。二、问题分析所谓高温作业 , 就是指操作场所的热源（炉子、明火、发热设备、灼热原料和大量蒸汽的设备) , 发出强烈的热辐射及对流热。为设计专用服装，将体

3、内温度控制在37C的假人放置在实验室的高温环境中，测量假人皮肤外侧的温度。问题一的分析：分析题意，可知由于是在两个坐标方向上存在温度梯度，分别是厚度和时间，传热虽然只发生在厚度这个方向上，但是温度的变化是随厚度和时间同时发生改变。假设系统在沿材料厚度变化的方向是无热源稳态热传导，沿厚度和时间两个方向建立微分方程模型，高温作业服装的温度分布可通过求解带有正确边界条件的导热方程来确定，借助傅里叶定律和热传导公式确定传热速率，从而得到温度分布，最后利用matlab软件求解方程，得到的结果导入excel表格。问题二的分析：利用第一问建立的微分方程进行求解，将其中的厚度设置为自变量进行求解，该问题的求解

4、的过程中，已知温度随着时间的变化，现在需要反向求解其厚度，并需要求得最优厚度，可以建立非线性优化模型，再运用matlab中遗传算法的工具箱，得到问题近似最优解。问题三的分析：与问题二相比多了一个决策变量，确定层和层的最优厚度，工作时间缩短为30分钟，可以建立多目标优化模型，层和层的厚度为目标函数，利用线性加权求和法，给两个层的厚度分权，将多目标模型转换为单目标模型，通过Lingo软件求解。三、模型假设1.忽略热辐射，仅考虑热传导和对流传热；2.假设系统在材料厚度方向上是一维、稳态、无热源和热导率为常数的导热；3.假设时间是温度变化的另一个维度；4.假设人体温度恒为37；四、符号说明定义符号环境

5、温度T,1层外表面温度T1层外表面温度T2层外表面温度T3皮肤外侧温度T,4材料层的比热c材料层的密度材料层的热传导率k导热速率qx热流密度q温度场T（x，y，z）时间t材料层厚度L环境的热传导率h1五、模型的建立与求解5.1模型一对高温防护图简化如图1所示。在图中，温度为T,1的热流体通过对流向温度为T1的衣服表面 = 1 * ROMAN I层传热，并靠导热传过 = 2 * ROMAN II层和 = 3 * ROMAN III层，再通过对流从温度为T3的 = 3 * ROMAN III面传至温度为T,4的冷流体。图1 高温防护图简化图5.1.1建立微分方程模型衣服各层的温度分布可通过求解带有

6、正确边界条件的导热方程3来确定。对稳态，无内热源和热汇的情况，相应的导热方程式1为q=-kT=-kiTx+jTy+kTz （5-1）式中，为三维倒三角算子；T（x，y，z）为标量温度场。确定防护服的温度分布，就是确定T（x，y，z）的温度场。由于物体本身所具有的热量不是内热源。故对于不存在内热源、物性为常数的二维稳态热传导方程1为 ct=xtx+yty+q （5-2）对于柱体发生的二维稳态导热，其数学描述为 2Tx2+2Ty2=1aT （5-3）x=0 , -Tx+1T=0 (5-4)x=L1 , -Tx+2T=0 (5-5)y=0 , -Ty+3T=0 (5-6)y=L2 , -Tx+4T=

7、0 (5-7)用分离变量法求解，设Tx，y，=X(x)Y(y)(),得到三个常微分方程d2Xdx2+2X=0 5-8 x=0 , -dXdx+1X=0 5-9x=L1 , dXdx+2X=0 5-10 d2Ydx2+2Y=0 5-11 y=0 , -dTdy+3Y=0 5-12 y=L2 , dYdx+4T=0 (5-13) dd+a2+2=0 5-14温度场的通解为Tx，y，=i=14j=1Ci，jXixYjyexp-ai2+j2 （5-15）5.2.2微分方程的求解该导热问题的求解实质上归结为在规定的初始条件及边界条件下求解。tx,0=t00 x （5-16）t(x,)x|x=0=0 （5

8、-17） ht,-t=tx,x|x= (5-18） 根据能量平衡法2，将上述式子重新整理可得 Ti-1,j+Ti,j+1+12Ti+1,j+Ti,j-1+hxkT-3+hxkTi,j=0 (5-19)利用matlab软件编程得到假人皮肤温度的曲线拟合方程：y=-4E-07x2+0.0029x+43.8506,曲线如图2所示，从图中可看到，随着时间的增大，假人皮肤外侧的温度先短时间内迅速增大，后逐渐趋于一稳定值。图2 假人皮肤温度曲线拟合图该问题所涉及的几何形状和边界条件不适合用分析法求解。本文选择用有限差分法1进行数值求解,具体步骤如下：(1)区域的离散或子区域的划分； (2)插值函数的选择；

9、(3)方程组的建立； (4)方程组的求解。把连续的定解区域用有限个离散点构成的网格来代替，如图3所示。图3 节点网格图这些离散点称作网格的节点，把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似，把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。利用ux,t关于t的向后差商： utujn-ujn-1t关于t的二阶中心差商：2ux2uj+1n-2ujn+uj-1n(x)2，对方程进行离散。离散

10、后的方程为：ujn-ujn-1t=a2uj+1n-2ujn+uj-1n(x)2 。令：r=a2t(x)2，即，1+2rujn-ruj+1n-ruj-1n=ujn-1。可化为矩阵形式：利用矩阵求逆法，通过matlab软件完成矩阵求逆。算出为未知温度。具体数据在Excel文件problem1.xlsx中。假人皮肤外侧的温度分布如图3所示。图3 假人皮肤外侧关于时间、厚度的温度分布图5.2模型二5.2.1建立非线性优化模型Min L=T,1-T,41h1A+LKA+ LKA+ L KA +1 KAcAT2-T3 (5-20)其中：qx=T,1-T31h1A+LKA+ LKA=T,1-T21h1A+L

11、KAs.t. T,447T,444且t44且t5tmax=30min (5-30) 5.3.2模型三的求解用目标达到法求解多目标规划的计算过程，可以通过调用Matlab软件系统优化工具箱中的fgoalattain函数实现。六、模型评价与推广6.1模型的评价6.1.1模型优点关于问题一，模型是通过微分方程知识推出的，而且分析之间关系连接密切，环环相扣，合理利用所给附件反映的已知条件，所以此模型最大的优点是它的逻辑准确性，符合实际。为以后研究高温作业专业服装设计打下了坚实的基础。模型二的计算充分利用了第一问的求解方法并结合遗传算法，并对本模型二中各参数影响极小的因素进行忽略，例如忽略了热辐射、热对

12、流等方面，使问题简单化，更加深刻的理解高温作业服装的原理，更加快捷的求解问题。模型三的求解原理与问题二的方法类似，前后的求解存在紧密联系，整个数学模型更具完整性。本文所建立的模型将所要解决的3个问题联系起来，运用最少的且合理的模型实现了求解，极大减少了运算量。6.1.2模型缺点但模型也存在不足，考虑的因素有些不足，在建立模型阶段忽略热辐射，皮肤表面蒸发等因素，实际生话中穿着高温作业服装存在这些情况，这样理想化的结果跟实际有一定的出入，实际会达不到完全理想化的结果。6.2模型的推广针对题目中所给的问题，不断改进模型，并对算法进行优化，可以推广进行一般的高温作业专用服装设计。同时，耐高温阻燃防护服

13、在保护消防员等应急救援人员的安危方面发挥着越来越重要的作用。国外在防护服面料及性能检测方面开展了大量的研究,而国内在自主研发防护服面料及对防护服的性能进行综合评价方面还远远落后于国外发达国家,因此进一步加强耐高温阻燃防护服的开发及性能检测的研究势在必行。本文结合实际问题，能有效得出相关参数并根据所差生误差最大的影响因素设计出了更加合理化的标定模板，为设计高温作业服装提供了很好的依据。参考文献1谢安来. 基于热传导方程正反问题的数值微分与数值实现D.东华理工大学,2016.2马强.关于热传导问题的研究J.电子测试,2017(18):104+87.3卢琳珍. 多层热防护服装的热传递模型及参数最优决

14、定D.浙江理工大学,2017.4黄青群.数学建模中优化模型的求解方法J.科技资讯,2016,14(24):139-140.5张武.MATLAB优化工具箱及其应用J.农业网络信息,2008(08):138-139+142.附件问题一程序：A=xlsread(1.xlsx);x=A(:,1);y=A(:,2);P=polyfit(x,y,2)xi=1:600:5400;yi=polyval(P,xi)plot(x,y,:o,xi,yi,-*)functionT=fitness(length)M=0.37;%热传导率cp=2100;%热容rou=862;%密度qw=50000;%导热速率a=M/ro

15、u/cp;c=qw/M;xspan=0length;tspan=05400;ngrid=200020;%空间网格数和时间网格数T,x,t=rechuandao(a,c,xspan,tspan,ngrid);%画图x,t=meshgrid(x,t);figure(1)mesh(x,t,T);xlabel(x)ylabel(t)zlabel(T)figure(2)plot(t(:,20),T(:,20);ylabel(t)zlabel(T)T=T(2000,20)T=abs(T(2000,20)-75)functionU,x,t=rechuandao(a,c,xspan,tspan,ngrid)%

16、热传导方程：n=ngrid(1);m=ngrid(2);h=range(xspan)/(m-1);x=linspace(xspan(1),xspan(2),m);k=range(tspan)/(n-1);t=linspace(tspan(1),tspan(2),n);r=a*k/h2;s=1-2*r;s2=c*h;U=zeros(ngrid);%赋初值U(1,:)=20;%差分计算forj=2:nfori=2:m-1U(j,i)=s*U(j-1,i)+r*(U(j-1,i-1)+U(j-1,i+1);endU(j,1)=50;U(j,m)=U(j,m-1);end问题二的程序：function

17、T=fitness(length)M=0.37;%热传导率cp=2100;%热容rou=862;%密度qw=50000;%导热速率a=M/rou/cp;c=qw/M;xspan=0length;tspan=03600;ngrid=200020;%空间网格数和时间网格数T,x,t=rechuandao(a,c,xspan,tspan,ngrid);x,t=meshgrid(x,t);T=T(2000,20)T=abs(T(2000,20)-65)functionU,x,t=rechuandao(a,c,xspan,tspan,ngrid)%热传导方程：n=ngrid(1);m=ngrid(2);h=range(xspan)/(m-1);x=linspace(xspan(1),xspan(2),m);k=range(tspan)/(n-1);t=linspace(tspan(1),tspan(2),n);r=a*k/h2;s=1-2*r;s2=c*h;U=zeros(ngrid);%赋初值U(1,:)=20;%差分计算forj=2:nfori=2:m-1U(j,i)=s*U(j-1,i)