一种有砟轨道大机捣稳组合作业模式决策方法，包括：采用离散元与多体动力学耦合算法，建立捣固稳定装置‑有砟道床耦合计算模型；现场采集待作业区段的道床支承刚度参数，将测试结果代入模型，使模型的初始道床支承状态与实测结果相拟合；根据现场作业条件，将多种大机捣稳组合作业模式作为不同工况输入至模型中，对大机捣稳组合作业影响下的有砟道床与轨枕的力学行为变化情况进行模拟；基于模型计算出不同大机捣稳组合模式作业下有砟道床状态指标参数；将计算结果代入至道床质量综合评价指数公式中，根据计算结果比选得出最佳大机捣稳组合作业模式。