DEM的静态力学解看起来很奇怪

我自己的DEM(离散元法)程序解决静态问题的方法似乎很奇怪。

我用我的程序来计算静压下悬臂梁的位移:

1.似乎我的程序显示收敛的静态问题，网格得到密集。随着网格密度的增大(这意味着颗粒变小)，自由尖端的位移(y方向)变小。

2.当质点较大时，自由尖位移大于理论解，但当质点足够小时，位移逐渐小于理论解。

3.计算解与理论解之间的差异约为4%(当粒子足够小时)