CEL方法计算静压桩遇到了问题

韦应侯爵

最近在做静压桩数值模拟，用的是CEL方法，土体采用的是摩尔-库伦模型，密度1.8t/立方米，弹模10000kpa，摩擦角30，凝聚力10kpa，桩采用的是各向同性体，底部设置了刃角。土体材料是欧拉体，桩是拉格朗日体。欧拉部分，上面留了一米高的VOID层。计算出来以后，提取桩结构上下两个截面的轴力，发现值特别小，不到一个kn，已经困扰我好多天了，不知道有没有大神遇到过类似的问题，真的迫切需要帮助。