Fused contactors and thermal magnetic circuit breakers are commonly applied protective devices in power distribution systems to protect the circuits when short-circuit faults occur. A power distribution system may contain various makes and models of protective devices, as a result, customizable simulation models for protective devices are demanded to effectively conduct system-level reliable analyses. To build the models, thermal energy-based data analysis methodologies are first applied to the protective devices’ physical properties, based on the manufacturer’s time/current data sheet. The models are further enhanced by integrating probability tools to simulate uncertainties in real-world application facts, for example, fortuity, variance, and failure rate. The customizable models are expected to aid the system-level reliability analysis, especially for the microgrid power systems.