Hyperbolic materials have attracted considerable interest for their unique open hyperbolic dispersion properties. These materials support high-momentum propagating modes and strong light confinement, leading to a wide range of applications including super-resolution technologies, negative refraction and long-life propagation. Even with these wonderful optical properties, hyperbolic materials, however, cause problems when applying perfectly matched layer (PML) boundary conditions in numerical simulation software such as COMSOL Multiphysics. Due to the unfit embedded attenuation function, the built-in PML of simulation software would result in a mass of reflections in the computational domain when the background medium is hyperbolic materials. Here, we take advantage of an imaginary coordinate mapping and the complex coordinate stretching of transformation optics theory to design a PML for biaxial hyperbolic materials, which avoids any reflections and can be tuned flexibly. The proposed recipe can provide antidote and new insights for hyperbolic material studies.