On the Localization of Checkerboarding in Multiaxial Stress Regions under SIMP Penalization

2026-06-09

AI Summary

The spatial localization mechanism underlying checkerboard artifacts in density-based topology optimization has long remained unclear. This study, through systematic numerical experiments combining the SIMP density method with linear finite element analysis, reveals for the first time that checkerboard patterns stably emerge only in multiaxial stress regions and vanish in uniaxial regions. The work demonstrates that this phenomenon arises because, after SIMP penalization suppresses intermediate densities, linear elements artificially fulfill multiaxial load-path requirements by adopting a checkerboard arrangement, thereby establishing a localized stiffness substitution mechanism. These findings provide a unified explanation for both the origin and spatial localization of checkerboard artifacts, clarifying the coupled interplay among global stress states, SIMP penalization, and element locking.

Abstract

Checkerboard patterns are a well-known numerical artifact in density-based topology optimization using the Solid Isotropic Material with Penalization (SIMP) method and linear finite elements. Existing explanations based on mixed-field incompatibility or locking-induced stiffness overestimation explain the artificial stiffness of checkerboard layouts but do not clarify their characteristic spatial localization. In this work, we show that checkerboard patterns systematically emerge in multiaxial load-transfer regions, whereas predominantly uniaxial stress regions remain checkerboard-free. Through systematic numerical investigations, we demonstrate that checkerboarding originates where continuous intermediate densities are mechanically favorable for multiaxial load transfer but are suppressed by SIMP penalization. Due to the characteristic behavior of linear elements, checkerboard layouts provide an artificially stiff discrete substitute for these penalized intermediate-density regions. In contrast, uniaxial load paths naturally favor continuous solid struts, rendering checkerboards mechanically disadvantageous. Our findings provide a unified mechanical interpretation of checkerboarding as the interplay between global stress states, SIMP penalization, and element-level locking, thereby explaining both its origin and the spatial localization.

Problem

• checkerboarding
• topology optimization
• SIMP
• multiaxial stress
• numerical artifact

Innovation

• checkerboarding
• SIMP penalization
• multiaxial stress
• topology optimization
• finite element locking