Summary
A new way of describing the density field in density‐based topology optimization is introduced. The new method uses finite elements constructed from Bernstein polynomials rather than the more common Lagrange polynomials. Use of the Bernstein finite elements allows higher‐order elements to be used in the density‐field interpolation without producing unrealistic density values, ie, values lower than zero or higher than one. Results on several test problems indicate that using the higher‐order Bernstein elements produces optimal designs with sharper estimates of the optimal boundary on coarse design meshes. However, higher‐order elements are also required in the structural analysis to prevent the appearance of unrealistic material distributions. The Bernstein element density interpolation can be combined with adaptive mesh refinement to further improve design accuracy even on design domains with complex geometry.