A series of large-scale molecular dynamics simulations have been performed to investigate the shock response of nanotwinned (NT) Cu, including shock-induced plasticity, strength behind the shock front, and spall behaviors. In this study, two configurations were investigated at an impact velocity of 600 m/s, i.e., the practical NT polycrystalline Cu with an average grain size of 10 nm and the simple NT single-crystalline Cu with an impact direction of [112]. In the NT polycrystalline Cu, the average flow stress behind the shock front first increases with decreasing twin-boundary spacing (TBS), reaching a maximum at a critical TBS, and then decreases as the TBS become even smaller. This trend of the average flow stress with decreasing TBS is due to two competitive dislocation activities under shock loading, with one being inclined to the twin boundaries (the dislocation-twin boundary intersecting) and the other parallel to the twin boundaries (detwinning with twin-boundary migration). Since voids always nucleate near the grain boundary (GB) junctions and then grow along the GBs to create spallation, no apparent correlation between the spall strength and TBS is observed in the NT polycrystalline Cu. However, the spall strengths of the NT single-crystalline Cu are found to increase with decreasing TBS. Two partial dislocation slips initiated from each twin boundary create voids at the intersections between the partial dislocation slips and twin boundaries. The smaller TBSs result in a larger number of twin boundaries and provide more nucleation sites for voids, requiring a higher tensile stress to create spallation in the NT single-crystalline Cu. These findings should provide insights for understanding the deformation physics of the NT metals subjected to shock loading.