A method of hydraulic fracturing a subsurface formation, the method comprising: using a three-dimensional finite element model implemented on one or more processors to simulate a deviated well with a wellbore casing (122, 314), a cement (130, 316) adjacent to the wellbore casing (122, 314), and a perforation cluster (140, 400) with at least two perforations (138, 402) through a side wall of the wellbore casing (122, 314), each of the at least two perforations (138, 402) extending through a rock of the subsurface formation via a respective perforation tunnel and having a different phase angle relative to a longitudinal axis of the wellbore casing (122, 314); solving the three-dimensional finite element model implemented on one or more processors to determine a breakdown pressure of the deviated well based on an amount of tensile damage of the perforation cluster (140, 400) induced by an applied pressure representing injected hydraulic fluid; drilling and completing a deviated well with a wellbore casing size, tubing size, and wellhead selected at least in part based on the determined breakdown pressure; and injecting hydraulic fluid into the deviated well at an injection pressure determined at least in part on the determined breakdown pressure to cause hydraulic fracturing of the rock of the subsurface formation, wherein the three-dimensional finite element model includes a concrete damage plasticity model representing constitutive behavior of the rock and a contact relationship between an interface of the wellbore casing (122, 314) and the cement (130, 316) adjacent to the wellbore casing (122, 314), the contact relationship modeling de-bonding of the interface of the wellbore casing (122, 314)and the cement (130, 316). using a three-dimensional finite element model implemented on one or more processors to simulate a deviated well with a wellbore casing (122, 314), a cement (130, 316) adjacent to the wellbore casing (122, 314), and a perforation cluster (140, 400) with at least two perforations (138, 402) through a side wall of the wellbore casing (122, 314), each of the at least two perforations (138, 402) extending through a rock of the subsurface formation via a respective perforation tunnel and having a different phase angle relative to a longitudinal axis of the wellbore casing (122, 314); solving the three-dimensional finite element model implemented on one or more processors to determine a breakdown pressure of the deviated well based on an amount of tensile damage of the perforation cluster (140, 400) induced by an applied pressure representing injected hydraulic fluid; drilling and completing a deviated well with a wellbore casing size, tubing size, and wellhead selected at least in part based on the determined breakdown pressure; and injecting hydraulic fluid into the deviated well at an injection pressure determined at least in part on the determined breakdown pressure to cause hydraulic fracturing of the rock of the subsurface formation, wherein the three-dimensional finite element model includes a concrete damage plasticity model representing constitutive behavior of the rock and a contact relationship between an interface of the wellbore casing (122, 314) and the cement (130, 316) adjacent to the wellbore casing (122, 314), the contact relationship modeling de-bonding of the interface of the wellbore casing (122, 314)and the cement (130, 316).