abaqus计算受重力作用下的土体表面重应力为什么不是零呢，这个是下面的个人一些想法，可以供小白参考。
举一个最简单例子，假设土体大小10X10X10米，材料密度2000kg/m3.弹性模量100Mpa，泊松比0.3，摩擦角30度，粘聚力30Kpa，只受重力作用，重力加速度取10。单元尺寸大小分别取0.5、1、2、5m。
计算地表竖向应力分别为0.5X104pa、1X104pa、2X104pa、5X104pa，可以看出，单元尺寸越小，地表单元的应力就越小，结果偏于更准确。这是因为单元的节点应力是单元积分点上通过积分推导出来的，重力是作用在每个单元的重心位置，该模型标准矩形，单元也规整，第一层每个单元的重心标高是单元网格尺寸的一半，第一层重心位置的应力就是密度X重力加速度X该层单元格重心深度，再通过有限元原理转化到每个单元的节点上，可想而知，要想地表网格节点尺寸为0，必须是单元网格大小足够小，接近于0，这就是为什么abaqus模拟岩土工程不准确的地方，不可能做到足够小，一般岩土工程的模型都是比较大的（几十米几百米几千米），模型越大网格尺寸会划分的很大，精确度也越低。
因篇幅问题不能全部显示，请点此查看更多更全内容