平面四节点等参单元matlab实现(13)

五、讨论比较

对所选取的节点的mises应力列表做对比： 1.沿半径方向（排列顺序为节点号从小到大） 表1 节点号

24 29 47 58 78 79 137

Ansys解 0.0388 2.2423 0.0969 0.0439 0.0521 0.0702 0.1100 程序解 0.0365 3.4201 0.0458 0.0415 0.0448 0.0673 0.1396 节点号

149

160

166

186

Ansys解 0.1444 0.2330 0.4766 0.5914 程序解 0.1076 0.1859 0.4403 3.2183

（注：单位为1.0e4Pa）

2沿圆周方向（排列顺序按单元从大到小） 表2 节点号 Ansys解 程序解

166 0.4765 0.4403

167 0.2350 0.3125

168 0.0983 0.2060

169 0.0885 0.1686

170 0.0620 0.2256

171 0.2032 0.6079