We have proposed a new design method to optimize a cooling pipe shape in the three-dimensional space for plastic injection molding die by applying the response surface approach and shape optimization technique. To evaluate the optimality of cooling pipe arrangement and shape, cooling rate and temperature distribution in the die, which are primal factors to cause residual stresses producing some defects in the molded plastic, are used as objective functions. Then the cooling pipe shape and position to minimize the two objective functions simultaneously by using optimization techniques has been obtained. The proposed method has been applied to some cooling pipe shape design of 3-dimensional problem to discuss the availability off the method.