This paper presents an alternative topology optimization method based on an efficient meshless smoothed particle hydrodynamics (SPH) algorithm. To currently calculate the objective compliance, the deficiencies in standard SPH method are eliminated by introducing corrective smoothed particle method and total Lagrangian formulation. The compliance is established relative to a designed density variable at each SPH particle which is updated by optimality criteria method. Topology optimization is realized by minimizing the compliance using a modified solid isotropic material with penalization approach. Some numerical examples of plane elastic structure are carried out and the results demonstrate the suitability and effectiveness of the proposed SPH method in the topology optimization problem.
727As suggested by Sigmund and Maute [4], non-FE based topology optimization methods should be explored to overcome the difficulties met in FE method. At present, meshless method without mesh connectivity has drawn much attention in various real-world engineering problems [13][14][15] and been developed increasingly for the structural topology optimization due to the simplicity and accuracy. As far as we know, the first research on topology optimization problem was performed by Kim [16] for two-dimensional continuum structures using element-free Galerkin (EFG) method; however, the checker-board phenomenon was not eliminated. EFG is also developed to solve the topology optimization problem undergoing geometrically nonlinear deformations, where the material distribution is described by a independent point-wise density interpolation method [17]. Meshless reproducing kernel method is successfully developed for geometrically nonlinear plane elasticity problems [18]. Li [19] has studied topology-optimization problems of elastic structures based on meshless local Petrov-Galerkin"mixed collocation" method. Luo [20] proposed a new structural topology optimization method using a dual-level point-wise density approximant and the meshless Galerkin weak-forms.Unlike the aforementioned meshless method with weak-forms, smoothed particle hydrodynamics (SPH) method discretizes directly the balance equation in strong form with collocation framework. Moreover, the background virtual cells independent of field nodes are required for Guass quadrature in EFG [16,17,20], reproducing kernel [18], and meshless local Petrov-Galerkin [19], but unnecessary in SPH [21]. These make SPH easier to be understood, simpler to be programmed and faster to be implemented, which contributes to the wild application of SPH in the region of computational fluid mechanics. However, the inconsistency and instability are the intrinsic problems in classical SPH method and limit the application in solid mechanics [21].Firstly, because of the boundary truncation of the support domain, SPH method cannot even ensure C 0 consistency, which decreases the approximation accuracy. Many researchers have devotes a lot to the promotion of the consistency and proposed corrective smoothe...