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

节点号 Ansys解 程序解 节点号 Ansys解 程序解 节点号 Ansys解 程序解

172 0.0700 0.2483 178 0.1078 0.1069 184 0.0569 0.2950

173 0.1019 0.2184 179 0.0617 0.0854 185 0.0770 0.2027

174 0.1291 0.1692 180 0.5409 1.5036

175 0.1331 0.1140 181 0.4426 0.4233

176 0.1252 0.2197 182 0.1214 0.1903

177 0.1067 0.1035 183 0.1341 0.1553

（注：单位为1.0e4Pa）

对比分析：在沿半径方向数据吻合的很好，误差基本上是很小的；在圆周方向数据基本吻合，但是在一些应力值很小的点存在较大的差别，这可能与ansys处理四节点单元节点应力与程序处理方式差异有关。