This article proposes the distribution similarity measure–based Gaussian mixtures model for the contact-state (CS) modelling in force-guided robotic assembly processes of flexible rubber parts. The wrench (Cartesian force and torque) signals of the manipulated object are captured for different states of the given assembly process. The distribution similarity measure–based Gaussian mixtures model CS modelling scheme is employed in modelling the captured wrench signals for different CSs. The proposed distribution similarity measure–based Gaussian mixtures model CS modelling scheme uses the Gaussian mixtures model in modelling the captured signals. The parameters of the Gaussian mixtures models are computed using expectation maximisation. The optimal number of Gaussian mixtures model components for each CS model is determined by considering the classification success rate as an index for the similarity measure between the distribution of the captured signals and the developed models. The optimal number of Gaussian mixtures model components corresponds to the highest classification success rate; hence, object elasticity variation would be accommodated by properly choosing the optimal number of Gaussian mixtures model components. The performance of the proposed distribution similarity measure–based Gaussian mixtures model CS modelling strategy is evaluated by a test stand composed of a KUKA lightweight robot doing peg-in-hole assembly processes for flexible rubber objects. Two rubber objects with different elasticity are considered for two experiments; in the first experiment, an elastic peg of 30 Shore A hardness is considered and that of the second experiment has hardness of 6 Shore A which is even softer than the one used in experiment 1. Employing the proposed distribution similarity measure–based Gaussian mixtures model CS modelling strategy excellent classification success rate was obtained for both experiments. However, more Gaussian mixtures model components are required for the softer one that gives a strong impression of the non-stationarity behaviour increment for softer materials. Comparison is performed with the available CS modelling schemes and the distribution similarity measure–based Gaussian mixtures model is shown to provide the best classification success rate performance with a reduced computational time.