Chapter 4 · Materials & the Environment

Brackets, Metals and Carbon

Which bracket, in which metal, for the least carbon? 379 crowd-designed jet-engine brackets, four aerospace alloys, 1 516 finite-element runs and their footprints.

1Which bracket, in which metal, for the least carbon?

The first three chapters chose a material for a fixed shape class: a tie, a beam, a panel. Real parts are shaped and specified together. Topology optimisation changes the shape, and the material changes what the shape must do. The TOSM benchmark (Topology Optimisation & Sustainable Metals) puts numbers on that coupling. It takes a large set of independently designed parts, simulates each one in several alloys, and attaches production and lifetime footprints to every result.

Every figure in this chapter is computed in your browser from the benchmark's own result files: 379 brackets × 4 alloys = 1 516 runs.

2The brackets

In 2013 GE ran an open challenge to design the lightest possible jet-engine bracket. Every entry had the same four bolt holes and the same interface point, and had to keep its peak stress below the yield stress of Ti-6Al-4V under four load cases: vertical, diagonal, horizontal and torsional. The SimJEB dataset collects 381 of those designs with clean meshes and finite-element models. Of these, 379 run cleanly in every alloy and form the benchmark. Their shapes range from chunky blocks to lattice-like organic forms; the source article includes a 3-D gallery of them.

The SimJEB decks are Altair OptiStruct models. A notebook (OptistructAuto) rewrites the material card of each deck for each alloy and runs linear static solves. Each run records the bracket's mass, its peak displacement and its peak von Mises stress over the four load cases.

3Four alloys and their footprints

The alloys and their per-kilogram production footprints come from the eco-design database of Yepes Llorente and Morlier.

Alloyρ (kg/m³)E (GPa)σy (MPa)CO₂ (kg/kg)Energy (MJ/kg)Water (L/kg)Cost ($/kg)

For carbon, the benchmark adds a use phase. A kilogram carried on an aircraft for its whole service life burns fuel, and Duriez et al. estimate that at 98.8 tonnes of CO₂ per kilogram of structure.

\text{CO}_{2,T} = M\cdot \text{CO}_{2,s} + 98.8\ \tfrac{\text{t}}{\text{kg}}\cdot M

Titanium's production footprint is 19.7 kg CO₂/kg, or 0.0197 t/kg: the use phase outweighs it by a factor of about 5 000. Minimising lifetime CO₂ is, to within a rounding error, minimising mass.

Figure 4.1Production against use phase. Lifetime CO₂ of the median-mass bracket in each alloy, split into production and use. Slide the factor to zero to see production on its own.

4What changes when you swap the alloy

The geometry is fixed and the analysis is linear, so swapping alloys changes outputs in predictable ways. Mass scales with density. Displacement scales roughly with 1/E. Peak stress hardly moves: with prescribed loads and one homogeneous material, linear elastic stresses do not depend on the modulus at all, and only weakly on Poisson's ratio. That is the source of the 1–3% shifts.

Figure 4.2Change relative to titanium. Median over the 379 brackets of each output's change when the same bracket is made in another alloy. Aluminium is about 37% lighter and cheaper in carbon, but deflects 56% more. Stainless steel and beryllium copper are nearly twice as heavy.

5Most designs are over-stressed

A design is admissible when its peak stress stays below the yield stress divided by a safety factor, \sigma \le \sigma_y/\text{SF}. Because stress depends on geometry and not on alloy, the alloy choice decides admissibility through \sigma_y alone.

Figure 4.3Utilisation by alloy. Cumulative share of brackets whose peak stress is below a given fraction of yield, σ/σy. The vertical line is the admissibility limit 1/SF; curves left of it count as passing.

Even in titanium, the alloy the challenge was set for, the median bracket peaks at about 1.6 times its yield stress. At SF = 1.25 only 22 of the 1 516 runs pass, 21 of them in titanium. No aluminium bracket passes at any safety factor. Part of this is the solver: peak von Mises stress at sharp re-entrant corners is mesh-dependent, and with SimJEB's own stress values 74 titanium brackets would pass instead of 21. The rest is material: every alloy here has a far lower yield stress than titanium, while the shapes were drawn for titanium.

6Finding the trade-off

The benchmark poses a multi-objective problem: minimise peak stress, lifetime CO₂, energy and water, subject to the stress constraint. The original notebooks solve it with NSGA-II through pymoo. With 1 516 discrete candidates there is no need for a genetic algorithm: the exact Pareto front can be found by enumeration, which the figure below does live.

Figure 4.4Pareto explorer. Every bracket–alloy run, with admissible runs in colour and the exact non-dominated front connected. Hover a point for its bracket number; click it to link the same bracket across all four alloys. Switch the constraint to the notebook's version to reproduce the original optimiser's result.

On the stress–CO₂ plane at SF = 1.25, the front holds five titanium brackets. Bracket 120 is the lightest admissible design, at 65.9 t of lifetime CO₂ and 666 MPa. Bracket 271 trades 32 t more CO₂ for a peak stress of 404 MPa, and bracket 180 buys only 5 MPa more margin for another 45 t.

7A note on the constraint sign

pymoo treats a constraint g(x) \le 0 as satisfied. The notebook codes the stress constraint as

g = -\sigma + \sigma_y/\text{SF},

which is non-positive exactly when \sigma \ge \sigma_y/\text{SF}, the opposite of the intent. The optimiser therefore searched only among over-stressed designs. Its saved result, ParetoFrontResults.xlsx, lists fourteen Al-2024 designs stressed between 2.4 and 57 times their yield limit, and re-running the notebook's setup (population 100, 500 generations, seed 1) reproduces them. Choose “Notebook (inverted)” in Figure 4.4: the exact enumeration over the inverted feasible set returns the same fourteen Al-2024 brackets, from bracket 180 at 2.4 times yield to bracket 373 at 57 times. The fix is one line:

G = stress - yield_stress / safety_factor

With the sign corrected, the front moves to titanium.

8Three objectives at once

Stress, CO₂ and displacement pull in three directions. Aluminium wins on CO₂ and loses on displacement; steel wins on displacement and loses on CO₂; titanium sits between them and is the only alloy that is admissible in practice. Pick displacement and CO₂ as the axes in Figure 4.4 and lower the safety factor to 1 to see the alloys separate into bands. The source article shows the three-objective fronts as rotatable 3-D plots.

9Where this goes next

Reproducing these figures

All figures read FinalResults.xlsx (1 516 runs) and Materials.xlsx from the TOSM repository, converted to data/tosm.js in this site. Lifetime CO₂, energy, water and cost are recomputed from mass and the per-kilogram values; they match the repository's columns to machine precision. Bracket models and renderings are from the SimJEB repository (MIT licence).