Topology optimization with linearized buckling criteria in 250 lines of Matlab

Federico Ferrari*, Ole Sigmund, James K. Guest

*Corresponding author for this work

    Research output: Contribution to journalJournal articleResearchpeer-review

    271 Downloads (Pure)


    We present a 250-line Matlab code for topology optimization for linearized buckling criteria. The code is conceived to handle stiffness, volume and buckling load factors (BLFs) either as the objective function or as constraints. We use the Kreisselmeier-Steinhauser aggregation function in order to reduce multiple objectives (viz. constraints) to a single, differentiable one. Then, the problem is sequentially approximated by using MMA-like expansions and an OC-like scheme is tailored to update the variables. The inspection of the stress stiffness matrix leads to a vectorized implementation for its efficient construction and for the sensitivity analysis of the BLFs. This, coupled with the efficiency improvements already presented by Ferrari and Sigmund (Struct Multidiscip Optim 62:2211–2228, 2020a), cuts all the computational bottlenecks associated with setting up the buckling analysis and allows buckling topology optimization problems of an interesting size to be solved on a laptop. The efficiency and flexibility of the code are demonstrated over a few structural design examples and some ideas are given for possible extensions.
    Original languageEnglish
    JournalStructural and Multidisciplinary Optimization
    Pages (from-to)3045–3066
    Publication statusPublished - 2021


    • Topology optimization
    • Matlab
    • Buckling optimization
    • Aggregation functions
    • Optimality criteria


    Dive into the research topics of 'Topology optimization with linearized buckling criteria in 250 lines of Matlab'. Together they form a unique fingerprint.

    Cite this