The SIMP (solid isotropic material with penalization) approach is perhaps the most popular density variable relaxation method in topology optimization. This method has been very successful in many applications, but the optimization solution convergence can be improved when new variables, not the direct density variables, are used as the design variables. In this work, we newly propose S-shape functions mapping the original density variables nonlinearly to new design variables. The main role of S-shape function is to push intermediate densities to either lower or upper bounds. In particular, this method works well with nonlinear mathematical programming methods. A method of feasible directions is chosen as a nonlinear mathematical programming method in order to show the effects of the S-shape scaling function on the solution convergence.