We review some features of topology optimization with a lower bound on the critical load factor, as computed by linearized buckling analysis. The change of the optimized design, the competition between stiffness and stability requirements and the activation of several buckling modes, depending on the value of such lower bound, are studied. We also discuss some specific issues which are of particular interest for this problem, as the use of non-conforming finite elements for the analysis, the use of inconsistent sensitivities in the optimization and the replacement of the single eigenvalue constraints with an aggregated measure. We discuss the influence of these practices on the optimization result, giving some recommendations.

Compact and efficient Matlab implementations of compliance Topology Optimization (TO) for 2D and 3D continua are given, consisting of 99 and 125 lines respectively. On discretizations ranging from 3•10 4 to 4.8 • 10 5 elements, the 2D version, named top99neo, shows speedups from 2.55 to 5.5 times compared to the well-known top88 code (Andreassen et al, 2011). The 3D version, named top3D125, is the most compact and efficient Matlab implementation for 3D TO to date, showing a speedup of 1.9 times compared to the code of Amir et al ( 2014), on a discretization with 2.2 • 10 5 elements. For both codes, improvements are due to much more efficient procedures for the assembly and implementation of filters and shortcuts in the design update step. The use of an acceleration strategy, yielding major cuts in the overall computational time, is also discussed, stressing its easy integration within the basic codes.

No abstract

This work presents a multilevel approach to large-scale topology optimization accounting for linearized buckling criteria. The method relies on the use of preconditioned iterative solvers for all the systems involved in the linear buckling and sensitivity analyses and on the approximation of buckling modes from a coarse discretization. The strategy shows three main benefits: first, the computational cost for the eigenvalue analyses is drastically cut. Second, artifacts due to local stress concentrations are alleviated when computing modes on the coarse scale. Third, the ability to select a reduced set of important global modes and filter out less important local ones. As a result, designs with improved buckling resistance can be generated with a computational cost little more than that of a corresponding compliance minimization problem solved for multiple loading cases. Examples of 2D and 3D structures, discretized by up to some millions of degrees of freedom in Matlab, are presented to show the effectiveness of the proposed method. Finally, a post-processing procedure is suggested in order to reinforce the optimized design against local buckling.

Summary The article presents an efficient solution method for structural topology optimization aimed at maximizing the fundamental frequency of vibration. Nowadays, this is still a challenging problem mainly because of the high computational cost required by spectral analyses. The proposed method relies on replacing the eigenvalue problem with a frequency response one, which can be tuned and efficiently solved by a multilevel procedure. Connections of the method with multigrid eigenvalue solvers are discussed in details. Several applications demonstrating more than 90% savings of the computational time are presented as well.