您的代码中有几个地方可以查看：

1. 查看$S^e_{1}$，局部刚度矩阵，和$T^e_{1}$，局部质量矩阵，用于节点编号的两种情况下的元素 [1]。– 局部矩阵和的创建方式可能存在错误（如评论中所述，一个常见错误是用于 CW 和 CCW 节点编号$S^e_{1}$$T^e_{1}$$A_1$
2. 确保元素 [1] 的编号不会影响和，其余未触及元素的局部刚度和质量矩阵。$S^e_{2,3,4}$$T^e_{2,3,4}$
3. 从元素 [2]、[3] 和 [4] 开始组装全局刚度矩阵然后，比较您的局部刚度矩阵如何添加到两个编号的全局刚度矩阵中。$S$$S^e_{1}$
4. 比较元素 [1] 的两个编号的 RHS。

我想我和你面临同样的问题。对于顺时针和逆时针本地编号，我有不同的结果。我正在使用三角形元素。我刚刚发现，只要我们正确地将符号放入局部刚度矩阵中，顺时针或逆时针局部编号应该具有相同的结果。我刚刚从参考资料中读到的规则是，如果局部编号方向与全局编号相反，则为局部编号刚度矩阵加上减号，反之亦然。有关详细信息，您可以阅读此参考资料https://drive.google.com/file/d/1Xd2gq8q6REr8DX1fumFxh8kYNaXG1I0o/view?usp=sharing