Abstract: The coefficients of strain energy density are incorporated in every finite element, which makes the expression of the globalized strain energy density constraint in the ICM method more reliable by enforcing the local stress constraint with high accuracy. Furthermore, the original coefficients of strain energy density are relaxed in order to accelerate the optimization iteration convergence. A power function formed as the ratio of the maximal structural von Mises stress and the allowable stress is applied to modify the constraint value of the globalized constraint to control the structural maximal stress or to accelerate the optimization iteration convergence. For the topology configuration to evolve in a reasonable direction, the continuous topological design variables are converted to discrete ones by using of a dynamic filter criterion computed automatically. Some examples show that the method proposed in this paper is reasonable and efficient. The optimal topology configurations are similar to the Michell truss, which is known as the theoretical optimal solution.