This paper investigates an integrated production and distribution scheduling problem in a hybrid flow shop with a batch delivery system. The goal is to schedule jobs and assign them to batches to minimize the total weighted number of tardy jobs and delivery costs, thereby improving both operational efficiency and customer satisfaction. To achieve this, a novel mixed-integer linear programming model is developed. Given its computational complexity for large instances, two metaheuristic algorithms, simulated annealing (SA) and a genetic algorithm (GA), are proposed and calibrated using the Taguchi method to enhance performance. Experimental results demonstrate that simulated annealing consistently outperforms the genetic algorithm in both solution quality and computational time. The results also demonstrate that the proposed approach reduces the overall cost of production and distribution by 10.45%.
