An integrated numerical model for hydride-induced failure of zirconium alloys is presented. The model solves the time-dependent and interconnected problems of temperature- and stress-directed diffusion of hydrogen, metal-hydride phase transformation, stress-directed hydride orientation, hydrogen-induced expansion and fracture within a two- dimensional finite element framework. The finite element method allows representation of arbitrary two-dimensional geometries and thermo-mechanical loading conditions. The model has a wide range of application, but is intended primarily for analyses of delayed hydride cracking in cladding tubes of light water reactor nuclear fuel rods.