材料计算与模拟2.ppt

第二章 材料计算示例 2.1 方程的求解与简单计算问题 2.2 线性回归问题 例2.1.2 热缺陷浓度的计算 热缺陷是由于热起伏引起的，并与温度有关。故在某一温度下，热缺陷的数目可以用热力学中自由能最小原理来进行计算。热缺陷浓度与温度的关系式为 式中：N表示完整的单质晶体的原子数目；n表示热振动形成的空位数；n/N表示热缺陷在总结点中所占分数，既热缺陷浓度。ΔGf分别代表空位形成自由能或填充缺陷形成自由能。计算ΔGf分别等于1eV, 2eV, 4eV, 6eV和8eV时，温度在100、1200和2000℃下的缺陷浓度。Boltzmann常数k = 1.38×10-23 J/K。1ev = 1.59×10-19 J。 例2.1.2 热缺陷浓度的计算 例2.1.2 热缺陷浓度的计算 例 2.1.3 自由能温度关系式 CaCO3 分解反应的自由能ΔG与温度的关系为: ΔG=186.08+10.7×10-3TlnT+4.187×10-6T2-5.23×102T-1-0.245T (1) 计算温度区间800～1400K范围内的ΔG。(间隔100K) (2) 绘出ΔG与温度的关系图。 (3) 计算ΔG=0的温度条件。 ? 例2.1.3 自由能温度关系式 例2.1.4 微分方程的数值解 常用函数有：ode23 – 精度10-3 , ode45 – 精度 10-6 格式： [T,Y] = solver(odefun,tspan,y0) [T,Y] = solver(odefun,tspan,y0,options) [T,Y] = solver(odefun,tspan,y0,options,p1,p2...) [T,Y,TE,YE,IE] = solver(odefun,tspan,y0,options) sol = solver(odefun,[t0 tf],y0...) where solver is one of ode45, ode23, ode113 odefun方程: 例如y’=f(t,y)或M(t,y)y’ = f(t,y) tspan 变量的间隔 y0 初始条件 例2.1.4：解方程y’ = -y2+x, y(0)=1 函数文件：(函数名：fode1) function dy = fode1(x,y) dy=-y^2+x 主程序： [x,y] = ode23(@fode1,[0,10],[1]); plot(x,y) xlabel(x); ylabel(y); 例2.1.5 解方程： 例2.1.5：Matlab程序 定义函数： function dy = fode2(x,y) dy = zeros(2,1); % a column vector dy(1) = y(2); dy(2) = 1000*(1 - y(1)^2)*y(2) - y(1); 主程序: [x,y] = ode23(@fode2,[0,5],[2,0]); plot(x,y(:,1)) xlabel(x); ylabel(y); 2.2 线性回归问题 例2.2.1 实验测量含有不同量α-氯代萘的一组线性聚乙烯试样的熔点，得到的数据如下： α- 氯代萘的体积分数φ1 0.00 0.06 0.16 0.32 0.52 0.75 0.95 聚乙烯的熔点Tm(K) 410.65 407.65 404.15 393.15 393.15 388.15 383.15 假定熔点Tm与低分子稀释剂的体积分数φ1的关系为 式中，ρu和ρ1分别是高分子和低分子稀释剂的密度。R气体常数。χ1是高分子和稀释剂的相互作用参数。ΔHu是每摩尔重复单元的熔融热。求回归系数。 2.2 线性回归问题 2.2 线性回归问题 2.2 线性回归问题 2.2 线性回归问题 * * 例2.1.1 求方程的解 y=x^3+10*x^2-2*sin(x)-50 分析：用matlab求解方程的符号解，所用函数为SOLVE(eqn1,eqn2,...,eqnN) SOLVE(eqn1,eqn2,...,eqnN,var1,var2,...,varN) SOLVE(eqn1,eqn2,...,eqnN,var1,var2,...varN) 求解非线性方程的数值解(最小二乘法)所用函数fsolve 函数文件格式(myfun为函数文件名) x = fsolve(@myfun,x0,option) 内联函数格式(fun为内联函数名) x = fs