An elliptical crack subjected to coupled phonon–phason loadings in a three-dimensional body of two-dimensional hexagonal quasicrystals is analytically investigated. Owing to the existence of the crack, the phonon and phason displacements are discontinuous along the crack face. The phonon and phason displacement discontinuities serve as the unknown variables in the generalized potential function method which are used to derive the boundary integral equations. These boundary integral equations governing Mode I, II, and III crack problems in two-dimensional hexagonal quasicrystals are expressed in integral differential form and hypersingular integral form, respectively. Closed-form exact solutions to the elliptical crack problems are first derived for two-dimensional hexagonal quasicrystals. The corresponding fracture parameters, including displacement discontinuities along the crack face and stress intensity factors, are presented considering all three crack cases of Modes I, II, and III. Analytical solutions for a penny-shaped crack, as a special case of the elliptical problem, are given. The obtained analytical solutions are graphically presented and numerically verified by the extended displacement discontinuities boundary element method.