6}=1w2 (4-6) 在单元中由虚位移引起的虚应变为, {E"}=[B]{G6"}° (4-7) 单元中的实际应力为 {G}°=[DB]{G}2 (4-8) 离散后的单元组合体的虚功方程为, }{F}=∑(6}°) B DIBS dxdydz (4-9) (F)=2(6))'J. [B][DILB]dxdyd-48 (4-10 K=.BDJB就是单元刚度矩阵 对于轴对称问题 [K]° [B] [DIB]rdrdEd0=27 [B] [DIB]rdrd= (4-11) 将(4-11)代入(4-10)可得 6}{F}=('}∑(GKrG)o (4-12) K]=∑(Kr(G)为整体刚度矩阵,得到方程组, [K]{G}={F (4-13)                      = * * * 2 * 2 * 1 * 1 * ... { } n n w u w u w u  （4-6） 在单元中由虚位移引起的虚应变为， e e { } [B]{ } * *  =  （4-7） 单元中的实际应力为， e e {} = [D][B]{} （4-8） 离散后的单元组合体的虚功方程为，  = = n i e T T e e T F B D B dxdydz 1 * * { } { } ({ } ) [ ] [ ][ ]{} （4-9） e n i T e T e T { } {F} ({ } ) [B] [D][B]dxdydz{ } 1 * *      = = （4-10） K B D B dxdydz T e e [ ] [ ] [ ][ ]  = 就是单元刚度矩阵。 对于轴对称问题， K B D B rdrdzd B D B rdrdz e T T [ ] [ ] [ ][ ] 2 [ ] [ ][ ] 2 0    =  =   （4-11） 将（4-11）代入（4-10）可得 { } { } { } ([ ] [ ] [ ]){ } * *  =    e T T T e F G K G （4-12） =  e T e [K] ([G] [K] ([G]) 为整体刚度矩阵，得到方程组， [K]{} = {F} （4-13）