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

Dominik Schillinger; Iulian Paunel; Jonathan Stollberg

arxiv: 2606.10509 · v1 · pith:UGWF5NB2new · submitted 2026-06-09 · 💻 cs.CE

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

Iulian Paunel , Jonathan Stollberg , Dominik Schillinger This is my paper

Pith reviewed 2026-06-27 11:34 UTC · model grok-4.3

classification 💻 cs.CE

keywords checkerboard patternsSIMP topology optimizationmultiaxial stressfinite element lockingdensity penalizationnumerical artifactsload transfer regions

0 comments

The pith

Checkerboard patterns in SIMP topology optimization localize to multiaxial stress regions due to penalization and element locking.

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

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.

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.

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.

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.

Where Pith is reading between the lines

These are 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.

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.

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.

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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

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

The new piece is the claim that checkerboard patterns appear specifically in multiaxial load-transfer zones while uniaxial paths stay free of them. The account is that intermediate densities would help with multiaxial transfer but get penalized by SIMP, so linear elements fall back on checkerboards as a stiff substitute; uniaxial cases favor solid struts instead. This ties the known artificial stiffness to a location pattern that earlier incompatibility or locking explanations left open.

The paper reports systematic numerical investigations to back this up. If those runs actually held mesh alignment, element aspect ratios, and boundary orientations fixed while only changing the local stress multiaxiality, the localization result would be useful. The abstract does not describe those controls, so the central interpretation still rests on unshown steps.

The reasoning itself is direct and avoids fitting parameters or self-reference. It connects global stress states, penalization, and element behavior without contradiction. No new entities are introduced.

This is for people who run or debug density-based topology optimization in structural mechanics. The work is coherent enough on its own terms to deserve referee time so the experiments can be checked for the missing isolation details. A review would likely focus on whether the pattern survives rotated loads on the same mesh or orthogonal cases on identical discretizations.

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)

[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.
[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

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
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

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).

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 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.

axioms (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.

pith-pipeline@v0.9.1-grok · 5717 in / 1464 out tokens · 31039 ms · 2026-06-27T11:34:21.709487+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

36 extracted references · 26 canonical work pages · 1 internal anchor

[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]

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

work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv 2018
[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]

(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]

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]

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]

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]

(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
[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
[11]

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

work page doi:10.1007/bf01743693 1995
[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
[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

work page doi:10.1007/s00158-021 2021
[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

work page doi:10.1002/nme.1064 2004
[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 optimization: existence of solutions and efficient implementation for minimum compliance p...

work page doi:10.1007/bf01279647 1996
[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
[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

work page doi:10.1002/nme.738 2003
[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

work page doi:10.1016/0045 1996
[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

work page doi:10.2307/2007276 1982
[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

work page doi:10.1007/s00158-010 2011
[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

work page doi:10.1002/nme.3072 2011
[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

work page doi:10.1016/s0045 1999
[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
[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

work page doi:10.1007/s00158-002 2002
[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

work page doi:10.1002/fld.1650010104 1981
[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
[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

work page doi:10.1007/s001580050176 2001
[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

work page doi:10.1007/s00158-013 2013
[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

work page doi:10.1007/bf01214002 1998
[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

work page doi:10.1515/rnam-2019 2019
[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

arXiv 2025
[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

work page doi:10.1007/s001580100129
[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

work page doi:10.1007/s001580050177
[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

work page doi:10.1007/s00158-008-0261-4 2009
[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

work page doi:10.1007/s00158 2015
[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
[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

work page doi:10.1007/s001580050179 2001

[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 internal anchor Pith review Pith/arXiv 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

[9] [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

[10] [11]

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

work page doi:10.1007/bf01743693 1995

[11] [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

[12] [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

work page doi:10.1007/s00158-021 2021

[13] [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

work page doi:10.1002/nme.1064 2004

[14] [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 optimization: existence of solutions and efficient implementation for minimum compliance p...

work page doi:10.1007/bf01279647 1996

[15] [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

[16] [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

work page doi:10.1002/nme.738 2003

[17] [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

work page doi:10.1016/0045 1996

[18] [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

work page doi:10.2307/2007276 1982

[19] [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

work page doi:10.1007/s00158-010 2011

[20] [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

work page doi:10.1002/nme.3072 2011

[21] [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

work page doi:10.1016/s0045 1999

[22] [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

[23] [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

work page doi:10.1007/s00158-002 2002

[24] [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

work page doi:10.1002/fld.1650010104 1981

[25] [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

[26] [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

work page doi:10.1007/s001580050176 2001

[27] [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

work page doi:10.1007/s00158-013 2013

[28] [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

work page doi:10.1007/bf01214002 1998

[29] [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

work page doi:10.1515/rnam-2019 2019

[30] [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

arXiv 2025

[31] [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

work page doi:10.1007/s001580100129

[32] [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

work page doi:10.1007/s001580050177

[33] [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

work page doi:10.1007/s00158-008-0261-4 2009

[34] [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

work page doi:10.1007/s00158 2015

[35] [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

[36] [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

work page doi:10.1007/s001580050179 2001