非光滑多体系统最大Lyapunov指数的计算方法

摘要

结合非光滑力学与非线性动力学的理论给出非光滑多体系统最大Lyapunov指数的数值计算方法.应用第一类Lagrange方法建立具有单边约束非光滑多体系统变拓扑结构动力学方程;应用线性互补理论与方法检测非光滑事件;应用混沌同步方法计算Lyapunov指数.在提高计算速度方面,提出了采用二分搜索法加快搜索速度以及采用调制参数方法加快系统达到同步状态的速度.通过算例说明本文算法的可行性与有效性.