PyTopo3D: A Python Framework for 3D SIMP-based Topology Optimization

2025-04-08

AI Summary

The Python scientific computing ecosystem lacks accessible, open-source tools for three-dimensional topology optimization (TO). Method: This paper introduces the first end-to-end, engineering-oriented 3D TO framework implemented entirely in Python. It combines the Solid Isotropic Material with Penalization (SIMP) method and the Optimality Criteria (OC) algorithm, and integrates STL geometry import/export, real-time 3D visualization, KD-tree-accelerated sensitivity filtering, and efficient sparse matrix operations. Optional parallel linear solvers (e.g., PETSc) are supported to enable complex geometry modeling and manufacturable STL output. Contribution/Results: To our knowledge, this is the first pure-Python implementation achieving millisecond-scale sensitivity updates, interactive convergence monitoring, and a fully automated 3D TO workflow. It significantly improves computational efficiency and usability while filling a critical gap in open-source Python-based 3D topology optimization tools.

Technology Category

• Search and Optimization: Sampling/Simulation-based Search
• Constraint Satisfaction and Optimization: Solvers and Tools
• Reasoning under Uncertainty: Stochastic Optimization

Application Category

• Systems and Infrastructure for Web, Mobile and WoT: Experiences and lessons learnt from Web-based algorithms and system deployments
• Graph Algorithms and Modeling for the Web: Efficient manipulation of static and dynamic Web-related graphs
• Responsible Web: Human-perceived consequences of algorithmic deployment on the web

Abstract

Three-dimensional topology optimization (TO) is a powerful technique in engineering design, but readily usable, open-source implementations remain limited within the popular Python scientific environment. This paper introduces PyTopo3D, a software framework developed to address this gap. PyTopo3D provides a feature-rich tool for 3D TO by implementing the well-established Solid Isotropic Material with Penalization (SIMP) method and an Optimality Criteria (OC) update scheme, adapted and significantly enhanced from the efficient MATLAB code by Liu and Tovar (2014). While building on proven methodology, PyTopo3D's primary contribution is its integration and extension within Python, leveraging sparse matrix operations, optional parallel solvers, and accelerated KD-Tree sensitivity filtering for performance. Crucially, it incorporates functionalities vital for practical engineering workflows, including the direct import of complex design domains and non-design obstacles via STL files, integrated 3D visualization of the optimization process, and direct STL export of optimized geometries for manufacturing or further analysis. PyTopo3D is presented as an accessible, performance-aware tool and citable reference designed to empower engineers, students, and researchers to more easily utilize 3D TO within their existing Python-based workflows.

Problem

• Lack of open-source Python tools for 3D topology optimization
• Need for efficient SIMP-based 3D TO with practical engineering features
• Limited integration of 3D TO workflows in Python environments

Innovation

• Python-based 3D SIMP topology optimization framework
• Sparse matrix operations and parallel solvers
• STL import/export and 3D visualization