This study aims to achieve the co-optimization of thrust force and thrust fluctuation using a long secondary double-sided linear flux switching permanent magnet motor (LSDLFSPM). Firstly, the motor model is constructed and derived using a theoretical approach. Subsequently, the motor parameters are subjected to sensitivity analysis using the Taguchi method to identify the significant influencing factors. Based on the screening results, the Response Surface Method (RSM) is employed to construct the test space and derive regression equations for thrust force and thrust fluctuation. The Multi-Objective Grasshopper Optimization Algorithm (MOGOA) is then utilized to iteratively optimize the regression equation for optimal parameter sizes. Finally, the optimized results are validated through finite element analysis (FEA) and compared with the original motor performance to demonstrate the effectiveness of the optimization approach proposed in this paper.