.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "examples/topology_optimization/topology_optimization_short_cantilever.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note Click :ref:`here ` to download the full example code .. rst-class:: sphx-glr-example-title .. _sphx_glr_examples_topology_optimization_topology_optimization_short_cantilever.py: Solve a 2D short cantilever topology optimization problem ========================================================= .. GENERATED FROM PYTHON SOURCE LINES 23-34 .. code-block:: default from __future__ import annotations import matplotlib.pyplot as plt from gemseo.api import configure_logger from gemseo.api import create_scenario from gemseo.problems.topo_opt.topopt_initialize import ( initialize_design_space_and_discipline_to, ) from matplotlib import colors configure_logger() .. rst-class:: sphx-glr-script-out .. code-block:: none .. GENERATED FROM PYTHON SOURCE LINES 35-38 Setup the topology optimization problem --------------------------------------- Define the target volume fraction: .. GENERATED FROM PYTHON SOURCE LINES 38-40 .. code-block:: default volume_fraction = 0.3 .. GENERATED FROM PYTHON SOURCE LINES 41-42 Define the problem type: .. GENERATED FROM PYTHON SOURCE LINES 42-44 .. code-block:: default problem_name = "Short_Cantilever" .. GENERATED FROM PYTHON SOURCE LINES 45-46 Define the number of elements in the x- and y- directions: .. GENERATED FROM PYTHON SOURCE LINES 46-49 .. code-block:: default n_x = 50 n_y = 25 .. GENERATED FROM PYTHON SOURCE LINES 50-51 Define the full material Young's modulus and Poisson's ratio: .. GENERATED FROM PYTHON SOURCE LINES 51-54 .. code-block:: default e0 = 1 nu = 0.3 .. GENERATED FROM PYTHON SOURCE LINES 55-56 Define the penalty of the SIMP approach: .. GENERATED FROM PYTHON SOURCE LINES 56-58 .. code-block:: default penalty = 3 .. GENERATED FROM PYTHON SOURCE LINES 59-60 Define the minimum member size in the solution: .. GENERATED FROM PYTHON SOURCE LINES 60-62 .. code-block:: default min_memeber_size = 1.5 .. GENERATED FROM PYTHON SOURCE LINES 63-64 Instantiate the :class:`.DesignSpace` and the disciplines: .. GENERATED FROM PYTHON SOURCE LINES 64-76 .. code-block:: default design_space, disciplines = initialize_design_space_and_discipline_to( problem=problem_name, n_x=n_x, n_y=n_y, e0=e0, nu=nu, penalty=penalty, min_member_size=min_memeber_size, vf0=volume_fraction, ) .. GENERATED FROM PYTHON SOURCE LINES 77-80 Solve the topology optimization problem --------------------------------------- Generate a :class:`.MDOScenario`: .. GENERATED FROM PYTHON SOURCE LINES 80-87 .. code-block:: default scenario = create_scenario( disciplines, formulation="DisciplinaryOpt", objective_name="compliance", design_space=design_space, ) .. GENERATED FROM PYTHON SOURCE LINES 88-89 Add the volume fraction constraint to the scenario: .. GENERATED FROM PYTHON SOURCE LINES 89-91 .. code-block:: default scenario.add_constraint("volume fraction", "ineq", value=volume_fraction) .. GENERATED FROM PYTHON SOURCE LINES 92-93 Generate the XDSM: .. GENERATED FROM PYTHON SOURCE LINES 93-95 .. code-block:: default scenario.xdsmize() .. rst-class:: sphx-glr-script-out .. code-block:: none INFO - 16:56:14: Generating HTML XDSM file in : xdsm.html .. GENERATED FROM PYTHON SOURCE LINES 96-97 Execute the scenario: .. GENERATED FROM PYTHON SOURCE LINES 97-99 .. code-block:: default scenario.execute(input_data={"max_iter": 200, "algo": "NLOPT_MMA"}) .. rst-class:: sphx-glr-script-out .. code-block:: none INFO - 16:56:14: INFO - 16:56:14: *** Start MDOScenario execution *** INFO - 16:56:14: MDOScenario INFO - 16:56:14: Disciplines: DensityFilter FininiteElementAnalysis MaterialModelInterpolation VolumeFraction INFO - 16:56:14: MDO formulation: DisciplinaryOpt INFO - 16:56:14: Optimization problem: INFO - 16:56:14: minimize compliance(x) INFO - 16:56:14: with respect to x INFO - 16:56:14: subject to constraints: INFO - 16:56:14: volume fraction(x) <= 0.3 INFO - 16:56:14: Solving optimization problem with algorithm NLOPT_MMA: INFO - 16:56:14: ... 0%| | 0/200 [00:00 .. GENERATED FROM PYTHON SOURCE LINES 112-113 .. image:: /_images/topology_optimization/Short_Cantilever_history.png .. GENERATED FROM PYTHON SOURCE LINES 116-117 Plot the solution: .. GENERATED FROM PYTHON SOURCE LINES 117-130 .. code-block:: default plt.ion() # Ensure that redrawing is possible fig, ax = plt.subplots() im = ax.imshow( -scenario.optimization_result.x_opt.reshape((n_x, n_y)).T, cmap="gray", interpolation="none", norm=colors.Normalize(vmin=-1, vmax=0), ) fig.show() im.set_array(-scenario.optimization_result.x_opt.reshape((n_x, n_y)).T) fig.canvas.draw() plt.savefig(problem_name + "_solution.png") .. image-sg:: /examples/topology_optimization/images/sphx_glr_topology_optimization_short_cantilever_001.png :alt: topology optimization short cantilever :srcset: /examples/topology_optimization/images/sphx_glr_topology_optimization_short_cantilever_001.png :class: sphx-glr-single-img .. GENERATED FROM PYTHON SOURCE LINES 131-132 .. image:: /_images/topology_optimization/Short_Cantilever_solution.png .. rst-class:: sphx-glr-timing **Total running time of the script:** ( 0 minutes 18.317 seconds) .. _sphx_glr_download_examples_topology_optimization_topology_optimization_short_cantilever.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: topology_optimization_short_cantilever.py ` .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: topology_optimization_short_cantilever.ipynb ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_