Distribution of averages of Ramanujan sums

The average value of a certain normalization of Ramanujan sums is determined in terms of Bernoulli numbers and odd values of the Riemann zeta function. The distribution of values and limiting behavior of such a normalization are then studied along subsets of Beurling type integers with positive density and sequences of moduli with constraints on the number of distinct prime factors.