Strain engineering is a powerful strategy for tuning the optical, electrical, vibrational properties of 2D nanomaterials. In this work, a four‐point bending apparatus is constructed to apply both compressive and tensile strain on 2D anisotropic black phosphorus flake. Further polarized Raman spectroscopy is used to study the vibrational modes of black phosphorus flakes under uniaxial strain applied along various crystalline orientations. Here, a strong anisotropic blue/redshift of A1g, B2g, and A2g modes is found under compressive/tensile strain, respectively. Interestingly, mode A1g exhibits the maximum/minimum shift while mode B2g and mode A2g present the minimum/maximum shift when the strain is applied along armchair/zigzag direction. Density functional theory calculations are carried out to investigate the anisotropic strain response mechanism, finding that the strain‐induced regulation of the P─P bond angle, bond length, and especially interlayer interaction has a giant influence on the Raman shift.