InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Agricultural and Biological Sciences » "Advances in Seed Biology", book edited by Jose C. Jimenez-Lopez, ISBN 978-953-51-3622-4, Print ISBN 978-953-51-3621-7, Published: December 6, 2017 under CC BY 3.0 license. © The Author(s).

Chapter 16

Application of Mechanics to Plant Seeds as a Granular or Particulate Material

By Michal Petrů and Ivan Mašín
DOI: 10.5772/intechopen.70838

Article top


Plant seeds as granular or particulate material.
Figure 1. Plant seeds as granular or particulate material.
Different shapes of seeds.
Figure 2. Different shapes of seeds.
Basic dimensions of seed with oval shape (L—length, W—Width, and T—Thickness).
Figure 3. Basic dimensions of seed with oval shape (L—length, W—Width, and T—Thickness).
Seed deformation at different forces.
Figure 4. Seed deformation at different forces.
Selected types of deformation diagrams for different plant seeds.
Figure 5. Selected types of deformation diagrams for different plant seeds.
Perzyna rheology model for single Jatropha curcas L. seed compression [21].
Figure 6. Perzyna rheology model for single Jatropha curcas L. seed compression [21].
Jatropha curcas L. seed (CAD model, FEM model, and simulation model).
Figure 7. Jatropha curcas L. seed (CAD model, FEM model, and simulation model).
Deformation of Jatropha curcas L. seed (real sample, FEM model, vector plot of FEM model).
Figure 8. Deformation of Jatropha curcas L. seed (real sample, FEM model, vector plot of FEM model).
Time response of Jatropha curcas L. seed deformation.
Figure 9. Time response of Jatropha curcas L. seed deformation.
Comparison of theory and FEM model of Jatropha curcas L. seed for compressive force.
Figure 10. Comparison of theory and FEM model of Jatropha curcas L. seed for compressive force.
Strain response: compressive stress, compressive energy, compressibility (individual seed).
Figure 11. Strain response: compressive stress, compressive energy, compressibility (individual seed).
Contact area aS of the seeds on their outer boundary.
Figure 12. Contact area aS of the seeds on their outer boundary.
Symmetry planes for modeling (left) and axial section for 2D model (right).
Figure 13. Symmetry planes for modeling (left) and axial section for 2D model (right).
Seed-system deformation in time.
Figure 14. Seed-system deformation in time.
Growth of the cracks in seeds during pressing (FEM model—above, experiment—bottom).
Figure 15. Growth of the cracks in seeds during pressing (FEM model—above, experiment—bottom).
Comparison of experiment and FEM model for compressive force.
Figure 16. Comparison of experiment and FEM model for compressive force.
Strain response: compressive stress, compressive energy, compressibility (system of seeds).
Figure 17. Strain response: compressive stress, compressive energy, compressibility (system of seeds).

Application of Mechanics to Plant Seeds as a Granular or Particulate Material

Michal Petrů and Ivan Mašín
Show details


The mechanical behavior of plant seeds as a granular or particulate material dramatically differs from the mechanical behavior of solid materials. This difference is caused by the possibility of partially autonomous movement and rotation of seeds, their mutual contacts, or due to the occurrence of the second fluid phase among the seeds at the stage of their moving or processing. For obtaining the economic effects from the seeds (energy, nutrients, livestock, etc.), seeds must often be subjected to mechanical treatment. In this context application of mechanics as science concerned with the behavior of physical bodies when subjected to forces or displacements is very important and needed. One of the goals of this chapter is therefore to provide an overview for readers who are not primarily concerned with mechanics but who are interested in the behavior of seeds in the context of biology, agriculture, and pharmacy or food industry. This chapter is therefore focused on both an overview of the principles of mechanics of granular or particulate materials and the presentation of experimental results particularly in the area of mechanical extraction of oil from seeds.

Keywords: seed, mechanics, granular material, particulate material, rheological model, FEM model, oil extraction

1. Introduction

Human life has become dependent on plants for the qualities and developments that they provide, which include agriculture, food production, and chemical industry. Plant seeds are one of the most important agricultural materials which affects billions of people. Seeds are the result of sexual reproduction in plants. Seeds are of immense biological and economic importance. They contain protein, carbohydrate, starch, and oil reserves that help in the early stages of growth and development in a plant. These reserves are what make many plant seeds important source for a large proportion of the world’s inhabitants. For obtaining the economic effects from the seeds, they must often be subjected to mechanical treatment (manipulation, conveying, separation, hulling, pressing, purification, packaging, etc.). That is why studies of mechanical behavior of seeds are nowadays an important field of science and engineering. The research of how seeds as granular or particulate material are deformed, cracked, or how they flow has big importance for industries as biotechnology, pharmacy, agriculture, or food processing and various nanotechnologies. Results of that research are very important for both proper design of machines or equipment and for technological processes optimization. Current research has shown that seeds mechanical behavior is changing over time within the context of their moisture or oil content, acting forces, geometrical parameters, process of crack formation, porosity, visco-elastic properties, or in connection with the changes in the mutual arrangement or reorganization of the seed layers. The resulting nonlinear mechanical behavior and complex motion or spatial orientation of individual seeds or seed layers are very difficult to study, analyze, describe, and predict. For these scientific and engineering tasks, it is necessary to combine knowledge of basic physical seed characteristics, trivial calculations, modeling with the help of simplified or rheological models, or modeling using numerical methods such as FEM (finite element method) and DEM (discrete element method) and experimentation combined with relevant measurements. Also, special methods utilizing the principle of FEM with variable geometry like FDMs (fictitious domain methods) and IB-BCE (immersed boundary—body conformal enrichment) can be applied. Almost all the physical processes can be solved using numerical methods. The difficulty lies in the time-consuming calculation of highly complex and nonlinear problems. In this case, the calculation time can only be reduced by using supercomputers. The following chapter is therefore divided into several sections dealing with descriptive seed properties, seed mechanics principles, rheological models, and advanced modeling of seeds behavior.

2. Geometrical and mechanical properties of seeds

Seeds can generally be seen as granular or particular material (Figure 1). The actual state of the seeds as a material depends on the time from their harvest, moisture/oil content, and on mechanical, thermal, chemical, or other effects induced by the other subjects or by environment. Granular material is composed of small, discrete entities as opposed to being continuous. The granules can be characterized as a mathematical set of macroscopic particles defined by their shape. Particulate material is substance composed of mutually contacting solid particles, or structural units, within the liquid and/or the gaseous phase [1].For various reasons, physical properties of seeds need to be studied. A physical property is any property of matter or energy that can be measured. It is an attribute of matter that can be observed or perceived. The physical properties of seeds can be divided into following categories: geometrical, mechanical, chemical, electrical, optical, and others. For mechanical engineering tasks, geometrical and mechanical properties are most important. Geometrical properties are those that can be derived from the geometry of a seed body. Seeds, however, come in a great variety of shapes and sizes. Seed size varies among species over a range of 10 orders of magnitude with extremes represented by orchids, e.g., Goodyera repens (weight 0.000002 g) and by the double coconut palm Lodoicea maldivica (weight 18,000–27,000 g) [2]. Seeds may be round, egg-shaped, triangular, long and slender, curved, coiled like a snail or irregular in shape (Figure 2). Some have a groove or depression with a ridge along the length; others may be flattened at one or both ends. Moreover, the seed shape may vary in immature or poorly developed seed. Since seed shapes are highly variable, their shapes are simplified (very often on an oval shape), and standardized geometrical parameters are used for the purpose of solving tasks from basic or advanced mechanics and simulation. This simplification is very important as a mean by which the size and shape of an irregular shaped seed can be easily described and quantified whether the seed is treated as an individual unit or as one that is a representative of many seeds in layer or bulk. Selected and often used geometrical properties are illustrated in Figure 3 and described in Table 1. Size of seeds determines the efficiency of processing and storage and the quality of semi-product or final product. For this reason, it is also necessary during mechanical problems solving to take into account the stochastic nature of geometrical properties, for example, distinguish size fractions of one type plant seed (most commonly divided into groups of small, medium and large seeds), frequency distribution curves, or coefficient of variation of individual geometrical property. Mechanical properties are physical properties that a material exhibits upon the application of forces. Examples of mechanical properties are the modulus of elasticity, tensile strength, elongation, hardness, and fatigue limit. Mechanical properties occur as a result of the physical properties inherent to each material and are determined through a series of standardized mechanical tests. Some of important mechanical properties of seeds are listed in Table 2.


Figure 1.

Plant seeds as granular or particulate material.


Figure 2.

Different shapes of seeds.


Figure 3.

Basic dimensions of seed with oval shape (L—length, W—Width, and T—Thickness).

Geometrical property of seedNomenclatureUnitFormula
Length of one seedLmm
Width of seedWmm
Thickness of one seedTmm
Arithmetic mean diameterDammDa=L+W+T/3(1)
Geometric mean diameterDgmmDg=LBT1/3(2)
Equivalent diameterϕmm4Sm/π1/2Smmeasured area(3)
Seed surface areaSmm2S=πDg2(4)
SphericitySϕe.g., Sϕ=Dg/L(5)
Volume of one seedVmm3e.g., V=4/3πLBT(6)

Table 1.

Selected geometrical properties of seeds.

Mechanical property of seedNomenclatureUnit
Bulk modulusKMPa
Coefficient of frictionμ
Compression energyWCJ mm−3
Consumed energy at rupture pointWRJ
Deformation coefficient of mechanical behaviorBDmm−1
Elasticity modulusEMPa
Energy used for ruptureERMPa
Rupture forceRFN
Seed oil dynamic viscosityηPa s
Shear modulusGMPa
Stress in the structureσSMPa

Table 2.

Selected mechanical properties of seeds.

Research on geometrical and mechanical properties has been focused on various types of seeds, such as rapeseed (Brassica napus L.) [2], Jatropha (Jatropha curcas L.) [4, 5], sunflower (Helianthus annuus L.) [6, 7], bean (e.g., Phaseolus vulgaris) [8], arigo seeds (Dacryodes edulis) [9], pine (Pinus pinea) [10], chia seeds (Salvia hispanical.) [11], wild legumes (e.g., Canavalia cathartica) [12], melon seeds (Cucumis melo L.) [13], tung seed (Aleurites Fordii) [14], sugarbeet seed (Beta vulgaris L.) [15], sesame (Sesamum indicum L.) [16], vetch seeds (Vicia sativa L.) [17], green soybean (Glycine max) [18], quinoa (Chenopodium album L.) [19], hemp (Cannabis sativa L.) [20], and so on. Geometrical and mechanical properties of selected seeds types are listed in Table 3.

SeedL (mm)W (mm)T (mm)Dg (mm)ER (Nmm)RF (N)μ (steel) (−)
Brassica napus (fraction: medium) [3]2.12 ± 0.1081.91 ± 0.09313.53 ± 2.649
Jatropha curcas L. (var. Kanlueang, nut) [4]21.02 ± 1.0311.97 ± 0.309.58 ± 0.2813.40 ± 0.3124.44 ± 19.95146.63 ± 14.82
Helianthus annuus L. (var. Morden, moisture 6.2%) [6]9.27 ± 0.684.78 ± 0.343.32 ± 0.275.39 ± 0.416)0.45
Phaseolus vulgaris (var. Hinis) [8]11.76 ± 0.778.85 ± 0.507.66 ± 0.589.26 ± 0.5363.79 ± 23.25 (X ax)145.88 ± 33.54 (X ax)0.227 ± 0.0038
Dacryodes edulis (moisture 10.3%) [9]19.00 ± 1.112.20 ± 0.8)10.10 ± 0.8)13.2 ± 1.4metal 0.25 ± 0.009

Table 3.

Selected geometrical and mechanical properties of seeds.

3. Mechanical behavior of single seed

When analyzing the mechanical behavior of individual seed, we can identify states similar to the mechanical behavior of the elastic body. Therefore, a model of an idealized flexible body can be used to study and describe the behavior of individual seeds. In this case, the forces and moments of the forces acting on the seed cause its state of strain accompanied by its deformation and fracture (Figures 4 and 5)


Figure 4.

Seed deformation at different forces.


Figure 5.

Selected types of deformation diagrams for different plant seeds.

Initially consider the compression of seed as a compression of linear elastic substance and constant side pressure coefficient K0|=const.. Then, according to Hook’s law:


where εSE,σSE are axial deformation and stress in the direction of compression during elastic deformation, μ is Poisson’s ratio, and ES is initial stiffness modulus.

During an elastic deformation, the structure of the substance does not change, and therefore, initial stiffness modulus will be constant ES|=const.. Deformation of single seed is a function of stress σSE which characterizes the intensity of structural changes. Then, Eq. (7) can be overwritten as:


where w0 is the initial seed width.

The stress in the structure of the sample σS (MPa) can be resolved into the sum of stresses consisting of the stress elastic component σSE and stress in the plastic component σSP, and which is characterized by oil viscosity component σSη. This is responsible for the damping of the structure. Friction component σSμ of the seed sample causes significant permanent deformation with geometrical and physical damage from certain value of strain. By using numerical analysis of the isothermal compression T=const. for the seed sample, it has been observed that the elastic component of stress a function which is defined approximately by the Hooke?s law [21]. This strain is characterized by a low oiliness point/low point of structure, where the isotropic properties of the structure can be considered. Values of low oiliness point for example tested samples of J. curcas L. are shown (Table 4). This area is characterized by a constant volume strain γ|εS ≤ 40 %  strain, which can be defined by Eq. (9), and in principle, it is same as Eq. (8).


where V is the volume (in an initial position: compression time t = 0), C is the compression pressure, and ψ the compressibility [22].

Maturity stageT (mm)S (mm2)εSLow (−)
Ripe10.0 ± 1.8180.1 ± 19.30.5
Over-ripe10.2 ± 1.4181.8 ± 19.50.4
Unripe9.2 ± 1.9172.2 ± 18.30.4

Table 4.

Dimensions of the seeds prior deformation tests and the coefficients characterizing the deformation of seeds (data in the table are mean ± SD).

S is the projected area of the seed, and ɛSlow is deformation of the seed (low oiliness point).

This can be visualized as a reversible elastic compression of the Hooke member, which is defined by the initial elastic modulus E|εS ≤ low point of structure and which can be supplemented by the damping factor of the Newtonian viscous member η|εS ≤ low point of structure. This behavior could be partly described by a generalized Maxwell rheological model (or generalized Maxwell model). The higher compression of the test structure from the point of oiliness or oil point passes the functional dependence of elastic stress in the elastic or plastic or visco-plastic components σSE=fεSlowpointofstructureσSP=fεSlowpointofstructure. This is reflected in breach of the yield stress of the test sample structure (Figure 5). Here, a comprehensive change in the volumetric strain γ|≥low point of structure was reflected and spread with the speed of plastic strain εSP which caused increased friction μ = μS and creep structure due to extrusion of the oil component from the structure. Functional dependency was reflected in a significant increase in compression strength with a parabolic character which is different to the elastic-plastic deformation of elastic structures. There is a whole range of analytical and empirical relations in the literature that extend the seed deformation description by porosity. Define that the density of the undried seed is ρS = mS/VS and the density of the dried seed at 105°C is ρD = mS/VD. The pore volume can be then defined by the porosity number:


where eS = ρS/ρD is the ratio by which the compression can be defined according to (9).

In principle, it is a semi-logarithmic Walker-Balsin’s law applicable to, for example, powder metallurgy:


where e0 is the initial porosity number at σS = σ0.

Let us assume a single seed is comprised under pressure σS1 to width w1 and subsequently comprised under pressure σS2 (σS1 < σS2) to width w2 < w1. Then, we can introduce Eq. (12), which (in the case σ0 is unit pressure) can be generalized to Eq. (13).


For FEM simulation of such anisotropic behavior, we can use equilibrium equations or applications of the constitutive equations of a rheology model with modified parameters like Maxwell model, Kelvin model, Saint-Venant model or Perzyna model. Acceptance of certain seed states as a particulate material leads to the use of rheology in describing its mechanical behavior. The science of rheology is important to the study of the flow behavior of solids suspended in fluids. Rheological models (visco-elastic models) used to characterize seeds are classified as non-Newtonian fluids. For these fluids, no constant of proportionality exists between shear stress and shear rate; their viscosity varies with changing shear rate. For seeds, it is advisable to use the Perzyna rheological model in particular (Figure 6).


Figure 6.

Perzyna rheology model for single Jatropha curcas L. seed compression [21].

The Perzyna model with visco-plastic parameters includes the strain rate, and its rheological behavior is applicable not only for the description of elastic–plastic material damage but also for the behavior of anisotropic structures of compressed seed samples of different crops (Figure 6) in the plastic strain state εSPlowpoint of structure. The Perzyna model derivation is based on studies of the stress limit description, which was defined by Bingham in 1920 to describe a plastic deformation of elastic materials. It is possible to use this model for quasi-static strain behind yield stress. Also, the model defines the time-dependent response of the structure cohesion, for example, the determination of oil extrusion initiation (so-called low oiliness point), immediate friction, and viscous and elastic parameters as functions or constants. Also, the Perzyna model can be used for modeling of the flow of soil such as clay. The velocity of stress σ̇S propagation during plastic deformation can be expressed by the constitutive Eq. (14), where the progress of plastic deformations is given by changing of volumetric strain γ, which can assume various values:


where σ̇S is stress rate of structure, ES initial elastic modulus, ε̇S is strain rate of structure, and ε̇SP is visco-plastic strain rate.

The plastic deformation of material structure denoted by the stress σS can be resolved according to Eq. (15) to the sum of the initial yield stress σ¯Sσ¯σSP and product of equivalent visco-plastic strain γSP and function of the strain hSP that describes following condition: if hSP<0, then a softening of the structure occurs hSP>0, and then a hardening occurs:


The numerical model can be subsequently designed for isothermal or thermal compression of different seed types. An example is the numerical isothermal study of three samples of Jatropha seeds (mature, immature, and precisely mature) with the same initial geometry ΩSripe=ΩSunripe=ΩSoverripet=0T=const.. Differently mature seeds show different deformation of the original geometry during compression t ≠ t0, which has been described in detail in [21]. The degree of compression of magnitude δ is also the compressibility function ψ according to Eq. (9). Seed dimensions for the CAD/FEM model for differently mature seeds are based on data listed in the Table 4. Modeled seeds had almost ellipsoidal shape with dimensions ~17 × 10 ± 1.8 × 9.3 mm, as shown in Figure 7.


Figure 7.

Jatropha curcas L. seed (CAD model, FEM model, and simulation model).

Compression experiments and simulations of J. curcas L. seeds were performed up to plastic deformation (Figure 8). The model showed the highest concentration of the main stress on the seed circumference for each level of maturity of tested seeds. This is because the greatest stress is concentrated in the center of the seed and consequently causes the extrusion of the structure from the center toward the periphery. This effect changed and damaged the smallest radius of the seed shape as seen from Figures 8 and 9. From this, it can be deduced that the seed is damaged by the internal pressure resulting from external compression forces.


Figure 8.

Deformation of Jatropha curcas L. seed (real sample, FEM model, vector plot of FEM model).


Figure 9.

Time response of Jatropha curcas L. seed deformation.

The energetic behavior of individual seeds during compression is shown in the following figures. Figure 10 shows the comparison of compression forces derived from mathematical equations and forces obtained by numerical models. It can be seen from Figure 10 that mathematical relationships (7)(13) do not allow to sufficiently describe the “point of fracture” as they are generally based on the empirical analysis and using graphs of standard mathematical functions (such as a parabola). Point of fracture can be better and more accurately determined from numerical models. Numerical and FEM models can be used to solve a whole range of engineering tasks. For example, using a FEM model, it is possible to determine the compressibility of mature, immature, and overmatured seeds (Figure 11) or to effectively assess the energy intensity of obtaining the oil from the seeds. The results of the experiment also showed differences in the values of the stress of individual seeds, which depend on their maturity. These findings can be used to find a proper seed position to control orientation of the seed in the press or to optimize the geometry of the press.


Figure 10.

Comparison of theory and FEM model of Jatropha curcas L. seed for compressive force.


Figure 11.

Strain response: compressive stress, compressive energy, compressibility (individual seed).

4. Mechanical behavior of seeds as a granular or particulate materials

Seeds as a granular or particulate material consists of individual seeds and adjacent voids filled with gas (air) or/and liquid (e.g., oil). Understanding the mechanical behavior of such kind of material, we must therefore take into account the geometrical presentation that describes spatial distribution of seeds and voids, orientation of seeds, and contacts of seeds, etc. Systems composed of seeds consist of mutually contacting phases, which can be solid (geometry, solid structure of seed, etc.), liquid phases (oil component of seed), and gas phases (gaseous environment). The liquid and gas phases fill the pores of the solid skeleton of the seed. In addition to these basic phases, there are mutual bonds between the seeds which, in terms of mechanical properties of the compressed seeds, are mainly frictional. When compressing seeds, we have to study frictional bonds seed-seed, seed-seeds, seeds-container wall, etc. To understand this complex process, we can start on the general theory of friction of solids. The two ideal solid particulate seeds, when moved relative to one another, can either slide or roll in an elastic state. Seeds in a plastic state can penetrate and break into each other. In solving mechanical problems, generally we have to distinguish dry seed surfaces, hydrodynamically lubricated surfaces, and limiting friction. The friction or contact bond is given by the normal stress σn on the geometric contact area ag (Figure 12). The normal stress of the contact points σnc on the surface area aS is then given by Eq. (16). Similarly, sliding friction τn can be described by Eq. (17).where τnc is the shear strength of contact point (assumed to be the same everywhere). From Eqs. (16) and (17) is Eq. (18):


Figure 12.

Contact area aS of the seeds on their outer boundary.


According to Eq. (18), we obtain the friction response between the seed surfaces. It follows that shear friction in the case of seeds also does not depend on the size of the surface and is directly proportional to the normal stress. The area of contact point plasticity allows a more accurate description of the plastic state of the stresses in it as well as for the emerging strain. The surface area ak is inversely proportional to the stress σnc according to Eq. (16) and at the same time increases the shear stress of the contact point [23]. This effect is particularly evident in elastic materials such as steel, in which case their friction coefficient ν increases considerably. In the case of structures such as seeds, this problem is limited by their skeletal fragility. In general, according to the friction adhesion theory, the intergranular friction coefficient is equal to the shear and normal strength of the contact point. Through this theory, we gain a function dependent on the amount of stress concentration at the boundaries between the seeds.

Therefore, the description of the mechanical behavior of the seeds as granular or particulate material must, for various reasons, be based on:

  • initial seed configuration (seeds geometry, seeds orientation, shape regularity, and ripeness);

  • the volume ratio of the individual phases in the press VSi|VA ≠ Vsi ≠ Voi (VA the volume of air, Vsi the volume of individual seeds, Voi the volume of the oil fraction of each seed);

  • the ratio of the initial elastic modules ESi|Eripe ≠ Eunripe ≠ Eoverripe (Eripe the ripe seed module, Eunripe the unripe seed module, Eoverripethe overripe seed module);

  • initial volumetric modulesGSi|Gripe ≠ Gunripe ≠ Goverripe; and

  • initial shear modulesKSi|Kripe ≠ Kunripe ≠ Koverripe.

As an example, we can consider a seed-system that will be compressed in an experimental container. Such a system will be loaded with axial stress σa and all-sided pressure σr. We will not consider a homogeneous system where the stress σa , σr would be totally consistent with the continuum, but we will consider a particle system where each part of the filled container (seeds, air, walls, and piston) contributes to σa , σr. Terzaghi [24] introduced in 1923 the construct of effective stress, which is a function of the total stress and the particle stress. The study states that the change in chamber pressure does not affect the dependency of (σa − σr) on axial strain of the system εa. Terzaghi also introduced the boundary condition that we can simplify the particular systems in the individual planes of symmetry (Figure 13). The greater the seed filling in the container, the more accurate the calculation will be, not only as a continuous system but also as a particle system. This claim can be demonstrated by the example of compressing a system of regularly arranged ellipsoidal shaped seeds (such as J. curcas L.). If we make a simple axial section of the seeds, we obtain a 2D arrangement of identical circles with a radius r, where adjacent circles touch each other (Figure 13).


Figure 13.

Symmetry planes for modeling (left) and axial section for 2D model (right).

During compression, the radius Δr will change. Radius Δr will be a function of contact forces that can be derived from the Hertz or Mindlin relationship. Basically, if the contact forces are known, deformation of the circles for axial strain and lateral strain can be calculated according to Eqs. (19) and (20).


It can be seen from Eqs. (19) and (20) that the size of the strain depends on the ratio of the main stresses σ1/σ3 ≡ σa/σr. There is a special case, where the contact forces of the ideal circular arrangement are statically determined (in case of hydrostatic tension, respectively, when σa = σr). These mathematical derivations of stress and deformation give us insights about mechanical changes of the compressed seed system. The Terzaghi’s principle is well described in [23], where a symmetrical numerical model was created to study the strain of the seed system during compression.

Figure 14 shows the transformation of seeds stored in a regular hexagonal arrangement where the initial seed filling αS will gradually increase during compression as it αS < αS1 < αSn|t0 < t1 < tn. Figure 14 also shows the reorganization of the seeds and their subsequent clustering and joining in the case of a plastic failure. Figure 15 is a comparison of the modeled individual seeds crack propagation with a real experiment. This gives us insight into the compressive behavior of individual seeds and thus the knowledge about the efficiency of the pressing process. By increasing the crack in the individual seeds, we also provide information on the movement and direction of the extruded oil.


Figure 14.

Seed-system deformation in time.


Figure 15.

Growth of the cracks in seeds during pressing (FEM model—above, experiment—bottom).

Consequently, the energy behavior of the seed system during compression can be evaluated, as shown in Figures 16 and 17. Figures compare numerical models and the experimental data (compressive forces and compressibility) for the different cylinder diameters. The results show a very good match that can be attributed to the simplified 2D symmetric model, which approximates empirically established relationships (19) and (20). A 2D seed-system modeling not only simplifies calculation and reduces time for computations but also allows for more accurate results compared to 3D modeling of single seed behavior (Figure 10). Appropriate simplification of mathematical models describing the mechanical behavior of seeds-system (edge conditions, symmetry, etc.) can provide us with results that are sufficient to comprehend complex processes and to optimize seed treatment technologies, compared to sophisticated models with disproportionate time-consuming calculation and interpretation. Mathematical models of seeds mechanical behavior are very important for different agricultural and engineering tasks (e.g., tillage mode, processing technology design, treatment control and optimization, procedures effectiveness increasing, seeds behavior prediction, etc.).


Figure 16.

Comparison of experiment and FEM model for compressive force.


Figure 17.

Strain response: compressive stress, compressive energy, compressibility (system of seeds).

As this topic goes beyond the scope of this publication, we refer readers to various works focused on various aspects of mechanic behavior of seeds, using different theories, models, and modeling methods:

  • modeling of seeds compression [21, 22, 23, 25, 31];

  • modeling of stress or force relaxation of seeds [30];

  • modeling of seeds movement and orientation [26];

  • various contact and friction models utilization [1, 29];

  • FEM application for seeds behavior modeling [21, 23]; and

  • DEM application for seeds behavior modeling [27, 28].

5. Conclusion

A detailed understanding of the mechanical behavior of plant seeds is currently an important science and technology challenge. Post-harvest geometrical and mechanical properties of the seeds or seeds mechanical behavior are important from the point of view of treatment process optimization and machines design. Basic dimensions and geometrical properties (length, width, thickness, mean diameters, seed surface area, sphericity, etc.) were described in this overview type of paper. Selected mechanical properties (e.g., compression energy, stress in the structure, coefficient of friction, rupture force, etc.) were also described. Mechanical behavior of seeds has been described and discussed in terms of both single seed and seeds-system. Selected experimental results from the research of mechanical behavior of J. curcas L. seeds were used to bring seeds behavior as a granular and particulate material closer to readers.


1 - Feda J. Mechanics of Particulate Materials: The Principles. Amsterdam, Oxford, New York: Elsevier Scientific Publishing Company; 1982
2 - Harper JL, Lovell PH, Moore K. The shapes and sizes of seeds. Annual Review of Ecology, Evolution, and Systematics. 1970;1:327-356
3 - Tanska M, Konopka I, Rotkiewicz D. Relationship of rapeseed strength properties to seed size, colour and coat fibre composition. Journal of the Science of Food and Agriculture. 2008;88:2186-2193
4 - Sirisomboon P, Kitchaiya P, Pholpho T, Mahuttanyavanitch W. Physical and mechanical properties of Jatropha curcas L. fruits, nuts and kernels. Biosystems Engineering. 2007;97:201-207
5 - Karaj S, Müller J. Determination of physical, mechanical and chemical properties of seeds and kernels of Jatropha curcas L. Industrial Crops and Products. 2010;32:129-138
6 - Gupta RK, Das SK. Physical properties of sunflower seeds. Journal of Agricultural Engineering Research. 1997;66:1-8
7 - Khodabakhshian R, Emadi B, Abbaspour Fard MH. Some engineering properties of sunflower seed and its kernel. Journal of Agriculture, Science and Technology. 2010;4:37-46
8 - Ozturk I, Kara M, Yildiz C, Ercisli S. Physico-mechanical seed properties of the common Turkish bean (Phaseolus vulgaris) cultivars ‘Hinis’ and ‘Ispir’. New Zealand Journal of Crop and Horticultural Science. 2009;37:41-50
9 - Davies RM. Some physical properties of arigo seeds. International Agrophysics. 2010;24:89-92
10 - Özgüven F, Vursavus K. Some physical, mechanical and aerodynamic properties of pine (Pinus pinea) nuts. Journal of Food Engineering. 2005;68:191-196
11 - Ixtainaa VY, Nolascoa SM, Tomás MC. Physical properties of chia (Salvia hispanica L.) seeds. Industrial Crops and Products. 2008;28:286-293
12 - Niveditha VR, Sridhar KR, Balasubramanian D. Physical and mechanical properties of seeds and kernels of Canavalia of coastal sand dunes. International Food Research Journal. 2013;20:1547-1554
13 - Mansouri A, Mirzabe AH, Ráufi A. Physical properties and mathematical modeling of melon (Cucumis melo L.) seeds and kernels. Journal of the Saudi Society of Agricultural Sciences. 2017;16:218-226
14 - Sharma V, Das L, Pradhan RC, Naik SN, Bhatnagar N, Kureel RS. Physical properties of tung seed: An industrial oil yielding crop. Industrial Crops and Products. 2011;33:440-444
15 - Dursun I, Tugrul KM, Dursun E. Some physical properties of sugarbeet seed. Journal of Stored Products Research. 2007;43:149-155
16 - Tunde-Akintunde TY, Akintunde BO. Some physical properties of sesame seed. Biosystems Engineering. 2004;88:127-129
17 - Yalcın I, Özarslan C. Physical properties of vetch seed. Biosystems Engineering. 2004;88:507-512
18 - Sirisomboon P, Pornchaloempong P, Teeranud Romphophak T. Physical properties of green soybean: Criteria for sorting. Journal of Food Engineering. 2007;79:18-22
19 - Vilche C, Gely M, Santalla E. Physical properties of quinoa seeds. Biosystems Engineering. 2003;86:59-65
20 - Sacilik K, Öztürk R, Keskin R. Some physical properties of hemp seed. Biosystems Engineering. 2003;86:191-198
21 - Petrů M, Novák O, Herák D, Simanjuntak S. Finite element method model of the mechanical behaviour of Jatropha curcas L. seed under compression loading. Biosystems Engineering. 2012;111:412-421
22 - Raji AO, Favier JF. Model for the deformation in agricultural and food particulate materials under bulk compressive loading using discrete element method II: theory, compression of oilseeds. Journal of Food Engineering. 2004;64:373-380
23 - Petru M, Novak O, Herak D, Masin I, Lepsik P, Hrabe P. Finite element method model of the mechanical behaviour of Jatropha curcas L. bulk seeds under compression loading: Study and 2D modelling of the damage to seeds. Biosystems Engineering. 2014;127:50-66
24 - Terzaghi K. Die Berechnung der Durchlassigkeitsziffer des Tones aus dem Verlauf der hydrodynamischen Spannungserscheinungen. Sitzungsber, Ak.Wiss., Math. - Naturwiss. Klasse, 132,3–4 (1923), 125–138.
25 - Hasseldine BPJ, Gao C, Collins JM, Jung HD, Jang TS, Song J, Li Y. Mechanical response of common millet (Panicum miliaceum) seeds under quasi-static compression: Experiments and modeling. Journal of the Mechanical Behavior of Biomedical Materials. 2017;73:102-113
26 - Koller AA, Taylor RK, Raun WB, Weckler PR, Buser MD. Modelling and validation of maize seed orientation by pushing. Biosystems Engineering. 2016;151:338-349
27 - Horabik J, Beczek M, Mazur R, Parafiniuk P, Ryzak M, Molenda M. Determination of the restitution coefficient of seeds and coefficients of visco-elastic Hertz contact models for DEM simulations. Biosystems Engineering. 2017;161:106-119
28 - Boac JM, Ambrose RPK, Casada ME, Maghirang RG, Maier DE. Applications of discrete element method in modeling of grain postharvest operations. Food Engineering Reviews. 2014;6:128-149
29 - Mohsenin NN. Physical Properties of Plant and Animal Materials. London: Gordon and Breach Science Publishers; 1986
30 - Herak D, Kabutey A, Choteborsky R, Petru M, Sigalingging R. Mathematical models describing the relaxation behaviour of Jatropha curcas L. bulk seeds under axial compression. Biosystems Engineering. 2015;131:77-83
31 - Masin I, Petru M. Complex approach to conceptual design of machine mechanically extracting oil from Jatropha curcas L. seeds for biomass-based fuel production. BioMed Research International. 2016:77-83