材料加工过程的数值模拟之温度场数值模拟.pptVIP

2-5不稳定导热的有限元解法二、有限元发的解题思想和步骤 建立温度的差值函数 对于三角形单元：T=f(Ti,Tj,Tm) T=NiTi+NjTj+NmTm 单元变分计算 2-5不稳定导热的有限元解法二、有限元发的解题思想和步骤 总体合成 得到线性方程组。 求解线性方程组 2-5不稳定导热的有限元解法三、内单元计算格式的建立 1、一维系统(略去课堂不讲) 1）模型： 2）泛函： 3)温度差值函数 2-5不稳定导热的有限元解法二、内单元计算格式的建立 4）单元变分计算 三、建立内节点差分方程/一维系统 1、模型： ??0，0?x?L 2、初始条件：T(x,0)=?(x) 3、边界条件：T(0, ?)=?1(?)， T(L, ?)=?2(?) 4、区域离散 距离步长：?x=xi-xi-1, xi =(i-1) ?x 时间步长： ? ?= ?n- ?n-1, ?n=n ? ? Tin=T(xi, ?n) 三、建立内节点差分方程/一维系统 5、有限差分方程建立 1）显示差分 点（i,n）的导热方程为： 三、建立内节点差分方程 三、建立内节点差分方程/一维系统 2)隐式差分格式 温度对距离的二阶偏微商是对应时刻n+1的，而温度对时间的一阶偏微商是对应时刻n的。差分方程为： 截断误差：O[? ?+（ ? x）2]，