Next Article in Journal
Weight Prediction of Landlly Pigs from Morphometric Traits in Different Age Classes Using ANN and Non-Linear Regression Models
Previous Article in Journal
Development of a General Prediction Model of Moisture Content in Maize Seeds Based on LW-NIR Hyperspectral Imaging
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An ANSYS/LS-DYNA Simulation and Experimental Study of Sectional Hob Type Laver Harvesting Device

1
College of Mechanical and Power Engineering, Dalian Ocean University, Dalian 116023, China
2
Marine Fishery Equipment Professional Technology Innovation Center of Liaoning Province, Dalian 116023, China
3
Key Laboratory of Environment Controlled Aquaculture, Ministry of Education, Dalian 116023, China
*
Author to whom correspondence should be addressed.
Agriculture 2023, 13(2), 361; https://doi.org/10.3390/agriculture13020361
Submission received: 20 December 2022 / Revised: 27 January 2023 / Accepted: 30 January 2023 / Published: 2 February 2023
(This article belongs to the Section Agricultural Technology)

Abstract

:
To solve the problems of low net harvesting rate, high loss rate, and uneven stubble height during the harvest of laver, the laver (Porphyra yezoensis) was selected as the research object, the analysis of the cultivation mode, biomechanical characteristics, harvesting trajectory and force of laver were carried out. A sectional hob type harvesting device was designed. A rigid-flexible coupling model related to the interaction between the cutting mechanism and the laver was constructed based on ANSYS/LS-DYNA. The Box–Behnken design method was used to simulate the effects of different structural parameters and process parameters on the force of laver cutting, and the bench test of the laver harvesting device was carried out. The simulation results showed that the four factors that significantly affect the force exerted on the laver during cutting in proper order were cutter revolving speed, knife extension length, knife inclination angle and forward velocity. When the combination of the forward velocity, the cutter revolving speed, the knife extension length and inclination angle was 0.77 m/s, 900 r/min, 40 mm, and 110°, respectively, the cutting force on laver was the smallest, which was 4.21 N. The bench test of harvesting performance showed that the cutter revolving speed has a significant impact on the recovery rate, and the forward velocity has a significant impact on the loss rate. When the harvesting speed ratio was λ4 (the cutter revolving speed was 900 r/min and the forward velocity was 0.77 m/s), the net harvesting rate and the loss rate were 97.45% and 3.38%, respectively, and the cutting proportion of laver can reach 77.5%. The results of the study provide a theoretical basis for the development of harvesting for laver.

1. Introduction

Laver (Porphyra yezoensis) is widely distributed in the world and is the main large economic algae cultivated in China [1,2]. China has a total yield of 2.22 × 105 t in 2020, with the economic output value reaching CNY 18.3 billion [3]. Harvesting is a crucial link in the aquaculture of laver. Although some mechanical devices have been applied to laver harvesting in China, such as ZS160-1.8 laver harvester developed by Lianyungang Agricultural Machinery Factory [4,5,6], most of these mechanical devices are designed and manufactured by imitating foreign products. As a result, laver harvesting is limited by a heavy labor burden, low production efficiency, low net recovery rate, and high loss rate [7,8,9]. In contrast, some researchers from Japan, South Korea and other countries have designed some advanced laver harvesting devices, and they also conducted explorations into the interaction between harvesting cutters and laver, as well as the relationship between the cutter life and production efficiency. For example, Yada et al. [10] revealed that the closer the cutting point is to the root of laver, the greater the cutting stress on laver and the better the harvesting effect. Kim et al. [11] found that the longest service life of the cutter can be achieved under the traveling speed of 0.53 m/s and the rotation speed of 700~800 r/min, and the harvesting efficiency can reach 1.15~1.26 kg/min, thus realizing high-quality harvesting of laver [12,13,14,15,16]. However, laver is fixed on a single point and grows on the net curtain, due to the lack of systematic research on the force analysis and fracture state of laver during harvesting, it is difficult to conduct in-depth research on single point support cutting of laver, and the problem of high harvesting damage rate has not been effectively solved.
For the single point support cutting problem of flexible crops such as laver, by establishing a rigid flexible coupling model for simulation research, this model can be employed to effectively reduce the research difficulty, shorten the research cycle, and decrease the development costs [17]. Bílek and Kovařík [18] constructed a rigid-flexible coupling model for determining the parameters of wood sawing. They performed a stress analysis and found that the maximum stress caused by the temperature gradient was about ten times that inferred from the centrifugal force. This stress state was basically consistent with that obtained in a finite element model. Meng et al. [19] established a coupling model related to the interaction between the circular saw blade and branches of Morus alba. They simulated and analyzed the changes of sawtooth number, rotation speed, feeding speed, and branch diameter in the stress and cutting force in the branch-cutting of process. The bench test results verified that the simulated value of the cutting force was consistent with the measured value. Zhao [20] constructed a rigid-flexible “picking mechanism-fruit tree” coupling model based on ANSYS/LS-DYNA to simulate the dynamics of the vibration system. The vibration response characteristics of Camellia oleifera trees were also explored under different excitation frequencies and holding heights. The results proved that the excitation frequency of 38 Hz and the holding height of 1300 mm constituted the optimal harvesting parameter combination. Dun et al. [21] simulated the sawtooth rotary cutting process for the straw of Glycine max. They selected the cutting linear speed, cutter inclination angle, and operation speed as test factors, and analyzed the primary and secondary factors affecting cutting resistance and the optimal cutting parameters. Currently, the application of rigid-flexible coupling model research mainly focuses on the interaction between knives and materials with fixed form such as plant stems or wood, and the research and application of the interaction between knives and crops with greater flexibility such as laver has not been reported. The harvesting of laver involves a complex and random high-speed penetrating collision process. As the cutter and laver contact in a discontinuous dynamic mode, it would be impossible to observe and analyze the interaction between the cutter and laver by conventional methods. According to the crop-cutting theory [22], cutting will cause less damage to laver, while the rupture under a tensile force will cause more damage. The fester caused by damage is the main factor affecting the yield of laver.
Therefore, laver was selected as the research object in this study. The parameter combination for avoiding the missing cutting or re-cutting of laver was explored by analyzing the relationship between the inclination angle of knife and the stress on laver. A rigid-flexible coupling model related to the interaction between the cutting mechanism and the laver was constructed based on ANSYS/LS-DYNA, to perform a harvesting simulation. Moreover, the influence of different structures and operating parameters on the cutting force on laver was analyzed. After a laver harvester prototype was designed and manufactured, the simulation results were verified by a bench test. These efforts may provide a theoretical foundation for the research and development of laver harvesters.

2. Structure and Working Principle of the Harvester

2.1. Physical Property and Harvesting Requirements of Laver

Laver is mainly cultivated on raft net, the cultivation facilities include net curtains and floating rafts. Specifically, net curtains provide the substrate for the growth of laver seedlings, while floating rafts serve as the brackets for hanging these net curtains. The specification of a net curtain for the cultivation of laver is 1.5 m × 12 m [1,23]. Floating rafts can be divided into three patterns according to different cultivation conditions, namely full floating raft, pillar raft, and semi-floating raft (Figure 1) [24,25].
After laver grows to a certain length and the algal quality becomes poor, harvesting contributes to decreased losses timely and increased laver yields. During the harvesting of laver, it is required to keep the re-production ability of the laver. The laver can be harvested when it grows to 15–20 cm, the stubble length is generally 5–7 cm [26,27,28], and the low damage rate of laver shall be ensured during harvesting. The main biological characteristics of laver [29] are shown in Table 1.

2.2. Structure and Principle of the Harvester

In this study, a laver harvester was designed based on the working environment and harvesting requirements. It was composed of a sectional hob mechanism, a stubble height adjustment mechanism, a net curtain support mechanism, a rack, and a hydraulic motor, as shown in Figure 2. Specifically, the sectional hob mechanism (length is 2000 mm) consisted of a knife shaft, 9 pairs of knives, and a connection support plate. The stubble height adjustment mechanism included an arc cover, a protective grid, and a rotary drum. The net curtain support mechanism contained a roller and a connector.
As shown in Figure 3, the laver is fixed on the net curtain at a single point. During the working of the laver harvester, the stubble length can be adjusted in advance by rotating the stubble height adjustment mechanism manually. The net curtain can be supported and flattened by the cooperation between the net curtain support mechanism and the arc cover. The laver can be cut by high-speed rotating knives during the forward movement of the loading device. Laver would fall into the collection box. The device controls the speed of the knife through the manual regulating flow valve, adjusts the speed to the required speed through the indication on the speed sensor, and realizes the emergency stop of the knife through the manual reversing valve.

3. Theoretical Analysis of the Harvesting Process

3.1. Knife Movement Trajectory Analysis

According to the working principle of the laver harvester, the movement trajectory of the knife top presents a trochoid curve combined by circular movement and horizontal movement [30]. The kinematic model of the knife top can be expressed as follows:
X = V0t + Rsinωt
Y = −Rcosωt
where V0 represents the traveling speed of the loading device (m/s), ω represents the angular speed of the cutter shaft (rad/s), t represents time (s), R represents the turning radius of the knife (mm).
The position of the knife top at any time can be expressed by Formulae (1) and (2). There are different shapes and characteristics of the movement trajectory with the changes in the traveling speed (V0), turning radius of the knife (R), and angular speed (ω).
As shown in Figure 3, the speed of the knife top in the horizontal and vertical directions is the derivative of displacement (X) and (Y) to time (t):
V X = d X d t = V 0 + R ω cos ω t
V Y = d Y d t = R ω sin ω t
As shown in Figure 4, due to the fact that the cutter can only contact laver in the solid line part during harvesting, only the movement trajectory of this part would be analyzed. The axle center of the cutter is point P5, and the intersection point O of the net curtain and the Y axis is the origin. When the anterior knife top I contacts the laver at point P1, the laver on the net curtain would be cut. Then, the knife top continues to move to the right side until the posterior knife top leaves the lowest pendulous point of the laver at point P2. The range from point P1 to point P2 indicates the area for cutting laver. Point P1 is the highest cutting point, and the rotation speed overlaps with the traveling speed, which induces the maximum cutting force on the laver. The distance from the net curtain to the cutting baseline represents the target stubble length (H2), and the laver length is H1. The movement trajectory of the posterior knife top II is the same as that of the anterior knife top I. Point P4 represents the position entering the cutting area, the missing cutting length of laver is H3, and the shaded area (P1 P3 P4) represents the missing cutting area.
To avoid missing cutting, reducing the missing cutting length (H3) can effectively reduce the missing cutting area (P1 P3 P4), H3 is related to the movement trajectory of the knife top and the cutting pitch. The cutting pitch (the distance from point P1 to point P3 in Figure 3) is the traveling distance of the harvester during the time interval for the cutting of laver by two adjacent cutters. If the time interval for the cutting of laver by two adjacent cutters is set as:
t s = 2 π 3 ω
the cutting pitch can be expressed as:
S = V 0 t s = 2 π V 0 3 ω
if the ratio of the rotation speed of the cutter to the traveling speed λ is set as:
λ = R ω V 0
the approximate calculation can be completed as per the following formula:
S 2 = a ω V 0 + R sin a
when the value of angle a is not large, it can be considered that sin a=a.
S 2 = a ( V 0 ω + R )
a = S 2 ( V 0 ω + R ) = 2 π R 3 λ 2 ( R λ + R ) = π 3 ( λ + 1 )
hence, the missing cut length can be simplified as:
H 3 = R ( 1 cos a ) = R ( 1 cos π 3 ( λ + 1 ) )
where V0 represents the traveling speed of the loading device (m/s), R represents the turning radius of the knife (mm).
Therefore, the traveling speed, the rotation speed of the cutter, and the turning radius of the knife directly affect the missing cutting length. When the turning radius is constant, the reducing speed ratio λ, the missing cutting length (H3) will decrease, and the missing cutting area.

3.2. Force Analysis in the Cutting Process

When the knife top cuts into the laver, the force is shown in Figure 5.
Laver is selected as the research object, the X axis represents the traveling direction, the Y axis represents the pendulous direction of laver, FX represents the horizontal component of the cutting force (F), and FY represents the vertical component of the cutting force (F), as expressed in Formulae (12) and (13).
FX = Fcosθ + TTrsinθr
FY = Fsinθ + GTrcosθr
where F represents the cutting force of the cutter on laver (N), T represents the horizontal thrust of the knife on laver (the force on laver when the device moves forward) (N), Tr represents the self-elasticity of laver (N), G represents the weight of laver below the action point (N), θ represents the cutting rotation angle (°), θr represents the deflection angle of laver (°).
The horizontal component (FX) is perpendicular to laver. One end of the laver is fixed, and the other end is unconstrained. If the weight of the laver is very small relative to the horizontal component, the laver may jump forward unsteadily, which may cause cutting failure. The weight (G) of the harvested laver will change with the movement of the action point of the knife. The closer to the root, the greater the weight (G) will be. Therefore, in order to cut off laver successfully, it is necessary to consider the weight of the laver below the action point.
The direction of the vertical component (FY) of laver is the same as that of the pendulous tension of laver. Some parts of the laver may be stretched and broken due to its soft texture, and hence there is a tensile force. The cutting stress and tensile stress on laver can be expressed as Formulae (14) and (15):
τ = F cos θ + T T r sin θ r A s
σ = F sin θ + T r cos θ r A s
where As represents the cross-sectional area of laver.
When τ > σ , laver would be cut off; when τ < σ , laver would be ruptured under a tensile force. If θr, F, and As in Formulae (14) and (15) are constant, τ and σ would be proportional to the cutting rotation angle (θ), and laver may be cut off or ruptured under a tensile force with the changes in the cutting rotation angle. When the cutting rotation angle is zero, namely when the cutting edge acts vertically on laver, τ would reach the maximum, and laver can be harvested only by the cutting stress. If the cutting rotation angle increases, τ would decrease; otherwise, σ would increase. However, when the cutting rotation angle is close to 90°, the cutter would not affect the harvesting of laver, it is impossible to complete the harvesting of laver. For rotary cutting, Li and Yada gave similar laws and conclusions [10,31]. In summary, when the cutting rotation angle decreases, the proportion of laver under a cutting force would increase; when the cutting rotation angle increases, the proportion of laver under a tensile force would increase.
The cutting rotation angle (θ) is determined by the knife position, knife structure parameters including extension length (B) and inclination angle (A), and operating parameters including rotation speed (D) and traveling speed (C). Under the same knife position and operating parameters, the cutting rotation angle (θ) is determined by cutter structure parameters. When the extension length (B) increases and the inclination angle (A) decreases, θ increases; otherwise, θ decreases. The single grid of the laver net curtain is 80 mm, and the length of the laver stubble is 50–70 mm. In comprehensive consideration of the structure size of the knife should not be too large and the stubble space requirement, the knife extension length is 35~45 mm. In order to make the top of knife better contact with the laver, the inclination angle is 100~120°. Based on the knife movement trajectory analysis results and the speed of the existing fishing boats, the traveling speed is determined to be 0.51~1.03 m/s, and the rotation speed is 700~1100 r/min.

4. Laver Cutting Simulation Study

A rigid-flexible coupling model related to the interaction between the cutting mechanism and the laver was constructed based on ANSYS/LS-DYNA. The laver model is endowed with flexible body features, and the knife is endowed with rigid body features. The laver model is calibrated by experiment and simulation, the error between the actual value and the simulated value of the fracture limit displacement is 3.75%, and the error between the actual value and the simulated value of the breaking force is 13%, which can be used for simulation research.

4.1. Model Establishment and Simulation Results

There are different lengths and multi-layer stacking in the laver. To realize simplification, the bending caused by stacking may be ignored, the laver model (length is 120 mm, width is 12 mm, thickness is 0.13 mm) was simplified to an inverted conical strip, and the model parameters [29] are shown in Table 2. According to the actual cultivation pattern of laver, the root of laver was fixed. The interference of external factors in the laver movement was also ignored. Meanwhile, to shorten the simulation time, the turning radius of the cutter (111 mm) remained unchanged. Further, the axial length of the cutter was shortened and simplified to a one section hob (length is 667 mm), and the vibration factors in actual work were also ignored. During model establishment, the knife shall be as close to the laver as possible. Moreover, an explicit dynamic analysis was performed on the finite element model of the cutter and laver in the harvesting process by the coupling solution method.
The model related to cutter and laver was established by Solidworks, it was imported into ANSYS/LS-DYNA. The rigid setting was selected for the cutter, while the flexible setting was selected for the laver with the same characteristics. After the laver slices were finely divided, the mesh size was 2 mm, and a total of 15,863 meshes were obtained (Figure 6).
The total simulation time was set as 0.012 s, the root of laver was fixed and restrained, and the linear traveling speed and rotation speed of the cutter were set as 0.77 m/s and 900 r/min, respectively. The simulated harvesting process is shown in Figure 7.
After the Body contact tracker command was added for the laver model slices, the stress curves of the laver model slices in the X and Y axes during harvesting were obtained, as shown in Figure 8a,b, respectively. The time required for cutting laver was 0.005 s. Under the action of the knife on laver, the force on laver in the X and Y axes reached a peak value and decreased rapidly with the rupture of laver. The peak value in the X and Y axes was 4.56 N and 1.52 N, respectively (no force was generated in the Z axis). The peak value of the resultant force on the laver reached 4.8 N after the force combination.

4.2. Response Surface Test Design

In this study, the rotation speed, traveling speed, inclination angle, and extension length of the knife were selected as influencing factors. The cutting force on the laver was selected as an evaluation index to perform a response surface test by the Box–Behnken method [31]. The codes of each factor level in the test are listed in Table 3.

4.3. Test Results and Analysis

The test scheme and results are listed in Table 4. A total of 29 groups of tests were performed, each group of tests was repeated three repeated times. The regression analysis and factor variance analysis were performed with Design-Expert 8.0.6. The significance of each factor on indexes was analyzed, and the regression model of the cutting force on laver was established.
Based on multiple regression fitting of the test data, a regression model established for the influence of various test factors on the cutting force (Y) on laver:
Y = 4.95 + 0.43B + 1.37D + 0.81AB + 0.77BC − 0.58C2
the variance analysis results of response surface tests (Table 5) indicated that there was a statistical significance in the model (p = 0.0006). The lack of fit of the regression model was not significant (p > 0.05), which indicated that the regression model had a high fitting degree, and the p-value of D was less than 0.01, which indicated that this term had a significant influence on the cutting force on the laver. The p-value of B, AB, BC, and C2 was less than 0.05, which indicated that the above four items had a significant influence on the cutting force on laver. Further, there was a quadratic relationship for the influence of related test factors on the response value. The four factors with the most significant influence on the cutting force on the laver can be ranked as the rotation speed (D), extension length (B), inclination angle (A), and traveling speed (C).
As shown in Figure 9a, the interaction between the extension length and the inclination angle of knife can exert a significant influence on the cutting force of laver. The extension length of knife has a larger influence than the inclination angle. The cutting force would increase with the increase in the extension length, and it would decrease first and then increase more significantly with the increase in the inclination angle. The main reason is that the extension length increases the knife torque, which leads to the increase in cutting force on laver. At the same time, the knife inclination angle increases, the cutting force would decrease first and then increase, and the increase is not obvious. When the extension length ranges from 35 mm to 40 mm and the inclination angle ranges from 105° to 115°, the cutting force would be relatively small. As shown in Figure 9b, the extension length exerts a more significant influence on the cutting force compared with the traveling speed. The cutting force increases with the extension length, and it would increase first and then decrease with the increase in the traveling speed. The main reason is that the traveling speed increases, the speed variation increases, and the impact force on the laver gradually increases. At the same time, with the increase in traveling speed, the cross-sectional area changes with the position of the cutting point of the knife. With the increase of traveling speed, the cutting point of the knife moves close to the root of the laver and then moves away from the root of the laver, and the cross-sectional area of the laver first increases and then decreases, resulting in an increase in the cutting force on laver first and then decreases. When the extension length ranges from 35 mm to 40 mm and the traveling speed ranges from 0.77 m/s to 1.03 m/s, the cutting force would be relatively small.

4.4. Parameter Optimization

Taking the minimum cutting force was selected as the optimization objective, the objective function can be expressed as follows:
{ m i n   Y ( A , B , C , D ) 100 ° A 120 ° 35   mm B 45   mm 0.51 m / s C 1.03   m / s 700   r / min D 1100   r / min
After optimization, the optimum parameter combination can be obtained. The minimum cutting force was 4.21 N under an inclination angle of 110° (A), an extension length of 40 mm (B), a traveling speed of 0.77 m/s (C), and a rotation speed of 900 r/min (D), thus achieving the optimal cutting effect. The optimal structure and operating parameter combination optimized by simulation tests provide a basis for the subsequent bench test.

5. Bench Test

5.1. Test Materials and Methods

The net curtain of laver was collected from Zhuanghe sea area (Dalian, China). The size of the net curtain was modified to 1.5 m × 1.2 m. The net curtain was soaked in seawater before the bench test.

5.1.1. Harvesting Bench Test

As shown in Figure 10, the test bench was mainly composed of a laver-harvesting device (3), a net curtain (1), and net curtain traveling simulation devices (2 and 4). The harvesting device remained fixed during the test. One side of the net curtain passed through the opposite roller mechanism, and the other side was connected with the hexagonal roller shaft through the upper part of the harvesting device. The rotation speed of the cutter was controlled by a throttle valve, and the net curtain was driven to move in translation by a hexagonal roller shaft, and speed adjusted by a frequency converter to simulate the relative motion of the cutter and the curtain.

5.1.2. Cutting Fracture Detection Test

To explore the rupture state of laver under different operating parameters of the harvester and determine the cutting and rupture proportions, a sampling inspection was performed. The rupture state was predicted by macroscopic observation and then judged by observing the cells at the rupture edge under a microscope. Three net curtains (the net curtain harvested by speed ratio λ3, λ4, and λ6) were selected for each group of tests, and 40 pieces of stubble laver were selected from each net curtain according to the five-point sampling method [32].
The test factors and levels are listed in Table 6. As listed in the table, the speed ratio (λ) refers to the ratio of the rotation speed to the traveling speed, and λ1–6 represents different speed ratios. During the bench test, the speed ratio test of each group was repeated three times.

5.2. Test Evaluation Index

According to the simulation results, the inclination angle of the cutter was 110°, and the extension length was 40 mm. The stubble height was adjusted to 50 mm in the bench test, and the harvested laver was collected and weighed after the test. Additionally, the traveling speed and rotation speed were selected as test factors, and the net recovery rate and loss rate were selected as evaluation indexes.
The evaluation of the net recovery rate is shown in Figure 11, which can be calculated as follows:
y 1 = n 1 n 1 + n 2 × 100 %
where n1 represents the weight of the harvested laver (kg), n2 represents the weight of the laver with an ineligible stubble length on the net curtain (kg).
The evaluation of the loss rate is shown in Figure 12. After the laver harvesting ship was simulated, the laver falling outside the specified area (outside the ship) was regarded as the loss amount, which can be calculated as follows:
y 2 = m 1 m 1 + m 2 × 100 %
where m1 represents the weight of the harvested lever outside the specified area (kg), m2 represents the weight of the harvested lever within the specified area (kg).
Determine the proportion of laver cutting and pulling according to the fracture morphology. Specifically, laver with a flat fracture can be classified into the cutting group (left side of Figure 13); while laver with a serrated fracture can be classified into the rupture group (middle part of Figure 13). As for laver with the two states (right side of Figure 13), the mesh in the figure can be selected as a classifying standard (each mesh was a square with a size of 5 × 5 mm). When the horizontal length of the fracture is larger than 2/3 of the whole length, it can be classified into the cutting group; otherwise, it can be classified into the rupture group.
The cell arrangement at the fracture of laver can be magnified 4 times or 10 times under the microscope with the assistance of TCapture 5.1.1 software. Specifically, laver with a relatively flat fracture and properly arranged marginal cells can be classified into the cutting group, while laver with an uneven fracture and cell fluid exudation can be classified into the rupture group.

5.3. Test Results and Analysis

5.3.1. Influence of the Rotation Speed on Harvesting Effects

As shown in Figure 14, under three tests (No.1, 2, and 3), the net recovery rate of laver increases first and then remains stable with the increase of the rotation speed. The loss rate increases with the increase of the rotation speed. After the rotation speed of the cutter increases from 800 r/min to 1000 r/min, the net recovery rate and loss rate of laver increases by 9.37% and 2.23%, respectively. With the increase in the rotation speed, the missing cutting area along the cutting trajectory decreases, and the relative cutting area and the net recovery rate increase. Due to the viscoelastic feature, laver may undergo plastic deformation before the fracture. With the increase in the rotation speed, the plastic deformation time would shorten, and the impact of cutters on laver would be enhanced. The higher splash weight of laver would also increase the loss rate. When the rotation speed of the cutter increases from 900 r/min to 1000 r/min, the net recovery rate only increases by 0.93%, which is not significant. In addition, the loss rate increases by 1.71%. This can be mainly explained that when the traveling speed is constant, the missing cutting area decreases and gradually remains stable after the rotation speed reaches a certain value. Meanwhile, the high rotation speed could exert a great impact on the laver, the throwing speed would also increase, and the splash weight outside the collection area would increase. Therefore, an appropriate rotation speed can improve the net recovery rate, avoid a larger loss rate, and reduce unnecessary energy consumption.

5.3.2. Influence of the Traveling Speed on Harvesting Effects

As shown in Figure 15, under three tests (No.2, 4, and 5), the net recovery rate gradually decreases with the increase in the traveling speed, while the loss rate increases sharply after a slow increase with the increase in the traveling speed. After the traveling speed increases from 0.51 m/s to 1.03 m/s, the net recovery rate decreases by 2.44%, and the loss rate increases by 4.02%. This can be mainly explained that the missing cutting area along the cutting trajectory increases with the increase in the traveling speed. When the cutting area decreases, the number of laver entering the cutting area increases per unit time. The cutting effect on the back row of the laver in the cutting area becomes worse, and incomplete cutting may be induced. Moreover, a higher traveling speed also contributes to a higher throwing speed. As a result, the net recovery rate is relatively reduced, and the loss rate is relatively increased.

5.3.3. Influence of Speed Ratio on the Proportion of Cutting and Pulling

The microscopic observation results related to fractures are shown in Figure 16. The laver harvested from cutting presents a flat fracture and properly arranged marginal cells (Figure 16a). The laver harvested from ruptures presents cell fluid exudation at the fracture and improperly arranged marginal cells (Figure 16b) The laver with a cutting-ruptured fracture only presents improperly arranged cells around the rupture (Figure 16c).
Figure 17 presents the statistical results of the microscopic examination of laver fractures under different speed ratios. When the speed ratios (λ3, λ4, and λ6) of the three groups changes from large to small, the proportion of harvested laver increases significantly. This may be explained that when the speed ratio is large, the downward vertical speed of the cutter is larger, which makes a great effective contribution to the tensile speed. Hence, the proportion of harvested laver is larger when the speed ratio is small. When the speed ratio is small, the forward horizontal speed of the cutter is larger, which has a great effective contribution to the cutting speed, thus realizing effective cutting. Hence, the proportion of harvested laver increases. After the net recovery rate and loss rate are compared under different operating parameters, the speed ratio (λ4) can be identified to ensure an appropriate net recovery rate and a relatively small loss rate, the laver harvester can achieve a favorable cutting effect.

6. Conclusions

In this study, a laver harvester was designed and manufactured. The relationship between the movement trajectory and the inclination angle of cutters and the stress on the laver was clarified. A rigid-flexible coupling model related to the interaction between the cutting mechanism and the laver was constructed. The Box–Behnken design method was used to simulate the effects of different structural parameters and process parameters on the force of laver cutting. In addition, the bench test and cutting fracture morphological microscope observation of laver harvesting device were carried out.
(1) The results indicated that the inclination angle of the knife was inversely proportional to the cutting stress and directly proportional to the tensile stress. When the inclination angle decreased, the cutting proportion increased; otherwise, the rupture proportion increased. The order of factors affecting the force of laver cutting was cutter revolving speed, knife extension length, knife inclination angle, and traveling velocity. It was confirmed that the cutting force on laver was 4.21 N under the optimal parameter combination with a traveling speed of 0.77 m/s, a rotation speed of 900 r/min, an extension length of 40 mm, and an inclination angle of 110°. The bench test under this parameter, the net recovery rate and the loss rate were 97.45% and 3.38%, respectively, the cutting proportion of laver can reach 77.5%. This result further validated that an excellent harvesting effect can be achieved under this parameter combination.
(2) According to the simulation results, the force required to harvest laver is 4.21~8.13 N, and it is recommended that the force in the actual harvesting process should not be lower than this range. At the same time, it is particularly important to match the traveling speed with the rotation speed of the device in the complex marine environment, in order to the harvesting process needs to have a low loss rate, damage rate, and a high net recovery rate. It is recommended that the traveling speed of the device and tool speed be 0.77 m/s and 900 r/min. This study provides reference for the research methods of flexible crops such as laver.

Author Contributions

Conceptualization, G.M., J.T. and W.L.; methodology, G.M., J.T., W.L., S.N., Z.S., C.Z., Z.S. and Q.Z.; investigation, G.M., J.T., W.L., S.N., H.L. and H.Z.; writing—original draft preparation, J.T.; writing—review and editing, G.M., X.L. and J.T.; supervision, G.Z. and G.M.; funding acquisition, X.L. and G.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Program of China, grant number 2019YFD0900701: the Scientific research project of Liaoning Provincial Department of Education, grant number LJKFZ20220241, the Dalian Science and Technology Innovation Fund Project, grant number 2021JJ12SN33, the Dalian Key Task Project in Scientific and Technological Research, grant number 2021JB11SN035.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The Engineering, Dalian Ocean University, for providing equipment and technical support in the process of the experiment. We are also grateful to our other colleagues such as Guojie Wang, Zhaoxin Jiang and Yusong Liu in the lab for their help during the experiment.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhu, J.Y.; Yan, X.H.; Ding, L.P. Color Atlas of Chinese Laver, 1st ed.; China Agriculture Press: Beijing, China, 2016; pp. 28–33. [Google Scholar]
  2. Zhou, W.; Hu, C.M.; Lu, Q.Q. Germplasm Innovation and Application of Pyropia yezoensis. J. Guangxi Acad. Sci. 2021, 37, 46–52. [Google Scholar]
  3. Agriculture and Rural Fisheries Administration. China Fishery Statistical Yearbook; China Agriculture Press: Beijing, China, 2021; pp. 23–50. [Google Scholar]
  4. Zhang, L.J.; Li, Y.C. ZS160-1.8 type laver harvester. Agric. Food Mach. 1994, 2, 33–34. [Google Scholar]
  5. Huang, X.P. Conveniently Pull Laver Harvester of Planting Frame. C.N. Patent 2,052,841,88U, 8 June 2016. [Google Scholar]
  6. Cai, Q.W. Seaweed Gathters in Clear Head All-in-One Machine. C.N. Patent 2,088,912,16U, 24 May 2019. [Google Scholar]
  7. Wang, Z.C.; Liu, T.Y.; Hu, Y.Z. Study on the Quality Status of Porphyra yezoensis in Jiangsu Province. J. Anhui Agric. Sci. 2018, 46, 207–209. [Google Scholar]
  8. Chen, X.Z.; Qin, X.; Luo, Q.J. Harvesting trials of cultivated Porphyra haitanensis. Sci. Fish Farming 2021, 6, 61–62. [Google Scholar]
  9. Chen, G.F.; Rao, D.Z. Application of mechanized harvesting technology for Porphyra haitanensis. Aquaculture 2013, 6, 37–38. [Google Scholar]
  10. Yada, S.; Chen, J.K.; Toda, M. Analysis of cutting mechanism of the roll type nori harvester. La Mer 1996, 34, 25–32. [Google Scholar]
  11. Kim, O.S.; Min, E.B.; Hwang, D.J. A study on the harvesting process and operating behavior of working ships for farming laver. J. Korean Soc. Fish. Ocean Technol. 2020, 56, 223–229. [Google Scholar] [CrossRef]
  12. Wang, X.N. Current status of research on Laver. China Food Saf. Mag. 2020, 14, 33. [Google Scholar]
  13. He, P.M.; Zhang, Z.Y.; Zhang, X.C. Seaweed Cultivation, 1st ed.; Science Press: Beijing, China, 2018; pp. 174–295. [Google Scholar]
  14. Wang, Y.D.; Zhu, Y.X. Study on Half-floating Cultivation Laver Harvester Solution for Beaches. Fish. Mod. 2004, 5, 29–30. [Google Scholar]
  15. Xu, P.; Zhang, X.C.; Wang, S.J. Reproduction and Development of Major Economic Marine Algae in China, 1st ed.; China Agriculture Press: Beijing, China, 2013; pp. 101–104. [Google Scholar]
  16. Yada, S.; Chen, J.K.; Toda, M. A study on cutting mechanism of the suction type nori harvester. La Mer 1995, 33, 177–184. [Google Scholar]
  17. Cao, W.B.; Sun, W.L.; Niu, C. Combed Safflower Picking Device Based on ANSYS/LS-DYNA. Trans. Chin. Soc. Agric. Mach. 2018, 49, 123–131. [Google Scholar]
  18. Bílek, O.; Kovařík, M. Stress analysis of the rotating circular saw blade. Solid State Phenom. 2017, 261, 259–266. [Google Scholar]
  19. Meng, Y.M.; Wei, J.D.; Wei, J. An ANSYS/LS-DYNA simulation and experimental study of circular saw blade cutting system of mulberry cutting machine. Comput. Electron. Agric. 2019, 157, 38–48. [Google Scholar] [CrossRef]
  20. Zhao, K.J. Research on Dynan1ics Characteristics of Rigid-F1exible Coupling of Camellia Oleifera Picking Mechanism-Fruit Tree. Master’s Thesis, Central South University of Forestry and Technology, Hunan, China, 2020. [Google Scholar]
  21. Dun, G.Q.; Yang, Y.Z.; Li, H.S. Working parameters optimization of soybean serration rotarycutter by ANSYS / LS-DYNA. J. Henan Agric. Univ. 2019, 53, 739–744. [Google Scholar]
  22. Wu, K.; Song, Y.P. Research progress analysis of crop stalk cutting theory and method. Trans. Chin. Soc. Agric. Mach. 2022, 53, 1–20. [Google Scholar]
  23. Zhejiang Marine Fisheries Research Institute. The Culture of Laver in Zhejiang, 1st ed.; Zhejiang Publishing United Group: Zhejiang, China, 1973; pp. 56–63. [Google Scholar]
  24. Zhao, C.N.; Wu, Y.L.; Cai, Z.Q. Experimental Study on Deep-water Rollover Type Cultivation of Porphyra yezoensis. J. Anhui Agric. Sci. 2020, 48, 107–109. [Google Scholar]
  25. Zhang, Y.J. Laver Culture, 1st ed.; Agriculture Press: Beijing, China, 1988; pp. 143–149. [Google Scholar]
  26. Zhang, X.W.; Chen, X.Z.; Chu, M.B. The technology of stem culture of Porphyra haitanensis. Sci. Fish Farming 2021, 3, 63–64. [Google Scholar]
  27. Yu, J. Laver harvester made in Daishan. Fish. Mod. 1982, 1, 61. [Google Scholar]
  28. Zhang, K.F.; Liu, Y.F. A comparative trial of two cultivation patterns of Porphyra haitanensis in Weitou Bay. J. Aquac. 2020, 10, 35–37. [Google Scholar]
  29. Lu, W.; Li, X.C.; Zhang, G.C.; Tang, J.H.; Ni, S.; Zhang, H.B.; Zhang, Q.; Zhai, Y.L.; Mu, G. Research on biomechanical properties of laver (Porphyra yezoensis Ueda) for mechanical harvesting and postharvest transportation. AgriEngineering 2022, 41, 48–66. [Google Scholar] [CrossRef]
  30. Li, B.F. Agromechanics, 1st ed.; China Agricultural Press: Beijing, China, 2018. [Google Scholar]
  31. Li, X.D. Research on Stalk Cutting Device of Corn and Stalk Combine Harvester. Master’s Thesis, Chinese Academy of Agricultural Mechanization Sciences Group Co., Ltd., Beijing, China, 2015. [Google Scholar]
  32. Zhang, X.G.; Du, L.X.; Jia, S.Z. Study on the control and analysis of Ips Subelongatus. Xiandai Nongye Keji 2016, 11, 157+164. [Google Scholar]
Figure 1. Laver cultivation method.(a) Full floating raft; (b) Pillar raft; (c) Semi-floating raft.
Figure 1. Laver cultivation method.(a) Full floating raft; (b) Pillar raft; (c) Semi-floating raft.
Agriculture 13 00361 g001
Figure 2. Structure diagram of laver harvester. 1. Arc cover 2. Connection support plate 3. Knife 4. Cutter shaft 5. Protective grid 6. Rotary drum 7. Connector 8. Roller 9. Rack.
Figure 2. Structure diagram of laver harvester. 1. Arc cover 2. Connection support plate 3. Knife 4. Cutter shaft 5. Protective grid 6. Rotary drum 7. Connector 8. Roller 9. Rack.
Agriculture 13 00361 g002
Figure 3. Working principle diagram of laver harvester. 1. Net curtain 2. Arc cover 3. Knife 4. Protective grid 5. Laver 6. Roller 7. Rack 8. Connection support plate 9. Rotation direction of cutters 10. Movement direction of the net curtain.
Figure 3. Working principle diagram of laver harvester. 1. Net curtain 2. Arc cover 3. Knife 4. Protective grid 5. Laver 6. Roller 7. Rack 8. Connection support plate 9. Rotation direction of cutters 10. Movement direction of the net curtain.
Agriculture 13 00361 g003
Figure 4. Knife movement trajectory diagram.
Figure 4. Knife movement trajectory diagram.
Agriculture 13 00361 g004
Figure 5. Force analysis of laver. A. Inclination angle B. Extension length C. Rotation speed D. Traveling speed.
Figure 5. Force analysis of laver. A. Inclination angle B. Extension length C. Rotation speed D. Traveling speed.
Agriculture 13 00361 g005
Figure 6. Simulation Model Meshing.
Figure 6. Simulation Model Meshing.
Agriculture 13 00361 g006
Figure 7. Simulated harvesting process.
Figure 7. Simulated harvesting process.
Agriculture 13 00361 g007
Figure 8. The force curve of laver in the X and Y directions. (a) The force curve of laver in the X direction; (b) The force curve of laver in the Y direction.
Figure 8. The force curve of laver in the X and Y directions. (a) The force curve of laver in the X direction; (b) The force curve of laver in the Y direction.
Agriculture 13 00361 g008
Figure 9. Response surface of the interaction of various factors on the cutting force of laver. (a) The influence of the inclination angle and extension length on the cutting force of laver; (b) The influence of the extension length and traveling speed on the cutting force of laver.
Figure 9. Response surface of the interaction of various factors on the cutting force of laver. (a) The influence of the inclination angle and extension length on the cutting force of laver; (b) The influence of the extension length and traveling speed on the cutting force of laver.
Agriculture 13 00361 g009
Figure 10. Schematic diagram of test bench construction. 1. Net curtain 2. Hexagonal roller shaft 3. Harvesting device 4. Opposite roller mechanism.
Figure 10. Schematic diagram of test bench construction. 1. Net curtain 2. Hexagonal roller shaft 3. Harvesting device 4. Opposite roller mechanism.
Agriculture 13 00361 g010
Figure 11. Schematic diagram of net recovery rate evaluation principle.
Figure 11. Schematic diagram of net recovery rate evaluation principle.
Agriculture 13 00361 g011
Figure 12. Schematic diagram of loss ratio evaluation.
Figure 12. Schematic diagram of loss ratio evaluation.
Agriculture 13 00361 g012
Figure 13. Macroscopical observation method.
Figure 13. Macroscopical observation method.
Agriculture 13 00361 g013
Figure 14. Technology parameters of prototype.
Figure 14. Technology parameters of prototype.
Agriculture 13 00361 g014
Figure 15. Technology parameters of prototype.
Figure 15. Technology parameters of prototype.
Agriculture 13 00361 g015
Figure 16. Electron microscope scan of sample incision. (a) Marginal cells of the cutting fracture laver; (b) Marginal cells of the ruptured fracture lave; (c) Marginal cells of the cutting-ruptured fracture lave.
Figure 16. Electron microscope scan of sample incision. (a) Marginal cells of the cutting fracture laver; (b) Marginal cells of the ruptured fracture lave; (c) Marginal cells of the cutting-ruptured fracture lave.
Agriculture 13 00361 g016
Figure 17. Proportion diagram of fracture state of laver.
Figure 17. Proportion diagram of fracture state of laver.
Agriculture 13 00361 g017
Table 1. The main physical parameters of laver.
Table 1. The main physical parameters of laver.
ParametersValues
Frond overall length (mm)78.3~236.65
Frond head width (mm)4.2~14.9
Frond middle width (mm)10~21.4
Frond root width (mm)8.8~19.6
Frond head thickness (mm)0.013~0.028
Frond middle thickness (mm)0.02~0.035
Frond root thickness (mm)0.016~0.033
Density (wet) (g∙cm−3)0.208~0.596
Density (dry) (g∙cm−3)0.571~1.684
Table 2. Material parameters of model.
Table 2. Material parameters of model.
ParameterLaver [29]Cutter
Poisson’s ratio 0.330.33
Shear modulus (Pa) 3.26 × 107-
Elastic modulus (Pa) 2.54 × 1062.00 × 1011
Tensile strength (MPa) 0.95-
Shear strength (MPa) 2.33-
Density (kg·m−3) 0.967850
Table 3. Level of test factors for response surface.
Table 3. Level of test factors for response surface.
Level Factor
Inclination Angle A (°) Extension Length B (mm)Traveling Speed C (m∙s−1)Rotation Speed D (r∙min−1)
−1100350.51700
0110400.77900
1120451.031100
Table 4. Level of test factors for response surface.
Table 4. Level of test factors for response surface.
Test No. Inclination Angle A (°) Extension Length B (mm) Traveling Speed C (m∙s−1) Rotation Speed D (r∙min−1)Cutting Force (N)
1110.0040.000.77900.004.96
2110.0035.000.51900.004.66
3110.0040.000.77900.005.41
4100.0045.000.77900.005.04
5110.0045.000.51900.004.42
6120.0045.000.77900.006.64
7100.0040.000.771100.006.00
8100.0040.000.51900.004.24
9110.0045.000.771100.005.48
10110.0035.000.771100.005.88
11110.0040.000.77900.005.16
12120.0040.000.771100.008.13
13110.0045.001.03900.006.14
14120.0040.000.77700.004.17
15110.0040.000.77900.005.00
16110.0035.001.03900.003.31
17110.0040.001.03700.002.79
18120.0040.001.03900.003.89
19110.0035.000.77700.003.19
20120.0040.000.51900.004.90
21110.0040.000.51700.003.16
22100.0040.001.03900.004.34
23110.0040.000.77900.004.21
24100.0040.000.77700.003.34
25120.0035.000.77900.003.77
26100.0035.000.77900.005.43
27110.0045.000.77700.003.63
28110.0040.001.031100.005.83
29110.0040.000.511100.005.45
Table 5. ANOVA results for response surface tests.
Table 5. ANOVA results for response surface tests.
IndexSourceSquaresdfMean SquareFp
YModel 34.86142.496.520.0006 **
A0.8110.812.110.1684
B2.1812.185.690.0317 *
C0.02310.0230.0610.8081
D22.66122.6659.30<0.0001 **
AB2.6612.666.950.0195 *
AC0.3110.310.810.3845
AD0.4210.421.110.3108
BC2.3612.366.170.0263 *
BD0.1810.180.460.5079
CD0.1410.140.370.5538
A20.4710.471.230.2858
B20.00210.0020.0060.9414
C22.1512.155.630.0326 *
D20.04810.0480.120.7294
Residual error 5.35140.38
Lack of fit 4.54100.452.260.2250
Error 0.8140.20Not significant
Sum 40.2128
Note: p < 0.01 (Extremely significant, **); p < 0.05 (Significant, *), the same below.
Table 6. Bench test factors and levels.
Table 6. Bench test factors and levels.
Test No.Rotation Speed (r∙min−1)Traveling Speed (m∙s−1)Speed Ratio (λ)
18000.51λ1
2900λ2
31000λ3
49000.77λ4
51.03λ5
68001.03λ6
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tang, J.; Li, X.; Zhang, G.; Lu, W.; Ni, S.; Sun, Z.; Li, H.; Zhao, C.; Zhang, H.; Zhang, Q.; et al. An ANSYS/LS-DYNA Simulation and Experimental Study of Sectional Hob Type Laver Harvesting Device. Agriculture 2023, 13, 361. https://doi.org/10.3390/agriculture13020361

AMA Style

Tang J, Li X, Zhang G, Lu W, Ni S, Sun Z, Li H, Zhao C, Zhang H, Zhang Q, et al. An ANSYS/LS-DYNA Simulation and Experimental Study of Sectional Hob Type Laver Harvesting Device. Agriculture. 2023; 13(2):361. https://doi.org/10.3390/agriculture13020361

Chicago/Turabian Style

Tang, Jiahong, Xiuchen Li, Guochen Zhang, Wei Lu, Shang Ni, Zhenyin Sun, Haidong Li, Cheng Zhao, Hanbing Zhang, Qian Zhang, and et al. 2023. "An ANSYS/LS-DYNA Simulation and Experimental Study of Sectional Hob Type Laver Harvesting Device" Agriculture 13, no. 2: 361. https://doi.org/10.3390/agriculture13020361

APA Style

Tang, J., Li, X., Zhang, G., Lu, W., Ni, S., Sun, Z., Li, H., Zhao, C., Zhang, H., Zhang, Q., & Mu, G. (2023). An ANSYS/LS-DYNA Simulation and Experimental Study of Sectional Hob Type Laver Harvesting Device. Agriculture, 13(2), 361. https://doi.org/10.3390/agriculture13020361

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop