铝酸钙力学性质的理论计算于艳超资料.ppt

Page ? * LOGO LOGO 指导教师：李美玲 答辩人：于艳超 2013年6月3日 目录 引言 密度泛函理论 结论 铝酸钙CaAl2O4力学性质的理论计算 引言 2密度泛函理论 2.2 CASTEP简介 CASTEP（Cambridge Sequential Total Energy Package）为一商业软件包，它始于剑桥大学凝聚态理论研究组开发的一系列程序。CASTEP可以进行总能量赝势用于解量子力学方程，以得到包含原子随意排列系统的电子状态，这些计算给出系统的基态能量和电子密度，允许计算与总能量相关的任何物理量（如晶格常数、弹性常数、几何结构等）。CASTEP足以应付在材料科学中发现的各种分子、固态及表面问题。现已在广泛的材料科学领域中得到应用。CASTEP软件包提供了界面友好的从头计算方法，是目前计算机模拟实验中最先进、最重要的方法之一。 3 铝酸钙CaAl2O4力学性质的理论计算 3.1晶体结构 本文应用密度泛函理论借助castep软件计算了能量最低的铝酸钙晶体结构，得到CaAl2O4晶格参数数据如表3-1。 实验值 理论值 Vo （?3） 144.62（8） 141.925 a （?） 5.556（2） 5.5118 b （?） 3.7627（8） 3.7352 c （?） 7.056（3） 7.0318 β 101.36（3） 101.377 密度/g/cm3 3.629 3.698 表3-1 计算的CaAl2O4晶格参数 从表中数据可看出，理论值与实验值十分接近。说明理论与实验吻合很好。由表3-1可看出计算得到的CaAl2O4晶格体积略小于实验值，这主要是由于计算材料不是单相引起的。 3.2 力学性质 由于铝酸钙为单斜晶系，独立的弹性常数有13个 分别为C11、C22、C33、C44、C55、C66、C12、C13、C15、C23、C25、C35、C46。表3-2给出了计算得到的具体数据。 3-2 铝酸钙的弹性常数表 Cij (Gpa） C11 397.58 C22 247.32 C33 290.63 C44 78.15 C55 69.60 C66 136.96 C12 94.58 C13 48.45 续表 C15 -48.51 C23 103.06 C25 0.25 C35 31.26 C46 -11.53 在多晶或宏观材料弹性常数的理论估算方面，Voigt、Reuss、Hill和Kroner等进行了大量的研究。在Voigt和Reuss模型中, 分别假设材料中的各晶粒处于等应变和等应力状态。显然,这样的假设和实际情况有很大的差别, Hill用极值原理证明,Voigt和Reuss模型的计算结果是弹性常数的上下限。Hill模型则将Voigt和Reuss模型的计算结果取一个简单的算术平均, 其结果和实际测定值较为接近。 在Voigt-Reuss-Hill近似下， 再根据单斜晶系力学稳定性判据：  由表3-2中的数据，计算得到以下数据：  满足力学稳定性判据，因而铝酸钙晶体材料力学性质是稳定的。 结论 本文采用基于密度泛函理论的第一性原理平面波赝势方法对铝酸钙结构进行了几何优化。此过程得到了能量最低的晶体结构，此时结构最稳定。在此基础上，对该材料的力学性质进行了计算，对取得的结果进行了分析，得到以下结论： （1）对铝酸钙的晶体结构进行了优化，得到了优化的晶体结构。晶格参数为a=5.5118 ?， b=3.7352 ?, C=7.0318 ? , =101.377，晶格体积V=141.925 ?3。与实验值（a=5.556 ?， b=3.7627 ?, C=7.056 ?, =101.36晶格体积V=144.62 ?3）相比略小。这主要是由于计算材料不是单相引起的。说明我们的计算合理，结果准确可靠。 （2）计算得到的铝酸钙晶体的体模量为153.64Gpa、杨模量162.38Gpa、切模量为66.28Gpa、泊松比为 =0.31。根据力学稳定性判据得知，单斜晶系的铝酸钙材料力学性质是稳定的。 由于目前对铝酸钙晶体的力学性质的研究比较少，从文献和资料中没有查到关于理想铝酸钙晶体体模量、杨模量、切模量等实验数值，因此本文所得数据可为铝酸钙晶体的晶体学及力学性能的研究做一些有价值的参考。 本文内容比较肤浅，请各位老师批评指正，谢谢！ * * * Page ? * LOGO LOGO * * *