关于外壳元素的UEL/VUEL

你好，我正在尝试编程一个四边形经典Mindlin平壳单元(对膜的完全积分，对横向剪切的弯曲和选择性简化积分)。

现在的情况是:
对于纯张力测试，它工作良好(因为只有膜部分工作)，
在纯弯曲测试中，与S4R相比，它的结果非常差，我发现这是因为每一步的增量旋转不同。

我隐约记得，从abaqus传递到vuel的增量旋转是用节点局部坐标表示的。但是我没有从手册或说明中找到任何关于用户外壳元素的信息。

1.有人成功地用vuel/uel编程了一个shell元素吗?旋转落点有什么特殊程序吗?
2.旋转是全局坐标还是局部坐标?

欢迎提出任何意见。