Summary
A bidirectional evolutionary structural optimization algorithm is presented, which employs integer linear programming to compute optimal solutions to topology optimization problems with the objective of mass minimization. The objective and constraint functions are linearized using Taylor's first‐order approximation, thereby allowing the method to handle all types of constraints without using Lagrange multipliers or sensitivity thresholds. A relaxation of the constraint targets is performed such that only small changes in topology are allowed during a single update, thus ensuring the existence of feasible solutions. A variety of problems are solved, demonstrating the ability of the method to easily handle a number of structural constraints, including compliance, stress, buckling, frequency, and displacement. This is followed by an example with multiple structural constraints and, finally, the method is demonstrated on a wing‐box, showing that topology optimization for mass minimization of real‐world structures can be considered using the proposed methodology.