计算机在材料学中综合作业.docx

《计算机在材料学中综合作业.docx》由会员分享，可在线阅读，更多相关《计算机在材料学中综合作业.docx（38页珍藏版）》请在三一办公上搜索。

1、计算机在材料学中的应用综合作业目 录一、传热学模型建立- 2 -（一）【模型一】平板焊接过程的二维温度场的数值模拟- 2 -（二）【模型二】立方体淬火过程三维非稳态温度场的有限差分计算- 5 -（三）【模型三】平板激光淬火过程的三维温度场数值模拟- 9 -（四）【模型四】长方体正火过程三维非稳态温度场的有限差分计算- 14 -（五）【模型五】端淬实验的数值模拟- 18 -三、计算机模拟计算结果处理- 26 -四、团队工作情况- 29 -五、成员投入程度评定排序- 30 -六、总结与体会- 31 -一、传热学模型建立（一）【模型一】平板焊接过程的二维温度场的数值模拟1.模型的建立二维模型的假设：

2、 1. 材料表面热物性参数不随温度变化。 2. 不考虑相变潜热。 3. 考虑工件的辐射与空气的对流换热。 4. 焊接热源能量分布为高斯分布。 5. 材料各向同性。 6. 工件为二维有限大物体。（Z方向厚度不计）2.热传导方程采用直角坐标系，则材料内部热传导方程为：式中，为材料热扩散率，为密度，为比热，为热传导率，T为温度，t为时间。材料的热物性参数，比热，热传导系数和密度均不随温度变化。3.边界条件上表面： 其它表面：式中，T是工件表面的温度，Ta是环境温度。n是其它表面的外法线方向。h是材料表面总的换热系数，包括空气对流和热辐射换热。h=+，为对流换热系数，为辐射换热系数。q(x,y,t)为

3、焊接热源能量分布函数P为热源有效功率，为热源半径，v为热源运动速度。4.有限差分方程由于对称性，取工件整体区域的一半求解。X方向节点数为N1，Y方向节点数为N2。X方向步长为x，Y方向步长为y，时间步长为t。热源沿X方向运动。内部节点差分方程：令，则内部节点差分方程变为：四条线边界节点的差分方程：四个角边界节点的差分方程：（二）【模型二】立方体淬火过程三维非稳态温度场的有限差分计算1.模型的建立钢的淬火是将钢加热到临界温度Ac3（亚共析钢）或Ac1（过共析钢）以上温度，保温一段时间，使之全部或部分奥氏体化，然后以大于临界冷却速度的冷速快冷到Ms以下（或Ms附近等温）进行马氏体（或贝氏体）转变的

4、热处理工艺，以大幅提高钢的刚性、硬度、耐磨性、疲劳强度以及韧性等，从而满足各种机械零件和工具的不同使用要求。也可以通过淬火满足某些特种钢材的铁磁性、耐蚀性等特殊的物理、化学性能。淬火过程是一个涉及到相变，热传导，热对流，热辐射的三维非稳态传热问题。为计算淬火过程的温度场，首先应建立淬火过程的传热学模型。本模型中采用以下几点假设。 1） 材料的热物性参数不随温度变化 2） 不考虑相变潜热3） 材料为各向同性4） 考虑工件辐射与空气对流换热。5） 工件为立方体2.热传导方程根据工件形状，采用三维直角坐标系，在工件内部的热传导方程为，式中，为材料的热扩散率，为热传导率，为材料密度，为比热，T为温度，

5、t为时间。 , 均不随温度变化。 3.边界条件 式中，T为工件表面温度，为淬火介质温度，n为其他表面的的外法线方向，h为淬火介质的对流换热系数。4.初始条件 初始时刻工件整体温度分布均匀，即，为一常数，取淬火加热温度，对35CrMo钢，取860。5.有限差分方程首先对求解区域离散化，求解区域为立方体，在X,Y,Z方向用网格来离散求解区域， X方向节点数为N1，Y方向节点数为N2，Z方向节点数为N3。X方向步长为X，Y方向步长为Y，Z方向步长为Z。时间步长为t。内部某节点编号为（i,j,k).左右相邻两节点为(i,j-1,k),(i,j+1,k),前后相邻两节点为(i.j,k-1),(i,j,k

6、+1).上下相邻两节点为（i,j-1,k),(i,j+1,k)。 内部节点差分方程为令，。 则内部节点差分方程变成了。 同理，六个面边界节点的差分方程分别为： 十二条棱上边界节点的差分方程：八个角边界节点的差分方程为：（三）【模型三】平板激光淬火过程的三维温度场数值模拟1.建立传热学模型图1是平板类工件激光热处理示意图。在平板类工件激光热处理过程的三维传热学模型中采用以下几点假设：（1）材料表面对激光的吸收系数不随温度变化；（2）材料的热物性参数随温度变化；（3）考虑相变潜热；（4）考虑工件的辐射与空气对流换热；（5）入射激光束能量分布为高斯分布（TEM00）；（6）激光功率恒定；（7）激光束

7、扫描速度不变；（8）材料各向同性；（9）工件为三维有限大物体；（10）工件初始温度恒定。 图1 平板类工件激光热处理示2.热传导方程P.S.如果激光束扫描的速度很快，则沿激光束扫描方向（Y方向）的传热可以忽略，这样三维传热问题就可简化为X方向和Z方向传热的二维传热问题，这就是平板类工件激光热处理过程的二维传热模型。在直角坐标系中，工件内部的热传导方程为：式中，为材料密度，Cp为材料比热，为热传导系数，T为温度，t为时间。3.边界条件边界条件为：上表面：其他表面：Q（x,y,t）是激光光斑能量分布函数。P为激光功率，A为吸收系数，R为激光光斑直径，V为激光光斑运动速度。Ta是环境温度，h是材料表

8、面总的换热系数，包括空气对流和热辐射换热。n是其他表面的外法线方向。4.初始条件初始条件为：初始时刻工件整体温度分布均匀，即，Ta为一常数。5.相变潜热的处理激光热处理加热时发生奥氏体相变，冷却时发生马氏体相变，其他相变不可能发生。目前对于在激光相变硬化等超快速加热条件下的奥氏体转变动力学尚不明白，因此，采用把奥氏体化的相变潜热折算成温度降低值的办法，把温度降低值等分，分步加到计算得到的温度场中。HA为奥氏体的相变潜热，n为一整数，一般取为10，或者其他合适的整数。可以用Koistinen-Marburger公式来计算钢铁材料中马氏体相变的转变量Fm为马氏体转变的体积百分数，Ms为马氏体开始转

9、变温度，T为节点温度，k和n为两个常数，与具体的钢种有关。每一步马氏体相变所放出的相变潜热可折算成温升，然后利用Tmn对温度场进行修正：Hm为马氏体转变的相变潜热。Tmn是每一步马氏体转变放出的相变潜热所引起的温升。6. 有限差分方程（激光热处理的三维差分方程）由偏微分方程和边界条件，初始条件构成传热的定解问题。用能量平衡法推导了这一问题的非均匀空间网络的三维显式差分方程。直角坐标系下差分方程： X方向节点数为101，y方向节点数为101，z方向节点数为51。X方向步长为X，y方向步长为Y，z方向步长为Z。 内部节点差分方程为令，。 则内部节点差分方程变成了。 六个面上节点的差分方程分别为：上

10、表面：下表面：左表面：右表面：前表面：后表面：棱上的节点的差分方程： 八个角上节点的差分方程为：（四）【模型四】长方体正火过程三维非稳态温度场的有限差分计算1、模型的建立钢铁材料在炉中加热到超过奥氏体转变点后，发生奥氏体转变，然后放到空气中冷却，就可得到含有珠光体的均匀组织，从而实现钢铁材料的正火处理。长方体是一种在工业界中常用的工件，他的正火过程是一个涉及到相变、热传导、热对流的三维非稳态传热过程，可用三维模型近似的描述长方体正火冷却的传热过程，在这个三维模型中采用以下几个假设：1、 工件为三维长方体2、 材料的热物性参数不随温度变化3、 不考虑相变潜热4、 考虑工件与空气的热对流换热5、

11、材料各向同性2、直角坐标系的建立及网格的划分 根据工件的形状，建立如图1所示的坐标系，将OX，OY分成100份，OZ分成50份，即节点为从（i,j,k)=(1,1,1)到（i,j,k)=(101,101,51)。此模型为长宽高等于10010050（cm）的钢铁材料，材料的密度为7.8g/ ,比热为0.5J/(g.k),导热系数为0.3。工件冷却过程为从边界向中心逐渐冷却。3、热传导方程 在工件内部的热传导方程为 ，式中， ，为材料的热扩散率，为热传导率，为材料密度，为比热，T为温度，t为时间。 , 均不随温度变化。 4、 初始条件初始时刻工件整体温度分布均匀，即 =，设为900C。设环境温度=

12、20C。5、边界条件外表面：式中，T为工件表面温度，为环境温度，n为外表面的外法线方向，h为材料表面的换热系数。h=+,为对流换热系数，为辐射换热系数。6、有限差分方程 X方向节点数为101，y方向节点数为101，z方向节点数为51。X方向步长为X，y方向步长为Y，z方向步长为Z。内部节点差分方程为令，。则内部节点差分方程变成了六个面上节点的差分方程分别为：上表面：下表面：左表面：右表面：前表面：后表面：十二条棱上的节点的差分方程：八个顶角上节点的差分方程为：（五）【模型五】端淬实验的数值模拟1、模型的建立端淬实验是将标准尺寸的端淬试样(25mm100mm)奥氏体化后，在专用设备上对其一端面喷

13、水冷却，后沿轴线方向的测出硬度-距水冷端距离的关系曲线的试验方法，是测定钢的淬透性的方法之一 。 由于模型轴向对称，计算时可以采用二维模型如图。 此次模拟实验的模型为标准尺寸（25mm100mm）的金属试样。先将式样加热到860并保温，使式样完全奥氏体化，然后将式样取出，放在端淬试验机上，在式样下端喷20的水120s。模型中采取如下假设：（1） 材料的热物性参数不随温度变化（2） 忽略相变潜热对温度场的影响（3） 不考虑轴向传热（4） 材料为各向同性材料（5） 考虑材料与空气的对流换热，采取第三类边界条件2、 网格划分 取式样延一个直径的剖面的右半部分建立如上右图的直角坐标系，将or方向分成1

14、25份，oz方向分成1000份。即从（i，k)=（1,1）到（i，k）=（N1，N2）=（126,1001）。3.热传导方程 材料内部的热传导方程是，式中为密度，为比热，为传导热，T为温度，t为时间。4.边界条件下边：；左边：0；其他边：式中，T是工件表面温度，Ta是环境温度（水温与环境温度相图）。n是其他表面的法线方向。是式样与冷却水的对流换热系数，是式样与空气的换热系数。材料的热物性参数，比热，热传导系数和密度均不随温度的变化而变化。5、 初始条件 端淬开始时，式样温度为860且温度均匀，完全奥氏体化。冷却水温度与环境温度始终为20。6、 有限差分方程 r方向的节点数为31，Z方向的节点数

15、为201，r方向的步长为，Z方向的步长为。 内部节点的差分方程为令， ，则内部节点的差分方程变为令 则：四个边上的节点的差分方程：上边：下边：左边：右边：四个角上节点的差分方程：右上角：右下角：左上角：左下角：二、计算机编程模拟程序设计程序如下：/*有限长圆柱体端淬*/*初始温度860*/*采用显式差分方程*/#include #include #define PI 3.14#define K 0.3 /导热系数#define CP 0.5 /比热#define P 7.8/密度#define H1 2.0/与淬火介质的换热系数#define H2 0.0003/与空气的换热系数#define

16、 N1 51#define N2 401double S,DS=0.001,S1; /时间与时间步长int T0=860,Ta=20; /初始温度double T1N1N2,T2N1N2;double FR,Fz,F1,F2; /傅里叶数double DR=0.025,DZ=0.025; /空间步长 void f_number() /傅里叶数的计算FR=K*DS/(P*CP*DR*DR);Fz=K*DS/(P*CP*DZ*DZ);F1=8*N1-12;F1=F1/(4*N1-5);F2=8*(N1-1);F2=F2/(4*N1-5);printf(F1=%ftF2=%fn,F1,F2);pri

17、ntf(FR=%ftFz=%fn,FR,Fz);void initial_temperature() /初始温度int i,j;for(i=0;iN1;i+)for(j=0;jN2;j+)T1ij=T0;void heat_conduction() /热传导方程int i,j;/*内部节点*/for(i=1;i(N1-1);i+)for(j=1;j(N2-1);j+)T2ij=T1ij+FR*(1.0-1.0/(2.0*i)*(T1i-1j-T1ij)+FR*(1.0+1.0/(2.0*i)*(T1i+1j-T1ij)+Fz*(T1ij-1-2*T1ij+T1ij+1);/*边上的节点*/*上

18、边*/for(i=1;i(N1-1);i+)T2iN2-1=T1iN2-1+FR*(1.0-1.0/(2.0*i)*(T1i-1N2-1-T1iN2-1)+FR*(1.0+1.0/(2.0*i)*(T1i+1N2-1-T1iN2-1)+2*Fz*(T1iN2-2-T1iN2-1)+2*Fz*DZ*H2/K*(Ta-T1iN2-1);/*下边*/for(i=1;i(N1-1);i+)T2i0=T1i0+FR*(1.0-1.0/(2.0*i)*(T1i-10-T1i0)+FR*(1.0+1.0/(2.0*i)*(T1i+10-T1i0)+2*Fz*(T1i1-T1i0)+2*Fz*DZ*H1/K*

19、(Ta-T1i0);/*左边*/for(j=1;j(N2-1);j+)T20j=T10j+Fz*(T10j-1-2*T10j+T10j+1)+4*FR*(T11j-T10j);/*右边*/for(j=1;j0;j=j-2) for(i=0;iN1;i=i+2) fprintf(fp6,%ft, T1ij); fprintf(fp6,n); fclose(fp6);int main()int i,j,m=0,p; printf(请输入淬火时间:); scanf(%lf,&S1);printf(n); f_number();printf(请判断傅里叶数，输入任意数字加回车继续n);scanf(%d

20、,&p);printf(结果写入文件中，请稍后n); initial_temperature();for(S=0;SS1;S+=DS) heat_conduction();if(m%1000=0) output_1();m=m+1;if(m%10000=0)output_2(); for(i=0;iN1;i+) for(j=0;jN2;j+) T1ij=T2ij;return 0;三、计算机模拟计算结果处理不同时刻端淬试样的温度场情况 50s 100s 150s 200s 250s 300s 温度标尺从以上图中可以看出，水冷端开始时温度降低很快，降低到100后开始减慢；试样上端开始时温度变化很

21、慢。整个的端淬过程中，水冷端温度迅速下降，低温区逐渐向上移动。端淬过程的时间温度图及分析材料的CCT曲线由cct曲线可知，材料的鼻尖温度大约为550，时间为1s。对比上面的端淬过程时间温度曲线可知，距水冷端2mm处1s时温度大约为550由此可知，材料可淬成马氏体的深度大约在2mm左右。将数值模拟结果导入Excel中，绘制如下曲线：四、团队工作情况本组组员共五人，分别是：林 鑫（201265052）；李开妍（201265108）；李 飞（201265068）；崔 岩（201265051）；范晓晴（201265031）。五个模型分别对应：【模型一】李 飞；【模型二】范晓晴；【模型三】林 鑫；【模型

22、四】李开妍；【模型五】崔 岩。经组内同学讨论，最终选取【模型五】进行编程模拟。组内分工如下：程序主要由崔岩同学完成设计并调试运行，李飞同学协助数据的后期处理，范晓晴、李开妍同学负责校对程序及结果，林鑫同学负责整合实验结果并整理排版通篇文档完成总结与体会。五、成员投入程度评定排序经组内全体同学投票表决，现将组内成员投入程度进行评定，排序如下：（ 1 ） 崔 岩 （ 2 ） 李开妍 范晓晴 （ 3 ） 林 鑫 李 飞 组员签字： 、 、 、 、 六、总结与体会通过本门课程的学习，组员对如何在材料学中应用计算机有了新的认识，同时也重新温习了C语言相关知识。以往的材料加工及处理工艺制定只能通过经验和大

23、量实验来进行，通过计算机编程进行数值模拟，大大提高了效率，节约了人力、物力、财力和研制周期。而且五人一组共同完成大作业的形式，增进了专业班级同学之间的交流与合作，有利于提高团队意识与集体的凝聚力。我们组的五名同学分工较为明确，各司其职，如此方使得本次大作业得以顺利完成。当然这其中也不乏坎坷与困难，比如建模型时的懵懵懂懂，设计C语言程序时的多次调试与修改，实验结果处理时的不尽如人意等等。程序语言的设计过程与实验结果后期处理尤为突出，C语言有一段时间不接触致使编写程序时问题不断；实验结果的后期处理是缺乏相关专业软件，使用常规软件处理的能力也有限。以上种种，有些问题通过小组内部的讨论与协商得到解决，但专业性太强的问题只能求助于老师和助教，在这里也感谢张老师的细心讲解与指导，大家的共同努力才能够较好又顺利地完成本次【计算机在材料学中的应用】综合作业。