Open Access
Issue
Mechanics & Industry
Volume 22, 2021
Article Number 16
Number of page(s) 8
DOI https://doi.org/10.1051/meca/2021017
Published online 09 March 2021

© D. Kang et al., Hosted by EDP Sciences 2021

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

With the great increase of our labor cost, there is a pressing need to possess the automated picking machine of camellia fruit for helping harvest in large scale farm [1]. Because of poor working environment, the conventional mechanism is hard to meet the demands [2]. Therefore, developing the picking device with high reliability and success rate of harvesting is indispensable, however, it strongly depends on the kinematic characteristics, as well as its weight [3].

In order to ensure the high precision movement of picking mechanism, Li et al. [4] presented an image recognition algorithm based on preference artificial immune net to meet the demand of morphological features recognition in picking robot machine-vision system. To develop the recognition rates of machine-vision system in camellia fruit picking robot, the approach of multi-feature integration utilizing preference aiNet was proposed [5], which could offer the accuracy of multi-feature integration exceeding 90% in both sunny and cloudy day. Gao et al. [6] invented a picking actor in oil-tea camellia fruit picking machine of tooth comb type, and the motion equation of executive mechanism was established by the D-H matrix transform. Rao et al. [7] designed a motor-driven picking actuator of camellia fruit with rotate rubber roller, and the critical conditions for dropping the camellia fruit off and the efficient adjustable range between upper and lower rubber roller of swinging strut were both determined. In order to improve working efficiency of motor-driven picking actuator of camellia fruit with rotate rubber roller, a hydraulic-driven camellia fruit picking machined was developed [8], and found that the significant factors affecting the damage rate of camellia bud were represented as the diameter of the rubber roller, the distance between top and bottom roller and the rotational speed of the rotating frame. Although those facilities were able to conduct on the picking function, low efficiency and manual auxiliary picking limited their engineering applications. Therefore, it needs to put forward a new mechanism.

In this paper, a novel multi-links-based picking machine was designed and produced firstly. The kinematic characteristics of executive body were calculated by theoretical analysis, and the simulated acceleration at pivotal position was compared with the experimental data. Then, based on the Kriging surrogate model, the optimization design of kinematic characteristics and lightweight for this prototype was both conducted. The optimal solution was obtained through the non-dominated sorting genetic algorithm II(NSGA-II).

2 Principle of proposed picking machine

Based on the forest field investigations and comparison of current picking machines [19], this paper presents a novel multi-links-based structure, as shown in Figure 1. In the picture, 1 is connection support, 2 is hydraulic motor, 3 is chassis, 4 is crank slider mechanism, 5 is hydraulic locking mechanism, 6 is tightening rod, and 7 is excavator boom. The main working process included following steps. Firstly, the excavator boom lifted the tightening rod over camellia tree, and then fell down and hooded the target. Secondly, the hydraulic locking mechanism generated the driving force to make the tightening rod narrow down the opening range until it was locked. Finally, the hydraulic motor rotated to drive the crank slider mechanism for reciprocating motion, which would frequently slap the branches and shake the camellia fruits. All power came from the hydraulic source of the excavator. According to the design requirements, the rotating speed of crank was 450 r/min, the radius of crank was 45 mm, the length of link rod was 360 mm, the thickness of crank was 18 mm, the radius of link rod was 15 mm, the weight of load was 150 kg and the weight of crank slider mechanism was 11.39 kg, which were listed in Table 1. Some of parametes were explained in Section 3.

thumbnail Fig. 1

Structure diagram of picking machine.

Table 1

Motion parameters of executive mechanism.

3 Estimation of kinematic characteristics

The kinematic characteristics of vibration-based picking machine were important to ensure the motion state. In order to understand the law of motion, the dynamic information of the pivotal crank slider mechanism was analysed by theoretical equation. Before constructing the motion equation, the kinematic diagram of mechanism was shown in Figure 2 to illustrate motion parts more clearly. In this picture, 1 is crank slider mechanism, 2 is hydraulic locking mechanism, 3 is tightening rod, and 4 is the target tree. As is known to all, the motion characteristics of slide block in crank slider system were the most important, and the motion schematic diagram of crank slider mechanism was also shown in Figure 3. In this picture, 1 is crank, and 2 is link rod.

According to the kinematic relation and the Newton's law, the displacement of slider could be expressed as following: x = r 0 cosθ + l r 0 2 sin 2 θ / 2 l (1)

where θ is rotating angle of crank.

The velocity of slidercould be described as following: v = w r 0 ( sinθ+ r 0 sin2 θ / 2 l ) (2)

The acceleration of slider could be determined as following: a = w 2 r 0 ( cosθ + r 0 cos2 θ / l ) (3)

Then, the change characteristics of those kinematic parameters in a circle of crank were shown in Figures 46, respectively. The maximum displacement was 405 mm, and the minimum value was 315 mm, while the origin was located at the centre of crank. The maximum velocity of slide block was 2.14 m/s, as shown in Figure 5. However, the maximum acceleration of slide block during the running period reached to be 87.44 m/s2, while the acceleration at the starting and ending position was almost −110 m/s2. The huge motion value was not suitable for this picking machine, and thus the optimization design of this executive mechanism was needed.

In order to verify the accuracy of calculated results, the acceleration testing was conducted on the prototype in a camellia fruit forest, as shown in Figure 7. The dynamic data acquisition instrument was TMR-211 produced by the Tokyo measurement research institute, and the sampling frequency was 1000 Hz, as shown in Figure 8. The arrangement position of accelerometer is shown in Figure 9, and its sensitivity was 100 mV/g. The testing time was set to be 5600 s, and then the acquired data is shown in Figure 10. The acceleration value should be multiplied by a scale factor g and it is 9.8. Then, the maximum tested acceleration was 90.11 m/s2 during the in-service time, which was quite close to the calculated one. Therefore, the calculating approach mentioned above was acceptable.

thumbnail Fig. 2

Schematic diagram of mechanism.

thumbnail Fig. 3

Motion schematic diagram of crank slider mechanism.

thumbnail Fig. 4

Displacement of slider in a circle.

thumbnail Fig. 5

Velocity of slider in a circle.

thumbnail Fig. 6

Acceleration of slider in a circle.

thumbnail Fig. 7

Acceleration testing of slider in a camellia fruit forest.

thumbnail Fig. 8

Layout of acceleration sensor on the slider.

thumbnail Fig. 9

Testing instrument.

thumbnail Fig. 10

Experimental acceleration data of slider.

4 Optimization design

4.1 Kriging surrogate model

Kriging method belongs to a kind of interpolation model, and is defined as the linear weighting result of known sample function response value as following [10,11]: y ( x ) = i = 1 n w ( i ) y ( i ) (4)where w (i) is weighting factor.

In order to determine the weighting factor, the unknown function was considered as Gaussian static stochastic process, and could be described as following [10,11]: Y ( x ) = β 0 + Z ( x ) (5)where β0 is unknown constant, and stands for the mathematical expectation of Y(x).Z(x) isstatic random process with zero mean and σ 2 variance.

By taking the correlations between random variables in design space into account, the covariance could be expressed as following [10,11]: Cov [ Z ( x ) , Z ( x ) ] = σ 2 R ( x , x ) (6)where R(x, x ) is correlation function.

In order to search for the optimal weighting factor w (i), the mean square error of y ( x ) has to be minimum value, and could be determined as following [10,11]: MSE[ y ( x ) ] = E [ ( w T Y S ) - Y ( x ) ) 2 ] (7)

Based on the Lagrange multiplier approach, the optimal weighting factor could be obtained from solving the following linear system of equations [10,11]: { j = 1 n w ( j ) R ( x ( i ) , x ( j ) ) + μ 2 σ 2 = R ( x ( i ) , x ) i = 1 n w ( i ) = 1 (8)where μ is Lagrange multiplier.

Then, equation (8) could be expressed as matrix representation [10,11]: [ R F F T 0 ] [ w μ ] = [ r 1 ] (9)

Substituting the solution of equation into equation (5), the estimating value of Kriging model could be described as following [10,11]: y ( x ) = [ r ( x ) 1 ] T [ R F F T 0 ] 1 [ y S 0 ] (10)

By solving the inversion based on the block matrix, it could be expressed as following [10,11]: y ( x ) = β 0 + r T ( x ) R 1 ( y s β 0 F ) (11)where β0 is equal to ( F T R 1 F ) 1 F T R 1 y S .

4.2 Determination of optimization objectives and variables

According to the calculated and tested results, the acceleration of slide block was higher than the allowable maximum one, which was around 60 m/s2. As mentioned above, the acceleration of slide block depended on the rotating speed of crank, radius of crank and length of link rod. Moreover, the dynamic behaviour also depended on the structuraldimension, however, the dimension size was relative with the weight as well. In order to improve the kinematic characteristics and limit the weight at the same time, based on the Kriging surrogate model, the optimization objectives could be directly described as following: min y 1 = f 1 ( x 1 , x 2 , x 3 ) min y 2 = f 2 ( x 2 , x 3 , x 4 , x 5 ) (12)where y1 is the maximum acceleration of slide block, and y2 is the weight of crank slidermechanism. x1 is the rotating speed of crank, and x1 ∈ [350, 550]. x2 is the radius of crank, and x2 ∈ [30, 60]. x3 is the length of link rod x3 ∈ [260, 460].x4 is the thickness of crank, and x4 ∈ [12, 24]. x5 is the radius of link rod x5 ∈ [10, 20].

In order to construct the relationship between optimization objective and three design variables, the optimal Latin hypercube sampling method was utilized to obtain the sample data [1214]. In this paper, 30 sample points were generated, and 6 of them were listed in Table 2. According to the equation (3), the corresponding maximum acceleration of slide block could be calculated, as well as the weight of crank slider mechanism.

Then, based on the sample points, the Pareto diagram of optimization objective y1 was shown in Figure 11. From this picture, the design variable x1 had the greatest influence on the optimization objective y1 for almost 50%, and following with the design variable x2 and item x1-x2. Most of terms had a positive influence on the target, except item x 2 2 and x 3 2 . The main effect of three design variables on the target was also shown in Figure 12, which illustrates that there was positive correlation between optimization target and three design variables. The interaction effect between design variable x1 and x2 is displayed in Figure 13. It could be clearly seen that there is almost parallel with each other, which implied that the impact on optimization objective from the design variables was independent.

Similarly, the Pareto diagram of optimization objective y2was shown in Figure 14. It could be clearly seen that the design variable x2 had the greatest influence on the optimization objective y2 for almost 40%, and following with the design variable x4 and item x2 − x4. Most of terms had a positive influence on the target, except item x 4 2 . Figure 15 shows that the main effect of four design variables produced on the target, which implied that there was positive correlation between optimization target and all design variables. The interaction effect between design variable x2 and x4 was shown in Figure 16. They were not parallel with each other, which meanthat the impact on optimization objective depended on each other for these two design variables.

Table 2

Sample data of design variables and corresponding response valued.

thumbnail Fig. 11

Pareto diagram of optimization objective y1.

thumbnail Fig. 12

Main effect for all design variables.

thumbnail Fig. 13

Interaction effect between x1 and x2.

thumbnail Fig. 14

Pareto diagram of optimization objective y2.

thumbnail Fig. 15

Main effect for all design variables.

thumbnail Fig. 16

Interaction effect between x1 and x2.

4.3 Optimization results

According to the references [1214], the relationship between optimization target and design variables obtained from the Kriging surrogate model was implicit, so that the explicit expression could not be given here. Before starting the optimization, the surrogate model was testified by comparing the predicted response values with the calculated ones for the other ten sample points, shown in Figures 17 and 18. It could be clearly seen that all points were almost located at the equivalent value line, which mean that the accuracy of the surrogate model was acceptable.

Then, the non-dominated sorting genetic algorithm II was utilized in this paper to search for the optimal solution [1214]. The number of population was taken as 50, and the number of iteration was regarded as 100. After 38 min, the optimal solution was acquired with 5000 iterations, listed in Table 3. All design variables were reduced, and the optimized maximum acceleration of slide block was 60.07 m/s2, which could totally meet the requirement. Compared with the tested maximum acceleration of slide block 90.11 m/s2, it was decreased by 31.30%. The weight of crank slider mechanism was reduced by 27.51% as well. It was changed from 11.39 to 8.26 kg.

thumbnail Fig. 17

Verification of Kriging surrogate model y1.

thumbnail Fig. 18

Verification of Kriging surrogate model y2.

Table 3

Optimized results.

5 Conclusion

In this paper, a novel multi-links-based picking machine used in camellia fruit forest was presented. The working principle and procedure of picking machine was explained. Then, based on the theoretical analysis, the motion equation of crank slider mechanism was constructed, and the kinematic characteristics were estimated by the calculated maximum displacement, velocity and acceleration of slide block. The calculated acceleration agreed well with the experimental data, but it was higher than the design requirement. In order to improve the kinematic characteristics and reduce weight of crank slider mechanism at the same time, the optimization design was conducted by utilizing the Kriging surrogate model. The relationship between optimization objectives and design variables was constructed and analysed. Through the non-dominated sorting genetic algorithm II, the optimal solution was found. The optimized maximum acceleration of slide block decreased by 31.30% with reduction of three design variables in different extent, and the weight of crank slider mechanism was reduced by 27.51% as well.

Nomenclature

r0 : Radius of crank (mm)

w : Rotating speed of crank (r/min)

l : Length of link rod (mm)

t : Thickness of crank (mm)

θ : Rotating angle of crank (°C)

s : Displacement of slide block (mm)

v : Velocity of slide block (m/s)

a : Acceleration of slide block (m/s2)

Acknowledgments

The authors gratefully thank for the support of Science and Technology Innovation Project of Hunan Province (Grant No.: 2018NK1030 and 2018NK2065) and Forestry Science and Technology Promotion Demonstration Project of Central Finance (Grant No.: 2018XT02) and Forestry Science and Technology Innovation Project of Hunan Province (Grant No.: XLK201910).

References

  1. Z.H. Yang, Y.Q. Dong, Problems and solutions of camellia industrial development in China, Journal of Jiangxi Agricultural University (Social Sciences Edition) 4 , 40–42 (2010) [Google Scholar]
  2. R. Polat, I. Gezer, M. Guner, Mechanical harvesting of pistachio nuts, Journal of Food Engineering 79 , 1131–1135 (2007) [Google Scholar]
  3. G.L. Blanco-Roldan, J.A. Gil-Ribes, K. Kouraba, Effects of trunk shaker duration and repetitions on removal efficiency for the harvesting of oil olives, Applied Engineering in Agriculture 25 , 329–334 (2009) [Google Scholar]
  4. L.J. Li, X. Li, Z.C. Gao, Camellia fruit image recognition based on preference artificial immune net, Transactions of the Chinese Society of Agricultural Machinery 43 , 209–213 (2012) [Google Scholar]
  5. X. Li, L.J. Li, Z.C. Gao, Image recognition of camellia fruit based on preference for aiNET multi-features integration, Transactions of the Chinese Society of Agricultural Engineering 28 , 133–137 (2012) [Google Scholar]
  6. Z.C. Gao, L.J. Li, X. Li, Development and test of picking actor in oil-tea camellia fruit picking machine of tooth comb type, Transactions of the Chinese Society of Agricultural Engineering 29 , 19–25 (2013) [Google Scholar]
  7. H.H. Rao, L.Y. Zhang, D.S. Huang, Design and test of motor-driven picking actuator of camellia fruit with rotate rubber roller, Transactions of the Chinese Society of Agricultural Machinery 49 , 115–121 (2018) [Google Scholar]
  8. H.H. Rao, D.S. Huang, Y.L. Wang, Design and experiment of hydraulic-driven camellia fruit picking machine, Transactions of the Chinese Society of Agricultural Machinery 50 , 133–139+147 (2019) [Google Scholar]
  9. S. Ye, L.J. Li, Z.C. Gao, Movement simulation and analysis of vibratory camellia fruit-harvesting machine based on Pro/E, Journal of Central South University of Forestry & Technology 33, 166–169 (2013) [Google Scholar]
  10. Z.H. Han, Kriging surrogate model and its application to design optimization: a review of recent progress, Acta Aeronautica et Astronautica Sinica 37 , 3197–3225 (2016) [Google Scholar]
  11. C.J. Mi, Z.Q. Gu, Y. Zhang, Frame weight and anti-fatigue co-optimization of a mining dump truck based on Kriging approximation model, Engineering Failure Analysis 66 , 99–109 (2016) [Google Scholar]
  12. Y.F. Tong, Y. He, Z.B. Gong, Research on genetic algorithm-based rapid design optimization, Mechanika 18 , 569–573 (2012) [Google Scholar]
  13. A. Povilionis, A. Bargelis, Structural optimization in product design process, Mechanika 81 , 66–70 (2010) [Google Scholar]
  14. C.J. Mi, W.T. Li, W.G. Wu, Fuzzy fatigue reliability analysis and optimization of A-type frame of electric wheel dump truck based on response surface method, Mechanika 25 , 44–51 (2019) [Google Scholar]

Cite this article as: D. Kang, Z.J. Chen, Y.H. Fan, C. Li, C. Mi, Y.H. Tang, Optimization on kinematic characteristics and lightweight of a camellia fruit picking machine based on the Kriging surrogate model, Mechanics & Industry 22, 16 (2021)

All Tables

Table 1

Motion parameters of executive mechanism.

Table 2

Sample data of design variables and corresponding response valued.

Table 3

Optimized results.

All Figures

thumbnail Fig. 1

Structure diagram of picking machine.

In the text
thumbnail Fig. 2

Schematic diagram of mechanism.

In the text
thumbnail Fig. 3

Motion schematic diagram of crank slider mechanism.

In the text
thumbnail Fig. 4

Displacement of slider in a circle.

In the text
thumbnail Fig. 5

Velocity of slider in a circle.

In the text
thumbnail Fig. 6

Acceleration of slider in a circle.

In the text
thumbnail Fig. 7

Acceleration testing of slider in a camellia fruit forest.

In the text
thumbnail Fig. 8

Layout of acceleration sensor on the slider.

In the text
thumbnail Fig. 9

Testing instrument.

In the text
thumbnail Fig. 10

Experimental acceleration data of slider.

In the text
thumbnail Fig. 11

Pareto diagram of optimization objective y1.

In the text
thumbnail Fig. 12

Main effect for all design variables.

In the text
thumbnail Fig. 13

Interaction effect between x1 and x2.

In the text
thumbnail Fig. 14

Pareto diagram of optimization objective y2.

In the text
thumbnail Fig. 15

Main effect for all design variables.

In the text
thumbnail Fig. 16

Interaction effect between x1 and x2.

In the text
thumbnail Fig. 17

Verification of Kriging surrogate model y1.

In the text
thumbnail Fig. 18

Verification of Kriging surrogate model y2.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.