The connectedness percolation threshold (ηc) and critical coordination number (Zc) of systems of penetrable spherocylinders characterized by a length polydispersity are studied by way of Monte Carlo simulations for several aspect ratio distributions. We find that (i) ηc is a nearly universal function of the weight-averaged aspect ratio, with an approximate inverse dependence that extends to aspect ratios that are well below the slender rod limit and (ii) that percolation of impenetrable spherocylinders displays a similar quasiuniversal behavior. For systems with a sufficiently high degree of polydispersity, we find that Zc can become smaller than unity, in analogy with observations reported for generalized and complex networks.PACS numbers: 64.60.ah, 61.46.Fg, 82.70.Dd Idealized elongated objects such as perfectly rigid cylinders, spherocylinders and prolate spheroids are prototypical models for a wide array of technologically relevant systems that include liquid crystals, nanocomposites based on filamentous fillers as well as fiber-reinforced materials. Percolation phenomena involving dramatic increases in, e.g., structural rigidity and electrical and thermal conductivities of composites with increasing filler loading are currently of particular interest [1]. These increases are caused by the formation of an infinite cluster of in some sense connected particles at the critical loading, i.e., the percolation threshold.It has been established by analytical [1-3] and numerical [4][5][6][7][8][9][10][11] studies that for dispersions of sufficiently elongated objects of identical size and shape, i.e, "monodisperse" objects, the geometric percolation threshold expressed in terms of the critical volume fraction of particles is inversely proportional to the aspect ratio of the filler particles. This property is exploited in the fabrication of conducting polymeric composites with very low conducting filler contents. Depending on the production processes of the composites, however, the filler particles almost invariably exhibit a pronounced polydispersity in both size and shape [12,13]. Although it represents a possible factor behind huge quantitative discrepancies between theory and experiments [14,15], such polydispersity has received relatively little attention in terms of theoretical modeling until fairly recently [1,16,17], Achievement of a theoretical understanding of how the continuum percolation of fibrous fillers is affected by polydispersity is thus key to the controlled design of a large class of composite materials for practical particle size and shape distributions.Recent analytical results obtained from integral equation methods [16] and from an heuristic mapping onto a generalized Bethe lattice [17] predict that in the slender rod limit, where the particles have asymptotically large values of the aspect ratio, the volume fraction at the percolation threshold is inversely proportional to the weight average L w = L 2 / L of the rod lengths, where the brackets imply number averages over the distribution of...