As a solid material between the crystal and the amorphous, the study on quasicrystals has become an important branch of condensed matter physics. Due to the special arrangement of atoms, quasicrystals own some desirable properties, such as low friction coefficient, low adhesion, high wear resistance and low porosity. Thus, quasicrystals are expected to be applied to the coating surfaces for engines, solar cells, nuclear fuel containers and heat converters. However, when the quasicrystals are used as coating material, it is very hard to simulate the coupling fields by the finite elements numerical methods because of its thin thickness and extreme stress gradient. This is the main reason why the structure of quasicrystal coating cannot be calculated accurately and stably by various numerical platform. A general solution method which can be used to solve this contact problem for a 1D hexagonal quasicrystal coating perfectly bonded to a transversely isotropic semi-infinite substrate under the point force is presented in this paper. The solutions of the Green’s function under the distributed load can be obtained through the superposition principle. The simulation results show that this method is correct and effective, which has high calculation accuracy and fast convergence speed. The phonon–phason coupling field and elastic field in the coating and semi-infinite substrate will be derived based on the axisymmetric general solution, and the complicated coupling field of quasicrystals in coating contact space is explicitly presented in terms of elementary functions. In addition, the relationship between the coating thickness or external force and the stress component is also obtained to solve practical problems in engineering applications. The solutions presented not only bear theoretical merits, but also can serve as benchmarks to clarify various approximate methods.