Direct resistance check by applying advanced numerical models is getting increasingly used for the design of steel slender plated structures. This method has to take into account the same uncertainties as traditional analytical design calculations and should ensure the Eurocode-based prescribed safety level. The application of the model factor gives the possibility to account for the model-related uncertainties. The current study focuses on the determination of the model factor for one specific failure mode, the patch loading resistance. Numerical model has been developed and validated based on laboratory test results. To evaluate the model uncertainties, physically possible modeling differences are introduced, and their effects are evaluated on the resistance. The final aim of the study is to determine the model factor for the analyzed girder type and failure mode based on statistical evaluation.