我正在用Python编写一段代码，用于3D结构的动态模拟。结构由不同方向的节点和元素组成。我正在尝试在全局坐标系中构建整个结构的刚度矩阵，并叠加单元刚度矩阵。我的问题是，对于每个元素，我需要一个从局部轴到全局轴的通用转换矩阵。我已经在互联网上查找过它，但大多数人都做了一些假设来使这种转换更容易（例如，局部 y 轴在全局 X-Y 平面中......我找到了计算变换矩阵的通用方法，但它需要局部坐标系中的第三个节点，我不知道应该如何找到它。
如果有人能在这里帮助我，我真的很感激。
谢谢
本文
的第9页可能会有所帮助： http://elvis.rowan.edu/~kay/papers/kinematics.pdf
首先计算每个元素的齐次变换矩阵（相对于全局坐标系）。然后将其与相应的刚度矩阵相乘，将其转换为全局刚度矩阵。
顺便说一句：http://scicomp.stackexchange.com 或 http://engineering.stackexchange.com 可能是讨论这个问题的更好社区。