摘要：

以三次非线性缓冲包装系统为研究对象，建立系统跌落冲击动力学方程，应用牛顿谐波平衡（NHB）法获得系统跌落冲击动力学响应的一阶、二阶近似解析解，并获得产品位移、加速度响应最大值以及跌落冲击时间等重要参数的解析表达式。与龙格-库塔（R-K）数值分析解和变分迭代（VIM）近似解比较，算例分析表明，牛顿谐波平衡法二阶近似解与龙格-库塔数值解最为接近，随跌落高度增加二阶近似解精度有所降低，但相对误差控制在2.5%以内，能满足工程设计需要。研究结果可为非线性包装系统跌落冲击响应分析提供一种新的高精度近似分析方法。

关键词: 振动与波, 牛顿谐波平衡, 非线性, 位移最大值, 加速度最大值, 跌落冲击时间