The paper presents a novel approach of constructing the state matrix of the multi-machine power system for SSR (subsynchronous resonance) analysis using the linearized equations of individual devices including electrical transmission network dynamics. The machine models in the local d-q reference frame are integrated with the network models in the common R-I reference frame by simply transforming their output equations into the R-I frame where the transformed output is used as the input to the network dynamics or vice versa. The salient feature of the formulation is that it allows for modular construction of various component models without rearranging the overall state space formulation. This universal SSR small signal stability program provides a flexible tool for systematic analyses of SSR small-signal stability impacts of both conventional devices such as generation systems and novel devices such as power electronic apparatus and their controllers. The paper also presents its application results to IEEE benchmark models.