In this paper, we propose a novel mathematical model for the integrated optimization for production planning, scheduling, and dynamic optimization of continuous manufacturing processes. With a novel approach to estimate the inventory cost, the integrated problem is first formulated as a mixed-integer dynamic optimization (MIDO) problem, which is then reformulated into a large scale mixed-integer nonlinear program (MINLP). By using metamodeling to characterize the detailed linking information between different decision layers, we develop an efficient solution method to decompose the integrated MINLP into a mixed-integer linear program (MILP) for an integrated planning and scheduling problem with flexible recipes, and a set of dynamic optimization problems for generating metamodels. We further apply a bilevel decomposition algorithm to improve the computational efficiency of solving the integrated planning and scheduling problem with flexible recipes. The proposed models and algorithms are demonstrated through case studies of methyl methacrylate (MMA) polymerization processes. The results show that the proposed methods reduce the computational time by more than 2 orders of magnitude compared with the approach of solving the entire integrated optimization problem directly.