iitrabhi / topo-fenics

This repository contains a simple 55 lines code to carry out large scale parallel Topology Optimization using FEniCS in python

topo-fenics

This repository contains a simple 55 line code written in python for 2D and 3D to carry out large-scale parallel Topology Optimization using the open-source package FEniCS.

Abstract

This paper presents a 55-line code written in python for 2D and 3D topology optimization based on the open-source finite element computing software (FEniCS), equipped with various finite element tools and solvers. PETSc is used as the linear algebra back-end, which results in significantly less computational time than standard python libraries. The code is designed based on the popular solid isotropic material with penalization (SIMP) methodology. Extensions to multiple load cases, different boundary conditions, and incorporation of passive elements are also presented. Thus, this implementation is the most compact implementation of SIMP based topology optimization for 3D as well as 2D problems.

Detailed explanation could be found here.

Citation

If this code is useful for your research, please cite our publication.

@article{gupta55lineCodeLargescale2020,
title = {A 55-Line Code for Large-Scale Parallel Topology Optimization in {{2D}} and {{3D}}},
author = {Gupta, Abhinav and Chowdhury, Rajib and Chakrabarti, Anupam and Rabczuk, Timon},
year = {2020},
doi = {https://doi.org/10.48550/arXiv.2012.08208},
journal = {arXiv preprint arXiv:2012.08208},
}

topo-fenics's Issues

Migrating to Fenics-x.

This repo runs on dolfin, however, many structural optimizations run on dolfinx.
So, we need to update this repository, to keep up with tech fashion.

Implementations for 3D and Parallel

Hi,

I just have come across your paper on "A 55-line code for large-scale parallel topology optimization in 2D and 3D".
It is really great to see this work for me, because I am a new student interested in this field.

In this github, it seems only the 2D implementation is uploaded.
So I am wondering if you plan to upload the 3D bridge problem in Figure 10 as well as the parallel implementation things using mpirun. If you upload those implementations, it definitely helpful for many students including me who really would like to jump into this field.

Thank you.