abaqus网格密度到底取多少最合适？ 70

我现在在做一个悬臂梁理论解和有限元解的一个对比。

悬臂梁长度200mm，截面是20mm的矩形，在顶面施加均布荷载1MPa。符合材料力学杆模型，理论解如下图

有限元中边界条件和荷载设置如下：

我用abaqus有限元建模，材料属性设置的是400MPa，塑性应变为0。采用C3D8R单元，采用不同的网格密度计算出来的弯矩引起的最大应力水平分别见下图所示，S33就是弯矩引起的正应力。

随着网格密度不断细化，网格密度为2.5mm的时候和理论解是最接近的，但是网格细化到1mm的时候，结果却比理论解大得多？请问这是为什么呢？我之前见有人说是网格越细，会导致刚度越大？不知道是不是和这个有关系？但是石亦平老师书里讲过网格越细，应当越接近理论解。这个不知道该怎么解释了。（可视化中显示的是节点应力，但是我在积分点中也查看了，结果和上面显示的最大值差不多）