An approach for simulation of light scattering from beetles exhibiting structural colors originating from periodic helicoidal structures is presented. Slight irregularities of the periodic structure in the exoskeleton of the beetles are considered as a major cause of light scattering. Two sources of scattering are taken into account: surface roughness and volume non-uniformity. The Kirchhoff approximation is applied to simulate the effect of surface roughness. To describe volume non-uniformity, the whole structure is modeled as a set of domains distributed in space in different orientations. Each domain is modeled as an ideal uniformly twisted uniaxial medium and differs from each other by the pitch. Distributions of the domain parameters are assumed to be Gaussian. The analysis is performed using the Mueller matrix formalism which, in addition to spectral and spatial characteristics, also provides polarization properties of the scattered light.