Skip to main content

Measuring the compressive modulus of elasticity of pith-filled plant stems



The compressional modulus of elasticity is an important mechanical property for understanding stalk lodging, but this property is rarely available for thin-walled plant stems such as maize and sorghum because excised tissue samples from these plants are highly susceptible to buckling. The purpose of this study was to develop a testing protocol that provides accurate and reliable measurements of the compressive modulus of elasticity of the rind of pith-filled plant stems. The general approach was to relying upon standard methods and practices as much as possible, while developing new techniques as necessary.


Two methods were developed for measuring the compressional modulus of elasticity of pith-filled node–node specimens. Both methods had an average repeatability of ± 4%. The use of natural plant morphology and architecture was used to avoid buckling failure. Both methods relied up on spherical compression platens to accommodate inaccuracies in sample preparation. The effect of sample position within the test fixture was quantified to ensure that sample placement did not introduce systematic errors.


Reliable measurements of the compressive modulus of elasticity of pith-filled plant stems can be performed using the testing protocols presented in this study. Recommendations for future studies were also provided.


The measurement of mechanical properties of plant stems helps in investigating early and late season stalk lodging [1]. But in spite of the economic significance of plants with thin-walled stems (e.g., maize, sorghum, wheat, etc.), few studies have investigated reliable methods for obtaining their mechanical properties under compressive loading. One of the most important mechanical properties is the modulus of elasticity, which provides a linear relation between stress and strain [2]. This mechanical property is essential for calculating stress states as well as physical deformation of a structure or a plant [3, 4].

The modulus of elasticity can be measured in a number of ways, including bending, tension, compression, vibration, and acoustic excitation tests. Bending tests have been used in a number of studies, including those focused on the mechanical properties of wood [5,6,7,8], sunflower stalks [9], sorghum stalks [10], wheat stems [11], and maize stalks [12, 13]. Bending tests are popular because they involve low loads, easily measurable deformation, and require little sample preparation. Bending tests can only be performed on test samples that are long and slender [14], and produce one estimate of the modulus of elasticity for each sample. As a result, this method produces rather poor spatial resolution for the modulus of elasticity. The accuracy of the modulus obtained by bending tests is also adversely affected by the nonlinear form of the bending equations, which tends to amplify measurement uncertainty.

Tensile testing is another common technique for obtaining the modulus of elasticity. This approach has been used to measure the modulus of elasticity of wood [7], excised rind sections of maize stems [15, 16], excised longitudinal sections of switchgrass stems [17], rice stems, and Arabidopsis stems [18]. However, sample preparation is more laborious as compared to bending and specimens must be gripped securely without inducing tissue damage. The gripping aspect of tensile test is often quite challenging.

Compression testing is very common in the wood literature [7, 19], but is not commonly used in the testing of thin-walled plant stems. This is because the plant rind tends to be highly susceptible to buckling deformation. Consequently, information on the compressive modulus of thin-walled plant stems is often not available.

Studies have reported that the tensile and compressive modulus of elasticity values can be different for lumber, wheat straw, and barley straw [6, 20]. This indicates that tensile testing alone may be insufficient for measuring the modulus of elasticity, and that bending tests (which induce both bending and compression) may yield modulus values which are unreliable. Techniques for measuring the compressive modulus of elasticity of plant stems are therefore needed.

For thin-walled plant stems, bending loads are primarily borne by longitudinal stresses in the rind tissue [21, 22], so this study focuses on the longitudinal modulus of elasticity.

The goal of this study was to develop a robust method for obtaining the compressive modulus of elasticity of the rind of pith-filled plant stems, and to study the factors that influence the accuracy and reliability of this method. For the sake of brevity, the abbreviated term “modulus of elasticity” will be used in place of the more precise term “longitudinal compressive modulus of elasticity” in the remainder of this paper.


Stalk samples

Dry maize stalks were used as test specimens in this study. Maize can be highly susceptible to late-season stalk lodging, which occurs due to compression-induced buckling of the rind [13]. Maize stalks were sampled from 2 replicates of four commercially available hybrids of dent corn (maize) seeded at 5 planting densities (119,000, 104,000, 89,000, 74,000, and 59,000 plants ha−1) [23]. Stalks were cut just above the ground and just above the ear node immediately before harvest. To prevent fungal growth, stalks were placed in forced-air dryers to reduce stalk moisture to approximately 10–15% moisture by weight, which closely mimics the state of stalks in the field just prior to harvest. To avoid confounding factors, only stalks found to be free of disease and pest damage were included in the study. One hundred (100) samples were selected for compression testing.

CT scanning

X-ray computed tomography was used to quantify cross-sectional areas of the rind and pith regions (Fig. 1). Stalks were scanned using an X5000 scanner (NorthStar Imaging, Rogers, MN, USA). The scanning process produced 2D cross-sectional images of the maize stalks. A customized computer program was used to extract the cross-sectional area of each stalk from the CT data. The scanning and morphology extraction are described in more detail in a previous study [23].

Fig. 1
figure 1

Transverse cross-section of a maize stalk as obtained by X-ray computed tomography. a X-ray CT image, b X-ray CT image overlaid with lines used to segment the image into rind and pith regions. Segmentation was performed with a custom computer algorithm [23]

Sample preparation

Technical standards have been developed for compression tests on metals [24], plastics [25], and biomaterials [26]. Each of these standards specifies that sample geometry is critical for accurate assessment of compressional stiffness. These standards require that samples are prepared with end faces that are planar and perpendicular to the loading axis (Fig. 2a). This insures that stresses applied during testing are evenly distributed throughout the specimen.

Fig. 2
figure 2

Compressive testing setup: a schematic diagram depicting geometric features of an ideal compression test; b a photograph of one specimen situated for testing

In this study, test specimens were cut from stalks with an abrasive saw (Bosch GCO2000, Gerlingen, Germany). The face of the rotating saw blade insured that end faces were planar. Because the rind is thickest just below the node line [13], specimens were cut just below each node, as shown in Fig. 3. This approach utilizes the natural architecture of the stalk to minimize the stresses applied to each end during testing. Specimens prepared in this manner were found to be very durable, thus enabling the performance of multiple tests on individual specimens without induced permanent tissue damage. The prepared specimens contained three distinct tissue regions, each of which differed in anatomy and geometry; internode tissue, elongation zone tissue, and a subapical primary elongation meristem region [27].

Fig. 3
figure 3

(Top) Maize rind thickness as a function of axial distance. (Bottom) X-ray computed tomography image of the corresponding maize stalk. Dashed lines indicate the locations where the rind is thickest. Samples were prepared by cutting near the dashed lines. A prepared sample is shown in Fig. 2b

Self-aligning compression platens are used in situations where the perpendicularity of sample end-faces is difficult to achieve. As a load is applied, these platens rotate until they are in alignment with the testing surface, thus accommodating any discrepancies in the angle of the end-face. Self-aligning platens (Cat No: S5722A, Instron Corp., Norwood, MA, USA) were therefore used at both ends of the specimens to accommodate any angular inaccuracies in the cutting process. Figure 2 provides a diagram and a photograph of a specimen situated for testing.

Compression testing equipment

Compression tests were performed using a universal testing machine (Instron 5965, Instron Corp., Norwood, MA, USA). Loads were measured with a 5 kN Instron load cell. Instrumentation control and data acquisition were managed with Instron software (Bluehill 3.0).

Two types of strain were measured for each sample in this study; overall strain (ε overall ) and local strain (ε local ). Overall strain was based on the total displacement of the universal testing machine, (i.e., the displacement between the two spherical platens) divided by the total initial length of the sample prior to loading. Local strain was measured using an Instron extensometer, which recorded the displacement of two points on the surface of the specimen (see Fig. 2). The extensometer had a reference length of 50 mm (Instron 2630 Series Dynamic Extensometer, Instron Corp., Norwood, MA, USA) (Fig. 2).

Compression testing procedure

When testing biological tissues, a preload and repeated application of load cycles is commonly used to bring the samples to a repeatable reference state [28]. This procedure is used to reduce measurement variability and is referred to as pre-conditioning [28,29,30,31]. The loading process is described below.

An initial load of 200 N was applied to each specimen. Five loading cycles were then applied. In each loading cycle, the load increased from 200 to 700 N and then returned to the 200 N initial state. The first cycle was used as a conditioning cycle. Only measurements from the latter four cycles were employed in the modulus of elasticity calculations. A strain rate of 0.1 mm/s and a sampling frequency of 33 Hz were used in this study. This rate is similar to that used in a previous report (0.0833 mm/s), where corn stalk specimens with a length to diameter ratio of 1:1 were tested [20]. Lower rates have been used in testing wheat/barley straw (0.04 mm/s) [20], lumber (0.005 mm/s) [7] and timber (0.042 mm/s) specimens [32]. Further investigation is needed in the future to elucidate the effect of strain rate on the compressive elastic moduli values of pith-filled plant stems.

Modulus of elasticity calculations

Compressive modulus is defined as the slope of a uniaxial stress–strain curve. Because the rind is the primary load-bearing tissue of the maize stalk [33], the compressional stress, σ, was obtained by dividing the applied force, F, by the cross-sectional area of the rind, A r (Eq. 1). Cross-sectional areas were measured 5 cm below the node.

$$\sigma = \frac{F}{{A_{r} }}$$

This approach neglects the structural contribution of the pith tissue, but allows the estimation of the rind stiffness from a single test. As will be shown in the results section, this assumption introduces relatively minor errors.

For small deformations, strain is obtained by dividing the change in length by the original length:

$$\varepsilon = \frac{\Delta L}{{L_{0} }} = \left( {\frac{{L_{f} - L_{0} }}{{L_{0} }}} \right)$$

The slope of the stress–strain curve, or compressive modulus, E, was calculated as follows:

$$E = \frac{\Delta \sigma }{\Delta \varepsilon } = \frac{{\sigma_{2} - \sigma_{1} }}{{\varepsilon_{2} - \varepsilon_{1} }} = \frac{{\left( {\frac{{F_{2} - F_{1} }}{{A_{r} }}} \right)}}{{\left( {\frac{{L_{2} - L_{1} }}{{L_{0} }}} \right)}} = \left( {\frac{\Delta F}{\Delta L}} \right) \frac{{L_{0} }}{{A_{r} }}$$

In this equation, ΔF and ΔL values in this study corresponded to the changes measured between F1 = 200 N and F2 = 700 N. Equation 3 was used to compute overall and local compressive modulus of each sample.

The above equations represent a standard approach to measuring the compressive modulus. Although the self-aligning platens accommodated non-perpendicular end-faces, the self-aligning nature of these platens in combination with the complex geometry of the stalk was found to cause circumferential variation in strain distribution within the specimens. To account for potential variations in stress, local strain measurements were obtained from 4 equally spaced positions around the circumference of each specimen, denoted by the angular position of each measurement: ε 0, ε 90, ε 180, and ε 270 (see Fig. 4). Equation 3 was used to calculate corresponding Compressive Modulus values (E 0, E 90, E 180, and E 270).

Fig. 4
figure 4

Top view of the self-aligning platen with a maize stalk cross-section and angular directions of strain measurement

These individual strain values were combined to obtain a more accurate assessment of the compressive modulus. Because strain is inversely related to the compressive modulus, special attention must be paid to the manner in which averaging is performed [34,35,36,37]. Instead of inserting individual strain values into Eq. 3, the strain values were first averaged to obtain a single average strain value (εlocal) representing the average cross-sectional strain:

$$\varepsilon_{local} = \frac{1}{4}\left( {\varepsilon_{0} + \varepsilon_{90} + \varepsilon_{180} + \varepsilon_{270} } \right)$$

This strain can be substituted into Eq. 3 as follows to obtain a final expression for the local compressive modulus, E local :

$$E_{local} = \frac{\Delta \sigma }{{\Delta \varepsilon_{local} }} = \frac{{\left( {\frac{{F_{2} - F_{1} }}{{A_{r} }}} \right)}}{{\frac{1}{4}\left\{ {\left( {\frac{{L_{2} - L_{1} }}{{L_{0} }}} \right)_{0} + \left( {\frac{{L_{2} - L_{1} }}{{L_{0} }}} \right)_{90} + \left( {\frac{{L_{2} - L_{1} }}{{L_{0} }}} \right)_{180} + \left( {\frac{{L_{2} - L_{1} }}{{L_{0} }}} \right)_{270} } \right\}}}$$

As before, the subscript indices 1 and 2 refer to the test conditions at loads of 200 and 700 N, respectively.

Assessing the contribution of pith tissue

After all specimens were tested, the contribution of pith tissue to overall stiffness was assessed by carefully drilling a hole of 5 mm in diameter through the nodal tissue at the end-face of each specimen. A common wood drill bit was used for this purpose. A round wood file was then used to gently abrade the pith tissue until only rind tissue remained. These hollow samples were then re-tested using the techniques described above.

Sensitivity of the compressive modulus to sample placement

The cross-sectional shape of the maize stalk is somewhat irregular (see Fig. 4). Placement of each specimen on the two self-aligning platens is therefore somewhat subjective. The sensitivity of the compressive modulus measurements to specimen placement was therefore assessed to determine if sample placement affected compressive modulus results.

These tests were performed by first placing a specimen at the apparent center of each self-aligning platen. After measuring the compressive modulus in the typical fashion, the specimen was shifted away from the center and the test was repeated. This process was repeated for shift distances of 2 mm and 4 shift directions (0°, 90°, 180°, and 270°). The compressive modulus was therefore measured at each of the 12 resulting shift locations. These measurements were balanced by 12 tests performed with the specimen in the center position. Testing alternated between centered and shifted positions to avoid potential bias caused by temporal effects.

Statistical analysis

Measurement repeatability

Repeatability of the compression test methodologies in this paper was performed according to standard procedures [38]. A set of 10 specimens were tested repeatedly according to the protocols described above. Each specimen was tested 5 times, and both methods for obtaining compressive modulus were used for each test. The standard deviation was used to quantify the test repeatability for each specimen.


Representative stress–strain curves of pith-filled maize samples

The stress–strain curves for both overall and local compressive modulus values were linear in nature. The loads used in this study typically resulted in strain values less than 0.5%. Representative curves are shown in Fig. 5, which illustrates that the stiffness measured via local strain measurements was generally higher than the overall stiffness. For all tests in this study, the coefficient of determination (R 2) between stress and strain was above 0.99.

Fig. 5
figure 5

Representative stress–strain curves for local and overall measurements. Slopes of each curve represent the respective compressive modulus values, E local , and E overall

Sensitivity of the compressive moduli to pith-filled sample placement

Preliminary testing revealed that compressive modulus values were sensitive to sample placement—but only when the sample was shifted more than 2 mm from the center of the platens. The authors’ experience in performing these tests is that a shift of more than 1 mm from the center is easily detectable to the human eye. To assess the influence of spatial position, 10 specimens were tested at the centered position and 2 mm shift positions. Each specimen was tested a total of 8 times: 4 times at the center location, and 4 tests with a 2 mm shift. Each of the 4 “shifted” tests involved shifting the specimen in a different direction, as shown in Fig. 4. The resulting data is shown in Fig. 6, which demonstrates that sample placement within ± 2 mm of the platen center had no significant effect on compressive modulus measurements.

Fig. 6
figure 6

Box plots illustrating the effect of spatial position on specimen placement. All data for a single specimen was normalized by the mean modulus value from the tests conducted at the centered position

Repeatability analysis for pith-filled samples

Recall that the test repeatability (i.e., test-to-test variation of a single specimen) was quantified by using the standard deviation for each of 10 samples. Both compression tests methods (using local strain or overall strain) were found to have a mean repeatability of 3.9%. Additional repeatability information is summarized in Table 1. The final column of Table 1 provides an upper bound on test-to-test variation at the 95% confidence level.

Table 1 Repeatability statistics obtained from repeated tests on a set of 10 specimens

Averaging local strains of pith-filled samples

Self-aligning platens induced slight circumferential variation in strain. This variation was captured by taking strain measurements at 4 circumferential locations for each specimen. To examine the effect of averaging circumferential strains, compressive modulus values were calculated by using 1, 2, 3, and 4 strain values. Equation 3 describes the calculation of compressive modulus for a single strain measurement while Eq. 4 describes the calculation process for four strain measurements. Similar expressions can be obtained for two and three strain measurements.

The effect of strain averaging is shown in Fig. 7. As expected, variation in the calculated compressive modulus decreased as the number of utilized strain measurement increased. This trend was evident at both the individual and group levels.

Fig. 7
figure 7

The effect of averaging local strain measurements around the circumference of the test specimens from 1, 2, 3, or 4 sides. Measuring strain from all sides reduces circumferential variation caused by structural asymmetry of the test specimens. Data in this chart is from 94 specimens, with 4 strain measurements per specimen (6 specimens were damaged during testing and therefore were excluded). Sample sizes reflect the number of different combinations for averaging strain measurements (e.g., given 4 strain measurements per specimen, there are 6 possible combinations when using groups of two, 4 combinations when using groups of three, etc.)

Neglecting the contribution of pith tissue

Several specimens were damaged either during the pith removal process or during subsequent testing, thus reducing the sample size for this portion of the study. The contribution of the pith tissue had a statistically significant effect on the compressive modulus values, which was found to be of approximately the same magnitude as specimen repeatability. The overall mean reduction in stiffness after pith removal was found to be approximately 4%. The variation in this effect was relatively high, which was likely due to the imprecise nature of the pith removal process. At a 95% confidence level, the mean effect of pith tissue on stiffness was found to be less than or equal to 6.3%. We therefore concluded that although the process of neglecting the pith tissue does introduce a consistent error, the magnitude of this error is not substantial. Table 2 provides a summary of the statistics related to pith removal. It is worth noting that the pith prevents failure due to buckling and may therefore significantly contribute to overall stalk strength but not stiffness [39].

Table 2 Statistical effects of pith removal

Local versus overall compressive moduli of pith-filled samples

We now examine differences between local and overall compressive moduli (E overall and E local ). These values were calculated for all specimens in this study. Recall that E overall is the stiffness of the entire specimen; whereas E local is the stiffness obtained near the center of each specimen (see Fig. 2).

Figure 8 provides distribution plots for all specimens tested in this study. The mean and standard deviation values for E overall and E local were (10.1 ± 1.5 and 12.8 ± 1.5 GPa, respectively). Figure 8 also provides distributions which were shifted downward by 4% to account for the effect of neglecting the pith tissue. Finally, Fig. 8 also provides comparisons to published data on the distribution of the compressive modulus values for dried wood from angiosperms and gymnosperms.

Fig. 8
figure 8

Overall versus local compressive modulus distributions for maize and the two major types of wood. The narrower, gray boxes indicate modulus that have been decreased by the average pith effect of 4%. Wood data from [40, 41]


This study involved the compressional testing of dry, non-diseased maize stalk segments consisting of two nodes and the intervening internode (see Fig. 1). Specimens were cut just below the node lines because tough nodal tissues and thicker rind in this region effectively distribute stresses, thus preventing premature tissue failures that can occur during compression testing when test specimens involve only internode tissues.

Certain challenges were encountered in this study. One of these was the difficulty of cutting two parallel end faces on maize stalk specimens, both of which should (according to compression testing standards) be perpendicular to the stalk axis. This challenge was addressed by using two self-aligning compression platens. However, this solution then generated a new challenge: the lack of structural symmetry induced circumferential variation in strain, thus necessitating the measurement of strain at multiple locations. These strains were averaged to obtain the compressive modulus of internodal tissues.

Accuracy, reliability, and test duration

Two different compressive modulus values were obtained for each specimen in this study: E overall and E local . The overall compressive modulus value is based on deformation that occurs throughout the entire specimen, including at the end faces, meristematic tissue, and internodal tissues. As such, the overall compressive modulus should be considered as an aggregate stiffness value, with tissue stiffness within the specimen varying above and below this value. The local modulus approach measures tissue strain in a region where tissue is regular and uniform and thus is likely more accurate. Deformation of the testing apparatus was negligible as compared to deformation in test specimens. The repeatability values of both tests were comparable.

The local compressive modulus values were higher than overall modulus values for every specimen in this study. Although spatial variation in stiffness was not the focus of this study, we believe that this is due to a lower tissue stiffness near each node and in the meristematic region [42]. More detailed studies will be necessary to confirm this. The calculation of the compressive modulus values was based on an assumption that the pith tissues have a negligible effect on stalk stiffness. The removal of pith tissue was found to decrease modulus values by an average of 4%. Thus, the values obtained for rind stiffness in this study are (on average) 4% higher than their true values.

As shown in Fig. 7, the reliability of local compressive modulus values improved as the number of circumferential sample points increased. However, unless multiple circumferential samples can be acquired simultaneously, each circumferential sample point increases the testing duration. Excluding sample preparation (which was similar for both test types), local modulus testing required approximately 10 min per specimen.

As shown in Fig. 7, the mean value is relatively insensitive to the number of circumferential strain values. However, a decrease in circumferential measurements also decreases test-to-test repeatability, thus artificially increasing the observed variation in the compressive modulus. The use of fewer circumferential measurements may be suitable in certain situations where the mean value is the primary objective.

If relative differences between plants is of primary concern, absolute accuracy may not be a primary concern. In such a case, the overall modulus may be a better choice. The overall modulus provides a single, average rind stiffness value for the entire specimen with reasonable reliability. Test duration for prepared samples was approximately 2 min per specimen.

Recommendations for future studies

One of the most important considerations when performing compression tests is the perpendicularity of end faces. This is a particular challenge when dealing with plant stems, which typically do not have straight edges that can be used as a reference. Spherical platens can be used to address this issue, and are recommended for future studies. If for some reason, spherical platens cannot be used, special attention should be paid to the preparation of end faces as well as the resulting load/deformation curves. An alternative approach is to embed each end of the sample in polymethyl methacrylate (PMMA) or some other kind of resin, a technique used in the testing of bone specimens [43, 44].

In the current study, rind thicknesses were obtained from X-ray computed tomography 2D images, but this approach requires special equipment and software. A more accessible technique is to obtain areas of rind and pith areas based on cross-sectional images obtained with a flatbed scanner [45].


Two methods were developed for measuring the compressive elastic modulus values of the rind of pith-filled plant stems such as maize. The two elastic modulus values were calculated using two different strain measurements. These test methodologies did not require that end faces were strictly parallel, and both methods produced consistent results (mean repeatability of 4%). Both methods utilized the natural shape of the plant stem to avoid stress concentrations and buckling failure which are common challenges when performing compression tests, especially with thin-walled specimens.

Both elastic moduli measurements presented in this study neglected the contribution of the pith tissue. This assumption had a mean effect of overestimating the rind stiffness by 4%, which was deemed to be acceptable for these purposes.

Each of these methods possesses unique advantages and disadvantages. The overall compressive modulus technique provides a single, average value for all rind tissue in the specimen, but can be obtained relatively quickly. In contrast, the measurement of local modulus required multiple strain measurements, thus requiring additional tests, but provided results which are likely more accurate.

The modulus of elasticity values reported in this study are relevant from the stalk-level down to scales of a few centimeters. At scales smaller than this, the cellular architecture of the stalk tissue should be considered. Finally, although these measurements were developed and tested for dry maize specimens, the methods and principles introduced in this study are likely applicable for other types of plant stems, such as sorghum, reed, bamboo, etc.


  1. Von Forell G, Robertson D, Lee SY, Cook DD. Preventing lodging in bioenergy crops: a biomechanical analysis of maize stalks suggests a new approach. J Exp Bot. 2015;66:4367–71.

    Article  Google Scholar 

  2. Beer FP, Russell Johnston E, DeWolf JT, Mazurek DF. Mechanics of materials. 6th ed. New York: Mc Graw Hill; 2012.

    Google Scholar 

  3. Boresi AP, Schmidt RJ. Advanced mechanics of materials. 6th ed. New York: Wiley; 2003.

    Google Scholar 

  4. Gurtin ME. The linear theory of elasticity. In: Truesdell C, editor. Linear theories of elasticity and thermoelasticity. Berlin: Springer; 1973. p. 1–295.

    Google Scholar 

  5. Buchanan AH. Bending strength of lumber. J Struct Eng ASCE. 1990;116:1213–29.

    Article  Google Scholar 

  6. Kin K, Shim K: Comparison between tensile and compressive Young’s modulus of structural size lumber. In: World conference on timber engineering. Riva del Garda, Italy, 20–24 June 2010.

  7. Kretschmann DE. The influence of juvenile wood content on shear parallel, compression, and tension perpendicular to grain strength and mode I fracture toughness of loblolly pine at various ring orientation. For Prod J. 2008;58:89–96.

    Google Scholar 

  8. Lindstrom H, Harris P, Nakada R. Methods for measuring stiffness of young trees. Holz Als Roh-Und Werkstoff. 2002;60:165–74.

    Article  Google Scholar 

  9. Ince A, Ugurluay S, Guzel E, Ozcan MT. Bending and shearing characteristics of sunflower stalk residue. Biosyst Eng. 2005;92:175–81.

    Article  Google Scholar 

  10. Bashford LL, Maranville JW, Weeks SA, Campbell R. Mechanical-properties affecting lodging of sorghum. Trans ASAE. 1976;19:962–6.

    Article  Google Scholar 

  11. Esehaghbeygi A, Hoseinzadeh B, Khazaei M, Masoumi A. Bending and shearing properties of wheat stem of alvand variety. World Appl Sci J. 2009;6:1028–32.

    Google Scholar 

  12. Robertson DJ, Smith SL, Cook DD. On measuring the bending strength of septate grass stems. Am J Bot. 2015;102:5–11.

    Article  PubMed  Google Scholar 

  13. Robertson DJ, Julias M, Gardunia BW, Barten T, Cook DD. Corn stalk lodging: a forensic engineering approach provides insights into failure patterns and mechanisms. Crop Sci. 2015;55:2833–41.

    Article  CAS  Google Scholar 

  14. Robertson D, Smith S, Gardunia B, Cook D. An improved method for accurate phenotyping of corn stalk strength. Crop Sci. 2014;54:2038–44.

    Article  Google Scholar 

  15. Zhang LX, Yang ZP, Zhang Q, Guo HL. Tensile properties of maize stalk rind. BioResources. 2016;11:6151–61.

    CAS  Google Scholar 

  16. Yu M, Igathinathane C, Hendrickson J, Sanderson M, Liebig M. Mechanical shear and tensile properties of selected biomass stems. Trans ASABE. 2014;57:1231–42.

    Google Scholar 

  17. Yu M, Womac AR, Igathinathane C, Ayers PD, Buschermohle MJ. Switchgrass ultimate stresses at typical biomass conditions available for processing. Biomass Bioenergy. 2006;30:214–9.

    Article  Google Scholar 

  18. Varanasi P, Katsnelson J, Larson DM, Sharma R, Sharma MK, Vega-Sánchez ME, Zemla M, Loque D, Ronald PC, Simmons BA, et al. Mechanical stress analysis as a method to understand the impact of genetically engineered rice and Arabidopsis plants. Ind Biotechnol. 2012;8:238–44.

    Article  Google Scholar 

  19. Young SA, Clancy P. Compression mechanical properties of wood at temperatures simulating fire conditions. Fire Mater. 2001;25:83–93.

    Article  CAS  Google Scholar 

  20. Wright CT, Pryfogle PA, Stevens NA, Steffler ED, Hess JR, Ulrich TH. Biomechanics of wheat/barley straw and corn stover. Appl Biochem Biotechnol. 2005;121:5–19.

    PubMed  Google Scholar 

  21. Robertson DJ, Julias M, Lee SY, Cook DD. Maize stalk lodging: morphological determinants of stalk strength. Crop Sci. 2017;57:926–34.

    Article  Google Scholar 

  22. Stubbs CJ, Baban NS, Robertson DJ, Al-Zube LA, Cook DD. Bending stress in plant stems: models and assumptions. In: Geitmann A, Gril J, editors. Plant biomechanics—from structure to function at multiple scales. Berlin: Springer; 2018.

  23. Robertson DJ, Lee SY, Julias M, Cook DD. Maize stalk lodging: flexural stiffness predicts strength. Crop Sci. 2016;56:1711–8.

    Article  Google Scholar 

  24. ASTM-E9: Standard test methods of compression testing of metallic materials at room temperature. ASTM International, West Conshohocken, PA. (2009). Accessed 1 Aug 2017.

  25. ASTM-D695: Standard test method for compressive properties of rigid. ASTM International, West Conshohocken, PA. (2015). Accessed 1 Aug 2017.

  26. ASTM-F2150: Standard guide for characterization and testing of biomaterial scaffolds used in tissue-engineered medical products. ASTM International, West Conshohocken, PA. (2013). Accessed 1 Aug 2017.

  27. Sachs RM. Stem elongation. Annu Rev Plant Physiol. 1965;16:73–96.

    Article  CAS  Google Scholar 

  28. Cheng SK, Clarke EC, Bilston LE. The effects of preconditioning strain on measured tissue properties. J Biomech. 2009;42:1360–2.

    Article  PubMed  Google Scholar 

  29. Bowman SM, Keaveny TM, Gibson LJ, Hayes WC, Mcmahon TA. Compressive creep-behavior of bovine trabecular bone. J Biomech. 1994;27:301–10.

    Article  CAS  PubMed  Google Scholar 

  30. Caler WE, Carter DR. Bone creep-fatigue damage accumulation. J Biomech. 1989;22:625–35.

    Article  CAS  PubMed  Google Scholar 

  31. Keaveny TM, Guo XE, Wachtel EF, Mcmahon TA, Hayes WC. Trabecular bone exhibits fully linear elastic behavior and yields at low strains. J Biomech. 1994;27:1127–36.

    Article  CAS  PubMed  Google Scholar 

  32. ASTM-D143: Standard test methods for small clear specimens of timber. ASTM International, West Conshohocken, PA. (2014). Accessed 1 Aug 2017.

  33. Maranville J, Clegg M: Morphological and physiological factors associated with stalk strength. In: Rosenberg G, editor. Sorghum root and stalk rots: a critical review. Proceedings of the consultative group discussion on research needs and strategies for control of sorghum root and stalk rot diseases, Bellagio, Italy: ICRISAT, Patancheru, India; 1984. p. 111–8.

  34. Westfall PH, Henning KS. Understanding advanced statistical methods. Boca Raton, FL: Taylor and Francis Group; 2013.

    Google Scholar 

  35. Robertson D, Cook D. Unrealistic statistics: how average constitutive coefficients can produce non-physical results. J Mech Behav Biomed Mater. 2014;40:234–9.

    Article  PubMed  Google Scholar 

  36. Robertson DJ, Cook DD. Hyperelasticity and the failure of averages. In: Kruis J, Tsompanakis Y, Topping BHV, editors. Proceedings of the fifteenth international conference on civil, structural and environmental engineering computing. Stirlingshire: Civil-Comp Press; 2015.

    Google Scholar 

  37. Cook DD, Robertson DJ. The generic modeling fallacy: average biomechanical models often produce non-average results! J Biomech. 2016;49:3609–15.

    Article  PubMed  Google Scholar 

  38. NIST-TN1297: Guidelines for evaluating and expressing the uncertainty of NIST measurement results. National Institute of Standards and Technology. (1994). Accessed 1 Aug 2017.

  39. Zuber MS, Colbert TR, Darrah LL. Effect of recurrent selection for crushing strength on several stalk components in maize. Crop Sci. 1980;20:711–7.

    Article  Google Scholar 

  40. Green DW, Winandy JE, Kretschmann DE. Mechanical properties of wood. In: Wood handbook: wood as an engineering material. Gen Tech Rep FPL-GTR-113. Madison, WI: USDA, Forest Services, Forest Products Laboratory; 1999. p 45.

  41. Kretschmann DE: Mechanical properties of wood. In: Wood handbook: wood as an engineering material. Gen Tech Rep FPL-GTR-113. Madison: WI: USDA, Forest Service, Forest Products Laboratory; 1999. p 41–4.

  42. Niklas KJ. Responses of hollow, septate stems to vibrations: biomechanical evidence that nodes can act mechanically as spring-like joints. Ann Bot. 1997;80:437–48.

    Article  Google Scholar 

  43. Keller TS, Liebschner MA. Tensile and compression testing of bone. In: Yuehuei HA, Robert AD, editors. Mechanical testing of bone and the bone-implant interface. Boca Raton: CRC Press; 1999. p. 181.

    Google Scholar 

  44. Untaroiu CD. A numerical investigation of mid-femoral injury tolerance in axial compression and bending loading. Int J Crashworthiness. 2010;15:83–92.

    Article  Google Scholar 

  45. Heckwolf S, Heckwolf M, Kaeppler SM, de Leon N, Spalding EP. Image analysis of anatomical traits in stalk transections of maize and other grasses. Plant Methods. 2015;11:26.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Authors’ contributions

LA, DR, and DC designed the research and wrote the manuscript. LA, DR, and DC developed the experimental procedure of the approach. LA, JE, and WS performed the experimental procedure. All authors read and approved the final manuscript.


We thank Monsanto Company, St. Louis, MO, USA for providing the maize stalk samples used in this study. This work was funded in part by the National Science Foundation (Award # 1400973), and the U.S. Department of Agriculture (Award # 2016-67012-24685).

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Consent for publication

No consent was required in this study.

Ethics approval and consent to participate

No human subjects or animals were used in this study.


This study was supported by the National Science Foundation, Arlington, VA, USA (Grant # 1400973) and the U.S. Department of Agriculture, Washington, DC, USA (Grant # 2016-67012-24685). The funding agencies did not have any role in designing the study or in collecting, analyzing, interpreting the data or in writing the manuscript.

Guidelines and legislation

The authors confirm following local UAE import regulations to get the stalk samples imported from the USA. No permissions and/or licenses for the study are required.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Douglas D. Cook.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Al-Zube, L.A., Robertson, D.J., Edwards, J.N. et al. Measuring the compressive modulus of elasticity of pith-filled plant stems. Plant Methods 13, 99 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: