柔性框架结构动力非线性分析的刚体准则法

大跨、高层等柔性结构,其动力响应往往表现出大位移、大转动等非线性特征。动力非线性问题的分析关键在于运动方程的高效稳定求解,以及单元大转动产生的结点力增量的有效计算。动力时程分析通常采用直接积分法,但对于强非线性动力问题,直接积分法难以兼顾计算精度与稳定性。该文基于几何非线性分析的刚体准则,针对杆件结构大转动小应变的非线性问题,提出了一种新型空间杆系结构动力非线性分析的刚体准则法。该方法采用满足刚体准则的空间非线性梁单元,结合HHT-α法求解结构运动方程,并将刚体准则植入动力增量方程的迭代求解过程以计算结点力增量。通过典型柔性框架算例结果表明,该文方法可以有效分析柔性框架结构的强动力非线性行为。与高精度单元相比,该文采用的单元刚度矩阵构造简明,计算过程简洁;与商业软件所用方法相比,单元数和迭代步少,精度高,适于工程应用。

Hilber-Hughes-Taylor-α法在接触约束多体系统动力学中的应用

以柔性梁在重力作用下绕转动铰做大范围定轴转动,并与刚性平面发生碰撞这一动力学过程为例,对Hilber-Hughes-Taylor(HHT-α)法在求解含接触约束的柔性多体系统动力学方程时的数值特性进行了研究.系统运动过程的全局动力学仿真由常微分方程组和微分-代数方程组的数值求解构成.柔性梁在无碰撞阶段系统动力学方程是一组常微分方程组.采用接触约束法模拟接触约束过程,系统的动力学方程为指标3的微分-代数方程组.采用HHT-α法对的该微分-代数方程组进行求解,并与Baumgarte违约修正法进行比较.分析了HHT-α法自由参数和违约修正常数对计算效率、动力学响应和系统机械能的影响,并对数值积分方法对模态截断数的敏感度以及速度约束和加速度约束的违约程度进行了分析.结果表明,违约修正常数对仿真结果影响非常明显,而HHT-α法的自由参数α对动力学响应的影响较小,从而避免了违约修正常数对数值积分结果的影响.HHT-α法的自由参数α可以消除碰撞高频模态的影响.