ONI/ aIN aN ON ]{T}a fo IN] OdQ2-CIN]q,dr (6-20) +,MNM]Td-后Nhd=0 公式(6-20)是n个联立的线性方程组,可以确定n个结点的温度T。按有限元格式将 (6-20)表示为, [K]{T}°={P} (6-21) 其中矩阵K]为单元的导热矩阵或称为温度刚度矩阵,{T}°为单元的结点温度向量,{P}° 称为单元的温度载荷向量或热载荷向量( Thermal load vector)。对于某个特定单元,单元导 热矩阵K]和温度载荷向量{P}°的元素分别为 N 12+ hN, N (6-22) P=Nq,dr+ Nht d+ n,Odr (6-23) 如果某个单元完全处于物体的内部 K,=l L)dQ ax a P=LNOdT 在整个物体上的加权积分方程是单元积分方程的和 aN aINIaN/ O[N ∑Non2-∑「M ∑MNN;dr-∑ C INThT, dr=0 根据单元结点的局部编号与整体编号的关系,直接求和得到整体刚度矩阵,整体方程组 为, [K]{T}={P}[ ] [ ]{ } [ ] 0 [ ] [ ] )]{ } [ ] ) ( [ ] ) ( [ ] ) ( [ ] [( 3 3 2 +  −  = −  −       +               h N N T d N hT d N Qd N q d T d y N y N x N x N f T e T e e e s T T e e y T x T e   （6-20） 公式（6-20）是 n 个联立的线性方程组，可以确定 n 个结点的温度 Ti 。按有限元格式将 （6-20）表示为， e e e [K] {T} = {P} （6-21） 其中矩阵[K]e 为单元的导热矩阵或称为温度刚度矩阵，{T}e 为单元的结点温度向量，{P}e 称为单元的温度载荷向量或热载荷向量（Thermal load vector）。对于某个特定单元，单元导 热矩阵[K]e和温度载荷向量{P}e的元素分别为，  +      +     =   d hN N d y N y N x N x N K i j e j i y j i x e ij 3 (  ) （6-22） =  +  +     P N q d N hT d NiQd e i f e i s e i 2 3 （6-23） 如果某个单元完全处于物体的内部，      +     =  d y N y N x N x N K j i y j i x e ij (  ) =   P NiQd e i 在整个物体上的加权积分方程是单元积分方程的和， [ ] [ ]{ } [ ] 0 [ ] [ ] )]{ } [ ] ) ( [ ] ) ( [ ] ) ( [ ] [( 3 3 2  −  = −  −  +      +                    h N N T d N hT d N Qd N q d T d y N y N x N x N f T e e T e e e e s T e T e e e y T x T e e   （6-24） 根据单元结点的局部编号与整体编号的关系，直接求和得到整体刚度矩阵，整体方程组 为， [K]{T} = {P}