Using the discrete element method (DEM) to simulate the rice machine transplanting operation is important for assessing the plant injury and optimizing the rice transplanter performance, while the DEM flexible model establishment that can accurately reflect the mechanical properties of the rice blanket seedling root blanket is an important foundation. Based on the root blanket's stratification and the root system structure's measurement and statistics, a new method for root blanket flexible modeling was proposed in this study. Firstly, the Hertz-Mindlin with bonding V2 contact model was used to establish substrate I (SI), substrate II (SII), substrate III (SIII), stemroot combination (SRC), and netted layer (NL) flexible models, respectively, and the model parameters were calibrated and determined by angle of repose (AOR), direct shear, and mechanical tests. The calibration results showed that the deviations of AOR simulated values for SI and SII were both less than 1.5 %, and the deviations of shear strength simulated values were both less than 4 %. Secondly, the shear characteristics of SI and SII were determined by direct shear test. The results showed that the physical and simulated shear stress-displacement relationship curves of SI and SII were basically the same; the hair roots mainly relied on the cohesive between them and the substrate to improve the substrate strength; the fitted lines of simulated shear strength and normal stress of SI and SII were in high agreement with these of the measured values; the deviations of the simulated cohesion and internal friction angle were both less than 5 %. After that, the Hertz Mindlin with JKR V2 contact model was used between SRC and substrate. The interfacial surface energy of the root blanket and the bonding parameters of SIII were calibrated by stem, half-SRC, and SRC pulling-out tests layer by layer. The calibration results showed that the deviation of the maximum pulling-out force of SRC was 5.83 %, verifying that the model could accurately simulate the intertwining effect of the crown roots. Finally, the flexible model of the root blanket was verified by cutting, curling, and tensile tests. The simulated test results were consistent with the trends of the physical test results; the deviations of the maximum cutting resistance of front cutting and side cutting were both within 8 %, the error percentage range of the marked points height was 0.35 % to 17.16 %, and the deviation of the maximum tensile force was 9.22 %, indicating the good feasibility of the modeling method and accuracy of the flexible model. The results of this study lay a foundation for the DEM simulation of the rice machine transplanting operation. They can also provide a reference for the numerical simulation of other multiplant root-soil complexes.