Determination of Strength Properties of Energy Plants on the Example of Miscanthus × Giganteus , Rosa Multiﬂora and Salix Viminalis

: Energy from biomass accounts for 70% of all renewables used for heat and electricity production. Such a signiﬁcant share of biomass determines the need for the investigation of their mechanical properties, as most of the lignocellulosic material requires cutting, chipping or milling before its utilization for energy purposes. Therefore, the knowledge about cutting resistance, bending sti ﬀ ness, and impact strength of the energy plants is very important. The values of these parameters are used in the proper selection of shredding machines and their elements, wrapping nets or determination of power demand during raw material conversion. This paper presents the results of research on the mechanical properties of selected energy plants. The scope of the research included three di ﬀ erent plant species: Miscanthus × giganteus , Rosa multiﬂora , and Salix viminalis , investigated in terms of cutting resistance, bending sti ﬀ ness and impact strength of stalks. The results showed that the average stalk cutting resistance for the rotation speed of 4200 RPM was 0.17 N · mm − 2 for the Miscanthus × giganteus , 0.15 N · mm − 2 for the Rosa multiﬂora and 0.2 N · mm − 2 for the Salix viminalis . Meanwhile, for a rotation speed of 3200 RPM, the cutting resistance amounted to 0.15 N · mm − 2 for Miscanthus × giganteus , 0.16 N · mm − 2 for Rosa multiﬂora and 0.18 N · mm − 2 for Salix viminalis . For the impact measurements, the Salix viminalis exceeded 40 J · mm − 2 of absorbed energy. Meanwhile, the average impact strength value for the Rosa multiﬂora was 0.53 J · mm − 2 and for the Miscanthus × giganteus was 0.22 J · mm − 2 . The bending sti ﬀ ness of Miscanthus × giganteus at an average modulus of 3.44 GPa was 1.1 N · m 2 for the basal zone, 0.78 N · m 2 for the central zone, and 0.72 N · m 2 of the apical zone. For the average Young’s modulus of 0.19 GPa, the bending sti ﬀ ness of the Rosa multiﬂora reached a value of 0.64 N · m 2 for the basal zone, 0.23 N · m 2 for the central zone, and 0.28 N · m 2 for the apical zone. The Salix viminalis , with an average modulus of elasticity of 0.23 GPa, achieved bending sti ﬀ ness in the basal zone of 0.99 N · m 2 , the central zone 0.33 N · m 2 , and the tip zone 0.38 N · m 2 . This research makes it possible to expand our knowledge in the ﬁeld of biomass processing and construction of agricultural machinery with higher processing e ﬃ ciency.


Introduction
According to the Global Bioenergy Statistics (2018), biomass accounts for 70% of the total primary energy supply of all renewables; moreover, it accounts for 95.9% of heat generation derived from renewables [1]. The use of biomass for energy purposes supports the reduction of carbon dioxide emissions and is in line with the strategy adopted by the European Union [2,3]. Biomass has already attracted attention from the fuel market in the context of controlling greenhouse gas emissions [4].
The European Union, with respect to environmental protection, encourages countries to introduce energy crops [15]. The harvested biomass coming from energy crops can be easily converted into a valuable form for the energy market. The most popular conversion methods are chipping, baling, briquetting or pelleting. The choice of method depends mainly on the market demand and needs of the final user/consumer. In the case of chipping, it is important to properly match both the cutting elements and the parameters of the chipping device. About 50% of energy is absorbed by forage harvester for shredding biomass material. The following plant parameters affect the grinding process: moisture, hardness, stiffness, compactness, structure, and coefficient of friction [16,17]. Additionally, the appropriate selection of biomass processing machines and methods can increase the efficiency of processing, reduce the risk of the failures, and reduce the consumption of energy and raw materials. Therefore, proper examination of the processed material allows the selection of the process, as well as the components and materials, that build the machine. For this purpose, it is necessary to determine the strength parameters of the grinding material [18][19][20][21][22]. Some of the most common methods for testing the strength properties of plants are bending stiffness, impact resistance and determination of destructive stresses, including shearing stresses [23].
Until now, the impact strength of energy plants has not been measured using the Charpy method, which is described in detail in the PN-EN ISO 148-1: 2017-02 standard [24]. Thanks to the Charpy test, it is possible to adjust the appropriate method of grinding biomass and possible elimination of some agricultural treatments that might not bring the expected results. Skubisz et al. [25], Woliński et al. [26], and Prasada [27] have measured stalk shear energy by the impact method using a knife blade as a tool that causes an impact in a Charpy's impact test stand. Moreover, the measurement method concerned individual samples and assumed the constant speed of the blade forced by its free fall. Furthermore, Lisowski et al. [28], Frączek et al. [29], Chen et al. [30] and Chattopadhyay et al. [31] studied the cutting energy of stalks of Sorghum, Hemps, Jerusalem artichokes, and Miscanthus × giganteus, accordingly. Yore et al. [32] determined the cutting energy for rice straw at a cutting speed of 500 strokes per minute. The authors performed the tests for one, two and five stems, and concluded that as the number of cut stalks (the cross-section) increases, the unit energy absorbed per cut decreases (concerning one stalk).
Igathinathane et al. [33] checked the force needed to cut the corn stalk depending on its orientation with the blade. Azadbakht [34] investigated the dependence of canola stalk cutting energy on moisture content and cutting height. However, the tests were performed only in quasi-static conditions. Nowakowski et al. [35] conducted research related to the modulus of elasticity of selected energetic plants, including Miscanthus × giganteus, Salix viminalis, and Rosa multiflora. However, no data concerning stiffness on the bending of the stalks, or the relationship between the measurement zones, diameter, and bending power during the bending test, were presented. Such measurements, for example, may allow the selection of a sufficiently durable net wrapping the plant bale. Kronbergs [22] examined the ultimate tensile, the shear strength, the modulus of elasticity, the shear modulus, and the energy demand for compacting and pressing of wheat stalks. Mathanker et al. [36] investigated the cutting energy of energy cane at different speeds and stalk orientation relative to the blade.
As a result, in the literature, there are limited data about the strength and cutting resistance at variable speeds of selected energy plants. Moreover, no information was found about energy crops with respect to their bending stiffness and the relations between the cross-sectional surface, the speed of quasi-static load, bending force and displacement. Therefore, for a more extensive plant survey that can be used to improve the current state of the agricultural technology, this work aimed to examine the resistance of selected energy plants to cut stalks at variable rotation speeds, as well as to determine their impact strength and rigidity for bending.

Materials Used in the Experiments
In the tests, the following annual biomass stalks from energy crops were used ( Figure 1): Miscanthus × giganteus, Rosa multiflora, and Salix viminalis. Salix viminalis is a perennial woody crop, Miscanthus × giganteus is a perennial herbaceous C4 crop, whereas Rosa multiflora of the Jatar cultivar is a perennial shrub that is usually harvested in two-year rotations. These plants are characterized by resistance to moderately cold climate and drought. Therefore, they are most common in the countries of Western-Northern Europe. In this climate, Rosa multiflora grows in the wild and has relatively low soil requirements. It tolerates drought and low temperatures, as well as being resistant to pathogens. Additionally, the high calorific value of this crop makes it suitable for utilization as an energy source in power plants. Miscanthus × giganteus and Salix viminalis have low soil requirements and are characterized by fast growth; therefore, they are used not only for energy purposes, but also as raw materials [37].
Energies 2019, 12, x; doi: FOR PEER REVIEW www.mdpi.com/journal/energies Nowakowski et al. [35] conducted research related to the modulus of elasticity of selected energetic plants, including Miscanthus x giganteus, Salix viminalis, and Rosa multiflora. However, no data concerning stiffness on the bending of the stalks, or the relationship between the measurement zones, diameter, and bending power during the bending test, were presented. Such measurements, for example, may allow the selection of a sufficiently durable net wrapping the plant bale. Kronbergs [22] examined the ultimate tensile, the shear strength, the modulus of elasticity, the shear modulus, and the energy demand for compacting and pressing of wheat stalks. Mathanker et al. [36] investigated the cutting energy of energy cane at different speeds and stalk orientation relative to the blade.
As a result, in the literature, there are limited data about the strength and cutting resistance at variable speeds of selected energy plants. Moreover, no information was found about energy crops with respect to their bending stiffness and the relations between the cross-sectional surface, the speed of quasi-static load, bending force and displacement. Therefore, for a more extensive plant survey that can be used to improve the current state of the agricultural technology, this work aimed to examine the resistance of selected energy plants to cut stalks at variable rotation speeds, as well as to determine their impact strength and rigidity for bending.

Materials Used in the Experiments
In the tests, the following annual biomass stalks from energy crops were used ( Figure 1): Miscanthus × giganteus, Rosa multiflora, and Salix viminalis. Salix viminalis is a perennial woody crop, Miscanthus × giganteus is a perennial herbaceous C4 crop, whereas Rosa multiflora of the Jatar cultivar is a perennial shrub that is usually harvested in two-year rotations. These plants are characterized by resistance to moderately cold climate and drought. Therefore, they are most common in the countries of Western-Northern Europe. In this climate, Rosa multiflora grows in the wild and has relatively low soil requirements. It tolerates drought and low temperatures, as well as being resistant to pathogens. Additionally, the high calorific value of this crop makes it suitable for utilization as an energy source in power plants. Miscanthus × giganteus and Salix viminalis have low soil requirements and are characterized by fast growth; therefore, they are used not only for energy purposes, but also as raw materials [37].  The biomass was obtained from the one-year cultivation of energy crops belonging to the Wroclaw University of Environmental and Life Sciences in March 2019.
The range of the sample diameters was as follows: 5.71-10.87 mm for Miscanthus × giganteus, 2.60-11.80 mm for Rosa multiflora and 3.10-13.25 mm for Salix viminalis. Due to the different strength properties of The range of the sample diameters was as follows: 5.71-10.87 mm for Miscanthus × giganteus, 2.60-11.80 mm for Rosa multiflora and 3.10-13.25 mm for Salix viminalis. Due to the different strength properties of individual plant parts, and the large differences in diameters along the plant, rigidity testing was performed after the division of the plant stalks into three measuring parts (zones): basal, middle and apical. The measuring zones were 0.3 m long. The basal zone was taken along the length from the cutting place to the designated measuring length, the middle zone was determined from the center of the plant after 0.15 meters towards the basal and apical zone, and the apical zone was taken from the top of the plant according to the determined measuring length.

Proximate Analysis of Biomass
To characterize the investigated biomass, proximate analysis was performed. The samples of a randomly collected portion of each biomass (1500 g) were collected, dried and ground in the laboratory mill LMN 100 (Testchem, Pszów, Poland) to the fraction below 1 mm. From the obtained pulverized material, the required samples (in accordance with the applied ISO Standards) for further analyzes were used, accordingly. The proximate analysis included parameters such as moisture content (MC), higher heating value (HHV), lower heating value (LHV), ash content (AC), and volatile matter content (VM).
The moisture content was determined according to PN-EN ISO 18134-2:2017-03E [38] using a laboratory moisture analyzer SARTORIUS MA150 (Sartorius, Goettingen, Germany). The higher heating value (HHV) was determined in a calorimetric bomb IKA C200 (IKA, Lucknow, India) in accordance with PN-EN ISO 18125:2017-07 [39]. Next, the lower heating value (LHV) was calculated as a function of the higher heating value (HHV) and moisture content in the biomass according to [40].
where HHV is the higher heating value (MJ·kg −1 ), MC is moisture content, and r is the latent heat of water vaporization (r = 2.

Procedure of the Cutting Resistance Test
The research device (Figure 2), for cutting resistance measurements, was designed and constructed in the Department of Fundamentals of Technology at the Wroclaw University of Environmental and Life Sciences. The cutting resistance tests were carried out for three different initial rotation speeds of the chopper, successively: 4200, 3200 and 2200 RPM. The power of the engine was 2.5 kW. A laser DT-2236 (Lutron Electronics Company, Coopersburg, Pennsylvania, United States) was used to measure the rotation speed with an accuracy of 1 RPM. Torque measurements were carried out using an 8652-6001 torque sensor (Burster GmbH & Co. KG, Gernsbach, Germany). The diameter of each sample was measured with a caliper DC-1 (Ovibell GmbH & Co. KG, Mulheim, Germany), and then the shear surface was determined. The cutting resistance was determined as follows: where F is the cutting force (N), τ is the torque during cutting, L is the position vector of force (66 mm).
where CR is the cutting resistance (N·mm −2 ), A is the cutting surface (mm 2 ). The cutting resistance was determined as follows: where F is the cutting force (N), τ is the torque during cutting, L is the position vector of force (66 mm).
where CR is the cutting resistance (N·mm −2 ), A is the cutting surface (mm 2 ).

Procedure of the Bending Stiffness Test
The bending test was performed according to the PN-EN 310:1994 standard [43]. Due to the differences in the diameter of individual parts of the material, the biomass stalks were divided into three measurement parts (zones): basal (1), middle (2) and apical (3), as shown in Figure 3. The cutting resistance was determined as follows: where F is the cutting force (N), τ is the torque during cutting, L is the position vector of force (66 mm).
where CR is the cutting resistance (N·mm −2 ), A is the cutting surface (mm 2 ).

Procedure of the Bending Stiffness Test
The bending test was performed according to the PN-EN 310:1994 standard [43]. Due to the differences in the diameter of individual parts of the material, the biomass stalks were divided into three measurement parts (zones): basal (1), middle (2) and apical (3), as shown in Figure 3. Before each measurement, the diameter of the sample was measured to determine the moment of inertia and the bending index. The prepared samples were subjected to a bending load using the Instron Before each measurement, the diameter of the sample was measured to determine the moment of inertia and the bending index. The prepared samples were subjected to a bending load using the Instron 5566 (Instron, Norwood, Massachusetts, United States) ( Figure 4) with a measuring range up to 1 kN. The measuring stand consisted of a head equipped with a sample loading attachment in its central point and a bottom plate, to which an element constituting two-point support for the material subjected to bending test was fixed. 5566 (Instron, Norwood, Massachusetts, United States) ( Figure 4) with a measuring range up to 1 kN. The measuring stand consisted of a head equipped with a sample loading attachment in its central point and a bottom plate, to which an element constituting two-point support for the material subjected to bending test was fixed. The load was controlled using different displacement speeds (1.8, 30, 100 and 500 mm·min −1 ) in five repetitions. Based on the analysis of the peak value of the bending force in the displacement function, the modulus of elasticity and breaking stress for stalks were determined using Formulas (4)-(8), accordingly: where: Mg-bending moment (N·m), P-force (N), L-spacing of supports (mm).  The load was controlled using different displacement speeds (1.8, 30, 100 and 500 mm·min −1 ) in five repetitions. Based on the analysis of the peak value of the bending force in the displacement function, the modulus of elasticity and breaking stress for stalks were determined using Formulas (4)-(8), accordingly: where: Mg-bending moment (N·m), P-force (N), L-spacing of supports (mm).
where: Ix-the moment of inertia (mm 4 ), d-diameter of the stalk (mm).

Procedure of the Impact Strength Test
The impact strength test was performed using the Charpy method. The measurement method was based on the PN-EN ISO 148-1: 2017-02 standard [24]. Due to the heterogeneity of investigated material and its unique dimensions, the biomass sample was adopted to the applied standard ( Figure 5), accordingly. The tests were carried out in five replications.
Potential energy absorption by biomass stalks was obtained from the Charpy impact test. The sample dimensions were determined using a caliper to calculate the fracture area. The following formula was used to determine the impact strength of the samples: where K is impact strength (J·cm −2 ), K 1 is energy absorbed by the sample (J) and A is fracture surface (cm 2 ).
Energies 2019, 12, x; doi: FOR PEER REVIEW www.mdpi.com/journal/energies its unique dimensions, the biomass sample was adopted to the applied standard ( Figure 5), accordingly.
The tests were carried out in five replications.  Potential energy absorption by biomass stalks was obtained from the Charpy impact test. The sample dimensions were determined using a caliper to calculate the fracture area. The following formula was used to determine the impact strength of the samples: where K is impact strength (J·cm −2 ), K1 is energy absorbed by the sample (J) and A is fracture surface (cm 2 ).

Biomass Characteristics
The results of proximate analysis (moisture and ash content, low heating value, and volatile matter) of investigated energy crops are presented in Table 2. The obtained values are very similar to the parameters found in the literature. Vassilev et al. [44] examined the ash content and volatile matter in Miscanthus × giganteus with a slightly higher result of ash content (4.80%) and VM of (79.00%). Lewandowski et al. [45] investigated the values of ash content in different seasons in five locations in Europe. Miscanthus × giganteus cultivated in Portugal, harvested in winter, had the most similar results in terms of ash content. Tursunov et al. [46] determined the same properties of Rosa multiflora (2.4% of ash content and 76.4% of VM). He also examined the LHV for Miscanthus × giganteus (16.9 MJ·kg −1 at the moisture content of 8.8%) and Rosa multiflora (17.6 MJ·kg −1 at moisture content 18.3%). In turn, Weger et al. [47] analyzed Salix viminalis, obtaining the following values: ash content (1.35%), VM (81.02%) and LHV (18.66 MJ·kg −1 ).
The strength properties of energy crops were measured for moisture content and the results are presented in Table 3, along with the division into plant zones. The moisture values are typical for investigated energy crops during harvesting process [10,11].  Tables 4 and 5 show the cutting resistances, strength, and power required for the shredding process at initial rotation speeds of 4200, 3200, 2200 RPM. Unfortunately, the tests revealed that 2200 RPM for the 2.5 kW power of the engine was too low to crush the material; therefore, the results are shown for an initial speed of 3200 and 4200 RPM only. Based on the analysis of Tables 4 and 5, it can be concluded that the highest cutting resistance was achieved by Salix viminalis, reaching values between 0.15 and 0.23 N·mm −2 , for which the cutting resistance increased at the moment of lowering the rotation speed ( Figure 6). In contrast to Salix viminalis, Miscanthus × giganteus and Rosa multiflora showed a reverse strength property, whose cutting resistances were 0.11-0.23 N·mm −2 and 0.12-0.17 N·mm −2 , respectively. In those cases, the cutting resistance to their shredding increased with the rotation speed. It can be explained by the low modulus of Salix viminalis with a relatively large diameter of stalks.   Chen et al. [30] studied the shear force in quasi-static conditions required for Hemp stalks, which amounted to 243 N. However, this value does not give complete information on the strength of the sample since the cross-section of the sample was not taken into account. Chattopadhyay et al. [31] examined the Sorghum stalks for the absorption of a part of the potential energy, thus determining the resistance in the form of energy falling on the shear surface. The results were 0.101-0.142 MJ·mm −2 , while Frączek et al. [29] investigated Topinambur, obtaining results between 0.014 and 0.023 MJ·mm −2 . However, both mentioned measurements concerned the energy consumption of the process (taking into account the cross-section of the sample absorbing the energy), and not the cutting resistance, as examined in this work. Therefore, the comparison of the results was not possible.

Cutting Resistance of the Biomass
Lisowski et al. [28] checked the resistance and cutting energy of maize stalks in internodes and nodes, obtaining values of 16.9-60.4 N·mm −2 and 10.93-26.82 J·mm −2 , depending on the measurement zone. These results showed a significant difference between the cutting resistance in the dynamic and static conditions. The cutting in the quasi-static conditions requires much more energy to cut the stalks compared to the dynamic conditions. Lisowski et al. indicated that the cutting resistance of corn stalks was one hundred times higher than the cutting resistance of Rosa multiflora or Miscanthus × giganteus at a rotation speed of 3200 RPM obtained in this study (0.15-0.19 N·mm −2 ). The cutting energy for rapeseed stalks and pods under dynamic conditions was investigated by Skubisz et al. [25]. The values of the cutting energy were 9.88·10 4 J·mm −2 for pods and 11.87·10 4 J·mm −2 for stalks. However, these values cannot be compared with the results of this study, because they take into account the energy only and not the force acting on the shear surface. Similar research was performed for buckwheat stalks by Woliński et al. [26], for which the cutting energy amounted to 0.33-0.39 J. Unfortunately, the authors did not include the shear surface of the investigated material, only the energy needed to cut the stem was considered, so the results of the test cannot be compared, as well.
Yore et al. [32] investigated the cutting energy for rice straw at the cutting speed of 500 strokes per minute. The results showed that as the number of stems (cutting surface) increases, the amount of energy needed to cut one stem decreases. The highest value of shearing energy for one rice straw stem was 0.22 J, while average shearing energy per stem was 0.15 J (for a five-stem test). The highest peak force was 8.3 N for one stem and 5 N per stem for normalized results. Unfortunately, as the energy and the shear force did not correspond to the cutting surface, the values cannot be compared with the results obtained in this study.
In turn, Igathinathane et al. [33] investigated the influence of the position of the corn stalks relative to the cutting blade. The results showed that the highest cutting force required the corn stalks to be oriented perpendicular to the blade (242 N), while for stalks set at an angle of 45 degrees, the force was 127 N, and for the parallel set only 46 N. In the case of these tests, the cutting cross-section was also not taken into account; moreover, the tests were carried out under quasi-static conditions, which prevented comparison of the results.
Based on the trend line (linear equation, in this case), it was observed that the increase in rotation speed causes a slight reduction in shearing forces (Figures 7 and 8). Such a slight difference means that the force needed for comminution is similar in both cases; however, a minimal increase in force when lowering the speed indicates that a further increase in force can be expected along with lowering the rotation speed, as it was detected using a speed of 2200, at which cutting resistance was not high enough to make it impossible to crumble the material. At a rotation speed of 4200 RPM, the force concerning the shear area increased more linearly. Additionally, the Miscanthus × giganteus at a rotation speed of 4200 RPM showed a similar relationship of force to the surface like for Rosa multiflora, while at a speed of 3200 RPM, Rosa multiflora showed a higher cutting resistance.       The increasing torque causes a reduction of the rotation speed under the load, which increases the shear force. The highest decrease in rotation speed was observed for Salix viminalis, characterized by the highest value of the shear resistance, while the lowest decrease was observed for Miscanthus × giganteus, characterized by the lowest shear resistance (Figures 9-11).
The increasing torque causes a reduction of the rotation speed under the load, which increases the shear force. The highest decrease in rotation speed was observed for Salix viminalis, characterized by the highest value of the shear resistance, while the lowest decrease was observed for Miscanthus × giganteus, characterized by the lowest shear resistance (Figure 9,10,11).        Additionally, the application of the statistical method ANOVA (Analysis of Variance) revealed that the surface (cross-section area) had the greatest influence on the cutting force for Miscanthus × giganteus and Rosa multiflora. In the case of Salix viminalis, a weak dependence of strength on the other parameters was observed. The details of this analysis (at the probability value p < 0.05) are presented in the Supplementary Materials (Table S1). In turn, for Miscanthus × giganteus, the analysis showed a dependence of cutting resistance on all invetigated parameters (Table S2 in Supplementary Materials).
Analysis of force variance and cutting resistance depending on the initial speed and decrease in rotation, presented in Supplementary Materials (Tables S3 and S4) showed the relationship between force and rotation speeds for Miscanthus × giganteus and Salix viminalis. However, the relationship between cutting resistance and rotation speeds was observed only for Miscanthus × giganteus. There was no correlation between the considered parameters for Salix viminalis.
The studies showed also that for a greater rotation speed (4200 RPM in this case) the power demand for the investigated energy plants chipping increased in relation to the lower rotation speed (3200 RPM in this case). The power demand can be calculated using the following formula: where: P is a power (W), n is a rotation speed (RPM), and M is a torque (N·m). For a rotation speed of 4200 RPM, the power needed to chop the material was 802 W for Miscanthus × giganteus, 806 W for Rosa multiflora and 896 W for Salix viminalis, while at 3200 RPM the power demand for Miscanthus × giganteus was 438 W for Rosa multiflora 482 W and 467 W for Salix viminalis. The differences between the power demanded for branch chipping at the given rotation speeds (3200 RPM and 4200 RPM) differed by 48% for Salix viminalis, 45% for Miscanthus × giganteus and 40% for Rosa multiflora, respectively. Thus, the rotation speed indirectly influences energy consumption during biomass processing (cutting, chipping, shredding). As a result, the rotation speed significantly affects also the energy efficiency of those processes, and its selection is an indispensable parameter that should not be overlooked during the selection of the equipment and settings of the operation conditions. During the test, it was observed that for Miscanthus × giganteus and Rosa multiflora, the cutting resistance was increasing along with rotation speed, in contrast to Salix viminalis. This relationship is due to Salix viminalis requiring higher blade rotation speeds because it is resistant to impact forces and its comminution is only due to the cutting forces. It is important information from a biomass processing point of view, indicating that the cutting efficiency, the energy consumption and the equipment used at a given rotation speed may not be effective for every type of energy plant.

Impact Strength Test of the Biomass
Based on the results of the tests carried out, it was found that the highest impact strength value was observed for Salix viminalis. The energy needed to break it went beyond the range of the measuring instrument and amounted to above 0.57 J·mm −2 . Tables 6 and 7 show the results of the Charpy impact test for Miscanthus × giganteus and Rosa multiflora. Based on the results, it can be concluded that Rosa multiflora has a higher impact strength than Miscanthus × giganteus. For Rosa multiflora, the energy absorbed exceeded twice the measuring range of the stand. Table 6. Impact strength of Miscanthus × giganteus. The differences in measurements were significant, and resulted from both the surface of the sample cross-section and the sampling zone, for the Miscanthus × giganteus standard deviation was 0.15 J·mm −2 and for the Rosa multiflora 0.06 J·mm −2 . The highest value of Rosa multiflora impact strength was 0.77 J·mm −2 , with an average impact strength 0.53 J·mm −2 , while the lowest value was 0.32 J·mm −2 . For Miscanthus × giganteus, the highest value amounted 0.30 J·mm −2 and the lowest one was 0.13 J·mm −2 . Moreover, applying ANOVA, a high dependence of impact strength (at the p < 0.05) on the surface was found for Miscanthus × giganteus, as shown in the Supplementary Materials (Table S5). Table 7. Impact strength of Rosa multiflora. The impact tests carried out led to the conclusion that Miscanthus × giganteus belongs to the group of plants characterized by a low impact strength, which indicates its fragility. Therefore, it is more suitable for impact crushing than cutting. The difference in impact strength in the tested plants proves that the grinding methods should be properly selected, depending on their mechanical properties. Grinding a plant with a low impact strength (Miscanthus × giganteus) using knives may be unnecessary and economically unprofitable, as its apparent cutting may result from its crushing. Then, it is unnecessary to use and constantly replace expensive blades, because a better solution would be to use a hammer mill.

Bending Stiffness of the Biomass
The dependence of Young's modulus, bending stress and bending stiffness on the stalk measurement zone are shown in Figures 12-14. In the case of Young's modulus (Figure 12), significantly higher values were observed for Miscanthus × giganteus. The greatest value was obtained in the middle zone of the stalk (4.42 GPa). In the basal and apical zones, Young's modulus amounted to 2.93 GPa and 2.97 GPa, respectively. The maximum achieved in the middle zone can be explained by the fact that Miscanthus × giganteus stalk is thick, but hollow in the basal zone, which weakens it in terms of the strength of the material against bending. The middle zone is filled with tissue in the whole volume and relatively thicker and stronger than the apical zone. In the case of Salix viminalis and Rosa multiflora, Young's modulus was much lower. Depending on the zone, for Salix viminalis, it was 0.13-0.29 GPa, and for Rosa multiflora, it was 0.15-0.22 GPa. This arises from the fact, that both plants have morphological features completely different from Miscanthus × giganteus. Salix viminalis and Rosa multiflora are flexible and contain much more reinforcing tissue cells responsible for their elasticity.    Differences were found in bending stresses between the zones. Rosa multiflora and Salix viminalis achieved the highest values of bending stress in the basal zone and the lowest in the apical zone. Meanwhile, Miscanthus × giganteus achieved the highest value of bending stresses in the central zone. The reason for this was that the central zone was characterized by the highest strength properties, because it was lignified enough to have a high concentration of reinforcing tissue, while it was not ripe enough to develop a hole in the stalk, weakening the stalk to forces acting along its cross-section. Differences were found in bending stresses between the zones. Rosa multiflora and Salix viminalis achieved the highest values of bending stress in the basal zone and the lowest in the apical zone. Meanwhile, Miscanthus × giganteus achieved the highest value of bending stresses in the central zone. The reason for this was that the central zone was characterized by the highest strength properties, because it was lignified enough to have a high concentration of reinforcing tissue, while it was not ripe enough to develop a hole in the stalk, weakening the stalk to forces acting along its cross-section.
Finally, the bending stiffness in three zones of the energy plants was investigated ( Figure 14). The highest bending stiffness was achieved by Miscanthus × giganteus. The values amounted as follow: 1.1 N·m 2 for the basal zone, 0.78 N·m 2 in the central zone and 0.72 N·m 2 for the apical zone.   The lowest value of bending stiffness was recorded for Rosa multiflora, for which the highest result was achieved by the basal zone (0.64 N·m 2 ) and the lowest by middle zone 0.23 N·m 2 , in the apical zone the value of the stiffness was 0.28 N·m 2 . Salix viminalis achieved the following results: for the basal zone, 0.99 N·m 2 ; middle zone, 0.33 N·m 2 ; and apical zone 0.38 N·m 2 . The basal zone was characterized by the greatest stiffness for all investigated energy plants. As the plant grows, the thickness of the stalk increases, reducing its flexibility, as well as its susceptibility during the baling process. On the other hand, too low a stiffness of the plant may have a negative impact on the cutting process; however, tests of cutting resistance have determined that at given stiffness the middle and apical zone of plants are susceptible to crushing by the knife method.
During the cutting resistance tests, the Miscanthus × giganteus was shredded due to the impact strength of the blade than because of the cutting force. Such a behavior of the investigated biomass is confirmed by the results of the bending stiffness, stress and Young's modulus tests. Miscanthus × giganteus was the most fragile plant characterized by the highest bending stiffness and much higher Young's modulus in relation to Rosa multiflora and Salix viminalis.
Additionally, the initial correlation between physical parameters was studied. Based on graphical maps in Figure 15, it can be concluded that the bending stiffness increased with the sample diameter. The speed of sample loading did not significantly affect the results of bending stiffness and bending stress. The apical zone was characterized by smaller stresses during sample loading, which was related to the bending rigidity decreasing along with the diameter.
Additionally, the initial correlation between physical parameters was studied. Based on graphical maps in Figure 15, it can be concluded that the bending stiffness increased with the sample diameter. The speed of sample loading did not significantly affect the results of bending stiffness and bending stress. The apical zone was characterized by smaller stresses during sample loading, which was related to the bending rigidity decreasing along with the diameter. The analysis of variance was also performed in the case of the tests of the bending stiffness (Table S6  and Table S7 in the supplementary materials). The aim was to determine the influence of the stalks' diameter, measuring zone and speed of sample loading on the bending stiffness of the investigated energy plants (stalks). According to the results shown in table S7, a significant dependence of bending stiffness on the diameter was observed for Rosa multiflora and Salix viminalis. For Miscanthus × giganteus, this correlation was not as relevant. The bending stiffness is also affected by speed and measuring zone (Table S7), but these parameters do not have an influence as great as that of the diameter of the plant (Table S6). The results seem to be satisfactory, but they need further verification, including tests with the use of other plants. The analysis of variance was also performed in the case of the tests of the bending stiffness (Tables S6 and S7 in the Supplementary Materials). The aim was to determine the influence of the stalks' diameter, measuring zone and speed of sample loading on the bending stiffness of the investigated energy plants (stalks). According to the results shown in Table S7, a significant dependence of bending stiffness on the diameter was observed for Rosa multiflora and Salix viminalis. For Miscanthus × giganteus, this correlation was not as relevant. The bending stiffness is also affected by speed and measuring zone (Table S7), but these parameters do not have an influence as great as that of the diameter of the plant (Table S6). The results seem to be satisfactory, but they need further verification, including tests with the use of other plants.

Conclusions
Knowledge regarding biomass mechanical properties is very important, as its potential utilization for energy purposes is still growing, especially in terms of energy crops and waste raw materials. The strength parameters influence machinery selection and other materials used during biomass processing (especially harvesting and shredding). The understanding of the biomass behavior contributes to the reduction of energy consumption, operation time and limitation of unexpected failures, and thus renewable fuel production costs.
The studies showed that energy plants have different mechanical properties. As a result, to optimize the process of further processing, it is necessary to use other technologies or machines. The mechanical tests with the use of Miscanthus × giganteus, Salix viminalis, and Rosa multiflora provided data supporting the decision making in proper planning of the harvested biomass treatment.
In terms of the processes' efficiency, it seems to be necessary to continue further measurements of the cutting efficiency at given rotation speeds using sieve analysis and to test the energy consumption by machinery using both methods of impact: grinding and cutting. Moreover, the expanding data regarding the strength properties of plants and the impact of various loads can be obtained using the finite element method on the designed discrete model of the plant. More studies are also required to validate the correlation model concerning the mechanical fuel properties, operation conditions of the cutting machinery and size of the energy plants.
Supplementary Materials: The following are available online at http://www.mdpi.com/1996-1073/12/19/3660/s1, Table S1: Correlation analysis for force dependent variable of surface and rotation speeds. Table S2: Correlation analysis for cutting resistance dependent variable of surface and rotation speeds. Table S3: Analysis of the force dependent variable from rotation speed and rotation decrease. Table S4: Analysis of the cutting resistance dependent variable from initial rotation speed and rotation decrease. Table S5: Multivariate analysis of variance for the impact strength dependent variable from the surface independent variable. Table S6: Correlation between stiffness, diameter and measuring speed for each of the plants concerned. Table S7: Bending stiffness variance analysis for speed and zone.
Author Contributions: M.S. and R.S. conceived, designed and performed the research; M.S. and A.D. analyzed the data, contributed materials analysis tools and wrote the paper.
Funding: This research received no external funding. The APC was funded by the Institute of Agricultural Engineering of the Wroclaw University of Environmental and Life Sciences (Statutory activity B010/0036/19).

Conflicts of Interest:
The author declares no conflict of interest.