abaqus是如何处理角位移的？

pacoyang

举刚体转动的一般情况的例子。
刚体参考点，在x,y,z轴（以此为4,5,6自由度）上施加角位移Φ1，Φ2，Φ3，三个角位移均为时间t的函数，即Φ1=Φ1(t),Φ2=Φ2(t),Φ2=Φ2(t)。t=0是Φ1=0，Φ2=0，Φ3=0。


刚体的一般转动，其角位移不是矢量，既不满足可加性，也不满足可交换性。所以有下面的疑问：
在t时刻，角位移分别为Φ1，Φ2，Φ3。在t+Δt时刻，角位移为Φ1+ΔΦ1，Φ2+ΔΦ2，Φ3+ΔΦ3。那么abaqus如何实现角位移增量（ΔΦ1，ΔΦ2，ΔΦ3）的。


猜想有两种方法：
一是近似为无限小位移，该位移具有矢量特性，就把这个角位移增量当做矢量。这是个近似做法（在理论上就是近似，是否等同于微分，不懂！！）。

二是将角位移增量转换为角速度，用角速度反推欧拉方程的三个角度变量(章动角，进动角，自转角),获得准确的结果。

到底abaqus是如何实现转动自由度的控制的？


我们都知道，定点刚体转动，以t=0时刻的初始构型为基础，在t时刻的构型，均可通过初始构型绕某个轴线，转动某个角度来获得。上面的问题等效于：已知角位移与时间的函数关系，如何求得这个轴线和这个角度？

louis_2022

好棒的问题，同问

louis_2022

我想咨询楼主一个问题。
就是abaqus输出的角位移的旋转顺序是怎样的？
我试了可能的6中方式，均不行。。。有点难受