Topology optimization of load-bearing capacity

The present work addresses the problem of maximizing a structure load-bearing capacity subject to given material strength properties and a material volume constraint. This problem can be viewed as an extension to limit analysis problems which consist in finding the maximum load capacity for a fixed geometry. We show that it is also closely linked to the problem of minimizing the total volume under the constraint of carrying a fixed loading. Formulating these topology optimization problems using a continuous field representing a fictitious material density yields convex optimization problems which can be solved efficiently using state-of-the-art solvers used for limit analysis problems. We further analyze these problems by discussing the choice of the material strength criterion, especially when considering materials with asymmetric tensile/compressive strengths. In particular, we advocate the use of a L1-Rankine criterion which tends to promote uniaxial stress fields as in truss-like structures. We show that the considered problem is equivalent to a constrained Michell truss problem. Finally, following the idea of the SIMP method, the obtained continuous topology is post-processed by an iterative procedure penalizing intermediate densities. Benchmark examples are first considered to illustrate the method overall efficiency while final examples focus more particularly on no-tension materials, illustrating how the method is able to reproduce known structural patterns of masonry-like structures. This paper is accompanied by a Python package based on the FEniCS finite-element software library.