单元刚度矩阵为 kr=2』[B[ DIB]rrde (4-28) 单元刚度矩阵的分块矩阵为, KK]=2[ [B] [D][B],rdrd-s (4-29) 由于几何矩阵中的元素不是常量,单元刚度矩阵需要通过积分得到,为简化计算,可以 用三角形单元形心位置的坐标r,=c代替[B矩阵中的变量r、z r +r+r 应变矩阵变成 B]=园 b b;+c1;c0 B (4-30) b 单元刚度矩阵的近似表达式为 [K=2mB]ID B 单元刚度矩阵的分块矩阵近似表达式为, 区]=2mBD[B [K] E(1-)「b+,,+A(b,+fb,)+Ac 2(1+p)(1-21)4 A(c, b,+crfs)+A,b,c (4-33) A,(b, cs +f, c,)+A,c, b C,+a,bb 44载荷移置 单元上的体力为{p},与平面问题相同,由虚功方程可以得到结点载荷, (R =2T[INI(p)rdrd= (4-34) 作用在单元上的面力为何},结点载荷为,单元刚度矩阵为， K B D B rdrdz e T [ ] 2 [ ] [ ][ ]  =  （4-28） 单元刚度矩阵的分块矩阵为， K B D B rdrdz s T rs r 2 [ ] [ ] [ ]  =  （4-29） 由于几何矩阵中的元素不是常量，单元刚度矩阵需要通过积分得到，为简化计算，可以 用三角形单元形心位置的坐标 c c r ,z 代替[B]矩阵中的变量 r、z。 ( ) 3 1 c i j m r = r + r + r ， ( ) 3 1 zc = zi + z j + zm 应变矩阵变成，  m  [B] = Bi Bj B               + + = i i i i i c c i i i c b c b c z r a b A B 0 0 0 2 1 [ ] （4-30） 单元刚度矩阵的近似表达式为： [K] 2 r [B] [D] [B] T c e =  （4-31） 单元刚度矩阵的分块矩阵近似表达式为，   s T K rs = 2rc [B] r [D] [B] （4-32）    + + +    + + + + + + + − − = r s r s r s r s r s r s r s r s c r s r s r s r s r s r s c c A b b A b c f c A c b A c b c f A b c b b f f A b f f b A c c A E r K 2 1 2 1 2 1 2 ( ) ( ) ( ) 2(1 )(1 2 ) (1 ) [ ]     (4-33) 4.4 载荷移置 单元上的体力为{p}，与平面问题相同，由虚功方程可以得到结点载荷，  R = N p rdrdz e T { } 2 [ ] { } (4-34) 作用在单元上的面力为 P ，结点载荷为