Summary: Based on rigorous operator theory, a general solution of the three-dimensional elasticity equations of equilibrium for 1D hexagonal quasicrystals is obtained. The solution is expressed in terms of four quasiharmonic functions, which is very simple and useful. The point phonon (phason) force solution of an infinite 1D hexagonal quasicrystal body is derived, all in terms of elementary functions. They can play an important role in numerical simulations such as boundary element method.