Objective: Most of effectiveness assessments of the widely-used Massage therapy were based on subjective routine clinical assessment tools, such as Visual Analogue Scale (VAS) score. However, few studies demonstrated the impact of massage on the Electroencephalograph (EEG) rhythm decoding of Motor imagery (MI) and motion execution (ME) with trunk left/right bending in patients with skeletal muscle pain. Method: We used the sample entropy (SampEn), permutation entropy (PermuEn), common spatial pattern (CSP) features, support vector machine (SVM) and logic regression (LR) classifiers. We also used the convolutional neural network (CNN) and attention-based bi-directional long short-term memory (BiLSTM) for classification. Results: The averaged SampEn and PermuEn values of alpha rhythm decreased in almost fourteen channels for five statuses (quiet, MI with left/right bending, ME with left/right bending). It indicated that massage alleviates the pain for the patients of skeletal pain. Furthermore, compared with the SVM and LR classifiers, the BiLSTM method achieved a better area under curve (AUC) of 0.89 for the classification of MI with trunk left/right bending before massage. The AUC became smaller after massage than that before massage for the classification of MI with trunk left/right bending using CNN and BiLSTM methods. The Permutation direct indicator (PDI) score showed the significant difference for patients in different statuses (before vs after massage, and MI vs ME). Conclusions: Massage not only affects the quiet status, but also affects the MI and ME. Clinical Impact: Massage therapy may affect a bit on the accuracy of MI with trunk left/right bending and it change the topography of MI and ME with trunk left/right bending for the patients with skeletal muscle pain.