The elastic strain energy potential for nonlinear fibre-reinforced materials is customarily obtained by superposition of the potentials of the matrix and of each family of fibres. Composites with statistically oriented fibres, such as biological tissues, can be seen as being reinforced by a continuous infinity of fibre families, the orientation of which can be represented by means of a probability density function defined on the unit sphere (i.e. the solid angle). In this case, the superposition procedure gives rise to an integral form of the elastic potential such that the deformation features in the integral, which therefore cannot be calculated a priori. As a consequence, an analytical use of this potential is impossible. In this paper, we implemented this integral form of the elastic potential into a numerical procedure that evaluates the potential, the stress and the elasticity tensor at each deformation step. The numerical integration over the unit sphere is performed by means of the method of spherical designs, in which the result of the integral is approximated by a suitable sum over a discrete subset of the unit sphere. As an example of application, we modelled the collagen fibre distribution in articular cartilage, and used it in simulating displacement-controlled tests: the unconfined compression of a cylindrical sample and the contact problem in the hip joint.