Numerical instabilities such as checkerboard patterns, gray-scale elements and mesh dependency always occur in topology optimization of continuum structures.Many optimization algorithms, including the method of using higher -order finite elements, perimeter control method and the filtering method, have been proposed and tested in optimization problems.