含摩擦与碰撞平面多刚体系统动力学线性互补算法

基于接触力学理论和线性互补问题的算法, 给出了一种含接触、碰撞以及库伦干摩擦, 同时具有理想定常约束(铰链约束) 和非定常约束(驱动约束) 的平面多刚体系统动力学的建模与数值计算方法. 将系统中的每个物体视为刚体, 但考虑物体接触点的局部变形, 将物体间的法向接触力表示成嵌入量与嵌入速度的非线性函数,其切向摩擦力采用库伦干摩擦模型. 利用摩擦余量和接触点的切向加速度等概念, 给出了摩擦定律的互补关系式; 并利用事件驱动法, 将接触点的黏滞-滑移状态切换的判断及黏滞状态下摩擦力的计算问题转化成线性互补问题的求解. 利用第一类拉格朗日方程和鲍姆加藤约束稳定化方法建立了系统的动力学方程, 由此可降低约束的漂移, 并可求解该系统的运动、法向接触力和切向摩擦力, 还可以求解理想铰链约束力和驱动约束力. 最后以一个类似夯机的平面多刚体系统为例, 分析了其动力学特性, 并说明了相关算法的有效性.      

车辆纵向非光滑多体动力学建模与数值算法研究

在建立车辆纵向多体系统的动力学模型中, 将车身与车轮视为刚体, 两者通过减振器链接; 将传动系统视为一个圆盘通过扭簧和阻尼器与驱动轮连接; 将车轮与路面间的接触力简化为法向约束力、摩擦力和滚阻力偶,其中摩擦力的模型采用库仑干摩擦模型. 采用笛卡尔坐标作为该系统的广义坐标用于描述该系统的位形, 给出系统单双边的约束方程, 应用第一类拉格朗日方法建立了系统的动力学方程. 由于摩擦与滚阻的非光滑性, 使得该系统动力学方程不连续. 为便于计算, 建立了车轮与路面接触点的相对切向加速度与摩擦力余量的互补条件、车轮角加速度与滚阻力偶余量的互补条件, 以及车轮轮心法向加速度与路面法向约束力的互补条件. 将接触—分离、黏滞—滑移的判断问题转化成线性互补问题的求解, 并给出了具有约束稳定化的基于事件驱动法的数值计算方法. 最后, 应用该方法对车辆纵向多体系统进行了仿真, 分析了输出扭矩、摩擦及滚阻系数对其动力学行为的影响.      

含摩擦滑移铰平面多刚体系统动力学的数值算法

基于LuGre摩擦模型和线性互补问题(LCP)的数值算法,给出了具有双边约束含摩擦滑移铰平面多体系统动力学的数值算法.首先,根据滑移铰的特点,当间隙充分小时,将其视为双边约束,给出了滑移铰中滑道作用于滑块上的法向接触力的互补关系;LuGre摩擦模型能有效地描述机械系统中的黏滞与滑移运动,将该模型用于描述滑块与滑道间的摩擦力.其次,结合Baumgarte约束稳定化方法,应用第一类Lagrange方程,建立了该多体系统的动力学方程,给出了Lagrange乘子与滑移铰中作用于滑块上的法向接触力的关系式.然后,将滑块与滑道间多种接触状态的判断以及作用于滑块上的法向接触力的计算转换为线性互补问题的求解,并用常微分方程的数值算法求解该多体系统的动力学方程.最后,通过数值仿真算例揭示了滑移铰中滑块的黏滞与滑移现象,以及滑块在滑道内的多种接触状态;另外,在文中分别采用Coulomb干摩擦模型和LuGre摩擦模型,对算例中的某些工况进行了数值仿真,并且分别用本文方法得到的数值仿真结果与已有方法得到的数值仿真结果对比,表明了本文给出的方法的有效性.      

非光滑多体系统动力学数值算法的研究进展

非光滑多体系统动力学数值计算方法是多体系统动力学研究的重要内容之一. 本文介绍了近年来含摩擦与碰撞的非光滑多体系统动力学数值算法方面的研究进展. 首先, 讨论了库仑摩擦模型和修正的库仑摩擦模型, 以及具有单边和双边约束的多体系统中法向约束力的特点. 其次, 回顾了基于连续模型和非连续模型的多体系统动力学方程的数值计算方法, 详细介绍了基于互补概念的非光滑多体系统动力学的事件驱动法和时间步进法, 分析比较了相关的数值算法. 最后, 指出了一些需要进一步研究的问题.     

空间可展机构非光滑力学模型和动力学研究

空间可展机构广泛应用于展开和支撑柔性太阳能帆板和航天工程领域中的有效载荷, 包括抛物面天线、平面相控阵雷达和合成孔径雷达等. 非光滑特性及其相应的动力学现象在空间可展机构的设计中有着非常重要的作用. 该文系统地综述了空间可展机构非光滑力学建模与非线性动力学的研究进展. 首先详细描述了含间隙铰链的接触碰撞力和摩擦力等非光滑特点;然后系统地介绍了含间隙机构的动力学建模方法、分析方法和参数设计;进一步简单介绍了含间隙铰链空间可展机构的非线性动力学特性, 如谐波共振、周期运动的稳定性和各类分岔等;最后提出了空间可展机构非光滑动力系统动力学、稳定性与控制中亟待解决的若干问题.     

圆弧足被动行走器非光滑动力学仿真研究

传统的被动行走器动力学研究方法是将被动 行走器模型视为双足轮流铰接于地面(或纯滚动)的倒立摆模型, 并应用角动量守恒假设计 算支撑腿转换时足地碰撞后的系统运动状态. 但该方法不能分析足地之间的滑动现 象、足地之间的非完全塑性碰撞现象和两足同时着地的情况. 采用基于线性互补 的非光滑动力学方法, 给出了一种基于中点格式积分一步的时间步进算法, 可以有效地分析 足地间的摩擦和碰撞现象. 应用该方法对被动行走器进行数值仿真, 发现了两种新的周期性 步态, 一种是带摩擦滑动现象的步态, 另一种是支撑足转换时发生多次碰撞现象的步态. 数 值仿真结果还表明, 较大的摩擦系数和较小的法向恢复系数更有利于被动行走器实现稳定周 期步态行走.     

链轮驱动系统的多刚体建模及分析

采用多刚体方法对链条链轮系统的驱动过程进行了建模和动力学分析．链节看作刚体,链轮看作等边刚性多边形,辊子看作两个链节之间的旋转副．啮合约束当作可变约束处理,当啮合力为正时,链节和链轮之间是固定约束,当力变负时,链节脱离链轮．利用龙格库塔法求解系统的拉格朗日方程．最后计算并分析了不同工况下的驱动过程,解释了系统的锁死现象,并进一步分析了两种避免锁死的处理方法．     

井筒内受压杆管后屈曲能量法分析与实验研究

目前不少理论力学教材都涉及刚体平面运动动力学方程的教学内容,但有些教材在叙述上有许多值得商榷之处.对这一重要的基本概念,本文提出作者的思考,供理论力学教学的教师与学过此内容的学生讨论.     

EXPERIMENTAL INVESTIGATION ON THE EVOLUTION OF INTERNAL SOLITARY WAVE OVER A RIDGE AND ITS ACTING FORCES ON A SUBMERGED SLENDER BODY

内孤立波沿山脊地形传播过程中的结构变化使得影响海洋结构物及水下航行器安全的不确定因素增加,用大型重力式分层流试验水槽系统研究内孤立波过山脊地形及其对细长潜体作用力特性,可有效提高对其危害性机理的认识.为此,利用电导率探头阵列,结合染色标识方法,测量内孤立波演变特性,同时利用三分量传感器测量水下细长体模型受力特性.研究结果表明：山脊地形显著改变下凹型内孤立波结构,表现为坡前波幅增大,坡顶背风波面抬升,坡后波长变长;内波结构变化直接影响细长体受力特性,表现为坡前向下作用力明显增强,坡顶向上作用力突显,坡后作用力持续时间变长;细长体受力特性变化影响其运动趋势,坡前下沉运动增强、朝山脊方向纵荡运动显著.     

冲压过程模拟的接触与摩擦关系数值方法

文中简介采用非经典算法求解板件冲压过程数值模拟的基本思想；就与模具间接触与摩擦关系的处理作重点论述，并给出算例。     

动、静摩擦对可接触型裂纹动态影响的一种算法

本文提出了一种分析考虑裂纹面摩擦影响的可接触型裂纹动力反应的时域显式计算方法，这一方法不但可以较好描述裂纹面之间碰撞过程，也能反映当裂纹面发生接触－滑移时，裂纹面间动、静摩擦力的影响．文中给出了使用这一方法的基本步骤，计算了无限大裂纹的动接触问题和直线型裂纹对入射Ｐ波的散射问题，初步验证了本文提出的方法．