如何用Abaqus隐式动力学画出李萨如图形？

李萨如图形是正交方向上满足一定频率比和相位差的两个简谐振动合成的规则、稳定的闭合曲线，在大学物理的电学实验课堂上，老师会让我们用示波器调出这种图形。

在一些机械振动或某些航天器的轨道中，也能看到它的影子，比如李萨如沙摆，沙子从运动的摆锤中漏出，绘制出规则的图形。

沙摆原理

沙摆在Y型线的作用下，绿色和红色两个正交方向上的摆长不同，频率也就不同，受到初始扰动后，摆锤轨迹可以合成为典型的李萨如图形。

李萨如沙摆的模拟可以用Abaqus显式动力学，考虑到我们的显式分析应用案例比较多，这期文章我们换一种方法，使用Abaqus隐式动力学来计算这个过程。

Abaqus隐式动力学使用隐式时间积分（Hilber-Hughes-Taylor算法、向后的欧拉算法）来计算系统的瞬态动力学或准静态响应，首先简单的看一下Abaqus中隐式动力学的几种应用方案：

• 瞬态保真（Transient fidelity）

不含接触模型的默认选项，涉及最小的系统能量耗散，比如卫星系统的分析，使用较小的时间增量来精确求解结构的振动响应，本文的沙摆采用这种方案。