The problem is to find the optimal shape of material that gives minimal compliance under constraints. Based on shape representation techniques, existing methods can be categorized into level set method, density field method, phase field method, etc. Optimization tools include gradient descent, optimality criteria method, method of moving asymptotes and neural networks. I will discuss these methods and present implementation of the phase field method to approach this problem.  Some preliminary results using threshold dynamics have also been obtained.