In a series of previous works, we established a novel approach to topology optimization for compliance minimization based on thermodynamic principles known from the field of material modeling. Hamilton's principle for dissipative processes directly yields a partial differential equation (referred to as the evolution equation) as an update scheme for the spatial distribution of density mass describing the topology. Consequently, no additional mathematical minimization algorithms are needed. In this article, we introduce a regularization scheme by penalization of the gradient of the spatial distribution of mass density. The parabolic evolution equation (owing to a similar structure to the transient heat-conduction equation) is solved most efficiently by an explicit time discretization. The Laplace operator is discretized via a Taylor series expansion yielding an operator matrix that is constant for the entire optimization process. This method shares some similarities to meshless methods and allows for an accurate application also on unstructured finite element meshes. The minimal size of the structure member can directly be controlled, a priori, by a numerical parameter introduced along with the regularization, similar to classical filter radii.
KEYWORDSmeshless Laplacian, parabolic PDE, regularization, thermodynamic topology optimization, unstructured meshes, variational material modeling
INTRODUCTIONOver recent decades, topology optimization has become a major field of interest in mechanical engineering. Approaches have been developed for different optimization objectives, constraints relating to physically feasible or manufacturable designs, or models for multiphysics problems. The most native approach is considered to be topology optimization with compliance minimization under a structure volume constraint for linear elastic materials, see, eg, the work of Bendsøe and Sigmund, 1 in which the elastic strain energy is subject to minimization and is referred to as (structural) compliance. Numerous solution approaches were developed of which the most popular are overviewed and compared in the work of Sigmund and Maute. 2 Usually, the topology has to be found within a given design space under mechanical boundary conditions, for which the mechanical problem is solved by a finite element approach. The design space contains material and void phases to identify the topology.This can be achieved either by tracing the phase borders or by introducing a density field that describes the material configuration in each point of the design space. The phase borders are subject to optimization in level-set approaches, 3,4 which Int J Numer Methods Eng. 2019;117:991-1017.wileyonlinelibrary.com/journal/nme