International audienceAn approach to model contact-friction interactions between beams within assemblies of fibers is presented in this paper in order to simulate the mechanical behaviour of entangled structures at the scale of individual fibers using the finite element method. The determination of contact elements associating pairs of material particles is based on the detection of proximity zones between beams and on the construction of intermediate geometries approximating the actual contact zone, and allowing to consider contact along zones of non-zero lengths. The penalty method for contact is improved by adjusting the penalty parameter for each contact zone, thus stabilizing contact algorithms and allowing to handle high numbers of co...