Electronics Guide

Topology Optimization and Lattice Design

Electronic products are rarely bare circuit boards. They ship inside chassis, brackets, cold plates, heat sinks, stiffeners, and shielding cans, and that structure often dominates mass, volume, and thermal performance. Topology optimization and lattice design are computational methods for shaping this structure. Instead of drawing a part and then analyzing it, the engineer defines a design space, the loads and temperatures it must carry, and a performance objective, and lets an optimizer decide where material belongs.

The two methods are complementary. Topology optimization determines the coarse distribution of material across a region: where the load paths run and where material can be removed entirely. Lattice design fills the regions that remain with engineered cellular architecture, trading a controlled amount of stiffness for large reductions in mass, or adding surface area where heat must be moved into a fluid. Both belong in the design-for-excellence toolkit because both make lifecycle trade-offs explicit at the moment the geometry is created, rather than after tooling exists.

Topology optimization reached industrial practice well before three-dimensional printing was common: automotive and aerospace engineers used it to shape cast, forged, and machined structures, with draft, parting-line, and symmetry constraints written into the problem statement. What additive manufacturing changed is the range of answers that can actually be built. The organic, load-path-following geometry an unconstrained optimizer produces is difficult or impossible to machine or mold, and lattices at useful scale are effectively an additive-only proposition. That coupling is also the discipline's central caution: an optimization posed without process limits returns a shape no available process can build, at a cost no program can absorb. The useful skill is not running the solver but posing the problem so that its answer is manufacturable, inspectable, and cheaper in total than the part it replaces.

What the Method Actually Computes

A topology optimization problem has four ingredients. The design domain is the volume the part may occupy, typically the largest envelope the surrounding assembly permits. Non-design regions are the parts of that volume the optimizer may not touch: bolt bosses, connector cutouts, bearing seats, mating flanges, and keep-out zones for cables and airflow. Boundary conditions define how the part is held and what it carries, including static loads, accelerations, thermal loads, and vibration inputs. Finally, an objective and a set of constraints state what "better" means.

The classic formulation minimizes compliance, the product of the applied loads and the displacements they produce, which for a linear elastic structure equals twice the stored strain energy. The minimization is subject to a constraint on the fraction of the domain that may remain solid. Minimizing compliance is equivalent to maximizing stiffness for the given loads. A volume fraction of 0.3 asks the optimizer for the stiffest arrangement that uses at most thirty percent of the available material. Other formulations invert the statement, minimizing mass subject to limits on displacement, stress, or natural frequency, which is usually closer to how a hardware requirement is actually written.

The domain is discretized into finite elements, and each element receives a design variable representing how much material it holds. The optimizer evaluates the structure, computes the sensitivity of the objective to each design variable, nudges material toward the elements that contribute most, and repeats. The loop stops when the design variables cease to change appreciably, when the objective flattens, or when a prescribed iteration limit is reached. Because sensitivities are obtained by the adjoint method, the cost of a full gradient is about one extra solve per response function no matter how many design variables exist, and the minimum-compliance objective is self-adjoint, so even that solve is unnecessary. That property is what makes problems with millions of design variables tractable.

Density Methods

The modern field dates to a 1988 paper by Martin Bendsøe and Noboru Kikuchi, which posed structural design as a material distribution problem and solved it by homogenization, treating each element as a microstructure whose hole size and orientation are the design variables. Homogenization in that original form is little used for design today, but its framing underlies everything built since, and homogenization itself returns below as the practical way to represent a lattice without meshing every strut.

The dominant approach now is the density method, usually implemented as Solid Isotropic Material with Penalization (SIMP), which Bendsøe proposed in 1989 as a simpler power-law interpolation and which Ole Sigmund, George Rozvany, and others developed extensively. Each element carries a pseudo-density between zero and one, and its stiffness is scaled by that density raised to a penalization exponent, commonly three. The exponent makes intermediate densities structurally inefficient: a half-dense element delivers about an eighth of the stiffness for half the material, so the optimizer is pushed toward a design that is nearly all solid or nearly all void. A small lower bound on stiffness keeps the finite element system from becoming singular where density approaches zero.

Density methods are robust, easy to implement, and available in the major commercial structural solvers. Their compact expression in teaching code is well known: Sigmund's ninety-nine-line MATLAB implementation, published in 2001, and the faster eighty-eight-line rewrite by Andreassen and colleagues in 2011 remain the standard introductions to the field.

Level-Set and Evolutionary Methods

Level-set methods represent the structure implicitly, as the region where a scalar field is positive, and evolve the boundary rather than element densities. The representation comes from Stanley Osher and James Sethian's 1988 work on propagating interfaces; Sethian and Wiegmann applied it to structural topology optimization in 2000, and Grégoire Allaire and coworkers developed the shape-derivative formulation in common use since. The result is a crisp boundary at every iteration with no gray transition zone, which suits problems where surface effects matter, such as stress concentration or convective heat transfer. The trade-off is that a conventional level-set method cannot nucleate new holes inside the material, so the initial design influences the outcome more strongly than in density methods. Hybrid schemes add the topological derivative, which measures the effect of introducing an infinitesimal hole at a point, to restore that ability.

Evolutionary structural optimization and its bidirectional variant, developed by Yi Min Xie and Grant Steven beginning in 1993, take a simpler approach: remove elements with low stress or low strain energy, and in the bidirectional form add elements back where they would help most. These methods are intuitive and inexpensive, but they lack the convergence guarantees of gradient-based formulations, and published counterexamples showed that a naive removal criterion can move a design away from the optimum rather than toward it. Later bidirectional formulations answered that criticism with sensitivity numbers, filtering, and averaging across iterations.

Generative Design Is Not the Same Thing

Commercial tools often market topology optimization under the label "generative design." The terms are not synonymous. Topology optimization solves one well-posed optimization problem and returns one converged answer. Generative design workflows wrap that solver in an outer loop that sweeps materials, manufacturing processes, and load cases, then presents a catalog of candidate designs for human selection. The underlying mathematics is the same; the difference is how many problems are posed and who chooses among the results.

Numerical Behavior and the Filters That Tame It

A naive implementation produces artifacts that look like real structure but are not. Sigmund and Petersson's 1998 survey of the subject named three failure modes that still recur, and every production tool now includes countermeasures.

Checkerboarding is the appearance of alternating solid and void elements in a grid pattern. It is a numerical artifact: low-order finite elements overestimate the stiffness of that pattern, so the optimizer finds it attractive. Mesh dependence means the optimizer returns more and thinner members as the mesh is refined, because a finer mesh permits finer structure. Without a length scale in the formulation there is no converged answer, only an answer that tracks the discretization. Local minima arise because the problem is non-convex; different starting points and different continuation schedules for the penalization exponent lead to different designs.

Filtering addresses the first two. Sensitivity filtering averages the computed sensitivities over a neighborhood of a chosen radius; density filtering applies the same averaging to the densities themselves. Either imposes a minimum length scale and suppresses checkerboards, at the cost of a band of intermediate density around every boundary. Heaviside projection schemes then push those intermediate values back toward zero or one so the final geometry is crisp. In practice the filter radius is the single most consequential setting an engineer chooses: it should be tied to the manufacturing process resolution, not to the mesh.

Stress-constrained problems deserve separate mention because they are substantially harder than compliance problems. Stress is a local quantity, so a strict formulation adds one constraint per element; stress is also singular in the low-density limit, producing an ill-behaved design space. Practical solvers use relaxation and constraint aggregation to make the problem tractable, and the results warrant more verification than a compliance-driven design. An optimizer that minimized compliance has said nothing about whether the part meets its stress allowable.

Lattice Structures

Where topology optimization decides whether material is present, lattice design decides what that material looks like at a finer scale. A lattice replaces solid material with a repeating cellular architecture, and its effective properties are set by the geometry of the unit cell and by relative density, the fraction of the cell volume that is solid.

Unit Cell Families

Strut-based lattices are built from beams meeting at nodes. The families are usually named after crystal structures: body-centered cubic, face-centered cubic, octet-truss, and their many reinforced variants. Triply periodic minimal surfaces (TPMS), including the gyroid, Schwarz primitive, and diamond surfaces, are defined by implicit equations rather than beams. They are smooth, self-supporting over a wide range of orientations, and free of the sharp nodal junctions where strut lattices concentrate stress. Each TPMS cell can be built in two ways: as a sheet, thickening the surface itself, or as a skeletal solid filling one of the two volumes the surface separates. The sheet form is generally reported as the stiffer of the two at equal relative density, while the skeletal form leaves a single connected void that is easier to flush. Stochastic foams mimic natural cellular solids and are easier to generate but harder to characterize and reproduce.

The most useful distinction is mechanical rather than geometric. In the framework developed by Lorna Gibson and Michael Ashby, and sharpened by Maxwell's counting rule for pin-jointed frames as Deshpande, Ashby, and Fleck applied it to cellular solids, a lattice is either stretch-dominated or bending-dominated. Maxwell's rule compares the number of struts with the number of degrees of freedom at the nodes: a cell with enough struts to be statically determinate carries load axially, while one without them relies on strut bending. In a stretch-dominated lattice such as the octet-truss, effective stiffness falls roughly in proportion to relative density. In a bending-dominated lattice such as a simple body-centered cubic cell, stiffness falls approximately with the square of relative density. The consequence compounds as the lattice gets lighter: halving relative density costs a stretch-dominated cell about half its stiffness but costs a bending-dominated cell about three quarters of its own, so the gap widens with every step toward lower density. Stretch-dominated cells are therefore the choice for stiffness-critical structure, while bending-dominated cells, which collapse at a long, flat plateau stress, are the choice for energy absorption in shock mounts and drop protection. Absolute stiffness figures for a given cell depend on the alloy, the strut aspect ratio, the node geometry, and the as-built accuracy of the print, so published coefficients should be treated as a starting point and confirmed by test.

Characteristic behavior of common lattice families
Family Dominant mechanism Approximate stiffness scaling Typical use in electronics hardware
Octet-truss and reinforced cubic Stretch Near-linear in relative density Chassis and bracket cores where stiffness per unit mass governs
Body-centered cubic and similar Bending Near-quadratic in relative density Shock and vibration isolation, crushable standoffs
Gyroid and other TPMS Mixed, shell-like Between the two, geometry dependent Cold plates and heat exchangers needing high surface area per volume
Stochastic foam Bending Near-quadratic, with wide scatter Damping and acoustic treatment where tolerance is loose

Grading and Multiscale Design

A uniform lattice wastes material because load is never uniform. Functionally graded lattices vary relative density, cell size, or cell type across the part, placing dense material along load paths and sparse material elsewhere. The usual workflow runs a topology optimization first, then maps the resulting density field onto a lattice grading rule, so that the optimizer's continuous answer becomes a physically realizable gradient of strut diameters.

Simulating every strut in a graded lattice is prohibitive for a full assembly. Homogenization solves this by characterizing the unit cell once, either analytically or by a detailed unit cell simulation, and representing large regions with an equivalent anisotropic continuum. The system-level model then runs at reasonable cost, and detailed models are reserved for critical regions identified by the coarse analysis. Because effective properties depend on the cell family and its orientation, the homogenized material is generally anisotropic, and treating it as isotropic is a common and costly modeling error.

Representing lattices in conventional boundary-representation CAD does not scale, because a modest part may contain hundreds of thousands of struts. Implicit modeling, which stores geometry as signed distance or field functions and evaluates surfaces only when needed, has become the standard approach in dedicated lattice tools, and it makes operations such as thickening, blending, and Boolean trimming inexpensive at any lattice density.

Applications in Electronic Hardware

Thermal Structures

Heat transfer problems admit the same treatment as structural ones. Minimizing an integrated measure of temperature in a conduction problem produces branching, tree-like material distributions that resemble the vascular structures found in biology, because both solve the problem of collecting a distributed source into a concentrated sink. Optimizing a convective heat sink or cold plate requires coupling the structural design variables to a fluid solve, most often by penalizing flow in solid regions with a Brinkman friction term, an approach established by Thomas Borrvall and Joakim Petersson in 2003. The result is a set of channels shaped by the actual flow field rather than by the extrusion constraints of a die.

The practical payoff in electronics is highest where power density is high and pumping power or mass is expensive: additively manufactured aluminum cold plates for traction inverters and power conversion, gyroid-cored heat exchangers where surface area per unit volume is the binding constraint, and conformal cooling that follows the contour of a device rather than a flat mounting plane. These designs must still be evaluated against conventional alternatives, since a skived or bonded-fin heat sink produced in volume is far cheaper than a printed one and frequently sufficient.

Structural and Dynamic Requirements

Avionics, launch, and vehicle hardware is usually governed by stiffness and by natural frequency rather than by static strength. A common requirement sets a minimum first natural frequency so that the assembly does not couple with the input vibration spectrum, and frequency maximization is a standard topology optimization objective. Optimized brackets, gimbal mounts, antenna backing structures, and instrument housings are the usual targets, and case studies report substantial mass reduction against machined baselines at equal or higher fundamental frequency. Individual percentages deserve scrutiny: most are published by the software vendor or the service bureau that produced the part, and the comparison is honest only when the baseline is a competently designed conventional part rather than the first one drawn.

Board-level structure benefits as well. Stiffeners and support ribs sized by optimization control printed circuit board deflection under shock and reduce solder joint strain, a direct contribution to reliability. Lattice-cored standoffs and mounts provide damping and controlled crush in drop events. In each case the design driver comes from an environmental specification rather than from a strength margin, which is precisely the kind of requirement gradient-based optimization handles well.

Electromagnetic and Photonic Structures

The same mathematics extends beyond mechanics. Topology optimization applied to Maxwell's equations underlies the inverse design of compact waveguide filters, antenna elements, and integrated photonic components, where an optimizer distributes dielectric or conductor to meet a spectral target. Devices reported in that literature are often far smaller than hand-designed equivalents of similar function, at the cost of geometry that is harder to interpret and to fabricate within lithographic design rules.

Designing Within Manufacturing Constraints

An optimizer honors only the constraints it is given. Because these methods are usually paired with additive manufacturing, the relevant limits come from the chosen process, and modern solvers accept many of them directly.

  • Minimum member size keeps features above process resolution and is imposed through the filter radius. Published design guidelines for laser powder bed fusion of metals put the smallest reliably built strut diameters in the range of a few tenths of a millimeter up to about half a millimeter, varying with alloy, strut inclination, and machine, and thin members are the first thing to distort.
  • Overhang angle constraints keep down-facing surfaces above the self-supporting limit, conventionally taken as forty-five degrees from horizontal for powder bed fusion, so that internal supports are unnecessary. The forty-five-degree figure is a rule of thumb rather than a physical boundary; tuned parameter sets and scan strategies build shallower overhangs, but the rule remains the safe default for a design that must transfer between machines. Self-supporting formulations build the constraint into the optimization rather than repairing it afterward.
  • Powder and support removal governs internal geometry. A closed internal lattice traps unfused powder that cannot be evacuated, so drainage paths and access ports must be part of the design, and supports inside a lattice cannot be reached for removal.
  • Build orientation affects strength, surface finish, distortion, and cost simultaneously. Powder bed parts are anisotropic, generally weaker and less ductile when loaded along the build direction, across the layer boundaries, so orientation must be fixed before the design is qualified rather than left to the build preparation technician.
  • Machining access matters because sealing faces, bearing bores, and threaded features are almost always finished conventionally. Those surfaces need stock allowance and tool clearance, which belong in the non-design region from the outset.
  • Extrusion, casting, and symmetry constraints apply when the target process is not additive. Topology optimization is entirely applicable to die-cast housings and extruded heat sinks provided the corresponding draft, parting line, and constant-cross-section constraints are enforced.

Standards work has caught up with practice, although the design documents offer guidance rather than acceptance criteria. The ISO/ASTM 52900 series supplies the vocabulary and framework for additive manufacturing, ISO/ASTM 52910 sets out requirements, guidelines, and recommendations for design generally, and ISO/ASTM 52911-1:2019 addresses design for laser-based powder bed fusion of metals specifically, with a companion part covering laser-based powder bed fusion of polymers. Programs in regulated industries should anchor their internal design rules to these documents rather than to vendor guidance alone, while recognizing that satisfying a recommendation in a design standard does not by itself discharge an airworthiness, medical, or customer qualification requirement.

From Optimized Result to Released Part

The raw output of a density-based optimization is a faceted, stair-stepped isosurface extracted from the element density field. It is not a model anyone should release. Three routes lead from that result to a controlled design.

Reinterpretation uses the optimizer's answer as guidance and rebuilds the part in parametric CAD with clean features. This yields the most maintainable model and frequently the cheapest part, because the engineer naturally simplifies toward what the shop can make; it also sacrifices some of the optimized performance. Direct reconstruction fits smooth surfaces to the optimized shape, preserving performance at the cost of a model that is difficult to edit later. Implicit or field-driven modeling keeps the geometry in the same representation used for the lattice, avoiding conversion entirely, and suits parts where the lattice is the design.

Whichever route is chosen, the resulting geometry must be re-analyzed as built, not as optimized. The optimizer's finite element model is a coarse idealization with filtered material properties; the released part must pass the same static, dynamic, thermal, and stress checks any other part would face, using a mesh suited to the final geometry. Skipping this step is the most common way an optimized part fails qualification.

Verification and Qualification

Optimized and lattice parts complicate inspection because most of the geometry is internal and much of it is not accessible to conventional metrology. X-ray computed tomography is the primary tool for verifying internal lattice geometry and detecting porosity, but resolution limits mean fine lattices are characterized statistically rather than strut by strut. Witness specimens built alongside production parts, on the same plate and in the same orientation, provide mechanical evidence that the process delivered the assumed properties.

Fatigue deserves specific attention. As-built additive surfaces are rough, and in a strut only a few tenths of a millimeter across, that roughness represents a large fraction of the cross-section and a severe stress concentrator. Nodes where struts intersect add geometric concentration on top of it. Consequently, lattice fatigue performance is often far below what bulk material data would suggest, and post-processing such as hot isostatic pressing to close internal porosity, or chemical and abrasive surface treatment, may be required to reach the assumed allowable. Design allowables for lattices should come from testing of representative specimens, not from handbook values for the alloy.

Where the Method Pays and Where It Does Not

Topology optimization and lattice design are not free improvements. They add analysis effort, specialized tooling and software, longer inspection cycles, and per-part costs that are often high compared with casting or machining. As a design-for-excellence decision, the question is where those costs are repaid.

The strongest cases share features: mass or volume carries a high shadow price, production quantities are low enough that tooling amortization does not favor conventional processes, part consolidation eliminates fasteners and joints along with their assembly labor and failure modes, and the thermal or dynamic requirement cannot be met by a conventional geometry at all. Satellite and airborne electronics, high-power conversion hardware with severe cooling demands, instrumentation, and motorsport hardware all fit this profile.

The weak cases are equally recognizable. High-volume consumer hardware is almost always better served by an injection-molded or die-cast housing, even a heavier one, because unit cost dominates. Parts whose geometry is fixed by interfaces rather than by loading have little design freedom for an optimizer to exploit. Applications requiring sealed pressure boundaries, cosmetic surfaces, or established qualification pedigrees face additional hurdles. And a part that is already lightly loaded has no meaningful mass to remove.

Common Pitfalls

  • Optimizing for a single load case, then discovering the part is weak in a direction that was never applied. Real hardware sees handling, shipping, and assembly loads that no functional requirement lists.
  • Minimizing compliance and assuming the stress requirement follows. It does not; stress must be constrained explicitly or verified afterward.
  • Choosing a filter radius from meshing convenience rather than from process capability, producing features the process cannot build.
  • Treating a homogenized lattice as an isotropic solid, which misestimates both stiffness and load distribution.
  • Designing internal lattices without powder evacuation paths, yielding a part that cannot be cleaned and therefore cannot be flown or shipped.
  • Using bulk-material fatigue allowables for thin as-built struts.
  • Comparing an optimized part against a poorly designed baseline. A well-executed conventional design is the honest benchmark.

Conclusion

Topology optimization determines where material belongs in a design space; lattice design determines what fills the space that remains. Together they convert structural and thermal requirements directly into geometry, and they reach performance that manual design rarely matches when mass, stiffness, natural frequency, or heat transfer per unit volume is the binding constraint. Density-based formulations such as SIMP dominate practice, with level-set and evolutionary methods serving particular problem classes, and filtering with an explicit length scale is what separates a converged design from a mesh-dependent artifact.

Applied as a design-for-excellence discipline, these methods succeed only when manufacturing, inspection, and cost constraints enter the problem statement rather than the review that follows it. Constrain minimum feature size and overhang angle from the start, plan powder and support removal, fix build orientation before qualification, re-analyze the final geometry rather than the optimizer's idealization, and derive lattice allowables from representative test specimens. Handled that way, the methods deliver parts that are lighter, cooler, and fewer in number. Handled otherwise, they deliver an elegant shape that cannot be built, cannot be inspected, and cannot be afforded.

Related Topics