This paper investigates commutation failure (CF) in line commutated HVdc converters within multi-infeed HVDC systems and thereby provides valuable insights for the operation and design of HVdc systems. CF susceptibility across various operational scenarios, fault types, and ac system representations at different frequencies is explored. An Electromagnetic Transients (EMT) model is constructed and parameter changes are applied using Monte-Carlo Simulation. The process is fully automated with the use of a controlling program written in Python, which varies parameter values and facilitates result retrieval and post-simulation analysis. A typical analysis results in several hundred thousands of simulation runs, requiring the procedure to be implemented on a parallel computing platform. The results show that two-phase faults are the most critical CF triggers, an aspect often overlooked in previous literature focused on three-phase or single-phase-to-ground faults. Furthermore, the frequency dependent characteristics of the ac network can result in very different CF susceptibilities.