6. 0 f10 [B;]= (下标轮换) (424) 2A0 Ci 由于在f是坐标r、z的函数,E分量在单元中不为常量,其它三个应变分量在单元 中仍为常量。 由轴对称问题的物理方程,得到弹性矩阵 1-1- E(1-4)1 (1+)(1-21) 0 0 0 2(1-) 令 2(1 )42,则弹性矩阵为, 1A1A10 [D]= E(-)4140 (4-26) (1+)(1-2)4141 0 000A 由弹性矩阵[D]和几何矩阵[B可以得到应力矩阵[S],并计算出单元内的应力分量, [S]=[D[B (4-27) S=Si Si smI b+f141 S;]=[DIB]= E(1-)4b+f (4-27) 2(1+1)(1-2u)4A1(b+f)c A2c A2b 下标轮换,可得到[S;1[S 由应力矩阵可知,除剪应力rx为常量,其它三个正应力分量都是r、z的函数。            = i i i i i i c b c f b A B 0 0 0 2 1 [ ] （下标轮换） （4-24） 由于在 i f 是坐标 r、z 的函数，   分量在单元中不为常量，其它三个应变分量在单元 中仍为常量。 由轴对称问题的物理方程，得到弹性矩阵，                       − − − − − − − − + − − = 2(1 ) 1 2 0 0 0 1 0 1 1 0 1 1 1 0 1 1 1 (1 )(1 2 ) (1 ) [ ]                 E  D （4-25） 令 1 1 = A −   ， 2 2(1 ) 1 2 = A − −   ，则弹性矩阵为，             + − − = 2 1 1 1 1 1 1 0 0 0 1 0 1 0 1 0 (1 )(1 2 ) (1 ) [ ] A A A A A A A E D    （4-26） 由弹性矩阵[D]和几何矩阵[B]可以得到应力矩阵[S]，并计算出单元内的应力分量， [S] = [D][B] （4-27） （4-28） [ ] [ m] S = Si S j S             + + + + − − = = i i i i i i i i i i i i i A c A b A b f c A b f A c b f A c A E S D B 2 2 1 1 1 1 2(1 )(1 2 ) ( ) (1 ) [ ] [ ][ ]    （4-27） 下标轮换，可得到 [ ], [ ] Si S j 。 由应力矩阵可知，除剪应力 zr  为常量，其它三个正应力分量都是 r、z 的函数