The correct computation of the Stress Intensity Factor (SIF) in ultrasonic Very-High-Cycle Fatigue (VHCF) loading conditions is a key issue when investigating the crack growth rate curve with pre-cracked specimens or when evaluating critical SIF values from fracture surfaces of failed specimens. Dynamic conditions related to the resonance of the vibrating specimen, contact nonlinearity between crack faces and stress singularity at the crack tip make the SIF computation difficult and cumbersome. Generally, numerical computation through Finite Element Models (FEMs) under non-linear dynamic conditions makes use of direct integration methods (implicit or explicit). However, in the high frequency regime of ultrasonic VHCF tests, the procedure may lead to an unacceptable computational time. In order to reduce the computational time, a hybrid procedure based on the Harmonic Balance Method (HBM) and on the Virtual Crack Closure Technique (VCCT) is originally presented and applied in this paper. The dynamic field parameters calculated with the HBM are used as input data for the computation of the SIF through the VCCT.