碰撞激发弹塑性波传播的动态子结构方法OALib.PDFVIP

第 卷 第 期 力 学 学 报 Vol. No. 年 月 Chinese Journal of Theoretical and Applied Mechanics ， 碰撞激发弹塑性波传播的动态子结构方法1) 2) 骞朋波，尹晓春 ，沈煜年，杨钧 （南京理工大学力学与工程科学系，南京 210094) 摘要：该文通过建立弹塑性碰撞动态子结构模型，推导了模态坐标下的控制方程，提出了模拟柔性结构碰撞激发 弹塑性波传播的动态子结构方法，并对其中的主模态的存在性和主模态截断的收敛性进行了证明。通过对柔性杆 纵向碰撞和柔性梁横向碰撞两个算例的计算，并将计算结果与理论解和三维动力有限方法计算结果进行了对比， 验证了该方法的数值收敛性和计算碰撞弹塑性波传播的有效性。 关键词：柔性体；碰撞；弹塑性波；动态子结构方法；数值方法 中图分类号： O347.4 ；O313.4 ；O313.7 文献标识码：A 文章编号： 引 言 收敛性。通过数值算例，并与理论解和有限元解的 柔性结构的碰撞响应行为一直受到工程界的关 对比，验证了该动态子结构方法研究碰撞弹塑性波 注[1-2] 。由于碰撞引起的局部弹塑性接触变形和弹性 传播的可行性。 波与塑性波的传播行为等具有强非线性特征，使得 1 动力学方程推导 研究主要依赖于数值和实验分析方法[2-3] 。动态子结 1.1 子结构模型 构方法在上世纪六十年代由 Hurty 等人提出[4]，通 如图 1 所示，考虑两个可能发生碰撞的柔性体 过单独分析子结构的模态和截断子结构主模态，大 B 和 B 。按照一定的规则将B 划分为 d 个子结构， 1 2 1 幅减缩计算自由度，在结构振动分析中得到了广泛 将 B2 划分为 N-d (Nd)个子结构，可能发生碰撞的 的运用。1990 年，Wu 和 Haug[5]首次将动态子结构 结点P 1 和 P2 分别位于第 d 个和 d+1 个子结构。 方法应用于柔性杆和梁的碰撞问题。之后，郭安萍 B 1 N [6] [7] 1 2 B2 N-1 等 做了进一步的研究。最近，沈煜年和尹晓春 ... d-1 ... d P 1 d+2 将该方法应用于柔性杆件碰撞弹性瞬态波传播的研 d+1 P2 究中。但是，到目前为止，动态子结构方法的研究