abaqus計算受重力作用下的土體表面重應力為什么不是零呢，這個是下面的個人一些想法，可以供小白參考。

    舉一個最簡單例子，假設土體大小10X10X10米，材料密度2000kg/m³.彈性模量100Mpa，泊松比0.3，摩擦角30度，粘聚力30Kpa，只受重力作用，重力加速度取10。單元尺寸大小分別取0.5、1、2、5m。

     計算地表豎向應力分別為0.5X10⁴pa、1X10⁴pa、2X10⁴pa、5X10⁴pa，可以看出，單元尺寸越小，地表單元的應力就越小，結果偏于更準確。這是因為重力是作用在每個單元的重心位置，該模型標準矩形，單元也規整，第一層每個單元的標高是單元網格尺寸的一半，第一層重心位置的應力就是密度X重力加速度X該層單元格重心深度，再通過有限元原理轉化到每個單元的節點上，可想而知，要想地表網格節點尺寸為0，必須是單元網格大小足夠小，接近于0，這就是為什么abaqus模擬巖土工程不準確的地方，不可能做到足夠小，一般巖土工程的模型都是比較大的（幾十米幾百米幾千米），模型越大網格尺寸會劃分的很大，精確度也越低。更多案例可以關注抖音abaquser。