Topology optimization

Last updated April 20, 2024

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, bio-chemical 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.

Problem statement

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

{\begin{aligned}&{\underset {\rho }{\operatorname {minimize} }}&&F=F(\mathbf {u(\rho ),\rho } )=\int _{\Omega }f(\mathbf {u(\rho ),\rho } )\mathrm {d} V\\&\operatorname {subject\;to} &&G_{0}(\rho )=\int _{\Omega }\rho \mathrm {d} V-V_{0}\leq 0\\&&&G_{j}(\mathbf {u} (\rho ),\rho )\leq 0{\text{ with }}j=1,...,m\end{aligned}}

The problem statement includes the following:

An objective function $F(\mathbf {u(\rho ),\rho } )$ . 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 $\rho (\mathbf {x} )$ . Material is either present, indicated by a 1, or absent, indicated by a 0. $\mathbf {u} =\mathbf {u} (\mathbf {\rho } )$ is a state field that satisfies a linear or nonlinear state equation depending on $\rho$ .
The design space $(\Omega )$ . This indicates the allowable volume within which the design can exist. Assembly and packaging requirements, 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.
$\scriptstyle m$ constraints $G_{j}(\mathbf {u} (\rho ),\rho )\leq 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 $\mathbf {u(\rho )}$ 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 material density of one indicates the presence of material, while zero 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 30000 variables have been reported.^[2]

Solving the problem with continuous variables

The earlier stated 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 zero and one. 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\;=\;E_{0}\,+\,\rho ^{p}(E_{1}-E_{0})$ . 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, $E_{0}$ , to make sure 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 software on the market. Most of them use topology optimization as a hint how the optimal design should look like, and manual geometry re-construction is required. There are a few solutions which produce optimal designs ready for Additive Manufacturing.

Examples

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 we would not want to spend a lot of money on the material. To obtain the total material utilized, 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 _{\rho }\;\int _{\Omega }{\frac {1}{2}}\mathbf {\sigma

subject to:

$\rho \,\in \,[0,\,1]$
$\int _{\Omega }\rho \,\mathrm {d} \Omega \;\leq \;V^{*}$
$\mathbf {\nabla } \cdot \mathbf {\sigma } \,+\,\mathbf {F} \;=\;{\mathbf {0} }$
$\mathbf {\sigma } \;=\;{\mathsf {C}}:\mathbf {\varepsilon }$

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

Mesh dependency—Mesh Dependency means that the design obtained on one mesh is not the one that will be obtained on another mesh. The features of the design become more intricate as the mesh gets refined.^[7]
Numerical instabilities—The selection of region in the form of a chess board.^[8]

Some techniques such as filtering based on image processing^[9] 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.^[10]

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-interaction and to take such effects into consideration is therefore critical in the design of many engineering applications. Topology optimisation for fluid structure interaction problems has been studied in e.g. references^[11]^[12]^[13] and.^[14] Design solutions solved for different Reynolds numbers are shown below. The design solutions depend on the fluid flow with indicate that the coupling between the fluid and the structure is resolved in the design problems.

Design solution and velocity field for Re=1

Fluid-Structure-Interaction-Topology-Optimization-2.png

Design solution and velocity field for Re=5

Fluid-structure-interaction-pressure-field-topology-optimization.png

Design solution and pressure field for Re=10

Fluid-structure-interaction-pressure-field-topology-optimization-4.png

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. Wall-flow-problem-topology-optimization-for-fluid-structure-interaction-problems.png — 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. 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. Fluid-structure-interaction-design-evolution.gif — Design evolution for a fluid-structure-interaction problem from reference. 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-sketch.png — 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. Topology-optimization-off-diagonal-design-evolution.gif — 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. Evolution-design solution.gif — 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-physic problem which concerns the interaction and coupling between electric and thermal energy in semi conducting 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.^[15] By spatially distributing two thermoelectric materials in a two dimensional design space with a topology optimisation methodology,^[16] it is possible to exceed performance of the constitutive thermoelectric materials for thermoelectric coolers and thermoelectric generators.^[17]

3F3D Form Follows Force 3D Printing

The current proliferation of 3D printer technology has allowed designers and engineers to use topology optimization techniques when designing new products. Topology optimization combined with 3D printing can result in less weight, improved structural performance and shortened design-to-manufacturing cycle. As the designs, while efficient, might not be realisable with more traditional manufacturing techniques.^{[ citation needed ]}

Internal contact

Internal contact can be included in topology optimization by applying the third medium contact method.^[18]^[19]^[20] 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.

Related Research Articles

Continuum mechanics is a branch of mechanics that deals with the deformation of and transmission of forces through materials modeled as a continuous medium rather than as discrete particles. The French mathematician Augustin-Louis Cauchy was the first to formulate such models in the 19th century.

<span class="mw-page-title-main">Navier–Stokes equations</span> Equations describing the motion of viscous fluid substances

The Navier–Stokes equations are partial differential equations which describe the motion of viscous fluid substances. They were named after French engineer and physicist Claude-Louis Navier and the Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842–1850 (Stokes).

The vorticity equation of fluid dynamics describes the evolution of the vorticity $ω$ of a particle of a fluid as it moves with its flow; that is, the local rotation of the fluid. The governing equation is:

The finite volume method (FVM) is a method for representing and evaluating partial differential equations in the form of algebraic equations. In the finite volume method, volume integrals in a partial differential equation that contain a divergence term are converted to surface integrals, using the divergence theorem. These terms are then evaluated as fluxes at the surfaces of each finite volume. Because the flux entering a given volume is identical to that leaving the adjacent volume, these methods are conservative. Another advantage of the finite volume method is that it is easily formulated to allow for unstructured meshes. The method is used in many computational fluid dynamics packages. "Finite volume" refers to the small volume surrounding each node point on a mesh.

In continuum mechanics, the Froude number is a dimensionless number defined as the ratio of the flow inertia to the external field. The Froude number is based on the speed–length ratio which he defined as:

The Rankine–Hugoniot conditions, also referred to as Rankine–Hugoniot jump conditions or Rankine–Hugoniot relations, describe the relationship between the states on both sides of a shock wave or a combustion wave in a one-dimensional flow in fluids or a one-dimensional deformation in solids. They are named in recognition of the work carried out by Scottish engineer and physicist William John Macquorn Rankine and French engineer Pierre Henri Hugoniot.

In fluid mechanics, or more generally continuum mechanics, incompressible flow refers to a flow in which the material density of each fluid parcel — an infinitesimal volume that moves with the flow velocity — is time-invariant. An equivalent statement that implies incompressible flow is that the divergence of the flow velocity is zero.

In fluid mechanics, the Taylor–Proudman theorem states that when a solid body is moved slowly within a fluid that is steadily rotated with a high angular velocity $, the fluid velocity will be uniform along any line parallel to the axis of rotation. must be large compared to the movement of the solid body in order to make the Coriolis force large compared to the acceleration terms.$

<span class="mw-page-title-main">Smoothed-particle hydrodynamics</span> Method of hydrodynamics simulation

Smoothed-particle hydrodynamics (SPH) is a computational method used for simulating the mechanics of continuum media, such as solid mechanics and fluid flows. It was developed by Gingold and Monaghan and Lucy in 1977, initially for astrophysical problems. It has been used in many fields of research, including astrophysics, ballistics, volcanology, and oceanography. It is a meshfree Lagrangian method, and the resolution of the method can easily be adjusted with respect to variables such as density.

Stokes flow, also named creeping flow or creeping motion, is a type of fluid flow where advective inertial forces are small compared with viscous forces. The Reynolds number is low, i.e. $. This is a typical situation in flows where the fluid velocities are very slow, the viscosities are very large, or the length-scales of the flow are very small. Creeping flow was first studied to understand lubrication. In nature, this type of flow occurs in the swimming of microorganisms and sperm. In technology, it occurs in paint, MEMS devices, and in the flow of viscous polymers generally.$

<span class="mw-page-title-main">Crash simulation</span> Virtual recreation of a destructive car crash

A crash simulation is a virtual recreation of a destructive crash test of a car or a highway guard rail system using a computer simulation in order to examine the level of safety of the car and its occupants. Crash simulations are used by automakers during computer-aided engineering (CAE) analysis for crashworthiness in the computer-aided design (CAD) process of modelling new cars. During a crash simulation, the kinetic energy, or energy of motion, that a vehicle has before the impact is transformed into deformation energy, mostly by plastic deformation (plasticity) of the car body material, at the end of the impact.

The Navier–Stokes existence and smoothness problem concerns the mathematical properties of solutions to the Navier–Stokes equations, a system of partial differential equations that describe the motion of a fluid in space. Solutions to the Navier–Stokes equations are used in many practical applications. However, theoretical understanding of the solutions to these equations is incomplete. In particular, solutions of the Navier–Stokes equations often include turbulence, which remains one of the greatest unsolved problems in physics, despite its immense importance in science and engineering.

Limited-memory BFGS is an optimization algorithm in the family of quasi-Newton methods that approximates the Broyden–Fletcher–Goldfarb–Shanno algorithm (BFGS) using a limited amount of computer memory. It is a popular algorithm for parameter estimation in machine learning. The algorithm's target problem is to minimize $over unconstrained values of the real-vector where is a differentiable scalar function.$

Robust optimization is a field of mathematical optimization theory that deals with optimization problems in which a certain measure of robustness is sought against uncertainty that can be represented as deterministic variability in the value of the parameters of the problem itself and/or its solution. It is related to, but often distinguished from, probabilistic optimization methods such as chance-constrained optimization.

In the numerical solution of partial differential equations, a topic in mathematics, the spectral element method (SEM) is a formulation of the finite element method (FEM) that uses high-degree piecewise polynomials as basis functions. The spectral element method was introduced in a 1984 paper by A. T. Patera. Although Patera is credited with development of the method, his work was a rediscovery of an existing method

The derivation of the Navier–Stokes equations as well as its application and formulation for different families of fluids, is an important exercise in fluid dynamics with applications in mechanical engineering, physics, chemistry, heat transfer, and electrical engineering. A proof explaining the properties and bounds of the equations, such as Navier–Stokes existence and smoothness, is one of the important unsolved problems in mathematics.

The Clausius–Duhem inequality is a way of expressing the second law of thermodynamics that is used in continuum mechanics. This inequality is particularly useful in determining whether the constitutive relation of a material is thermodynamically allowable.

JCMsuite is a finite element analysis software package for the simulation and analysis of electromagnetic waves, elasticity and heat conduction. It also allows a mutual coupling between its optical, heat conduction and continuum mechanics solvers. The software is mainly applied for the analysis and optimization of nanooptical and microoptical systems. Its applications in research and development projects include dimensional metrology systems, photolithographic systems, photonic crystal fibers, VCSELs, Quantum-Dot emitters, light trapping in solar cells, and plasmonic systems. The design tasks can be embedded into the high-level scripting languages MATLAB and Python, enabling a scripting of design setups in order to define parameter dependent problems or to run parameter scans.

Ole Sigmund is a Danish Professor in Mechanical Engineering who has made fundamental contributions to the field of topology optimization, including microstructure design, nano optics, photonic crystals, Matlab code, acoustics, and fluids. In 2003 he co-authored the highly cited book "Topology Optimization: Theory, Methods and Applications" with Martin P. Bendsøe. His research group was the first to achieve giga-resolution topology optimization, making it for the first time possible to optimize an entire Boeing 777 wing structure.

The third medium contact(TMC) is an implicit formulation for contact mechanics. Contacting bodies are embedded in a highly compliant medium (the third medium), which becomes increasingly stiff under compression. The stiffening of the third medium allows tractions to be transferred between the contacting bodies when the third medium between the bodies is compressed. In itself, the method is inexact; however, in contrast to most other contact methods, the third medium approach is continuous and differentiable, which makes it applicable to applications such as topology optimization.

References

↑ Sigmund, Ole; Maute, Kurt (2013). "Topology optimization approaches". Structural and Multidisciplinary Optimization. 48 (6): 1031–1055. doi:10.1007/s00158-013-0978-6. S2CID 124426387.
↑ Beckers, M. (1999). "Topology optimization using a dual method with discrete variables" (PDF). Structural Optimization. 17: 14–24. doi:10.1007/BF01197709. S2CID 122845784.
↑ Bendsøe, M. P. (1989). "Optimal shape design as a material distribution problem". Structural Optimization. 1 (4): 193–202. doi:10.1007/BF01650949. S2CID 18253872.
↑ , a monograph of the subject.
↑ Bendsøe, M. P.; Sigmund, O. (1999). "Material interpolation schemes in topology optimization" (PDF). Archive of Applied Mechanics. 69 (9–10): 635–654. Bibcode:1999AAM....69..635B. doi:10.1007/s004190050248. S2CID 11368603.
↑ 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
↑ Allaire, Grégoire; Henrot, Antoine (May 2001). "On some recent advances in shape optimization". Comptes Rendus de l'Académie des Sciences. Series IIB - Mechanics. 329 (5). Elsevier: 383–396. Bibcode:2001CRASB.329..383A. doi:10.1016/S1620-7742(01)01349-6. ISSN 1620-7742 . Retrieved 2021-09-12.
↑ Shukla, Avinash; Misra, Anadi; Kumar, Sunil (September 2013). "Checkerboard Problem in Finite Element Based Topology Optimization". International Journal of Advances in Engineering & Technology. 6 (4). CiteSeer: 1769–1774. CiteSeerX 10.1.1.670.6771 . ISSN 2231-1963 . Retrieved 2022-02-14.
↑ Bourdin, Blaise (2001-03-30). "Filters in topology optimization". International Journal for Numerical Methods in Engineering. 50 (9). Wiley: 2143–2158. Bibcode:2001IJNME..50.2143B. doi:10.1002/nme.116. ISSN 1097-0207. S2CID 38860291 . Retrieved 2020-08-02.
↑ Sigmund, Ole; Maute, Kurt (October 2012). "Sensitivity filtering from a continuum mechanics perspective". Structural and Multidisciplinary Optimization. 46 (4). Springer: 471–475. doi:10.1007/s00158-012-0814-4. ISSN 1615-1488. S2CID 253680268 . Retrieved 2021-06-17.
↑ 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. Bibcode:2010IJNME..82..591Y. doi:10.1002/nme.2777. S2CID 122993997.
↑ 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.
↑ 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. S2CID 124632210.
1 2 Lundgaard, Christian; Alexandersen, Joe; Zhou, Mingdong; Andreasen, Casper Schousboe; Sigmund, Ole (2018). "Revisiting density-based topology optimization for fluid-structure-interaction problems" (PDF). Structural and Multidisciplinary Optimization. 58 (3): 969–995. doi:10.1007/s00158-018-1940-4. S2CID 125798826.
↑ Rowe, David Michael. Thermoelectrics handbook: macro to nano. CRC press, 2005.
↑ Lundgaard, Christian; Sigmund, Ole (2018). "A density-based topology optimization methodology for thermoelectric energy conversion problems" (PDF). Structural and Multidisciplinary Optimization. 57 (4): 1427–1442. doi:10.1007/s00158-018-1919-1. S2CID 126031362.
↑ Lundgaard, Christian; Sigmund, Ole; Bjørk, Rasmus (2018). "Topology Optimization of Segmented Thermoelectric Generators". Journal of Electronic Materials. 47 (12): 6959–6971. Bibcode:2018JEMat..47.6959L. doi:10.1007/s11664-018-6606-x. S2CID 105113187.
1 2 Frederiksen, Andreas Henrik; Sigmund, Ole; Poulios, Konstantinos (2023-10-07). "Topology optimization of self-contacting structures". Computational Mechanics. 73 (4): 967–981. arXiv: 2305.06750 . Bibcode:2023CompM..73..967F. doi:10.1007/s00466-023-02396-7. ISSN 1432-0924.
↑ Bluhm, Gore Lukas; Sigmund, Ole; Poulios, Konstantinos (2021-03-04). "Internal contact modeling for finite strain topology optimization". Computational Mechanics. 67 (4): 1099–1114. arXiv: 2010.14277 . Bibcode:2021CompM..67.1099B. doi:10.1007/s00466-021-01974-x. ISSN 0178-7675. S2CID 225076340.
↑ 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. Bibcode:2013CompM..52..837W. doi:10.1007/s00466-013-0848-5. ISSN 0178-7675. S2CID 254032357.

External links

Topology optimization animations

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Sigmund, Ole; Maute, Kurt (2013). "Topology optimization approaches". Structural and Multidisciplinary Optimization. 48 (6): 1031–1055. doi:10.1007/s00158-013-0978-6. S2CID 124426387.

[2] Beckers, M. (1999). "Topology optimization using a dual method with discrete variables" (PDF). Structural Optimization. 17: 14–24. doi:10.1007/BF01197709. S2CID 122845784.

[3] Bendsøe, M. P. (1989). "Optimal shape design as a material distribution problem". Structural Optimization. 1 (4): 193–202. doi:10.1007/BF01650949. S2CID 18253872.

[book-4] , a monograph of the subject.

[5] Bendsøe, M. P.; Sigmund, O. (1999). "Material interpolation schemes in topology optimization" (PDF). Archive of Applied Mechanics. 69 (9–10): 635–654. Bibcode:1999AAM....69..635B. doi:10.1007/s004190050248. S2CID 11368603.

[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] Allaire, Grégoire; Henrot, Antoine (May 2001). "On some recent advances in shape optimization". Comptes Rendus de l'Académie des Sciences. Series IIB - Mechanics. 329 (5). Elsevier: 383–396. Bibcode:2001CRASB.329..383A. doi:10.1016/S1620-7742(01)01349-6. ISSN 1620-7742 . Retrieved 2021-09-12.

[8] Shukla, Avinash; Misra, Anadi; Kumar, Sunil (September 2013). "Checkerboard Problem in Finite Element Based Topology Optimization". International Journal of Advances in Engineering & Technology. 6 (4). CiteSeer: 1769–1774. CiteSeerX 10.1.1.670.6771 . ISSN 2231-1963 . Retrieved 2022-02-14.

[9] Bourdin, Blaise (2001-03-30). "Filters in topology optimization". International Journal for Numerical Methods in Engineering. 50 (9). Wiley: 2143–2158. Bibcode:2001IJNME..50.2143B. doi:10.1002/nme.116. ISSN 1097-0207. S2CID 38860291 . Retrieved 2020-08-02.

[10] Sigmund, Ole; Maute, Kurt (October 2012). "Sensitivity filtering from a continuum mechanics perspective". Structural and Multidisciplinary Optimization. 46 (4). Springer: 471–475. doi:10.1007/s00158-012-0814-4. ISSN 1615-1488. S2CID 253680268 . Retrieved 2021-06-17.

[11] 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. Bibcode:2010IJNME..82..591Y. doi:10.1002/nme.2777. S2CID 122993997.

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

[13] 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. S2CID 124632210.

[Lundgaard_FSI-14] 1 2 Lundgaard, Christian; Alexandersen, Joe; Zhou, Mingdong; Andreasen, Casper Schousboe; Sigmund, Ole (2018). "Revisiting density-based topology optimization for fluid-structure-interaction problems" (PDF). Structural and Multidisciplinary Optimization. 58 (3): 969–995. doi:10.1007/s00158-018-1940-4. S2CID 125798826.

[15] Rowe, David Michael. Thermoelectrics handbook: macro to nano. CRC press, 2005.

[16] Lundgaard, Christian; Sigmund, Ole (2018). "A density-based topology optimization methodology for thermoelectric energy conversion problems" (PDF). Structural and Multidisciplinary Optimization. 57 (4): 1427–1442. doi:10.1007/s00158-018-1919-1. S2CID 126031362.

[17] Lundgaard, Christian; Sigmund, Ole; Bjørk, Rasmus (2018). "Topology Optimization of Segmented Thermoelectric Generators". Journal of Electronic Materials. 47 (12): 6959–6971. Bibcode:2018JEMat..47.6959L. doi:10.1007/s11664-018-6606-x. S2CID 105113187.

[:0-18] 1 2 Frederiksen, Andreas Henrik; Sigmund, Ole; Poulios, Konstantinos (2023-10-07). "Topology optimization of self-contacting structures". Computational Mechanics. 73 (4): 967–981. arXiv: 2305.06750 . Bibcode:2023CompM..73..967F. doi:10.1007/s00466-023-02396-7. ISSN 1432-0924.

[19] Bluhm, Gore Lukas; Sigmund, Ole; Poulios, Konstantinos (2021-03-04). "Internal contact modeling for finite strain topology optimization". Computational Mechanics. 67 (4): 1099–1114. arXiv: 2010.14277 . Bibcode:2021CompM..67.1099B. doi:10.1007/s00466-021-01974-x. ISSN 0178-7675. S2CID 225076340.

[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. Bibcode:2013CompM..52..837W. doi:10.1007/s00466-013-0848-5. ISSN 0178-7675. S2CID 254032357.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]

[14]

[15]

[16]

[17]

[18]

[19]

[20]