计算材料学PPT幻灯片课件.ppt

计算材料学Computationalmaterialsscience,张晖材料科学与工程学院Email:huizhang;Mobile1,课程内容,1计算材料学概述2应用Fortran语言实现计算材料学在介观尺度上材料学的一些应用的具体实例。（蒙特卡洛方法、元胞自动机等方法）3MaterialsStudioThermo-Calc;Jmatpro热物性计算，ANSYS有限元模拟等软件的简单介绍和应用。,2,本堂课内容：计算材料学概论本课程研究对象和内容学习要求Fortran语言语法结构简单回顾,3,科学计算各学科（工科和文科）广泛的应用,4,计算模拟及设计在材料研究中的地位,5,计算模拟及设计在材料研究中的地位,6,计算模拟及设计在材料研究中的地位,7,计算模拟与设计在材料研究中的地位,8,计算材料学的概念,计算材料学是近20年里发展起来的一门新兴交叉学科。它综合了凝聚态物理、材料物理学、理论化学、材料力学和工程力学、计算机算法等多个相关学科。本学科的目的是利用现代高速计算机，模拟材料的各种物理化学性质，深入理解材料从微观到宏观多个尺度的各类现象与特征，并对于材料的结构和物性进行理论预言，从而达到设计新材料的目的。本学科目前尚无统一的称呼，计算材料学（Computationalmaterialsscience），计算材料模拟（Computationalmaterialssimulation）、计算材料设计（Computationalmaterialsdesign）“计算分析与模型化”(ComputerAnalysisandModeling)等称谓往往同时使用。,9,计算材料学内容,计算材料学主要包括两个方面的内容：一方面是计算模拟，即从实验数据出发，通过建立数学模型及数值计算，模拟实际过程；另一方面是材料的计算机设计，即直接通过理论模型和计算，预测或设计材料结构与性能。前者使材料研究不是停留在实验结果和定性的讨论上，而是使特定材料体系的实验结果上升为一般的、定量的理论，后者则使材料的研究与开发更具方向性、前瞻性，有助于原始性创新，可以大大提高研究效率。因此，计算材料学是连接材料学理论与实验的桥梁。,10,计算机模拟的概率和步骤,11,计算材料学作用和意义,材料的组成、结构、性能、服役性能是材料研究的四大要素，传统的材料研究以实验室研究为主，是一门实验科学。但是，随着对材料性能的要求不断的提高，材料学研究对象的空间尺度在不断变小，只对微米级的显微结构进行研究不能揭示材料性能的本质，纳米结构、原子像已成为材料研究的内容，对功能材料甚至要研究到电子层次。因此，材料研究越来越依赖于高端的测试技术，研究难度和成本也越来越高。另外，服役性能在材料研究中越来越受到重视，服役性能的研究就是要研究材料与服役环境的相互作用及其对材料性能的影响。随着材料应用环境的日益复杂化，材料服役性能的实验室研究也变得越来越困难。,12,计算材料学作用和意义,总之，仅仅依靠实验室的实验来进行材料研究已难以满足现代新材料研究和发展的要求。然而计算机模拟技术可以根据有关的基本理论，在计算机虚拟环境下从纳观、微观、介观、宏观尺度对材料进行多层次研究，也可以模拟超高温、超高压等极端环境下的材料服役性能，模拟材料在服役条件下的性能演变规律、失效机理，进而实现材料服役性能的改善和材料设计。因此，在现代材料学领域中，计算机“实验”已成为与实验室的实验具有同样重要地位的研究手段，而且随着计算材料学的不断发展，它的作用会越来越大。就实际工程而言，应用数值近似方法进行预测计算，可以有效的减少在优化材料和设计新工艺方面所必须进行的大量试验。材料模拟和材料制备工艺的进步，极大的促进了新产品的优化和开发。,13,计算材料设计的一些成功应用,1.设计新型材料和器件i.高性能磁光记录材料：Tb/Bi/FeCo与Tb/Pb/FeCo超晶ii.超硬材料：C3N4（硬度可以媲美金刚石？）iii.新型锂电池阴极材料：LixCoO2的替代品，Al替代Co？2.预言晶体结构（e.g.，针对70种合金，120晶体进行10000个第一性原理能量计算，六个月）3.计算材料相图4.获得实验难以实现的极端条件下（如高温、高压）的材料结构与物性,14,计算材料学范围和层次,计算材料学主要有两种分类方法：一是按理论模型和方法分类，二是按材料计算的特征空间尺寸(Characteristicspacescale)分类。材料的性能在很大程度上取决于材料的微结构，材料的用途不同，决定其性能的微结构尺度会有很大的差别。例如，对结构材料来说，影响其力学性能的结构尺度在微米以上，而对于电、光、磁等功能材料来说可能要小到纳米，甚至是电子结构。因此，计算材料学的研究对象的特征空间尺度从埃到米。,15,按材料计算的特征空间尺寸分类,根据研究对象的空间和时间尺度不同，材料计算的方法也有很大差别：研究材料的电子结构的方法有基于密度泛函理论、量子力学及固体物理理论的第一原理；介观层次上对体系的模拟近年来有较快的发展，如合金中的相变微观组织演化过程可采用相场动力学或原胞自动化方法，这些方法使人们能够定量地描述不同过程中的组织变化的动力学规律，探索不同因素对微观组织形成的作用；宏观层次上的计算模拟常常采用有限元和有限差分方法，这些方法已经被广泛用语解决材料工程的实际问题，可为实际工艺的设计提供定量化的指导。对于不同的过程其发生的时间尺度也是迥然不同的，相应需要采用不同的模拟方法。对于许多材料的性质，常常由几个层次的结构来决定，因而近年来将不同方法结合起来的多尺度方法受到广泛的重视。,16,计算材料设计的层次和典型模拟方法,17,18,本学期学习重点,常用的计算方法包括：蒙特卡洛方法:是在简单的理论准则基础上，采用反复随即抽样的方法，解决复杂系统的问题。其实质是一种概率和统计的问题。元胞自动机方法:生命游戏，是一种建立在离散的时间和空间上的动力学系统。散布在规则格网（LatticeGrid）中的每一元胞（Cell）取有限的离散状态，遵循同样的作用规则，依据确定的局部规则作同步更新。大量元胞通过简单的相互作用而构成动态系统的演化。,19,第一性原理法,根据原子核和电子相互作用的原理及其基本运动规律，运用量子力学原理，从具体要求出发，经过近似处理后直接求解薛定谔方程的算法，习惯上称为第一性原理法（abinitio）。但由于其计算非常复杂，对计算机要求很高，物理模型是将被研究对象的物理图像理想化和简单化，因为真实的物体运动的行为太复杂，即使是宏观物体，企图得到绝对严格的描述其运动规律也是不可能的，而微观世界的运动更复杂，所以在第一性原理计算中都是合理地忽略影响运动的次要因素，保留主要的因素。第一性原理方法由于直接基于基本的物理原理而不依赖于经验参数，从量子力学基本原理导出的计算方法直接描述了微观粒子的运动规律，所以是精确的，因而具有很强的预测性，在未来合成材料之前先预测其可能的性质，因而对材料的设计具有很强的指导意义，近年其应用得到迅速发展，如金属中合金化效应的预测、金属间化合物中合金原子占据位置的预测、缺陷复合体的电子结构与性质的预测等，但由于其计算中考虑了电子的自由度，其运算量极大，所能研究的体系的尺度很小；,20,本课程的预期目标,开阔眼界，对计算材料科学和材料设计的基本方法和原理有初步的了解。从计算模拟的角度加深对晶体形核和长大的认识。对于将来可能从事实验研究的，在必要时，能够在本人所从事的领域进行一定的计算模拟与设计.对于将来可能从事材料理论和计算方向的，能够在本课程基础上深入学习，顺利进入本领域开展工作。掌握一些对今后科研工作有用的知识,21,主要参考书,1.M.P.AllenandD.J.Tildesley,ComputerSimulationofLiquids,(OxfordUniversityPress,1987).2.J.M.Thijssen,ComputationalPhysics,(CambridgeUniversityPress,1999).3.K.Ohno,K.Esfarjani,K.Kawazoe,ComputationalMaterialsScience,(Springer,1999).4.M.Springborg,MethodofElectronic-StructureCalculations,(Wiley,2000).5.E.Kaxiras,AtomicandElectronicStructureofSolids,(CambridgeUniversityPress,2003).6.R.M.Martin,ElectronicStructure,(CambridgeUniversityPress,2004).7.熊家炯主编，材料设计，（天津大学出版社，2000）8.荷Frenkelread(*,*)MyGuessif(MyGuess=Number)exitif(MyGuessNumber)thenwrite(*,*)Toohigh.Tryagain.elsewrite(*,*)Toolow.Tryagain.endifenddowrite(*,*)Youarelucky.Itis,Numberend,36,数组,输入10个整数，并按输入时的逆序输出，每行5个。,PROGRAMNUMBERINTEGERDIMENSION(1:10)AINTEGERIDOI=1,10READ(*,*)A(I)ENDDOWRITE(*,*)(A(11-I),I=1,10)ENDPROGRAMNUMBER,37,1.数组名INTEGERA(1:2,1:3)READ(*,*),A相当于READ*,A(1,1),A(2,1),A(1,2),A(2,2),A(1,3),A(2,3)2.隐含DO循环WRITE(*,*)(A(I,J),J=1,3),I=1,2)READ(*,*)(A(I,J)，I=1,2)，J=1,3)相当于READ(*,*),A3.DO结构加隐含DO循环最常用方法,DOI=1,2READ(*,*),(A(I,J),J=1,3)ENDDO,DOI=1,2WRITE(*,*)(A(I,J),J=1,3)ENDDO,38,根据所操作的文件是否在内存中可以把文件分为内部文件和外部文件。当把内存中的数据记录到到磁盘的文件中或输入输出到其他外部设备如打印机、显示器、键盘上时，被称为外部文件。连接着一个外部文件的设备描述必须是整型表达式或是星号(*)，其整型表达式的取值范围在-32768到32767之间。例：OPEN(UNIT=10,FILE=f:IF1.DAT)WRITE(10，(A)howareyou?例：OPEN(7,FILE=E:HYZ.DAT)WRITE(7,*)ACLOSE(7)在Fortran中有4个预定义的外部文件(设备)：,外部文件,39,利用随机过程产生10个处于区间10,99上的随机整数。,PROGRAMRANDOM_ARRAYINTEGERIINTEGERN_RAN(1:10)REALX_RANDOI=1,10CALLRANDOM_NUMBER(X_RAN)N_RAN(I)=INT(X_RAN*90+10)ENDDOWRITE(*,*)(N_RAN(I),I=1,10)ENDPROGRAMRANDOM_ARRAY,40,属性说明,PARAMETER属性PARAMETER属性也称常数名属性。被说明对象一旦附加了PARAMETER属性，就不再是变量名，而是一个常数的名字，它的形式虽与变量名形式一样，但在程序中不能改变值，只能当常数使用。例如：INTEGER，PARAMETER:K=20REAL(8),PARAMETER:PI=3.141592654,K_PAI=K*PI,Light_Speed=2.99654E10常数名被说明后，在程序中就不可再更改值，如要更改，则系统提示出错信息。这一性质，有助于防止程序中误改不可改的常数。使用常数名可以加强可读性，也有助于程序维护。设原编程序是解四阶方程以后要扩充为能解五阶方程，修改时必须把所有反映阶数，写在程序各处的4改成5，这就难免修改时有逃漏，甚至把不该改的4也误改为5，造成程序混乱。使用常数名就可避