# A Phenomenological Thermal-Mechanical Viscoelastic Constitutive Modeling for Polypropylene Wood Composites

Authors: Xiongqi Peng; Hongling Yin; Jun Chen; Xue Liu
Journal: Advances in Materials Science and Engineering
 (2012)

---

## Abstract

This paper presents a phenomenological thermal-mechanical viscoelastic constitutive modeling for polypropylene wood composites. Polypropylene (PP) wood composite specimens are compressed at strain rates from 10−4 to 10−2 s−1 and at temperature of 90∘C, 130∘C, and 170∘C, respectively. The mechanical responses are shown to be sensitive both to strain rate and to temperature. Based on the Maxwell viscoelastic model, a nonlinear thermal-mechanical viscoelastic constitutive model is developed for the PP wood composite by decoupling the effect of temperature with that of the strain rate. Corresponding viscoelastic parameters are obtained through curve fitting with experimental data. Then the model is used to simulate thermal compression of the PP wood composite. The predicted theoretical results coincide quite well with experimental data. The proposed constitutive model is then applied to the thermoforming simulation of an automobile interior part with the PP wood composites.


---

## Body

## 1. Introduction

Wood plastic composites (WPCs) are a new class of composite materials which combine the characteristics of plastic and wood fibers. In appearance, WPC sheets are similar to wood sheets, but the low stiffness of plastics makes the modulus of composites significantly lower than that of solid woods, and thus brings in a greatly improved formability. This makes it feasible to form WPC sheets into products with complex curvatures. WPCs have many additional advantages such as low cost, lightweight, and renewability. The increasing environmental concern has created enormous opportunities for introducing this kind of nature fiber plastic composites into high-volume consumer products such as automobiles, especially for car interior parts [1–4].Nowadays, automobile interior parts made from WPCs are primarily manufactured by the process of thermoforming which is well suited for mass production with a cost-effective way. In thermoforming processes, stress relaxation and viscoelastic properties [5–8] of WPCs are very important. Temperature and pressure conditions in the forming process should be determined carefully taking into account the thermoviscoelastic property of WPCs. High temperature (about 170°C) makes wood plastic composites show a strong nonlinear viscoelastic material behavior during forming. Under an inadequate forming condition, WPC specimen may be damaged or a high residual stress may present in the final formed parts. Besides, inappropriate processing parameters may incur quality problems such as serious thinning and insufficient filling defects. Many research efforts have been devoted to understanding the WPC thermoforming mechanism for solving process problems and optimum forming process design, either by experiment approaches [9, 10] or with numerical simulation techniques [11–16].Thermoforming of WPCs is a thermal-mechanical coupled process which involves complex phenomenon such as nonlinear material behavior with large deformation, heat transfer, viscosity, and flow. Numerical simulations of WPC thermoforming process could provide a powerful tool for materials selection, tool design, and process optimization. Material models for WPCs constitute the cornerstone for numerical simulation of WPC thermoforming process. To obtain an accurate and meaningful simulation result, it is crucial to have an appropriate material model for WPCs, which is capable of reflecting the material behavior of WPCs rendering in thermoforming. By using cellular solid theories, Wolcott et al. [17] characterized the nonlinear compression behavior of wood composites under various temperature and moisture contents. Oudjehane et al. [18] developed an analytical nonlinear continuum mechanics based model for predicting the consolidation of wood flake mats taking into account the effect of sheet thickness and press closing rate. By combining the homogenization method and a network mechanics model, Stålne and Gustafsson [19] proposed a three-dimensional model for characterizing the stiffness and hygroexpansion of fiber composite materials. Pooler and Smith [20] carried out a series of creep and recovery tests on High Density Polyethylene (HDPE) wood composites and then used a Prony series to describe the nonlinear viscoelastic of the wood-plastic composites. Damage effects were taken into account by using an effective stress. Bogren et al. [21] measured the dynamic stress transfer between fiber reinforcement and matrix ground materials to investigate the influence of fiber-matrix interface properties on the macroscopic mechanical properties of wood-fiber reinforced polylactide. A viscoelastic micromechanical model based on perfect interfacial stress transfer was taken for a comparison. Pramanick and Sain [22] developed a generic creep prediction model to describe the nonlinear viscoelastic behavior of HDPE wood composites. Erchiqui et al. [23] used a thermodynamical approach to express external work in terms of enclosed gas volume to investigate the effect of latter on the mechanical properties on wood-plastic composites. The Lodge model was selected for characterizing the viscoelastic behavior.In this study, the thermoviscoelasticity of PP wood composites is characterized by unidirectional thermocompression tests with various strain rates from 10−4 to 10−2 s−1 and under temperature 90°C, 130°C, and 170°C, respectively. A nonlinear thermomechanical viscoelastic constitutive model for PP wood fiber composite is developed with a phenomenological approach by using the Maxwell model. The effect of temperature is decoupled with that of the strain rate to greatly facilitate the material characterization process for the PP wood composite. Corresponding viscoelastic parameters of the PP wood composite are obtained through curve fitting with experimental results. Model validation is implemented by comparing numerical stresses with experimental data. The developed thermomechanical viscoelastic constitutive model is then applied to the thermoforming simulation of an automobile interior part with the PP wood fiber composite.


## 2. Thermocompression Test and Results

The compression test is set up according to the ASTM D6641(M) standard. A PP wood fiber composite sheet with thickness of 2.57 mm is cut into10mm×10mm specimens. A CMT4000 series electronic universal testing machine with an oven and a digital temperature controller is used to carry out the unidirectional thermocompression tests for the PP wood fiber composite. Figure 1 shows the experimental set-up along with some test samples. Since the thermoforming of the PP wood composites is performed above the glass transition temperature (Tg) of PP which is about 80°C, the compression tests are accordingly performed under different temperature (90°C, 130°C, and 170°C), which are all above the Tg of PP, for the purpose of characterizing the mechanical behavior of the WPC during forming. Three different loading rates (0.1 mm/min, 0.5 mm/min, and 5 mm/min) are chosen according to two reasons. First, a very small rate (0.1 mm/min) is chosen for a quasistatic test. The other two loading rates are set to reflect the effect of strain rates on the material response as well as to cover the forming speed range of the wood composite. The tests are approximated as constant strain rate. Four specimens are prepared for each temperature and loading rate combination case.Figure 1
Thermo-compression test for PP wood composites.The load-displacement curves for each thermocompression test are recorded and then averaged and converted to true stress-strain curves. Figure2 shows the true stress-strain curves under different temperature and different strain rate. The trends of the true stress-strain curves are similar. Initially, it exhibits linear elastic deformation in a small strain range. When the strain increases to some value, it exhibits nonlinear behavior. There is no apparent yield point in the whole thermocompression stress-strain curves.True stress-strain curves under compression with different strain rate (s−1).
(a)
(b)
(c)The stress-strain curves in Figure2 indicated that the mechanical responses of PP wood fiber composites are sensitive to strain rate and temperature as well. With the temperature increasing, the maximum stress decreases. In addition, increasing strain rate has an equivalent effect on the stress as reducing the temperature. Figure 2 also shows that the effect of temperature on the stress is more prominent than the strain rate.


## 3. Viscoelastic Constitutive Model for PP Wood Composite

In PP wood fiber composites, the PP resin matrix shows viscoelastic characteristics. Viscoelastic behavior in wood is not as prominent as it is in polymer. Hence, for simplification of determining model constants, it is assumed that wood fiber renders elastic properties. The Maxwell viscoelastic model is selected to characterize the comprehensive mechanical properties of the PP wood fiber composites. Figure3 shows the schematic of the Maxwell viscoelastic model. In Figure 3, E0, E1, and η are the to-be-determined viscoelastic parameters.Figure 3
The generalized Maxwell viscoelastic model for PP wood composites.In Figure3, the stress of the spring is equal to that of the damper in the serial connection part. The total strain ε in the serial connection part is the summation of the spring strain ε1 and the damper strain ε2, that is,
(1)ε=ε1+ε2,ε˙=σ˙1E1+σ1η,
where E1 is the equivalent Young’s modulus of PP with a unit of Pa, η is a viscosity with a unit of N-s/m2. Combining  leads to
(2)σ˙1+σ1θ=E1ε˙,
where θ=η/E1 is the relaxation time of the PP wood composite. The stress σ1 can be obtained by convolution integral as
(3)σ1=∫-∞tE1exp(-t-τθ)ε˙(τ)dτ.
In large deformation, the stress σ2 of the spring in the parallel part is usually nonlinear. Based on the pattern of the stress-strain curves shown in Figure 2, σ2 can be expressed in an exponential format as
(4)σ2=E0εn,
where E0 is the Young’s modulus of the wood fiber with a unit of Pa, n is a material constant. The total stress of the model is then the summation of σ1 and σ2(5)σ=E0εn+∫-∞tE1exp(-t-τθ)ε˙(τ)dτ,
where ε˙ stands for strain rate. Experimental studies [24, 25] on wood fiber behavior indicated that the mechanical properties of wood fiber have insignificant variation under different temperature and testing rate. Hence, it is assumed that the elastic response of the wood fiber (σ2) is independent of temperature and strain rate, and the elastic behavior of the PP resin (the spring in the serial connection part) is affected by temperature only. As a result, the elastic parameters E0 and n are constants independent of temperature and strain rate, and the parameter E1 is a function of temperature only. Taking advantage of these assumptions and simplifications and implementing least square fitting on each individual experimental stress-strain curve, the material parameters in the proposed nonlinear thermal-mechanical viscoelastic constitutive model for the PP wood composite can be obtained. The elastic material constants E0 and n are
(6)E0=3.37MPa,n=1.73.
Other material parameters are listed in Table 1.Table 1
Viscoelastic parameters of PP wood composites under different temperature.
ε
˙ (s−1)
90°C
130°C
170°C
E
1/MPa
θ/s
E
1/MPa
θ/s
E
1/MPa
θ/s
0.0006
320.1
92.42
183.1
94.66
58.5
91.5
0.0030
21.63
21.67
22.56
0.0303
2.66
2.63
3As shown in Figure4, the viscoelastic modulus E1 linearly changes with temperature for the testing temperature levels, which cover the forming temperature. Hence, E1 (MPa) can be obtained by curve fitting as
(7)E1=-3.27T+612.3,
where T denotes temperature with a unit of degrees of Celsius. The relaxation time θ under various strain rate and temperature is shown in Figure 5. As indicated by Table 1 and Figure 5, the relaxation time θ is insensitive to temperature and can be assumed as a function of strain rate only. By curving fitting in Figure 5, the relaxation time θ can be expressed as
(8)θ=0.1138(ε˙)0.9044.
In summary, the nonlinear thermal-mechanical viscoelastic constitutive model for the PP wood composite can be expressed as
(9)σ=3.37ε1.73+∫-∞t(-3.27T+612.3)×exp(-(t-τ)0.1138ε˙0.9044)ε˙(τ)dτ.Figure 4
Relations between viscoelastic modulusE1 and temperature.Figure 5
Relations between relaxation timeθ and strain rate.


## 4. Model Validation and ThermoformingSimulation

The developed nonlinear thermoviscoelastic constitutive model for the PP wood composite is implemented as a UMAT in LS-Dyna. The previous thermocompression tests are simulated with the developed viscoelastic model. Figure6 shows the comparison of the numerical stresses predicted from the viscoelastic model with corresponding experimental data for each test. In Figure 6, the solid lines denote numerical stresses and the discrete data points represent experimental stresses. As can be seen from Figure 6, in general, the numerical stresses are in a good agreement with experimental data, with an error less than 10%.Figure 6
Comparison of numerical stress of thermocompression tests with experimental data.To further validate the developed nonlinear thermoviscoelastic constitutive model, thermocompression tests are carried out under strain rate 0.012 s−1 with temperature 90°C, 130°C, and 170°C, respectively. Simulation results based on the developed constitutive model are compared with corresponding experimental data in Figure 7. The solid lines in Figure 7 represent numerical stresses and the discrete data points are experimental stresses under strain rate 0.012 s−1 with various temperature. As presented in Figure 7, reasonable agreement is obtained in a sense of overall. It demonstrates that the proposed constitutive model can accurately and effectively characterize the nonlinear thermal-mechanical viscoelastic material behavior of the PP wood composite.Figure 7
Comparison of simulation results with experimental data at strain rate 0.0121.Finally, the thermoforming simulation of a car interior part of the PP WPC is implemented with the viscoelastic material model. The forming system and its FEM model are shown in Figure8. Experimental study on thermoforming of the car interior part will be carried out in future work. The WPC composite sheet size is 1750mm×810mm with a thickness of 2.75 mm. The sheet is heated up to a temperature above the glass transition temperature of PP, then pressed and cooled to room temperature to form the final part.Figure 8
Thermal forming setup.Figure9 presents the thickness distribution contour of the formed part under an initial temperature of 130°C. As shown in Figure 9, there is a significant thickness reduction at some locations with small features. The maximum thickness-reducing ratio reaches to about 70%, which will result in serious defects for the formed part. The thickness distribution contour of the formed part under an initial temperature of 170°C is shown in Figure 10. The maximum thickness-reducing ratio improves to about 38%. This demonstrates that the heating temperature is a critical processing parameter for the thermoforming of the PP wood composite sheet.Figure 9
Thickness distribution with forming temperate of 130°C.Figure 10
Thickness distribution with forming temperate of 170°C.


## 5. Conclusions

Based on thermocompression test results of PP wood fiber composites under different temperature and strain rate, a thermal-mechanical viscoelastic model is developed for the PP wood composite with a phenomenological approach by using the Maxwell model. The material parameters in the model are either simplified as constants or assumed to be temperature-related or strain-rate-related only. With those reasonable simplifications and assumptions, the effect of temperature is decoupled with that of the strain rate in the material modelling. This greatly facilitates the material characterizing process. The developed viscoelastic model is validated by comparing numerical results with experimental data under a strain rate differing from those used in the material characterizing process. This study provides a basis for the thermoforming simulation and processing design for the PP wood fiber composite.



---

*Source: 793617-2012-10-14.xml*