(计算机在材料科学和 与工程中的应用)第四章材料科学和 与工过程中典型物理场的数值模拟.pptVIP

Matlab PDE二维不稳态焊接热传导求解 3 网格设置 初值和误差设置 单击Solve菜单中Parameters…选项，打开Solve Parameters对话框，输入Time为0:0.5:5，u(t0)为20,其他不变。 解方程 单击 ，或者加密网格，单击 。 单击 。 简单传热学计算机分析 7)整理数据 单击Mesh-Export Mesh…输出 p e t 的数值，单击Solve-Export Solution…输出u. 回到Matlab主窗口执行下面两条命令： u1=[p,u(:,7)] %将节点坐标和其在3s时的温度组成新矩阵 u2=sortrows(u1,3) %将u1按温度值大小升序排列。 u1=[p,u(:,4)] %将节点坐标和其在1.5s时的温度组成新矩阵 u2=sortrows(u1,3) %将u1按温度值大小升序排列。 Matlab PDE二维不稳态焊接热传导求解 3 1.5s时温度场 简单传热学计算机分析 Matlab PDE二维不稳态焊接热传导求解 3 3s时温度场 简单传热学计算机分析 4.2 浓度场的数学模型及求解 4.2 浓度场的数学模型及求解 任何不均质的材料，在—定的热力学条件下，都将趋向于均匀化。 譬如，通过扩散退火可以改善因凝固带来的成分不均匀性．这是在合金中分布不均匀的溶质原子从高浓度区域向低浓度区域运动(扩散)的结果。所以固态中的扩散本质是在扩散力(浓度、电场、应力场等的梯度)作用下，原于定向、宏观的迁移。这种迁移运动的结果是使系统的化学自由能下降。 材料的扩散现象在工程中广泛存在，如压力加工时的动态恢复再结晶，双全属板的生产、焊接过程，热处理中的相变，化学热处理以及粉末冶金的烧结等。 扩散理论的研究主要方面之—是宏观规律的研究，它重点讨论扩散物质的浓度分布与时间的关系，根据不同条件建立一系列的扩散方程，并技其边界条件求解。用计算机数值计算方法代替传统的、复杂的数学物理方程对浓度场问题进行研究已成为发展的趋势。 4.2 浓度场的数学模型及求解 上坡扩散 上坡扩散是沿着浓度升高的方向进行扩散，即由低浓度向高浓度方向扩散，使浓度发生两级分化。 例如，奥氏体向珠光体转变的过程中，碳原子由浓度较低的奥氏体向浓度较高的渗碳体扩散，就是上坡扩散。 例如，将碳含量相近的碳钢（Wc=0.441%）与硅钢（Wc=0.478%）对焊在一起，在1050℃进行长时间(13d)加热后，焊接面处两侧的碳浓度分布如图。硅的存在推高了碳的活度和化学位，从而驱使碳从含硅的一侧向另一侧上坡扩散。 发生浓度低的向浓度高的方向扩散，产生成分的偏聚而不是成分的均匀化。 硅钢侧 碳钢侧 4.2 浓度场的数学模型及求解 菲克第一定律 稳定浓度场模型如图4-21和4-22所示 ，x轴上两单位面积1（A）和2（B），间距dx，面上原子浓度为CA、CB， 则平面1到平面2上原子数n1=C1dx，平面2到平面1上原子数n2=C2dx，若原子平均跳动频率f, dτ时间内跳离平面1的原子数为n1f·dτ，跳离平面2的原子数为n2fdτ。 4.2 浓度场的数学模型及求解 菲克第一定律(一定时间内，浓度不随时间变化dC/dτ=0)，单位时间内通过垂直于扩散方向的单位截面积的扩散物质流量（扩散通量）与该面积处的浓度梯度成正比。定义：组分i每单位时间通过单位面积的质量传输量正比于浓度梯度。定义式： D—扩散系数。 负号表示质量传输的方向与浓度梯度的方向相反。J—扩散通量，g/cm2*s 式中负号表明扩散通量的方向与浓度梯度方向相反。可见，只要存在浓度梯度，就会引起原子的扩散，物体的扩散系数单位：m2/s，物理意义：单位传质量相当于单位浓度梯度下的扩散传质通量。影响因素：物体的种类，物体的结构，温度、压力等。 4.2 浓度场的数学模型及求解 菲克第二定律 解决溶质浓度随时间变化的情况，即dc/dt≠0。 如图：单元模型 则单元体积中溶质积累速率为： 两个相距dx垂直x轴的平面组成的微体积，J1、J2为进入、流出两平面间的扩散通量，扩散中浓度变化为： 由Fick第一定律： 4.2 浓度场的数学模型及求解 即第二个面的扩散通量为第一个面注入的溶质与在这一段距离内溶质浓度变化引起的扩散通量之和。若D不随浓度变化，则 菲克第二定律在三维直角坐标系下的形式： 4.2 浓度场的数学模型及求解 渗碳浓度场的有限差求解 渗碳是最常用的一种化学热处理工艺，工件渗碳后表面具有高硬度和耐磨性，心部具有良好的强韧性。工件渗碳后的碳浓度分布决定了工件的性能，求解菲克第二定律偏微分方程，能预测