Topology optimization

Short description: Mathematical method for optimizing material layout under given conditions

Topology optimization is a mathematical method that optimizes material layout within a given design space, for a given set of loads, boundary conditions, and constraints with the goal of maximizing the performance of the system. Topology optimization is different from shape optimization and sizing optimization in the sense that the design can attain any shape within the design space, instead of dealing with predefined configurations.

The conventional topology optimization formulation uses a finite element method (FEM) to evaluate the design performance. The design is optimized using either gradient-based mathematical-programming techniques such as the optimality criteria algorithm and the method of moving asymptotes or non-gradient-based algorithms such as genetic algorithms.

Topology optimization has a wide range of applications in aerospace, mechanical, biochemical, and civil engineering. Currently, engineers mostly use topology optimization at the concept level of a design process. Due to the free forms that naturally occur, the result is often difficult to manufacture. For that reason, the result emerging from topology optimization is often fine-tuned for manufacturability. Adding constraints to the formulation in order to increase the manufacturability is an active field of research. In some cases, results from topology optimization can be directly manufactured using additive manufacturing; topology optimization is thus a key part of design for additive manufacturing.

## Contents

* 1 Problem statement
* 2 Implementation methodologies
* 2.1 Solving with discrete/binary variables
* 2.2 Solving the problem with continuous variables
* 2.3 Commercial software
* 3 Examples
* 3.1 Structural compliance
* 3.2 Multiphysics problems
* 3.2.1 Fluid-structure interaction
* 3.2.2 Thermoelectric energy conversion
* 3.3 3F3D: Form Follows Force 3D Printing
* 3.4 Internal contact
* 4 See also
* 5 References
* 6 Further reading
* 7 External links

## Problem statement

A topology-optimization problem can be written in the general form of an optimization problem as:

minimizeρF=F(𝐮(ρ),ρ)=∫Ωf(𝐮(ρ),ρ)dVsubjecttoG0(ρ)=∫ΩρdV−V0≤0Gj(𝐮(ρ),ρ)≤0 with j=1,...,m.

The problem statement includes the following:

* An objective function F(𝐮(ρ),ρ). This function represents the quantity that is being minimized for best performance. The most common objective function is compliance, where minimizing compliance leads to maximizing the stiffness of a structure.
* The material distribution as a problem variable. This is described by the density of the material at each location, ρ(𝐱). Material is either present, indicated by a 1, or absent, indicated by a 0. 𝐮=𝐮(ρ) is a state field that satisfies a linear or nonlinear state equation depending on ρ.
* The design space (Ω). This indicates the allowable volume within which the design can exist. Assembly and packaging requirements and human and tool accessibility are some of the factors that need to be considered in identifying this space. With the definition of the design space, regions or components in the model that cannot be modified during the course of the optimization are considered as non-design regions.
* m constraints Gj(𝐮(ρ),ρ)≤0 a characteristic that the solution must satisfy. Examples are the maximum amount of material to be distributed (volume constraint) or maximum stress values.

Evaluating 𝐮(ρ) often includes solving a differential equation. This is most commonly done using the finite-element method, since these equations do not have a known analytical solution.

## Implementation methodologies

There are various implementation methodologies that have been used to solve topology-optimization problems.

### Solving with discrete/binary variables

Solving topology-optimization problems in a discrete sense is done by discretizing the design domain into finite elements. The material densities inside these elements are then treated as the problem variables. In this case, a material density of 1 indicates the presence of material, while 0 indicates an absence of material. Owing to the attainable topological complexity of the design being dependent on the number of elements, a large number is preferred. Large numbers of finite elements increases the attainable topological complexity, but come at a cost. Firstly, solving the FEM system becomes more expensive. Secondly, algorithms that can handle a large number (several thousands of elements is not uncommon) of discrete variables with multiple constraints are unavailable. Moreover, they are impractically sensitive to parameter variations.[1] In literature, problems with up to 30,000 variables have been reported.[2]

### Solving the problem with continuous variables

The aforementioned complexities with solving topology optimization problems using binary variables has caused the community to search for other options. One is the modelling of the densities with continuous variables. The material densities can now also attain values between 0 and 1. Gradient-based algorithms that handle large amounts of continuous variables and multiple constraints are available. But the material properties have to be modelled in a continuous setting. This is done through interpolation. One of the most implemented interpolation methodologies is the Solid Isotropic Material with Penalisation method (SIMP).[3][4] This interpolation is essentially a power law: E=E0+ρp(E1−E0). It interpolates the Young's modulus of the material to the scalar selection field. The value of the penalisation parameter p is generally taken between [1,3]. This has been shown to confirm the micro-structure of the materials.[5] In the SIMP method, a lower bound on the Young's modulus is added, E0, to make sure that the derivatives of the objective function are non-zero when the density becomes zero. The higher the penalisation factor, the more SIMP penalises the algorithm in the use of non-binary densities. Unfortunately, the penalisation parameter also introduces non-convexities.[6]

### Commercial software

There are several commercial topology-optimization softwares on the market. Most of them use topology optimization as a hint to how the optimal design should look, and manual geometry re-construction is required. There are a few solutions which produce optimal designs ready for additive manufacturing.[7]

## Examples

File:Checkerboards in Topology Optimization.tif File:Topology Optimization with filtereing.tif

Topology optimization of a compliance problem

### Structural compliance

A stiff structure is one that has the least possible displacement when given certain set of boundary conditions. A global measure of the displacements is the strain energy (also called compliance) of the structure under the prescribed boundary conditions. The lower the strain energy, the higher the stiffness of the structure. So, the objective function of the problem is to minimize the strain energy.

On a broad level, one can visualize that the more the material, the less the deflection, as there will be more material to resist the loads. So, the optimization requires an opposing constraint, the volume constraint. This is in reality a cost factor, as one would not want to spend a lot of money on the material. To obtain the total material used, an integration of the selection field over the volume can be done.

Finally, the elasticity-governing differential equations are plugged in so as to get the final problem statement:

minρ∫Ω12σ:εdΩ

subject to:

* ρ∈[0,1]
* ∫ΩρdΩ≤V*
* ∇⋅σ+𝐅=𝟎
* σ=C:ε

But, a straightforward implementation in the finite-element framework of such a problem is still infeasible due to issues such as:

1. Mesh dependency—The design obtained on one mesh can be very different from that obtained on another mesh. The features of the design become more intricate as the mesh gets refined.[8]
2. Numerical instabilities—A small change to an input parameter can produce a large change in the computed solution.[9]

Some techniques such as filtering based on image processing[10] are currently being used to alleviate some of these issues. Although it seemed like this was purely a heuristic approach for a long time, theoretical connections to nonlocal elasticity have been made to support the physical sense of these methods.[11]

### Multiphysics problems

#### Fluid-structure interaction

Fluid-structure interaction is a strongly coupled phenomenon and concerns the interaction between a stationary or moving fluid and an elastic structure. Many engineering applications and natural phenomena are subject to fluid-structure interactions, and taking such effects into consideration is therefore critical in the design of many engineering applications. Topology optimization for fluid-structure interaction problems has been studied in e.g. references,[12][13][14] and.[15] Design solutions solved for different Reynolds numbers are shown below. The design solutions depend on the fluid flow, which indicates that the coupling between the fluid and the structure is resolved in the design problems.

Design solution and velocity field for Re=1

Design solution and velocity field for Re=5

Design solution and pressure field for Re=10

Design solution and pressure field for Re=40

Design solutions for different Reynolds number for a wall inserted in a channel with a moving fluid.

Sketch of the well-known wall problem. The objective of the design problem is to minimize the structural compliance. Design evolution for a fluid-structure-interaction problem from reference.[15] The objective of the design problem is to minimize the structural compliance. The fluid-structure-interaction problem is modelled with Navier-Cauchy and Navier-Stokes equations.

#### Thermoelectric energy conversion

A sketch of the design problem. The aim of the design problem is to spatially distribute two materials, Material A and Material B, to maximise a performance measure such as cooling power or electric power output Design evolution for an off-diagonal thermoelectric generator. The design solution of an optimisation problem solved for electric power output. The performance of the device has been optimised by distributing Skutterudite (yellow) and bismuth telluride (blue) with a density-based topology optimisation methodology. The aim of the optimisation problem is to maximise the electric power output of the thermoelectric generator. Design evolution for a thermoelectric cooler. The aim of the design problem is to maximise the cooling power of the thermoelectric cooler.

Thermoelectricity is a multi-physics phenomenon which concerns the interaction and coupling between electric and thermal energy in semiconducting materials. Thermoelectric energy conversion can be described by two separately identified effects: the Seebeck effect and the Peltier effect. The Seebeck effect concerns the conversion of thermal energy into electric energy, and the Peltier effect concerns the conversion of electric energy into thermal energy.[16] By spatially distributing two thermoelectric materials in a two-dimensional design space with a topology-optimization methodology,[17] it is possible to exceed the performance of the constitutive thermoelectric materials for thermoelectric coolers and thermoelectric generators.[18]

### 3F3D: Form Follows Force 3D Printing

### Internal contact

File:3D Topology Optimization with Internal Contact.webm

Design development and deformation of self-engaging hooks resulting from topology optimization of a contact problem using the TMC method.[19]

Internal contact can be included in topology optimization by applying the third medium contact method.[20][21] The third medium contact (TMC) method is an implicit contact formulation that is continuous and differentiable. This makes TMC suitable for use with gradient-based approaches to topology optimization. Monolithic[22] as well as staggered approaches,[19][23] which are more common in topology optimization, have been used to create various designs with internal contact. Recently, thermal contact has been included in the TMC topology optimization framework.[24]

## See also

* Generative design

## References

1. ↑ Sigmund, Ole; Maute, Kurt (2013). "Topology optimization approaches". Structural and Multidisciplinary Optimization 48 (6): 1031–1055. doi:10.1007/s00158-013-0978-6.
2. ↑ Beckers, M. (1999). "Topology optimization using a dual method with discrete variables". Structural Optimization 17: 14–24. doi:10.1007/BF01197709. http://empslocal.ex.ac.uk/people/staff/reverson/uploads/MoodSwings/beckers.pdf.
3. ↑ Bendsøe, M. P. (1989). "Optimal shape design as a material distribution problem". Structural Optimization 1 (4): 193–202. doi:10.1007/BF01650949.
4. ↑ [1], a monograph of the subject.
5. ↑ Bendsøe, M. P.; Sigmund, O. (1999). "Material interpolation schemes in topology optimization". Archive of Applied Mechanics 69 (9–10): 635–654. doi:10.1007/s004190050248. Bibcode: 1999AAM....69..635B. http://www.giref.ulaval.ca/~deteix/bois/documents_references/bendsoe1999.pdf.
6. ↑ van Dijk, NP. Langelaar, M. van Keulen, F. Critical study of design parameterization in topology optimization; The influence of design parameterization on local minima.. 2nd International Conference on Engineering Optimization, 2010
7. ↑ "Topological optimization software for 3d printing". 23 September 2021. https://www.3dnatives.com/en/topological-optimization-software-for-3d-printing-230920214/.
8. ↑ Allaire, Grégoire; Henrot, Antoine (May 2001). "On some recent advances in shape optimization" (in en). Comptes Rendus de l'Académie des Sciences. Series IIB - Mechanics (Elsevier) 329 (5): 383–396. doi:10.1016/S1620-7742(01)01349-6. ISSN 1620-7742. Bibcode: 2001CRASB.329..383A. https://linkinghub.elsevier.com/retrieve/pii/S1620774201013496. Retrieved 2021-09-12.
9. ↑ Shukla, Avinash; Misra, Anadi; Kumar, Sunil (September 2013). "Checkerboard Problem in Finite Element Based Topology Optimization" (in en). International Journal of Advances in Engineering & Technology (CiteSeer) 6 (4): 1769–1774. ISSN 2231-1963. http://citeseerx.ist.psu.edu/viewdoc/download;jsessionid=5BECA0B81B176391016659D59276F7FE?doi=10.1.1.670.6771&rep=rep1&type=pdf. Retrieved 2022-02-14.
10. ↑ Bourdin, Blaise (2001-03-30). "Filters in topology optimization" (in en). International Journal for Numerical Methods in Engineering (Wiley) 50 (9): 2143–2158. doi:10.1002/nme.116. ISSN 1097-0207. Bibcode: 2001IJNME..50.2143B. http://doi.wiley.com/10.1002/nme.116. Retrieved 2020-08-02.
11. ↑ Sigmund, Ole; Maute, Kurt (October 2012). "Sensitivity filtering from a continuum mechanics perspective". Structural and Multidisciplinary Optimization (Springer) 46 (4): 471–475. doi:10.1007/s00158-012-0814-4. ISSN 1615-1488. http://link.springer.com/10.1007/s00158-012-0814-4. Retrieved 2021-06-17.
12. ↑ Yoon, Gil Ho (2010). "Topology optimization for stationary fluid-structure interaction problems using a new monolithic formulation". International Journal for Numerical Methods in Engineering 82 (5): 591–616. doi:10.1002/nme.2777. Bibcode: 2010IJNME..82..591Y.
13. ↑ Picelli, R.; Vicente, W.M.; Pavanello, R. (2017). "Evolutionary topology optimization for structural compliance minimization considering design-dependent FSI loads". Finite Elements in Analysis and Design 135: 44–55. doi:10.1016/j.finel.2017.07.005.
14. ↑ Jenkins, Nicholas; Maute, Kurt (2016). "An immersed boundary approach for shape and topology optimization of stationary fluid-structure interaction problems". Structural and Multidisciplinary Optimization 54 (5): 1191–1208. doi:10.1007/s00158-016-1467-5.
15. ↑ 15.0 15.1 Lundgaard, Christian; Alexandersen, Joe; Zhou, Mingdong; Andreasen, Casper Schousboe; Sigmund, Ole (2018). "Revisiting density-based topology optimization for fluid-structure-interaction problems". Structural and Multidisciplinary Optimization 58 (3): 969–995. doi:10.1007/s00158-018-1940-4. https://backend.orbit.dtu.dk/ws/files/163153999/grayscale_Lundgaard_C._Alexandersen_J._Zhou_M._Andreasen_C._S._Sigmund_O_2018_.pdf.
16. ↑ Rowe, David Michael. Thermoelectrics handbook: macro to nano. CRC press, 2005.
17. ↑ Lundgaard, Christian; Sigmund, Ole (2018). "A density-based topology optimization methodology for thermoelectric energy conversion problems". Structural and Multidisciplinary Optimization 57 (4): 1427–1442. doi:10.1007/s00158-018-1919-1. https://backend.orbit.dtu.dk/ws/files/163153924/grayscale_Lundgaard_C._Sigmund_O_2018_.pdf.
18. ↑ Lundgaard, Christian; Sigmund, Ole; Bjørk, Rasmus (2018). "Topology Optimization of Segmented Thermoelectric Generators". Journal of Electronic Materials 47 (12): 6959–6971. doi:10.1007/s11664-018-6606-x. Bibcode: 2018JEMat..47.6959L. https://www.researchgate.net/publication/323143969.
19. ↑ 19.0 19.1 Frederiksen, Andreas Henrik; Sigmund, Ole; Poulios, Konstantinos (2023-10-07). "Topology optimization of self-contacting structures" (in en). Computational Mechanics 73 (4): 967–981. doi:10.1007/s00466-023-02396-7. ISSN 1432-0924. Bibcode: 2023CompM..73..967F.
20. ↑ Wriggers, P.; Schröder, J.; Schwarz, A. (2013-03-30). "A finite element method for contact using a third medium". Computational Mechanics 52 (4): 837–847. doi:10.1007/s00466-013-0848-5. ISSN 0178-7675. Bibcode: 2013CompM..52..837W.
21. ↑ Frederiksen, Andreas H.; Rokoš, Ondřej; Poulios, Konstantinos; Sigmund, Ole; Geers, Marc G. D. (2024-12-01). "Adding friction to Third Medium Contact: A crystal plasticity inspired approach". Computer Methods in Applied Mechanics and Engineering 432. doi:10.1016/j.cma.2024.117412. ISSN 0045-7825. Bibcode: 2024CMAME.43217412F.
22. ↑ Bluhm, Gore Lukas; Sigmund, Ole; Poulios, Konstantinos (2021-03-04). "Internal contact modeling for finite strain topology optimization". Computational Mechanics 67 (4): 1099–1114. doi:10.1007/s00466-021-01974-x. ISSN 0178-7675. Bibcode: 2021CompM..67.1099B.
23. ↑ Frederiksen, Andreas H.; Dalklint, Anna; Sigmund, Ole; Poulios, Konstantinos (2025-03-01). "Improved third medium formulation for 3D topology optimization with contact". Computer Methods in Applied Mechanics and Engineering 436. doi:10.1016/j.cma.2024.117595. ISSN 0045-7825. Bibcode: 2025CMAME.43617595F.
24. ↑ Dalklint, Anna; Alexandersen, Joe; Frederiksen, Andreas Henrik; Poulios, Konstantinos; Sigmund, Ole (2025). "Topology Optimization of Contact-Aided Thermo-Mechanical Regulators" (in en). International Journal for Numerical Methods in Engineering 126 (2). doi:10.1002/nme.7661. ISSN 1097-0207. https://onlinelibrary.wiley.com/doi/full/10.1002/nme.7661.

## Further reading

* Pedersen, Claus B. W.; Allinger, Peter (2006). "Industrial Implementation and Applications of Topology Optimization and Future Needs". IUTAM Symposium on Topological Design Optimization of Structures, Machines and Materials. Solid Mechanics and Its Applications. 137. Springer. pp. 229–238. doi:10.1007/1-4020-4752-5_23. ISBN 978-1-4020-4729-9.
* Schramm, Uwe; Zhou, Ming (2006). "Recent Developments in the Commercial Implementation of Topology Optimization". IUTAM Symposium on Topological Design Optimization of Structures, Machines and Materials. Solid Mechanics and Its Applications. 137. Springer. pp. 239–248. doi:10.1007/1-4020-4752-5_24. ISBN 978-1-4020-4729-9.
* Mahdavi, A.; Balaji, R.; Frecker, M.; Mockensturm, E. M. (2006). "Topology optimization of 2D continua for minimum compliance using parallel computing". Structural and Multidisciplinary Optimization 32 (2): 121–132. doi:10.1007/s00158-006-0006-1.
* Leiva, Juan; Watson, Brian; Kosaka, Iku (1999). "Modern structural optimization concepts applied to topology optimization". 40th Structures, Structural Dynamics, and Materials Conference and Exhibit. American Institute of Aeronautics and Astronautics. doi:10.2514/6.1999-1388.

## External links

* Topology optimization animations

0.00

(0 votes)

Original source: https://en.wikipedia.org/wiki/Topology optimization. Read more

|

Retrieved from "https://handwiki.org/wiki/index.php?title=Topology_optimization&oldid=4205427"