1.一种接触界面的仿真优化方法，应用于碰撞体的运动状态仿真，其特征在于，所述方法包括：初始化输入参数，所述输入参数包括碰撞体接触界面之间的摩擦系数 ；

运动状态求解：根据初始化的 求解碰撞体之间的相互摩擦力 ，以及由 引起的摩擦力冲量 和摩擦导致的切向加速度 ，利用摩擦力冲量 和切向加速度 更新碰撞体的位置和速度；

摩擦系数优化：利用损失函数优化摩擦系数 ，获得摩擦系数 ；

利用摩擦系数 再次进行运动状态求解和摩擦系数优化，获得摩擦系数  ，利用摩擦系数 再次进行运动状态求解和摩擦系数优化，获得摩擦系数 ……如此循环迭代，直至所述损失函数收敛于最小值；

其中，所述损失函数包括输出损失函数 ：

；

对于各次迭代的摩擦系数满足：

其中， ，为学习率，用于控制更新的步长； 表示第一碰撞体实际停下的位置与预期位置的欧式距离， 表示第二碰撞体实际停下的位置与预期位置的欧式距离； 表示在 时刻碰撞体的状态变量，包括碰撞体在 时刻的位置 和速度 ，，表示整个控制时间区间的终止时间，表示迭代次数。

2.根据权利要求1所述的接触界面的仿真优化方法，其特征在于：所述损失函数还包括输入损失函数以及总损失函数，输入损失函数为：其中，为权重常数， 表示在连续时刻的外部控制输入， ， 表示在时刻碰撞体的状态变量，包括碰撞体在时刻的位置 和速度 ；总损失函数为：总损失函数 的约束条 件为： ， ， ，

；

其中， 表示求最小值； 表示系统状态的时间导数， ； 是系统的动力学方程，描述了在外部控制输入 作用下，状态 的演变； 表示发生第 次碰撞的时间， 表示第 次碰撞发生后碰撞体的状态， 表示第 次碰撞前碰撞体的状态，函数 表示碰撞对碰撞体状态的影响， 用于描述碰撞之后碰撞体状态的变化； 是碰撞检测函数，当 满足 时，表示即将在 时刻发生碰撞；

将所述总损失函数和约束条件进行离散化，获得离散化的总损失函数：}；

以及离散化约束条件： ；

其中， 为离散化时间步长， ； 表示系统在第N个时间步长内的状态；表示系统在第 个时间步长内的状态， 表示在第 个时间步长内的外部控制输入， 表示状态更新函数。

3.根据权利要求2所述的接触界面的仿真优化方法，其特征在于，所述摩擦系数优化还包括：新增损失函数，以优化碰撞体与地面之间的摩擦系数，其中，新增损失函数为：是第 个碰撞体在仿真中最后一步的速度，为权重常数，优化目标为损失函数收敛于其最小值 ；

新增另一损失函数，以优化碰撞体与墙体之间的摩擦系数，该新增损失函数为：其中，表示仿真过程中共发生了 次碰撞事件， 表示发生第 次碰撞后碰撞体的实际速度变化量， 表示发生第 次碰撞后对碰撞体预期的速度变化量,优化目标为损失函数 收敛于其最小值 。

4.根据权利要求3所述的接触界面的仿真优化方法，其特征在于，利用初始化后的摩擦系数 进行运动状态求解包括以下过程：利用连续碰撞检测CCD计算待仿碰撞体的穿透时间 ；

向前步进一个穿透时间 ，以计算待仿碰撞体的碰撞相关参数，所述碰撞相关参数包括碰撞位置 、碰撞速度 和碰撞方向向量 ；

根据所述碰撞相关参数计算待仿碰撞体所受的推力冲量 ；

根据所述推力冲量 以及初始化后的摩擦系数 计算待仿碰撞体之间的相互摩擦力 ，并确定摩擦力 的方向；

计算由摩擦力 对所述待仿碰撞体所造成的冲量 ；

摩擦力冲量 和摩擦导致的加速度 被用于更新待仿碰撞体的位置和速度，这个过程在每个时间步长中重复执行。

5.根据权利要求4所述的接触界面的仿真优化方法，其特征在于：所述穿透时间 根据下式计算：其中，d为穿透深度， 为第二碰撞体在第 个时间步长时的穿透速度， 为第一碰撞体在第 个时间步长时的穿透速度，表示穿透方向；穿透深度d根据下式计算：其中， 表示第二碰撞体在第 个时间步长时的穿透位置向量， 表示第一碰撞体在第 个时间步长时的穿透位置向量， 表示求向量的模；第一碰撞体和第二碰撞体均为球体，表示第一碰撞体的半径，表示第二碰撞体的半径；

其中，穿透方向 根据下式计算：

；

对于第一碰撞体 ，其穿透位置向量 通过辛欧拉积分方式进行计算，如下：其中， 为第一碰撞体在第 个时间步长时的位置向量， 为第一碰撞体在第个时间步长时的穿透速度， 表示一个时间步长； 表示第一碰撞体在第 个时间步长时的速度， 表示第一碰撞体在第 个时间步长时的加速度；

对于第二碰撞体，其穿透位置向量 通过辛欧拉积分方式进行计算，如下：其中， 为第二碰撞体在第 个时间步长时的位置向量， 为第二碰撞体在第个时间步长时的穿透速度； 表示第二碰撞体在第 个时间步长时的速度， 表示第二碰撞体在第 个时间步长时的加速度。

6.根据权利要求5所述的接触界面的仿真优化方法，其特征在于，所述推力冲量 根据下式计算：其中，为弹性系数， 是第一碰撞体与第二碰撞体之间的相对速度在碰撞方向上的分量， 为第一碰撞体的质量， 为第二碰撞体的质量，碰撞方向向量 根据下式计算：其中， 表示第二碰撞体在第 个时间步长时的碰撞位置向量， 表示第一碰撞体在第 个时间步长时的碰撞位置向量；

其中，碰撞位置向量 根据下式计算：

其中， 为第二碰撞体在第 个时间步长时的碰撞速度；

其中，碰撞位置向量 根据下式计算：

其中， 为第一碰撞体在第 个时间步长时的碰撞速度。

7.根据权利要求6所述的接触界面的仿真优化方法，其特征在于：所述碰撞体之间的相互摩擦力 根据下式计算：其中， 为法向压力，根据下式计算：

其中，为法向冲量，根据下式计算：

；

确定摩擦力 的方向为相对移动速度的相反方向：；

其中，“”表示向量点乘，（ ）表示第一碰撞体与第二碰撞体的相对移动速度， 为切线方向向量，表示为：其中， 为向量 的 向分量，为向量 的 向分量；

摩擦力冲量 根据下式计算：

为 的符号函数。

8.根据权利要求7所述的接触界面的仿真优化方法，其特征在于，所述摩擦导致的加速度 根据下式计算：其中， 表示在 时刻对第 个碰撞体的摩擦冲量， 表示第 个碰撞体的质量，所述 包括碰撞体之间的摩擦冲量 。

9.根据权利要求8所述的接触界面的仿真优化方法，其特征在于，所述方法还包括：计算碰撞体与地面之间的摩擦力大小 ：

当碰撞体的速度小于预定速度 时，采用上述公式（1）计算摩擦力大小 ；当碰撞体的速度大于等于预定速度 时，采用上述公式（2）计算摩擦力大小 ；其中，为静摩擦系数， 为动摩擦系数， ，为碰撞体的质量，为重力加速度，预定速度 用于决定静摩擦和动摩擦的分界线；

计算摩擦力方向向量：

计算碰撞体与地面之间的摩擦矢量 ：

计算由 引起的摩擦冲量 ：

利用摩擦冲量 再次更新物体的位置和速度，这个过程在每个时间步长中重复执行。

10.根据权利要求9所述的接触界面的仿真优化方法，其特征在于，所述方法还包括：计算碰撞体与阻挡物之间的摩擦力 ，其中，所述阻挡物包括墙体，摩擦力 根据下式计算：其中， 为碰撞体与墙体之间的摩擦系数， 表示碰撞体与墙体碰撞时作用在墙体上的力的法向分量， 表示摩擦力的方向，是沿X轴的单位矢量， 和 的乘积决定了摩擦力在X方向上的作用效果；

法向分量 根据下式计算：

其中， 表示物体与墙体碰撞产生的瞬时冲量； 表示碰撞体与墙体之间造成的穿透时间；

根据下式计算：

其中， 是碰撞体在时刻的速度向量在 方向上的分量；

计算由 引起的摩擦冲量 ：

利用摩擦冲量 再次更新物体的位置和速度，这个过程在每个时间步长中重复执行。