扩散方程解和合金元素分布均匀性.pdf

扩散方程解和合金元素分布均匀性 王崇琳 (中国科学院金属研究所沈阳 110016) 摘要：本文讨论了扩散微分方程和几种条件下的积分解，研究了用辛普生积分法计算误差函数 Fe中扩散分布计算结果，在通常烧结温度下，Mo和Ni等扩散距离仅1—3 lam，因此，若加入合金元 素粉末，其尺度应控制于此值，以保证扩散均匀性，讨论了合金分布均匀性对性能的影响。 关键词： 扩散；误差函数；合金粉末 在粉末冶金零件中需加入合金元素，以提高其机械性能，这些元素分布均匀性对性能有较大的影 响，人们希望根据文献数据预测选定烧结条件下的分布状况，从而优化工艺；另一方面，随着电镜等 微观分析技术的发展，人们能得到合金分布的实测结果，希望由此推出该合金的扩散数据。为此，本 文尝试在研究扩散数学理论基础上，用计算机解法建立这种联系。 1、扩散方程及其通解 1．1扩散微分方程 金属和合金的原子在其晶格中不断振动，当温度升高时，原子或空位会脱离其平衡位置。若物系 中存在温度场，则热量从高温向低温流动，这是物体内能量流动，即热传导现象；若物系中存在浓度 A．Fick采用类似的数学方程表示了扩散规律。即Fick定律【l】。 描述某元素扩散通量-，‘和其浓度梯度gradC‘关系是Fick第一定律，其通式为 J‘=一DgradC‘ (1) 温度的关系符合Arrhenius方程【2】 D=Aexp(-E／RT) (2) 式中：D，扩散系数，Cm2SeC～；A，频率常数，cm2sec～；E，扩散激活能，kJmol一；R，气体 常数，8．31434Jk_1tool-1【2】或1．986calk_1mol～；T，绝对温度，K。在文献【2—3】中有各种类型扩散的E 和A值，从而可以方便的计算出不同温度下的扩散系数D值。 描述各处浓度随时间变化关系的是Fick第二定律，其通式为 兰二=div(D．gradC‘l (3) 、 ’ df 常数，在一维直角坐标系中，Fick第二定律可以简写为 等=。等 ㈣ 至于用极坐标系表示在圆柱和圆球中的扩散方程可参见文献[1】。 1．2几种扩散源的通解和误差函数 人们最感兴趣的是经过扩散后的元素分布，它将影响材料的性能。根据上述微分方程，对不同扩 散源进行数学计算，可得到扩散一定时间后各点浓度变化的方程式Ⅲ。 (1)若有一平藤扩散源(Plane source)，该物质总量为M，如图1向一维X方向扩散，则有 e。似 (5) c㈧=研M initial (2)若有一半无限平厩扩散源(Extended distribution)，其初始浓度为co，向一维X方向扩散，如 图2，在经历时间f后，其浓度分布可以用上例由P点(孝=石)积分至无穷，经数学变换，7=∥2阍， 得 14dtdf (6) c(彬)2面Co…T[e_‘2 2》√2,P／-西却 数学上有一种函数称为误差函数(errorfunction)，其表达式如下 P哦={p却 (7) 1／．20 ． 差函数后，式(6)可以