Multimethods of simulation and experiment have been performed to investigate the interaction between glycidyl azide polymer (GAP) matrix and insensitive energetic plasticizers N-butyl-N-(2-nitroxy-ethyl)nitramine (Bu-NENA) and bis(2,2-dinitropropyl)formal/acetal (BDNPF/A). To start with, the blending energy distribution and Huggins parameters have been calculated, indicating fine miscibility between the GAP matrix and both plasticizers. The solubility parameter and binding energies show better compatibility between Bu-NENA and the GAP matrix than BDNPF/A, owing to stronger interactions. The interaction mechanism includes both hydrogen bonds and van der Waals forces. The low field NMR physical cross-link density and dynamic rheological behaviors imply larger disentanglement effect of Bu-NENA in the GAP matrix. The dynamic mechanical performance of elastomers show lower glass transition temperature of GAP/Bu-NENA blends, as supportive proof of stronger interactions between the GAP matrix and Bu-NENA in comparison with BDNPF/A.