This study presents an optimization strategy developed for the stretch-blow molding process. The method is based on a coupling between the Nelder-Mead optimization algorithm and finite element (FE) simulations of the forming process developed using ABAQUS 1 1 1 . FE simulations were validated using in situ tests and measurements performed on 18.5 g-50 cl polyethylene terephthalate bottles. To achieve that, the boundary conditions were carefully measured for both the infrared heating and the blowing stages. The temperature distribution of the perform was predicted using a 3D finite-volume software, and then applied as an initial condition into FE simulations. In addition, a thermodynamic model was used to predict the air pressure applied inside the preform, taking into account the relationship between the internal air pressure and the enclosed volume of the preform, i.e., the fluid-structure interaction. It was shown that the model adequately predicts both the blowing kinematics and the thickness distributions of the bottle. In the second step, this model was combined to an optimization loop to automatically compute the best preform temperature distribution, providing a uniform thickness for the bottle.