In order to study the damage of polymer bonded explosive（PBX） under dynamic shock loads， the constitutive curve of the mechanical damage of PBX under different strain rates was obtained by using the Split-Hopkinson pressure bar（SHPB） test device. And the constitutive parameters in the Z-W-T constitutive model with damage variables were fitted in sections. Then， based on the fitting results， the finite element theory， elastoplastic mechanics and ABAQUS/VUMAT， the subroutine of PBX containing damage is completed， and is verified by the finite element simulation. The model is used for finite element engineering examples. Results show that the verification results of finite element samples are in good agreement with the experimental results， and the correlation degree of the results is higher than 0.95.The stress cloud diagram and damage variable cloud diagram for fragment penetration example reflect the change and possible damage location of PBX under shock loads.The evolution process of damage cloud diagram reflects the correlation between damage evolution and strain rate effect of PBX.