Discrete Element Modelling Method and Parameter Calibration of Mussel Based on Bonding V2 Model

To address the inefficiency and high labor intensity associated with traditional manual mussel seedling unloading, this study proposes an automated traction-rope mussel unloading machine. This study focuses on the thick-shelled mussel (Mytilus coruscus) as the research subject. Furthermore, the key mussel unloading processes were simulated using the EDEM software to analyze mechanical interactions during detachment. A breakable mussel discrete element model was developed, and its Bonding V2 model parameters were systematically calibrated. Using the ultimate crushing displacement (2.25 mm) and ultimate crushing load (552 N) as response variables, the model was optimized through a sequential experimental design comprising Plackett–Burman screening, the steepest ascent method, and the Box–Behnken response surface methodology. The results demonstrate that the optimal parameter combination consists of unit area normal stiffness (2.48 × 1011 N/m3), unit area tangential stiffness (3.80 × 108 N/m3), critical normal stress (3.15 × 106 Pa), critical tangential stress (2.90 × 107 Pa), and the contact radius (1.60 mm). The model’s accuracy was validated through integrated discrete element simulations and prototype testing. The equipment achieves an exceptionally low mussel damage rate of only 1.2%, effectively meeting the operational requirements for mussel unloading. This study provides both theoretical foundations and practical insights for the design of mechanized mussel unloading systems in China.