PRINT sel pile local，axis下x‘，y'，z'表示pilesel的局部坐标，
如第一行，x'：(5.175e-8,-2.561e-7,1.000) 表明桩单元局部坐标的x轴在全局坐标系中的向量为(5.175e-8,-2.561e-7,1.000)，
y'，z'同样的道理
从你图中看，要查看绕y轴的弯矩，而y'平行于y轴，所以应查看my
一般来说，这种竖直桩，局部坐标系应该很精确，（0，0，1），（0，1，0），（1，0，0），你这个为什么有很小的偏差，不知道什么原因
边界条件施加方法是对的，1/2对称处也只约束边界面法向位移即可（可能有特殊情况，视情况而定）
mz出现那种情况是否跟横支撑、腰梁的连接有关，
另外模拟挡土桩，可以试试beam单元，pile单元参数较多而且很多不好确定（个人看法）
多做点试算工作，看看英文手册。
