We consider a well-studied generalization of the maximum clique problem which is defined as follows. Given a graph G on n vertices and a fixed parameter $$d\ge 1$$ , in the maximum diameter-bounded subgraph problem (MaxDBS for short) the goal is to find a (vertex) maximum subgraph of G of diameter at most d. For $$d=1$$ , this problem is equivalent to the maximum clique problem and thus it is NP-hard to approximate it within a factor $$n^{1-\epsilon }$$ , for any $$\epsilon >0$$ . Moreover, it is known that, for any $$d\ge 2$$ , it is NP-hard to approximate MaxDBS within a factor $$n^{1/2-\epsilon }$$ , for any $$\epsilon >0$$ . In this paper we focus on MaxDBS for the class of unit disk graphs. We provide a polynomial-time constant-factor approximation algorithm for the problem. The approximation ratio of our algorithm does not depend on the diameter d. Even though the algorithm itself is simple, its analysis is rather involved. We combine tools from the theory of hypergraphs with bounded VC-dimension, k-quasi planar graphs, fractional Helly theorems, and several geometric properties of unit disk graphs.