Next Article in Journal
Can Visual Aesthetic Components and Acceptance Be Traced Back to Forest Structure?
Previous Article in Journal
Studies into Fungal Decay of Wood in Ground Contact—Part 2: Development of a Dose–Response Model to Predict Decay Rate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Finite Element Method Simulations and Experiments of Detachments of Lycium barbarum L.

1
College of Mechanical and Electronic Engineering, Northwest A&F University, Yangling 712100, China
2
Graduate School of Bioagricultural Sciences, Nagoya University, Furo-Cho, Chikusa, Nagoya 464-8601, Japan
*
Authors to whom correspondence should be addressed.
These authors equally contributed to this work.
Forests 2021, 12(6), 699; https://doi.org/10.3390/f12060699
Submission received: 20 April 2021 / Revised: 26 May 2021 / Accepted: 27 May 2021 / Published: 28 May 2021
(This article belongs to the Section Wood Science and Forest Products)

Abstract

:
When harvesting Lycium barbarum L., excess amounts of detachments of the half-ripe fruit, unripe fruit, flowers, and leaves significantly affect the yield and adversely affect the subsequent processing, such as drying and grading. Finite element method (FEM) simulations and experiments of detachments were performed to harvest more ripe fruit and less half-ripe fruit, unripe fruit, flowers, and leaves. Three-dimensional (3D) models of the ripe fruit, half-ripe fruit, unripe fruit, flowers, leaves, fruit calyxes (flower calyx), fruit stems (flower stem), and branches were constructed using a 3D scanner, and material mechanics models of the above parts were established based on physical tests with universal testing machines. Detachment simulations and experiments of the ripe fruit, half-ripe fruit, unripe fruit, flowers, and leaves were performed to determine the detachment mechanisms and sequences. The detachment forces of each set of two parts were obtained. The field experiments showed that the detachment force between the fruit and calyx of ripe fruit was the lowest value of these forces, and only the ripe fruit was the first to detach from the calyx when harvesting. The results provided data support on the mechanics properties of wood and the optimization basis for the harvesting method of L. barbarum.

1. Introduction

Lycium barbarum L. is a solanaceae Lycium deciduous shrub with high added value, and the fruit has high profitability [1,2,3,4]. Its ripe fruit after being dried was widely favored by more and more people due to the ingredients, which are beneficial to health, such as a group of L. barbarum polysaccharides, which have neuroprotective, anti-oxidant, and immunomodulation functions [3,5]. However, the price of the fruit remains high because L. barbarum is harvested manually [4,6], and the harvesting cost accounts for half of the total cost of production, currently. With the increase of the planting area of L. barbarum, the problem of labor shortage for the fruit harvesting is becoming increasingly prominent [6]. Thus, it is urgent to design and develop an L. barbarum harvester [7,8,9,10]. Scholars have developed some prototype harvesters of L. barbarum [11,12,13,14]. It indicated that half-ripe and unripe fruit were easy to be harvested by mistake [6,11]. Meanwhile, excess amounts of detachments of the half-ripe fruit, unripe fruit, and flowers substantially affect the yield. Hence, it is essential to analyze the detachment mechanisms of L. barbarum to achieve the high-efficiency and low-loss harvesting performance. Furthermore, there is also little research on the wood mechanics properties of L. barbarum. This leads to a lack of the key parameters in the analysis and affects the accuracy of the results. In addition, impurities, such as the half-ripe fruit, unripe fruit, flowers, and leaves, significantly affect the subsequent processing, such as the drying and grading [15]. Therefore, it is crucial to investigate the detachment sequences of L. barbarum to harvest more ripe fruit and less half-ripe fruit, unripe fruit, flowers, and leaves. The essence of detachment forces was the strength of plant attachments. The branch-to-stem diameter ratio affected the strength of plant attachments [16,17]. In addition, the picking patterns have a significant effect on the detachment forces [18,19]. Thus, dividing the ripe fruit, half-ripe fruit, unripe fruit, flowers, and leaves of L. barbarum into different plant parts according to the materials was more accurate to analyze the detachment mechanisms and sequences of L. barbarum. With the development of the reverse engineering technology, it was easy to obtain accurate three-dimensional (3D) models. Modeling each plant part using the reverse engineering technology was more conducive to simulation analyses accurately. Additionally, the material mechanics parameters of different plant parts were different. It was necessary to establish material mechanics models through physical tests [20,21]. On this basis, it was expected to access the detachment mechanisms and sequences of L. barbarum by performing finite element method (FEM) simulations and experiments.
In this study, to harvest more ripe fruit and less half-ripe fruit, unripe fruit, flowers, and leaves, 3D models of the ripe fruit, half-ripe fruit, unripe fruit, flower, leaf, fruit calyxes (flower calyx), fruit stems (flower stem), and branches were constructed using a 3D scanner, and material mechanics models of the above parts were established based on physical tests with universal testing machines. FEM simulations and experiments of detachments of the ripe fruit, half-ripe fruit, unripe fruit, flowers, and leaves were performed to determine the detachment mechanisms and sequences of L. barbarum. The detachment forces of each set of two parts were also obtained. Such findings provided data support on the mechanics properties of wood and the optimization basis for the harvesting method of L. barbarum.

2. Materials and Methods

2.1. Reverse Engineering Modeling Using the 3D Scanner

The more realistic the 3D models are, the more accurate the results of subsequent detachment simulations based on FEM will be. Reverse engineering modeling was used to obtain accurate 3D models using the 3D scanner (type: HandySCAN 700; precision: 0.03 mm; manufactured by Creaform Shanghai Co., Ltd., Shanghai, China). The material mechanics properties differ for different plant parts, affecting the results of detachment simulations and experiments. The ripe fruit, half-ripe fruit, and unripe fruit could be divided into four parts (i.e., the fruit, fruit calyx, fruit stem, and branch), as shown in Figure 1a–c. As shown in Figure 1d, the flower could be divided into four parts (i.e., the flower, flower calyx, flower stem, and branch). As shown in Figure 1e, the leaf could be divided into two parts (i.e., the leaf and branch). The creations of the 3D models of the ripe fruit, half-ripe fruit, unripe fruit, flowers, leaves, fruit calyxes (flower calyx), fruit stems (flower stem), and branches included calibrating, configuring, attaching marking points, scanning the orientation point, scanning the sample, meshing, and rendering in the software (type: Vxelements 3D software version-6.1; manufactured by Creaform Inc., Lévis, Canada). To better display the 3D models, their surfaces were rendered with the color blue. There were some defects in the surfaces due to the shade and illumination during the scanning process. The defects were corrected based on the morphology method and rendered with the color yellow in the 3D models.

2.2. Physical Tests Using Universal Testing Machines

The density, elastic modulus, and Poisson ratio are essential parameters for the detachment simulations. The densities of the plant parts were calculated by measuring the masses and volumes of the corresponding parts. Two electronic scales (type: Professional digital jewelry scale 8028-series; measurement range: 0–100 g; precision: 0.001 g; manufactured by Shenzhen Diheng Electronics Co., Ltd., Shenzhen, China and type: SA-120; measurement range: 0–120 g; precision: 0.0002 g; manufactured by Shinko Denshi Co., Ltd., Tokyo, Japan) were used to measure the masses. Two electronic vernier calipers (type: AIRAJ second-generation product; measurement range: 0–300 mm; precision: 0.01 mm; manufactured by Qingdao Yigou Hardware Tools Co., Ltd., Qingdao, China and type: CD-15APX; measurement range: 0–150 mm; precision: 0.01 mm; manufactured by Mitutoyo Corporation, Kanagawa, Japan) were used to obtain the measurements. The volumes could be calculated according to the obtained measurements. The Poisson ratios were obtained from relevant references. However, the elastic moduli should be calculated from the stress–strain curves obtained in the physical tests. Hence, two electronic universal testing machines (type: DDL10; max testing force: 10 kN; manufactured by Sinotest Equipment Co., Ltd., Changchun, China and type: AG-1; max testing force: 100 kN; manufactured by Shimadzu Corporation, Kyoto, Japan) were used to perform all physical tests. The cutting surfaces of the samples were polished before the tests. Meanwhile, in order to get more accurate stress–strain curves, compression or tensile tests should be selected according to the samples. Therefore, physical tests of the ripe fruit, half-ripe fruit, unripe fruit, flower, leaf, fruit calyxes (flower calyx), fruit stems (flower stem), and branches were conducted.

2.3. Establishing the Material Mechanics Models

2.3.1. Densities

The more realistic the hypothetical shapes are, the more accurate the densities are when calculating the volumes of the plant parts. Since the shapes of the ripe fruit, half-ripe fruit, and unripe fruit are similar to an ellipsoid, the volumes of the ripe fruit, half-ripe fruit, and unripe fruit were calculated using the following Equation (1) of an ellipsoid.
V = 4 3 π a b c
where V is the volume of the sample, mm3; a, b, and c are half of the length of each axis of the ellipsoid-shaped sample, mm.
The flower is cone-shaped; thus, the volume of the flower was calculated using the following Equation (2) of a cone.
V = 1 3 π R 2 h
where R is the bottom radius of the cone-shaped sample, mm; h is the height of the cone-shaped sample, mm.
The leaf is cuboid-shaped; thus, the volume of the leaf was calculated using the following Equation (3) of a cuboid.
V = a b h
where a is the length of the cuboid-shaped sample, mm; b is the width of the cuboid-shaped sample, mm; h is the height of the cuboid-shaped sample, mm.
The calyxes of the ripe fruit, half-ripe fruit, unripe fruit, and flowers are cone-shaped; thus, the volumes of these parts were calculated using the Equation (2) of a cone. The stems of the ripe fruit, half-ripe fruit, unripe fruit, and flowers and the branches of the ripe fruit, half-ripe fruit, unripe fruit, flowers, and leaves are similar to a cylinder. The volumes of these parts were calculated using the following Equation (4) of a cylinder.
V = π R 2 h
where R is the bottom radius of the cylinder-shaped sample, mm; h is the height of the cylinder-shaped sample, mm.
The densities of all plant parts were calculated using the following Equation (5). The results are listed in Table 1.
ρ = m V
where ρ is the density of the sample, t/mm3; m is the mass of the sample, t.

2.3.2. Elastic Moduli

The stresses and strains were calculated based on the following Equations (6) and (7) and the references [3,22,23]:
σ = F A 0
ε = Δ l l 0
where σ is the stress of the sample, MPa; F is the load of the test, N; A0 is the cross-sectional area of the sample, mm2; ε is the strain of the sample, mm·mm−1; ∆l is the change in the length of the sample, mm; l0 is the nominal length of the sample, mm.
The elastic modulus is defined as the following Equation (8) based on the references [3,22,23]:
E = σ ε
where E is the elastic modulus of the sample, MPa.
The load of the test and the change in the length of the sample were obtained from tests with the universal testing machines. The nominal length of the sample was measured using the electronic vernier calipers. However, the cross-sectional area of the sample should be determined based on the sample shapes. The mean values of the largest and smallest ellipse area of the ripe fruit, half-ripe fruit, and unripe fruit were considered the cross-sectional area of the sample. The mean value of the largest and smallest round area of the flower was considered the cross-sectional area of the sample. The mean value of the largest and smallest rectangle area of the leaf was considered the cross-sectional area of the sample. The mean values of the largest and smallest round area of the calyxes of the ripe fruit, half-ripe fruit, unripe fruit, and flowers were considered the cross-sectional area of the sample. The mean values of the largest and smallest round area of the stems of the ripe fruit, half-ripe fruit, unripe fruit, and flowers and the branches of the ripe fruit, half-ripe fruit, unripe fruit, flowers, and leaves were considered the cross-sectional area of the sample. Therefore, the elastic moduli of all plant parts were calculated based on the stress–strain curves and are listed in Table 1.

2.3.3. Poisson Ratio

The Poisson ratio is the ratio of the absolute value of the transverse positive strain to the axial positive strain when a material is under a unidirectional compression or tension test, also known as a transverse deformation coefficient, which is an elastic constant reflecting the transverse deformation of a material [24]. According to the previous studies, the Poisson ratios of many materials have been measured. Since the Poisson ratio has little effect on the simulation result, the usual practice is to assume the Poisson ratio of the material according to the Poisson ratio of the similar material obtained previously. Therefore, the Poisson ratios of all plant parts were assumed and obtained from the references [3,25,26,27,28,29,30] and are listed in Table 1.

2.4. Detachment Simulations Based on FEM

After obtaining the material mechanics parameters, the FEM simulations results could be analyzed to access the detachment mechanisms of the ripe fruit, half-ripe fruit, unripe fruit, flowers, and leaves theoretically. The simulations were conducted using the Abaqus software. The 3D models were inputted into the Part module. The material mechanics parameters (i.e., the density, elastic modulus, and Poisson ratio) of the plant parts (Table 1) were inputted into the Property module. The velocity was chosen in the predefined field type, and the definition was translational only in the Load module. Each branch was selected as the embedded option, and the mass centers of the above parts were assigned a vertical downward velocity in this module. The 3D models of the ripe fruit, half-ripe fruit, unripe fruit, flowers, and leaves were meshed in the Mesh module, as shown in Figure 2a, Figure 3a, Figure 4a, Figure 5a, and Figure 6a, respectively. Many grids were used to make up the models. Every grid was green, and the boundaries between the grids were black.

3. Results

3.1. Analyses of Detachment Simulations

The analyses and post-processing of the ripe fruit, half-ripe fruit, unripe fruit, flowers, and leaves were performed in the Visualization module. The simulation results of the ripe fruit, half-ripe fruit, unripe fruit, flowers, and leaves are shown in Figure 2b, Figure 3b, Figure 4b, Figure 5b, and Figure 6b, respectively. The magnitude was the relative displacement comparing to the initial state and selected as the display option of simulation results. The greater the magnitude value, the greater the relative displacement. The relative displacements of each part during the detachment progress were displayed in real time. The first fracture of the ripe fruit occurred between the fruit and calyx. That of the half-ripe fruit occurred between the stem and branch mostly. That of the unripe fruit occurred between the stem and branch. That of the flower occurred between the flower and calyx mostly. That of the leaf occurred between the leaf and branch.

3.2. Verification Using Field Experiments

Field experiments were conducted to verify the accuracy of the detachment simulations results of the ripe fruit, half-ripe fruit, unripe fruit, flowers, and leaves. A digital display tensile force meter (type: HP-30; measurement range: 0–30 N; precision: 0.01 N; manufactured by Yueqing Aidebao Instrument Co., Ltd., Yueqing, China) was used to measure the detachment forces. Each experiment was repeated ten times to eliminate random errors. As shown in Table 2, the detachment forces of each two parts were obtained, and the unit was N. The mean values and standard deviations were also obtained. The standard deviations reflected the degree of dispersion of a dataset. According to the standard deviations, the mean values of each force were valid, and the force value was evenly distributed for the same detachment force. To better display the detachment sequences, representative samples with the similar time dimension were selected for the comparison. The detachment experiments results of the ripe fruit, half-ripe fruit, and unripe fruit are shown in Figure 7 and that of the flowers and leaves are shown in Figure 8.

4. Discussion

Based on the field experiments results, as shown in Figure 7; Figure 8, the detachment mechanisms and sequences of the ripe fruit, half-ripe fruit, unripe fruit, flowers, and leaves were determined. To be specific, the ripe fruit fractured first between the fruit and calyx, second between the stem and branch, and the stem fractured last. Mostly, the half-ripe fruit fractured first between the stem and branch, second between the fruit and calyx, and the stem fractured last. The unripe fruit fractured first between the stem and branch, second between the fruit and calyx, and the stem fractured last. The flower fractured first between the flower and calyx mostly, the stem fractured first in some cases, and it fractured first between the stem and branch least often. The leaves fractured between the leaf and branch. The field experiments results were consistent with the simulations results. This also showed that the established material mechanics models of the plant parts were accurate.
As shown in Table 2, it indicated that the detachment force between the fruit and calyx of ripe fruit was the lowest value of these forces. The detachment forces of the ripe fruit, half-ripe fruit, unripe fruit, flowers, and leaves were also distributed stably. Therefore, only the ripe fruit was the first to detach from the calyx when harvesting. Scholars who have designed the L. barbarum harvesters also proposed that impurities, such as unripe fruit and flowers, affected the comprehensive harvesting performance [1,2,4,6,11,15]. Based on the obtained results, it was expected to design and develop a L. barbarum harvester, which could harvest more ripe fruit and less half-ripe fruit, unripe fruit, flowers, and leaves. When designing the harvester, the component forces of the harvester mechanism should be integrated to form the resultant force, firstly. Then, the resultant force was compared with the detachment forces. As long as the resultant force was greater than the detachment force between the fruit and calyx of ripe fruit and less than the other detachment forces, it was expected to achieve the high-efficiency and low-loss harvesting performance.
Furthermore, codominant stems were easy to split apart compared with branches that were small relative to stem size [16]. It illustrated that the detachment forces of the ripe fruit, half-ripe fruit, unripe fruit, flowers, and leaves were different. The strength of attachment was the major factor that affected the detachment force [16,17]. Hu et al. proposed that studies of the detachment force were a basis to design a fruit harvester. It revealed that the detachment forces of different forms were different [31]. It was consistent with the obtained results in this study. The results in this study provided data support on mechanics properties of wood and the optimization basis for the harvesting method of L. barbarum. In addition, it indicated that the tension was the dominant factor during the detachment process [32]. Based on the tension, the detachment mechanisms and sequences of L. barbarum were determined, and the detachment forces of each two parts were also obtained in this study. It also showed that harvesting patterns were also critical for the efficient harvesting from references [18,19,32]. Especially, there were four picking methods, i.e., the horizontal pull, vertical tension, bending, and twisting [32]. Therefore, based on the obtained results in this study, further study is required to determine which harvesting patterns (i.e., the pull, bend, and torsion) are more efficient for harvesting L. barbarum.

5. Conclusions

In this study, to harvest more ripe fruit and less half-ripe fruit, unripe fruit, flowers, and leaves, FEM simulations and experiments of detachments of L. barbarum were conducted to determine the detachment mechanisms and sequences. Firstly, 3D models of the ripe fruit, half-ripe fruit, unripe fruit, flowers, leaves, fruit calyxes (flower calyx), fruit stems (flower stem), and branches were constructed using the 3D scanner. Then, material mechanics models of the above parts were established based on physical tests with universal testing machines. Lastly, detachment simulations and experiments of the ripe fruit, half-ripe fruit, unripe fruit, flowers, and leaves were performed to determine the detachment mechanisms and sequences. The detachment forces of each set of two parts were also obtained. The field experiments showed that the detachment force between the fruit and calyx of ripe fruit was the lowest value of these forces, and only the ripe fruit was the first to detach from the calyx when harvesting. According to the obtained results, it was expected to design and develop a L. barbarum harvester, which could harvest more ripe fruit and less half-ripe fruit, unripe fruit, flowers, and leaves. When designing the harvester, the component forces of the harvester mechanism should be integrated to form the resultant force, firstly. Then, the resultant force was compared with the detachment forces. As long as the resultant force was greater than the detachment force between the fruit and calyx of ripe fruit and less than the other detachment forces, it was expected to achieve the high-efficiency and low-loss harvesting performance. The results could provide data support on the mechanics properties of wood and the optimization basis for the harvesting method of L. barbarum. Further study is required to determine which harvesting patterns (i.e., the pull, bend, and torsion) are more efficient for harvesting L. barbarum based on the knowledge earned in this study.

Author Contributions

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

Funding

This research and APC were funded by the National Key Research and Development Program of China, grant number 2018YFD0701102. The authors also appreciate the financial support provided by the China Scholarship Council.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xu, L.; Chen, J.; Wu, G.; Yuan, Q.; Ma, S.; Yu, C.; Duan, Z.; Xing, J.; Liu, X. Design and operating parameter optimization of comb brush vibratory harvesting device for wolfberry. Trans. CSAE 2018, 34, 75–82. [Google Scholar] [CrossRef]
  2. Zhang, Z.; Xiao, H.; Ding, W.; Mei, S. Machanism simulation analysis and prototype experiment of Lycium barbarum harvest by vibration mode. Trans. CSAE 2015, 31, 20–28. [Google Scholar] [CrossRef]
  3. Zhao, J.; Sugirbay, A.; Chen, Y.; Zhang, S.; Liu, F.; Bu, L.; Chen, Y.; Wang, Z.; Chen, J. FEM explicit dynamics simulation and NIR hyperspectral reflectance imaging for determination of impact bruises of Lycium barbarum L. Postharvest Biol. Technol. 2019, 155, 102–110. [Google Scholar] [CrossRef]
  4. Wang, Y.; Chen, Y.; Han, B.; Chen, J. Research on laws of wolfberry dropping based on high-speed camera. J. Agric. Mech. Res. 2018, 40, 166–170. [Google Scholar] [CrossRef]
  5. Amagase, H.; Farnsworth, N.R. A review of botanical characteristics, phytochemistry, clinical relevance in efficacy and safety of Lycium barbarum fruit (Goji). Food Res. Int. 2011, 44, 1702–1717. [Google Scholar] [CrossRef]
  6. Chen, J.; Zhao, J.; Chen, Y.; Bu, L.; Hu, G.; Zhang, E. Design and experiment on vibrating and comb brushing harvester for Lycium barbarum. Trans. CSAM 2019, 50, 152–161. [Google Scholar] [CrossRef]
  7. So, J.D. Vibration characteristics of boxthorn (Lycium chinense Mill) branches. Appl. Eng. Agric. 2001, 17, 755–760. [Google Scholar] [CrossRef]
  8. So, J.D. Vibratory harvesting machine for boxthorn (Lycium chinense Mill) berries. Trans. ASAE 2003, 46, 211–221. [Google Scholar] [CrossRef]
  9. Zhang, W.; Li, Z.; Tan, Y.; Li, W. Optimal design and experiment on variable pacing combing brush picking device for Lycium barbarum. Trans. CSAM 2018, 49, 83–90. [Google Scholar] [CrossRef]
  10. Zhang, W.; Zhang, M.; Zhang, J.; Li, W. Design and experiment of vibrating wolfberry harvester. Trans. CSAM 2018, 49, 97–102. [Google Scholar] [CrossRef]
  11. Zhao, J.; Chen, Y.; Wang, Y.; Chen, J. Experimental research on parameter optimization of portable vibrating and harvesting device of Chinese wolfberry. J. Agric. Mech. Res. 2019, 41, 176–182. [Google Scholar] [CrossRef]
  12. Wang, Y. Research on Key Technology of Wolfberry Vibration Harvest. Master’s Thesis, Northwest A&F University, Yangling, China, 2018. [Google Scholar]
  13. Zhou, B.; He, J. Design of simulate hand wolfberry picking machine. Trans. CSAE 2010, 26 (Supp. 1), 13–17. [Google Scholar] [CrossRef]
  14. He, M.; Kan, Z.; Li, C.; Wang, L.; Yang, L.; Wang, Z. Mechanism analysis and experiment on vibration harvesting of wolfberry. Trans. CSAE 2017, 33, 47–53. [Google Scholar] [CrossRef]
  15. Zhao, J.; Sugirbay, A.; Liu, F.; Chen, Y.; Hu, G.; Zhang, E.; Chen, J. Parameter optimization of winnowing equipment for machineharvested Lycium barbarum L. Span. J. Agric. Res. 2019, 17, e0203. [Google Scholar] [CrossRef]
  16. Gilman, E.F. Branch-to-stem diameter ratio affects strength of attachment. J. Arboric. 2003, 29, 291–293. [Google Scholar]
  17. Kane, B.; Farrell, R.; Zedaker, S.M.; Loferski, J.R.; Smith, D.W. Failure mode and prediction of the strength of branch attachments. Arboric. Urban For. 2008, 34, 308–316. Available online: https://www.researchgate.net/profile/Brian-Kane/publication/238757178_Failure_Mode_and_Prediction_of_the_Strength_of_Branch_Attachments/links/00b7d536b533e75fc5000000/Failure-Mode-and-Prediction-of-the-Strength-of-Branch-Attachments.pdf (accessed on 25 May 2021).
  18. Li, J.; Karkee, M.; Zhang, Q.; Xiao, K.; Feng, T. Characterizing apple picking patterns for robotic harvesting. Comput. Electron. Agric. 2016, 127, 633–640. [Google Scholar] [CrossRef] [Green Version]
  19. Davidson, J.; Silwal, A.; Karkee, M.; Mo, C.; Zhang, Q. Hand-picking dynamic analysis for undersensed robotic apple harvesting. Trans. ASABE 2016, 59, 745–758. [Google Scholar] [CrossRef]
  20. Slater, D.; Ennos, A.R. Determining the mechanical properties of hazel forks by testing their component parts. Trees 2013, 27, 1514–1524. [Google Scholar] [CrossRef]
  21. Viet, D.D.; Ma, T.; Inagaki, T.; Kim, N.T.; Chi, N.Q.; Tsuchikawa, S. Physical and Mechanical Properties of Fast Growing Polyploid Acacia Hybrids (A. auriculiformis × A. mangium) from Vietnam. Forests 2020, 11, 717. [Google Scholar] [CrossRef]
  22. Celik, H.K. Determination of bruise susceptibility of pears (Ankara variety) to impact load by means of FEM-based explicit dynamics simulation. Postharvest Biol. Tec. 2017, 128, 83–97. [Google Scholar] [CrossRef]
  23. Shen, C.; Chen, Q.; Li, X.; Zhang, B.; Huang, J.; Tian, K. Test and analysis of axial compressive mechanical properties for ramie stalk. Acta Agric. Univ. Zhejiangensis 2016, 28, 688–692. [Google Scholar] [CrossRef]
  24. The Poisson’s Ratio. Available online: https://baike.baidu.com/item/%E6%B3%8A%E6%9D%BE%E6%AF%94/5920115?fr=aladdin (accessed on 25 May 2021).
  25. Shen, C.; Li, X.; Tian, K.; Zhang, B.; Huang, J.; Chen, Q. Experimental analysis on mechanical model of ramie stalk. Trans. CSAE 2015, 31, 26–33. [Google Scholar] [CrossRef]
  26. Wang, Y. Abaqus Analysis User’s Guide: Materials Volume, 1st ed.; China Machine Press: Beijing, China, 2018; pp. 27–28. [Google Scholar]
  27. Shi, N. Peeling Method of Cotton Stalk and Development of Rubbing Peeling Machine. Ph.D. Thesis, Northwest A&F University, Yangling, China, 2018. [Google Scholar]
  28. Wu, J.; Huang, Y.; Wang, Y.; Wang, W. Study on axial compression properties of cotton stalks. J. Agric. Mech. Res. 2004, 26, 148–149, 152. [Google Scholar] [CrossRef]
  29. Wang, Y.; Su, P. Research of mechanics characteristics for cotton stalks compression. J. Agric. Mech. Res. 2006, 28, 171–172. [Google Scholar] [CrossRef]
  30. The State General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China; China National Standardization Administration Commission. Test Method for Flexural Properties of Sandwich Constructions; National standards of the People’s Republic of China: Beijing, China, 2005; pp. 1–9. [Google Scholar]
  31. Hu, G.; Bu, L.; Chen, J. Simulation to determination of significant parameters on apple stress for combing harvesting in trellis trained trees. Sci. Hortic. 2020, 274, 109654. [Google Scholar] [CrossRef]
  32. Bu, L.; Hu, G.; Chen, C.; Sugirbay, A.; Chen, J. Experimental and simulation analysis of optimum picking patterns for robotic apple harvesting. Sci. Hortic. 2020, 261, 108937. [Google Scholar] [CrossRef]
Figure 1. The photos and 3D models of the ripe fruit (a), half-ripe fruit (b), unripe fruit (c), flowers (d), and leaves (e).
Figure 1. The photos and 3D models of the ripe fruit (a), half-ripe fruit (b), unripe fruit (c), flowers (d), and leaves (e).
Forests 12 00699 g001
Figure 2. The meshed model (a) and the result (b) of the ripe fruit in the detachment simulation.
Figure 2. The meshed model (a) and the result (b) of the ripe fruit in the detachment simulation.
Forests 12 00699 g002
Figure 3. The meshed model (a) and the result (b) of the half-ripe fruit in the detachment simulation.
Figure 3. The meshed model (a) and the result (b) of the half-ripe fruit in the detachment simulation.
Forests 12 00699 g003
Figure 4. The meshed model (a) and the result (b) of the unripe fruit in the detachment simulation.
Figure 4. The meshed model (a) and the result (b) of the unripe fruit in the detachment simulation.
Forests 12 00699 g004
Figure 5. The meshed model (a) and the result (b) of the flower in the detachment simulation.
Figure 5. The meshed model (a) and the result (b) of the flower in the detachment simulation.
Forests 12 00699 g005
Figure 6. The meshed model (a) and the result (b) of the leaf in the detachment simulation.
Figure 6. The meshed model (a) and the result (b) of the leaf in the detachment simulation.
Forests 12 00699 g006
Figure 7. The detachment experiments results of the ripe fruit, half-ripe fruit, and unripe fruit.
Figure 7. The detachment experiments results of the ripe fruit, half-ripe fruit, and unripe fruit.
Forests 12 00699 g007
Figure 8. The detachment experiments results of the flowers and leaves.
Figure 8. The detachment experiments results of the flowers and leaves.
Forests 12 00699 g008
Table 1. The material mechanics parameters of the plant parts.
Table 1. The material mechanics parameters of the plant parts.
ItemDensity (t/mm3)Elastic Modulus (MPa)Poisson Ratio
The ripe fruit9.3339 × 10−101.3440 × 10−10.4000
The half-ripe fruit9.6742 × 10−104.0970 × 10−10.4000
The unripe fruit8.5822 × 10−104.46940.4000
The flower3.1461 × 10−111.5500 × 10−10.4500
The leaf5.3178 × 10−104.20900.4500
The calyx of ripe fruit3.3684 × 10−113.5500 × 10−20.4200
The calyx of half-ripe fruit4.3356 × 10−111.6100 × 10−20.4200
The calyx of unripe fruit1.5180 × 10−108.8500 × 10−20.4200
The calyx of flower8.3132 × 10−111.2900 × 10−20.4200
The stem of ripe fruit2.4618 × 10−98.21100.3800
The stem of half-ripe fruit2.1556 × 10−91.13220.3800
The stem of unripe fruit2.2060 × 10−94.94440.3800
The stem of flower1.5778 × 10−93.26110.3800
The branch of ripe fruit4.1148 × 10−104.2564 × 1020.3000
The branch of half-ripe fruit1.6544 × 10−91.9119 × 1020.3000
The branch of unripe fruit1.0988 × 10−92.5738 × 1020.3000
The branch of flower1.2096 × 10−95.6472 × 1020.3000
The branch of leaf9.9668 × 10−104.5531 × 1020.3000
Table 2. The detachment forces of L. barbarum.
Table 2. The detachment forces of L. barbarum.
ItemMaximum ValueMinimum ValueMean ValueStandard Deviation
The first detachment force of ripe fruit (R1)1.260.320.620.33
The second detachment force of ripe fruit (R2)2.540.571.210.58
The first detachment force of half-ripe fruit (H1)3.350.161.080.86
The second detachment force of half-ripe fruit (H2)3.620.701.510.82
The first detachment force of unripe fruit (U1)2.150.381.130.52
The second detachment force of unripe fruit (U2)3.161.322.100.55
The first detachment force of flowers (F1)1.100.490.750.20
The first detachment force of leaves (L1)3.931.472.940.85
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, J.; Ma, T.; Inagaki, T.; Chen, Q.; Gao, Z.; Sun, L.; Cai, H.; Chen, C.; Li, C.; Zhang, S.; et al. Finite Element Method Simulations and Experiments of Detachments of Lycium barbarum L. Forests 2021, 12, 699. https://doi.org/10.3390/f12060699

AMA Style

Zhao J, Ma T, Inagaki T, Chen Q, Gao Z, Sun L, Cai H, Chen C, Li C, Zhang S, et al. Finite Element Method Simulations and Experiments of Detachments of Lycium barbarum L. Forests. 2021; 12(6):699. https://doi.org/10.3390/f12060699

Chicago/Turabian Style

Zhao, Jian, Te Ma, Tetsuya Inagaki, Qingyu Chen, Zening Gao, Lijuan Sun, Haoxuan Cai, Chao Chen, Chuanlin Li, Shixia Zhang, and et al. 2021. "Finite Element Method Simulations and Experiments of Detachments of Lycium barbarum L." Forests 12, no. 6: 699. https://doi.org/10.3390/f12060699

APA Style

Zhao, J., Ma, T., Inagaki, T., Chen, Q., Gao, Z., Sun, L., Cai, H., Chen, C., Li, C., Zhang, S., Tsuchikawa, S., & Chen, J. (2021). Finite Element Method Simulations and Experiments of Detachments of Lycium barbarum L. Forests, 12(6), 699. https://doi.org/10.3390/f12060699

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