Understanding the dynamic behavior of soil layers under seismic loading is pivotal for accurate seismic design and risk assessment. This study conducts a two-dimensional equivalent linear site response analysis using the Scaled Boundary Finite Element Method with Rayleigh damping to enhanced modeling accuracy. SBFEM combines the advantages of finite and boundary element methods, offering high efficiency in simulating wave propagation and stress concentrations in semi-infinite domains. A MATLAB i