The objective of this study is to develop an integrated computer-aided engineering (CAE) optimization system that can quickly and intelligently determine the optimal process conditions for injection molding. This study employs support vector regression (SVR) to establish the surrogate model based on executions of three-dimensional (3D) simulation for a selected dataset using the latin hypercube sampling (LHS) technique. Once the surrogate model can satisfactorily capture the characteristics of simulations with much less computing resources, a hybrid optimization genetic algorithm (GA) or a multi-objective optimization GA is then used to evaluate the surrogate model to search the global optimal solutions for the single or multiple objectives, respectively. The performance and capabilities of other surrogate modeling approaches, such as polynomial regression (PR) and artificial neural network (ANN), are also investigated in terms of accuracy, robustness, efficiency, and requirements for training samples. Experimental validations and applications of this work for process optimization of a special box mold and a precision optical lens are presented.