材料计算与模拟2讲解.pptVIP

第二章 材料计算示例 ? 2.1 方程的求解与简单计算问题 ? 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 = fsolve(fun,x0,option); fplot(x^3+10*x^2-2*sin(x)-50,[-15 5]),grid on -15 -10 -5 0 5 -1200 -1000 -800 -600 -400 -200 0 200 400 f=inline(x^3+10*x^2-2*sin(x)- 50,x); fplot(f,[-15 5]); grid on fsolve(f,-10) fsolve(f,-1) fsolve(f,1) 或： fsolve(x^3+10*x^2-2*sin(x)-50,1) fsolve(x^3+10*x^2-2*sin(x)-50,-1) fsolve(x^3+10*x^2-2*sin(x)-50,-10) 或： fsolve(@ff,1) fsolve(@ff,-1) fsolve(@ff,-10) 其中 ff 为函数文件，格式为： function f=ff(x) f=x^3+10*x^2-2*sin(x)-50; 例 2.1.2 热缺陷浓度的计算 热缺陷是由于热起伏引起的，并与温度有关。故在某一温度下， 热缺陷的数目可以用热力学中自由能最小原理来进行计算。热缺 陷浓度与温度的关系式为 式中： N 表示完整的单质晶体的原子数目； n 表示热振动形成的空 位数； n/N 表示热缺陷在总结点中所占分数，既热缺陷浓度。 ΔG f 分别代表空位形成自由能或填充缺陷形成自由能。计算 ΔG f 分别等 于 1eV, 2eV, 4eV, 6eV 和 8eV 时，温度在 100 、 1200 和 2000 ℃下的 缺陷浓度。 Boltzmann 常数 k = 1.38 × 10 -23 J/K 。 1ev = 1.59 × 10 - 19 J 。 ? ? ? ? ? ? ? ? ? ? ? kT G N n f 2 exp 例 2.1.2 热缺陷浓度的计算 例 2.1.2 热缺陷浓度的计算 例 2.1.3 自由能温度关系式 CaCO 3 分解反应的自由能 ΔG 与温度的关系为 : ΔG=186.08+10.7 × 10 -3 TlnT+4.187 × 10 -6 T 2 - 5.23 × 10 2 T -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 = -y 2