Calibration of the Discrete Element Method Parameters in Living Juvenile Manila Clam ( Ruditapes philippinarum ) and Seeding Veriﬁcation

: The Manila Clam is an important economic shellﬁsh in China’s seafood industry. In order to improve the design of juvenile Manila Clam seeding equipment, a juvenile clam discrete element method (DEM) particle shape was established, which is based on 3D scanning and EDEM software. The DEM contact parameters of clam-stainless steel, and clam-acrylic were calibrated by combining direct measurements and test simulations (slope sliding and dropping). Then, clam DEM simulation and realistic seeding tests were carried out on a seeding wheel at different rotational speeds. The accuracy of the calibrated clam DEM model was evaluated in a clam seeding veriﬁcation test by comparing the average error of the variation coefﬁcient between the realistic and simulated seeding tests. The results showed that: (a) the static friction coefﬁcients of clam-acrylic and clam-stainless steel were 0.31 and 0.23, respectively; (b) the restitution coefﬁcients of clam-clam, clam-acrylic, and clam-stainless steel were 0.32, 0.48, and 0.32, respectively. Furthermore, the results of the static repose angle from response surface tests showed that when the contact wall was acrylic, the coefﬁcient rolling friction and static friction of clam-clam were 0.17 and 1.12, respectively, and the coefﬁcient rolling friction of clam-acrylic was 0.20. When the contact wall was formed of stainless steel, the coefﬁcient rolling friction and static friction of clam-clam were 0.33 and 1.25, respectively, and the coefﬁcient rolling friction of clam-stainless steel was 0.20. The results of the veriﬁcation test showed that the average error between the realistic and simulated value was <5.00%. Following up from these results, the clam DEM model was applied in a clam seeding simulation.


Introduction
The Manila Clam (Ruditapes philippinarum) is one of the most important farmed shellfish in China and has high economic benefit. China's Manila Clam production exceeded 4 × 10 6 t in 2019 [1]. The seeding quality of juvenile Manila Clams (clam) is an important factor affecting production. Currently, artificial seeding is the foremost clam seeding method and has the problems of high labor intensity, uneven seeding, and considerable product damage [2]. The clams' mechanized seeding can reduce farming costs and increase economic returns. However, there is no mechanized clam seeding equipment on the market at the moment. Therefore, the development of mechanized seeding equipment is imperative.
If designing clam seeding equipment by traditional empirical design methods, there are disadvantages such as time consumption and high testing cost, and it is difficult to accurately analyze the complex dynamic behavior of clam during the mechanized seeding process. The discrete element method (DEM) can accurately simulate the mechanized seeding of clams and better reflect the state of particles in complex motion, with the advantage of improving design efficiency and cost savings, and DEM has been widely used in the research and design of agricultural equipment [3][4][5][6][7][8]. However, when DEM is used to simulate the seeding process, the accurate DEM particle model (clam DEM shape and simulation contact parameters) is the key factor that can improve the accuracy of the simulation results and is thereby conducive to the structural design of the clam seeding equipment and the actual seeding effect. Due to the morphological and physical trait differences between a clam and the DEM particle model, there will inevitably be some errors between the realistic contact parameters and the simulated contact parameters. Hence, the simulated contact parameters need calibrating [9][10][11].
Several studies on the establishment of an agricultural granular material DEM model and the acquisition of simulation contact parameters have been undertaken. There are two methods to obtain the simulation contact parameters, namely a direct measuring approach and a DEM simulation calibration approach [12]. The contact parameters of some granular materials have been measured directly, e.g., flax seeds [13], corn, and olives [14]. However, the accuracy of the granular material DEM model is low due to the shape differences between reality and simulation, and the subsequent large measuring error [12,15]. The DEM simulation calibration approach, a laboratory-scale experiment, is performed with a granular material, which is then simulated using DEM. The DEM input simulation contact parameters are then repetitively changed until the results are similar to the actual experiment's outputs [16]. The DEM simulation calibration approach has the advantage of a reduced error and greater accuracy, and it has been widely used in agricultural engineering. Therefore, several studies have calibrated the simulated contact parameters of agricultural granular materials through DEM simulation calibration. The calibration objects include corn [17][18][19][20], dried cherry fruit [21], rice [3,22,23], wheat [11], soybean [24,25], minipotatoes [9,26], rape stalks [27], peanuts [28], and yams [29]. Although DEM simulation calibration has been widely used in the calibration field, there are still some deficiencies requiring improvements, such as few studies on the simulation calibration of the contact parameters of heterogeneous and live materials, e.g., clam. Not all parameters can be obtained by simulation calibration approach, etc. Therefore, the combination of experimental design method and simulation method was used for predicting the unavailable parameters in this study, so as to achieve the purpose of consistency between simulation effect and actual effect.
To provide more accurate reference to the design and optimization of the clam seeding equipment, the living clams were taken as research objects. A clam DEM particle shape was established, as well as calibration of the simulation contact parameters between clam and various materials. The accuracy of the clam DEM model was verified with a clam seeding test.

Materials
In this paper, the calibration object was a 15-20 mm long clam used for seeding in the shellfish farm area of Dandong City, Liaoning Province, China. According to the seeding equipment design requirements, 5 mm thick acrylic and stainless steel plates were selected as the contact wall materials. The initial set of material parameters for acrylic (AC), stainless steel (SS), and clam are shown in Table 1.

Clam DEM Model
Due to the oval and irregular clam shape (Figure 1a), the clam three-dimensional (3D) shape ( Figure 1b) was difficult to draw accurately, and so it was created by reverse engineering with a 3D scanner (Einscan-S, Shining3D, CHN) and Geomagic Studio software (Geomagic Studio 2016, Geomagic, Morrisville, NC, USA). Then the clam DEM shape (Figure 1c), based on the clam 3D shape, was filled with different size particles using the multi-sphere method (MSM) in EDEM software (EDEM 2018, DEM Solutions Ltd., Edinburgh, UK) [22]. The Hertz-Mindlin no-slip contact model, one of the most commonly employed contact models, was selected in EDEM [22,23,32]. The simulation contact parameters input in EDEM between the living clam and different contact materials were calibrated using the DEM simulation calibration approach. Due to the oval and irregular clam shape (Figure 1a), the clam threeshape (Figure 1b) was difficult to draw accurately, and so it was created neering with a 3D scanner (Einscan-S, Shining3D, CHN) and Geomagic (Geomagic Studio 2016, Geomagic, Morrisville, NC, USA). Then the clam ure 1c), based on the clam 3D shape, was filled with different size particle sphere method (MSM) in EDEM software (EDEM 2018, DEM Solutions UK) [22]. The Hertz-Mindlin no-slip contact model, one of the most com contact models, was selected in EDEM [22,23,32]. The simulation contact in EDEM between the living clam and different contact materials were the DEM simulation calibration approach.

The Static Friction Coefficient (μs)
The particle-wall static friction coefficient (μs-pw) was determined u plate method [14,33]. The clam and digital inclinometer (DMI 410, Rio were placed on the inclined plate. One end of the inclined plate was fix end ascended gradually with the lifting platform. When the clam began t value of the inclination angle (β) was recorded by digital inclinometer a was calculated (Figure 2a). Furthermore, the value of the particle-part coefficient (μs-pp) was difficult to measure directly and the μs-pp was deter 2.4.2. The particle-wall static friction coefficient (µ s-pw ) was determined using the inclined plate method [14,33]. The clam and digital inclinometer (DMI 410, Rion, Beijing, CHN) were placed on the inclined plate. One end of the inclined plate was fixed, and the other end ascended gradually with the lifting platform. When the clam began to slide down, the value of the inclination angle (β) was recorded by digital inclinometer and the µ s-pw value was calculated ( Figure 2a). Furthermore, the value of the particle-particle static friction coefficient (µ s-pp ) was difficult to measure directly and the µ s-pp was determined in Section 2.4.2.
value of the inclination angle (β) was recorded by digital i was calculated (Figure 2a). Furthermore, the value of the coefficient (μs-pp) was difficult to measure directly and the 2.4.2.

The Restitution Coefficient (e)
The restitution coefficient (e) was determined through drop tests [14]. SS, AC, and clam (the clams had been closely adhered to the SS plate) were used as the bottom plate, and the graph paper in 5 × 5 mm form was placed vertically as the experimental background to assist in measuring the height of rebounding. The clam fell freely from a height (H 0 ) of 150 mm and initially rebounded to a maximum height (H 1 ) [32,34]. The drop process was recorded with a high-speed video camera (Fastcamsa5, Photron, Tokyo, Japan) running at 1024 × 1024 pixels@1000 fps (Figure 3a). Then the video was converted into images using Matlab (Matlab 2014, Mathworks, Natick, MA, USA) to measure H 1 . The value of e was calculated by Equation (1).

The Restitution Coefficient (e)
The restitution coefficient (e) was determined through drop tests [14 clam (the clams had been closely adhered to the SS plate) were used as the and the graph paper in 5 × 5 mm form was placed vertically as the expe ground to assist in measuring the height of rebounding. The clam fell freely (H0) of 150 mm and initially rebounded to a maximum height (H1) [32,34]. cess was recorded with a high-speed video camera (Fastcamsa5, Photron, running at 1024 × 1024 pixels@1000 fps (Figure 3a). Then the video was c images using Matlab (Matlab 2014, Mathworks, Natick, MA, USA) to me value of e was calculated by Equation (1).
Before impact Impact After impact Before impact Impact After impact  The static repose angle reflects the internal friction and scattering ch granular materials. The clam static repose angle was measured by a liftin [15]. A hollow cylinder (R = 50 mm, H = 350 mm) was placed vertically plate, and 500 clams were poured into it. When the 500 clams in the hollow in a stable state, the cylinder was elevated at a speed of 0.1 m/s [35]. The 5 steadily deposited due to gravity, and the static repose angle of the four X−, Y+, Y− was measured with a digital inclinometer ( Figure 4a).

The Static Repose Angle (θ)
The static repose angle reflects the internal friction and scattering characteristics of granular materials. The clam static repose angle was measured by a lifting cylinder test [15]. A hollow cylinder (R = 50 mm, H = 350 mm) was placed vertically on the bottom plate, and 500 clams were poured into it. When the 500 clams in the hollow cylinder were in a stable state, the cylinder was elevated at a speed of 0.1 m/s [35]. The 500 clams were steadily deposited due to gravity, and the static repose angle of the four directions: X+, X−, Y+, Y− was measured with a digital inclinometer ( Figure 4a). plate, and 500 clams were poured into it. When the 500 clams in in a stable state, the cylinder was elevated at a speed of 0.1 m/s steadily deposited due to gravity, and the static repose angle o X−, Y+, Y− was measured with a digital inclinometer (Figure 4a).

Calibration of the Fixed Contact Parameters
Calibration of the simulation contact parameters were carrie the DEM simulation calibration results consistent with the direct simulation test was run as shown in Figures 2b, 3b and 4b. The (simulation inclination angle (β') and simulation rebound height EDEM analysis. The simulation contact parameters range was pr rect measurements and extending the range at both ends. The si eters range was divided into 5 groups, each group with a step were respectively input into the DEM simulation calibration te measured contact parameters in Sections 2.3.1 and 2.3.2 (e, μs-pw) w

Calibration of the Fixed Contact Parameters
Calibration of the simulation contact parameters were carried out by EDEM. To keep the DEM simulation calibration results consistent with the direct measurements, the DEM simulation test was run as shown in Figures 2b, 3b and 4b. The DEM simulation results (simulation inclination angle (β') and simulation rebound height (H 1 )) were measured by EDEM analysis. The simulation contact parameters range was predicted by taking the direct measurements and extending the range at both ends. The simulation contact parameters range was divided into 5 groups, each group with a step length of 0.05, and they were respectively input into the DEM simulation calibration test [10]. Only one of the measured contact parameters in Section 2.4.1 and Section 2.4.2 (e, µ s-pw ) was taken as the variable, and the other simulation contact parameters (the simulation contact parameters without influence on the calibration test results) were set to 0 in each DEM simulation calibration test. For example, only the particle-wall simulation static friction coefficient (µ s-pw ) had an impact on the β' results in the DEM simulation calibration test, so the remaining contact parameters, apart from µ s-pw , were set to 0. The DEM simulation test data were entered into Original software (Original2019, Origin Lab, Northampton, MA, USA) to obtain the quadratic polynomial fitting curve and equation. The accuracy of the quadratic polynomial fitting equation was evaluated by the R 2 value. Subsequently, the results of direct measurement parameters (β and H 1 ) were brought into the quadratic polynomial fitting equation as the target values to calculate and calibrate the simulation contact parameters. The calibrated simulation contact parameters were re-entered into the EDEM software, and the DEM simulation test was repeated to acquire the results. The relative error (δ) between the DEM simulation test results and direct measurements was calculated to verify the accuracy of the simulation contact parameters.

Response Surface Simulation Test
As the particle-wall (µ r-pw ) and particle-particle (µ r-pp ) rolling friction coefficients and the particle-particle static friction coefficient (µ s-pp ) were difficult to measure directly, the simulation contact parameters (particle-wall simulation rolling friction coefficient (µ r-pw ), particle-particle simulation static friction coefficient (µ s-pp ), particle-particle simulation rolling friction coefficient (µ r-pp )) could not be obtained by the method in Section 2.4.1. Therefore, response surface methodology (RSM) was used to find the best combination of simulation contact parameters. Box-Behnken Design (BBD), an RSM method, was applied. The calibrated simulated contact parameters in Section 2.4.1 were input into EDEM for a clam stacking simulation pre-test, and the value range of the other uncalibrated simulation contact parameters (µ r-pp-ss , µ s-pp-ss , µ r-pw-ss , µ r-pp-ac , µ s-pp-ac , µ r-pw-ac ) was estimated. Then a 3 factor and 3 level (high (+1), medium (0), and low (−1)) BBD test was applied for the numerical design [21], and the simulation repose angles of Clam-SS (θ' ss ) and Clam-AC (θ' ac ) were evaluated as the response index.
The analysis of variance (ANOVA) was applied to assess the main effects and interactions. A multiple regression equation, without non-significant factors, was obtained based on ANOVA results of BBD tests. The directly measured clam static repose angle was placed as the target value into the multiple regression equation, and the optimal combinations of the other simulation contact parameters (µ r-pw , µ s-pp , µ r-pp ) under different contact parameters were calculated. All the calibrated simulation contact parameters were input into the EDEM software for repeating the simulation lifting cylinder test, and the simulation static repose angle was measured. Additionally, by calculating the clam static repose angle relative error (δ a ) between the measured and simulated values, the accuracy of the clam DEM model was verified.

Clam DEM Model Seeding Verification Test
The accuracy of the calibrated clam DEM model was evaluated in a clam seeding verification test by comparing the average error of the variation coefficient between the direct and simulated seeding tests. Taking the seeding wheel rotation speed as the variable, a clam seeding test of the five groups was conceived. The clam feeding speed of the blanking hopper was 0.2 kg s −1 , and the conveyor belt speed (i.e., seeding equipment advancing speed) was 0.75 m s −1 in the test. In each test, the weight of clams discharged onto conveyor belts (i.e., clam seeding area) by seeding wheels on five continuous grids was measured. The clam seeding equipment, designed and used in direct seeding tests, consisted of a blanking hopper, a seeding tray, and a seeding wheel. The seeding tray and wheel were manufactured through 3D printing, and the inner wall of the seeding wheel adhered to AC plates (Figure 5a). The 3D seeding equipment model, clam DEM model, and the AC plate were imported into EDEM software for the seeding simulation test (Figure 5b). The variation coefficient was calculated according to Equations (2)-(4). The testing protocol of the DEM simulation and realistic seeding tests were the same: where C V was the variation coefficient, %; S was the standard deviation, g; n was the grid number; M 0 was the average feed weight collected by the grid in the collection domain, g; M i was the feed weight collected in i computing grid, g.
ing 2021, 3 and the AC plate were imported into EDEM software for the seeding ure 5b). The variation coefficient was calculated according to Equation protocol of the DEM simulation and realistic seeding tests were the s where CV was the variation coefficient, %; S was the standard deviati number; M0 was the average feed weight collected by the grid in the Mi was the feed weight collected in i computing grid, g.

Statistical Analysis
Each direct measurement group and DEM simulation calibratio 10 times to calculate the mean value. The response surface test result Design Expert software (Design expert 10, Stat-Ease, Minneapolis, M nificant impact level of 95% (p < 0.05), and a highly significant imp 0.01).

Statistical Analysis
Each direct measurement group and DEM simulation calibration test were repeated 10 times to calculate the mean value. The response surface test results were calculated by Design Expert software (Design expert 10, Stat-Ease, Minneapolis, MN, USA), with a significant impact level of 95% (p < 0.05), and a highly significant impact level of 99% (p < 0.01).

The Static Friction Coefficient (µ s )
The direct measurement results of β and µ s-pw between clam and SS, and AC are shown in Table 2. Due to a µ s-pw-ss and µ s-pw-ac of 0.26 and 0.34, respectively, the resulting µ s-pw prediction range was 0.20-0.40. The quadratic polynomial fitting curve and equation y s-ss , y s-ac based on the DEM simulation test were fitted and are shown in Figure 6a. The coordinates obtained by virtual line marking in Figure 6a are the simulation contact parameters (µ s-pw ) and their test results (β') in the DEM simulation calibration tests.    Table 2, the directly measured βac and μs-pw-ac were 1.31 times the βss and μs-pw-ss, and a similar result occurred in the DEM simulation test. Due As shown in Table 2, the directly measured β ac and µ s-pw-ac were 1.31 times more than the β ss and µ s-pw-ss , and a similar result occurred in the DEM simulation test. Due to the AC surface roughness, the static friction between clam and AC was sizeable. Therefore, in the clam seeding equipment design, AC would be the most appropriate material when increasing friction was required. For example, to avoid the clams dropping too quickly and causing congestion at the outlet of the dropping hopper, it would be necessary to add a high friction guide plate in the blanking hopper; hence, AC would be the recommended material. Conversely, SS would be the obvious choice when friction needs to be minimized.
The simulation value µ s-pw was significantly smaller than µ s-pw (as shown in Table 2). The reason may be that the clam DEM model was composed of smooth particles; however, clams have a growth line on their surface that grows with increasing clam size. As the clam growth lines gradually deepened, the surface roughness and the static friction coefficient (µ s ) also gradually increased, which concurs with previous studies that suggest that the larger the clam size, the greater the static friction coefficient is [31]-hence why the µ s-pw was larger than µ s-pw , as also in the calibration case of panax notoginseng seeds [36].
Furthermore, the value of the inclination angle relative error (δ a ) between β' and β included: δ a-ss = 4.0%, δ a-ac = 1.8%, respectively. It indicated that the simulation inclination angle was consistent with the direct measurements. The DEM calibration test results of the static friction coefficient were accurate, which could then be applied to the clam simulation study.

The Coefficient of Restitution (e)
The coefficients of restitution (e) for SS, AC, and Clam are shown in Table 3. According to the calculated restitution coefficient results: e pw-ss = 0.30, e pp = 0.21, e pw-ac = 0.41, the predicted range of the simulation restitution coefficient was 0.25-0.55. Figure 6b shows the relationship between the simulation restitution coefficient (e') and rebound height (H' 1 ) in the fitting curve and equation y h-ss , y h-ac, y h-pp . The coordinates obtained by virtual line marking in Figure 6b are the simulation contact parameters (e') and their test results (H 1 ') in the DEM simulation test. As illustrated in Table 3, the e varied in direct measurement, and the e pp < e pw-ss < e pw-ac . The lower the restitution coefficient, the lower rebound height is. In the process of clam seeding, when clams fall into the groove on the seeding wheel through the bottom of the blanking hopper, they rebound. When the rebound height of the clams is higher than the depth of the groove, the clams are crushed by the seeding wheel and the seeding traying. To avoid this affecting, SS, with a small restitution coefficient, should be selected as the surface contact material of the seeding wheel. AC with a high restitution coefficient would be an appropriate contact material if the impact of rebound height was negligible and the clam breakage rate could be reduced by modifying the equipment structure.
Additionally, the e' was greater than in Table 3; specifically, the e' pw-ac was 0.48, which is 17.1% larger than the e' pw-ac . This may be because the center of gravity of the clam is different from that of its DEM model. The clam is composed of an external shell, internal flesh, and a small amount of water, which are heterogeneous granular materials. Due to the different shapes and water content between clams, the center of gravity of each clam also varies. Therefore, when each clam lands on the bottom plate, the impact position and rebound height are different. However, the clam DEM model in the simulated drop test was filled with solid homogeneous granular materials, and the gravity center and impact position were more fixed than the living clam. Therefore, the direct measurement rebound height was significantly lower than in the DEM simulation test result.
The rebound height relative error (δ H1 ) between H' 1 and H 1 was: 1.7%, 1.7%, 2.1%, respectively. The DEM simulation test result was similar to the direct measurement, which could effectively replace the realistic drop test.

Response Surface Simulation Test and ANOVA
The results of the directly measured static repose angles of Clam-SS (θ ss ) and Clam-AC (θ ac ) were θ ss = 31.75 • , θ ac = 38.07 • . The range of the simulation contact parameters was predicted by a clam stacking simulation pre-test. With an SS wall, the simulation rolling coefficient of Clam-Clam (µ r-pp ) was in the range of 0.14−0.22, the simulation statics coefficient of Clam-Clam (µ s-pp ) was in the range of 1.04−1.12, and the rolling coefficient of Clam-SS (µ r-pw-ss ) was in the range of 0.14−0.22. The simulation contact parameter range for an AC wall was also predicted. The factors and levels from the response surface simulation test are shown in Table 4. Table 4. Factors and levels.

SS Factors AC Factors
µ r-pp-ss µ s-pp-ss µ r-pw-ss µ r-pp-ac µ s-pp-ac µ r-pw-ac In this study, 17 experiments were carried out to find the best combination of simulation contact parameters and to study the effect of the µ r-pp , µ s-pp , and µ r-pw on the clam simulation static repose angle, based on the BBD method [37]. The corresponding simulation results are shown in Table 5. Table 5. Experimental scheme and response results based on the BBD method.

Runs
Factors Response In ANOVA, the p-value represents the significance of the factors. The F-value represents the primary and secondary order of influence that the factors had on the response. The larger the F-value, the stronger the influence on the response was. The ANOVA results for the quadratic polynomial model are shown in Table 6. With an SS bottom plate, the simulation contact parameters: µ r-pp , µ s-pp , µ r-pw , µ r-pp µ s-pp , µ r-pp µ r-pw , µ s-pp 2 , and µ r-pp 2 showed highly significant influence (p < 0.01), whereas µ s-pp µ r-pw and µ r-pw 2 showed insignificant influence. The influence order of the factors was µ s-pp > µ s-pp 2 > µ r-pp > µ r-pp µ s-pp > µ r-pp 2 > µ r-pp µ r-pw > µ r-pw > µ r-pw 2 > µ s-pp µ r-pw . With an AC bottom plate, the simulation contact parameters: µ r-pp , µ r-pw 2 showed highly significant influence and µ s-pp , µ s-pp µ r-pw , µ s-pp 2 showed significant influence (p < 0.05), whereas µ r-pp µ r-pw and µ r-pw 2 showed insignificant influence. The influence order of the factors was µ r-pp > µ r-pw 2 > µ s-pp > µ s-pp 2 > µ s-pp µ r-pw > µ r-pp µ s-pp > µ r-pp 2 > µ r-pw > µ r-pp µ r-pw . Note: * shows that the item is significant (p < 0.05); ** shows that the item is extremely significant (p < 0.01).
In both SS and AC regression models, the parameters such as the lack of fit p value, the coefficient of variation (CV), determination coefficient (R s 2 ), correction determination coefficient (Adj−R s 2 ), and the Adeq−Precision demonstrated good predictability with the multiple regression equation (Equations (5) and (6)).
The static repose angle in the stacking test was determined as θ ss < θ ac by comparing the direct measurement AC and SS results. This may be because the roughness of the AC surface is greater than that of smoother SS. The larger the µ s-pw-ac , the larger the θ ac with the AC bottom plate is. It could be inferred that the µ s-pw , µ s-pw , µ r-pp , and µ s-pp have a significant effect on the static repose angle of clams, as illustrated in Table 6. Consequently, it could be inferred that the clam contact parameters µ s-pw , µ r-pp , µ s-pp had a significant effect on the application of the clam DEM model, and the response surface simulation test conclusions were similar to the published result [11]. Therefore, the influence of µ s-pw , µ r-pp , and µ s-pp on the seeding effect should be fully considered when designing and optimizing clam seeding equipment.

Clam Seeding Verification Test
As shown in Figure 7, as the rotation speed of the seeding wheel increases to 30 r min −1 , the δ t value and the realistic and simulated values of the coefficient of variation decrease; however, when this rotation speed increases from 30 r min −1 to 45 r min −1 , the said values start to increase. When the seeding wheel speed was around 30 r min −1 , all the variation coefficients were at their lowest values. The smaller the variation coefficient, the more uniform the distribution of clams is. Therefore, an optimal working rotation speed of around 30 r min −1 , which is the best matching relationship with the feeding speed and the equipment advancing speed in this study, for clam seeding equipment is recommended. With a slower seeding wheel rotation speed, the seeding continuity was reduced, consequently increasing the clam seeding variation coefficient. On the contrary, when the rotation speed of the seeding wheel was higher, the groove of the seeding wheel directly below the clam blanking hopper was not full; the groove had already rotated away. Therefore, the clam seeding variation coefficient and test error were large. When the seeding wheel rotation speed was 30 r min −1 , the feeding speed and advancing speed were more appropriate. Moreover, the verification results were also helpful in the design and optimization of the clam seeding equipment and could effectively help to understand the motion state of clams for future designers, whilst providing a technical reference.
variation coefficients were at their lowest values. The smaller the varia more uniform the distribution of clams is. Therefore, an optimal work of around 30 r min −1 , which is the best matching relationship with the the equipment advancing speed in this study, for clam seeding equ mended. With a slower seeding wheel rotation speed, the seeding conti consequently increasing the clam seeding variation coefficient. On the rotation speed of the seeding wheel was higher, the groove of the seed below the clam blanking hopper was not full; the groove had already ro fore, the clam seeding variation coefficient and test error were large. wheel rotation speed was 30 r min −1 , the feeding speed and advancing appropriate. Moreover, the verification results were also helpful in th mization of the clam seeding equipment and could effectively help to u tion state of clams for future designers, whilst providing a technical re As shown in Figure 7, the minimum and maximum values of the v error (δt) were 9.65%, and 1.60%, and the average value of δt was showed that the clam seeding realistic and simulated tests were highly be inferred that the accuracy of the clam DEM model was high. Theref ing simulation could effectively simulate actual clam seeding.

Conclusions
In this study, based on the clam DEM particle shape, the clam parameters and a clam seeding verification test were calibrated and The conclusions are as follows: (1) The clam simulation contact parameters of different contact mat were determined by a DEM simulation calibration test. When SS w terial, the simulation contact parameters in the calibrated set inclu e'pw-ss = 0.38, e'pp = 0.32, μ′r-pp-ss = 0.33, μ′s-pp-ss = 1.25, μ′r-pw-ss = 0.34. S μ′s-pw-ac =0.31, e'pw-ac = 0.48, e'pp = 0.32, μ′r-pp-ac = 0.20, μ′s-pp-ac = 1.12, μ′r-(2) The accuracy of the simulation contact parameters (μ′s-pp, μ′r-pp, μ′s As shown in Figure 7, the minimum and maximum values of the variation coefficient error (δ t ) were 9.65%, and 1.60%, and the average value of δ t was 4.98%. The results showed that the clam seeding realistic and simulated tests were highly similar, so it could be inferred that the accuracy of the clam DEM model was high. Therefore, the clam seeding simulation could effectively simulate actual clam seeding.

Conclusions
In this study, based on the clam DEM particle shape, the clam simulation contact parameters and a clam seeding verification test were calibrated and analyzed by DEM. The conclusions are as follows: (1) The clam simulation contact parameters of different contact materials (SS and AC) were determined by a DEM simulation calibration test. When SS was the contact material, the simulation contact parameters in the calibrated set included: µ s-pw-ss = 0.22, e' pw-ss = 0.38, e' pp = 0.32, µ r-pp-ss = 0.33, µ s-pp-ss = 1.25, µ r-pw-ss = 0.34. Similarly, with AC, µ s-pw-ac =0.31, e' pw-ac = 0.48, e' pp = 0.32, µ r-pp-ac = 0.20, µ s-pp-ac = 1.12, µ r-pw-ac = 0.17. (2) The accuracy of the simulation contact parameters (µ s-pp , µ r-pp , µ s-pw ) was simultaneously verified, and they were found to have a significant effect on the clam DEM model. (3) Through comparing the results of the clam realistic and simulated seeding tests, it was verified that the accuracy of the clam DEM model was high, and therefore it could be applied for clam seeding DEM simulation and equipment design in the future.