In this study, an optimal shape design was performed on a fillet model which is subject to surface traction through minimizing the maximum stress of the fillet. A 2-dimensional quadratic isoparametirc element with 8 nodes was used in stress analysis for finite element method, and Hooke-Jeeves direct search algorithm was adopted for optimi- zation. From the resulting optimal shape, it was found that the maxium von Mises stress on the boundary of fillet was reduced by 36%, compared to other paper in which the cross sectional area of fillet was minimized. In conclusion, a real optimal fillet shape could be obtained in the viewpoint of yielding, and more pratical fillet design could be accomplished.