:orphan:

..
   Copyright 2021 IRT Saint Exupéry, https://www.irt-saintexupery.com

   This work is licensed under the Creative Commons Attribution-ShareAlike 4.0
   International License. To view a copy of this license, visit
   http://creativecommons.org/licenses/by-sa/4.0/ or send a letter to Creative
   Commons, PO Box 1866, Mountain View, CA 94042, USA.

Topology optimization
=====================

In these examples, 2D topology optimization problems are solved thanks to GEMSEO.

Topology optimization aims at finding the "best" material layout withing a given design space.
It is especially useful in preliminary design phases when one only knows a component of the operating conditions and wants to
find good design candidates achieving stiffer structures within a given mass budget.
The formulation adopted in these examples is the Solid Isotropic Material with Penalization (SIMP) [1]_ approach.

The examples are the same from MATLAB and Python implementations described in [2]_ and [3]_ .

Problem Formulation
-------------------
Given a 2D design space with loads and boundary conditions,
let us take the example of an MBB structure:

.. image:: /_images/topology_optimization/TO_design_space_MBB.png

First, the solid design space is meshed with 2D bi-linear squared finite elements.
All the examples proposed here consider rectangular domains.
This means that one only needs to define the number of elements in horizontal (x) and vertical (y) direction.


.. image:: /_images/topology_optimization/TO_mesh_MBB.png

In the above figure it is possible to find the element and degree of freedoms numbering convention adopted for a
4x3 finite element mesh. What gives a very large design freedom to topology optimization is that a design variable
:math:`x \in \{0,1\}^N` is associated with each finite element. These variables are equal to 0 when the finite element is
void and equal to 1 when the finite element is filled with solid material.
In order to use the convergence rate of the gradient-based optimization solvers, the
design variable are relaxed :math:`x \in [0,1]^N`. To enforce a discrete solution at convergence,
the intermediate valuess of the design variables are penalized using the SIMP approach that introduces a power low relationship between
the local material density and the Young's modulus. In order to avoid numerical difficulties such as mesh-dependent solutions and
checkerboard patterns, a density filtering technique is implemented.
In these examples, topology optimization is employed to
minimize structural compliance subjected to a mass budget or equivalently a volume fraction target:

.. math:: \min_{x \in [0,1]^N}{F \cdot U(x)}
.. math:: s.t.
.. math:: \frac{1}{N}\sum_{i=1}^N{x_i}\leq \overline{V}
.. math:: K(x)U(x) = F

where :math:`F` is the load vector,
:math:`U` is the displacement vector,
:math:`N` is the number of elements,
:math:`\overline{V}` is the allowable volume fraction and :math:`K` is the stiffness matrix.



.. [1] Bendsøe, M. P. (1989). Optimal shape design as a material distribution problem. Structural optimization, 1(4), 193-202.
.. [2] Sigmund, O. (2001). A 99 line topology optimization code written in Matlab. Structural and multidisciplinary optimization, 21(2), 120-127.
.. [3] Andreassen, E., Clausen, A., Schevenels, M., Lazarov, B. S., & Sigmund, O. (2011). Efficient topology optimization in MATLAB using 88 lines of code. Structural and Multidisciplinary Optimization, 43(1), 1-16.

.. toctree::
   :glob:
   :maxdepth: 1
   :hidden:

   **



.. raw:: html

    <div class="sphx-glr-thumbnails">


.. raw:: html

    <div class="sphx-glr-thumbcontainer" tooltip="Solve a 2D L-shape topology optimization problem">

.. only:: html

  .. image:: /examples/topology_optimization/images/thumb/sphx_glr_plot_topology_optimization_L_shape_thumb.png
    :alt:

  :ref:`sphx_glr_examples_topology_optimization_plot_topology_optimization_L_shape.py`

.. raw:: html

      <div class="sphx-glr-thumbnail-title">Solve a 2D L-shape topology optimization problem</div>
    </div>


.. raw:: html

    <div class="sphx-glr-thumbcontainer" tooltip="Solve a 2D MBB topology optimization problem">

.. only:: html

  .. image:: /examples/topology_optimization/images/thumb/sphx_glr_plot_topology_optimization_MBB_thumb.png
    :alt:

  :ref:`sphx_glr_examples_topology_optimization_plot_topology_optimization_MBB.py`

.. raw:: html

      <div class="sphx-glr-thumbnail-title">Solve a 2D MBB topology optimization problem</div>
    </div>


.. raw:: html

    <div class="sphx-glr-thumbcontainer" tooltip="Solve a 2D short cantilever topology optimization problem">

.. only:: html

  .. image:: /examples/topology_optimization/images/thumb/sphx_glr_plot_topology_optimization_short_cantilever_thumb.png
    :alt:

  :ref:`sphx_glr_examples_topology_optimization_plot_topology_optimization_short_cantilever.py`

.. raw:: html

      <div class="sphx-glr-thumbnail-title">Solve a 2D short cantilever topology optimization problem</div>
    </div>


.. raw:: html

    </div>


.. toctree::
   :hidden:

   /examples/topology_optimization/plot_topology_optimization_L_shape
   /examples/topology_optimization/plot_topology_optimization_MBB
   /examples/topology_optimization/plot_topology_optimization_short_cantilever


.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-gallery

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download all examples in Python source code: topology_optimization_python.zip </examples/topology_optimization/topology_optimization_python.zip>`

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download all examples in Jupyter notebooks: topology_optimization_jupyter.zip </examples/topology_optimization/topology_optimization_jupyter.zip>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_