Next Article in Journal
Structural Investigation of Poly(ethylene furanoate) Polymorphs
Next Article in Special Issue
Use of Irradiated Polymers after Their Lifetime Period
Previous Article in Journal
60Co γ-ray Irradiation Crosslinking of Chitosan/Graphene Oxide Composite Film: Swelling, Thermal Stability, Mechanical, and Antibacterial Properties
Previous Article in Special Issue
The Effect of Irradiation on Mechanical and Thermal Properties of Selected Types of Polymers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computer Modeling for Single-Screw Extrusion of Wood–Plastic Composites

by
Krzysztof Wilczyński
*,
Kamila Buziak
,
Krzysztof J. Wilczyński
,
Adrian Lewandowski
and
Andrzej Nastaj
Polymer Processing Department, Faculty of Production Engineering, Warsaw University of Technology, 02-524 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Polymers 2018, 10(3), 295; https://doi.org/10.3390/polym10030295
Submission received: 9 February 2018 / Revised: 3 March 2018 / Accepted: 6 March 2018 / Published: 9 March 2018
(This article belongs to the Special Issue Model-Based Polymer Processing)

Abstract

:
Experimental and theoretical studies have been performed on the single-screw extrusion of wood–plastics composites. Experimental research has been carried out on the flow and melting of polypropylene (PP)-based composites with different wood flour (WF) content in the single-screw extruder. Based on these experimental observations, elementary models of the process phenomena have been proposed and a global model of the process has been developed. This global computer model includes the models of solid conveying, melting dependent on the wood flour content, melt flow in the screw, and melt flow in the die. 3-D non-Newtonian finite element method (FEM) screw pumping characteristics have been applied to model the melt flow in the metering section of the screw. The model predicts the extrusion output, pressure and temperature profiles, melting profile, and power consumption. The model has been validated experimentally. An effect of material slip on the extruder operation has been discussed including both slipping in the screw/barrel surfaces and in the extruding die.

Graphical Abstract

1. Introduction

Wood–plastic composites (WPC) are made up of thermoplastics and wood flour or fiber. Both are polymer-based materials, but synthetic and natural, respectively. Thermoplastics are high-molecular-weight materials forming a continuous matrix that surrounds the wood component. These matrix materials are usually low-cost common thermoplastics that flow easily when heated, which translates to processing flexibility when wood flour or fiber is composed with them. Thermoplastics shrink and swell, but absorb only a small amount of moisture, and can be an efficient barrier to moisture penetrating into the composite. Wood flour or fiber is composed of natural polymers such as lignin, cellulose, etc., and has very different properties from the synthetic polymers. Wood is less expensive and stronger than synthetic polymers, and it may be a useful filler or reinforcement. The relatively low density of these natural fillers is also an important advantage. Wood does not shrink and swell substantially with temperature, but easily absorbs moisture which leads to biodegradation when not protected.
WPC exhibit high durability, relatively high stiffness and strength, and low price compared with other competitive materials. They are resistant to water and weather; and can be used for various outdoor applications. The WPC market has increased considerably in the last few years which is mainly driven by the growing building and construction market and by the automative industry [1]. The most commonly used thermoplastics for WPC are high-density polyethylene (HDPE), polypropylene (PP), and polyvinyl chloride (PVC) [1,2,3,4,5,6].
Worldwide WPC production has increased from 2.43 million tonnes in 2012 to 3.83 million tonnes in 2015, and it keeps growing. Although North America was the world’s production leader in 2012 (1.1 million t), ahead of China (900,000 t) and Europe (260,000 t), it was expected that China (with 1.8 million t) would have overtaken North America (1.4 million t) in 2015. European production was to increase by around 10% per year and reach 350,000 tonnes in 2015 [1].
Currently, due to extensive use of WPC in the form of profiles, the basic technology for processing these materials is extrusion. Extrusion is the most important technology in the polymer processing industry. It is broadly used for profile production (pipes, films, sheets, etc.), as well as for specialty polymer processing operations, such as compounding, pelletizing, mixing, devolatilization, chemical reactions, etc.
Extrusion of WPC differs considerably from plastics extrusion; this results from the different thermal and rheological properties of these highly filled materials, different structure, etc. Very limited research has been carried out on the rheology and extrusion of wood–plastic composites. Fundamental research has been performed by Xiao and Tzoganakis [7,8,9,10] as well as by Li and Wolcott [11,12,13]. Important contributions to this subject have been also delivered by Oksman Niska and Sain [5], Vlachopoulos et al. [14,15,16], and Zolfaghari et al. [17].
Now, it is known that WPC are pseudoplastic and viscoelastic materials. The viscosity of these materials decreases with an increase of shear rate and temperature; and increases with an increase of filler content. A decrease of elasticity with an increase of filler content has been also observed (Li and Wolcott [11,12,13], Mohanty et al. [3], Klyosov [4], Oksman Niska and Sain [5]). WPC exhibit yield stress and wall slip during extrusion, as reported by Xiao and Tzoganakis [7], Li and Wolcott [11], and Hristov et al. [14]. The wall slip velocity is dependent on the wood filler content and shear rate. With an increase of shear rate the slip velocity increases, which leads to the plug flow (Vlachopoulos and Hristov [15]). Also, an increase in wood content promotes the plug flow (Zolfaghari et al. [17]).
Recently, an interesting study has been performed by Laufer et al. [18] for WPC with different matrices and different types and contents of wood fillers. It has been concluded that shear thinning behavior of the WPC melts can be well described by the power law, and the consistency index of the power law well correlates with the wood volume fraction. The concept of equivalent inner shear rate has been developed, and on this basis the shear thinning behavior of WPC can be estimated with good accuracy, taking into account the wood volume fraction.
Designing the extrusion process of wood–plastic composites requires some knowledge about the flow mechanism of these highly filled materials in the extruder and in the die. Unfortunately, current theories do not predict the extrusion process course of these composites, e.g., material melting, pressure, and temperature distribution [8,10,16].
A few melting studies have been performed by Xiao and Tzoganakis [7,8,10] for extrusion of high-density polyethylene (HDPE) composite, and experimental data has been compared to simulations using commercial software. However, the predictions of pressure were not consistent with those measured.
Recently, Wilczyński et al. [19] have carried out extensive experimental studies on extrusion of a polypropylene (PP) composite. It was observed that material transport and melting were strongly dependent on the material composition and process parameters. The contiguous solids melting mechanism (CSM), known as a Tadmor mechanism [20], was not seen for composites of higher wood flour (WF) content. However, it was seen for composites of lower WF content of up to 50%. Several steps of the process were distinguished: first, compacting of the material in the solid conveying section close to the active flight and spreading to the rest of material; then, melting strongly dependent on the material composition; and finally shear or plug flow with possible slipping at the walls in the melt conveying section.
Up to the present, there is not any global model of WPC extrusion that would include solid conveying, melting, and melt flow. The state of the art of the extrusion of filled polymers has been recently discussed by Ariffin and Ahmad [21], and that of modeling of polymer extrusion by Wilczyński et al. [22]. In this paper, modeling of single-screw extrusion of wood–plastic composites is discussed, and an original computer model of the process is proposed and validated experimentally.
The basis for modeling the extrusion process is a knowledge of a material flow mechanism in the machine, i.e., solid transport, melting, and melt flow. The melt flow is quite well understood, but only for viscous flows. Viscoelastic flows are much less well understood. The flow of filled plastics that usually exhibit yield stress and wall slip is also poorly understood. Thus, global extrusion models are generally limited to viscous flows. To produce global models for wall-slipping materials it is necessary to develop models both for a plasticating unit and for an extrusion die [23,24,25,26].
Solids conveying and melting are much less understood and investigated, especially for advanced polymeric materials such as polyblends and composites. It is important to note that a melting model is always fundamental for the development of the global model of a process.
Melting in single-screw extruders has been widely studied and modeled. The first observations were reported by Maddock [27] and Street [28]. Their studies were based on the rapid cooling of an extruder, removing the screw and the polymer strip that was cross sectioned. This technique is called the “screw pulling-out technique”. In the melting section, a melt layer was seen along the barrel, and a melt pool close to the active flight. These observations were confirmed by many other researchers, e.g., by Wilczyński et al. [29] (Figure 1).
Tadmor [30,31] and Tadmor and Klein [20] first developed a model for the melting of polymers in single-screw extruders, and finally the first global model of the single-screw extrusion process. According to this model, a melt layer is formed between the barrel and solid, and molten polymer is scraped off by the transverse flow and accumulates at the active flight of the screw. The solid bed is gradually decreased by the combined effect of heating from the barrel and viscous dissipation in the melt. Several melting models were later developed and have been recently discussed by Altinkaynak et al. [32], and global extrusion models have been built up, e.g., by Agur and Vlachopoulos [33], Vincelette et al. [34], Potente et al. [35], and Wilczyński [36].
Contrary to the melting in single-screw extruders, the studies on melting in twin-screw extruders, both co-rotating and counter-rotating, appeared much later in the literature. Several melting models have been built up for co-rotating extruders, e.g., by Bawiskar and White [37,38], Sakai [39], Potente and Melish [40], Gogos et al. [41], and Vergnes et al. [42,43], and a few global models have been developed, e.g., by Bawiskar and White [44], Vergnes et al. [45] and Potente et al. [46].
Limited studies have been performed on melting in counter-rotating extruders, e.g., by Wilczyński and White [47], and Wilczyński et al. [48]. A melting model and the first global models have been developed by Wilczyński and White [49], Wilczyński et al. [50,51], and Jiang et al. [52].
The literature on starve-fed single-screw extrusion is extremely limited. Recent experimental studies performed by Wilczyński et al. [29] revealed the mechanism of melting in starve-fed single-screw extruders, and then the authors developed a mathematical model of melting of polymers in these machines [53], and finally a global model of the process [54].
Summarizing, reviews on the modeling of polymer extrusion have been recently presented by Wilczyński et al. [22], Malik et al. [24], Ilinca and Hetu [55], and Teixeira et al. [56]. A few basic monographs have also taken up this problem, e.g., by White and Potente [57], Tadmor and Gogos [58], and Rauwendaal [59].

2. Experimental

An instrumented single-screw extruder was used in the study, equipped with a barrel of diameter D = 45 mm and a length/diameter ratio (L/D) = 27. The all-purpose three-sectional screw was applied. It has feeding, compression, and metering sections with length/diameter ratios equal to (L/D)F = 10.78, (L/D)C = 7.11, and (L/D)M = 9.11. The compression ratio, which is defined as a ratio of the channel depth (HF) in the feeding section to the channel depth in the metering section (HM), i.e., CR = HF/HM, was equal to CR = 2.66 (HF = 8 mm, HM = 3 mm). The die for extruding rods of diameter D = 5 mm was applied. Pressure, temperature, and power consumption were measured in the study. A “screw pulling-out technique” was applied to investigate the material flow and melting mechanism. This technique is described elsewhere [29]. The screw speed was set at N = 20 rpm, N = 50 rpm, N = 80 rpm, and the temperature at TI = 160 °C, TII = 180 °C, TIII = 190 °C, TIV = 190 °C in the four barrel sections (of equal length) and TDIE = 180 °C in the die. All these parameters have been collected in Table 1 and Table 2.
Various polypropylene (PP)-based composites were used in the study (manufactured by Beologic). These were WPC of WF/PP compositions equal to 60/40, 50/50, 40/60, and 30/70. We used (i) WPC 60/40 PP ext 044/0-400, (ii) WPC 50/50 PP copo inj 4, (iii) WPC 70/30 PPC ext 101/0-400, and PP copo inj 2. The composites with WF/PP content equal to 40/60 and 30/70 were obtained by blending WPC 70/30 with PP in a single-screw extruder using a granulating die. A small portion of the material (3%) was colored with 2% Ramafin-Orange CEK071 dye-stuff to distinguish between solid and melt. The material was dried at 80 °C for 4–6 h.
Polypropylene has density 0.9 g/cm3, melt flow index MFI = 25 g/10 min (230 °C, 2.16 kg), and a melting point of about 160 °C. The WPC composites have density 1.1–1.3 g/cm3, bulk density 0.4–0.6 g/cm3, and mass flow rate index MFR = 3 g/10 min (200 °C, 5 kg).
Viscosity curves were determined for WPC 50/50 PP copo inj 4. Rheometric tests were carried out using a high-pressure capillary rheometer Rheograph 6000 (from Goettfert, Buchen, Germany) operating on the principle of constant shear rate. The measurements were performed in the temperature range 175–195 °C and the shear rate range 5–3500 s−1. Rabinowitsch and Bagley corrections were applied in rheological calculations. The viscosity characteristics are depicted in Figure 2. The viscosity curves present a typical non-Newtonian pseudoplastic behavior with decreasing viscosity with an increase of shear rate.
The mass flow rate index (MFR) was additionally determined for this composite using a Melt Indexer MI-2 (from Goettfert). The tests were carried out at 190 °C, and two loads were applied—High Load (HL = 21.6 kg) and Low Load (LL = 10 kg). We obtained MFRHL= 38 g/10 min, and MFRLL = 4 g/10 min.

3. Results

The main goal of experimental investigations was to learn the mechanism of flow and melting of the studied composites as the physical basis for modeling. Sample results of the tests are presented in Figure 3, Figure 4, Figure 5 and Figure 6.
It was observed that in the vicinity of the hopper the material forms a loose structure which gradually thickens starting from the active flight of the screw, and the material starts to melt after being compacted. The effect of wood flour (WF) content on the material structure and flow in the screw channel is shown in Figure 3.
It is clearly visible that with a lower WF content there is a classical melting mechanism with accumulation of the molten material (melt pool) at the active flight of the screw. With higher WF content this mechanism does not occur. This can be explained by the small amount of molten material being not enough to form the melt pool, and also by insufficient pressure in this area which is necessary to initiate forming the melt pool. It can be predicted that, in this case, the molten material penetrates into the unmolten layer and melting has a one-dimensional character, proceeding from the barrel to the core of the screw. A similar mechanism was observed when extruding plastics in powdered form [60], as well as in the case of our previous tests [19]. After melting was completed, shear melt flow or plug flow with slipping at the walls occurred. It was observed that an increase in wood content promoted existence of the plug flow.
The influence of the screw rotational speed on the material flow is shown in Figure 4. It can be seen here that with an increase of the rotational speed of the screw, melting begins later, i.e., further away from the hopper, and closer to the die. Melting also seems to be slower since the residence time of the material in the extruder is shorter.
The classical mechanism of melting of wood–plastic composites (WPC) in the single-screw extruder is shown in Figure 5, in the side view of the screw and in the cross section.
Figure 6 presents cross-sectional views of composites with different wood flour (WF) content (in the melting section) compared with cross-sectional views of neat polypropylene (PP). The similarity is obvious, although with an increase in the WF content the melt pool decreases.
During the tests, the pressure and flow rate were measured; these were compared with the simulation data and are presented in Table 3. The influence of the rotational screw speed on the material flow rate (i.e., the extrusion throughput) and extrusion pressure (i.e., the die pressure) was quite obvious. As the screw speed increased, the flow rate and the die pressure increased. However, this dependence was not proportional but rather weaker.

4. Modeling

On the basis of the performed research, a mechanism for flow and melting of wood composites in a single-screw extruder was proposed; this mechanism is depicted in Figure 7.
An extrusion process is the result of interaction between the extruder and the die. Extrusion modeling must take this interaction into account. It is not possible to consider the phenomena occurring in the extruder in isolation from the flow of material in the die. The flow of material in the die is relatively simple and well described (for viscous flows), although not in the case of the wood composites tested in this work. The flow of material in the extruder is much more complex, since it involves transporting the material in the solid state, melting, and flow of the molten material. The description of this flow is difficult even in the case of neat polymers, and especially difficult and very poorly understood in the case of composite materials.
In the conventional extrusion process (without dosing), the flow rate of the material is not determined by the extruder operator, but it results from the cooperation of the extruder–die system. Conditions of this cooperation are determined by the operating point of the extruder, which determines the extrusion throughput (polymer flow rate), and the extrusion pressure (die pressure). The basis of the calculation algorithm is to solve the problem of determining the flow rate and pressure distribution in the plasticating unit and the die. This problem can only be solved using the iterative calculation procedure in which the compliance of the pressure increase in the plasticating unit is tested with the pressure drop in the die. This is classical approach presented, e.g., in [36].
The global extrusion model is a series connection of elementary models describing the material flow in the plasticating unit of the extruder; that is, models describing transport of solid material, melting, and flow of molten material, and then the melt flow in the die.
Since the flow rate during extrusion is constant, the process can be considered as a series of separate elementary segments having the same mass flow rate (G), and locally constant material and process parameters. This is the lumped parameter approach which is sufficient for engineering computations. The process parameters, e.g., pressure, temperaturę, or solid content, at the beginning of each segment are equal to the process parameters at the end of the previous segment, i.e.,
pi_in (z) = pi-1_out (z)
where pi_in (z) is the process parameter (e.g., pressure or temperature) at the beginning of the ith segment, pi-1_out (z) is the process parameter at the end of the (i − 1)th segment, and z is the location of the segment Δz along the screw channel length.
The global extrusion model presented here contains some elements of our previous developments in the range of solid conveying [36], and introduces new solutions for melting and melt flow (applying 3-D screw pumping characteristics).
In summary, the model consists of the following elementary models:
Solid conveying, where delay in melting is allowed;
Melting 1-D or 2-D, depending on the wood flour content;
Melt conveying in the plasticating unit using 3-D finite element method (FEM)-based screw pumping characteristics;
Melt flow in the die 1-D or 3-D FEM.
The scheme of computations is depicted in Figure 8. The computations start assuming some initial value of the flow rate, e.g., equal to the drag flow rate. Then the process is simulated with dependence on the wood flour content (WF) according to 1-D or 2-D melting model. Next, 3-D screw pumping characteristics are applied for modeling of the melt conveyance, and the pressure at the screw exit is finally computed. This pressure is compared with the pressure drop in the die (in general computed as 1-D flow), and convergence is checked. Depending on the result (negative or positive) the flow rate is changed (increased or decreased) and computations are repeated until convergence is reached.
A one-dimensional melting process (1-D model) can be described on the basis of the heat balance on the interface of the solid and melt, and of the mass balance in the melt and solid. In order to make such a balance, one should determine the velocity and temperature profile in the melt film and the temperature profile in the solid.
According to the heat balance, the heat flux used to melt the material is equal to the difference between the heat flux from the film into the melt/solid interface and the heat flux from the interface into the solid bed. This balance does not differ from the heat balance of the two-dimensional classic approach (2-D model) [36,58].
However, the mass balance in the melt and solid differs since the solid bed width is constant in this case and equal to the screw channel width. The material is melted and the volume of melt increases; however, the melt pool at the active flight of the screw is not formed. An increase in melt volume, i.e., in film thickness, decreases the dissipation since the shear rate decreases when the film thickness increases. This causes slower melting.
Global modeling of the extrusion process requires iterative, multiple calculations (sometimes hundreds of iterations). Therefore, the use of time-consuming calculations using the finite element method (FEM) is not used here. However, you can develop the so-called flow characteristics based on FEM calculations, and then implement them in the form of regression models to the global process model. A good accuracy of the calculation is then obtained with reasonable calculation time.
These flow characteristics, i.e., screw pumping characteristics, are defined as functions of dimensionless flow rate and dimensionless pressure gradient:
Q* = f (Δp*)
where Q* is the dimensionless flow rate, and Δp* is the dimensionless pressure gradient.
This concept has been recently applied by the authors for starve-fed and flood-fed single-screw extrusion with conventional screws [54] and with mixing screws [61,62], as well as for counter-rotating twin-screw extrusion [63,64].
In this study, we model single-screw extrusion with conventional screws, and in this case the screw pumping characteristics are defined as
Q* = 2Q/WHVbz,
Δp* = Hn+1sin φ/6mVbzn·(Δpc/Lf),
where Q is the flow rate, W is the screw channel width, H is the screw channel depth, Vbz = πDbNcos φ is the z-component of the circumferential velocity at the barrel, Db is the barrel diameter, N is the screw speed, φ is the helix angle, Δpc is the pressure change, Lf is the screw length of the pressure change, m is the consistency coefficient, and n is the power-law exponent.
These characteristics are depicted in Figure 9a for different values of power-law exponent.

5. Simulations

Simulations were performed for the range of experimental material and processing data. These parameters are collected in Table 1 and Table 2.
The viscous flow properties were modeled using the Klein rheological equation [65]:
ln η = A 0 + A 1 ln γ ˙ + A 11 ln 2 γ ˙ + A 12 ln γ ˙ T + A 2 T + A 22 T 2
where η is the viscosity, Pa·s; γ ˙ is the shear rate, s−1; T is the temperaturę, °C; and A0, A1, A11, A12, A2, A22 are the model parameters (A0 = 12.469780638, A1 = −0.8345507, A11 = −0.017832191, A12 = 0.001331159, A2 = −0.008413991, A22 = −0.000025745).
Simulation results (dimensionless) are depicted in Figure 10, Figure 11, Figure 12 and Figure 13. Overall extrusion characteristics are depicted in Figure 10. These include pressure and temperature profiles, melting profile, i.e., solid bed profile (SBP), and power consumption profile, as well as the fill factor profile for starve-fed extrusion which determines the extent of filling of the screw.
We simulated the process using two different models of melting—the classical two-dimensional model and a one-dimensional model. The results are depicted in Figure 11, Figure 12 and Figure 13. It is clearly seen that melting according to the 1-D model is substantially slower; however, the impact on the pressure profile is not significant.
We validated pressure computations, and the results of validation are depicted in Figure 14. The pressure profile is well predicted; however, in general, it is overestimated.

6. Modeling of Slip Effects

Wood–plastic composites are wall-slipping materials and this fact has not been considered in modeling in the previous sections. This may be one of the causes of discrepancies between simulations and experimentations. Detailed study on the problem of global modeling of the extrusion process with slip effects will be presented elsewhere [66]. In this paper, some crucial aspects of this issue are presented.
ANSYS Polyflow CFD software [67] has been used for simulations of the extrusion process with slip effects, and the generalized Navier’s law has been applied for slip modeling. This law can be expressed in the following way:
f s = F slip ( V wall V s ) | V s V wall | e slip 1
where fs is the shear stress, Vs is the tangential velocity of the fluid, Vwall is the tangential velocity of the wall, and Fslip and eslip are the material parameters (Fslip = 0 denotes full slip conditions, i.e., fs = 0, and eslip = 1 denotes linearity of the law).
The constitutive equation of power-law fluid has been used in this study for flow modeling, and is expressed as
τ   = m γ ˙ n
where τ is the shear stress, γ ˙ is the shear rate, m is the consistency coefficient and n is the power-law exponent.
The parameters of the model were m = 30,111.9 Pa·sn and flow index n = 0.249.
Fully three-dimensional non-Newtonian simulations have been performed of the flow with slip effects both in the extruder (on the screw and barrel surfaces) and in the die. This allowed global study of operation of the extrusion process of wall-slipping polymers. Examples are depicted in Figure 15 and Figure 16. It is clearly seen that the velocity profile drastically changes and pressure substantially drops when slipping increases, both in the die and in the extruder.
Screw pumping characteristics at different slip conditions (with slip and without slip) are shown in Figure 9 for Newtonian and non-Newtonian cases. Decreasing the power-law exponent, which increases the non-Newtonian character of the flow, importantly reduces the screw pumping capacity (for positive pressure gradient). It is also clearly seen that slipping decreases the pressure gradient both when positive or negative.
Extrusion is a process of pushing of the polymer through a forming die, and can be considered as a series connection of the extruder and the die. Since the die exerts a resistance to flow, pressure is required to push the polymer through the die. One should know that the die pressure is caused by the die, and not by the extruder. The extruder simply produces pressure to push the polymer through the die. An extrusion process can be described by the extruder operating characteristics, which are determined by the screw characteristics and the die characteristics (Figure 17). Actual extrusion operating conditions are defined by the extruder operating point which results from an intersection of the screw characteristics and the die characteristics. This point determines the flow rate (throughput) and the extrusion pressure (die pressure).
Figure 17 shows the effect of slip conditions determined by the Fslip parameter of Navier’s law (6) on the extruder operating point. The starting point (1) has been obtained for low slip at the screw and no slip at the die. When slip at the die is allowed, and the slip at the screw increases, the operating point moves to lower flow rate and lower pressure, i.e., to the point (2). It can be concluded that slip at the screw and at the die have an important impact on the location of the extruder operating point—that is, on the flow rate and extrusion pressure. When slip at the screw increases, the flow rate and pressure decrease; when slip at the die increases, the flow rate increases and the pressure decreases.

7. Conclusions

Experimental research has shown that flow and melting of polypropylene (PP)-based wood composites are considerably different from those of neat polypropylene. In the initial part of the screw, the processed material forms a loose structure that gradually thickens and, after compaction, begins to melt. Melting depends on the content of wood flour (WF) in the material. In a simplification, it can be concluded that in the case of lower WF content classical 2-D melting takes place, while in the case of higher WF content 1-D melting is observed without a melt pool at the active flight of the screw. In this case, an increase in the melt volume, i.e., in the film thickness, decreases the dissipation since the shear rate decreases with an increase of the film thickness, and this causes slower melting.
A global computer model of single-screw extrusion of wood–plastic composites has been developed that describes solid conveying, melting dependent on the wood flour content, melt flow in the screw, and melt flow in the die. The model predicts the extrusion output, pressure and temperature profiles, melting profile, and power consumption. The model has been successfully validated experimentally; however, the pressure was overpredicted.
It has been concluded that slip at the screw and the die have an important impact on extruder operation. When slip at the screw/barrel increases, the extrusion output and pressure decrease, while when slip at the die increases, the extrusion output increases and the pressure decreases. We suggest an additional modeling study taking into account possible yield stress of the processed material.

Acknowledgments

The authors would like to acknowledge support from National Science Center, Poland (DEC-2015/19/B/ST8/00948).

Author Contributions

Krzysztof Wilczyński wrote the paper, Kamila Buziak elaborated the results of experimentation and simulation, Krzysztof J. Wilczyński performed rheological studies, Adrian Lewandowski simulated the 3-D flow using Polyflow software, Andrzej Nastaj implemented new model elements into the code of the global program. The authors performed experiments and modeled the process.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nova-Institut GmbH. New Market and Trend Report Published. Available online: http://news.bio-based.eu/new-market-and-trend-report-published (accessed on 21 November 2017).
  2. Bledzki, A.K.; Reihmane, S.; Gassan, J. Thermoplastics Reinforced with Wood Fillers: A Literature Review. Polym. Plast. Technol. Eng. 1998, 37, 451–468. [Google Scholar] [CrossRef]
  3. Mohanty, A.K.; Misra, M.; Drzal, L.T. Natural Fibers, Biopolymers, and Biocomposites; Taylor & Francis: Boca Raton, FL, USA, 2005; ISBN 978-084-931-741-5. [Google Scholar]
  4. Klyosov, A.A. Wood-Plastic Composites; John Wiley & Sons: Hoboken, NJ, USA, 2007; ISBN 978-047-014-891-4. [Google Scholar]
  5. Oksman Niska, K.; Sain, M. Wood-Polymer Composites; Woodhead Publishing Ltd.: Cambridge, UK, 2008; ISBN 978-184-569-272-8. [Google Scholar]
  6. Zajchowski, S.; Ryszkowska, J. Kompozyty polimerowo-drzewne—Charakterystyka ogólna oraz ich otrzymywanie z materiałów odpadowych. Polimery 2009, 10, 674–682. [Google Scholar]
  7. Xiao, K.; Tzoganakis, C. Rheological Properties and Their Influence on Extrusion Characteristics of HDPE-Wood Composite Resins. In SPE ANTEK; Society of Plastics Engineers: Houston, TX, USA, 2002; Volume 39, pp. 252–257. [Google Scholar]
  8. Xiao, K.; Tzoganakis, C. Extrusion Characteristic of HDPE-Wood Composites. In SPE ANTEK; Society of Plastics Engineers: Houston, TX, USA, 2003; Volume 39, pp. 197–201. [Google Scholar]
  9. Xiao, K.; Tzoganakis, C. Rheological Properties of HDPE-Wood Composites. In SPE ANTEC; Society of Plastics Engineers: Houston, TX, USA, 2003; Volume 39, pp. 975–979. [Google Scholar]
  10. Xiao, K.; Tzoganakis, C. Extrusion of HDPE-Wood Blends. In SPE ANTEK; Society of Plastics Engineers: Houston, TX, USA, 2004; Volume 40, pp. 448–453. [Google Scholar]
  11. Li, T.Q.; Wolcott, M.P. Rheology of wood plastics melt. Part 1. Capillary rheometry of HDPE filled with maple. Polym. Eng. Sci. 2005, 45, 549–559. [Google Scholar] [CrossRef]
  12. Li, T.Q.; Wolcott, M.P. Rheology of wood plastics melt, part 2: Effects of lubricating systems in HDPE/maple composites. Polym. Eng. Sci. 2006, 46, 464–473. [Google Scholar] [CrossRef]
  13. Li, T.Q.; Wolcott, M.P. Rheology of wood plastics melt, part 3: Nonlinear nature of the flow. Polym. Eng. Sci. 2006, 46, 114–121. [Google Scholar] [CrossRef]
  14. Hristov, V.; Takács, E.; Vlachopoulos, J. Surface tearing and wall slip phenomena in extrusion of highly filled HDPE/wood flour composites. Polym. Eng. Sci. 2006, 46, 1204–1214. [Google Scholar] [CrossRef]
  15. Vlachopoulos, J.; Hristov, V. Rheology of Wood Polymer Composites. In Proceedings of the Regional PPS Meeting, Yamagata, Japan, 2–6 July 2006. [Google Scholar]
  16. Santi, C.R.; Hage, E.; Vlachopoulos, J.; Correa, C.A. Rheology and Processing of HDPE/Wood Flour Composites. Int. Polym. Process. 2009, 24, 346–353. [Google Scholar] [CrossRef]
  17. Zolfaghari, A.; Behravesh, A.H.; Shakouri, E.; Soury, E. Flow balancing in die design of wood flour/HDPE composite extrusion profiles with consideration of rheological effect. Polym. Eng. Sci. 2010, 50, 543–549. [Google Scholar] [CrossRef]
  18. Laufer, N.; Hansmann, H.; Koch, M. Rheological Characterisation of the Flow Behaviour of Wood Plastic Composites in Consideration of Different Volume Fractions of Wood. J. Phys. Conf. Ser. 2017, 790, 012017. [Google Scholar] [CrossRef]
  19. Wilczyński, K.; Nastaj, A.; Lewandowski, A.; Wilczyński, K.J.; Buziak, K. Experimental Study for Extrusion of Polypropylene/Wood Flour Composites. Int. Polym. Process. 2015, 30, 113–120. [Google Scholar] [CrossRef]
  20. Tadmor, Z.; Klein, I. Engineering Principles of Plasticating Extrusion; Van Nostrand Reinhold Co.: New York, NY, USA, 1970; ISBN 978-044-215-635-0. [Google Scholar]
  21. Ariffin, A.; Ahmad, M.S.B. Review: Single Screw Extruder in Particulate Filler Composite. Polym. Plast. Technol. Eng. 2011, 50, 395–403. [Google Scholar] [CrossRef]
  22. Wilczyński, K.; Nastaj, A.; Lewandowski, A.; Wilczyński, K.J. Multipurpose Computer Model for Screw Processing of Plastics. Polym. Plast. Technol. Eng. 2012, 51, 626–633. [Google Scholar] [CrossRef]
  23. Potente, H.; Ridder, H.; Cunha, R.V. Global Concept for Describing and Investigation of Wall Slip Effects in the Extrusion Process. Macromol. Mater. Eng. 2002, 287, 836–842. [Google Scholar] [CrossRef]
  24. Malik, M.; Kalyon, D.M.; Golba, J.C. Simulation of Co-Rotating Twin Screw Extrusion Process Subject to Pressure-Dependent Wall Slip at Barrel and Screw Surfaces: 3D FEM Analysis for Combinations of Forward- and Reverse-Conveying Screw Elements. Int. Polym. Process. 2014, 29, 51–62. [Google Scholar] [CrossRef]
  25. Gooneie, A.; Schuschnigg, S.; Duretek, I.; Holzer, C. Numerical simulations of the flow of wood polypropylene composites with wall slipping in a profile die: The significance of material data. In AIP Conference Proceedings 1664; American Institute of Physics: New York, NY, USA, 2015; p. 050014. [Google Scholar]
  26. Duretek, I.; Schuschnigg, S.; Gooneie, A.; Langecker, G.R.; Holzer, C. Rheological properties of wood polymer composites and their role in extrusion. J. Phys. Conf. Ser. 2015, 602, 012014. [Google Scholar] [CrossRef]
  27. Maddock, B.H. A Visual Analysis of Flow and Mixing in Extruder Screws. Soc. Plast. Eng. 1959, 15, 383–389. [Google Scholar]
  28. Street, L.F. Plastifying Extrusion. Intern. Plast. Eng. 1959, 1, 289–296. [Google Scholar]
  29. Wilczyński, K.; Lewandowski, A.; Wilczyński, K.J. Experimental study for starve-fed single screw extrusion of thermoplastics. Polym. Eng. Sci. 2012, 52, 1258–1270. [Google Scholar] [CrossRef]
  30. Tadmor, Z. Fundamentals of plasticating extrusion. I. A theoretical model for melting. Polym. Eng. Sci. 1966, 6, 185–190. [Google Scholar] [CrossRef]
  31. Tadmor, Z.; Duvdevani, I.; Klein, I. Melting in plasticating extuders theory and experiments. Polym. Eng. Sci. 1967, 7, 198–217. [Google Scholar] [CrossRef]
  32. Altınkaynak, A.; Gupta, M.; Spalding, M.A.; Crabtree, S.L. Melting in a Single Screw Extruder: Experiments and 3D Finite Element Simulations. Int. Polym. Process. 2011, 26, 182–196. [Google Scholar] [CrossRef]
  33. Agur, E.E.; Vlachopoulos, J. Numerical simulation of a single-screw plasticating extruder. Polym. Eng. Sci. 1982, 22, 1084–1094. [Google Scholar] [CrossRef]
  34. Vincelette, A.R.; Guerrero, C.S.; Carreau, P.J.; Lafleur, P.G. A Model for Single-screw Plasticating Extruders. Int. Polym. Process. 1989, 4, 232–241. [Google Scholar] [CrossRef]
  35. Potente, H.; Hanhart, W.; Schöppner, V. Potential Applications for Computer-aided Extruder Design. Int. Polym. Process. 1993, 8, 335–344. [Google Scholar] [CrossRef]
  36. Wilczyński, K. Single Screw Extrusion Model for Plasticating Extruders. Polym. Plast. Technol. Eng. 1999, 38, 581–608. [Google Scholar] [CrossRef]
  37. Bawiskar, S.; White, J.L. Solids Conveying and Melting in a Starve Fed Self-wiping Co-rotating Twin Screw Extruder. Int. Polym. Process. 1995, 10, 105–110. [Google Scholar] [CrossRef]
  38. Bawiskar, S.; White, J.L. Melting model for modular self wiping co-rotating twin screw extruders. Polym. Eng. Sci. 1998, 38, 727–740. [Google Scholar] [CrossRef]
  39. Sakai, T. The development of on-line techniques and novel processing systems for the monitoring and handling of the evolution of microstructure in nonreactive and reactive polymer systems. Adv. Polym. Technol. 1995, 14, 277–290. [Google Scholar] [CrossRef]
  40. Potente, H.; Melisch, U. Theoretical and Experimental Investigations of the Melting of Pellets in Co-Rotating Twin-Screw Extruders. Int. Polym. Process. 1996, 11, 101–108. [Google Scholar] [CrossRef]
  41. Gogos, C.C.; Tadmor, Z.; Kim, M.H. Melting phenomena and mechanisms in polymer processing equipment. Adv. Polym. Technol. 1998, 17, 285–305. [Google Scholar] [CrossRef]
  42. Vergnes, B.; Delacour, M.L.; Souveton, G.; Bouvier, J.M. A Study of Polymer Melting in a Co-Rotating Twin Screw Extruder. In Proceedings of the International Polymer Processing Society Annual Meeting PPS-15’s, Hertogenbosch, The Netherlands, 31 May–4 June 1999. [Google Scholar]
  43. Vergnes, B.; Souveton, G.; Delacour, M.L.; Ainser, A. Experimental and Theoretical Study of Polymer Melting in a Co-rotating Twin Screw Extruder. Int. Polym. Process. 2001, 16, 351–362. [Google Scholar] [CrossRef]
  44. Bawiskar, S.; White, J.L. A Composite Model for Solid Conveying, Melting, Pressure and Fill Factor Profiles in Modular Co. -Rotating Twin Screw Extruders. Int. Polym. Process. 1997, 12, 331–340. [Google Scholar] [CrossRef]
  45. Vergnes, B.; Valle, G.D.; Delamare, L. A global computer software for polymer flows in corotating twin screw extruders. Polym. Eng. Sci. 1998, 38, 1781–1792. [Google Scholar] [CrossRef]
  46. Potente, H.; Bastian, M.; Flecke, J. Design of a compounding extruder by means of the SIGMA simulation software. Adv. Polym. Technol. 1999, 18, 147–170. [Google Scholar] [CrossRef]
  47. Wilczyński, K.; White, J.L. Experimental Study of Melting in an Intermeshing Counter-Rotating Twin Screw Extruder. Int. Polym. Process. 2001, 16, 257–262. [Google Scholar] [CrossRef]
  48. Wilczyński, K.; Lewandowski, A.; Wilczyński, K.J. Experimental study of melting of LDPE/PS polyblend in an intermeshing counter-rotating twin screw extruder. Polym. Eng. Sci. 2012, 52, 449–458. [Google Scholar] [CrossRef]
  49. Wilczyński, K.; White, J.L. Melting model for intermeshing counter-rotating twin-screw extruders. Polym. Eng. Sci. 2003, 43, 1715–1726. [Google Scholar] [CrossRef]
  50. Wilczyński, K.; Jiang, Q.; White, J.L. A Composite Model for Melting, Pressure and Fill Factor Profiles in a Metered Fed Closely Intermeshing Counter-rotating Twin Screw Extruder. Int. Polym. Process. 2007, 22, 198–203. [Google Scholar] [CrossRef]
  51. Wilczyński, K.; Nastaj, A.; Lewandowski, A.; Wilczyński, K.J. Modelowanie przepływu tworzyw w procesie wytłaczania dwuślimakowego przeciwbieżnego, Cz.2. Badania symulacyjne i doświadczalne. Polimery 2011, 56, 45–50. [Google Scholar]
  52. Jiang, Q.; White, J.L.; Yang, J. A Global Model for Closely Intermeshing Counter-rotating Twin Screw Extruders with Flood Feeding. Int. Polym. Process. 2010, 25, 223–235. [Google Scholar] [CrossRef]
  53. Wilczyński, K.; Nastaj, A.; Wilczyński, K.J. Melting Model for Starve Fed Single Screw Extrusion of Thermoplastics. Int. Polym. Process. 2013, 28, 34–42. [Google Scholar] [CrossRef]
  54. Wilczyński, K.J.; Nastaj, A.; Lewandowski, A.; Wilczyński, K. A composite model for starve fed single screw extrusion of thermoplastics. Polym. Eng. Sci. 2014, 54, 2362–2374. [Google Scholar] [CrossRef]
  55. Ilinca, F.; Hétu, J.F. Three-dimensional Finite Element Solution of the Flow in Single and Twin-Screw Extruders. Int. Polym. Process. 2010, 25, 275–286. [Google Scholar] [CrossRef]
  56. Teixeira, C.; Gaspar-Cunha, A.; Covas, J.A. Flow and Heat Transfer Along the Length of a Co-rotating Twin Screw Extruder. Polym. Plast. Technol. Eng. 2012, 51, 1567–1577. [Google Scholar] [CrossRef]
  57. White, J.L.; Potente, H. Screw Extrusion; Carl Hanser Verlag: Munich, Germany, 2003; ISBN 344-619-624-2. [Google Scholar]
  58. Tadmor, Z.; Gogos, C. Principles of Polymer Processing; Wiley & Sons: New York, NY, USA, 2006; ISBN 047-138-770-3. [Google Scholar]
  59. Rauwendaal, C. Polymer Extrusion; Hanser: Munich, Germany, 2014; ISBN 978-156-990-516-6. [Google Scholar]
  60. Campbell, G.A.; Spalding, M.A. One Dimensional Melting in Single-Screw Extruders. In SPE ANTEC; Society of Plastics Engineers: Houston, TX, USA, 2011; Volume 47, pp. 1367–1373. [Google Scholar]
  61. Wilczyński, K.J.; Lewandowski, A.; Nastaj, A.; Wilczyński, K. Modeling for Starve Fed/Flood Fed Mixing Single-Screw Extruders. Int. Polym. Process. 2016, 31, 82–91. [Google Scholar] [CrossRef]
  62. Wilczyński, K.J.; Lewandowski, A.; Nastaj, A.; Wilczyński, K. A Global Model for Starve-Fed Nonconventional Single-Screw Extrusion of Thermoplastics. Adv. Polym. Technol. 2017, 36, 23–35. [Google Scholar] [CrossRef]
  63. Wilczyński, K.; Lewandowski, A. Study on the Polymer Melt Flow in a Closely Intermeshing Counter-Rotating Twin Screw Extruder. Int. Polym. Process. 2014, 29, 649–659. [Google Scholar] [CrossRef]
  64. Lewandowski, A.; Wilczyński, K.J.; Nastaj, A.; Wilczyński, K. A composite model for an intermeshing counter-rotating twin-screw extruder and its experimental verification. Polym. Eng. Sci. 2015, 55, 2838–2848. [Google Scholar] [CrossRef]
  65. Wilczyński, K.; Nastaj, A.; Buziak, K. Modelowanie procesu wytlaczania jednoślimakowego polimerowych kompozytów drzewnych przy dozowanym zasilaniu wytłaczarki. Polimery 2018, 63. in press. [Google Scholar]
  66. Lewandowski, A.; Wilczyński, K. Global Modeling for Single Screw Extrusion with Slip Effects. Int. Polym. Process. 2018, 33. in press. [Google Scholar]
  67. Ansys Polyflow v. 17. Available online: www.ansys.com (accessed on 7 March 2018).
Figure 1. Classical melting mechanism of polymers in the single-screw extruder on an example of extrusion of polypropylene (PP). Reprinted with permission from [29]. Copyright Wiley 2012.
Figure 1. Classical melting mechanism of polymers in the single-screw extruder on an example of extrusion of polypropylene (PP). Reprinted with permission from [29]. Copyright Wiley 2012.
Polymers 10 00295 g001
Figure 2. Viscosity curves of wood–plastic composite (WPC) 50/50 polypropylene (PP) copo inj 4.
Figure 2. Viscosity curves of wood–plastic composite (WPC) 50/50 polypropylene (PP) copo inj 4.
Polymers 10 00295 g002
Figure 3. The effect of wood flour (WF) content on the material flow in the screw channel (samples of the material of different wood content removed from the screw, obtained by extrusion at the screw speed N = 50 rpm): (a) WF = 30%, (b) WF = 40%, (c) WF = 50%, (d) WF = 60%.
Figure 3. The effect of wood flour (WF) content on the material flow in the screw channel (samples of the material of different wood content removed from the screw, obtained by extrusion at the screw speed N = 50 rpm): (a) WF = 30%, (b) WF = 40%, (c) WF = 50%, (d) WF = 60%.
Polymers 10 00295 g003
Figure 4. The effect of screw speed on the material flow in the screw channel (samples of the material of 40% WF content removed from the screw): (a) N = 20 rpm, (b) N = 50 rpm, (c) N = 80 rpm.
Figure 4. The effect of screw speed on the material flow in the screw channel (samples of the material of 40% WF content removed from the screw): (a) N = 20 rpm, (b) N = 50 rpm, (c) N = 80 rpm.
Polymers 10 00295 g004
Figure 5. Mechanism of melting of wood–plastic composite (WPC) in the single-screw extruder: (a) a view from the screw side (WF = 40%), (b) cross sections (WF = 30%); A—molten material, B—solid material.
Figure 5. Mechanism of melting of wood–plastic composite (WPC) in the single-screw extruder: (a) a view from the screw side (WF = 40%), (b) cross sections (WF = 30%); A—molten material, B—solid material.
Polymers 10 00295 g005
Figure 6. The effect of wood flour (WF) content (30–60%) on the melting mechanism (presented as cross sections) of wood–plastic composites (WPC) in the single-screw extruder compared with [29] melting of neat polypropylene (PP). Reprinted with permission from [29]. Copyright Wiley 2012.
Figure 6. The effect of wood flour (WF) content (30–60%) on the melting mechanism (presented as cross sections) of wood–plastic composites (WPC) in the single-screw extruder compared with [29] melting of neat polypropylene (PP). Reprinted with permission from [29]. Copyright Wiley 2012.
Polymers 10 00295 g006
Figure 7. Geometrical model of melting mechanism of wood–plastic composites (WPC) in the single-screw extruder with dependence on wood flour (WF) content: 1—barrel, 2—screw, A—molten material, B—solid material.
Figure 7. Geometrical model of melting mechanism of wood–plastic composites (WPC) in the single-screw extruder with dependence on wood flour (WF) content: 1—barrel, 2—screw, A—molten material, B—solid material.
Polymers 10 00295 g007
Figure 8. Scheme of computations: Q—flow rate, Q0—initial flow rate, p—pressure, Δpscrew—total pressure increase in the screw, Δpdie—total pressure drop in the die, T—temperaturę, SBP—melting profile (Solid Bed Profile), P—power consumption.
Figure 8. Scheme of computations: Q—flow rate, Q0—initial flow rate, p—pressure, Δpscrew—total pressure increase in the screw, Δpdie—total pressure drop in the die, T—temperaturę, SBP—melting profile (Solid Bed Profile), P—power consumption.
Polymers 10 00295 g008
Figure 9. Screw pumping characteristics: (a) without slip, (b) with slip (see Section 6).
Figure 9. Screw pumping characteristics: (a) without slip, (b) with slip (see Section 6).
Polymers 10 00295 g009
Figure 10. Dimensionless overall extrusion process characteristics, N = 80 rpm.
Figure 10. Dimensionless overall extrusion process characteristics, N = 80 rpm.
Polymers 10 00295 g010
Figure 11. Solid bed profile and pressure profile for 1-D and 2-D simulations, N = 20 rpm.
Figure 11. Solid bed profile and pressure profile for 1-D and 2-D simulations, N = 20 rpm.
Polymers 10 00295 g011
Figure 12. Solid bed profile and pressure profile for 1-D and 2-D simulations, N = 50 rpm.
Figure 12. Solid bed profile and pressure profile for 1-D and 2-D simulations, N = 50 rpm.
Polymers 10 00295 g012
Figure 13. Solid bed profile and pressure profile for 1-D and 2-D simulations, N = 80 rpm.
Figure 13. Solid bed profile and pressure profile for 1-D and 2-D simulations, N = 80 rpm.
Polymers 10 00295 g013
Figure 14. Validation of pressure computation.
Figure 14. Validation of pressure computation.
Polymers 10 00295 g014
Figure 15. Die flow simulations.
Figure 15. Die flow simulations.
Polymers 10 00295 g015
Figure 16. Screw flow simulations.
Figure 16. Screw flow simulations.
Polymers 10 00295 g016
Figure 17. Extruder operating characteristics: Q—flow rate, Δp—pressure change, S1, S2—screw characteristics, D1, D2—die characteristics.
Figure 17. Extruder operating characteristics: Q—flow rate, Δp—pressure change, S1, S2—screw characteristics, D1, D2—die characteristics.
Polymers 10 00295 g017
Table 1. Extruder screw geometry and dimensions.
Table 1. Extruder screw geometry and dimensions.
Square-Pitch Screw
Feed/compression/metering sections10.78/7.11/9.11 turns
Barrel inside diameter, Db45 mm
Screw pitch45 mm
Channel depth in feed section, HF8 mm
Channel depth in metering section, HM3 mm
Flight width5 mm
Table 2. Composite properties and processing variables.
Table 2. Composite properties and processing variables.
Extruder Operating Conditions
Barrel and die temperature profile160/180/190/190/180 °C
Frequency of screw rotation20, 50, and 80 rpm
Material and rheological properties
Density
  - bulk, ρbulk500 kg/m3
  - solid, ρs 1080 kg/m3
  - melt, ρm930 kg/m3
Polymer–barrel friction factor 0.25
Polymer–screw friction factor0.15
Heat capacity
  - melt, Cpm1850 J/(kgK)
  - solid, Cps1160 J/(kgK)
Thermal conductivity—melt, km0.243 W/(mK)
Heat of fusion, λ133,850 J/kg
Melting temperature, Tm125 °C
Viscosity, η (see Figure 2, Equation (5))
  - A012.469780638
  - A1−0.834550700
  - A11−0.017832191
  - A120.001331159
  - A2−0.008413991
  - A22−0.000025745
where η [Pa∙s], γ̇ [s−1], T [°C]
Table 3. Predicted versus measured process parameters.
Table 3. Predicted versus measured process parameters.
ExperimentSimulation (See Section 5)
Screw Speed
rpm
Throughput
kg/h
Die Pressure
MPa
Throughput
kg/h
Die Pressure
MPa
207.9106.610.7
5017.41116.812.1
8025.81327.512.6

Share and Cite

MDPI and ACS Style

Wilczyński, K.; Buziak, K.; Wilczyński, K.J.; Lewandowski, A.; Nastaj, A. Computer Modeling for Single-Screw Extrusion of Wood–Plastic Composites. Polymers 2018, 10, 295. https://doi.org/10.3390/polym10030295

AMA Style

Wilczyński K, Buziak K, Wilczyński KJ, Lewandowski A, Nastaj A. Computer Modeling for Single-Screw Extrusion of Wood–Plastic Composites. Polymers. 2018; 10(3):295. https://doi.org/10.3390/polym10030295

Chicago/Turabian Style

Wilczyński, Krzysztof, Kamila Buziak, Krzysztof J. Wilczyński, Adrian Lewandowski, and Andrzej Nastaj. 2018. "Computer Modeling for Single-Screw Extrusion of Wood–Plastic Composites" Polymers 10, no. 3: 295. https://doi.org/10.3390/polym10030295

APA Style

Wilczyński, K., Buziak, K., Wilczyński, K. J., Lewandowski, A., & Nastaj, A. (2018). Computer Modeling for Single-Screw Extrusion of Wood–Plastic Composites. Polymers, 10(3), 295. https://doi.org/10.3390/polym10030295

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop