3维有限元后处理，求不标高的弯矩，代码

我现在在做 有限元分析后处理，模型如下图，是桩土相互作用的有限元分析模型，有以下两个问题
1、在ABAUS中，计算桩身弯矩的原理与公式是什么？是如何计算桩身不同标高处的弯矩的？
2、我一直桩身各节点编号，节点坐标，节点力；各单元编号，单元正应力，剪应力，材料编号，请问如何计算桩身不同标高处的弯矩？

需求如下：
1、根据我给的网格文件（节点编号、节点坐标、单元编号、单元材料、单元类型等）、计算结果文件，写个程序，来提取不同标高处的弯矩与转角时程曲线
2、可以自行输入不同标高位置
3、可以根据单元材料属性，选择性输出（只输出材料为钢材---桩单元的弯矩，转角）
4、代码提出的弯矩、转角可以与ABAQUS对应结果
5、可以输入不同单元编号、输出单元孔压比
6、弯矩、转角、孔压比的公式与理论我可以提供，如果有更严谨专业的，也更好