The discrete element model of Mactra veneriformis currently employs an oversimplified multi-sphere approach using EDEM’s Hertz–Mindlin model, assuming uniform shell–flesh mechanical properties. This study developed an advanced dual-layer flexible bonding model through comprehensive biomechanical testing. Mechanical properties and shell morphology were experimentally characterized to inform model development. Parameter optimization combined free-fall experiments with Plackett–Burman screening, steepest ascent method, and Box–Behnken RSM, yielding optimal contact parameters: flesh–flesh stiffness (X1) = 3.64 × 1011 N/m3, shell–flesh interface (X3) = 1.48×1013 N/m3, shell–shell tangential stiffness (X6) = 3.23 × 1012 N/m3, and normal strength (X7) = 8.35 × 106 Pa. Validation showed only 4.89% deviation between simulated and actual drop tests, with hydraulic impact tests confirming excellent model accuracy. The developed model accurately predicts mechanical behavior and shell fracture patterns during harvesting operations. This research provides a validated numerical tool for optimizing clam cultivation and harvesting equipment design, offering significant potential to reduce shell damage while improving harvesting efficiency in bivalve aquaculture systems.
Xu et al. (Mon,) studied this question.