The forage crop Caragana korshinskii Kom. is of high quality, and the biomechanical properties of its plant system are of great significance for the development of harvesting equipment and the comprehensive utilisation of crop resources. However, the extant research on the biomechanical properties of Caragana korshinskii Kom. is inadequate to enhance and refine the theoretical techniques for mechanised harvesting. In this study, we established a discrete element model of CKS based on the Hertz-Mindlin bonding contact model. By combining physical experiments and numerical simulations, we calibrated and validated the intrinsic and contact parameters. The Plackett-Burman design test was employed to identify the significant factors influencing bending force, and the optimal parameter combination for these factors was determined through response surface analysis. When the shear stiffness per unit area was 3.56×109 Pa, the bonded disk scale was 0.93 mm, the normal stiffness per unit area was 9.68×109 Pa, the normal strength was 5.62×107 Pa, the shear strength was 4.27×107 Pa, the discrete element numerical simulation results for three-point bending, radial compression, axial tension, and shear fracture exhibited a maximum failure force error of 3.32%, 4.37%, 4.87% and 3.74% in comparison to the physical experiments. In the cutting experiments, a smaller radial angle between the tool edge and the stem resulted in less damage to the cutting section, which was beneficial for the smoothness of the stubble after harvesting and the subsequent growth of the stem. The discrepancy in cutting force between the physical and numerical simulations was 3.89%, and the F-x (force versus displacement) trend was consistent. The multi-angle experimental validation demonstrated that the discrete element model of CKS is an accurate representation of the real biomechanical properties of CKS. The findings offer valuable insights into the mechanisms underlying crop-machine interactions.