We develop an Fe-C-H interatomic potential based on the modified embedded-atom method (MEAM) formalism based on density functional theory to enable large-scale modular dynamics simulations of carbon steel and hydrogen.