数学建模在材料科学中的应用举例.doc

《数学建模在材料科学中的应用举例.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、数学建模在材料科学中的应用举例现代科学技术发展的一个重要特征是各门科学技术与数学的结合越来越紧密。数学的应用使科学技术日益精确化、定量化，科学的数学化已成为当代科学发展的一个重要趋势。数学模型是数学科学连接其他非数学学科的中介和桥梁，它从定量的角度对实际问题进行数学描述，是对实际问题进行理论分析和科学研究的有力工具。数学建模是一种具有创新性的科学方法，它将现实问题简化，抽象为一个数学问题或数学模型，然后采用适当的数学方法求解，进而对现实问题进行定量分析和研究，最终达到解决实际问题的目的。计算机技术的发展为数学模型的建立和求解

2、提供了新的舞台，极大地推动了数学向其他技术科学的渗透。 材料科学作为21世纪的重要基础科学之一，同样离不开数学。通过建立适当的数学模型对实际问题进行研究，已成为材料科学研究和应用的重要手段之一。从材料的合成、加工、性能表征到材料的应用都可以建立相应的数学模型。有关材料科学的许多研究论文都涉及到了数学模型的建立和求解，甚至产生了一门新的边缘学科——计算材料学（Computational Materials Science），正是这些数学手段才使材料研究脱离了原来的试错法（Trial or Error）研究，真正成为一门科学。以

3、下给出一些与材料科学有关的具体建模实例。例1：金属中空位形成能建模研究1）建模准备金属中空位研究的重要性，研究空位缺陷的形成能。高能粒子对材料性能的影响，尤其是反应堆的金属材料在高能粒子的辐射作用下，性质如何变化，如何保证其安全运行？固体受辐射后产生三种效应：电离、蜕变和离位(产生空位和填隙粒子)，其中离位是金属中最主要的辐照效应。空位形成的模型1－形成空位2－蒸发3－形成填隙原子2）建模假设与模型构造　a.金属材料为晶体，晶体为面心立方结构，原子间的交互作用限于最近邻；b.空位的形成能定义为：从晶体内部取出一个原子放到晶体

4、表面所需的能量；c.为不显著影响晶体表面的形状，取出的晶体放在晶体表面的台阶。3）模型求解和模型分析面心立方体（配位数为12）取出离子要割断12个键，而在表面台阶处放置一个原子，要形成6个鍵，因此净效应为割断6个鍵，其能量净效应等于晶体的结合能。结论：空位形成能与结合能间有密切的关系。4）模型检验结合能愈大，熔点愈高，则空位形成能愈大。结论基本符合实验事实。实际空位形成能只有结合能的1/2－1/4　　－－未考虑金属鍵的特征和空位周围原子的位移。需重新建模。Ø从内部取出正离子的情况u从晶体点阵中取出一个正离子，设所带正电荷均匀

5、散布于晶体中以抵消它的价电子，使整个晶体仍保持电中性，形成穴位静电效应。u对于单价金属，空位的静电效应，即空位的附加电荷为-Ze(Z=-1)，引起导带电子的屏蔽效应。达到平衡后，空位周围只保留局部的干扰电势(Vp)。(1-1)u设导带的电子浓度为n，静电能的增加等于nVp的体积积分：式中：r表示积分元到空位中心的距离R表示积分区域的半径。(1-2)根据电子屏蔽模型：Vp与Z应满足关系式中：N0(Em)表示晶体导带在费米能级Em处的能态密度。从式(1-1)和(1-2)中消去积分，可得：(1-3)(1-4)对于自由电子：(1-5

6、)这里C是一个常数，将式(1-4)、(1-5)代入(1-3)，则：(1-6)Ø正离子放在表面台阶由于自由电子气的膨胀，造成费米能的下降。由式(1-5)，即：(1-7)式中：N为晶体中的原子数，V为晶体的总体积。若体积膨胀了ΔV，费米能级的变化ΔEm。对式(1-7)两边取微分(1-8)(1-9)对式(1-7)两边取微分，得：由模型可知，ΔV等于一个原子体积，ΔV/V就等于1/N;而自由电子的平均动能为(3/5)EM,，因此总的费米能变化为：因此，在此模型下，即考虑取出原子后的内部能量变化E1和放置在台阶上的能量变化E2后，E1

7、和E2相加就等于空位的形成能：重新建模结果：E=4/15Em空位形成能（实验值与理论值的对照）金属Uf/eV(淬火)Uf/eV（正电子湮灭）Uf/eV（中位值）Uf/eV（计算值）Cu1.271.291.280.8-1.0Ag1.101.161.130.6-0.92Au0.940.970.950.6-0.77u若考虑到空位周围的原子略有松驰，可能降低能量，因而有第三项E3，比如对铜的估计值为-0.3eV。这样，E1、E2、E3三项的迭加就近似等于实验值。金属Uf/eV(淬火)Uf/eV（正电子湮灭）Uf/eV（最佳值）Uf/

8、eV（计算值）Cu1.271.291.280.8-1.0+0.3Ag1.101.161.130.6-0.92Au0.940.970.950.6-0.775）模型应用　　根据模型可求解空位浓度：　式中：c为平衡状态下的空位深度；Uf为空位形成能；Sf为形成一个空位改变了周围原子振动所引起的扰