Materials and Design

Transcription

1 Materials and Design 32 (20) Contents lists available at ScienceDirect Materials and Design journal homepage: Energy absorption properties of multi-layered corrugated paperboard in various ambient humidities Zhi-Wei Wang a,c,, Yu-Ping E a,b,c a Packaging Engineering Institute, Jinan University, 206 Qianshan Road, Zhuhai 59070, PR China b College of Materials and Textiles, Zhejiang Sci-Tech University, Xiasha, Hangzhou 3008, PR China c Key Laboratory of Product Packaging and Logistics of Guangdong Higher Education Institutes, Jinan University, 206 Qianshan Road, Zhuhai 59070, PR China article info abstract Article history: Received 27 August 200 Accepted 8 January 20 Available online 8 February 20 Keywords: B. Sandwich structures E. Environmental performance H. Material selection charts This paper develops a mathematical model to depict the energy absorption properties of multi-layered corrugated paperboard (MLCP) in various ambient humidities. It is a piecewise function to model the energy absorptions corresponding to three deformation stages of MLCP (elastic stage, plateau stage and densification stage) separately. Simple formulas are derived for each stage which relating the energy absorption capacity of MLCP to the thickness-to-flute pitch ratio (t c =k) of corrugated-core cell, the mechanical properties of corrugated medium tested under a controlled atmosphere [23 C and 50% relative humidity (RH)], and the RH. The theoretical energy absorption curves are then compared with experimental ones and good agreements are achieved for MLCP with wide range ratios of t c =k in various ambient humidities. Results of this research can be applied in the optimum design and material selection of cushioning packaging with multi-layered corrugated paperboard. Ó 20 Elsevier Ltd. All rights reserved.. Introduction Corresponding author at: Packaging Engineering Institute, Jinan University, 206 Qianshan Road, Zhuhai 59070, PR China. Tel./fax: address: wangzw@jnu.edu.cn (Z.-W. Wang). Multi-layered corrugated paperboard (MLCP, as shown in Fig. a) is a structure formed by gluing a specified number of paralleled single-wall corrugated paperboards, which are sandwich structures comprised of one layer of corrugated medium and two layers of linerboard (as illustrated in Fig. b). In recent years, MLCP has been gaining increasing attention as a replacement for polymeric materials in protective packaging of product. These benefited mainly from the environmentally friendly quality and the ability to absorb impact loading by converting it into plastic deformation energy. Therefore, the evaluation of cushioning and energy absorption property of MLCP is essential to engineering applications. Energy absorption diagram as an effective tool to characterise the energy absorption capacity of cushioning materials has been reported being successfully applied to plastic foam and aluminium alloy material [ 6]. Wang [7] then employed this tool to examine the energy absorption property of MLCP under an atmosphere of 20 C in temperature and 65% in relative humidity (RH). In fact, MLCP is a paper-based packaging material and its energy absorption property may be influenced by the ambient humidity besides structural parameters. However, the existed literature in discussing the cushioning and energy absorption properties of MLCP were mostly focus on the structural factors [8 ], neglecting the impact of the relative humidity. Thus, this study aimed to investigate the effect of ambient humidity on the energy absorption characteristics of MLCP with various thickness-to-flute pitch ratios and to construct the relationship between the energy absorption of MLCP and ambient humidity, as well as the structural parameters thereof. For this purpose, the research work done by Wang and E [2] can serve as a reference for our study. 2. Modelling process 2.. Simplified constitutive model of MLCP The energy absorption property of MLCP under flatwise compression (Fig. a, loading in T-direction) was investigated in this paper. A corrugated-core unit was shown in the enlarged drawing in Fig. a, where t c denotes the thickness of corrugated medium and k signifies the flute pitch, h and h c are the flute height and the thickness of single-wall corrugated board, respectively. The thickness of linerboard is negligible compared with the total thickness of corrugated board, thus h c approximately equals to h. The stress strain curve of MLCP under flatwise compression exhibits four deformation stages (as can be seen in Fig. 2a), that is, elastic stage, first buckling stage, sub-buckling stage and densification stage. Elastic stage can be modelled by a perfect elastic material while the densification stage can be simplified as a vertical line. Compressive curve in first buckling stage and sub-buckling /$ - see front matter Ó 20 Elsevier Ltd. All rights reserved. doi:0.06/j.matdes

2 Z.-W. Wang, Y.-P. E / Materials and Design 32 (20) Nomenclature e pc r p0 E p e 0c e Dc r pd E E 0 E sc energy absorption per unit volume stress at arbitrary point on stress strain curve strain at arbitrary point on stress strain curve average crushing stress in first buckling stage slope of the curve in sub-buckling stage strain at the end of elastic stage densification strain of MLCP stress corresponding to the densification strain elastic modulus of MLCP tested under flatwise compression in arbitrary humidity environment elastic modulus of MLCP tested under flatwise compression in an environment of 23 C and 50% RH elastic modulus of corrugated medium tested in arbitrary humidity environment elastic modulus of corrugated medium tested in an environment of 23 C and 50% RH r ysc r 0 ysc = = t c k t c =k h h c n c yield strength of corrugated medium tested in arbitrary humidity environment yield strength of corrugated medium tested in an environment of 23 C and 50% RH normalised energy absorption per unit volume normalised stress thickness of corrugated medium flute pitch thickness-to-flute pitch ratio of corrugated-core cell flute height thickness of single-wall corrugated paperboard number of single-wall corrugated boards in a multi-layered corrugated paperboard Fig.. Structural diagram of (a) Multi-layered corrugated paperboard and (b) Single-wall corrugated board. stage is a wave-like uplift curve, but in the simplified material model, these two stages were collectively referred to as plateau stage. Although it was not a plateau in the strict sense, we borrowed this term in contrast to the deformation stages of paper honeycomb [2]. Constitutive relation of this stage can be modelled by curvilinear equation analogy to that of plastic foam [7,]. However, in order to facilitate the research, wave-like uplift curve in this stage was simplified as a line with the slope of E p. The

3 3478 Z.-W. Wang, Y.-P. E / Materials and Design 32 (20) Fig. 2. Typical stress strain curve of MLCP under flatwise compression. (a) Original stress strain curve (I. Elastic stage; II. First buckling stage; III. Sub-buckling stage; IV. Densification stage); (b) Simplified material model. simplified material model of MLCP was shown in Fig. 2b. Thus, the stress strain relationship of MLCP can be depicted in a piecewise function 8 >< ¼ E e pc ; 0 6 e pc 6 e 0c ¼ r p0 þ E p ðe pc e 0c Þ; e 0c 6 e pc 6 e Dc ðþ >: P r p0 þ E p ðe Dc e 0c Þ; e pc ¼ e Dc where and e pc are the stress and strain at arbitrary point on stress strain curve, respectively; r p0 denotes the average crushing stress of the first buckling; E signifies the elastic modulus of MLCP; e 0c is the strain at the end of elastic stage, and e Dc signifies the densification strain. Energy absorption per unit volume of MLCP can be obtained and compared by integration of the area under the simplified and the original stress strain curve prior to the onset of densification (shaded part in the inset of Fig. 3). As illustrated in Fig. 3, the data lied closer to the unit-slope straight line, indicating that the simplified material model proposed here is feasible and accurate in predicting the energy absorption of MLCP Model of energy absorption curve Energy absorption in elastic stage Energy absorption per unit volume of MLCP in this stage can be calculated by the following formula Z epc ¼ de pc ¼ r 2 pc 2 E : ð2þ 0 Normalised by, the elastic modulus of corrugated medium tested under a controlled environment (23 C, 50% RH), then we obtain ¼ 2 E! 2 ; ð3þ

4 Z.-W. Wang, Y.-P. E / Materials and Design 32 (20) y=x Experiments of energy absorption per unit volume WC /J/cm Stress σ pc /MPa σ p0 Energy absorption E p per unit volume Predictions of energy absorption per unit volume /J/cm 3 ε 0c Strain ε pc ε Dc Fig. 3. Comparison of the total energy absorption per unit volume between the simplified stress strain relationships and the experimental ones. where E is the elastic modulus of MLCP in arbitrary RH. Defining the elastic modulus of MLCP tested under the controlled environment aforementioned as E 0, and then assuming the relative elastic modulus of MLCP (the elastic modulus ratio of MLCP tested in arbitrary atmosphere to that under the controlled atmosphere) changed with the ambient humidity as follows: E E 0 ¼ f ðrhþ; where f (RH) is to be determined by experiments. The elastic modulus ratio of MLCP to corrugated medium varies proportionally to the cube of the thickness-to-flute pitch ratio [8], that is: E 0 3 k 0 t c ; ð5þ k Where k 0 is a constant. After substituting Eqs. (4) and (5) into Eq. (3), we have! 2 2 f ðrhþ k 0 ðt : ð6þ c=kþ 3 Let C C ¼ 2 C C f ðrhþ. Eq. (6) can be simplified as: k 0 ðtc=kþ3 ð4þ! 2 : ð7þ As can be seen from above equation, in elastic stage, the relationship between the normalised energy absorption per unit volume ð = Þ and the normalised stress ð=þ for MLCP with different t c =k ratios under various ambient humidities is a family of inclined lines with the same slop of 2 in logarithmic coordinates, but their intercepts depend upon the t c =k value of corrugated-core cell and the ambient humidity Energy absorption in plateau stage MLCP absorbs a predominant portion of kinetic energy in this stage. By neglecting the energy absorption in elastic stage, the energy absorption capacity can be approximately evaluated by ¼ r p0ð r p0 Þ þ ð r p0 Þ 2 ¼ r 2 pc r2 p0 ; ð8þ E p 2 E p 2 E p where r p0 and E p can be computed by following equations [3] 2 t r P0 ¼ D C r c t ysc ¼ r c yscg k ; ð9þ k E p ¼ 2:249r ysc expð :0448kÞ ¼r ysc g 2 ðkþ; ð0þ where g ð tc Þ¼D k Cð tc k Þ2, and g 2 ðkþ ¼2:249 expð :0448kÞ. D C is a constant related to the configuration of the corrugated flute, and equals to 8.923, and.3289 corresponding to the MLCP with A-flute, B-flute and C-flute, respectively; r ysc signifies the yield strength of corrugated medium. Denoting r 0 ysc as the yield strength of corrugated medium tested under the controlled environment, then supposing the relative yield strength of corrugated medium (the yield strength ratios of corrugated medium measured in arbitrary environment to that under the controlled atmosphere) varied with the RH as follows: r ysc ¼ f r 0 2 ðrhþ; ysc ðþ where f 2 (RH) is an undetermined function. Substituting Eqs. (9) () into Eq. (8) and then normalised by, the following is obtained: ¼ 2 r 0 ysc f 2 ðrhþ g 2 ðkþ! 2 r 0 ysc 2 g 2 ðt c=kþ f 2 ðrhþ: g 2 ðkþ ð2þ Let C 2C ¼ 2 r, and C 0 f ysc 2 ðrhþ g 2 ðkþ 3C ¼ r 0 ysc g 2 2 E 0 ðtc=kþ f g sc 2 ðkþ 2 ðrhþ. Then, Eq. (2) can be rewritten as:! 2 ¼ C 2C C 3C: ð3þ Therefore, in plateau stage, the energy absorption curve of MLCP with different t c =k ratios in various ambient humidities is a family of parabolas in linear coordinations, while a train of similar curves in logarithmic coordinates.

5 3480 Z.-W. Wang, Y.-P. E / Materials and Design 32 (20) Energy absorption in densification stage The energy absorbed in densification stage equals to that at the end of plateau stage, which can be predicted by substituting e Dc into Eq. (3)! 2 r pd ¼ C 2C C 3C: ð4þ The stress corresponding to the densification strain () can be assessed by: r pd ¼ r p0 þ E p ðe Dc e 0c Þ: ð5þ Substituting Eq. (5) into Eq. (4), incorporating Eqs. (9) and (0), and ignoring the negligible t c =k small elastic strain e 0c, we get! r ysc C 2C f 2 2 ðrhþ g t c þ g k 2 ðkþe Dc C 3C: ð6þ Let C 4C ¼ C 2C ð r0 ysc Þ 2 f E 0 2 2ðRHÞðg ð tc Þþg k 2ðkÞe Dc Þ 2 C 3C. Then, the normalised energy absorption per unit volume in densification stage can sc be represented as: C 4C ð7þ Thus, the energy absorption curve in densification stage is a family of horizontal lines; their levels are determined by the configuration of corrugated-core cell, the mechanical characteristics of corrugated medium and the RH. Combining Eqs. (7), (3), and (7), then we establish the energy absorption model for MLCP with any ratio in arbitrary ambient humidities: 8 >< >: C C ð rp Þ 2 ; E 0 s ¼ C 2C ð rp Þ 2 C E 0 3C ; s C 4C ; where 8 C C ¼ 2 C 2C ¼ 2 >< >: C 3C ¼ 2 f ðrhþ r 0 ysc r 0 ysc k 0 ðtc=kþ3 f 2 ðrhþ g 2 ðkþ r p0 g 2 ðtc=kþ f g 2 ðkþ 2 ðrhþ 6 r p0 6 6 r pd P r pd C 4C ¼ C 2C ð r0 ysc Þ 2 f E 0 2 2ðRHÞðg ð tc Þþg k 2ðkÞe Dc Þ 2 C 3C sc ð8þ : ð9þ Above two equations connect the energy absorption capacity of MLCP with the thickness-to-flute pitch ratio of corrugated-core cell, the mechanical characteristics of corrugated medium and the RH. Therefore, we can predict the energy absorption capacity for a given MLCP applied in known humidity environments. 3. Experiment materials and method 3.. Materials Both corrugated medium and single-wall corrugated board were supplied by the professional manufacturing industry of Jackson Packaging Co. Ltd. (Dongguan, China). The corrugated medium was made from recycled fibre, and grammages were 05, 20 and 45 g/m 2 (the corresponding thickness are, respectively, 0.5, 0.9 and 0.22 mm). Single-wall corrugated board was composed of two layers of linerboard with 50 g/m 2 Kraft paper and one layer of fluting made from corrugated mediums aforesaid, and they were glued together by starch adhesives along the outsides of the peaks and flute valleys. Three common used types of corrugated-core, A- flute, B-flute and C-flute, were employed in this research, and they differed from each other in the flute pitch and the flute height (as shown in Fig. b). Multi-layered corrugated paperboard was obtained manually by gluing a designated number of identical single-wall corrugated boards. The inner-liner of single-wall corrugated board was pasted to the outer-liner of single-wall corrugated board next, and the flute peaks or valleys of each layer should keep on a straight line. The adhesive was PVA emulsion (white latex). Totally eight kinds of MLCP were used here, and the specifications of each kind of MLCP specimen were listed in Table. In the designation system of specimens, the number in middle of the symbol denotes the basis weight of corrugated medium (3, 5 and 7 stand for the corrugated medium with basis weight of 05, 20 and 45 g/m 2, respectively), letter A on either side of these numbers signifies the linerboard with the basis weight of 50 g/m 2, and letters A, B and C in parenthesis indicate the flute type. In order to investigate the influence of the t c =k ratio and the ambient humidity on the energy absorption properties of MLCP, the number of single-wall corrugated board for each MLCP specimen was designated as eight layers. Therefore, the total thicknesses of MLCPs with A-flute, B-flute and C-flute approximately equal to 40, 25 and 3 mm, correspondingly Experiment methods Experiment equipments The test equipments involved were a temperature humidity programmable controller (GDJS-00, Zhongya Test Experiment Co. Ltd., Huaian, China) and two universal material testing machines with different full-scale load range (CMT 8502, used for testing the mechanical properties of corrugated medium, full-scale load range of 500 N, MTS Systems Co. Ltd., Shenzhen, China; and WDW-0C, used for testing the compression characteristics of MLCP, full-scale load range of 0 kn, Hualong Test Instruments Co. Ltd., Shanghai, China) Experiment standards ISO Paper, board and pulps: standard atmosphere for conditioning and testing and procedure for monitoring the atmosphere and conditioning of samples [4]; GB/T Plastic-determination of tensile propertiespart : general principles [5]; GB/T Testing method of static compression for package cushioning materials [6]; GB/T Test method for flatwise compression properties of sandwich constructions or cores [7]. Table Specifications and characteristics of MLCP specimens (n c = 8). Specimen Basic weight of linerboard (g/m 2 ) Basis weight of corrugated medium (g/ m 2 ) Thickness of corrugated medium t c (mm) Flute type A3A(A) A 40 A3A(B) B 25 A3A(C) C 3 A5A(A) A 40 A5A(B) B 25 A5A(C) C 3 A7A(B) B 25 A7A(C) C 3 Thickness of MLCP T (mm)

6 Z.-W. Wang, Y.-P. E / Materials and Design 32 (20) Experiment conditions All specimens were pre-conditioned in the temperature humidity programmable controller to reach their own equilibrium moisture content in each experimental RH, i.e., 40%, 50%, 65%, 75%, 85% and 95%, at fixed temperature of 23 C for 48 h prior to tests. The corrugated medium was tested to determine its elastic modulus and yield strength under different environmental conditions. All specimens were cut into the same area of mm 2, with a gauge length of 65 mm. The long edge of specimen was paralleled to the machine direction (MD) of corrugated medium. The tests were conducted under a constant displacement velocity of ± 0.5 mm/min, and every test was conducted in six repeats. Linear regression of stress strain curve between strains of 0.05 and 0.25 was defined as the elastic modulus of corrugated medium, and the stress corresponding to strain of 0.3 was taken as the yield strength, as recommended by China National Standard GB/T MLCP was tested under quasi-static compression under different environmental conditions to determine their energy absorption characteristics as response to ambient humidities. The uniaxial compression characteristics of MLCP may be affected by the dimension effect. Some researchers have studied the dimension effect in metal foam and metal honeycomb [8 20], the results advised that the boundary-layer effects will vanish when the specimens larger than seven cell sizes. Thus for corrugated-core cell, the side length of specimen should not less than 60 mm. By considering the dimension of the compression apparatus we employed in this study, the dimension of specimen was g/m 2 20g/m 2 45g/m Fig. 4. Relationship between the relative yield strength of corrugated medium and the ambient humidities..2.0 Experiments Fitting curve Fig. 5. Relationship between the relative elastic modulus of MLCP and the ambient humidities.

7 3482 Z.-W. Wang, Y.-P. E / Materials and Design 32 (20) Table 2 Densification strain of MLCP in different humidities (23 C, n c = 8). Specimen Thickness-to-flute pitch ratio t c=k Densification strain e Dc 40% RH 50% RH 65% RH 75% RH 85% RH 95% RH Average A3A(A) A3A(C) A5A(A) A5A(C) A3A(B) A7A(C) A5A(B) A7A(B) settled to be mm 2 finally. The compressive load was applied with a constant displacement velocity of 2 ± 3 mm/ min. Every compression test was effectively conducted in five repeats, and each measurement was completed within 5 min after the specimen being removed from the temperature humidity programmable controller. Linear regression of stress vs. strain data located in the elastic stage was defined as the elastic modulus of MLCP. Fig. 6. Theoretical energy absorption curves of (a) MLCP with the same thickness-to-flute pitch ratio under various relative humidities; and (b) MLCP with different thickness-to-flute pitch ratios under a controlled atmosphere (23 C, 50%RH).

8 Z.-W. Wang, Y.-P. E / Materials and Design 32 (20) Results and discussion 4.. Mechanical characteristics of corrugated medium and MLCP According to the experimental results, the elastic modulus and the yield strength of corrugated medium tested under the controlled environment were [3] 05 g=m 2 : ¼ 2:943 GPa; r0 ysc ¼ 8:7 MPa; 20 g=m 2 : ¼ :904 GPa; r0 ysc ¼ 5:73MPa; 45 g=m 2 : ¼ 2:582 GPa; r0 ysc ¼ 7:349 MPa: ð20þ The yield strength of corrugated medium in arbitrary environment was measured by experiments [3]. Then, the relative yield strength of corrugated medium was computed by using Eq. () and plotted against the RH in Fig. 4, and the empirical expression of the function f 2 (RH) was acquired by curve fitting, :00 f 2 ðrhþ ¼ þ expð RH 92:038Þ ð2þ :206 The elastic modulus of MLCP in arbitrary and the controlled environment aforementioned were tested by experiments [3]. Then, the relative elastic modulus of MLCP was computed by using Eq. (4) and plotted against the RH in Fig. 5. By applying curve fitting to these data, the expression of the function f (RH) was obtained: 0:982 f ðrhþ ¼ : ð22þ þ exp RH 86:888 8:058 According to Eq. (5), the elastic modulus ratio of MLCP to corrugated medium varies proportionally to the cube of the thickness-to-flute Fig. 7. Energy absorption diagram of MLCP based on the effect of relative humidity (solid line stand for the contour with the same value of t c=k while dashed line represent the contour with constant level of RH). Fig. 8. Comparison of the energy absorption curves of MLCP between the theoretical one and the experimental one (t c=k = , 23 C, 50%RH).

9 3484 Z.-W. Wang, Y.-P. E / Materials and Design 32 (20) pitch ratio [8]. Thus, a linear approximation dependency between the ratio of E 0 to and the cube of t c=k value was obtained by experiments [3], E 0 3 : ¼ 2:908 t c k ð23þ The densification strain is defined as the onset of densification and taken to be the turning point from sub-buckling stage to densification stage. The densification strains of MLCP in various relative humidities were listed in Table 2. As can be seen from these data, the densification strain was inert to relative humidity, thus the experiment results under six levels of RH can be averaged and taken as the densification strain of this kind of MLCP. The average densification strains of MLCP with different configurations were inversely related to t c =k except for two small fluctuations, and their relationship can be represented as: e Dc ¼ 0:80 2:39 t c k : ð24þ 4.2. Energy absorption diagram of MLCP Two steps were involved to acquire the energy absorption curves for a given MLCP. Firstly, parameters of C C, C 2C, C 3C and C 4C in Eq. (9) have to be computed with the help of the known expressions of f (RH) and f 2 (RH) and the known values of t c =k, as well as the RH; and then substituting the calculated C C, C 2C, C 3C and C 4C into Eq. (8), the relationship between the normalised energy absorption in each stage and the normalised stress can be determined. Secondly, the normalised energy absorption per unit volume is plotted against the normalised stress for a range of thickness-to-flute pitch ratios in logarithmic coordinates by using MAT- LAB Ò software (version 7. R4, The Math Works Inc., USA), and then the energy absorption curves were constructed. Energy absorption curves obtained from the theoretical model aforesaid were plotted to illustrate the energy absorption properties of MLCP with the specified thickness-to-flute pitch ratio of under different surrounding environments (respectively at a b Fig. 9. Energy absorption curves obtained from experiments. (a) MLCP with the same thickness-to-flute pitch ratio under various relative humidities; (b) MLCP with different thickness-to-flute pitch ratios under a controlled atmosphere (23 C, 50% RH).

10 Z.-W. Wang, Y.-P. E / Materials and Design 32 (20) %, 50%, 65%, 75%, 85% and 95% in RH with the temperature fixed at 23 C) in Fig. 6a, and to illustrate the energy absorption properties of MLCP with different values of t c =k under the controlled environment in Fig. 6b. Inflexions of these curves signify the optimum points for energy absorption properties of MLCP. As illustrated in Fig. 6a, the inflexions of energy absorption curves for MLCP with the same t c =k ratio shift to the lower-left corner with the increasing of RH, suggesting a weakening of load carrying property and energy absorption capacity of MLCP. Fig. 6b indicates that the inflexions of energy absorption curves for MLCP used under the same relative humidity of 50% move to the higher-right corner as the increasing of the t c =k ratio. The energy absorption diagram of MLCP illustrated in Fig. 7 is constructed by connecting the inflexions of energy absorption curves with the same t c =k value and the inflexions of energy absorption curves having constant RH level respectively. As illustrated in this figure, energy absorption properties of MLCP increase proportionally to the t c =k ratios on the contours with constant RH levels, while inversely proportional to the ambient humidities on the contours with the same t c =k ratios. Eight kinds of MLCPs (thickness-to-flute pitch ratios are 0.085, , , , 0.029, , and , respectively) were compressed under six levels of relative humidity (40%, 50%, 65%, 75%, 85% and 95% in RH, respectively, with the temperature fixed at 23 C). Therefore, eight contours with the same value of and six contours having constant RH level were obtained from experiments. However, in order to make the energy absorption diagram of MLCP clear and distinguishable, only eight contours were shown in Fig. 7 to represent energy absorption properties of MLCPs with four t c =k ratios of 0.085, , and under four levels of RH (i.e., 40%, 65%, 85% and 95%). By using these contours, the optimal points of energy absorption for MLCP with arbitrary value of t c =k in prescribed humidity environment can be acquired by means of interpolation Verification of model Take the MLCP with the thickness-to-flute pitch ratio of for example, theoretical and experimental energy absorption curves were plotted and compared in Fig. 8. It can be seen that the theoretical energy absorption curve correlated closely with the experimental one particularly in plateau stage, indicating the energy absorption model established here can well simulate the energy absorption performance of MLCP. The theoretical curve located at the upper left of the experimental one in elastic stage, and lowered than the experimental one in densification stage. These differences were mainly caused by the simplification of material model of MLCP. Energy absorption curves of MLCP with different t c =k ratios tested under various ambient humidities (40%, 50%, 65%, 75%, 85% and 95%, respectively) were plotted in Fig. 9. These curves were converted from the original stress strain curves by integrating technique using the MATLAB Ò program, and then compared with the theoretical energy absorption curves in Fig. 6. Results show that the theoretical model can reflect the influence of relative humidity and thickness-to-flute pitch ratio on the energy absorption capacity of MLCP precisely. Furthermore, both theoretical and experimental energy absorption curves have proved that the energy absorption properties of MLCP decrease with the increasing of RH, and increase with the increasing of the thickness-to flute pitch ratio. 5. Conclusions pitch ratios under various ambient humidities simulating the actual logistics environment. The model was normalised by the elastic modulus of corrugated medium tested under the controlled atmosphere, which makes the diagram universal to MLCP made of different base materials. A correlation between the theoretical and experimental energy absorption curves indicated that the theoretical model can be used to evaluate the energy absorption properties of MLCP in various ambient humidities. Both theoretical and experimental energy absorption curves show that the energy absorption properties of MLCP decrease with the increasing of RH, and increase with the increasing of the thickness-to-flute pitch ratio. The model presented here can facilitate the construction of energy absorption curves of MLCP; therefore, we can quickly seek the optimal energy absorption points of MLCPs with different thickness-to-flute pitch ratios in various humidity conditions. Acknowledgements The authors acknowledge the financial support of this research by Project supported by National Natural Science Foundation of China. References [] Maiti SK, Gibson LJ, Ashby MF. Deformation and energy absorption diagrams for cellular solids. Acta Metal 984;32(): [2] Prakash O, Sang H, Embury JD. Structure and properties of AL-SiC foam. Mat Sci Eng A Struct 995;99(2): [3] Gibson LJ, Ashby MF. Cellular solids: structure and properties. second ed. Cambridge: Cambridge University Press; 997. [4] Avalle M, Belingardi G, Montanini R. Characterization of polymeric structural foams under compressive impact loading by means of energy-absorption diagram. Int J Impact 200;25(5): [5] Lu GX, Yu TX. Energy absorption of structures and materials. Woodhead Publishing Ltd. and CRC Press LLC; [6] Wang ZH. Studies on the dynamic mechanical properties and energy absorption of aluminum alloy foams. Doctoral dissertation of Taiyuan University of Technology. Taiyuan, China; 2005 [in Chinese]. [7] Wang DM. Energy absorption diagram of multi-layer corrugated boards. J Wuhan Univ Technol (Mater Sci Ed) 200;25(2):58 6. [8] Wang DM. Cushioning property and characteristic studies on honeycomb paperboards and corrugated paperboards. Doctoral Dissertation of Jiangnan University. Wuxi; 2007 [in Chinese]. [9] Pan XZ. Study on the structural properties of corrugated cardboard composites. Master s Dissertation of Nanjing Forestry University. Nanjing; 2007 [in Chinese]. [0] Naganathan P, He J, Kirkpatrick J. The effect of compression of enclosed air on the cushioning properties of corrugated fiberboard. Packag Technol Sci 999;2(2):8 9. [] Minett M. A study of air flow effects on the cushioning characteristics of multi layered pre-compressed fiberboard. Doctoral dissertation of Victoria University. Melbourne, Australia; [2] Wang Z-W, E YP. Mathematical modelling of energy absorption property for paper honeycomb in various ambient humidities. Mater Des 200;3(9): [3] E YP. Influence of relative humidity and strain rate on the energy absorption properties of paper-based cushioning materials. Doctoral Dissertation of Jiangnan University. Wuxi; 200 [in Chinese]. [4] ISO paper, board and pulps: standard atmosphere for conditioning and testing and procedure for monitoring the atmosphere and conditioning of samples. International Organization for Standardization; 990. [5] GB/T Plastic-determination of tensile properties-part : general principles. Beijing: Standards Press of China; 2006 [in Chinese]. [6] GB/T Testing method of static compression for package cushioning materials. Beijing: Standards Press of China; 2008 (in Chinese). [7] GB/T Test method for flatwise compression properties of sandwich constructions or cores. Beijing: Standards Press of China; 2005 [in Chinese]. [8] Andrews EW, Gioux G, Onck PR, et al. Size effects in ductile cellular solids. Part II: experimental results. Int J Mech Sci 200;43(3):70 3. [9] Hayes AM, Wang AJ, Dempsey BM, et al. Mechanics of linear cellular alloys. Mech Mater 2004;36(8): [20] Onck PR. Scale effects in cellular metals. MRS Bull 2003;28(4): A mathematical model was developed to predict the energy absorption properties of MLCPs with different thickness-to-flute