Efficiently managing multiple process parameters is critical for achieving optimal performance in additive manufacturing. This study investigates the relationship between eight key parameters in fused deposition modeling (FDM) and their impact on responses like average surface roughness (Ra), tensile strength (TS), and flexural strength (FS) of carbon fiber-reinforced polyamide 12 (PA 12-CF) material. The study integrates response surface methodology (RSM), grey relational analysis (GRA), and grey wolf optimization (GWO) to achieve this goal. A total of 51 experiments were planned using a definitive screening design (DSD) based on response RSM. The printing process parameters, including layer thickness, infill density, and build orientation, significantly affect Ra, TS, and FS. GRA combines responses into a single measure, grey relational grade (GRG), and a regression model is developed. GWO is then employed to optimize GRG across parameters. Comparison with GRA-optimized parameters demonstrates GWO’s ability to discover refined solutions, reducing average surface roughness to 4.63 μm and increasing tensile strength and flexural strength to 88.5 MPa and 103.12 MPa, respectively. Practical implications highlight the significance of GWO in industrial settings, where optimized parameters lead to reduced costs and improved product quality. This integrated approach offers a systematic methodology for optimizing FDM processes, ensuring robustness and efficiency in additive manufacturing applications.