第二章_407606549.pdfVIP

第二章计算材料学基础 2.1 材料模拟的基本思路 • 定义一系列自变量和因变量 • 建立数学模型 • 解数学模型 自变量和因变量 • 自变量：一般是指对空间和时间的量化 • 因变量：一般是指态变量（如温度，位移） 这些变量的选择要基于满足所研究材料性质 的计算精度要求 数学模型 由两部分组成 • 状态方程 – 用于描述由给定态变量定义的材料性质（如胡克定律， 相互作用原子势函数）， – 状态方程与路径无关 • 演化方程： – 用于描述态变量作为自变量的函数的变化情况(如牛顿 运动方程） – 当微结构处于热力学非平衡态时，决定演化方程的因 素与路径有关。 这些方程通常都是微分方程，通常要把这些方程 变成差分方程的形式，并确定求解问题的初始条 件和边界条件 数学模型 • 量子力学模型和经典模型 – 材料的物理性能：常常需要量子力学模型解 – 材料的微观结构问题：常常可以用经典模型解 • 离散模型和连续模型 – 电子和晶格缺陷问题：离散模型（atomistic simulation) – 介观/宏观层次模拟：连续模型 解数学模型 • 密度泛函方法/能带计算/分子轨道理论（量子） • 分子动力学/蒙特卡洛方法 （经典/量子，离散） • 有限元法（经典，连续） 计算材料学对数学模型公式化的基本步骤 1. 定义自变量 – 例如时间和空间 2. 定义因变量 – 例如温度，浓度、位移等 3. 建立运动学方程 – 亦即在不考虑实际作用力时，确定描述质点坐标变化 的函数关系。例如，在一定约束条件下，建立根据位 移梯度计算应变的方程。 4. 建立状态方程 – 亦即从因变量的取值出发，确定描述材料实际状态且 与路径无关的函数 计算材料学对数学模型公式化的基本步骤 5. 建立演化方程 – 亦即从因变量的变化，给出描微结构演化的 且与路径有关的函数 6. 相关物理参数的确定 7. 边界条件和初始条件 8. 确定用于求解由步骤1－7建立的联立方程 组的数值算法或解析算法 2.2计算材料学的基本原理 (参考罗伯书，第四章) r r  r,t  对于一个质量为m，坐标为 且在势场U( )中运动的粒子，若令其波函数为  ，则它的运动可由含时间的薛定谔方程表示 2  r,t  2     r,t U(r ) r,t i 2m     t (2-1-1) 式中  h 2 2 2 2 ，    。 2 x2 y2 z2 r 如果势场U( ) 与时间无关，令  E t  (2-1-2)  r,t  (r )exp i k   k   r