计算材料学PPT幻灯片课件

《计算材料学PPT幻灯片课件》由会员分享，可在线阅读，更多相关《计算材料学PPT幻灯片课件（42页珍藏版）》请在金锄头文库上搜索。

1、计算材料学Computationalmaterialsscience 张晖材料科学与工程学院Email huizhang Mobile 13865606861 1 课程内容 1计算材料学概述2应用Fortran语言实现计算材料学在介观尺度上材料学的一些应用的具体实例 蒙特卡洛方法 元胞自动机等方法 3MaterialsStudioThermo Calc Jmatpro热物性计算 ANSYS有限元模拟等软件的简单介绍和应用 2 本堂课内容 计算材料学概论本课程研究对象和内容学习要求Fortran语言语法结构简单回顾 3 科学计算 各学科 工科和文科 广泛的应用 4 计算模拟及设计在材料研究中的地

2、位 5 计算模拟及设计在材料研究中的地位 6 计算模拟及设计在材料研究中的地位 7 计算模拟与设计在材料研究中的地位 8 计算材料学的概念 计算材料学是近20年里发展起来的一门新兴交叉学科 它综合了凝聚态物理 材料物理学 理论化学 材料力学和工程力学 计算机算法等多个相关学科 本学科的目的是利用现代高速计算机 模拟材料的各种物理化学性质 深入理解材料从微观到宏观多个尺度的各类现象与特征 并对于材料的结构和物性进行理论预言 从而达到设计新材料的目的 本学科目前尚无统一的称呼 计算材料学 Computationalmaterialsscience 计算材料模拟 Computationalmater

3、ialssimulation 计算材料设计 Computationalmaterialsdesign 计算分析与模型化 ComputerAnalysisandModeling 等称谓往往同时使用 9 计算材料学内容 计算材料学主要包括两个方面的内容 一方面是计算模拟 即从实验数据出发 通过建立数学模型及数值计算 模拟实际过程 另一方面是材料的计算机设计 即直接通过理论模型和计算 预测或设计材料结构与性能 前者使材料研究不是停留在实验结果和定性的讨论上 而是使特定材料体系的实验结果上升为一般的 定量的理论 后者则使材料的研究与开发更具方向性 前瞻性 有助于原始性创新 可以大大提高研究效率 因此

4、计算材料学是连接材料学理论与实验的桥梁 10 计算机模拟的概率和步骤 11 计算材料学作用和意义 材料的组成 结构 性能 服役性能是材料研究的四大要素 传统的材料研究以实验室研究为主 是一门实验科学 但是 随着对材料性能的要求不断的提高 材料学研究对象的空间尺度在不断变小 只对微米级的显微结构进行研究不能揭示材料性能的本质 纳米结构 原子像已成为材料研究的内容 对功能材料甚至要研究到电子层次 因此 材料研究越来越依赖于高端的测试技术 研究难度和成本也越来越高 另外 服役性能在材料研究中越来越受到重视 服役性能的研究就是要研究材料与服役环境的相互作用及其对材料性能的影响 随着材料应用环境的日益复

5、杂化 材料服役性能的实验室研究也变得越来越困难 12 计算材料学作用和意义 总之 仅仅依靠实验室的实验来进行材料研究已难以满足现代新材料研究和发展的要求 然而计算机模拟技术可以根据有关的基本理论 在计算机虚拟环境下从纳观 微观 介观 宏观尺度对材料进行多层次研究 也可以模拟超高温 超高压等极端环境下的材料服役性能 模拟材料在服役条件下的性能演变规律 失效机理 进而实现材料服役性能的改善和材料设计 因此 在现代材料学领域中 计算机 实验 已成为与实验室的实验具有同样重要地位的研究手段 而且随着计算材料学的不断发展 它的作用会越来越大 就实际工程而言 应用数值近似方法进行预测计算 可以有效的减少在

6、优化材料和设计新工艺方面所必须进行的大量试验 材料模拟和材料制备工艺的进步 极大的促进了新产品的优化和开发 13 计算材料设计的一些成功应用 1 设计新型材料和器件i 高性能磁光记录材料 Tb Bi FeCo与Tb Pb FeCo超晶ii 超硬材料 C3N4 硬度可以媲美金刚石 iii 新型锂电池阴极材料 LixCoO2的替代品 Al替代Co 2 预言晶体结构 e g 针对70种合金 120晶体进行10000个第一性原理能量计算 六个月 3 计算材料相图4 获得实验难以实现的极端条件下 如高温 高压 的材料结构与物性 14 计算材料学范围和层次 计算材料学主要有两种分类方法 一是按理论模型和方

7、法分类 二是按材料计算的特征空间尺寸 Characteristicspacescale 分类 材料的性能在很大程度上取决于材料的微结构 材料的用途不同 决定其性能的微结构尺度会有很大的差别 例如 对结构材料来说 影响其力学性能的结构尺度在微米以上 而对于电 光 磁等功能材料来说可能要小到纳米 甚至是电子结构 因此 计算材料学的研究对象的特征空间尺度从埃到米 15 按材料计算的特征空间尺寸分类 根据研究对象的空间和时间尺度不同 材料计算的方法也有很大差别 研究材料的电子结构的方法有基于密度泛函理论 量子力学及固体物理理论的第一原理 介观层次上对体系的模拟近年来有较快的发展 如合金中的相变微观组织

8、演化过程可采用相场动力学或原胞自动化方法 这些方法使人们能够定量地描述不同过程中的组织变化的动力学规律 探索不同因素对微观组织形成的作用 宏观层次上的计算模拟常常采用有限元和有限差分方法 这些方法已经被广泛用语解决材料工程的实际问题 可为实际工艺的设计提供定量化的指导 对于不同的过程其发生的时间尺度也是迥然不同的 相应需要采用不同的模拟方法 对于许多材料的性质 常常由几个层次的结构来决定 因而近年来将不同方法结合起来的多尺度方法受到广泛的重视 16 计算材料设计的层次和典型模拟方法 17 18 本学期学习重点 常用的计算方法包括 蒙特卡洛方法 是在简单的理论准则基础上 采用反复随即抽样的方法

9、解决复杂系统的问题 其实质是一种概率和统计的问题 元胞自动机方法 生命游戏 是一种建立在离散的时间和空间上的动力学系统 散布在规则格网 LatticeGrid 中的每一元胞 Cell 取有限的离散状态 遵循同样的作用规则 依据确定的局部规则作同步更新 大量元胞通过简单的相互作用而构成动态系统的演化 19 第一性原理法 根据原子核和电子相互作用的原理及其基本运动规律 运用量子力学原理 从具体要求出发 经过近似处理后直接求解薛定谔方程的算法 习惯上称为第一性原理法 abinitio 但由于其计算非常复杂 对计算机要求很高 物理模型是将被研究对象的物理图像理想化和简单化 因为真实的物体运动的行为太复

10、杂 即使是宏观物体 企图得到绝对严格的描述其运动规律也是不可能的 而微观世界的运动更复杂 所以在第一性原理计算中都是合理地忽略影响运动的次要因素 保留主要的因素 第一性原理方法由于直接基于基本的物理原理而不依赖于经验参数 从量子力学基本原理导出的计算方法直接描述了微观粒子的运动规律 所以是精确的 因而具有很强的预测性 在未来合成材料之前先预测其可能的性质 因而对材料的设计具有很强的指导意义 近年其应用得到迅速发展 如金属中合金化效应的预测 金属间化合物中合金原子占据位置的预测 缺陷复合体的电子结构与性质的预测等 但由于其计算中考虑了电子的自由度 其运算量极大 所能研究的体系的尺度很小 20 本

11、课程的预期目标 开阔眼界 对计算材料科学和材料设计的基本方法和原理有初步的了解 从计算模拟的角度加深对晶体形核和长大的认识 对于将来可能从事实验研究的 在必要时 能够在本人所从事的领域进行一定的计算模拟与设计 对于将来可能从事材料理论和计算方向的 能够在本课程基础上深入学习 顺利进入本领域开展工作 掌握一些对今后科研工作有用的知识 21 主要参考书 1 M P AllenandD J Tildesley ComputerSimulationofLiquids OxfordUniversityPress 1987 2 J M Thijssen ComputationalPhysics Cambr

12、idgeUniversityPress 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 CambridgeUniversityPr

13、ess 2004 7 熊家炯主编 材料设计 天津大学出版社 2000 8 荷 Frenkel Smit著 分子模拟 从算法到应用 化学工业出版社 2002 9 德 D 罗伯编著 计算材料学 化学工业出版社 2002 10 吴兴惠 项金种编著 现代材料计算与设计教程 电子工业出版社 2002 11 陈舜麟编著 计算材料科学 化学工业出版社 2005 12 马文淦编著 计算物理学 科学出版社 2005 22 Fortran语言语法结构简单回顾 1Fortran的历史第一代Fortran语言是在1954年提出来的 称为FortranI 23 1 2 1编程实例 C 求两种平均値 PROGRAMExa

14、mple 1 1REALa b av1 av2READ a bav1 a b 2av2 sqrt a b WRITE av1 av2END程序中第1行是注释行 对程序起说明作用 F77注释行是以 C 或 作为该行第一个字符的 F90可在任意一行末以 开始作为注释符 第2行是主程序名 第3行是变量类型定义 第4行是输入语句 第5 6行是赋值部分 第7行是打印输出语句 最后是程序结束 24 常量和变量类型 a 常量常量是指其值始终不变的一些量 整型 实型 双精度和复型常量是算数型常量 也为常数 整型 实型 复型 逻辑型 字符型 数组 b 变量变量是指在程序运行期间其值是可以变化的量 系统为程序中的

15、每一个变量开辟一个存储单元 用来存放变量的值 常量是分为类型的 而变量是用来存放常量的 因此变量也相应地区分为整型变量INTEGER 实型变量REAL 双精度变量DOUBLEPRECISION 复型变量COMPLEX 逻辑型变量LOGICAL 字符型变量CHARACTER 在程序中应当说明哪些变量是整型变量 哪些变量是实型变量 变量在内存中所占的字节数和数据存储形式与相应类型的常数相同 例如 实型变量一般占4个字节 按指数形式存放 在程序中规定变量的类型可以用以下几种方法 隐含约定 Fortran规定 凡以字母I J K L M N六个字母开头的变量名 如无另外说明则为整型变量 以其它字母开头

16、的变量为实型变量 可以将这个隐含约定称为 I N规则 表示用I到N之间的字母开头的变量为整型 例如 I J IMAX NUMBER LINE JOB Kl为整型变量 而A Bl COUNT AMOUNT TOTAL BOOK为实型变量 25 内在函数 26 1 4 4赋值语句f 执行时的类型转换一个算术赋值语句中的被赋值的变量 V 和表达式 e 的类型可以相同 也可以不相同 F90作如下规定 如果变量V与表达式e的类型相同 则直接进行赋值 如果类型不同 则应先进行表达式的求值 然后将该表达式的值转换为被赋值变量的类型 如 I 3 5 2 1表达式3 5 2 1的值为7 35 实型 而变量I为整型 因此I的值等于7 又如 T 3 5 7表达式的值为2 整型 由于T为实型 故系统先将整数2转换成实数2 0 再赋给变量T T的值为2 0 当赋值号两侧的类型不同时 往往会产生程序设计者事先预想不到的结果 所以在编写程序时 应尽可能使两侧保持同类型 27 关系表达式 a 关系运算符关系运算符就是关系比较符 Fortran中有六个关系运算符 28 选择结构 2 3 1块IF构造a 块IF的组成和执