基于绝对节点坐标法的弹性线方法研究

相比于浮动坐标系法,绝对节点坐标法(absolute nodal coordinate formulation, ANCF)在处理柔性体非线性大变形问题上具有显著优势, ANCF将单元节点坐标定义在全局坐标系下,采用斜率矢量代替节点转角坐标,具有常数质量阵,不存在科氏离心力等优点,然而弹性力阵为非线性项,其求解将比较耗时且占用资源.据此,在弹性力求解方法中,引入弹性线方法 (elastic line method, ELM),该方法将格林–拉格朗日应变张量定义在中心线上,采用曲率公式来定义弯曲应变,转角公式来定义扭转应变.同时采用有限元法对三维柔性梁位移场进行离散,求解梁单元常数质量阵、广义刚度阵、广义力阵,进而得到单元的动力学方程,通过转换矩阵得到三维梁的动力学方程.接着从理论上指出连续介质力学方法 (continuum mechanics method, CMM)和弹性线方法在求解弹性力上的不同点,并编制动力学仿真软件.最后分别采用连续介质力学方法和弹性线方法对柔性单摆以及履带式车辆的动力学问题进行仿真分析,结果表明:弹性线方法能在保证精度的前提下有效提高计算效率.     

考虑端部质点的中心刚体-楔形梁系统动力特性研究

对带有端部质点的中心刚体-楔形梁系统的动力学特性进行了研究．楔形梁为Euler Bernoulli 梁,其宽度与高度随着长度方向改变．采用假设模态法对楔形梁的变形场进行离散,在考虑梁的横向和纵向变形以及由于横向弯曲而引起的纵向缩短量的条件下,运用第二类拉格朗日方程推导出系统的动力学方程,并编制动力学仿真软件．通过仿真算例对系统的动力学特性进行研究,结果表明：同等条件下梁的形状与有无端部质点对系统的刚柔耦合动力学响应以及系统频率的影响十分明显．因此采用楔形梁结构进行仿真更具实际意义且仿真结果更加精确．     

基于绝对节点坐标法的弹性线方法研究

相比于浮动坐标系法, 绝对节点坐标法(absolute nodal coordinateformulation, ANCF)在处理柔性体非线性大变形问题上具有显著优势,ANCF将单元节点坐标定义在全局坐标系下,采用斜率矢量代替节点转角坐标, 具有常数质量阵,不存在科氏离心力等优点, 然而弹性力阵为非线性项,其求解将比较耗时且占用资源. 据此, 在弹性力求解方法中,引入弹性线方法(elastic line method, ELM),该方法将格林--拉格朗日应变张量定义在中心线上,采用曲率公式来定义弯曲应变, 转角公式来定义扭转应变.同时采用有限元法对三维柔性梁位移场进行离散,求解梁单元常数质量阵、广义刚度阵、广义力阵,进而得到单元的动力学方程, 通过转换矩阵得到三维梁的动力学方程.接着从理论上指出连续介质力学方法(continuum mechanics method,CMM)和弹性线方法在求解弹性力上的不同点, 并编制动力学仿真软件.最后分别采用连续介质力学方法和弹性线方法对柔性单摆以及履带式车辆的动力学问题进行仿真分析,结果表明:弹性线方法能在保证精度的前提下有效提高计算效率.      

基于LuGre摩擦模型的接触约束法旋转柔性梁斜碰撞研究

基于接触约束法和LuGre摩擦模型对在重力场作用下作大范围旋转运动的柔性梁系统和斜坡发生含摩擦斜碰撞的动力学问题进行研究. 首先运用刚柔耦合的多体系统动力学理论对大范围运动的柔性梁进行离散化和动力学建模, 在碰撞时采用冲量动量法求出跳跃速度, 其次在法向上引入接触约束求解出碰撞力, 在切向上采用LuGre摩擦模型分两种方式求解摩擦力, 第一种是在滑动时摩擦力由摩擦系数和碰撞力计算得出, 黏滞状态下引入切向约束计算拉格朗日乘子反应实际摩擦力, 根据黏滞/滑动切换判断计算出碰撞过程摩擦力(与Coulomb摩擦模型计算摩擦力一致); 第二种根据LuGre摩擦模型摩擦系数和法向碰撞力计算其摩擦力, 从而在碰撞时无需黏滞/滑动切换, 采用相同的摩擦力计算公式. 通过与Coulomb摩擦模型对比发现, LuGre摩擦模型描述碰撞切向摩擦过程更精确, LuGre摩擦模型黏滞时建立约束方程和碰撞采用统一的摩擦力公式这两种建模方式描述的斜碰撞动力学特性没有区别, 进而说明采用法向接触约束和LuGre摩擦模型具有满足碰撞非嵌入情况、避免黏滞/滑动切换、描述摩擦力相对准确的优势.