Unterschiede
Hier werden die Unterschiede zwischen zwei Versionen angezeigt.
| Nächste Überarbeitung | Vorhergehende Überarbeitung | ||
| en:grundlagenhandbuch:drehmoment_und_antriebsleistung:enthalpiemodell [2025/05/15 08:22] – angelegt deppe2 | en:grundlagenhandbuch:drehmoment_und_antriebsleistung:enthalpiemodell [2026/02/07 16:04] (aktuell) – [Validation] deppe2 | ||
|---|---|---|---|
| Zeile 1: | Zeile 1: | ||
| ====== Enthalpy Model ====== | ====== Enthalpy Model ====== | ||
| - | *[[en: | + | ===== Power calculation via enthalpy ===== |
| - | | + | |
| - | *[[en:Grundlagenhandbuch:Drehmoment und Antriebsleistung:Enthalpiemodell: | + | In order to guarantee the three goals of improving accuracy, reducing complexity and considering all process zones in the new power model, a model approach was developed using the first law of thermodynamics. To understand the modelling of the model, the enthalpy input into the polymer is shown formulaically in the three process zones, solids conveying zone, melting zone and melt zone. |
| - | *[[en: | + | |
| - | | + | ===== Solids conveying zone ===== |
| - | *[[en:Grundlagenhandbuch: | + | |
| - | *[[en:Grundlagenhandbuch: | + | First of all, a specific enthalpy curve of a semi-crystalline plastic is considered (Picture 1). At the crystallite melting temperature, |
| - | *[[en:Grundlagenhandbuch: | + | |
| - | *[[en:Grundlagenhandbuch:Drehmoment und Antriebsleistung:Enthalpiemodell: | + | {{ :en:grundlagenhandbuch:drehmoment_und_antriebsleistung:en_sigma150_dlg_grundlagenhandbuch_drehmoment_und_antriebsleistung_009.png? |
| - | *[[en:Grundlagenhandbuch: | + | |
| + | **Picture 1:** Specific enthalpy curve of a partially crystalline polymer | ||
| + | |||
| + | The specific enthalpy change at a supporting point compared to the filling point of the material results from the multiplication of the solid enthalpy by a percentage factor which is dependent on the filling temperature, | ||
| + | |||
| + | $$\Delta h_{FF} = \Delta h_f \cdot \frac{T_{FS} - T_E}{T_K - 0}\tag{Equation 0-1}$$ | ||
| + | |||
| + | For the driving power of an element area the difference of the specific enthalpy change between two supporting points is calculated with equation 0-2. | ||
| + | |||
| + | $$\Delta \Delta h_{FF} = \Delta h_{FF_n} - \Delta h_{FF_{n-1}}\tag{Equation 0-2}$$ | ||
| + | |||
| + | Using the calculated specific enthalpy difference, the pressure difference and the introduced heat flow, the required drive power required between two support points in the solids transport zone can be calculated from this (equation 0-3). | ||
| + | |||
| + | $$P_{D/FF} = \dot{m} \ast (\Delta \Delta h_{FF} + \frac{\Delta p}{\rho}) - \dot{Q}\tag{Equation 0-1}$$ | ||
| + | |||
| + | ===== Melting zone and melt conveying zone ===== | ||
| + | |||
| + | The melting zone and the melt conveying zone differ from the solids conveying zone in that the material is two-phase. Consequently, | ||
| + | |||
| + | For the specific enthalpy change of a pure heating of the melt, the following applies in relation to the crystallite melting temperature through the integration of the specific heat capacity via the temperature equation 0-4. | ||
| + | |||
| + | $$\Delta h_{RS} = cp_0 \ast (T_M - T_K) + \frac{m_{cp}}{2} \ast (T_M^2 - T_K^2)\tag{Equation 0-4}$$ | ||
| + | |||
| + | In addition to the enthalpy, energy was already introduced into the molten material by heating the melt in order to heat the solid up to the crystallite melting temperature. Formally, the specific enthalpy introduced up to the melting point can be determined according to equation 0-1 with equation 0-5. | ||
| + | |||
| + | $$\Delta h_{FA} = \Delta h_a + \Delta h_f \ast \frac{T_K - T_E}{T_K - 0}\tag{Equation 0-5}$$ | ||
| + | |||
| + | From the addition of the equations 0-4 and 0-5, equation 0-6 results for the specific enthalpy increase from the starting point of the filling temperature to the current melt temperature for the melt range: | ||
| + | |||
| + | $$\Delta h_S = \Delta h_a + \Delta h_f \ast \frac{T_K - T_E}{T_K - 0} + cp_0 \ast (T_M - T_K) + \frac{cp_m}{2} \ast (T_M^2 - T_K^2)\tag{Equation 0-6}$$ | ||
| + | |||
| + | The solid bed in the melting zone, on the other hand, undergoes a different specific enthalpy increase. This can be determined with equation 0-7, following equation 0-1 in the solids conveying area. | ||
| + | |||
| + | $$\Delta h_{FA} = \Delta h_f \ast \frac{T_{FS} - T_E}{T_K - 0}\tag{Equation 0-7}$$ | ||
| + | |||
| + | The enthalpy increases of the solid bed and the melt zone must now be weighted and added up on the basis of the degree of melting in order to obtain the total enthalpy change at a supporting point in the melting zone and the melt conveying zone. For weighting, the total mass flow is multiplied by the respective percentage of the solid or melt. By adding the two specific enthalpy changes, the total enthalpy increase at one support point can be calculated. The differentiation between melting zone and melt area is made by the degree of melting. This is calculated in advance in SIGMA and can also be displayed in visual form. | ||
| + | |||
| + | In the pure melt conveying zone, for example, the melting degree is equal to one, so that the second term of equation 0-8 is omitted. | ||
| + | |||
| + | $$\Delta h_{AS/S} = m_{asv} \ast \Delta h_S + (1 - m_{asv}) \ast \Delta h_{FA}\tag{Equation 0-8}$$ | ||
| + | |||
| + | The required driving power of an element area can be determined with the available results via the difference of the specific enthalpy increase of two supporting points (equation 0-8) via equation 0-9. | ||
| + | |||
| + | $$\Delta \Delta h_{AS/S} = \Delta h_{AS/S_n} - \Delta h_{AS/ | ||
| + | |||
| + | $$P_{D_{AS/ | ||
| + | |||
| + | ===== Power calculation of a single-stage compounding process ===== | ||
| + | |||
| + | In the single-stage compounding process, two polymers are metered into the hopper and then compounded. This means that the components are dosed in a certain ratio in the hopper and simultaneously plasticized. As a result, both materials have the same temperature at one point along the extrusion process. However, the polymers differ in their melting behaviour. At the same temperature, | ||
| + | |||
| + | In the solids transport zone, the enthalpy level of the individual polymers can be calculated using equation 0-1. These are then weighted with equation 0-11 and added together to obtain the enthalpy level at a supporting point. | ||
| + | |||
| + | $$\Delta h_{FF_{C1}} = \frac{\dot{m}_1}{\dot{m}_{SS}} \ast \Delta h_{FF1} + \frac{\dot{m}_2}{\dot{m}_{SS}} \ast \Delta h_{FF2} + \cdots + \frac{\dot{m}_n}{\dot{m}_{SS}} \ast \Delta h_{FFn}\tag{Equation 0-11}$$ | ||
| + | |||
| + | Subsequently, | ||
| + | |||
| + | The same method is used in the melting zone and the melt conveying zone. The enthalpy levels of the polymers are considered separately in the first step and then added together. Equations 0-6 and 0-7 of the basic model provide the enthalpy levels of the corresponding melt and solid content for the individual polymers. In equation 0-8 the proportions are weighted with the enthalpy level. For a process with several polymers that melt simultaneously, | ||
| + | |||
| + | $$\Delta h_{AS/ | ||
| + | |||
| + | Subsequently, | ||
| + | |||
| + | ===== Power calculation for a two-stage compounding process ===== | ||
| + | |||
| + | At the beginning of the process, first one material is plasticized and later in the direction of extrusion two or more materials are plasticized simultaneously. As a result, different numbers of polymers are present at different support points in the extrusion process, which determine the enthalpy level. For this reason, the process must be regarded as a separate application. For the modelling of a two-stage process, it is again possible to fall back on the models already created. Before the calculation, | ||
| + | |||
| + | ===== Power model for the compounding of fillers ===== | ||
| + | |||
| + | During compounding, | ||
| + | |||
| + | |||
| + | ===== Power model for a melt extruder ===== | ||
| + | |||
| + | In the melt extruder, the specific increase in enthalpy due to melting of the solid is eliminated. The melt conveyed in the extruder only undergoes a change in enthalpy due to a temperature variation. For this reason, the calculation of the enthalpy increase up to the melting point can be neglected and the reference point for the enthalpy change is not the filling temperature but the crystallite melting temperature. With later differentiation of the enthalpy changes this portion would be shortened out again anyway. The specific enthalpy change of a melt at a support point in the melt extruder in relation to the crystallite melting temperature can be calculated using Equation 0-6 with Equation 0-13. | ||
| + | |||
| + | $$\Delta h_{SE} = cp_0 \ast (T_M - T_K) + \frac{cp_m}{2} \ast (T_M^2 - T_K^2)\tag{Equation 0-13}$$ | ||
| + | |||
| + | For the driving power, the difference of the specific enthalpies is calculated again (equation 0-14): | ||
| + | |||
| + | $$\Delta \Delta h_{SE} = \Delta h_{SE_n} - \Delta h_{SE_{n-1}}\tag{Equation 0-14}$$ | ||
| + | |||
| + | Equation 0-15 applies to the drive power: | ||
| + | |||
| + | $$P_{D_{SE}} = \dot{m} \ast (\Delta \Delta h_{SE} + \frac{\Delta p}{\rho}) - \dot{Q}\tag{Equation 0-15}$$ | ||
| + | |||
| + | |||
| + | ===== Power model for a melt extruder with several polymers ===== | ||
| + | |||
| + | In the case of a twin-screw extruder, which functions as a melt extruder and is equipped with a melt consisting of several polymers, the model must be modified for a melt extruder. The procedure is similar to that for compounding. The individual enthalpy levels of the components are determined in a first step and then weighted with the aid of the mass flows. Formally, equation 0-16 is obtained. | ||
| + | |||
| + | $$\Delta h_{SE_M} = \frac{\dot{m}_1}{\dot{m}_{ges}} \left(cp_{01} \ast (T_M - T_K) + \frac{cp_{m1}}{2} \ast (T_M^2 - T_K^2)\right) + \frac{\dot{m}_2}{\dot{m}_{ges}} \left(cp_{02} \ast (T_M - T_K) + \frac{cp_{m2}}{2} \ast (T_M^2 - T_K^2)\right) $$ | ||
| + | $$+ \cdots + \frac{\dot{m}_n}{\dot{m}_{ges}} \left(cp_{0n} \ast (T_M - T_K) + \frac{cp_{mn}}{2} \ast (T_M^2 - T_K^2)\right)\tag{Equation 0-16}$$ | ||
| + | |||
| + | Subsequently, | ||
| + | |||
| + | ===== Total drive power ===== | ||
| + | |||
| + | With the relationships shown, the individual required drive powers between two support points or the individual elements can be determined and visualized in SIGMA. Finally, the individual calculated element powers must be added up to obtain the total power. The sole consideration or difference of the initial and final enthalpy to determine the power leads to errors in various cases. This is the case, for example, if the melt temperature drops once during the extrusion process due to fillers and the melt is then heated again. The subsequent heating of the melt must be done by new energy from outside. This requires a further input of power, which would be neglected if the initial and final enthalpy were simply compared. | ||
| + | |||
| + | ===== Validation ===== | ||
| + | |||
| + | The new power model was implemented and verified in SIGMA. Identical to other models, the new power model was compared with the values of the experimental investigations to validate the model. The deviations of different process points and material combinations are shown in the following picture. | ||
| + | |||
| + | {{ :en:grundlagenhandbuch:drehmoment_und_antriebsleistung:en_sigma150_dlg_grundlagenhandbuch_drehmoment_und_antriebsleistung_0010.png? | ||
| + | |||
| + | **Figure:** Comparison of experimental investigations and simulations | ||