ABSTRACT This work presents a novel approach for the solution of multi-material topology optimization problems of elastic structures coupling the Direct Meshless Local Petrov–Galerkin (DMLPG) method with the Gradual Bi-direction Evolutionary Structural Optimization (gBESO) method. In recent years, meshless methods have been used in several engineering fields, showing some advantages compared to mesh-based methods. DMLPG is a truly meshless method, which achieves results with good accuracy and computational efficiency, avoiding the numerical integration of complicated shape functions by adopting low-degree polynomials. In the proposed algorithm, DMLPG is used to obtain smooth nodal displacements, strains and stresses, and gBESO updates the structural geometry based on design sensitivity values of the compliance while it gradually increases the Young's modulus of the material. The optimization problem minimizes the compliance objective function subject to volume constraints. Numerical examples demonstrate the applicability and validity of the technique.