In this study, a viscoplastic FEM simulation of the stretch blow molding process is presented. The material＇s behavior at elevated temperature is assumed to be sh ain-rate-dependent, and a non-Newtonian creeping material model is employed to specify the strain-rate-sensitive characteristic of the material. A degenerated quadrilateral shell element, which enables proper evaluation of the bending effect, is employed. The mold surface is described by triangular surfaces, and the contact phenomenon is modeled as a node-to-mesh contact. In using the code developed, the stretch blow molding process is taken as an example and discussed qualitatively. First, the process under three boundary conditions is simulated and examined, and then the condition leading to the wrinkle phenomenon, which is encountered in the simulation, is discussed.