Test Trials and Analysis of Pod-Shattering Characteristics of Harvested Rapeseed Silique

: In order to reduce the silique shattering loss of the rapeseed mechanical harvesting process, based on the state of force on the silique during the rapeseed harvesting reel branch stage, Ningza 1810, Zhenyou 8, and Fengyou 306 were used as research objects, and the experimental research on the factors affecting rapeseed silique shattering was carried out using the swing impact method. The experimental analysis showed that rapeseed varieties, silique moisture content, silique growth position, collision material, impact speed, force position, and other factors had signiﬁcant effects on silique shattering. The impact velocity was less than 1.5 m · s − 1 , the difference in the effect of each factor on pod shattering was not signiﬁcant, and it was not easy to shatter when the moisture content of the rapeseed silique was higher. The impact resistance of the front side of rapeseed was two to four times that of the bonding surface of rapeseed petals, the shattering rate of the top rapeseed silique was twice that of the bottom siliques, and when siliques were supported, they were more likely to shatter under external forces than when they were unsupported. The experimental study of the mechanical properties of rapeseed siliques was carried out using the impending fracture method; the experimental analyses showed that the support position and force position of the silique, the loading speed, and the growth position of the silique had a signiﬁcant effect on the mechanical properties of the silique. The maximum cracking force was higher and the bending strength was stronger when the body of the silique was supported; the range of the maximum cracking force was 3.05 N to 4.16 N, and the bending strength range was 8.48 MPa to 11.57 MPa. The maximum cracking force and bending strength of the silique were stronger when the front side of the silique petal was pressurized than when the bonding surface of the petal was pressurized. Based on Pearson’s correlation and grey correlation analysis, the morphological characteristics of rapeseed siliques were ranked in order of their inﬂuence on the performance of siliques in terms of the angle between the silique and stalk, stalk diameter, petal thickness, beak length, silique thickness, silique width, and silique length. This study can be used as a reference for the design and optimization of the rapeseed harvesting reel branch mechanism and the selection of machine-harvestable rapeseed varieties.


Introduction
Rapeseed is an ancient cultivated crop, planted all over the world, with the planting area and yield second only to soybean in oil crops. China is one of the countries with the longest history of rapeseed origin and cultivation, with a perennial planting area of about 100 million acres, and the area and production account for about 20% of the global total [1][2][3][4]. A high harvesting loss of the rapeseed machine is prominent, with the loss of rapeseed siliques to pod shattering caused by the action of reels or plant pulling each other during harvesting accounting for a high percentage. Losses due to silique shattering have been reported in the literature to be 8-12% of the total seed yield, and shattering loss can

Sample Material
The rapeseed samples used in the experiment were collected from the rapeseed fields planted by Hejiaqiao Agricultural Machinery Professional Cooperative in Daitou Town, Liyang City, Jiangsu Province, from 1 to 5 June 2022, and the rapeseed varieties were Zhenyou 8, Ningza 1818, and Fengyou 306. Five plants with uniform growth and the same color and luster and free from pests and diseases were randomly selected from the field of each rapeseed variety. We cut the siliques from the main branches of rapeseed and the apical, middle, and bottom ends of the branches, without separating the petals and retaining the stalks, stored them in zip-lock bags by variety and category to reduce water loss, and marked them with marking pens.
In order to compare the differences in shattering performance between different rapeseed varieties, all three rapeseed varieties were selected for testing at the harvesting stage of the rapeseed combine, and the water contents of the seeds and siliques were chosen to be similar to each other. The moisture contents of the samples measured before the test were as follows: moisture content of Ningza 1818 siliques was 16.31%, moisture content of seeds was 18.85%, moisture content of FengYou 306 siliques was 17.15%, moisture content of seeds was 19.72%, moisture content of ZhenYou 8 siliques was 17.52%, and moisture content of seeds was 19.43%.
The morphological and geometrical characteristics of rapeseed silique are shown in Figure 1. Vernier caliper (precision 0.02 mm) was used to measure the length, beak length, silique width, silique thickness, and pericarp thickness of rapeseed siliques, and 20 siliques were measured to calculate the maximum value, minimum value, and average value. The width of the silique is the distance between the two sides of the middle of the pericarp of the silique, and the thickness of the silique is the distance between the two flaps of the plane where the pseudoseptum of the silique is located.
of the rapeseed combine, and the water contents of the seeds and siliques were chos be similar to each other. The moisture contents of the samples measured before the were as follows: moisture content of Ningza 1818 siliques was 16.31%, moisture con of seeds was 18.85%, moisture content of FengYou 306 siliques was 17.15%, moisture tent of seeds was 19.72%, moisture content of ZhenYou 8 siliques was 17.52%, and m ture content of seeds was 19.43%.
The morphological and geometrical characteristics of rapeseed silique are show Figure 1. Vernier caliper (precision 0.02 mm) was used to measure the length, beak len silique width, silique thickness, and pericarp thickness of rapeseed siliques, and 2 liques were measured to calculate the maximum value, minimum value, and ave value. The width of the silique is the distance between the two sides of the middle o pericarp of the silique, and the thickness of the silique is the distance between the flaps of the plane where the pseudoseptum of the silique is located. The geometrical parameters of the three rapeseed varieties are shown in Table 1 geometrical dimension discrepancy of different rapeseed varieties varies consider with the length of the silique of Ningza 1818 being the longest, with a mean value of mm, the width of the silique of Zhenyou 8 being the largest, with a mean value of mm, and the thickness of the silique of Ningza 1818 being the largest, with a mean v of 3.37 mm.  The geometrical parameters of the three rapeseed varieties are shown in Table 1. The geometrical dimension discrepancy of different rapeseed varieties varies considerably, with the length of the silique of Ningza 1818 being the longest, with a mean value of 68.6 mm, the width of the silique of Zhenyou 8 being the largest, with a mean value of 5.23 mm, and the thickness of the silique of Ningza 1818 being the largest, with a mean value of 3.37 mm. The swing impact method was adopted to simulate the collision impact of reel operation process on the siliques during rapeseed harvesting, and the collision impact test was carried out on the single-constrained and supported-constrained siliques to analyze the effects of the different collision speeds, collision positions, and other parameters on the pod shattering of rapeseed siliques in the harvesting time. Due to the small silique stalk and support friction, external load force can easily cause it to slide down, and therefore, medical tape was wrapped 3 to 4 times around the silique stalk as part of the enhancement of the Appl. Sci. 2023, 13, 9369 4 of 19 friction between the silique stalk and the fixture, as well as to prevent fixture clamping and fixing the process of the silique stalk caused by the damage.

Test Equipment and Methods
Based on the self-designed swing silique collision impact test device, as shown in Figure 2, the stalk of the rapeseed silique was fixed with a fixture, and by releasing the impact parts and counterweights of different heights, different collision impact speeds were realized to hit the silique, and the number of siliques shattered was recorded in the test. It was found in a previous preliminary test that if twenty rapeseed siliques were used in each group of tests, more accurate results could be obtained. ation process on the siliques during rapeseed harvesting, and the collision impact test carried out on the single-constrained and supported-constrained siliques to analyz effects of the different collision speeds, collision positions, and other parameters on pod shattering of rapeseed siliques in the harvesting time. Due to the small silique and support friction, external load force can easily cause it to slide down, and there medical tape was wrapped 3 to 4 times around the silique stalk as part of the enhancem of the friction between the silique stalk and the fixture, as well as to prevent fixture cla ing and fixing the process of the silique stalk caused by the damage.
Based on the self-designed swing silique collision impact test device, as show Figure 2, the stalk of the rapeseed silique was fixed with a fixture, and by releasing impact parts and counterweights of different heights, different collision impact sp were realized to hit the silique, and the number of siliques shattered was recorded in test. It was found in a previous preliminary test that if twenty rapeseed siliques were in each group of tests, more accurate results could be obtained. The number of pod-shattering siliques was counted using the following criteria petals were cracked more than 2 mm from the stalk, the petals were separated from pseudoseptum on at least one side, or the rapeseed could be seen flying away and fa off. The silique shattering rate was calculated as follows: where p -silique shattering rate, x -number of siliques ruptured in the experiment n -total number of rapeseed siliques used in the experiment.

Impending Fracturing Test
A universal material testing machine was used to carry out the mechanical prope of rapeseed silique bending test, to study the silique in the plant in different growth tions, silique force positions, support modes, loading speeds, etc., and silique mecha properties of the direct relationship between the siliques. The rapeseed silique impen fracturing test was carried out in the panicle crop harvesting laboratory of Nanjing A cultural Mechanisation Research Institute (NIAM), and the test equipment used wa WDW-10 microcomputer-controlled electronic universal testing machine, which has a range of 0 to 1 kN, a maximum traveling stroke of 700 mm, and a displacement sensor force sensor accuracy of ±0.1%.
During the test, we first set the downward loading speed of the universal test te apparatus to 15 mm·min −1 in the computer software and then placed the silique sam on the support of the three-point bending test fixture in accordance with the test req ments to ensure that the probe of the testing machine was located at the midpoint o The number of pod-shattering siliques was counted using the following criteria: the petals were cracked more than 2 mm from the stalk, the petals were separated from the pseudoseptum on at least one side, or the rapeseed could be seen flying away and falling off. The silique shattering rate was calculated as follows: where p-silique shattering rate, x-number of siliques ruptured in the experiment, and n-total number of rapeseed siliques used in the experiment.

Impending Fracturing Test
A universal material testing machine was used to carry out the mechanical properties of rapeseed silique bending test, to study the silique in the plant in different growth positions, silique force positions, support modes, loading speeds, etc., and silique mechanical properties of the direct relationship between the siliques. The rapeseed silique impending fracturing test was carried out in the panicle crop harvesting laboratory of Nanjing Agricultural Mechanisation Research Institute (NIAM), and the test equipment used was the WDW-10 microcomputer-controlled electronic universal testing machine, which has a test range of 0 to 1 kN, a maximum traveling stroke of 700 mm, and a displacement sensor and force sensor accuracy of ±0.1%.
During the test, we first set the downward loading speed of the universal test tensile apparatus to 15 mm·min −1 in the computer software and then placed the silique samples on the support of the three-point bending test fixture in accordance with the test requirements to ensure that the probe of the testing machine was located at the midpoint of the silique. We adjusted the position of the tensile tester with the sensor so that it was at the critical point of contact with the silique and started the preload < 0.5 N, and the span L of the support was adjusted according to the test factors. Before the start of each group of tests, one needs to zero the displacement and load parameters to avoid human error. We manually clicked the run button, the universal test tensile tester began to run slowly downward, and the software showed the bending deformation and load force change curve which slowly increased. At the end of the test, the change curve data in the software were exported, and the change curve of the compression load with bending deformation of rapeseed silique was plotted using MATLAB R2018b software.
The rapeseed variety, silique stress surface, load loading speed, and silique support position were selected as the test factors, and the silique cracking force and the bending strength were used as the test evaluation indexes. Silique cracking force refers to the load force corresponding to the highest point at the straight line segment presenting linear elastic deformation with the increase in bending deformation, which is also the peak force of the sudden drop fluctuation of the pressure load when the first crack or petal crack appears in the silique.
The bending strength of the silique was calculated as follows: The cross-section of the rapeseed silique has an approximately elliptical circular shape, the direction of the silique width is defined as the x-axis, and its moment of inertia to the x-axis is where a 1 -large cross-section long semi-axis length, that is, half of the width of the silique, mm; b 1 -large cross-section short semi-axis length, that is, half of the thickness of the silique, mm; a 2 -small cross-section long semi-axis length of the small section, that is, half of the width of the silique minus the thickness of the petal, mm; and b 1 -length of the short semi-axis of the small section, that is, half of the thickness of the silique minus the thickness of the petal, mm.
Silique width direction bending test process of the maximum bending moment is where F x−max -maximum bending force, N; and L-distance between two supports, mm. Bending section modulus in the direction of the width of the silique: Bending strength in the direction of the width of the silique: Similarly, the silique thickness direction is defined as the y-axis, and its moment of inertia to the y-axis is The maximum bending moment during the bending test in the silique thickness direction is where F y−max -maximum bending force, N; and L-distance between the two supports, mm. The bending section modulus in the direction of the thickness of the silique is The flexural strength of the silique is Appl. Sci. 2023, 13, 9369 6 of 19

Swing Impact Test Results and Analysis
The growth state of rapeseed siliques on the plant can be regarded as a cantilever beam structure, which deforms and shakes freely under the action of the collision impact force, and the length of its stalk directly affects the state of silique rupturing and shattering. It was found in a previous preliminary test that the position of the stalk was 2 mm from the connection point, and there was a significant difference between the silique shattering with different collision impact speeds. In order to differentiate the effects of different factors on silique shattering, the position of the stalk was selected to be 2 mm from the connection point as a prerequisite for the test of the swing impact test.

Swing Impact Test Results and Analysis
The growth state of rapeseed siliques on the plant can be regarded as a can beam structure, which deforms and shakes freely under the action of the collision force, and the length of its stalk directly affects the state of silique rupturing and s ing. It was found in a previous preliminary test that the position of the stalk was from the connection point, and there was a significant difference between the siliqu tering with different collision impact speeds. In order to differentiate the effects of ent factors on silique shattering, the position of the stalk was selected to be 2 mm fr connection point as a prerequisite for the test of the swing impact test.  In swing impact tests with the SUS304 material, the shattering rate was 4 Zhenyou 8 and 45% for Fengyou 306. With the REEK nylon material, the shatteri of Zhenyou 8 was 30%, and the shattering rate of Fengyou 306 was 40%. With d types of collision materials, the shattering rate of Fengyou 306 was higher than Zhenyou 8. The differences in shattering caused by different material types were a nificant, and the shattering rates of the same rapeseed varieties with REEK nylon m were smaller than those of SUS304 materials, indicating that both rapeseed variet collision material types had a significant effect on the shattering of rapeseed siliqu

Impact Velocity and Stress Surface Effects on Silique Shattering
During rapeseed harvesting operations, reel teeth contact with siliques instant duces collision impact force because rapeseed silique shells consist of two shell-like and two linear petals. In order to analyze the impact of the rapeseed varieties and s of different stress surfaces by the collision of external forces on silique shattering, different collision impact speed on the Ningza 1818 and Zhenyou 8 rapeseed va from the right ahead and the bonding face of the silique petals processed, an imp In swing impact tests with the SUS304 material, the shattering rate was 40% for Zhenyou 8 and 45% for Fengyou 306. With the REEK nylon material, the shattering rate of Zhenyou 8 was 30%, and the shattering rate of Fengyou 306 was 40%. With different types of collision materials, the shattering rate of Fengyou 306 was higher than that of Zhenyou 8. The differences in shattering caused by different material types were also significant, and the shattering rates of the same rapeseed varieties with REEK nylon materials were smaller than those of SUS304 materials, indicating that both rapeseed varieties and collision material types had a significant effect on the shattering of rapeseed siliques.

Impact Velocity and Stress Surface Effects on Silique Shattering
During rapeseed harvesting operations, reel teeth contact with siliques instantly produces collision impact force because rapeseed silique shells consist of two shell-like petals and two linear petals. In order to analyze the impact of the rapeseed varieties and siliques of different stress surfaces by the collision of external forces on silique shattering, with a different collision impact speed on the Ningza 1818 and Zhenyou 8 rapeseed varieties, from the right ahead and the bonding face of the silique petals processed, an impact test was conducted; the point of impact is the central part of the rapeseed silique, and the results of the silique shattering rate test are shown in Table 2. According to the test data shown in Table 2, when the other factors were the same, Ningza 1818 and Zhenyou 8 showed that when subjected to the same collision impact, collision with the bonding surface of the petal compared to the right ahead of the petal is more likely to lead to silique shattering, and with the increasing speed of the collision impact, the silique shattering rate with the increase in the speed of the collision impact is increased. Table 2 also shows that the increase in impact velocity will cause the silique shattering rate to increase, but the increase does not have a linear proportional relationship: when the collision impact velocity is small, with different collision stress surfaces of siliques, the silique shattering rate difference is small, and when the impact velocity increases, with different collision stress surfaces of siliques, the silique shattering rate difference is obvious. Between different varieties with the collision impact speed greater than 2.5 m·s −1 , the silique shattering rate difference is also obvious, and when the collision impact speed is lower, the silique shattering rate difference is not obvious.

Effect of Moisture Content on Shattering in Siliques
In order to analyze the effect of moisture content on silique shattering, with Ningza 1818 as the research object, in the rapeseed maturity period, the silique at the middle of the main branch was cut. We selected the moisture content of the three levels. The bonding surface of the rapeseed silique petal was used as the collision force surface. The collision impact speed was 3.5 m·s −1 . The collision impact material was SUS304. Using the collision of the siliques at the middle part, the measured experimental data of the rapeseed silique shattering rate are shown in Table 3. From the data in Table 3, it can be seen that the shattering rate of rapeseed siliques with different moisture content is also different; the shattering rate of siliques decreases gradually with the increase in moisture content, and the shattering rate of siliques with 20.61% moisture content is 50% lower than the shattering rate of siliques with 9.24% moisture content. Rapeseed siliques are not easy to shatter under the condition of high moisture content.

Influence of Silique Growth Position on Silique Shattering
In order to analyze the shattering influence of the different growth positions of rapeseed siliques on the plant, with Ningza 1818 rapeseed varieties as the research object and the rapeseed silique petal bonding surface as the collision force surface, the collision impact speed was 3.5 m·s −1 , and the collision impact test of the middle of the silique was carried out at the top, middle, and bottom of the main rapeseed branch. The test results are shown in Table 4. Table 4. Experimental data of different growth positions of silique on shattering rate.

Bottom of Main Branch
Shattering rate/% 50 40 25 The shattering rates of siliques at the top, middle, and bottom of the main rapeseed branches were 50%, 40%, and 25%, respectively. The shattering rates of siliques at the top and middle of the branch were close to each other but differed significantly from those at the bottom of the branch. Under the same external force, the shattering rate of the silique from the top to the bottom of the main rapeseed branch gradually decreased. The reason is that during the rapeseed maturity period, the maturity degree of the silique at the top, middle, and bottom of the same plant is not consistent, and siliques with low maturity are strong and are not easy to shatter, while the moisture content of rapeseed siliques with relatively high maturity is on the low side, the brittleness is more obvious when subjected to the same collision impact force, and plastic deformation leads to silique shattering.

The Effect of Supporting Role on Silique Shattering
Rapeseed has many branches that are intertwined with each other, and each silique layer also overlaps. In the process of rapeseed harvesting reeling and branching operations, when reeling rapeseed plants, due to cross pulling and other reasons, it can be easily seen that the siliques are blocked by other plants to support the role of the siliques, so the study of the impact of the external supporting role on silique shattering to reduce silique shattering in reeling and branching has reference value.
Zhenyou 8 rapeseed was selected as the research object, the silique at the middle position of the rapeseed plant was selected, the clamping position was 2 mm away from the carpopodium, the collision impact material was SUS304, the collision impact speed was 3.5 m·s −1 , the silique force position was the middle part of the fruit flap bonding surface, in order to simulate the different blocking positions, three distances (0 mm, 5 mm, 10 mm) were set at the horizontal position of the silique, and a ∅3.86 mm rapeseed branch was used as blocking material; the test is shown in Figure 4. the rapeseed silique petal bonding surface as the collision force surface, the collision impact speed was 3.5 m·s −1 , and the collision impact test of the middle of the silique was carried out at the top, middle, and bottom of the main rapeseed branch. The test results are shown in Table 4. Table 4. Experimental data of different growth positions of silique on shattering rate.

Middle of Main Branch Bottom of Main Branch
Shattering rate/% 50 40 25 The shattering rates of siliques at the top, middle, and bottom of the main rapeseed branches were 50%, 40%, and 25%, respectively. The shattering rates of siliques at the top and middle of the branch were close to each other but differed significantly from those at the bottom of the branch. Under the same external force, the shattering rate of the silique from the top to the bottom of the main rapeseed branch gradually decreased. The reason is that during the rapeseed maturity period, the maturity degree of the silique at the top, middle, and bottom of the same plant is not consistent, and siliques with low maturity are strong and are not easy to shatter, while the moisture content of rapeseed siliques with relatively high maturity is on the low side, the brittleness is more obvious when subjected to the same collision impact force, and plastic deformation leads to silique shattering.

The Effect of Supporting Role on Silique Shattering
Rapeseed has many branches that are intertwined with each other, and each silique layer also overlaps. In the process of rapeseed harvesting reeling and branching operations, when reeling rapeseed plants, due to cross pulling and other reasons, it can be easily seen that the siliques are blocked by other plants to support the role of the siliques, so the study of the impact of the external supporting role on silique shattering to reduce silique shattering in reeling and branching has reference value.
Zhenyou 8 rapeseed was selected as the research object, the silique at the middle position of the rapeseed plant was selected, the clamping position was 2 mm away from the carpopodium, the collision impact material was SUS304, the collision impact speed was 3.5 m·s −1 , the silique force position was the middle part of the fruit flap bonding surface, in order to simulate the different blocking positions, three distances (0 mm, 5 mm, 10 mm) were set at the horizontal position of the silique, and a ∅3.86 mm rapeseed branch was used as blocking material; the test is shown in Figure 4. Experimental data of the silique shattering rate are shown in Table 5. Experimental data of the silique shattering rate are shown in Table 5. From the data in Table 5, it can be seen that when the moisture content of rapeseed siliques was 16.32%, the silique shattering rate was 40% under the unsupported effect. When there was support, the shattering rate was 75~100%, and the shattering rate was significantly higher than the shattering rate under unsupported conditions, and the effect of support distance on shattering was also extremely significant, with the shattering rate being 100% when the support point was 0 mm away from the horizontal distance of the top of the silique. When there was support, the reason for the high shattering rate is mainly that the upper end of the silique is blocked by another object; at this time, the state of the silique can be regarded as a simply supported beam, and the attachment state of the silique itself on the plant can be regarded as a cantilever beam structure. Under the action of the same external force, the center of the simply supported beam is more likely to be bent and broken. Under the action of two-end support, it is more likely to lead to silique shattering when the rapeseed silique is subjected to external force.

Effect of Support Position on Mechanical Properties of Silique
Taking Zhenyou 8 siliques as the research object, the siliques at the middle of the main rapeseed branch were selected, and the loading speed was 15 mm·min −1 . Rapeseed silique support positions were divided into the following three categories: 1 carpopodium and horn beak placed at both ends of the support, and rapeseed silique petals are completely in the state of stagnation; 2 from the horn beak closest to the end of the petals placed at one end of the support, and the rapeseed carpopodium is at the other end of the support; and 3 the rapeseed silique body is placed at both ends of the support. The corresponding spans L of the two supports were 60 mm, 55 mm, and 40 mm, respectively, with the loading position at the center of the support and the stress position at the front of the silique petals, as shown in Figure 5. From the data in Table 5, it can be seen that when the moisture content of rapeseed siliques was 16.32%, the silique shattering rate was 40% under the unsupported effect. When there was support, the shattering rate was 75~100%, and the shattering rate was significantly higher than the shattering rate under unsupported conditions, and the effect of support distance on shattering was also extremely significant, with the shattering rate being 100% when the support point was 0 mm away from the horizontal distance of the top of the silique. When there was support, the reason for the high shattering rate is mainly that the upper end of the silique is blocked by another object; at this time, the state of the silique can be regarded as a simply supported beam, and the attachment state of the silique itself on the plant can be regarded as a cantilever beam structure. Under the action of the same external force, the center of the simply supported beam is more likely to be bent and broken. Under the action of two-end support, it is more likely to lead to silique shattering when the rapeseed silique is subjected to external force.

Effect of Support Position on Mechanical Properties of Silique
Taking Zhenyou 8 siliques as the research object, the siliques at the middle of the main rapeseed branch were selected, and the loading speed was 15 mm·min −1 . Rapeseed silique support positions were divided into the following three categories: ① carpopodium and horn beak placed at both ends of the support, and rapeseed silique petals are completely in the state of stagnation; ② from the horn beak closest to the end of the petals placed at one end of the support, and the rapeseed carpopodium is at the other end of the support; and ③ the rapeseed silique body is placed at both ends of the support. The corresponding spans L of the two supports were 60 mm, 55 mm, and 40 mm, respectively, with the loading position at the center of the support and the stress position at the front of the silique petals, as shown in Figure 5.  The moisture content of rapeseed siliques was measured to be 14.8%, and that of seeds was 16.31% at the time of the test. The variation curves of the load applied to the siliques at different support positions with bending deformation are shown in Figure 6. The moisture content of rapeseed siliques was measured to be 14.8%, and that of seeds was 16.31% at the time of the test. The variation curves of the load applied to the siliques at different support positions with bending deformation are shown in Figure 6.
In the initial stage, the load increases with the increase in the bending deformation of the silique, which is approximately linear, and when the peak value appears for the first time, the silique begins to yield, and the bonding of the petal and the pseudoseptum begins to crack. At this time, it is the maximum yield load of the silique, and it is also the maximum cracking force of the silique. With the increase in the displacement, the load is essentially stable, the silique cracks are expanding, and the two petals of the silique begin to slowly undergo bending deformation. Finally, as the displacement continues to increase, the load gradually decreases and fluctuates around a certain value, the petals separate, and the seeds fall off in this process. When the body of the rapeseed silique is supported, the load curve shows two peaks, and the load of the second peak is much higher than the first one, which is because the yield critical value of the silique petals and false septum is reached at the second peak. The experimental data on the maximum cracking force and flexural strength of the silique based on the support position of the rapeseed silique are shown in Table 6. In the initial stage, the load increases with the increase in the bending deformation of the silique, which is approximately linear, and when the peak value appears for the first time, the silique begins to yield, and the bonding of the petal and the pseudoseptum begins to crack. At this time, it is the maximum yield load of the silique, and it is also the maximum cracking force of the silique. With the increase in the displacement, the load is essentially stable, the silique cracks are expanding, and the two petals of the silique begin to slowly undergo bending deformation. Finally, as the displacement continues to increase, the load gradually decreases and fluctuates around a certain value, the petals separate, and the seeds fall off in this process. When the body of the rapeseed silique is supported, the load curve shows two peaks, and the load of the second peak is much higher than the first one, which is because the yield critical value of the silique petals and false septum is reached at the second peak. The experimental data on the maximum cracking force and flexural strength of the silique based on the support position of the rapeseed silique are shown in Table 6. The mean value of the maximum cracking force of the silique when it was at the position of support at the carpopodium and beak was 1.46 N, the standard deviation was 0.23, and the coefficient of variation was 15.76%; the mean value of the flexural strength of the silique was 6.09 MPa, the standard deviation was 0.95, and the coefficient of variation was 15.61%. The mean value of the maximum cracking force of the silique when it was at the position of the support of the carpopodium and the petals was 2.13 N, the standard deviation was 0.18, and the coefficient of variation was 8.44%; the mean value of the bending strength of the silique was 8.15 MPa, the standard deviation was 0.7, and the coefficient of variation was 8.56%. The mean value of the maximum cracking force of the silique at the position of silique body support was 3.44 N, the standard deviation was 0.44, and the coefficient of variation was 12.86%; the mean value of the flexural strength of the  The mean value of the maximum cracking force of the silique when it was at the position of support at the carpopodium and beak was 1.46 N, the standard deviation was 0.23, and the coefficient of variation was 15.76%; the mean value of the flexural strength of the silique was 6.09 MPa, the standard deviation was 0.95, and the coefficient of variation was 15.61%. The mean value of the maximum cracking force of the silique when it was at the position of the support of the carpopodium and the petals was 2.13 N, the standard deviation was 0.18, and the coefficient of variation was 8.44%; the mean value of the bending strength of the silique was 8.15 MPa, the standard deviation was 0.7, and the coefficient of variation was 8.56%. The mean value of the maximum cracking force of the silique at the position of silique body support was 3.44 N, the standard deviation was 0.44, and the coefficient of variation was 12.86%; the mean value of the flexural strength of the silique was 9.51 MPa, the standard deviation was 1.24, and the coefficient of variation was 13.03%.
Silique support position had a significant effect on both the maximum cracking force and bending strength; when supported at the carpopodium and beak, the load was very easy to bend and deform, so the measured maximum cracking force and bending strength were small, with the range of the maximum cracking force of 1.11 N to 1.73 N and the bending strength range of 4.62 MPa to 7.21 MPa. When the silique body was supported at both ends of the support point, the bending capacity and bending strength of the silique became larger, and the maximum cracking force required to shatter was the largest, with the range of the maximum cracking force of 3.05 N to 4.16 N and the bending strength range of 8.48 MPa to 11.57 MPa.

Effect of Stress on the Mechanical Properties of Siliques
In order to analyze the mechanical properties of the different stress surfaces of rapeseed siliques under pressure conditions, three rapeseed varieties, Zhenyou 8, Ningza 1818, and Fengyou 306, were loaded uniformly from the front side of the petals and the bonding surface at a rate of 5 mm·min −1 . Due to the length of different rapeseed siliques and the influence of the different support positions on the mechanical properties of the siliques, in order to check if the results were comparable, silique carpopodium and beak ends were placed on the support for the test conditions, according to the length of the rapeseed silique, which was adjusted to support the span of the seat. The Zhenyou 8 corresponding span was 40 mm, the Ningza 1818 corresponding span was 53 mm, and the Fengyou 306 corresponding span was 50 mm. The test is shown in Figure 7, and the measured maximum cracking force and flexural strength of the different stress surfaces of siliques are shown in Table 7.
became larger, and the maximum cracking force required to shatter was the largest, the range of the maximum cracking force of 3.05 N to 4.16 N and the bending stre range of 8.48 MPa to 11.57 MPa.

Effect of Stress on the Mechanical Properties of Siliques
In order to analyze the mechanical properties of the different stress surfaces of seed siliques under pressure conditions, three rapeseed varieties, Zhenyou 8, Ningza and Fengyou 306, were loaded uniformly from the front side of the petals and the bon surface at a rate of 5 mm·min −1 . Due to the length of different rapeseed siliques an influence of the different support positions on the mechanical properties of the sili in order to check if the results were comparable, silique carpopodium and beak ends placed on the support for the test conditions, according to the length of the rapese lique, which was adjusted to support the span of the seat. The Zhenyou 8 correspon span was 40 mm, the Ningza 1818 corresponding span was 53 mm, and the Fengyo corresponding span was 50 mm. The test is shown in Figure 7, and the measured m mum cracking force and flexural strength of the different stress surfaces of silique shown in Table 7.
(1) petal front (2) petal bonding surface  From the experimental data in Table 7, it can be seen that for different rapeseed eties, the maximum cracking force and flexural strength of siliques were different results of the homogeneity test of variance that was carried out on the data are show Table 7, the maximum cracking force of siliques was p = 0.082 > 0.05, the flexural stre  From the experimental data in Table 7, it can be seen that for different rapeseed varieties, the maximum cracking force and flexural strength of siliques were different. The results of the homogeneity test of variance that was carried out on the data are shown in Table 7, the maximum cracking force of siliques was p = 0.082 > 0.05, the flexural strength was p = 0.412 > 0.05, and the data show a normal distribution, which meets the requirement of analysis of variance. The analysis of variance of the maximum cracking force of siliques under the different rapeseed varieties and different stress surfaces impending fracturing test is shown in Table 8. Note: *** indicate significant correlation at 0.05 level.
As shown in Table 8, the effects of rapeseed varieties (p < 0.001) and the different stress surfaces of siliques (p < 0.001) on the maximum cracking force of siliques were highly significant. Using the SNK method for multiple comparisons of rapeseed varieties [30], the maximum cracking force of Zhenyou 8 was significantly different from the two varieties of Ningza 1818 and Fengyou 306, and the difference between the maximum cracking force of Fengyou 306 and Ningza 1818 was non-significant. For the same rapeseed variety, the maximum cracking force of the petal bonding surface was less than that of the petal front surface, and the difference was significant.
The analysis of variance of the flexural strength of siliques under rapeseed varieties and different stress surfaces is shown in Table 9. The influence of rapeseed varieties with different stress surfaces on the flexural strength of the silique was highly significant. Using the SNK method for multiple comparisons of rapeseed varieties, the flexural strength of Ningza 1818 was significantly different from the two varieties of Ningza 1818 and Fengyou 306, the difference between the flexural strength of Fengyou 306 and Zhenyou 8 was non-significant, and the flexural strength of the petal bonding surface was less than that of the front side of the petals, and the difference was significant.

Effect of Loading Speed and Silique Growth Position on Mechanical Properties of Siliques
Zhenyou 8 was used as the research object, and 12 siliques each were selected from the top, middle, and bottom of the main rapeseed branch. The previous analysis showed that the stability of the silique carpopodium and beak support was poor, and the silique was easier to rupture and shatter when the bonding surface of the petals was subjected to force. The silique body was selected as the support position and the direction of the front of the silique petals to be loaded at the loading speeds of 5 mm·min −1 , 15 mm·min −1 , and 25 mm·min −1 , respectively. Before the start of the test, the water content of the silique was measured to be 14.8%, the water content of the seed was 16.31%, and the test data are shown in Table 10. From the data in Table 10, it can be seen that the maximum cracking force and flexural strength required for silique shattering at the middle of the rapeseed plant were the greatest, and the maximum cracking force and flexural strength of the silique at the middle of the rapeseed plant were also greater than those of the bottom and top siliques at the same loading speed. The average values of the maximum cracking force and flexural strength of siliques were maximum at a loading speed of 15 mm·min −1 .
The homogeneity test of variance for maximum cracking force (p = 0.304 > 0.05) and flexural strength (p = 0.282 > 0.05) was performed on the test results under the distribution position of siliques and different loading speed treatments, respectively, and the data showed normal distribution, which satisfied the requirement of analysis of variance. The analysis of variance of the maximum cracking force of siliques under different loading speeds and silique growth positions is shown in Table 11. Table 11. Analysis of variance of silique position and loading speed on the maximum cracking force of siliques.

Source of Variation
Type III Sum of Squares The p-values for both locations of silique growth and loading speed were less than 0.01, indicating that the effect of the two factors on the maximum cracking force of rapeseed siliques was extremely significant within 99% of the confidence interval. Multiple comparisons of the maximum cracking force of siliques under different distribution positions of siliques were carried out using the SNK method. The maximum cracking force of siliques at the middle of the main rapeseed branch was significantly different from that of siliques in other positions, and the difference in the maximum peak force of shattering between siliques at the bottom and top of the main rapeseed branch was not significant. Multiple comparisons of load loading speeds showed that all three levels of loading speeds had significant differences in the maximum cracking force, and the effects on the maximum cracking force of silique shattering were ranked as follows: medium speed (15 mm·min −1 ) > high speed (25 mm·min −1 ) > low speed (5 mm·min −1 ).

Degree of Freedom
The analysis of variance of the flexural strength of siliques under different loading speeds and silique growth positions is shown in Table 12. The model significance (p < 0.001) of the flexural strength of siliques indicated that the silique distribution location (p = 0.008 < 0.01) and the load loading speed (p < 0.001) had a highly significant effect on flexural strength. Multiple comparisons of the flexural strength of siliques under different distribution positions of siliques were carried out using the SNK method. The flexural strength of siliques at the middle of the main rapeseed branch was significantly different from that of siliques in other positions, and the difference in flexural strength between siliques at the bottom and top of the main rapeseed branch was not significant. Multiple comparisons of load loading speeds were made, and it was seen that there was a significant difference between the three loading speeds on the flexural strength of the silique, with the effects on the flexural strength of the silique ranked as follows: medium speed (15 mm·min −1 ) > high speed (25 mm·min −1 ) > low speed (5 mm·min −1 ).

Correlation Analysis between Silique Morphological Characters with Silique Shattering
Taking Ningza 1818 rapeseed siliques as the research object, the Vernier caliper (precision 0.02 mm) was used to measure seven silique traits, such as silique length, beak length, silique middle width, silique middle thickness, pericarp thickness, the diameter of the carpopodium, the angle between the carpopodium and the silique, etc., and the statistical data of the silique traits of Ningza 1818 rapeseed siliques were measured and are shown in Table 13. Using the swing impact test method, the dynamic impact test was carried out at the middle position of the frontal part of the petal silique, with a collision impact speed of 3.5 m·s −1 , the collision impact material was SUS304, and the silique shattering corresponding to each morphological feature of the silique was recorded. Pearson's correlation analysis was used to study the correlation between 90 rapeseed siliques' shattering and silique morphological characteristics. Silique shattering situations were divided into two cases, with 0 representing no silique shattering and 1 representing silique shattering; it was a dichotomous variable and was analyzed by point-two series correlation analysis. The results of correlation analyses between silique shattering and morphological characters are shown in Table 14. Note: * and ** indicate significant correlation at 0.05 level (two-tailed) and highly significant correlation at 0.01 level (two-tailed), respectively.
As shown in Table 14, there were different degrees of correlation between the silique shattering situation and the morphological characteristics of the silique itself. There was a highly significant positive correlation between silique length, beak length, and silique width, a highly significant positive correlation between silique thickness and petal thickness and silique width, a significant positive correlation between silique thickness and beak length and silique width and pericarp thickness, and a highly significant positive correlation between stalk diameter and silique width, silique thickness, and petal thickness, and a significant positive correlation between petiole diameter and silique length. Pearson's correlation coefficients between the various morphological characters of siliques revealed that longer siliques had a greater chance of having a longer beak, a greater width, and a thicker stalk diameter, while thicker siliques had a greater chance of having a wider width, a longer beak, a thicker flap, and a thicker stalk diameter.
A point-two series correlation analysis between silique shattering and the morphological characteristics of siliques, as measured by the swing impact test, showed that there was a significant positive correlation between shattering and beak length (point-two series correlation coefficient of 0.226) and a highly significant positive correlation with stalk diameter (point-two series correlation coefficient of 0.594, which was moderately correlated), provided that all other influencing factors were constant.

Grey Correlation Analysis between Silique Shattering with Silique Morphological Characters
In order to further investigate the correlation between the silique morphological characteristics and whether or not the siliques were shattered, grey correlation analysis was carried out using IBM SPSS 20 software for seven morphological parameters of the silique (length of the silique/mm, length of the beak/mm, width of the middle part of the silique/mm, thickness of the middle part of the silique/mm, thickness of the pericarp/mm, diameter of the silique stalk/mm, and the angle of the silique stalk to the silique/ • ), as well as for the data of 90 experiments. The correlation between the seven morphological parameters and shattering was also investigated by considering whether the siliques were shattered or not as the "reference value" (parent series). When using grey correlation analysis, the resolution coefficient was taken as 0.5, the value of the correlation coefficient was calculated by combining it with the formula of the correlation coefficient, and finally, the value of correlation was calculated, and the results of the analysis are shown in Table 15. Using grey correlation analysis, it was concluded that the influence of silique morphological characteristics on the silique shattering situation was in the following order: the angle of the silique stalk to the silique > carpopodium diameter > pericarp thickness > beak length > silique thickness > silique broad surface > silique length. The morphological characteristics of rapeseed siliques can be used as the appearance index for screening shatter-resistant rapeseed and provide a certain reference for selecting and breeding rapeseed varieties with strong shatter resistance to machine harvesting.

Discussion
Rapeseed siliques are susceptible to pod shattering which makes them extremely unfavorable to mechanized harvesting, limiting the further improvement of production efficiency, so based on the actual force state of rapeseed siliques during machine harvesting, carrying out the analysis of the shatter resistance of different varieties of rapeseed siliques, the selection of rapeseed varieties suitable for machine harvesting [31], and the optimization of the structure of the reel and separate mechanism and working parameters, to reduce the loss of rapeseed mechanical harvesting, has very important theoretical and practical significance.
In this paper, we chose the main promoted rapeseed varieties in the Chinese Yangtze River, China, Zhenyou 8, Ningza 1818, and Fengyou 306, and used the swinging impact method and the impending fracturing method to carry out the rapeseed pod-shattering and fracturing tests, respectively. For the same variety of rapeseed siliques, the top, middle, and bottom of the rapeseed silique of the main branch were subdivided. In the swing impact and impending fracturing tests, respectively, according to the biological characteristics of rapeseed siliques, the stress surface was divided into the front surface and the bonding surface of rapeseed siliques.
Under the same stress conditions, the pod-shattering rate of REEK nylon material was reduced by 5% to 10%. REEK nylon material can be used to reduce pod shattering in the process of reeling and separating in rapeseed harvesting. The impact resistance of the front side of rapeseed was two to four times that of the bonding surface of rapeseed petals, the maximum cracking force of the bonding surface of the rapeseed petals was smaller than that of the front side of the rapeseed petals in impending fracturing, the difference was significant, the maximum cracking force was small, which also indicated that the bonding surface of the rapeseed petals was easy to crack and shatter, and the two tests verified that the bonding surface of the rapeseed petals was more likely to shatter when impacted.
Under the same stress conditions, in the swing impact test, the shattering rate of the top rapeseed silique was twice that of the bottom siliques, and the silique at the bottom of the main branch was not easy to shatter, followed by the middle siliques, and at the top, siliques were most easy to shatter. The impending fracturing test also showed that the maximum cracking force and the bending strength of rapeseed siliques in different parts of the main branch differed significantly, but the impending fracturing test indicated that the maximum cracking force and bending strength of rapeseed siliques in the middle part of the main branch were not enough to shatter, siliques in the middle of the main branch were the highest, and the results differed from those of the swing impact test method, which might be due to the influence of the loading speed and the support position on the results of the test; the next step could be to carry out a correlation analysis and analyze the data to find out the reasons for the difference between the two methods.

Conclusions
(1) Based on the swing silique impact test method, according to the silique force state in the process of rapeseed harvesting and reeling, we carried out the silique crash impact test of the Ningza 1810, Zhenyou 8, and Fengyou 306 oilseed rape. Rapeseed varieties, the moisture content of siliques, silique growth position, collision materials, impact speed, force position, and other factors had a significant effect on silique shattering. Collision impact speed was lower than 1.5 m·s −1 , and the effect of various factors on silique shattering was not significant. Rapeseed siliques with a higher moisture content were not easy to shatter, and the front side of the silique petal was more resistant to external impacts than the bonding surface of the petals. The impact resistance of the front side of rapeseed was two to four times that of the bonding surface of rapeseed petals. Rapeseed siliques growing at the bottom of the plant were superior to those at the top and center of the plant in terms of resistance to external impacts, the shattering rate of the top rapeseed silique was twice that of the bottom siliques, and when siliques were supported, they were more likely to shatter under external forces compared to unsupported ones. (2) Based on the universal material testing machine, an experimental study on the mechanical properties of rapeseed siliques was carried out using the impending fracturing method, and the silique support position, force position, loading speed, and silique growth position had significant effects on the mechanical properties of siliques. The maximum cracking force and flexural strength were higher when the fruit body of the silique was supported, the range of the maximum cracking force was 3.05 N to 4.16 N, and the bending strength range was 8.48 MPa to 11.57 MPa. The maximum cracking force and flexural strength of silique frontal surfaces were higher than those of the bonding surfaces, and the rapeseed silique at the middle of the branch was stronger than those at the top and bottom of the branch in terms of the maximum cracking force and flexural strength. (3) The correlation between the morphological characteristics of silique and the characteristics of silique shattering was tested by using the swing impact test method, based on Pearson's correlation analysis and grey correlation analysis. The order of magnitude of the effect of morphological characteristics of rapeseed siliques on silique shattering performance was as follows: the angle of the silique stalk to the silique, diameter of stalk, thickness of pericarp, length of beak, thickness of silique, broad surface of silique, and length of silique.  Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the authors.