This paper presents a new algorithm to construct the modular network model for SSR analysis by simply applying KCL to each node and KVL to all branches connected to the node sequentially. This method has advantages that the model can be derived directly from the system data for transient stability study and turbine/generator shaft model, the resulted model in the form of augmented state matrix is very sparse, and thus efficient SSR study of a large scale system becomes possible. The proposed algorithm is verified with the IEEE First and Second Benchmark models.