This report is a part of the larger project of non-linear dynamics investigation of three coupled physical pendulums with damping and with arbitrary situated barriers, and externally driven. The set of differential equations and the set of algebraic inequalities (representing a barrier) governing the motion of three coupled rods are presented in the non-dimensional form. The system of governing equations is integrated between two successive impacts, and the discontinuity points are detected (by halving time step until a required precision is obtained). In each impact time, the state of the system is transformed using the extended restitution coefficient rule. The theory of Aizerman and Gantmakher is used to calculate the fundamental solution matrices in the analyzed system exhibiting discontinuities. The fundamental matrices are used during calculation of Lyapunov exponents, during stability analysis of periodic solutions (Floquet multipliers) and in shooting method applied to detect and trace periodic orbits. Some examples for three coupled identical rods with horizontal barrier are reported.