This study focuses on a virtual synchronous machine (VSM) based on voltage source converters to mimic the behavior of synchronous machines (SMs) and improve the damping ratio of the power system. The VSM model is simplified according to some assumptions (neglecting the speed variation and the stator transients) to allow for the possibility of dealing with low-frequency oscillation in large-scale systems with many VSMs. Furthermore, a virtual power system stabilizer (VPSS) structure is proposed and tuned using a method based on a linearized power system dynamic model. The linear and nonlinear analyses examine the stability of two modified versions of a 16-machine power system in which, in the first case, partial classical machines are replaced by VSMs, and in the second case, all SMs are replaced by VSMs. The simulation results of the case studies validate the efficiency of the proposed control strategy. Index Terms-Synchronous machine (SM), small-signal stability , transient stability, virtual power system stabilizer (VPSS), virtual synchronous machine (VSM). NOMENCLATURE A. Operators Δ Perturbation-variation operator p Differential operator d/dt s Laplace operator ° Equilibrium script T Transposition script B. Variables δ Angular position of rotor with respect to synchronous rotating frame ω, ω s Angular speed and synchronous speed ψ abc =[ψ a , ψ b , ψ c ] T Vector of abc instantaneous stator flux ψ dq0 = [ψ d , ψ q , ψ 0 ] T Vector of dq0 instantaneous stator flux ζ Damping ratio 1/K Integrator gain A, B, C, and D Matrices of state, input, output, and feed-forward E q Synchronous internal voltage mimicked by virtual synchronous machine (VSM) E B Infinite bus voltage e abc =[e a , e b , e c ] T Vector of abc instantaneous stator voltage e dq0 =[e d , e q , e 0 ] T Vector of dq0 instantaneous stator voltage e fd , i fd , R fd , ψ fd Field voltage, current, resistance, and flux H, D p Inertia constant and damping constant K A , T A Time constant and gain of automatic voltage regulator (AVR) K VPSS Gain of virture power system stabilizer (VPSS) i abc =[i a , i b , i c ] T Vector of abc instantaneous stator current i dq0 =[i d , i q , i 0 ] T Vector of dq0 instantaneous stator current L aa0 , L ab0 Proper and mutual inductances of stator L abc Matrix of abc stator inductance L abc-r (θ) Vector of abc stator-rotor mutual induc-tance L ad , L aq d-and q-axis mutual inductances owing to flux that links rotor circuits L dq0 Matrix of dq0 stator inductance L dq0-r (θ), L r-dq0 (θ) Vectors of dq0 stator-rotor and rotor-stator mutual inductance L d , L q d-and q-axis stator self-inductances L ffd , L afd Proper inductance of field circuit and mutual inductance of stator-field circuit L l Leakage inductance owing to flux that does not link any rotor circuit P θ Park transformation matrix from abc variables to dq0 variables that keeps the same Manuscript