Pith. sign in

REVIEW 2 major objections 36 references

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

T0 review · 2 major / 0 minor · reviewed 2026-06-27 · grok-4.3

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Checkerboard patterns in SIMP topology optimization localize to multiaxial stress regions due to penalization and element locking.

desk verdict The paper gives a mechanical story for why checkerboards localize to multiaxial regions under SIMP, but the numerical evidence for isolating stress state from mesh and boundary effects is not shown in enough detail to confirm the claim. read the letter →

arxiv 2606.10509 v1 pith:UGWF5NB2 submitted 2026-06-09 cs.CE

Iulian Paunel , Jonathan Stollberg , Dominik Schillinger This is my paper · ORCID
classification cs.CE
keywords checkerboardpatternsSIMPtopologyoptimizationmultiaxialstressfiniteelementlockingdensitypenalizationnumericalartifactsloadtransferregions
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper establishes that checkerboard patterns systematically emerge in multiaxial load-transfer regions while uniaxial stress regions stay free of them. Continuous intermediate densities would efficiently carry loads in several directions at once, but SIMP penalization suppresses those densities in favor of solid or void states. Linear finite elements then form checkerboard layouts that artificially overestimate stiffness through locking, serving as a discrete substitute. This does not occur along uniaxial paths, where solid struts are naturally preferred without the artifact. The work supplies a mechanical account of both the origin and the spatial localization of the patterns.

What carries the argument

The localization mechanism arising from the combination of multiaxial stress states, SIMP penalization of intermediate densities, and locking-induced artificial stiffness in linear finite elements.

What would settle it

A numerical experiment showing checkerboard patterns throughout a purely uniaxial stress field or their complete absence in a confirmed multiaxial region under identical optimization settings would falsify the localization claim.

Watch

Extended reading notes

Core claim

Checkerboard patterns originate where continuous intermediate densities are mechanically favorable for multiaxial load transfer but are suppressed by SIMP penalization. Linear elements provide an artificially stiff discrete substitute for these penalized regions through their locking behavior. In contrast, uniaxial load paths favor continuous solid struts, making checkerboards mechanically disadvantageous. This supplies a unified interpretation of checkerboarding as the interplay between global stress states, SIMP penalization, and element-level locking.

Load-bearing premise

That the observed checkerboard localization results from the interplay of multiaxial stresses, SIMP penalization, and linear element locking rather than from mesh alignment or other numerical factors.

Editorial extensions

If this is right

  • Checkerboarding is expected at junctions, corners, or bends where stress directions change.
  • Straight axial members or truss-like structures will remain checkerboard-free.
  • Higher-order elements may reduce the artifact but will not remove the underlying mechanical preference for checkerboards in multiaxial zones.
  • The phenomenon scales with the global stress distribution rather than depending solely on local element properties.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Optimization settings such as penalization power may need adjustment depending on whether the problem is dominated by multiaxial or uniaxial load transfer.
  • Mesh refinement alone is unlikely to eliminate the patterns if the mechanical favorability for intermediate densities persists.
  • Alternative density interpolation schemes could be evaluated specifically in multiaxial subdomains to isolate the contribution of penalization.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 0 minor

Summary. The paper claims that checkerboard patterns in density-based topology optimization with the SIMP method and linear finite elements systematically localize to multiaxial load-transfer regions (while uniaxial regions remain free of them) because intermediate densities are mechanically favorable for multiaxial load transfer but suppressed by penalization, with checkerboard layouts serving as an artificially stiff discrete substitute due to linear-element behavior. This is demonstrated via systematic numerical investigations and framed as a unified mechanical account involving global stress states, SIMP penalization, and element-level locking.

Significance. If the central claim holds after verification, the work would supply a useful mechanical account of checkerboard localization that extends prior explanations focused on incompatibility or locking. It could inform targeted regularization or element-selection strategies in topology optimization under complex loading, particularly if the numerical cases cleanly isolate stress-state effects.

major comments (2)
  1. [Abstract] Abstract: The localization claim rests on 'systematic numerical investigations' that are said to isolate global stress multiaxiality, yet the abstract supplies no description of controls that hold mesh alignment, element aspect ratios, boundary-condition orientations, and solver tolerances fixed while varying only local stress state (e.g., rotated loads on identical meshes). This isolation is load-bearing for the asserted causal mechanism.
  2. [Numerical investigations] Numerical investigations (presumably §3–4): Without explicit reporting of exclusion criteria, data sets, and fixed parameters across the load cases, it is impossible to confirm that the observed spatial pattern is produced by the proposed stress-state/SIMP/element interplay rather than by the very mesh or BC confounders flagged in the weakest assumption.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments, which identify opportunities to strengthen the clarity and reproducibility of our numerical evidence. We address each major comment below and will incorporate revisions to improve the manuscript.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The localization claim rests on 'systematic numerical investigations' that are said to isolate global stress multiaxiality, yet the abstract supplies no description of controls that hold mesh alignment, element aspect ratios, boundary-condition orientations, and solver tolerances fixed while varying only local stress state (e.g., rotated loads on identical meshes). This isolation is load-bearing for the asserted causal mechanism.

    Authors: We agree that the abstract would benefit from explicitly summarizing the controls used to isolate stress-state effects. In the revised version we will add a concise clause describing the fixed parameters (mesh alignment, element aspect ratios, boundary-condition orientations, and solver tolerances) and the use of rotated loads on identical meshes to vary only the local stress state. This change directly addresses the concern that the isolation is load-bearing for the claimed mechanism. revision: yes

  2. Referee: [Numerical investigations] Numerical investigations (presumably §3–4): Without explicit reporting of exclusion criteria, data sets, and fixed parameters across the load cases, it is impossible to confirm that the observed spatial pattern is produced by the proposed stress-state/SIMP/element interplay rather than by the very mesh or BC confounders flagged in the weakest assumption.

    Authors: We acknowledge that a more systematic and explicit presentation of the numerical protocol is needed. While §§3–4 describe the load cases, we will add a dedicated summary (table or subsection) that lists the fixed parameters held constant across all cases, the exclusion criteria for load cases, the data sets employed, and the precise manner in which only the stress state is varied. This revision will allow readers to verify that the observed localization arises from the stress-state/SIMP/element interplay rather than mesh or boundary-condition confounders. revision: yes

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: interpretation rests on numerical observations, not self-referential definitions or fitted predictions

full rationale

The paper advances a mechanical interpretation of checkerboard localization based on systematic numerical investigations of SIMP-penalized topology optimization with linear elements. The central claim—that checkerboards emerge preferentially in multiaxial load-transfer regions because intermediate densities are mechanically favorable there but penalized, while uniaxial regions favor solid struts—follows from observed spatial patterns in the computed designs rather than from any equation that reduces to its own inputs by construction. No self-definitional loops, fitted parameters renamed as predictions, or load-bearing self-citations appear in the provided abstract or description. The derivation chain is therefore self-contained against external benchmarks (the optimization runs themselves).

Assumptions & free parameters 0 free parameters · 2 assumptions · 0 invented entities

The claim depends on domain assumptions about mechanical favorability of intermediate densities under multiaxial versus uniaxial stress and on the locking behavior of linear elements; no free parameters or invented entities are introduced in the abstract.

assumptions (2)
  • domain assumption Continuous intermediate densities are mechanically favorable for multiaxial load transfer but suppressed by SIMP penalization.
    Invoked directly in the abstract to explain why checkerboards serve as substitutes only in multiaxial regions.
  • domain assumption Linear finite elements exhibit characteristic behavior that makes checkerboard layouts artificially stiff.
    Cited as the element-level mechanism enabling the discrete substitute in penalized regions.

how reviews work

0 comments
Cite this review

Pith. "Pith review of On the Localization of Checkerboarding in Multiaxial Stress Regions under SIMP Penalization." pith.science (2026). https://pith.science/paper/UGWF5NB2

@misc{pith2026260610509,
  author       = {Pith},
  title        = {Pith review of: On the Localization of Checkerboarding in Multiaxial Stress Regions under SIMP Penalization},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/UGWF5NB2}},
  note         = {Machine review of arXiv:2606.10509}
}
read the original 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.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

36 extracted references · 27 canonical work pages

  1. [1]

    Ambrosio, L., & Buttazzo, G. (1993). An optimal design problem with perimeter penalization. Calc. Var. Partial Differ. Equ.,1(1), 55–69, https://doi.org/10.1007/BF02163264

    work page doi:10.1007/bf02163264 1993

  2. [2]

    Kara, L. (2018). 3d topology optimization using convolutional neural networks.arXiv preprint, , https://doi.org/10.48550/arXiv .1808.07440

    work page Pith review arXiv doi:10.48550/arxiv 2018

  3. [3]

    Structural Optimization , volume = 1, number = 4, pages =

    Bathe, K.-J. (2014).Finite element procedures (2nd ed.). Watertown, MA: Klaus-J¨ urgen Bathe. Bendsøe, M.P. (1989). Optimal shape design as a material distribution problem.Struct. Mul- tidiscip. Optim.,1(4), 193–202, https:// doi.org/10.1007/BF01650949 Bendsøe, M.P., & Kikuchi, N. (1988). Generat- ing optimal topologies in structural design using a homoge...

    work page doi:10.1007/bf01650949 2014

  4. [4]

    (2013).Mixed finite element methods and applications(1st ed.)

    Boffi, D., Brezzi, F., Fortin, M. (2013).Mixed finite element methods and applications(1st ed.). Berlin, Heidelberg: Springer

    2013

  5. [5]

    Borrvall, T. (2001). Topology optimization of elas- tic continua using restriction.Arch. Com- put. Methods Eng.,8(4), 351–385, https:// doi.org/10.1007/BF02743737

    work page doi:10.1007/bf02743737 2001

  6. [6]

    Bourdin, B. (2001). Filters in topology opti- mization.Int. J. Numer. Methods Eng., 50(9), 2143–2158, https://doi.org/10.1002/ nme.1160

    2001

  7. [7]

    Brackett, D., Ashcroft, I., Hague, R. (2011). Topology optimization for additive manufac- turing.2011 International Solid Freeform Fabrication Symposium, , https://doi.org/ 10.26153/tsw/15300

    work page doi:10.26153/tsw/15300 2011

  8. [8]

    (1991).Mixed and hybrid finite element methods(1st ed.)

    Brezzi, F., & Fortin, M. (1991).Mixed and hybrid finite element methods(1st ed.). New York, NY: Springer

    1991

Show all 36 references
  1. [9]

    Bruggi, M. (2009). Generating strut-and-tie pat- terns for reinforced concrete structures using topology optimization.Comput. Struct., 87(23-24), 1483–1495, https://doi.org/10 .1016/j.compstruc.2009.06.003

    2009

  2. [11]

    Diaz, A., & Sigmund, O. (1995). Checkerboard patterns in layout optimization.Struct. Multidiscip. Optim.,10(1), 40–45, https:// doi.org/10.1007/BF01743693

    1995 doi

  3. [12]

    (2003).Finite ele- ment methods for flow problems(1st ed.)

    Donea, J., & Huerta, A. (2003).Finite ele- ment methods for flow problems(1st ed.)

    2003

  4. [13]

    Gangwar, T., & Schillinger, D. (2021). Con- current material and structure optimization of multiphase hierarchical systems within a continuum micromechanics framework. Struct. Multidiscip. Optim.,64(3), 1175– 1197, https://doi.org/10.1007/s00158-021 -02907-1

    2021 doi

  5. [14]

    Guest, J.K., Pr´ evost, J.H., Belytschko, T. (2004). Achieving minimum length scale in topology optimization using nodal design variables and projection functions.Int. J. Numer. Methods Eng.,61(2), 238–254, https:// doi.org/10.1002/nme.1064

    2004 doi

  6. [15]

    Haber, R.B., Jog, C.S., Bendsøe, M.P. (1996). A new approach to variable-topology shape design using a constraint on perimeter. Struct. Multidiscip. Optim.,11(1), 1–12, https://doi.org/10.1007/BF01279647 H¨ agg, L., & Wadbro, E. (2017). Nonlinear filters in topology optimizati...

    1996 doi

  7. [16]

    (1987).The finite element method: linear static and dynamic finite ele- ment analysis(1st ed.)

    Hughes, T.J.R. (1987).The finite element method: linear static and dynamic finite ele- ment analysis(1st ed.). Englewood Cliffs, NJ: Prentice-Hall

    1987

  8. [17]

    Park, C., Kim, M.-N. (2003). Checkerboard- free topology optimization using non- conforming finite elements.Int. J. Numer. Methods Eng.,57(12), 1717–1735, https:// doi.org/10.1002/nme.738

    2003 doi

  9. [18]

    Jog, C.S., & Haber, R.B. (1996). Stability of finite element models for distributed-parameter 21 optimization and topology design.Com- put. Methods Appl. Mech. Eng.,130(3- 4), 203–226, https://doi.org/10.1016/0045 -7825(95)00928-0

    1996 doi

  10. [19]

    Johnson, C., & Pitk¨ aranta, J. (1982). Analy- sis of some mixed finite element methods related to reduced integration.Math. Com- put.,38(158), 375–400, https://doi.org/ 10.2307/2007276

    1982 doi

  11. [20]

    Nomura, T., Kondoh, T., Nishiwaki, S. (2011). Heaviside projection based topology optimization by a pde-filtered scalar func- tion.Struct. Multidiscip. Optim.,44(1), 19–24, https://doi.org/10.1007/s00158-010 -0562

    2011 doi

  12. [21]

    Lazarov, B.S., & Sigmund, O. (2011). Filters in topology optimization based on helmholtz- type differential equations.Int. J. Numer. Methods Eng.,86(6), 765–781, https://doi .org/10.1002/nme.3072

    2011 doi

  13. [23]

    Petersson, J. (1999). Some convergence results in perimeter-controlled topology optimization. Comput. Methods Appl. Mech. Eng.,171(1- 2), 123–140, https://doi.org/10.1016/S0045 -7825(98)00248-5

    1999 doi

  14. [24]

    Petersson, J., & Sigmund, O. (1998). Slope constrained topology optimization. Int. J. Numer. Methods Eng.,41(8), 1417–1434, https://doi.org/10.1002/ (SICI)1097-0207(19980430)41:8⟨1417:: AID-NME344⟩3.0.CO;2-N

    1998

  15. [25]

    Poulsen, T.A. (2002). A simple scheme to pre- vent checkerboard patterns and one-node connected hinges in topology optimization. Struct. Multidiscip. Optim.,24(5), 396– 399, https://doi.org/10.1007/s00158-002 -0251-x

    2002 doi

  16. [26]

    Sani, R.L., Gresho, P.M., Lee, R.L., Griffiths, D.F. (1981). The cause and cure (?) of the spurious pressures generated by certain fem solutions of the incompressible navier-stokes equations: Part 1.Int. J. Numer. Meth- ods Fluids,1(1), 17–43, https://doi.org/ 10.1002/fld.1650010104

    1981 doi

  17. [27]

    (1994).Design of material structures using topology optimization(Unpublished doctoral dissertation)

    Sigmund, O. (1994).Design of material structures using topology optimization(Unpublished doctoral dissertation). Technical University of Denmark, Lyngby

    1994

  18. [28]

    Sigmund, O. (2001). A 99 line topology optimiza- tion code written in matlab.Struct. Mul- tidiscip. Optim.,21(2), 120–127, https:// doi.org/10.1007/s001580050176

    2001 doi

  19. [30]

    Sigmund, O., & Maute, K. (2013). Topology opti- mization approaches: A comparative review. Struct. Multidiscip. Optim.,48(6), 1031– 1055, https://doi.org/10.1007/s00158-013 -0978-6

    2013 doi

  20. [31]

    Sigmund, O., & Petersson, J. (1998). Numeri- cal instabilities in topology optimization: a survey on procedures dealing with checker- boards, mesh-dependencies and local min- ima.Struct. Multidiscip. Optim.,16(1), 68– 75, https://doi.org/10.1007/BF01214002

    1998 doi

  21. [32]

    Sosnovik, I., & Oseledets, I. (2019). Neural net- works for topology optimization.Russian J. 22 Numer. Anal. Math. Modelling,34(4), 215– 223, https://doi.org/10.1515/rnam-2019 -0018

    2019 doi

  22. [33]

    Schillinger, D. (2025). Multiscale topol- ogy optimization of functionally graded lattice structures based on physics- augmented neural network material models. Comput. Methods Appl. Mech. Eng., 438, 117808, https://doi.org/10.1016/ j.cma.2025.117808

    2025

  23. [34]

    Stolpe, M., & Svanberg, K. (2001a). An alterna- tive interpolation scheme for minimum com- pliance topology optimization.Struct. Mul- tidiscip. Optim.,22(2), 116–124, https:// doi.org/10.1007/s001580100129

    doi

  24. [35]

    Stolpe, M., & Svanberg, K. (2001b). On the trajectories of penalization methods for topology optimization.Struct. Multidiscip. Optim.,21(2), 128–139, https://doi.org/ 10.1007/s001580050177

    doi

  25. [36]

    Talischi, C., Paulino, G.H., Le, C.H. (2009). Hon- eycomb wachspress finite elements for struc- tural topology optimization.Struct. Multi- discip. Optim.,37(6), 569–583, https:// doi.org/10.1007/s00158-008-0261-4

    2009 doi

  26. [37]

    Wadbro, E., & H¨ agg, L. (2015). On quasi- arithmetic mean based filters and their fast evaluation for large-scale topology optimiza- tion.Struct. Multidiscip. Optim.,52(5), 879–888, https://doi.org/10.1007/s00158 -015-1273-5

    2015 doi

  27. [38]

    Xu, S., Cai, Y., Cheng, G. (2010). Volume preserv- ing nonlinear density filter based on heavi- side functions.Struct. Multidiscip. Optim., 41(4), 495–505, https://doi.org/10.1007/ s00158-009-0452-7

    2010

  28. [39]

    Zhou, M., Shyy, Y., Thomas, H. (2001). Checker- board and minimum member size control in topology optimization.Struct. Multidiscip. Optim.,21(2), 152–158, https://doi.org/ 10.1007/s001580050179 23

    2001 doi

Pith tools

Reviewed June 27, 2026 · model on record in the stance chip above.