x(N]2-)-x ON OT 假定单元内的导热系数为常数, N b 4故=4bb园a4 (6-27) b: b, b, b,bm IT. b,b, b, b,bm rT b, bm ba ITm Nx(4)=N) dT、2NyaT ∫,aNr 4421°c,c,cmT Cj 4A C, Cm C, m Cm Tm 单元的刚度矩阵为, 6: b6 b 6 C C:c [K]=4 bb, b2 c, c, c c,cm bb bb b2 CCCc 显然,单元的导热矩阵是对称的 如果单元的内部热源为常数,由内部热源产生的温度载荷项为, IN] OdA=ON=CAI (6-29) 由 Green公式可得 (N]42)+-(N,]dA ax a (6-30) f(NTA, an,+[N]'a, aTm ds 方便起见,把换热边界统一表示为第三类换热边界,x T x N x x T N x x T x N T x T x T     −     =     [ ] [ ] ( ) ([ ]  )  假定单元内的导热系数为常数，                     = −                     = −     −   m j i i m j m m i j j j m i i j i m x m j i i j m m j i A x T x A T T T b b b b b b b b b b b b b b b A dA T T T b b b b b b A dA x T x N 2 2 2 2 4 [ ] 4 [ ] 1    （6-27） y T y N y T N y y T y N T y y T y T     −     =     [ ] [ ] ( ) ([ ]  )                      = −                     = −     −   m j i i m j m m i j j j m i i j i m x m j i i j m m j i A x T y A T T T c c c c c c c c c c c c c c c A dA T T T c c c c c c A dA y T y N 2 2 2 2 4 [ ] 4 [ ] 1    （6-28） 单元的刚度矩阵为，           +           = 2 2 2 2 2 2 4 4 [ ] i m j m m i j j j m i i j i m y i m m m i j j j m i i j i j e x c c c c c c c c c c c c c c c b b b b b b b b b b b b b b b K j   显然，单元的导热矩阵是对称的。 如果单元的内部热源为常数，由内部热源产生的温度载荷项为，           =           =   1 1 1 3 [ ] QA dA N N N N QdA Q m j i A T A （6-29） 由 Green 公式可得 n dS y T n N x T N dA y T N x y T N x y y T x x T s y T x T A ([ ] [ ] ) [ ([ ] ) ([ ] ]   +   =     +           （6-30） 方便起见，把换热边界统一表示为第三类换热边界