(计算机在材料科学与工程中的应用)第四章材料科学与工过程中典型物理场的数值模拟.ppt

7）整理数据 单击Mesh-Export Mesh…输出 p e t 的数值，单击Solve-Export Solution…输出u. 回到Matlab主窗口执行下面两条命令： u1=[p‘,u(:,14401)] ！将节点坐标和其在14400s（即4h）时的碳浓度组成新矩阵 u2=sortrows(u1,2) ！将u1按y值大小排列。 8）求解 观察u2可发现，y=0.0004时，C=0.0016,并可看出渗碳厚度≤7.09×10-4m。 4.2 浓度场的数学模型及求解 用PDE工具箱求解浓度分布 对一10cm×10cm×5cm的碳含量为wc＝0.8%长方体钢材进行渗硼，渗硼温度为850℃。已知850℃时硼在T8中的扩散系数D=1.10×10-13 m2/s，分析渗层厚度随时间（0～9000s）变化的情况。 4.2 浓度场的数学模型及求解 用PDE工具箱求解浓度分布 2 4.2 浓度场的数学模型及求解 用PDE求解扩散方程 3 to solve the diffusion equation using the Forward-Time /Centered-Space (FTCS) scheme. Ref: Alejandro Garcia, Numerical Methods for Physics, 2e, 2002. /matlabcentral/fileexchange/2270 Run:dftcs.m Enter time step: 0.0001 Enter the number of grid points: 51 Solution is expected to be stable 4.2 浓度场的数学模型及求解 用PDE求解扩散方程 3 Enter time step: 0.0001 Enter the number of grid points: 51 Solution is expected to be stable 4.2 浓度场的数学模型及求解 用PDE求解扩散方程 3 to solve the neutron diffusion equation (中子扩散方程 )using the Forward-Time /Centered-Space (FTCS) scheme. Ref: Alejandro Garcia, Numerical Methods for Physics, 2e, 2002. /matlabcentral/fileexchange/2270 Run: neutrn.m Enter time step: 0.0005 Enter the number of grid points: 61 Enter system length: 2 = System length is subcritical Solution is expected to be stable Enter number of time steps: 12000 4.2 浓度场的数学模型及求解 用PDE求解扩散方程 3 Enter time step: 0.0005;Enter the number of grid points: 61 Enter system length: 2 = System length is subcritical Solution is expected to be stable Enter number of time steps: 12000 4.2 浓度场的数学模型及求解 * * * * 薄板焊接中移动热源为例,取焊件的一半作为模型进行离散化(图4-3).电弧起始点为O点,此后以速度v沿y轴移动,经过τ时间后到这O’点，此时由于热源引起的热能分布为： * 薄板焊接中移动热源为例,取焊件的一半作为模型进行离散化(图4-3).电弧起始点为O点,此后以速度v沿y轴移动,经过τ时间后到这O’点，此时由于热源引起的热能分布为： * 薄板焊接中移动热源为例,取焊件的一半作为模型进行离散化(图4-3).电弧起始点为O点,此后以速度v沿y轴移动,经过τ时间后到这O’点，此时由于热源引起的热能分布为： * 薄板焊接中移动热源为例,取焊件的一半作为模型进行离散化(图4-3).电弧起始点为O点,此后以速度v沿y轴移动,经过τ时间后到这O’点，此时由于热源引起的热能分布为： * 薄板焊接中移动热源为例,取焊件的一半作为模型进行离散化(图4-3).电弧起始点为O点,此后以速度v沿y轴移动,经过τ时间后到这O’点，此时由于热源引起的热能分布为： * 薄板焊接中移动热源为例,取焊件的一半作为模型进行离散化(图4-3).电弧