The meshless local Petrov–Galerkin (MLPG) method is employed for anisotropic transient thermoelasticity analysis of 2D decagonal quasicrystals (QCs) subjected to transient thermal and mechanical shock loadings. The wave type model and the elasto-hydrodynamic model are applied to derive the phonon and phason governing equations, respectively. The temperature affects only the phonon field. To find the temperature distributions on the assumed 2D domain, the anisotropic heat conduction problem is solved using the MLPG method. Also, the MLPG method is successfully employed to obtain the transient behaviors of both phonon and phason displacements by solving the governing equations in local integral equations (LIEs) forms. Making use a unit step function as the test functions in the local weak-form of governing equations, we derived the local integral equations (LIEs) considered on small subdomains identical with support domains of test functions around each node. The radial basis functions are used for approximation of the spatial variation of field variables. The Laplace-transform technique is utilized to discretize the time variations.