本发明公开了一种六辊冷轧机的轧制稳定性预测方法，涉及轧制过程自动化生产技术领域。该方法考虑了入口油膜挤压效应，将轧辊垂向振动速度引入油膜厚度计算公式，获得动态入口油膜厚度，并结合粗糙度分布假设，计算变形区摩擦应力分布随时间的变化情况；考虑了轧辊垂向振动速度的卡尔曼微分方程推导，并带入变形区摩擦应力分布，计算动态轧制力及由轧辊垂向振动引发的轧制力波动量；根据轧辊、轧件和牌坊间的受力关系，建立轧机系统的垂向振动动力学方程，然后采用Newmark‑Beta法求解，并以轧辊垂向位移作为判断轧机稳定性的依据，若轧辊位移曲线收敛，则轧机稳定，若轧辊位移曲线发散，则轧机不稳定。该方法能够更精准地预测出轧制过程中的轧机稳定性。