含碰撞的平面摩擦系统半解析半数值算法研究

利用含平面库仑摩擦力系统的解析解结果,通过碰撞边界条件和时间的隐式关系,获得精确的碰撞时刻,由此确定从当前时刻到碰撞发生时的步长及平面上的碰撞位置。结合滑移/黏滞转换时位置和时刻的解析解,引入布尔变量B定义系统的状态,建立一种求解含碰撞的平面摩擦系统响应的数值算法。研究结果表明:该算法能够求解滑移/黏滞转换、碰撞以及平面滑移运动存在方向剧烈变化拐点时系统的响应,并有效地提高了系统响应的精度。