Abstract: The topology optimization model of continuum structures with stress’ and displacement’s constraints is formulated, applying ICM (Independent Continuous and Mapping) method of high accuracy approximation, in which the exponential type function is adopted as fast filter function. Stress constraints are globalized by use of structural distortional strain energy density. Displacement constraint is explicited by use of unit virtual load method. The constraint limit values are normalized to solve the mismatch of the order of two kinds of constraints. The parameter selection of fast filter functions of exponential type with different properties is given. The examples of single loading case and multiple loading cases have shown that the ICM method of high accuracy approximation is available and efficient to solve the topology optimization of continuum structures with multiple constraint types.