Hybrid simulation is a method used to investigate the dynamic response of a system subjected to a realistic loading scenario. The system under consideration is divided into multiple individual loading-rate-sensitive substructures, out of which one or more are tested physically, whereas the remaining are simulated numerically. The coupling of all substructures forms the so-called hybrid model. Although hybrid simulation has been extensively used across various engineering disciplines, it is often the case that the hybrid model and related excitation is conceived as deterministic. However, associated uncertainties are present whilst simulation deviation due to their presence could be significant. To this regard, global sensitivity analysis based on Sobol' indices can be used to determine the sensitivity of the hybrid model response due to the presence of the associated uncertainties. Nonetheless, estimation of the Sobol' sensitivity indices requires unaffordable amount of hybrid simulation evaluations. Therefore, surrogate modeling techniques are used to alleviate this burden. In this paper, three different surrogate modeling methods are examined, namely polynomial chaos expansion, Kriging and polynomial chaos Kriging. A case study encompassing a virtual hybrid model is employed and hybrid model response quantities of interest are selected. Their respective surrogates are developed using all three aforementioned techniques. The Sobol' indices obtained utilizing each examined surrogate are compared with each other and results highlight potential deviations.