Perturbation method for solving elastic three-dimensional (3D) problems for 3D icosahedral quasicrystals is proposed. Considering an infinite 3D icosahedral quasicrystal weakened by a circular crack, we obtain the uniformly valid asymptotic solutions up to O(R2) for the mode I loading, where R is the elastic constant of phonon-phason coupling.