1. Introduction
Wheat is one of the earliest crops cultivated by humans and can provide over 20% of the protein required by the human body, equivalent to the combined protein content of meat, eggs, and dairy products, highlighting its high nutritional value [
1]. China is the world’s largest producer and consumer of wheat. In 2023, the wheat cultivation area was 23,627 thousand hectares, with a total output of 136.59 million tons and a total trade value exceeding CNY 319.6 billion. The production and processing of wheat play a significant role in the national economy [
2]. Wheat must be processed into flour before it can be consumed, with the processing primarily divided into two main stages: cleaning and milling. The efficiency of each step in the milling process is crucial to the quality of the flour. Traditionally, milling involves grinding and breaking the cleaned wheat kernels directly, followed by separating the wheat bran from the flour through sieving, brushing, and clearing, which often results in flour with a relatively high bran content and lower quality [
3,
4,
5]. To address these issues, the peeling and milling process has been increasingly adopted in high-quality flour production. This process involves removing some of the outer layers from cleaned wheat kernels before sieving, brushing, and clearing to separate any remaining bran from the flour, resulting in a higher quality product [
6,
7,
8,
9]. This method not only enhances flour quality but also reduces production costs and improves efficiency [
10].
In the peeling powder-making process, a peeling machine is used to remove the surface skin from wheat kernels. The structure and working principle of the peeling machine are illustrated in
Figure 1. During operation, wheat kernels undergo continuous extrusion inside the peeling chamber. This continuous pressure causes gentle grinding between the wheat and the peeling plate, as well as between individual kernels, allowing for the separation of the skin layer (bran) without damaging the kernels. The bran is expelled through a screen, while the dehulled wheat kernels are discharged from the machine for subsequent processing [
11,
12]. However, if the continuous extrusion pressure is too high, prolonged exposure can cause excessive damage to the wheat kernels, leading to breakage when the accumulated damage surpasses the plastic limit. After crushing, large-sized wheat kernels may have difficulty shedding the skin layer, resulting in increased bran content in the final product and failing to achieve the desired processing outcome. Additionally, small-sized kernels, which are crushed into powder, are discharged along with the bran, causing significant material loss. Therefore, the damage to and crushing of wheat kernels under continuous extrusion are crucial factors affecting the efficacy of the peeling and milling process. The structural and operational parameters of the peeling machine are key factors influencing the extrusion pressure inside the chamber. Due to the unclear damage mechanisms of wheat kernels under continuous extrusion conditions, the design and parameter optimization of peeling machines lack a reliable basis. This leads to issues such as incomplete peeling or excessive crushing and material loss during processing [
13,
14]. Thus, understanding the damage characteristics of wheat under continuous extrusion conditions is essential for the rational design of peeling machine structures and operational parameters. This knowledge will contribute to improving processing quality, reducing costs, and enhancing efficiency.
Research on the damage characteristics of wheat kernels provides valuable theoretical insights for improving wheat processing technologies, optimizing equipment parameters, reducing processing losses, and enhancing the quality of wheat products. Early studies in this area primarily focused on identifying the types of damage wheat kernels sustain. For instance, Chen ZP et al. investigated the resistance of wheat kernels to external forces through compression and friction tests. Their findings indicated that wheat kernels exhibit better resistance to abrasion than to extrusion [
15]. Since wheat kernels are more prone to compression damage, subsequent research concentrated on factors affecting their compressive damage characteristics. Cheng X et al., for example, examined the influence of moisture content on the compression and crushing behavior of wheat kernels and proposed a model relating the compressive force to moisture content [
16]. Zhang K et al. analyzed the compressive damage characteristics across different wheat varieties and concluded that the internal structure, shape, and size of the kernels play significant roles in determining their crushing load [
17]. Similarly, Barrer GN et al. demonstrated, through uniaxial compression tests, that kernel stiffness, fracture force, and deformation are negatively correlated with moisture content [
18]. These studies underscore the significant impact of factors such as moisture content and variety on the damage characteristics of wheat kernels. However, they do not provide an in-depth analysis of the entire damage process. To bridge this gap, later research shifted to investigating the damage evolution of wheat kernels and the changes in kernel characteristics throughout each damage stage. For example, Omarov A et al. studied the deformation and crushing behavior of wheat kernels using compression tests, and by analyzing force-deformation curves, they found that elliptical wheat kernels exhibit quasi-elastic and elastic–plastic deformation behaviors under applied force [
19]. Ponce-García N et al. further identified that the damage process of wheat kernels includes a distinct elastic–plastic deformation phase, which is influenced by the kernel variety and moisture content. Kernels with a higher moisture content and softer texture displayed greater plasticity, while those with lower moisture content were more brittle [
20,
21]. Fernandes LS et al. also confirmed through mechanical testing that the compressive force decreases as the moisture content increases, while the proportional deformation modulus rises as the moisture content decreases, highlighting the relationship between water content and kernel deformation [
22]. Although these studies have advanced our understanding of wheat kernel damage characteristics, their practical relevance to optimizing dehulling parameters remains limited. This is primarily because the experiments involved single compression tests, which do not accurately reflect the continuous extrusion and grinding that wheat kernels undergo during the dehulling process. Additionally, there has been no theoretical modeling of the elastic–plastic phase in wheat kernel compression.
In the actual wheat peeling and milling process, wheat kernel damage occurs under continuous loading, manifesting as progressive damage. Therefore, studying the continuous damage characteristics of wheat kernels is of great importance for optimizing the parameter regulation in the actual wheat peeling process. To date, no research has been conducted specifically on the continuous damage characteristics of wheat kernels. However, several scholars have studied the continuous damage processes and modeled them for other materials, such as rapeseed, brown rice, and walnuts. For instance, Yuan Jiacheng et al. applied continuous damage theory to establish a compression damage evolution equation for rapeseed kernels, analyzing the cumulative compression damage process through experimental validation [
23]. Similarly, Liu X et al. used continuous damage theory to develop a continuous damage model for brown rice kernels, verifying the model through experiments and identifying the critical elastic–plastic deformation during compression as the initial state when the damage occurs [
24]. Man XL et al. established and validated a damage model for walnut kernels, analyzing the effects of particle size and impact energy on the continuous damage process, and determining the relationship between the damage accumulation coefficient and material size [
25]. These studies have successfully constructed continuous damage models for various kernels using continuous damage theory, which have been further verified and analyzed by scholars. For example, Tavares LM et al. investigated the influence of particle size and shape on continuous damage parameters and verified the accuracy of a continuum damage model for particle rupture under continuous loading conditions [
26]. Han T et al. analyzed the fatigue strength of individual particles under continuous compression by establishing an ideal fatigue model, showing that particle fatigue strength decreases with increasing compression cycles before rupture occurs [
27]. In summary, the existing research demonstrated that the continuous damage process of material particles results from the accumulation of single compression damages, and damage parameters change significantly as the load accumulates. While continuous damage theory has been widely applied to materials such as rapeseed, brown rice, and walnuts, no study has yet focused on establishing a continuous damage model for wheat kernels. Due to differences in species, shape, size, and other factors, the continuous damage processes and models for different materials vary, making models developed for other materials not entirely applicable to wheat kernels. Therefore, it is necessary to develop a specific continuous damage model tailored to wheat kernels.
This study addressed the aforementioned practical issues by investigating the continuous extrusion damage characteristics of wheat kernels. Through the analysis of damage accumulation and parameter changes during the continuous extrusion process, we developed a continuous damage model for wheat kernels. This model provides critical damage parameters, offering theoretical insights for optimizing the structure and operational parameters of wheat peeling machines and improving the effectiveness of peeling and processing technology.
2. Construction of a Continuous Damage Model for Wheat Kernels
There has been substantial research and analysis on the continuous damage characteristics of granular materials. For instance, Tavares used continuous damage theory and Hertzian contact theory to develop a theoretical model for continuous compression damage to spherical particles [
26]. This model has been widely applied in studies of the continuous damage characteristics of conventional spherical granular materials, as shown in Equation (1):
where
P is the compression load;
d is the particle diameter;
is the effective stiffness; and
δ is the compression deformation.
The effective stiffness
of the material properties of the particles is
where
D is the amount of compression damage of the particles and
k is the initial stiffness.
The damage evolution equation for compression is defined as
where
δc is the maximum deformation of particles during compression crushing and
γ is the damage index.
The Tavares continuous damage model indicates that external compression loads cause deformation in particles, leading to damage. However, according to Hertzian contact theory [
28,
29], if the particle contact deformation remains within the elastic deformation limit, the particles are in a state of complete elastic contact and will revert to their original form upon unloading, thus avoiding damage. This limitation suggests that the Tavares model has constraints in certain scenarios. Furthermore, the research by Omarov, A et al. revealed that ellipsoidal wheat kernels exhibit quasi-elastic and elastic–plastic deformation behaviors under external loads [
19]. Therefore, the Tavares continuous damage model may not be fully applicable to wheat kernels due to their ellipsoidal shape and elastic–plastic properties. In this study, we developed a continuous damage model for wheat kernels, building on continuous damage theory and the Tavares model, and making necessary adjustments to account for the specific material properties and compression damage characteristics of wheat kernels.
2.1. Correction of Material Properties of Wheat Kernels
Ellipsoidal wheat kernels undergo both elastic and plastic deformation when subjected to successive external loads. During n successive compressions, the wheat kernels accumulate damage with each compression, resulting in changes to the elastic modulus after each compression, as illustrated in
Figure 2.
Based on the actual performance behavior of these wheat kernels, the modulus of elasticity En* after the nth compression was used in this study to replace the initial stiffness k in Equation (2), thereby correcting the initial material property parameters.
The maximum radius of curvature
Rn, the minimum radius of curvature
, and the directional angle
θn of the contact area when the wheat kernels are subjected to the nth compressive load were calculated as follows:
where
Ln,
Bn, and
Hn are the length, width, and thickness of the wheat kernels after the nth compression, respectively.
The value of
cosθn can be obtained from the contact relationship tables in the ASAE S368.4 DEC2000 (R2017) standard [
30] to calculate the dimensionless coefficient
Kn. Using the Hertz contact theory formula for ellipsoids, the apparent modulus of elasticity
En of the wheat kernels after the nth compression can be calculated as follows:
where
F is the normal load;
μ is the Poisson’s ratio, and the wheat Poisson’s ratio takes a value of 0.42 [
16]; and α is the compression deformation.
The modulus of elasticity
En* of wheat kernels after the nth compression is
2.2. Correction for the Amount of Compression Damage in Wheat Kernels
When wheat kernels are in the elastic compression stage, the deformation is reversible and no kernel damage occurs. Once the compression deformation exceeds the elastic limit, the wheat undergoes plastic deformation and begins to exhibit damage. When the compression deformation reaches the maximum crushing deformation, the wheat is fully damaged and crushed. Assuming the damage of wheat kernels in the nth compression is
Dn*, and the elastic moduli of the wheat kernels in the (n − 1)th and nth compressions are
En−1* and
En*, respectively, these can be corrected using Formula (2):
The relationship between the modulus of elasticity
En* of the wheat kernels after the nth compression and the initial modulus of elasticity
E0* of the wheat kernels is
where
α is the compression deformation of the wheat kernels;
αs is the yield limit deformation; and
αc is the maximum crushing deformation of the wheat kernels.
From Equations (9) and (10), it can be observed that after n consecutive compressions:
If the compression deformation α is less than the yield limit deformation αs, then the damage variable Dn* = 0, indicating no damage to the kernels.
If the compression deformation α exceeds the yield limit deformation αs, then Dn* > 0, and the kernels sustain damage.
If the compression deformation α reaches the maximum crushing deformation αc, then Dn* = 1, indicating that the kernels are completely damaged and crushed.
Based on the Tavares continuous damage model, Equation (3) is corrected to determine the amount of compression damage sustained by wheat kernels after n consecutive compressions as follows:
2.3. Continuous Damage Modeling of Wheat Kernels
Assuming that wheat kernels are approximately homogeneous ellipsoids exhibiting both elastic and elastic–plastic deformation behaviors under external loads, this study constructed a continuous extrusion damage model for wheat kernels. To elucidate the relationship between changes in normal load and compression deformations during the continuous extrusion process, this study referred to the ASAE S368.4 DEC2000 (R2017) standard and the elastic–plastic material model developed by Brizmer V et al. based on Hertz theory.
The equivalent radius
Ren of the wheat kernels at the nth compression is
By substituting Equations (8) and (12) into the equation for determining the shape variable along the load direction in Hertz’s contact theory [
26], the mechanical compression equation for wheat kernels during the elastic compression stage of the continuous compression process can be obtained as follows:
By referring to the elastic–plastic material model developed by Brizmer V et al. [
31], along with the formula for calculating the elastic modulus of wheat kernels during the continuous damage process (Equation (9)), and the formula for the mechanical compression equation of wheat kernels in the elastic compression stage (Equation (13)) established in this paper, the mathematical model relating the normal load
F and the compression shape variable
α for each compression during the continuous damage process can be determined as follows:
4. Results and Discussion
4.1. Analysis of Continuous Compression Test Results
The results of the continuous compression test for the Ningmai 22 kernels are summarized in
Table 1. The initial deformation intervals for small, medium, and large kernels were 0 mm to 0.35 mm. Within this range, applying 30 consecutive mechanical compressions did not lead to kernel fragmentation. However, when the deformation intervals were increased to 0.35 mm to 0.60 mm, 0.35 mm to 0.65 mm, and 0.35 mm to 0.75 mm, continuous application of the compression load resulted in kernel fragmentation. Notably, the number of consecutive compressions required to cause fragmentation decreased with an increasing deformation interval. For deformation amounts exceeding 0.60 mm, 0.65 mm, and 0.75 mm, a single compression load was sufficient to directly cause fragmentation. As examples, the variation curves of the normal load with the compression set for small kernels (0.50 mm), medium kernels (0.55 mm), and large kernels (0.60 mm) during the entire process of successive damage are illustrated in
Figure 7a,
Figure 7b, and
Figure 7c, respectively.
The results of the continuous compression test for the Jichun 1 kernels are summarized in
Table 1. The initial deformation intervals for small, medium, and large kernels were 0 mm to 0.45 mm, 0 mm to 0.45 mm, and 0 mm to 0.50 mm, respectively. Within these ranges, applying 30 consecutive mechanical compressions did not result in kernel fragmentation. However, when the deformation intervals were increased to 0.45 mm to 0.70 mm, 0.45 mm to 0.80 mm, and 0.50 mm to 0.85 mm, the continuous application of compression loads led to kernel fragmentation. The number of consecutive compressions required to cause fragmentation decreased as the deformation interval increased. For deformation amounts exceeding 0.70 mm, 0.80 mm, and 0.85 mm, a single compression load was sufficient to directly cause particle crushing. For example, the curves showing the variation in the normal load with the compression deformation variable for small kernels with deformation variables of 0.60 mm and 0.75 mm, medium kernels with deformation variables of 0.65 mm and 0.85 mm, and large kernels with deformation variables of 0.70 mm and 0.90 mm are illustrated in
Figure 7d,
Figure 7e, and
Figure 7f, respectively.
The reasons for the above results are as follows: the compression and deformation process of wheat kernels generally occurs in three stages: the elastic compression stage, the plastic compression stage, and the intensive destruction stage. In the initial stage of deformation, when the compression remains within the elastic–plastic threshold, the wheat kernels undergo elastic compression. At this stage, the internal structure of the kernels remains unchanged under mechanical compression, and they retain their elasticity. Upon removal of the external load, the wheat kernels naturally return to their original shape without sustaining damage.
When the deformation caused by compression exceeds the elastic–plastic threshold, the kernels enter the plastic compression stage. In this stage, continued mechanical compression induces microscopic defects within the kernels. As the number of compression cycles increases, these defects expand and merge, ultimately leading to kernel fracture. The greater the compression deformation at this stage, the faster the defects propagate, resulting in fewer compression cycles before the kernels break.
Finally, when the deformation surpasses the maximum crushing limit, the wheat kernels reach the stage of intensified destruction. At this point, a single compression is sufficient to cause rapid initiation, expansion, and fusion of internal damage, which becomes visible on the surface as the kernels are crushed.
Combining the above analysis with the results of the study by Yuan Jiacheng et al. [
23] on the continuous damage crushing process of rapeseed kernels, which share similar elastic–plastic properties with wheat kernels, it is evident that the critical value of elastic–plastic deformation and the maximum crushing deformation are key parameters for analyzing the behavior of kernels during each stage of the continuous damage process. The performance of elastic–plastic kernels across the elastic compression, plastic compression, and intensive destruction stages remains consistent. However, due to variations in kernel type, size, moisture content, and other intrinsic attributes, the critical value of elastic–plastic deformation and the maximum crushing deformation differ significantly. Consequently, the size of the deformation zones in each of these three stages also varies.
4.2. Analysis of Continuous Compression Crushing Process of Wheat Kernels
The parameters of the compression crushing results for the wheat kernels are detailed in
Table 2 and
Table 3. The successive damage
F–
α curves for each wheat variety and size revealed that the overall relationship between the normal load and compression deformation was consistent across the different samples. Specifically, the curves exhibited a similar trend, with the slope changing once the compression deformation reaches the elastic–plastic deformation threshold. Additionally, as the number of consecutive compressions increased, the maximum normal load observed during kernel unloading decreased continuously.
When Ningmai 22 and Jichun 1 were subjected to continuous compression, the values of the equivalent radius (Ren) and contact stiffness (Kn) increased with the number of compressions. This phenomenon occurs because continuous axial mechanical compression causes the deformation of the wheat kernels to exceed the critical elastic–plastic threshold. As a result, internal damage accumulates within the kernels, which is reflected in a numerical decrease in the elastic modulus (En*). As the axial thickness of the kernels gradually decreases, their ability to recover from deformation diminishes. Simultaneously, the contact area during compression increases, leading to a corresponding rise in the equivalent radius (Ren) and contact stiffness (Kn).
Figure 8 illustrates the entire process of continuous damage and crushing of wheat kernels. Throughout this process, the wheat kernels undergo repeated compression, leading to various stages of damage [
24,
33]. The damage process can be categorized into the following stages:
- (1)
Small Deformation Stage: The initial compression causes minimal deformation with no significant damage.
- (2)
Deformation and Strengthening Stage: As the compression continues, the kernels deform and undergo strengthening due to the applied pressure.
- (3)
Crack Generation Stage: Micro-cracks start to form within the kernels as the deformation exceeds the elastic limit.
- (4)
Crack Expansion Stage: The existing cracks begin to expand under the continued pressure.
- (5)
Crack Accumulation Stage: Cracks accumulate and grow, leading to significant structural damage.
- (6)
Fusion Crushing Stage: The accumulation of cracks and continued pressure result in the complete crushing and fragmentation of the kernels.
To determine the critical value of elastic–plastic deformation (αs) for wheat kernels of various varieties and sizes, the steps were as follows:
- (1)
Calculate the Elastic Modulus: Use the previously described method to compute the modulus of elasticity (En*) of the wheat kernels for each compression. This involves measuring the modulus of elasticity after each compression event during the continuous compression process.
- (2)
Obtain the Mechanical Compression Function Curve: With the calculated elastic modulus values, substitute these into Equation (14) to generate the mechanical compression function curve. This curve represents the relationship between the normal load and the compression deformation for each compression stage.
- (3)
Determine the Critical Elastic–Plastic Deformation (αs): Compare the mechanical compression function curve with the continuous damage F–α curve obtained from experimental tests. The intersection point of these two curves represents the critical value of elastic–plastic deformation (αs) for the wheat kernels.
- (4)
Compile Material Property Parameters: Gather and record the initial elastic modulus and other relevant material property parameters for the wheat kernels under continuous compression, as shown in
Table 3.
From Equation (13), the relationship between
F and
α3/2 is linear. To simplify the compression characteristic parameters of the continuous damage process of wheat kernels, the coefficient
M is used for their definition:
By substituting the compression characteristic parameters from the continuous compression process of wheat kernels into Equation (15), the average mechanical compression function equation coefficients
M for Ningmai 22 small, medium, and large kernels were 263.75, 289.71, and 311.17, respectively, while for Jichun 1 small, medium, and large kernels, the coefficients were 213.05, 226.93, and 263.75, respectively. Consequently, the mechanical compression equations for the wheat kernels were as follows:
Based on the above mechanical compression function curves of wheat kernels and the continuous damage
F–
α curves obtained from the tests, the compression and crushing data for each variety and size of wheat kernels were summarized, as shown in
Table 4.
4.3. Continuous Compression Characterization
Based on the results of the continuous compression test on the wheat kernels, it was observed that if the deformation caused by the compression load is below the critical threshold for elastic–plastic deformation, continued loading does not result in kernel crushing. However, if the deformation exceeds this critical threshold but remains below the maximum crushing deformation, the continuous load will lead to kernel crushing. In this scenario, the number of compression cycles and the degree of kernel deformation are inversely related. Conversely, if the deformation surpasses the maximum crushing deformation, a single load will cause immediate kernel crushing. Additionally, if the number of consecutive loads exceeds the maximum crushing deformation, a single load will directly lead to kernel crushing.
Taking the small kernels of Ningmai 22 and Jichun 1 as examples, tests were conducted with compression deformation variables of 0.50 mm and 0.60 mm and compression cycles of 5 and 6, respectively (
Figure 3a,d). The parameters obtained from the continuous compression tests are summarized in
Table 5 and
Table 6.
As shown in
Table 5 and
Table 6, the mean values of the elastic–plastic critical deformation for the Ningmai 22 and Jichun 1 kernels in each successive compression were 0.37 mm and 0.45 mm, respectively. These values align with the damage limit ranges presented in
Table 1, specifically 0.35 mm to 0.40 mm and 0.45 mm to 0.50 mm. The elastic–plastic critical force in each compression cycle exhibited an increasing trend as the number of compressions increased, while the critical deformation remained relatively stable.
Figure 7a,d and
Table 5 and
Table 6 demonstrate that the
Fs–
αs3/2 ratio and the equivalent radius (
Ren) of the wheat kernels tended to increase with successive compressions, whereas the integrated modulus of elasticity (
En*) showed a declining trend. According to the continuous damage theory and the wheat kernel model described in the previous section, when the compression deformation surpasses the critical value of elastic–plastic deformation, irreversible damage occurs in the wheat kernels. With continued compression, this damage accumulates, leading to an increase in kernel deformation and a corresponding deterioration in mechanical properties. This degradation is reflected numerically by a significant increase in the equivalent radius (
Ren) and a decrease in the integrated elastic modulus (
En*), with the rate of increase in Ren notably outpacing the rate of decrease in
En* [
23,
24,
25].
Analyzing the above results of the internal microscopic changes in the kernels, it is clear that wheat kernels are not uniform isotropic bodies; they contain micro-defects, such as pore spaces. As the external load and compression increase, these micro-defects expand and merge, causing damage to the original pore structure. The pores shrink and become denser, resulting in a gradual increase in the elastic–plastic critical force of the kernels as the number of compressions rises. Since the elastic–plastic critical deformation is an inherent property of the kernels—dependent only on factors like variety, size, and moisture content, and unaffected by external conditions—it did not show significant changes in the test results.
Referring to Yuan Jiacheng et al.’s study on the continuous damage characteristics of rapeseed kernels, it was observed that during the nine consecutive compression cycles of Huayouza 62 rapeseed kernels, the elastic–plastic critical force gradually increased from the initial 7.94 N to the final crushing force of 8.76 N. Throughout these compression cycles, the elastic–plastic critical deformation remained stable at 0.33 mm, without any significant variation. While wheat kernels and rapeseed kernels differ in shape, size, and variety, both can be regarded as non-homogeneous elastic–plastic bodies. Thus, the comprehensive elastic modulus and equivalent radius of wheat kernels obtained in this experiment align closely with the continuous damage progression observed in rapeseed kernels, despite noticeable differences in the parameters measured in the experiments of Yuan Jiacheng and others [
23]. In addition, Man Xiaolan et al. also used force–deformation curve analysis to describe the changes in equivalent radius and integrated elastic modulus in walnuts subjected to a continuous damage process, and obtained the size of the critical value parameter that determines the fracture mechanism of walnuts. The change rule of the above experimental parameters is basically consistent with the previous analysis [
25].
Considering that smaller wheat kernels are more susceptible to breakage during processing, and based on the results of the continuous compression tests presented in
Table 5 and
Table 6, the following parameters are recommended for wheat peeling machines.
For processing Ningmai 22 kernels, the peeling pressure should be set to less than 62.48 N, which is the average elastic–plastic critical pressure of the kernels. Additionally, the deformation amount should be kept below 0.37 mm, the average elastic–plastic critical deformation for these kernels.
For processing Jichun 1 kernels, the peeling pressure should be maintained at less than 58.76 N, the average elastic–plastic critical pressure of these kernels. Similarly, the deformation amount should be restricted to below 0.45 mm, the average elastic–plastic critical deformation for Jichun 1 kernels.
The variation curves of the integrated modulus of elasticity (
En*) and the equivalent radius (
Ren) for Ningmai 22 and Jichun 1 under continuous compression are illustrated in
Figure 9.
Figure 9 illustrates that the integrated modulus of elasticity (
En*) of the wheat kernels of all varieties decreased progressively with an increasing number of compressions (
n). Specifically, Ningmai 22 and Jichun 1, which underwent 5 and 6 compressions, respectively, experienced crushing. At the final compression, the integrated modulus of elasticity (
En*) of Ningmai 22 and Jichun 1 was reduced to 3/5 and 1/3 of the initial value, respectively. Conversely, the equivalent radius (
Ren) of the wheat kernels of all varieties increased with the number of compressions, following a nearly quadratic growth curve. At the final compression, the equivalent radius (
Ren) of Ningmai 22 and Jichun 1 increased to 5/2 and 15/2 times their initial values, respectively. These results are consistent with the previous theoretical analysis and support the notion that the continuous compression damage evolution of wheat grains conforms to the linear cumulative damage criterion proposed by continuous damage theory.
4.4. Continuous Damage Model Validation
To further validate the accuracy of the linear cumulative damage theory for wheat kernels proposed in this study, the developed continuous damage model was subjected to validation.
Substituting the results from the continuous compression crushing tests of the wheat kernels, as detailed in
Table 5 and
Table 6, into Equation (9), the damage variables (
D) were calculated as follows: for Ningmai 22 kernels, the values were 0.027, 0.042, 0.088, and 0.279; for Jichun 1 kernels, the values were 0.015, 0.037, 0.117, 0.164, and 0.529. By incorporating these results into the continuous damage model for wheat kernels, as described by Equation (15), the
F–
α theoretical curves for Ningmai 22 and Jichun 1 at each compression stage in the damage process were obtained. These theoretical compression curves were then compared and fitted against the actual compression curves to assess the accuracy of the continuous damage model developed in this study, as illustrated in
Figure 10 and
Figure 11.
The results shown in
Figure 10 and
Figure 11 indicate that the elastic–plastic critical deformation variable (α
s) of Ningmai 22 and Jichun 1 stabilized at around 0.35 mm and 0.45 mm, respectively, during each compression cycle. This behavior aligns with the trend observed in the previous analysis of the continuous damage of wheat kernels. Additionally, the overall slope of the actual compression curves exhibited a noticeable decreasing trend, which corresponds well with the theoretical analysis of the comprehensive elastic modulus (
En*) and equivalent radius (
Ren).
To further compare the actual compression curves of Ningmai 22 and Jichun 1 with their theoretical counterparts, the maximum deviations between the actual and theoretical curves were calculated. The ratio of these deviations to the maximum normal load during crushing was determined, and the goodness-of-fit (R²) of the compression curves was also evaluated. For Ningmai 22, the maximum deviations of the actual normal load from the theoretical load during each compression were 2.14 N, 4.52 N, 2.93 N, 3.36 N, and 2.17 N, respectively. For Jichun 1, the maximum deviations were 2.31 N, 4.43 N, 3.31 N, 4.98 N, 3.46 N, and 2.98 N, respectively. In all cases, the deviations accounted for less than 5% of the crushing normal load, and the R² values for the curve fitting exceeded 0.9.
In conclusion, the continuous damage model for wheat kernels developed in this study demonstrated a certain degree of accuracy.