我現在在做一個懸臂梁理論解和有限元解的一個對比。

懸臂梁長度200mm，截面是20mm的矩形，在頂面施加均布荷載1MPa。符合材料力學桿模型，理論解如下圖

有限元中邊界條件和荷載設置如下：

我用abaqus有限元建模，材料屬性設置的是400MPa，塑性應變為0。采用C3D8R單元，采用不同的網格密度計算出來的彎矩引起的最大應力水平分別見下圖所示，S33就是彎矩引起的正應力。

隨著網格密度不斷細化，網格密度為2.5mm的時候和理論解是最接近的，但是網格細化到1mm的時候，結果卻比理論解大得多？請問這是為什么呢？我之前見有人說是網格越細，會導致剛度越大？不知道是不是和這個有關系？但是石亦平老師書里講過網格越細，應當越接近理論解。這個不知道該怎么解釋了。（可視化中顯示的是節點應力，但是我在積分點中也查看了，結果和上面顯示的最大值差不多）