Next Article in Journal
Synthesis and Properties of Thermally Self-Healing PET Based Linear Polyurethane Containing Diels–Alder Bonds
Next Article in Special Issue
Homogenisation of the Local Thermal Conductivity in Injection-Moulded Short Fibre Reinforced Composites
Previous Article in Journal
The Improvement of Sumatran Elephant (Elephas maximussumatranus) Dung Particleboard Characteristics Using Bamboo Layering
Previous Article in Special Issue
Statistical Analysis of the Mechanical Behavior of High-Performance Polymers: Weibull’s or Gaussian Distributions?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Enhanced Simulation of Infrared Heating of Thermoplastic Composites Prior to Forming under Consideration of Anisotropic Thermal Conductivity and Deconsolidation by Means of Novel Physical Material Models

by
Manuel Längauer
1,*,
Gernot Zitzenbacher
1,
Hannes Stadler
2 and
Christoph Hochenauer
3
1
School of Engineering, University of Applied Sciences Upper Austria, Stelzhamerstraße 23, 4600 Wels, Austria
2
Transfercenter fuer Kunststofftechnik GmbH, Franz-Fritsch-Straße 11, 4600 Wels, Austria
3
Institute of Thermal Engineering, Graz University of Technology, Inffeldgasse 25 b, 8010 Graz, Austria
*
Author to whom correspondence should be addressed.
Polymers 2022, 14(16), 3331; https://doi.org/10.3390/polym14163331
Submission received: 14 July 2022 / Revised: 10 August 2022 / Accepted: 13 August 2022 / Published: 16 August 2022

Abstract

:
In recent years, thermoplastic composites have found their place in large business sectors and are in direct rivalry to thermoset matrix composites. In order to ensure efficient and lean processes, process modeling gains ever-growing attention. This work shows the computational fluid dynamics (CFD)-modeling of a typical heating step in a thermoforming process of a thermoplastic composite sheet. When heating thermoplastic composites, the heat conduction proceeds anisotropic, and the sheets are subject to thermal deconsolidation when heated above the melting temperature of the polymer matrix adding to the anisotropic effect. These effects are neglected in known process models and this study shows the first successful attempt at introducing them into CFD-modeling of the heating of thermoplastic composite sheets. Thus, the simulation requires temperature dependent values for the anisotropic thermal conductivity and the coefficient of linear thermal expansion, which are calculated with novel physical models which were developed solely for this cause. This alters the behavior of an isotropic CFD-model and allows the successful validation via laboratory experiments using glass fiber reinforced polypropylene (PP/GF) sheets with embedded thermocouples to check the internal temperature distribution when the sheet is heated to the designated forming temperature in a composite thermoforming press. The incorporation of this newly developed process model reduces the error in the core temperature prediction from close to 70 °C to 3 °C at the forming temperature.

1. Introduction

Manufacturing of parts from thermoplastic matrix composites and the process modeling thereof has been in the scientific spotlight in the last couple of years. In particular, their high strength, simple processability, and high wear resistance make them very attractive [1,2,3]. The lightweight industry has a huge demand for recyclable, high performance materials to replace metal and thermoset materials. To meet those demands, high standards are imposed on material quality and lean, efficient processes. Many parts are manufactured in stamp-forming processes which are a kind of thermoforming process, which in turn is adopted from thermoplastic processing [4,5,6,7,8]. Therein semifinished sheets are transported to a heating station—mostly infrared heating—to be heated above the transition temperature of the polymeric matrix to make the material formable [9]. This is accompanied by thermal deconsolidation and the volume of the sheet increases significantly, which again coincides with a decrease in the thermal conductivity [10,11,12,13,14,15,16]. When the desired temperature is reached, the sheets are transported to a two-sided mold and formed into a three-dimensional part, cooled down, and demolded.
While much effort has been made on this matter already, especially the heating step in thermoforming processes is still a topic that requires significant work in both academic and industrial prospects. This step is vital, both for the process efficiency since it is the limiting factor when considering cycle time, and also for the quality of the finished product as the temperature and its distribution predetermines the forming behavior of the sheet.
Models for the anisotropic thermal conductivity are available and can be used for heating simulations. Ever wider spreading access to CFD methods has simplified the modeling of the heating behavior and has led to a rise in attempts thereof.
Brogan and Monaghan conducted extensive research on the heating behavior of carbon fiber reinforced polyether ether ketone (PEEK/CF) and proposed an attempt to calculate the temperature distribution when using quartz heaters [17].
One of the first visualized finite element simulations was performed by Johnson and Pickett [18] who investigated double curvature forming of carbon fiber reinforced polyetherimide (PEI/CF). In the same year, Sweeney et al. [19] conducted a study on the heating of PEEK/CF considering different heater setups and composite layups. They concluded that heater-to-composite distance had the most significant effect on the resulting temperature distribution.
Hwang and Hwang [20] used a right angle die to study the forming behavior of carbon fiber reinforced polyamide 6 (PA 6/CF) parts at different heating setups of the laminate. They studied the defects of overheated and underheated parts with micrographs and load tests and concluded that the forming temperature had a larger effect than holding time or heater distance.
Hsiao and Kikuchi (PEEK/CF), Abbassi et al. (carbon fiber reinforced polyphenylene sulfide (PPS/CF)), Chen et al. (PPS/CF) as well as Harrison, Gomes, and Curado-Correia (PP/GF) published numerical studies for the forming of already heated sheets [21,22,23,24]. This approach was also chosen by Stamopoulos and Di Ilio [25] and D’Emilia et al. [26] who both studied the forming of semi-spherical preheated composite parts with the use of numerical and analytical methods.
McCool et al. also worked with PPS/CF and described thermal deconsolidation in the heating stage when reporting the effect of forming and mold temperature on the properties of the formed sheets [27].
Nardi and Sinke [28] calculated the center temperature of a composite sheet analytically with the knowledge of the heater and the initial sheet temperature, the specific heat capacity, the density, and an isotropic thermal conductivity parameter. They considered the whole process and made estimations of the mechanical behavior of the final parts.
This work is a follow-up to preceding studies, in which a model for the anisotropic thermal conductivity of neat consolidated composite sheets was formulated [29]. This model was then refined for thermal deconsolidation, which occurs at temperatures greater than the main transition temperature of the polymeric matrix (melting temperature for PP) [30].
As shown above, most authors presume isotropic properties or isothermal conditions, which is not the case in industrial processes. This study highlights the importance of including temperature dependent anisotropic thermal conductivity and deconsolidation when simulating the heating step in thermoforming of thermoplastic composites with the help of CFD methods. For this reason, heating experiments in a laboratory thermoforming press are performed with varying heater power, heater-sheet distance, sheet thicknesses, and heater setups. Moreover, a CFD model of the process is created, and three different sets of material properties are used as input data. Starting with a very simple isothermal and isotropic model and increasing the complexity to the point at which the anisotropic thermal conductivity, the anisotropic coefficient of linear thermal expansion and the specific heat capacity at constant pressure are temperature dependent. Including these data in a CFD-model for the heating of thermoplastic composites is a novelty to the field. It is shown that the temperature distribution in the sheets can only be predicted at a satisfactory level when the highest level of complexity is used. Simple, isotropic models are proven to be inaccurate and unusable.

2. Materials and Methods

In order to examine the infrared heating behavior of thermoplastic composite sheets, it was necessary to manufacture sheets with embedded thermocouples. The chosen material were pre-consolidated PP/GF sheets type Tepex dynalite 104 RG600(4)/47% with a thickness of 2 mm from LANXESS Deutschland GmbH, Cologne, Germany. Those sheets are made up of 2/2 twill-woven glass fibers, embedded in a polypropylene matrix system.
Those sheets were stacked in a parallel plate mold (350 mm × 250 mm) with varying numbers of thermocouples in between the layers and/or on the top and bottom surfaces. Since the sheets were later cut to four smaller samples, the thermocouples were placed correspondingly, as shown in Figure 1. The samples (II) in the figure are used to double-check the results from heating experiments with samples (I) and therefore only exhibit one embedded thermocouple. Furthermore, several samples with only 2 stacked sheets (Figure 1a(III)) were manufactured.
The manufacturing of sheets containing thermocouples at a known position was essential for this study and yet, this was the most difficult task since the mold is almost perfectly sealed when closed which led to many broken thermocouple-cables and unusable composites sheets. The sealing is necessary to keep the matrix from being pressed out of the mold and to keep the fiber volume fraction constant.
The pressing of the thick sheets, inspired by the works of Kiss et al. [31,32,33] was performed in a two-stage process using a heating press Wickert WLP 80/4/3 (Wickert Maschinenbau GmbH, Landau in der Pfalz, Germany) and a cooling press Höfer H10 (Höfer Presstechnik GmbH, Taiskirchen, Austria) which are connected by a mold transfer shuttle system.
Temperatures were monitored during every impregnation trial to guarantee equal processing conditions with the thin-film thermocouples type-K, 402–716 (TC Mess-und Regeltechnik GmbH, Mönchengladbach, Germany) and the heating and cooling processes were logged with a Picotech TC-08 logger (Pico Technology, St. Neots, United Kingdom). The following steps were necessary to fabricate the samples:
  • Heating the hot press to 300 °C;
  • Insertion of the mold containing the stack and applying a pressure of 2 bar;
  • Once the thermocouples read a temperature of over 165 °C, the pressure was increased to 10 bar and held until a temperature of 190 °C is reached;
  • The hot press was opened and the mold was transferred to the cooling press;
  • The mold was cooled to a temperature below 100 °C at a pressure of 5 bar;
  • The sheets were demolded using ejector pins.
The sheets were then cut to receive the final parts with 175 mm × 120 mm.
The heating experiments were conducted using a thermoforming press type LZT-OK-220-L (Langzauner GmbH, Lambrechten, Austria). This press was equipped with double-sided infrared heating utilizing a total of 18 IR-heaters type Krelus MINI G14-25 M (Leister Technologies AG, Sarnen, Switzerland) with a maximum power of 2.5 kW each. The heaters were vertically adjustable in a range of 140 mm to 440 mm from each other and the maximum power can be adjusted individually. Moreover, the press utilizes a shuttle system in which the semifinished parts were clamped in order to transfer the sheets from the press to the heating station and back to the press when the forming temperature was reached. This work deals solely with the heating step and forming was not considered. The setup is depicted in Figure 2.
The sheets were first clamped into the shuttle frame and the thermocouples were connected to the logger and a semiautomatic process was started:
(1)
Transfer the sheet to the heating station;
(2)
Activate infrared heating at a predefined power level for a set amount of time;
(3)
Shut off heaters and transfer the sheet to the molding station;
(4)
Activate convectional cooling.
The temperature of the infrared heaters was determined using an infrared camera system type Optris Xi 400 (Optrix GmbH, Holzkirchen, Germany). One single heater was activated at different power settings and the mean surface temperature was recorded until a steady state was reached.

3. CFD Modeling

The Model was created in CATIA V5 2012 (Dassault Systèmes, Vélizy-Villacoublay, France) and mimics the heating station of the thermoforming machine. Figure 2 shows the used model featuring 18 individual heaters of 250 by 250 mm, the clamping frame and the thermoplastic composite sheet.
This model is then exported to ANSYS Workbench 2021 R2 (Ansys Inc., Canonsburg, PA, USA). This tool supports a wide range of customized materials and heating mechanisms. For this work, three different thermal material behaviors were studied. They differ in their linear thermal expansion, specific heat capacity and thermal conductivity. For further reference, the naming and general physical behavior are summarized in Table 1.
The calculation of the anisotropic thermal conductivity was presented in detail in preceding studies of Längauer et al. [29,30]. The basic principle is to define a unit cell of the composite sheet with an edge length a and a roving width b from which the thermal material properties are calculated. The fibers are considered to be located in a “fiber-rich” layer with the thickness s1, that is unharmed by thermal deconsolidation. The “matrix-rich” layer with the thickness s2 is subject to deconsolidation and allows parallelly and serially conducting voids to form above the transition temperature. The parallelly conducting voids are located within the “matrix-rich” layer and therefore increase its thickness with rising temperature. The serially conducting voids form a layer of their own with the thickness sv. The total void volume is calculated from the temperature dependent fiber volume fraction φ ( T ) that is derived numerically. Basically, the total volume of the unit cell V T ( T ) is increasing with rising temperature
V T ( T ) = a ( 1 + α x Δ T ) a ( 1 + α y Δ T ) s T ( 1 + α z Δ T ) ,
where α x , y , z are the coefficients of linear thermal expansion in the directions of space, V T ( T ) is the temperature dependent total volume of the unit cell, s T is the total unit cell thickness and Δ T is the temperature difference. The thermal expansion in both fiber directions was assumed to be constant and limited to that of glass fibers [31]. The coefficient of linear thermal expansion in transversal direction α z for the anisotropic and deconsolidating case is therefore calculated as
α z = 1 Δ T ( V T ( T ) ( a ( 1 + α F Δ T ) ) 2 s T 1 )   ,
With
V T ( T ) = V F ( T ) φ ( T )
Here, V F ( T ) is the temperature dependent fiber volume which is determined as
V F ( T ) = m A a 2 ρ F ( 1 + α F Δ T ) 3 ,
where m A is the area density of the textile and ρ F is the fiber density.
The exact input data for the transversal coefficient of linear thermal expansion is sketched in Figure 3. The isotropic case considers an isothermal coefficient in all directions of space. This value is provided by the material datasheet. The transversal coefficient of linear expansion for the anisotropic material was set to 8.25 × 10−5 K−1, which corresponds to the mean value of the coefficients of a typical PP grade and glass fibers [34,35]. The transversal coefficient of linear expansion of the anisotropic and deconsolidating material according to Equation (2) exhibits a rapid increase at the melting temperature of the polymeric matrix and an almost steady decline from there. This of course has an impact on the density function that is required for the thermal simulation.
Following [29], the thermal conductivity in the fiber direction k for the anisotropic case was calculated as
k = ( A F , 2 k F + A P , 2 k P ) 1 + a b b · ( A F , 2 k F + A P , 2 k P ) ( A F , 1 k F + A P , 1 k P ) 1 ( s 1 + s 2 ) · b ,
where kP and kF are the thermal conductivity of the polymer and the fiber material, respectively and AF and AP are the corresponding areas of fiber and polymer, respectively, in each layer (with thickness s1 and s2), as the indices 1 and 2 imply [29].
In the preceding study [30], it was concluded that the fiber material is the dominant factor for heat conduction and that thermal deconsolidation cannot break the fiber network. For this reason, the thermal conductivity in the fiber direction is constant over the deconsolidation temperature resulting in the same values for k for the anisotropic and the anisotropic and deconsolidating case.
According to [30] for the transversal direction, the thermal conductivity k for the anisotropic case is calculated as
k = s 1 + s 2   s 2 · k P 1 + s 1 s 2 · k P a 2 A P 1 k P + A F k F ,
where AP and AF are again the corresponding areas of polymer and fiber, respectively, and the index 1 stands for the “fiber-rich” layer. When thermal deconsolidation takes place, a new void layer is formed, and the “matrix-rich” layer expands leading to a temperature dependent transversal thermal conductivity k ( T )
k ( T ) = k v   k p a r a l l e l ( s 1 + s 2 + V v   T T 0 a 2 + s v ) k v ( s 1 + s 2 + V v   T T 0 a 2 ) + k p a r a l l e l s v ,
where Vv is the initial void volume, T0 is the initial temperature, kv is the thermal conductivity of the voids and k p a r a l l e l is the thermal conductivity in the original unit cell, considering the corresponding areas of voids, fibers and matrix but neglecting the newly created void layer.
k p a r a l l e l = ( A P 2 k P + A v k v ) ( A F k F + A P 1 k P ) ( s 1 + s 2 + V v   T T 0 a 2 ) ( A P 2 k P + A v k v ) a 2 s 1 + a 2 ( A F k F + A P 1 k P ) ( s 2 + V v   T T 0 a 2 ) .
The thickness of the void layer s v is derived from the subtraction of the void volume increase in layer 2 given by the general gas equation from the total expansion calculated from the temperature dependent fiber volume fraction and the initial fiber volume fraction φ 0
s v = m A ρ f ( 1 φ ( T ) 1 φ 0 ) V v   T T 0 a 2
Figure 4 shows the input data for the thermal conductivity. As one can see, PP/GF isotropic uses constant values for all directions. The anisotropic case is a bit more advanced and uses values from Equations (5) and (6) for the anisotropic thermal conductivity, distinguishing between the fiber and the transversal direction. The input data which is also used for anisotropic and deconsolidating is provided in Table 2. Equation (7) is utilized to receive the transversal thermal conductivity for the anisotropic and deconsolidating case in which there is a steady decline until reaching the melting point of the matrix material where a drop in the thermal conductivity is visible due to the rapid void growth and creation. From there the transversal thermal conductivity is declining in a linear fashion again.
The temperature dependent specific heat capacity was previously tested and is derived from [29]. The isothermal value for the specific heat capacity at constant pressure of the isotropic material was set to 1.01 J·g−1·K−1, since this is the rounded measuring value at 30 °C.
The resulting values were provided to the engineering data for each material. In a following step, the transient heating evaluation was performed using surface to surface radiation from the heaters to the sheet (emissivity of one). Furthermore, convection in the heating station was considered at all sheet surfaces. The heat transfer coefficient was calculated according to [9] and is between 3.6 W·m−2·K−1 and 10.6 W·m−2·K−1 in the relevant temperature range of 30 °C to 300 °C. The meshing of the heaters and the aluminum clamping frame was uncritical, so a maximum element size of 5 cm and auto fitting was used (Figure 5a). In order to minimize the skewness in critical areas, face sizing with an element size of 2 mm was used. This is clearly visible at the clamps of the frame and the reflectors of the heaters (red arrows in Figure 5a). The sheet itself was meshed with body sizing, cartesian elements and quadratic element order with a maximum edge length of 2 mm to provide sufficient data (Figure 5b). Moreover, a comparison with an element size of 1 mm was made (Figure 5c).
Table 3 contains the information about the different components of the model and their mesh statistics. The skewness is supposed to be smaller than 0.9, the element quality as well as the aspect ratio should be close to 1. At the heaters and the clamping frame, which are both not as relevant for the simulation as the sheet, the average values for those parameters are in an acceptable range. The sheet was meshed twice with differing element count. Yet, Table 3 shows that the mesh statistics are almost the same for both cases with only minor deviations.
The heating was simulated for at least 90 s, and as a result the total temperature distribution of the part and the temperature of the measuring spots on and in the sheet was received for each condition. This was accomplished by placing temperature probes in the model at local coordinate systems.
The temperature of the heating elements was calculated according to Equation (4) for each heater setting. Figure 6 shows the results of the temperature readings of the heating tests of the single heater in the laboratory thermoforming press at different conditions over time. The infrared camera records one reading for each pixel. Yet, in the simulation a mean value for the heater temperature is used for the infrared radiation heating which is mimicked by averaging over a large area on the heater surface—even the unheated crosspiece in the center.
A logarithmic fit with the parameters C and D was applied to calculate the heater temperature TH over the time t
T H = C ln ( t ) + D
Table 4 shows the used fit parameters for various heater settings.

4. Results

The heating tests at the laboratory thermoforming station revealed a very inhomogeneous temperature distribution over the sheet thickness but also at the sheet surfaces. This can also be seen in Figure 7, where a comparison of the simulated sheet surface between isotropic (a) and anisotropic and deconsolidating (b) material behavior is made. The isotropic case reacts in a more homogenous way to the temperature increase but appears overheated in almost all areas (T > 300 °C). Only in the areas in and surrounding the clamps the temperature is lower. The anisotropic and deconsolidating case shows signs of overheating only in the outer areas, where the heat cannot be transported into the sheet in a sufficient manner. Figure 7c shows a picture of a heated sheet where the signs of overheating (unwetted, brittle, white fiber rovings) at the edges are very prominent (dashed region), and, also, the still consolidated regions at the clamps can clearly be located (dotted region). The center of the sheet exhibits typical signs of deconsolidation.
In order to validate the meshing result, the core temperature curves for the utilized fine and coarse mesh at the same heating conditions (4 mm sheet thickness, 80% heater power, and 150 mm sheet to heater distance) are plotted in Figure 8. The plots are congruent for the whole monitored heating time with one minor deviation at the melting point, which is highlighted in the figure as a zoomed plot. Since this discrepancy is insignificant, the simulations are commenced with the coarse mesh.
Figure 9 shows the measured core temperature of composite samples with 4 mm thickness and a sheet to heater distance of 150 mm on each side dependent on the heating time for different heater power settings and the corresponding model values derived from CFD-analysis. As expected, the sheet’s inner plies are heated faster with higher power. The isotropic material shows an almost linear heating behavior, that is clearly contradicted by the real course of the curve. The anisotropic case and the anisotropic and deconsolidating case are fairly close to each other, yet the major difference is the behavior after the transition temperature at 165 °C when thermal deconsolidation takes place. From there, the slope of the curve of the experiment is much closer to the model one, which considers thermal deconsolidation. The deviations at the end of the experiment are 7.6 °C (3.4% relative error), 15.9 °C (7% relative error), and 90.7 °C (40% relative error) at 80% power for the anisotropic and deconsolidating, anisotropic, and isotropic case, respectively. At 100% heater power, the surface of the composite is overheated quickly, which is why the experimental temperature in the core does not even reach temperatures above 200 °C before the test has to be aborted.
Figure 10 shows the relative errors of the heating tests dependent on the heater distance setting for samples of 4 mm thickness, comparing them directly to the model values for all viewed material properties. Again, the most reliable results are derived from the anisotropic and deconsolidating case with lower deviations than the other material property cases for both heater distance settings. It is noticeable that at both heater distances, the error reaches a local maximum after a very short heating period. With time, the error drops to values of 0% for the anisotropic and the anisotropic and deconsolidating case when approaching the melting temperature. From there, the errors show a small incline, reaching up to 2.9% at 150 mm and 8.6% at 200 mm heater distance for the anisotropic and deconsolidating case, and 7.4% at 150 mm and 13.7% at 200 mm heater distance for the anisotropic case. The error in the isotropic case is close to 50% for both heater distance settings.
Figure 11 shows the temperature distribution at heating times of 20 s and 140 s from heating tests of the 8 mm material when only one single heater is activated on the top and the bottom heater bank with a power setting of 80% and a distance of 150 mm with the corresponding model values and the experimental data. This sheet was equipped with five thermocouples, that all passed the heating test unharmed. For this case, the whole cross-sectional temperature distribution can be investigated.
After 20 s of heating time, the temperature in the core of the sheet has not changed dramatically, and all of the models used deliver results within the accepted margin. With longer heating times, the isotropic case once again overshoots, especially in the center layers of the sheet. The anisotropic and the anisotropic and deconsolidating cases perform very well. The difference is not as clear, since the melting temperature has not been reached in the core after 140 s heating time. On the outside, the values are almost matching until the melting point is reached. In the experiments, this is also where thermal deconsolidation takes place first. The thermocouples tend to lose the perfect contact to the sheets, and the experimental values are not completely trustworthy any longer.
It is, also, notable that the measured temperature in the experiments is not symmetrical over the sheet thickness. The top surface shows higher readings than the bottom surface, even though the sheet to heater distance is the same. The reason for this behavior could be ascending heated air that is trapped under the upper heater bank which provides additional thermal input.

5. Discussion

The work shows the numerical and experimental heating analysis of thermoplastic composite sheets in a typical heating step of a thermoforming application. It can clearly be demonstrated by the experiments, that the temperature distribution through thickness and on the surface of the composite sheets is not uniform but rather strongly inhomogeneous. The only model that delivers satisfactory results, considers the anisotropic thermal conductivity, the anisotropic coefficient of linear thermal expansion, and the specific heat capacity at constant pressure to be temperature dependent.
Without those considerations, the models prove to be very inaccurate. Especially isotropic modeling leads to massive errors (40–50%) compared to the experiments. The consequence when using such models for process simulations would be significant under-heating of the parts. The sheets would be formed with a core temperature of up to 25 °C below the melting point, which will lead to inferior mechanical and optical properties.
When anisotropy is considered but deconsolidation is neglected, the simulation matches the experiments in a satisfactory way when the temperature of the sheet is below the transition temperature of the polymeric matrix. The knowledge of the temperature dependent specific heat capacity is a great asset as well, as it leads to a kink in the heating curve, that is also found in real experiments. From the point where the transition temperature is reached, a good estimate of the core temperature of the composite can only be made when using material models that consider thermal deconsolidation.
The combination of thermal deconsolidation and anisotropic thermal conductivity in CFD methods is new to the field and shows excellent agreement with actual experimental results.
In almost every numerical thermoforming study cited in the introduction, a uniform temperature distribution in the sheet is presumed before the forming step [19,20,22,23,24]. The knowledge of the real temperature distribution after heating could lead to significantly different results. Due to the inaccuracy of the heating modeling, not only the temperature distribution of the sheets is affected, in the long run this also influences cycle time calculations and over- and underheating issues and therefore the mechanical properties of the parts. In [28], an isotropic approach was used, which also matched the experimental data. This was enabled by numerical fitting of a parameter to calculate the heat transfer coefficient to equate the temperature predictions of the model and the results of the experiment.
Judging from practical experience, PP/GF sheets are well formable at core temperatures of 200 °C. In industrial processes, core temperature tests are not available which is why, mostly, a trial-and-error principle is used. Figure 12 shows experimental sheet core temperatures at the time each model predicts 200 °C core temperature for 200 mm sheet to heater distance and 100% heater power and for 150 mm sheet to heater distance and 80% heater power. Furthermore, the figure shows the heating time that is considered for each case. The experimental heating time for those cases are 70 s and 87 s for a distance of 200 mm and 100% power and for 150 mm distance and 80% power, respectively. The red dashed line marks the transition temperature of the PP matrix. If these models were used in industrial processes, the sheets would not be formable for the isotropic case. The anisotropic material is a borderline case at the first setting and would be formable at the second one. Only the anisotropic and deconsolidating case is within the forming window of the PP matrix for both settings. The deviation in the first case looks large, yet the heating time is just 2 s shorter than the true experimental value to reach a temperature of 200 °C in the core.
This work yields the possibility to determine the temperature and volume of a composite sheet in every processing step. The future scope is to use the full potential of the model presented in this study to simulate the entire thermoforming process from the heating step to transportation and forming of the sheets to cooling and demolding of the final part.

6. Conclusions

This work shows how the use of novel, state-of-the-art material models for the temperature dependent anisotropic thermal conductivity and the anisotropic coefficient of linear thermal expansion can influence the results of CFD-modelling of thermoforming applications for thermoplastic composite materials in comparison to much simpler models. Since these materials must fulfil the highest standards in terms of mechanical behavior and aesthetics, precise process modeling is demanded. The numerical simulation is validated with experiments in which PP/GF sheets are heated in a laboratory thermoforming machine using different sheet to heater distances, heater power settings, and sheet thickness.
It can be shown that the anisotropic effects of thermal deconsolidation strongly influence the thermal conductivity and the linear coefficient of thermal expansion. This behavior can solely be simulated using the presented temperature dependent models, as the resulting comparisons between simulation and experimental data prove. The error in the heating time prediction can be lowered to under 3% for thick composite sheets. Using simpler models leads to unacceptable errors and inefficient processes. The heating time estimation herein delivers an error of over 25%, compared to experimental data. This work thus presents a powerful tool for thermoforming simulations and will help to improve processes and products.

Author Contributions

Conceptualization, M.L., G.Z. and C.H.; methodology, M.L. and G.Z.; formal analysis, M.L.; investigation, M.L. and H.S.; writing—original draft preparation, M.L.; writing—review and editing, M.L., G.Z., C.H. and H.S.; visualization, M.L.; supervision, C.H. and G.Z.; project administration, G.Z.; funding acquisition, G.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the federal government of Upper Austria and the European Union, grant number Wi-2015-132734 for the project “ProFVK”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Li, N.; Yang, X.; Bao, F.; Pan, Y.; Wang, C.; Chen, B.; Zong, L.; Liu, C.; Wang, J.; Jian, X. Improved Mechanical Properties of Copoly(Phthalazinone Ether Sulphone)s Composites Reinforced by Multiscale Carbon Fibre/Graphene Oxide Reinforcements: A Step Closer to Industrial Production. Polymers 2019, 11, 237. [Google Scholar] [CrossRef] [PubMed]
  2. Wang, B.; Li, N.; Cheng, S.; Hu, F.; Li, G.; Guo, H.; Liu, C.; Chen, Y.; Jian, X. Thermal conductivity and mechanical properties enhancement of CF/PPBESK thermoplastic composites by introducing graphene. Polym. Compos. 2022, 43, 2736–2745. [Google Scholar] [CrossRef]
  3. Wang, B.; Li, N.; Cheng, S.; Hu, F.; Bao, Q.; Hao, H.; Liu, C.; Chen, Y.; Jian, X. Enhance the thermal conductivity and mechanical properties of CF/PPBESK thermoplastic composites by growth ZnO nanowires with tunable length and diameter on prepreg. Polymer 2022, 244, 124662. [Google Scholar] [CrossRef]
  4. Vaidya, U.K.; Chawla, K.K. Processing of fibre reinforced thermoplastic composites. Int. Mater. Rev. 2008, 53, 185–218. [Google Scholar] [CrossRef]
  5. Pühringer, J.F.; Zitzenbacher, G.; Spreitzer, C. Study of Heat Absorption in the Thermoforming Process for Transparent and Filled Polystyrene. In Proceedings of the 27th Annual Meeting of the Polymer Processing Society, Marrakesh, Morocco, 10–14 May 2011. [Google Scholar]
  6. Nowacki, J.; Neitzel, M. Thermoforming of Reinforced Thermoplastic Stiffened Structure. Polym. Compos. 2000, 21, 531–538. [Google Scholar] [CrossRef]
  7. Manson, J.A.E.; Wakeman, M.D.; Bernet, N. Composite processing and manufacturing-an overview. Compr. Compos. Mater. 2000, 2, 577–607. [Google Scholar]
  8. Pühringer, J.F.; Zitzenbacher, G.; Spreitzer, C. Study of heat absorption in thermoforming for transparent and filled polystyrene. Int. Polym. Processing 2013, 28, 14–23. [Google Scholar] [CrossRef]
  9. Längauer, M.; Zitzenbacher, G.; Burgstaller, C.; Hochenauer, C. Enhanced Infrared Heating of Thermoplastic Composite Sheets for Thermoforming Processes. Int. Polym. Process. 2021, 36, 35–43. [Google Scholar] [CrossRef]
  10. Wolfrath, J.; Michaud, V.; Manson, J.A.E. Deconsolidation in glass mat thermoplastic composites: Analysis of the mechanisms. Compos. Part A 2005, 36, 1608–1616. [Google Scholar] [CrossRef]
  11. Gröschel, C.; Drummer, D. The Influence of Moisture and Laminate Setup on the De-Consolidation Behavior of PA6/GF Thermoplastic Matrix Composites. Int. Polym. Process. 2014, 29, 660–668. [Google Scholar] [CrossRef]
  12. Ye, L.; Lu, M.; Mai, Y.-W. Thermal de-consolidation of thermoplastic matrix composites—I. Growth of voids. Compos. Sci. Technol. 2002, 62, 2121–2130. [Google Scholar] [CrossRef]
  13. Lu, M.; Ye, L.; Mai, Y.-W. Thermal de-consolidation of thermoplastic matrix composites—II.“Migration” of voids and “re-consolidation”. Compos. Sci. Technol. 2004, 64, 191–202. [Google Scholar] [CrossRef]
  14. Ye, L.; Chen, Z.-R.; Lu, M.; Hou, M. De-consolidation and re-consolidation in CF/PPS thermoplastic matrix composites. Compos. Part A 2005, 36, 915–922. [Google Scholar] [CrossRef]
  15. Shi, H.; Fernandez-Villegas, I.; Bersee, H.E.N. Modelling of heat transfer and consolidation for thermoplastic composites resistance welding. In Proceedings of the 18th International Conference on Composites Materials, Jeju, Korea, 21–26 August 2011; pp. 1–6. [Google Scholar]
  16. Shi, H.; Fernandez-Villegas, I.; Bersee, H.E.N. Analysis of void formation in thermoplastic composites during resistance welding. J. Thermoplast. Compos. Mater. 2016, 30, 1654–1674. [Google Scholar] [CrossRef]
  17. Brogan, M.T.; Monaghan, P.F. Thermal simulation of quartz tube infra-red heaters used in the processing of thermoplastic composites. Compos. Part A 1996, 27, 301–306. [Google Scholar] [CrossRef]
  18. Johnson, A.F.; Picket, A.K. Numerical simulation of the forming process in long fibre reinforced thermoplastics. Trans. Eng. Sci. 1996, 10, 233–242. [Google Scholar]
  19. Sweeney, G.J.; Monaghan, P.F.; Brogan, M.T.; Cassidy, S.F. Reduction of infra-red heating cycle time in processing of thermoplastic composites using computer modelling. Compos. Manuf. 1995, 6, 255–262. [Google Scholar] [CrossRef]
  20. Hwang, S.-F.; Hwang, K.-J. Stamp forming of locally heated thermoplastic composites. Compos. Part A 2002, 33, 669–676. [Google Scholar] [CrossRef]
  21. Hsiao, S.-W.; Kikuchi, N. Numerical Analyses of Deep Drawing Process for Thermoplastic Composites Laminates. J. Eng. Mater. Technol. 1997, 119, 314–318. [Google Scholar] [CrossRef]
  22. Abbassi, F.; Elfaleh, I.; Mistou, S.; Zghal, A.; Fazzini, M.; Djilali, T. Experimental and numerical investigastions of a thermoplastic composites thermoforming. Struct. Control Health Monit. 2011, 18, 769–780. [Google Scholar] [CrossRef]
  23. Chen, Q.; Boisse, P.; Park, C.H.; Saouab, A.; Breard, J. Intra/inter-ply shear behaviors of continuous fiber reinforced thermoplastic composites in thermoforming processes. Compos. Struct. 2011, 93, 1692–1703. [Google Scholar] [CrossRef]
  24. Harrison, P.; Gomes, R.; Curado-Correia, N. Press forming a 0/90 cross-ply advanced thermoplastic composite using the double-dome benchmark geometry. Compos. Part A 2013, 54, 56–69. [Google Scholar] [CrossRef]
  25. Stamopoulos, A.; Di Ilio, A. Numerical and experimental analysis of the thermoforming process parameters of semi-spherical glass fibre thermoplastic parts. Procedia CIRP 2021, 99, 420–425. [Google Scholar] [CrossRef]
  26. D’Emilia, G.; Gaspari, A.; Natale, E.; Stamopoulos, A.; Di Ilio, A. Experimental and numerical analysis of the defects induced by the thermoforming process on woven textile thermoplastic composites. Eng. Fail. Anal. 2022, 135, 106093. [Google Scholar] [CrossRef]
  27. McCool, R.; Murphy, A.; Wilson, R.; Jiang, Z.; Price, M. Thermoforming of Continuous Fibre Reinforced Thermoplastic Composites. AIP Conf. Proc. 2011, 1353, 901–906. [Google Scholar]
  28. Nardi, D.; Sinke, J. Design analysis for thermoforming of thermoplastic composites: Prediction and machine learning-based optimization. Compos. Part C 2021, 5, 100126. [Google Scholar] [CrossRef]
  29. Längauer, M.; Brunnthaller, F.; Zitzenbacher, G.; Burgstaller, C.; Hochenauer, C. Modeling of the anisotropic thermal conductivity of fabrics embedded in a thermoplastic matrix system. Polym. Compos. 2021, 42, 2050–2060. [Google Scholar] [CrossRef]
  30. Längauer, M.; Zitzenbacher, G.; Heupl, S.; Plank, B.; Burgstaller, C.; Hochenauer, C. Influence of thermal deconsolidation on the anisotropic thermal conductivity of glass fiber reinforced, pre-consolidated polypropylene sheets used for thermoforming applications. Polym. Compos. 2022, 43, 2264–2275. [Google Scholar] [CrossRef]
  31. Kiss, P.; Stadlbauer, W.; Burgstaller, C.; Stadler, H.; Fehringer, S.; Haeuserer, F.; Archodoulaki, V.M. Development of high-performance glass fibre-polypropylene composite laminates: Effect of fibre sizing type and coupling agent concentration on mechanical properties. Compos. Part A Appl. Sci. Manuf. 2020, 138, 106110. [Google Scholar] [CrossRef]
  32. Kiss, P.; Schoefer, J.; Stadlbauer, W.; Burgstaller, C.; Archodoulaki, V.M. An experimental study of glass fibre roving sizings and yarn finishes in high-performance GF-PA6 and GF-PPS composite laminates. Compos. Part B Eng. 2021, 204, 108487. [Google Scholar] [CrossRef]
  33. Kiss, P.; Glinz, J.; Stadlbauer, W.; Burgstaller, C.; Archodoulaki, V.M. The effect of thermally desized carbon fibre reinforcement on the flexural and impact properties of PA6, PPS and PEEK composite laminates: A comparative study. Compos. Part B Eng. 2021, 215, 108844. [Google Scholar] [CrossRef]
  34. Domininghaus, H.; Elsner, P.; Eyerer, P.; Hirth, T. Kunststoffe—Eigenschaften und Anwendungen, 7th ed.; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 2008. [Google Scholar]
  35. Terzvergil, A.; Lassila, L.V.J.; Vallittu, P.K. The effect of fiber orientation on the thermal expansion coefficients of fiber-reinforced composites. Dent. Mater. 2003, 19, 471–477. [Google Scholar] [CrossRef]
Figure 1. (a) Manufactured sheet (left) and different sample cross-sections with embedded thermocouples (right). The areas are cross-hatched to enable telling apart different layers of stacked sheets. (I) a sample with 5 thermocouples, (II) a sample with one central thermocouple, and (III) a sample made of two sheets with one central thermocouple. (b) depicts a sheet of type (II).
Figure 1. (a) Manufactured sheet (left) and different sample cross-sections with embedded thermocouples (right). The areas are cross-hatched to enable telling apart different layers of stacked sheets. (I) a sample with 5 thermocouples, (II) a sample with one central thermocouple, and (III) a sample made of two sheets with one central thermocouple. (b) depicts a sheet of type (II).
Polymers 14 03331 g001
Figure 2. Model of the heating station of the thermoforming device (top) and a photograph of the real heating station (bottom), where (1) is the upper heater bank, (2) is the thermoplastic composite sample, (3) is the clamping frame, and (4) is the lower heater bank.
Figure 2. Model of the heating station of the thermoforming device (top) and a photograph of the real heating station (bottom), where (1) is the upper heater bank, (2) is the thermoplastic composite sample, (3) is the clamping frame, and (4) is the lower heater bank.
Polymers 14 03331 g002
Figure 3. Transversal coefficient of linear thermal expansion α z for anisotropic and deconsolidating composites and isotropic materials with isothermal behavior.
Figure 3. Transversal coefficient of linear thermal expansion α z for anisotropic and deconsolidating composites and isotropic materials with isothermal behavior.
Polymers 14 03331 g003
Figure 4. Temperature dependent thermal conductivity used in the CFD models for each specific material behavior.
Figure 4. Temperature dependent thermal conductivity used in the CFD models for each specific material behavior.
Polymers 14 03331 g004
Figure 5. (a) Meshed model in ANSYS Workbench 2021 R2 (the red arrows indicate areas with finer face meshing) with a detailed view of the meshed, clamped thermoplastic composite sheet in (b) and a finer meshed model for the sheet with a suppressed mesh at the clamping frame (c).
Figure 5. (a) Meshed model in ANSYS Workbench 2021 R2 (the red arrows indicate areas with finer face meshing) with a detailed view of the meshed, clamped thermoplastic composite sheet in (b) and a finer meshed model for the sheet with a suppressed mesh at the clamping frame (c).
Polymers 14 03331 g005
Figure 6. (a) Infrared camera image of a single heater. The white frame shows the region of interest, that was selected to determine the mean heater temperature. (b) Fitted heater temperature over the heating time with varying heater settings including one experimental series over time and the maximum measured mean heater temperature for each condition.
Figure 6. (a) Infrared camera image of a single heater. The white frame shows the region of interest, that was selected to determine the mean heater temperature. (b) Fitted heater temperature over the heating time with varying heater settings including one experimental series over time and the maximum measured mean heater temperature for each condition.
Polymers 14 03331 g006
Figure 7. Surface temperature distribution after 100 s of heating with a power of 80% at a distance of 150 mm in case of an isotropic sheet of 4 mm thickness (a) and an anisotropic and deconsolidating sheet of 4 mm thickness (b) in comparison to a real composite sheet that has been subjected to these heating conditions (c). The dashed area indicates signs of overheating whereas the dotted area is still consolidated.
Figure 7. Surface temperature distribution after 100 s of heating with a power of 80% at a distance of 150 mm in case of an isotropic sheet of 4 mm thickness (a) and an anisotropic and deconsolidating sheet of 4 mm thickness (b) in comparison to a real composite sheet that has been subjected to these heating conditions (c). The dashed area indicates signs of overheating whereas the dotted area is still consolidated.
Polymers 14 03331 g007
Figure 8. Comparison of two simulations of the heating of a sheet with a thickness of 4 mm at 150 mm sheet to heater distance and 80% heater power for two different meshes with a detailed view (dashed area) of the region of the melting point of the polypropylene matrix.
Figure 8. Comparison of two simulations of the heating of a sheet with a thickness of 4 mm at 150 mm sheet to heater distance and 80% heater power for two different meshes with a detailed view (dashed area) of the region of the melting point of the polypropylene matrix.
Polymers 14 03331 g008
Figure 9. Result of heating tests of a 4 mm PP-GF sheet including the values derived from CFD-modeling. (a) Shows the values for 150 mm sheet to heater distance and 80% heater power and (b) shows the values for 150 mm sheet to heater distance and 100% heater power.
Figure 9. Result of heating tests of a 4 mm PP-GF sheet including the values derived from CFD-modeling. (a) Shows the values for 150 mm sheet to heater distance and 80% heater power and (b) shows the values for 150 mm sheet to heater distance and 100% heater power.
Polymers 14 03331 g009
Figure 10. Relative errors of the CFD-model cases compared to the experimental values when heating a 4 mm PP/GF sheet at 100% heater power and 150 mm sheet to heater distance (a) and 100 % heater power and 200 mm sheet to heater distance (b).
Figure 10. Relative errors of the CFD-model cases compared to the experimental values when heating a 4 mm PP/GF sheet at 100% heater power and 150 mm sheet to heater distance (a) and 100 % heater power and 200 mm sheet to heater distance (b).
Polymers 14 03331 g010
Figure 11. Temperature distribution through the thickness of a PP/GF sheet with a thickness of 8 mm that is heated with one active top and one active bottom heater running at 80% power at a distance of 100 mm at two heating times.
Figure 11. Temperature distribution through the thickness of a PP/GF sheet with a thickness of 8 mm that is heated with one active top and one active bottom heater running at 80% power at a distance of 100 mm at two heating times.
Polymers 14 03331 g011
Figure 12. Experimental sheet core temperature at the time each model predicts 200 °C core temperature for two different settings (left: heater power of 100% at a distance of 200 mm; right: heater power of 80% at a distance of 150 mm) and a sheet of 4 mm thickness. The numbers above the bars are the experimental heating time at which the simulated data was retrieved. The red dashed line indicates the melting point of the PP matrix.
Figure 12. Experimental sheet core temperature at the time each model predicts 200 °C core temperature for two different settings (left: heater power of 100% at a distance of 200 mm; right: heater power of 80% at a distance of 150 mm) and a sheet of 4 mm thickness. The numbers above the bars are the experimental heating time at which the simulated data was retrieved. The red dashed line indicates the melting point of the PP matrix.
Polymers 14 03331 g012
Table 1. Studied material behaviors and the corresponding material properties used for modeling.
Table 1. Studied material behaviors and the corresponding material properties used for modeling.
NameLinear Thermal ExpansionSpecific Heat CapacityThermal Conductivity
isotropicisothermal and isotropicisothermalisotropic and isothermal
anisotropicisothermal and anisotropictemperature dependentanisotropic
anisotropic and deconsolidatinganisotropic and temperature dependenttemperature dependentanisotropic and subject to thermal deconsolidation
Table 2. Material data used for the modeling of the anisotropic thermal conductivity and the deconsolidation volume [29,30].
Table 2. Material data used for the modeling of the anisotropic thermal conductivity and the deconsolidation volume [29,30].
ParameterDescriptionValue
φ Initial fiber volume fraction0.47
kvThermal conductivity (voids)0.026 W m−1 K−1
kFThermal conductivity (fibers)1.03 W m−1 K−1
kPThermal conductivity (polymer)0.20 W m−1 K−1
aUnit cell edge length4.1 mm
bRoving width4.0 mm
ρFFiber density2.5 g cm−3
αFCoefficient of linear thermal expansion of the fibers5 × 10−6 K−1
αPCoefficient of linear thermal expansion of the polymer1.6 × 10−4 K−1
mAArea density of the textile600 g m−2
vm0Initial matrix volume fraction0.53
mPower law index, dependent on fiber configuration10.2
pPressure at deconsolidation101,000 Pa
p0Initial pressure202,000 Pa
T0Initial temperature298.15 K
φ 0 fiber volume fraction in the relaxed state0.3
CcEfMaterial constant multiplied by the fiber bending modulus350 Pa
Table 3. Mesh statistics of the CFD model for the heaters and the clamping frame and two different element counts at the sheet (Avg. = Average value; σ = standard deviation).
Table 3. Mesh statistics of the CFD model for the heaters and the clamping frame and two different element counts at the sheet (Avg. = Average value; σ = standard deviation).
MeshElement
Count
SkewnessElement
Quality
Aspect Ratio
minmaxAvg.σminmaxAvg.σminmaxAvg.σ
Heaters and
Clamping
Frame
fine2.4 × 1067 × 10−40.840.260.110.20.990.810.091.28.81.910.42
Sheetfine84,00010−1010−1010−10011110−71110
coarse10,56010−1010−1010−1000.990.990.999 × 10−81.0061.0061.00610−7
Table 4. Fit parameters to calculate the heater temperature for various heater settings.
Table 4. Fit parameters to calculate the heater temperature for various heater settings.
Heater Setting (%)Fit Parameter CFit Parameter D
2028.0783.51
3035.04111.41
4041.73138.19
5047.50161.30
6051.59177.66
7057.13199.85
8060.63213.83
9064.57229.63
10068.17244.05
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Längauer, M.; Zitzenbacher, G.; Stadler, H.; Hochenauer, C. Enhanced Simulation of Infrared Heating of Thermoplastic Composites Prior to Forming under Consideration of Anisotropic Thermal Conductivity and Deconsolidation by Means of Novel Physical Material Models. Polymers 2022, 14, 3331. https://doi.org/10.3390/polym14163331

AMA Style

Längauer M, Zitzenbacher G, Stadler H, Hochenauer C. Enhanced Simulation of Infrared Heating of Thermoplastic Composites Prior to Forming under Consideration of Anisotropic Thermal Conductivity and Deconsolidation by Means of Novel Physical Material Models. Polymers. 2022; 14(16):3331. https://doi.org/10.3390/polym14163331

Chicago/Turabian Style

Längauer, Manuel, Gernot Zitzenbacher, Hannes Stadler, and Christoph Hochenauer. 2022. "Enhanced Simulation of Infrared Heating of Thermoplastic Composites Prior to Forming under Consideration of Anisotropic Thermal Conductivity and Deconsolidation by Means of Novel Physical Material Models" Polymers 14, no. 16: 3331. https://doi.org/10.3390/polym14163331

APA Style

Längauer, M., Zitzenbacher, G., Stadler, H., & Hochenauer, C. (2022). Enhanced Simulation of Infrared Heating of Thermoplastic Composites Prior to Forming under Consideration of Anisotropic Thermal Conductivity and Deconsolidation by Means of Novel Physical Material Models. Polymers, 14(16), 3331. https://doi.org/10.3390/polym14163331

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