In the battery pack assembly, it is essential to ensure that the cell-to-busbar joints can be produced with high quality and with minimal impact on the individual battery cells. This study examines the influence of process parameters on the joint quality for nickel-plated copper and steel plates, laser welded in an overlap configuration. Artificial neural network-based meta models, trained on numerical results from computational fluid dynamics simulations of the laser welding process, are used to predict and evaluate the joint quality. A set of optimized process parameters is identified, in order to simultaneously maximize the interface width for the joints, and minimize the formation of undercuts and in-process temperatures. In an meta model-based multi-objective optimization approach, the non-dominated sorting genetic algorithm II (NSGA-II) is used to efficiently search for trade-off solutions and the meta models are used for objective approximation. As a result, the objective evaluation time is decreased from around 9 h, when evaluated directly from numerical simulations, to only tenths of a second. From the Pareto-optimal front of trade-off solutions, three optimal solutions are selected for validation. The selected solutions are validated through laser welding experiments and numerical simulations, resulting in joints with large interface widths and low in-process temperatures without a full penetration.