Experimental approach to modeling of the plasticizing operation in the hot plate welding process

The paper discusses the topic of butt welding of polyurethane drive belts by the hot plate method in the context of modeling the process of this technological operation. Based on the analysis of the butt welding process, a series of studies of the thermomechanical properties of the material from which the belt is made has been planned. The results will be used for mathematical modeling of the welding process, and in particular its most important phase: the plasticizing operation. On this basis, the study of the compression of cylindrical specimens taken from the belt has been performed at two different speeds. Their result is the relationship between the compressive stress σc and the modulus of longitudinal elasticity Ec at compression and: deformation εc, temperature value T, as well as the compressive velocity vc. In the next step, dynamic viscosity η of the belt material was determined based on the results of dynamic thermomechanical analysis. The research work culminated in the attempts to plasticize the material on a hot plate, in conditions similar to the process of industrial welding. These studies were performed at different speeds vpl, resulting in the correlation between the force required for plasticizing Fpl and the value of the speed of the belt end vpl relative to the hot plate heated to a temperature Tp. The obtained results will be used to formulate a mathematical model of plasticizing the material, based on the selected mechanical deformation models.


Source of the research problem
Endless belts with a circular cross section and a diameter of several millimeters, made of thermoplastic elastomers, are widely used in the driving and transporting systems in industrial machines [1], similar to perforated transporting belts with a rectangular cross section [2]. That is why they are often offered by leading belt manufacturers [3]. The process of their production is two-stage [4]. In the first stage, the belt band with a length of several hundred meters is produced by continuous extrusion. This action ends with its winding on the reel, which facilitates transport and storage. The final step in the production of a finished drive or conveyor belt is to cut to the appropriate length and make a permanent connection of its ends, to obtain a tight belt band with a closed circumference and uniform geometry along the entire length [5].
In many cases, the connection of the ends of the belt is carried out by butt welding using the hot plate method. This method is known as one of the cheapest and the easiest ways of joining thermoplastic materials, both for general industrial purposes [6] and in some specialized industries, for example in medical devices manufacturing [7]. This process is usually carried out manually, using dedicated, simple tools [3]. Unfortunately, this operation usually does not provide a satisfactory repeatability of the geometric dimensions of the product after the connection, as well as high-quality weld. This is mainly due to the manual application of the force needed to plasticize and join the ends of the belt, which does not ensure the repeatability of the conditions for carrying out this operation [8]. To avoid this phenomenon, there is a need for automating this process and a requirement to introduce special control methods. In industrial practice, for example, these ideas can be realized using the non-contact heating operation [9] or by restriction of the displacement of welded elements during their movement [10]. There is also a third solution-controlling only the displacement of welded elements. This can be realized in the easiest way in industrial practice.
In response to the demand of the industry for the automation of the process of butt welding of drive belts, the authors have proposed a design solution for an automatic belt welding machine. The correct configuration of the control system requires the development and verification of technological parameters of the welding process. This is essential for: • properly describing physical phenomena which occur during welding [1], which will allow performing the process in the most effective way possible, • properly controlling the mechanical equipment in the designed machine [4], • properly planning the course of technological operations accompanying welding, e.g., cutting the belt to the appropriate length or cutting the flash [5].
In this context, research work has been undertaken on butt joint welding of drive belts, in particular on the operation of plasticizing the material. Its ultimate goal is to develop a mathematical model of this process.

Current state of knowledge about hot plate welding
Research on butt welding of plastics has already been carried out many times in various aspects. Watson [11] investigated the process of welding plastics in the form of tanks made of polypropylene, polystyrene or polyphenylene oxide. They showed the correlation between the strength of the joint and the technological parameters of the bonding process in an experimental way. Results of their research provided information about the significant positive influence of heating time on the strength of the joint, especially for polypropylene. In addition, influences of heating pressure, heating displacement and consolidation pressure are ambiguous, and strongly depend on the material. Polyphenylene oxide can be characterized by the most stable response for changing the process parameters. Gehde [12] conducted research on the welding of perpendicular beams made of different types of polypropylene in the context of the analysis of joint defects formed during welding with different parameters. During studying the conclusions of their work, we can learn that bond failure mechanisms strongly depend not only on the type of material but also on the welding parameters. Particularly noteworthy is the fact that increasing the joining pressure could give a strength-increasing effect, but the excessive value of this parameter can cause a crack in the joint cross section. Stokes conducted advanced research on plastics welding, which was based mainly on determining the correlation between the relative strength of the welded joint and the process parameters. He conducted research on polycarbonate [10]. As the results of this work, we can find a lot of information about weld strength and failure strain depending on the hot plate temperature, heating and sealing time, and state of sample preparation status. The research was conducted using the hybrid approach: pressure-displacement-controlled welding, using mechanical stops to restrict maximal displacement of material during operation. Research results are a kind of database of welding parameters for rigid polycarbonate. In other works, this researcher obtains similar results for: polycarbonate, polyetherimide, and poly(butylene terephthalate) [13]. Using the same methodology, Stokes obtained a kind of database of weld strength and failure strain for a lot of standard polymer materials (e.g., ABS, PA, PEI, etc.) [14]. It is worth noticing that these studies were carried out with an innovative method of process control-pressure-controlled with a limitation of the movement of the welded element. The second group of his main achievements concerned the experimental comparison of the results of the hot plate welding and vibration welding of rigid thermoplastics: polyamides [15], ABS [16], and poly(vinyl chloride) [17]. Concluding these three papers which are also kinds of databases, we can say that for these materials, using vibration welding we can obtain joints with higher relative strength (where average strength can achieve almost 100% of bulk strength), than using the hot plate welding (for which joint strength can achieve even more than 100%, but average values are much lower). Potente also research about the hot plate welding process, but in this case, his work mainly concerns mathematical modeling and process optimization, using thermodynamic and rheological models. In [18] authors showed derivation of the mathematical model for one-dimensional hot plate welding operation of rigid, thermoplastic material. Their model was validated by the experimental procedure. This work shows a very important issue related to hot plated welding: the complex nature of thermal interaction between hot plate, welded material, and environment. First, there is conduction between the welded material and the hot plate and also between the welded sample and grips used for holding the sample. For the second, there is convection between the hot plate, sample, and the environment. The presence of this phenomenon is due to the atmospheric environment. At last, there is radiation from the sample and the hot plate. Researchers proved that these actions have influence on final temperature distribution during the process. This parameter in turn has a great influence on the process and final joint strength. These works were continued in [19] with the added stress calculations. This author also carried out research on the formation of the effective geometry of the contact surfaces in joined elements before the connection process [20], in which it was proved that shaping of the head surfaces of joined elements does not give any positive effect on the final weld strength. Almost the same results were obtained by Liu [21], who gained an almost constant correlation between relative weld strength and interface shape. Using the same relations as in the model derivation [18], Potente also worked on scaling-up thermodynamic and rheological laws to simplify the hot plate welding model. As results, there is a model for rigid materials which uses non-dimensional coefficients [22]. Potente also conducted research about optimization of the hot plate welding process in different aspects: using computer-aided design tools [23] and for use in self-optimizing heated-tool welder [24]. In cooperation with Schneiders, he also worked on high-speed hot plate welding using a designed machine [25]. All of his impressive work can be characterized by three main features: • he proved that considering all of the thermal phenomena (conduction, convection, and radiation) is important for proper modeling of this process, • all obtained models are used for welding rigid plastics (ABS, PC, PMMA, PE-HD), • due to this fact, the general principle of the welding operation, which is modeled for this applications, is based on a pressure-controlled mechanism, without controlling the displacement.
Nonhof [26] and Cocard [27], similarly to Potente worked on optimizing the welding process for pressure-controlled welding. Poopat conducted a lot of work connected with the hot plate welding process, referring to contact [28] as well as contactless [29] ways to carry out this technological operation. From his experimental considerations, the characteristics of joint strength vs. melt layer thickness, and melt layer thickness vs. heating time and temperature can be obtained. These papers can be treated as a database of sorts for welding HDPE using pressure-controlled hot plate welding. Mokhtarzadeh also carried out a number of research activities in the area of butt welding of plastics, including analysis of welding and comparison of other types of welding with the hot plate method, in terms of quality and strength of the joint, for two similar materials [30] as well as for two different materials [31]. Riahi conducted research on plasticizing polyethylene [32] also using analytical and numerical modeling. The main goal of his research was to estimate the strength of the welds in the function of welding time and the hot plate temperature. Oliveira conducted research on plasticizing polypropylene in order to estimate the size and structure of the flash, depending on the process parameters, both for monolithic [8] and for glass-reinforced polypropylene [33]. A lot of researchers conducted work concerning one of the most important operations from the hot plate welding cycle-material plasticizing. Wood [34] attempted to model the operation of plasticizing a rigid polymer using flow continuity equations and a classical thermodynamic approach. Yoo [35] obtained the model of plasticizing the rigid pipes using the finite element method including Galerkin's solutions. Ezekoye and Nieh using the theory of polymer chain diffusion [36] and thermodynamic similarity numbers [37], gained simplified models of plasticizing and prediction of the weld strength. Lee [38] experimentally obtained the correlation between the degree of plasticity of the material and the process parameters, using some specific methods, e.g., differential scanning calorimetry or thermogravimetric analysis. During literature analysis, it can be noticed that the heat exchange during the hot plate welding process is a very important issue. Savija [39] and Kim [40] analyzed the contact thermal resistance between two surfaces, especially for hot plate welding purposes. Myers [41] and Poslinski [42] investigated the heat exchange between a hot plate and a heated polymer in the context of hot plate welding using the heat exchange laws.
Analysis of the current state of knowledge about the technological process of hot plate welding indicates that the research work undertaken so far relates to butt welding of classic plastics with thermoplastic properties and relatively high rigidity at standard environment temperature (polyamide, polycarbonate, polyethylene and polypropylene). In addition, the cases of butt welding cited in the literature mostly concern elements with a relatively complex crosssectional shape, as a result of which their susceptibility to buckling is low. Therefore, it is common to model the butt welding process as a set of compression operations of the welded material under variable temperature conditions. These phenomena occur in the first stages of the process, when the material is plasticized, and then when the melted surfaces are pressed together. The natural result of this approach to the process description is to analyze the force (or pressure) of material pressure in the various process phases. This approach is common in the works of many authors , regardless of the deviations in the final, practical implementation of the process, which consist, for example, in the use of displacement limiters of the plasticized element [10,[13][14][15][16][17]. This translates directly into industrial applications where it is common to use force or pressure adjustments to maintain fixed, well-defined values of these parameters in the various process stages. much better solution in this case seems to be to control the movement and the speed of the belt end in the various stages of the welding process. What is meant by this is maintaining a constant value of velocity of the ends of the belt in the various stages of the process and implementing the specified displacements by the working units of the device.
A prototype of an automated welding machine using this process control method has been developed [43] taking into account technical and utilitarian conditions in industrial practice (Fig. 1). In consequence, the main features of the designed machine are as follows: • controlling the length of the belt at each of the process stages to obtain the precisely defined final length of the belt, • using mostly the electrical components, rather than pneumatic drives, • maximal simplification of the construction, using as few measurement components as possible.
Due to these assumptions, the main principle of the designed automatic welder consists of controlling the position and velocity of the movable belt gripper (2a), which has one end of the belt (1) clamped. This part is driven by a power screw mechanism and it is responsible for making the required displacement of the end of the belt and pressure to Fig. 1 The prototype of the automated hot plate welder which was introduced to industrial production of the belts; a the whole device during the belt dosage process, b the focused view to the welding unit during plasticizing of the belt: 1-belt, 2a-movable belt gripper, 2b-fixed belt gripper, 3-hot plate the hot plate (3). From the other side, the hot plate (3) is also driven by the same mechanism, but with a velocity that is two times lower. So, in this case, simultaneous movement of the movable gripper and the hot plate, with different velocities, can be obtained. This phenomenon allows obtaining pressure from the hot plate to the second end of the belt (1) which is clamped by a fixed belt gripper (2b). Therefore, simultaneous pressure between the ends of the belt and hot plate during the plasticizing process can be obtained. This action is controlled only by displacement or velocity. The same situation is during the pressing stage when the hot plate is retracted.
This approach makes it possible to obtain the main functional features as per the customer requirements: • obtaining the precise value of the length of belt due to displacement control, • simplifying the structure of the control system and device design, as only one drive unit is used also for hot plate and belt displacement, • exclusion of the possibility of the lack of precision in controlling the welding process, which could occur during controlling of the small-value forces in case of the pressure-controlled system. This phenomenon can occur due to some value of measurement system errors, which can have a big percentage influence in the case of small forces value, • elimination of the stick-slip effect, which could occur in the presumed pneumatic drive system, in case of small inlet pressure (in case of small required forces).
In this case, in contrast to the classic approach [19,44], where the welding process is divided into 5 phases characterized by the fact that the individual operations include control of the pressing force (or pressure), a new method of process handling based on speed and displacement control was proposed for the description of belt welding, while the overall process was simplified to 4 phases (Fig. 2).
During the first stage-plasticizing of the belt (Fig. 2a), the hot plate (3) is pressed against the end of the belt (1b) mounted in the fixed jaw (2b), as a result of the plate (3) moving at controlled velocity v pl towards this part of the belt. At the same time, the other end of the belt (1a), mounted in the movable jaw (2a), is pressed against the hot plate (3) from the other side, as a result of its velocity relative to the fixed end (1a) being 2 × v pl , while this value relative to the hot plate (3) is v pl . This ensures that both ends of the belt are evenly pressed against the hot plate during the plasticizing process. However, it should be noted that the value of the plasticizing force F pl is a variable value, not known and not regulated. Its value depends on the type of material and actual temperature state, which has influence on the obtained force F pl value.
The application of the flat front surface of the belt (1a and 1b) to the hot plate (3) causes a non-linear increase in its temperature along the axis [45]. In this phase, the flash is formed, as a result of early plasticizing and melting of the belt (1a and 1b) on the hot plate (3). Partial melting of the surface causes the surface of the belt ends (1a and 1b) to match the flat surface of the plate (3). This converts the localized heat conduction between the hot plate and the vertices of the belt end surface irregularities into a contact conduction present on the entire surface. This improves the efficiency of the heating process and has a positive effect on the quality of the connection, since the parallel alignment of the interconnected belt surfaces is ensured (Fig. 2a) [1].
When plasticizing the end of the belt (1a and 1b) on a hot plate (3), physical phenomena associated with heat exchange can be observed: conduction, convection and radiation [18,19,27,29,[43][44][45]. Taking into account the individual heat fluxes, one can distinguish ( Fig. 2a In this case, the dominant modes of heat exchange that have an obvious effect on the increase in belt temperature are: contact conduction (Q p3-1 ) between the hot plate (3) and the belt (1a and 1b) and thermal conduction (Q p1 ) inside the belt material (1a and 1b). The thermal radiation (Q r3 i Q r1 ) is not as significant as the conduction (Q p3-1 i Q p1 ) because the process temperature is relatively low and usually does not exceed 300 °C [3,43]. Therefore, it is possible to distinguish two characteristic geometric dimensions: • length p at which the temperature of the belt in its central axis exceeds the welding temperature T w (the temperature at which the material is plasticized and ready to form a permanent joint), • the length h at which the temperature in the central axis of the belt exceeds the environment temperature T 0 .
In the phase of removing the heating device from the welding area (Fig. 2b), the hot plate (3 in Fig. 2a) being removed from the area between the ends of the belt (1a and 1b). At this stage, the heating process of the belt ends is completed, and at the same time, the effect of cooling them as a result of the influence of the gas surrounding them is noticeable. It is important that the temperature of their flat surfaces does not decrease to a level lower than the welding temperature of T w . Therefore, taking into account a certain non-zero time needed to remove the hot plate (3), the temperature of the hot plate T p , is slightly higher than the required welding temperature T w .
In the pressing phase (Fig. 2c), the movable end of the belt (1a) is pressed against the fixed one (1b) at velocity v j to obtain a connection. As a result of the fact that the only Fig. 2 Diagram of the hot plate welding process of the drive belt divided into 4 phases and marked by the expected temperature distribution in the belt axis: 1a-movable end of the belt, 1b-fixed end of the belt, 2a-movable belt holder, 2b-fixed belt holder, 3-hot plate; v pl velocity of belt material plasticization, v j velocity of pressure during joining, v c velocity of pressure during cooling, F pl pressure force during belt plasticization, F j pressure force when connecting the belt ends, F c pressure force during cooling of the connection, T p temperature of the hot plate, T w welding temperature, T 0 environ-ment temperature, p belt plasticizing distance, h section heated to a temperature higher than environment temperature; Q p3-1 heat conduction between the hot plate and the belt, Q p1-2 heat conduction between the belt and the holder, Q p1 heat conduction inside of the material, Q r1 heat transmitted to the environment from the belt by radiation, Q r3 heat transmitted to the environment from the hot plate by radiation, Q c1 heat transmitted to the environment from the belt by convection, Q c3 heat transmitted to the environment from the hot plate by convection parameter controlled in this case is the velocity v j , the value of the compressive force F j is unknown and unregulated. The course of this phase is extremely important for the correctness of the entire welding process, since the interactions between the ends of the belt begin there: • chemical reactions between polymer macromolecules, • mechanical interactions between macromolecules, consisting in their splicing, which facilitates the occurrence of chemical reactions [6].
In the connection cooling phase (Fig. 2d) the movable end of the belt (1a) does not move relative to the fixed one (1b), velocity of the compression v c ≅ 0. However, it should be noted that the belt ends (1a and 1b) are pressed together by the residual force F c remaining after the compression process (Fig. 2c). During this phase (Fig. 2d), the connection cools down, as a result of heat exchange with the environment. Ideally, the connection should be kept in this state for a period of time allowing its temperature to be lowered to T 0 . The main assumptions of this phase are: • the temperature of the belt and the connection is equalized due to thermal conductivity in the material structure (heat Q p1 ), • chemical reactions and mechanical interactions between macromolecules continue, • cooling of the connection occurs as a result of heat exchange with the environment by convection (Q c1 ) and radiation (Q r1 ).
After cooling, the connected belt can be further processed, for example by removing the flash.
Compared to the classic description of the welding process [19,44], the modified welding process is characterized by the fact that in the removal phase of the heating device ( Fig. 2b), the belt ends (1a and 1b) are not far away from the hot plate (3). It is retracted during its contact with the belt material, which provides improved energy efficiency. The plate is covered with the material that provides a low coefficient of friction in combination with the plasticized material. Additionally, during cooling of the connection (Fig. 2d) velocity of movement of the movable end of the belt (1a) v c ≅ 0. Therefore, the external force causing compression during cooling (F c in Fig. 2d) is also close to zero. This phenomenon is purposeful, taking into account the high susceptibility of the belt to buckling-as a result of which a not completely cooled joint can be damaged.
The analysis of the hot plate welding process points out that from the solid mechanic's point of view, the main physical phenomenon that occurs during process flow is compression of the material at elevated, but variable temperature. In addition, it should be noted that one of the most important stages of the process is the plasticizing of the belt on the hot plate. The course of this stage determines the course of the rest of the welding process.
The dependence of the controlled velocity can be showed in the process cyclogram (Fig. 3). Contrary to the standard cyclograms from the classical approach to hot plate welding [19,44], the cycle was simplified to four phases. Using this assumption, a change of only one parameter is required, which allows authors to simplify the control system and constructional solutions of the machine.
A review of the state of knowledge of the described welding process also provides information on the lack of a clear approach to modeling the welding process in terms of the constant movement velocity during these technological operations. In addition, in all analyzed publications, the complex model of the hot plate welding operation, in the sense of dependency between process parameters, e.g.: is not present. The research carried out so far focuses on plasticizing under conditions of a certain-adjustable pressing force, rather than displacement, and they are focused only on chosen aspects (e.g., temperature distribution, stress Fig. 3 Cyclogram of the hot plate welding process in designed welder: I-plasticizing, II-switch over, III-pressing, IV-cooling of the weld, v pl velocity of belt material plasticization, v j velocity of pressure during joining cracking, contact thermal resistance, etc.). It is, therefore, advisable to develop a mathematical model, in particular of the process of belt plasticizing under conditions of adjustable velocity and displacement.

Material conditions
Drive belts with a circular cross section are usually made of thermoplastic elastomers [1,3]. They belong to a group of materials that can be subjected to very large deformations and after being unloaded they return to their original shapes [46]. One of the most commonly used materials in drive belts is a thermoplastic elastomer based on polyurethane. For further analysis, the material with the trade mark TPU C85A [3], which is a case study in the analysis of hot plate welding of drive belts, was adopted.
In general, polyurethanes have a segmented structure (Fig. 4). Their hydrocarbon chain consists of alternating flexible segments (2): methylene, ester or ether, as well as rigid segments (1): urea, urethane or aromatic. Rigid segments give the material high values of strength parameters. Flexible ones, on the other hand, greatly affect its ability to deform. The segments do not mix with each other, forming a two-phase, heterogeneous structure ( Fig. 4) [47].
The macromolecule chains of the polyurethane thermoplastic elastomer are reticulated by means of physical hydrogen and van der Waals forces, which undergo reversible decomposition at elevated temperature. At the basic temperature of use, the material exhibits a highly elastic state, where reversible deformations up to several hundred percent caused by relatively small loads result from the possibility of stretching the chains of macromolecules under the influence of load. Such properties of the structural network of this type of polyurethane cause that on a macroscopic scale this material exhibits properties from the intersection of elastomeric and thermoplastic plastics, which allows its repeated thermal processing, including extrusion, injection or welding [48].
A number of studies have already been carried out on polyurethane thermoplastic elastomers in order to determine their mechanical characteristics. Qi [49] and Boyce [50] conducted research on these materials to formulate their constitutive model, taking into account the elastic hysteresis loop during cyclic deformation. Similar work was carried out by Kukla [51], Drozdov [52] and Wang [53]. This work is based on the analysis of elastic strain energy, which is typical for models of hyperelastic materials, but does not take into account the deformation at elevated temperature. Diani [54] described the elastic deformation of thermoplastic elastomers in a novel way, but also based on models of hyperelastic deformations. Das [55] developed a constitutive model based on the temperatures of two characteristic thermodynamic points for rubber-like materials. Eberlein [56] described the properties of the polyurethane thermoplastic elastomer TPU with a model of a non-linearly elastic-plastic material, but only in the range of standard environment temperature. Henze [57] and Yildrim [58], focused in turn on the description of the mechanical properties of TPU, in the context of modifying its microstructure. However, little information can be found on the description of the properties of the final product-drive belts. Krawiec [59] conducted research on the distribution of temperature influence on the thermal distribution of polyurethane V-belts, but these studies do not provide information on their mechanical behavior under compression conditions at elevated temperature.
In view of the unsatisfactory results of the review of the state of knowledge, the authors have undertaken research activities to determine the selected thermomechanical properties of the welded belt material and its reaction to heating. Based on their results, the following conclusions can be drawn in terms of the phenomena occurring during hot plate welding, and in particular the process of plasticizing the material: • the material does not show a clear melting point, its behavior at elevated temperature is similar to that of amorphous bodies. As a result, it is not possible to indicate a clear plasticizing temperature [60], • heating the belt with a hot plate is a transient process [45], • the material exhibits strongly non-linear characteristics of the stress response to forcing in the form of longitudinal or transverse deformation [61,62].
These conclusions are the basis for modeling the process of belt plasticizing on the hot plate in order to estimate the effective displacement velocity of the belt end relative to the hot plate. With this information, it is reasonable to undertake research activities and refine what is known about the thermomechanical properties of this material in the context of modeling the developed welding methodology.

An approach to modeling the plasticizing process
Developing the mathematical model of this operation is required to develop a set of parameters for the control Fig. 4 Example structure of polyurethanes [47]: 1-rigid segment, 2-flexible segment, 3-crosslinking system. To prepare the mathematical model of the plasticizing process, it is necessary to take into consideration the following facts: • the welded material has specific properties, which is shown by the features from the borderline of elastomers and thermoplastics [46,47]. That is why during the preparation of the mathematical model it is necessary to take into account both elastic and viscous behavior during the compression. It should be noted that these features are also strictly connected with plasticizing velocity, • compression during the plasticizing operation is carried out at an elevated temperature, where the temperature distribution is a function of technological parameters. This phenomenon causes local changes in the properties of the material. What is more, the values of process temperature exceeds the phase change temperature, so the influence of the temperature on the behavior of the material can be very high. Its observable effect is that the mechanical response for axial forces (which are present during plasticizing) can vary very gradually during the process flow.
To consider these facts during the modeling of the plasticizing process, the authors decided to use the combination of basic material models: Newton's and Hooke's [63][64][65]. This approach allows for considering viscoelastic behavior of the material. Using these elementary components, the Kelvin-Voigt model can be built. Contrary to the simpler situations (e.g., during modeling the properties of material which works at room temperature), for the plasticizing operation, the authors decided to prepare a model with controlled parameters (Fig. 5).
The most characteristic feature of this proposed approach to model the plasticizing parameters is the fact, that the elastic parameter (modulus of proportionality E) and the viscous parameter (dynamic viscosity η) will depend on the temperature of the material. Additionally, the ratio between these two properties, which can be described by an additional parameter k, also depends on the temperature. In the authors' opinion, this approach can imitate the material behavior at changed and elevated temperature in a good way (where material at room temperature shows the domination of the elastic properties, but at melting temperature the viscous one dominates). The assumed model can be described by the following formula (1): where is the stress, E is the modulus of proportionality, is the strain, is the dynamic viscosity, ̇ is the strain rate, k is the coefficient of proportionality between the elastic and the viscous component, which can be described by the function (2): where T is the temperature, x is the distance from the hot plate, measured along the axis of the belt (Fig. 2).
To properly describe the planned plasticizing model, considering the current state of knowledge (especially the issue of material features), it is necessary to make a series of thermomechanical tests of the properties connected with the plasticizing operation.
As a result, the following tests were planned: • the compression test of the belts at elevated temperature to obtain the proportionality modulus E c , • the dynamic thermomechanical analysis for evaluation of the dynamic η viscosity parameter, which makes it possible to fulfill and calculate the mathematical model, and • the designation of the plasticizing force F pl as a function of time t or displacement s at a constant plasticizing velocity v pl , which makes it possible to obtain the verifiable parameter for the model derivation.

Belt compression tests at elevated temperature
First of all, axial compression of cylindrical samples of the drive belt material was tested at two values of the compression velocity v c at a variable temperature T. The tests were carried out on the MTS Insight 50 strength testing machine, with a built-in climatic chamber (Fig. 6).
Compression tests were performed to determine the compressive stresses σ c in the function of deformation ε c , depending on the temperature T and the compressive velocity v c . Prior to testing, the samples were conditioned in a climatic chamber, in which a compression test was then  (Table 1), the values of which were selected on the basis of standards PN-EN ISO 604 and PN-80/C-04246. 5 testing cycles were performed for all parameters. The graphs with the correlation between the compressive stress σ c and the deformation ε c were recorded.
A compression plate was used for testing (Fig. 7), equipped with a cylindrical notch with a diameter of 25 mm and a depth of 1 mm, which prevented the slippage of the specimens, without hampering the deformation of their side surfaces. The front surfaces of the specimens were wetted with silicone grease designed to operate at temperatures up to 200 °C. This maintains the axial stress state of the specimens during compression.

Determination of dynamic viscosity of the belt material
The dynamic thermomechanical analysis of DMTA was then tested on the Anton Paar MCR302 rotary rheometer. For testing, specimens were taken from the belt with a circular cross section and a diameter of d = 18 mm, subjected to machining to the form of rectangular beams with dimensions 53.50 ± 0.31 × 10.21 ± 0.20 × 3.74 ± 0.33 mm.
The testing method consisted in introducing the orthogonal beams in controlled oscillations, by inflicting torsioninducing deformations, under varying temperature conditions. The temperature T change cycle included a single heating. A similar methodology was used commonly in research work concerning thermoplastic materials properties. One example is research on self-healing on polyurethane elastomers [66], which depends on manipulation on hydrogen bonds. The DMTA methodology was used for the verification of thermomechanical properties of the material, especially for glass transition temperature. Another example of using the DMTA test for the evaluation of polymeric properties was conducted for modified thermoplastic material [67] when damping factor from DMTA test was used for evaluation of the changes obtained by modification. Test results, especially components connected with a modulus Fig. 6 The workstation for testing thermomechanical properties of drive belts with MTS Insight 50 strength testing machine and climatic chamber during testing: 1-strength testing machine, 2-climatic chamber, 3-climatic chamber controller, 4-compression grips, 5temperature sensor, 6-specimen before test of elasticity G and damping factor δ allows designating the dynamic viscosity [65]. Its value will be necessary to develop the planned model of plasticizing the belts during hot plate welding. The study had the following assumptions: • temperature range from − 30 °C to + 180 °C, recorded continuously, • velocity of temperature change |∇T| = +5 • C∕ min |nabla |∇T| = +5 • C∕ min, • oscillation amplitude A = 0.02%, • oscillation frequency f = 1 Hz.
The measuring system of the device recorded the response of the sample to the forced oscillations and the current temperature conditions. Based on the obtained results, the dynamic viscosity η at torsional loads was determined.

Empirical verification of the plasticizing process
The last of the performed tests consisted of axial plasticizing of the belt end on a hot plate, under conditions that simulate industrial welding of belts. The main assumption of the tests was a constant velocity of the belt movement relative to the heating device v pl , over the entire range of plasticizing time t, with the tests being carried out at its different values. The tests were performed on the MTS Insight 50 kN strength testing machine (Fig. 8), whose jaw (1) and handle (2) pressed the end of the belt (3) against the heating plate (4). This element (4) was heated using Multi-TC, the 70 W electric welding machine for drive belts, manufactured by BEHABelt. The plate (4) had its lower surface adjoining to the thermal insulation pad (5), made of a ceramic plate.
Each time a calibration (I) was performed before the measurement (Fig. 9), which consisted in the fact that a sample of the belt (1) with a height h (the free end of which extends from the grip to a distance h 1 ) and a diameter d, mounted in the grip (2), was moved with a constant velocity v i towards the thermal insulation plate (5) with a constant thickness t = 0.5 mm, until physical contact (detected by the strength measurement system of the testing machine). The plate (5) was adjacent to the hot plate (3) lying on the ceramic pad (4). The traverse movement was then stopped and the plate (5) removed.
Next, the actual testing (II) was started, in which the belt was pressed against the hot plate at the working velocity v pl , while simultaneously measuring: the displacement s of the traverse of the testing machine, plasticizing time t, the pressing force F pl and the temperature value of the hot plate T p . 50 testing cycles were made for: 5 velocities of belt to plate pressure v pl and two specimen diameters d. Plasticizing was carried out in a range far exceeding the actual deformation during welding-up to 50% of the length of the free end h 1 . The parameters of the performed tests are listed in Table 2 (see Fig. 9).
The introduction of the calibration phase (I) made it possible to compensate for the inaccuracy of the specimen height dimension h. The use of a thermal insulating spacer limited the convection heating of the specimen in the calibration phase before the actual testing was started.
Analyzing the course of the compression curves (σ c -ε c ) as a function of the temperature T, the expected effect of a decrease in the compressive stress σ c was observed with an increase in the specimen temperature T. Regardless of the  temperature T, in the studied deformation range (ε max = 30%), the compression curve (σ cε c ) exhibits a nonlinear, degressive character, regardless of the compression velocity v c . However, it should be noted that in the initial characteristic range (σ c -ε c ), at a strain ε c not exceeding 5%, a linear segment can be distinguished, indicating a proportional change in the compressive stresses σ c with the strain ε c . This phenomenon occurs independent of the temperature T of the compressed specimen and the compressive velocity v c .

Determination of material viscosity
Based on the DMTA analysis, the average values of the individual components of the transverse modulus of elasticity were determined: the storage modulus G ′ avg and the loss modulus G ′′ avg , as well as the loss factor tgδ avg depending on the temperature T of the material (Fig. 11).

Plasticizing material on the hot plate
The result of the tests of plasticizing the belt on the hot plate belt are the characteristics of the plasticizing force F pl , depending on the displacement of the belt end s and plasticizing time t (Fig. 12).
Analyzing the obtained results, it can be observed that the characteristics of the plasticizing force value F pl , depending on the displacement of the belt end s or time T, demonstrate a variable course, whereas the character of the variability is repeated for all tests. In the first part of the characteristics (I), the course of the plasticizing force F pl is almost constant, and its value is relatively small in comparison with the value Fig. 9 Diagram of the workstation for pressing the belt to the hot plate: I-calibration phase, II-actual testing phase; 1belt specimen, 2-grip, 3-hot plate, 4-thermal insulation pad, 5-thermal insulation spacer, v i initial calibration speed, v pl testing velocity  (Fig. 12). This phenomenon is associated eliminating the slack between the plasticized specimen (1 in Fig. 9), and the hot plate (3 in Fig. 9), after removing the thermal insulation spacer (5 in Fig. 9). Therefore, it should be assumed that this phase of the material plasticizing operation is not relevant in the analysis of the process, and as such can be omitted in further consideration, unlike the other two distinguishable phases (II and III) related to changes in strength depending on the plasticizing parameters (Fig. 12, Tables 3, 4). Table 3 shows the average characteristics of the process of plasticizing a belt with a diameter d = 12 mm on a hot plate depending on the velocity of the operation v pl , while Table 4 shows for a belt with a diameter d = 18 mm. In both cases, the first phase (I) related to eliminating the slack at the testing workstation was omitted.

Discussion
The described results of the experimental studies on the thermomechanical properties of the welded material can be used to estimate the effective plasticizing velocity v pl of the belt material depending on its geometric dimensions. To this end, it is necessary to carry out the analysis of the obtained results, with a view to establish the correlation between the studied parameters and the process parameters, which will be used in further research work to formulate and validate the mathematical model of the belt plasticizing process.
Compression characteristics of the belt material (Fig. 10), regardless of temperature T and velocity v c can be approximated by the 4th degree polynomial function in the following form: where a 1 ….a 5 are constant factors, determined empirically.
Assuming that a 5 = 0 (starting point of the characteristic σ c = 0 MPa and ε c = 0), it is possible to determine the approximated characteristics, with a correlation coefficient R 2 > 0.9991. They were used to compare the compressive stress σ c obtained with the same strain ε c over the entire range of compression curves (ε max = 30%), for the two studied compression velocities v c1 = 1 mm/min and v c10 = 10 mm/ min. A comparative analysis of these values provided information on the relative difference of compressive stresses δσ c for both compressive velocities v c , which was calculated from the following dependency: where σ c1 is the compressive stress at velocity v c1 = 1 mm/ min, σ c10 is the compressive stress at velocity v c10 = 10 mm/ min. Its results are shown in Fig. 13, in the form of the correlation of the compressive stress difference δσ c between the values obtained for both compressive velocities v c , at different temperature values of the specimen T.
A change in the compressive velocity v c results in a difference in compressive stresses σ c , not exceeding 10% over the specimen temperature range T < 175 °C, with compressive stresses σ c assuming higher values for the compressive velocity v c = 10 mm/min (Fig. 10). Thus, it can be concluded that this difference is negligible and in this temperature range, the strength of the material does not depend on the deformation velocity v c .
This characteristic is somewhat different (Fig. 13) for the maximum specimen temperature (T = 175 °C), where compressive stresses σ c at a higher compressive velocity v c are higher by about 30%. This indicates that at higher temperature T, the viscosity of the test material begins to prevail over the properties of elasticity. This results in the dependence of the macroscopic mechanical properties of the material on the deflection velocity, as a result of which the effect of the obtained compressive stresses σ c on the value of the compressive velocity v c becomes observable [63].
Approximating the start curve (0 < ε c < 5%) of the material compression characteristics (Fig. 10) using a linear function, using the least squares method, it is possible to determine the proportionality modulus at compression E c , with the correlation coefficient R 2 > 0.994 (Fig. 14), from the dependency:   Based on the obtained results, the correlation between the modulus of elasticity E c and the temperature T at compression was determined (Fig. 15), for two compression velocities v c .
Analyzing the presented characteristics, it can be concluded that the determined value of the modulus E c decreases with an increase in temperature T, regardless of the compression velocity v c , which is the expected effect. The relative difference between the values of the modulus of elasticity at δE c for both compression velocities can be determined from the dependency: where E c1 is the modulus at velocity v c1 = 1 mm/min, E c10 is the modulus at velocity v c10 = 10 mm/min. As in the case of compressive stress values σ c (Fig. 15), for temperature T < 175 °C, the differences between the determined values of the proportionality modulus at compression δE c are small and do not exceed 5% in this case. An increase in specimen temperature T results in an increased role of the viscous properties of the material, which is demonstrated by an increase in the difference in the value of the proportionality modulus δE c to about 45%, with a change in the deformation velocity v c . The dependency of the proportionality modulus E c averaged over two compression velocities to specimen temperature T can be approximated using the least squares method with a correlation coefficient R 2 > 0.979 to the following dependency: where a 1 ….a 3 are constant factors. This results in the correlation between the averaged proportionality modulus at compression E c_avg_EMP and temperature T (Fig. 16).
The compression tests of the belt material were carried out at a maximum temperature T = 175 °C. Above this limit, the belt specimens are self-destructed as a result of clear material flow, without compressive force. Taking into account the fact that the process of hot plate welding takes place at higher temperatures, it is necessary to know the proportionality modulus at compression E c in a wider temperature range T. Due to the fact that with the assumed test methodology it is impossible to determine the modulus of proportionality E c and compressive Fig. 13 The correlation of the average percentage difference δσ c between the compressive stresses σ c for the two compressive velocities v c1 = 1 mm/min and v c10 = 10 mm/min to the strain ε c , at different temperature values T of the compressed specimen  Fig. 15 The correlation between the longitudinal modulus of elasticity E c and the difference between the determined moduli δE c and the specimen temperature T for two compression velocities v c1 = 1 mm/ min and v c10 = 10 mm/min stresses σ c empirically at temperature T exceeding 175 °C, the obtained variability function of the modulus E c relative to specimen temperature has been extrapolated for a temperature range T extended to 300 °C. By performing this operation for the averaged proportionality modulus at compression E c_avg_EMP , the expected value of the proportionality modulus at compression E c_avg_EXT is obtained in the temperature range T from 25 to 300 °C (Fig. 16).
As can be inferred from the analysis of the mechanical properties during compression of the belt material at variable temperature T, the values of the basic mechanical properties do not depend on the compression velocity v c at temperature T not exceeding approximately 175 °C. This observation allows us to conclude that for this viscoelastic material, in this temperature range, the dominant effect of loading it is deformation, mainly of elastic nature. The viscous behavior of the material is revealed at a higher temperature.
The performed dynamic thermomechanical analysis study consisted in essence of oscillatory twisting of the specimen, at a constant maximum oscillation amplitude A and a constant frequency f. The obtained results (Fig. 11), however, are not satisfactory from the point of view of the analysis of the operation of plasticizing the belt material on the hot plate. This is due to the fact that the characteristics of the constituents of the transverse modulus of elasticity G and the loss factor tgδ, do not directly provide useful information in the context of compressing the material at elevated temperature. It is necessary to carry out an analysis of the material viscosity η, whose value will play an important function in further studies.
Analyzing the test conditions of the dynamic thermomechanical analysis, the maximum amplitude of the dimensionless transverse deformation γ 0 can be assumed as [63,64,68]: Therefore, the change in deformation γ(t) over time is described by the dependency: and change in the deformation velocity γ'(t): where t is the test time, ω is the circular oscillation frequency, described by the dependency: Taking into account the specified oscillation parameters, it is possible to determine the course of deformation γ and the deformation velocity γ' as a function of time t. On this basis, it is possible to determine the maximum deformation velocity γ'(t), obtained during the DMTA test, whose maximum absolute value is approximately 1.25 × 10 -3 s −1 (Fig. 17).
This value is crucial for the analysis of material behavior. The analysis of these velocities is used in modeling the process of molten polymer processing [63]. Taking into account the standard, typical characteristic curves of polymer flow, it can be concluded that the considered range of deformation velocity lies in the first region of the Newtonian flow of polymer material, for which the value of this velocity does not usually exceed 10 1 s −1 [63,64,69]. This means that in this range of deformation velocity γ', it is possible to take into account the laws concerning the Newtonian fluid, that is the behavior of a viscous polymer material such that any stresses τ (σ) depend linearly on the deformation velocity γ'. This determines the conclusion that the total modulus of elasticity G (or E) during cyclic deformations does not depend on the amplitude of deformations γ 0 for the viscoelastic body, = 2 × × f . Fig. 16 The correlation between the averaged modulus of proportionality E c and the difference between the determined moduli δE c and the specimen temperature T for two compression velocities v c = 1 mm/min and v c = 10 mm/min Fig. 17 Characteristics of deformation variability γ and deformation velocity γ' depending on the duration of the DMTA test which this polymer undoubtedly is [64]. Therefore, the studied material can be treated as a viscoelastic body complying with the Newtonian flow law, within the studied range of deformation velocity γ'.
Therefore, further analysis of the obtained results is warranted (Fig. 11) in the context of determining the dynamic viscosity of the material η, based on the components of the transverse modulus of elasticity G [65].
On this basis, it is possible to determine the total averaged, experimentally determined dynamic viscosity η EXP_AVG and its variability together with temperature T (Fig. 18), using the dependency [63,65]: Since the maximum temperature T during the DMTA test was limited to 180 °C, and the welding process takes place at a temperature of up to 300 °C [4], there is a need to extend the range of temperature T for which the dynamic viscosity η will be determined. In this case, for materials with a welldefined microstructure (crystalline, amorphous or semi-crystalline body), extrapolation is typically used for the range of variation of any mechanical parameter (e.g., viscosity η) together with temperature T, using known empirical laws [37,42,63,70]. For partially crystalline plastics, this is the Williams-Landel-Ferry law (WLF) in the following form: where η 0 is the viscosity at reference temperature T 0 , C 1 and C 2 is the are the material constants. In the case of amorphous plastics, Arrhenius law is usually used in the following form: where R is the gas constant, E is the activation energy of the material. Unfortunately, in this case, both the solid material equations WLF (13) and the activation energy of the material (14) are unknown parameters, as is the actual structure of the material. Therefore, in order to extend the range of known dynamic viscosity η, the empirical data were approximated using the least squares method with a correlation coefficient R 2 > 0.9975. Next, the curve was extrapolated to determine the model viscosity η MOD_AVG in the extended temperature range T (Fig. 18). Therefore, the dependency of this parameter to temperature T, takes the following form: where a 1 , a 5 , b 1 , b 2 are constant factors, determined empirically.
It is important, however, that the physical phenomena occurring in plasticizing the material in the process of hot plate welding are completely different from those occurring in the DMTA test. The main difference is the method of deformation of the material, which during welding involves compression, rather than torsion. Therefore, the results of testing the dynamic viscosity η of the material in the transverse direction are not sufficient. It is necessary to determine the viscosity in the longitudinal direction. For this purpose, the obtained deformation velocities should be analyzed (Fig. 19).
Analyzing the deformation of the belt material during plasticizing, it should be noted that the deformation is not uniform over the entire length of the compressed belt, which distinguishes this situation from the model compression process. In this case, the deformation of the material occurs on the section h 1 located between the belt grip and the hot plate. Taking into account the fact that the belt material deforms unilaterally, as a result of a decrease in strength under the influence of elevated temperature (heating occurs unilaterally), it can be assumed that the part of the belt with temperature T = T 0 is perfectly rigid, in comparison with the part whose temperature exceeds T 0 (comparing in particular to the section of the material heated to temperature T w ). Therefore, it can be assumed that during plasticizing only the section of the belt of length h is subject to deformation. Thus, the deformation to which the belt is subjected will take the value: The following dependency is also true:

Fig. 18
Dynamic viscosity variation characteristics η: η EXP_AVG averaged, experimentally determined, η MOD_AVG model, based on extrapolating the curve therefore, the deformation velocity can be described by the following dependency: However, it should be noted that the delivery of heat to the belt material by means of the hot plate (Fig. 2) causing the temperature of the material to rise, will result in different susceptibility to deformation. As a result, the distribution of the deformation velocity ε' along the belt axis will not be homogeneous, and the expected effect is its relation to the temperature distribution T(x) (Fig. 19). In addition, the value of the distance h also depends on the temperature distribution T(x). Without knowledge of this function, it is impossible to determine precisely the distribution of the deformation velocity ε´. However, taking into account the performed tests of plasticizing the material, they occur at a certain maximum speed: the maximum deformation velocity will be: Regardless of the actual distribution of the plasticizing velocity ε´, its value cannot be greater. Therefore, during belt plasticizing under consideration, deformation velocities remain in the Newtonian viscosity range, i.e., below 10 1 s −1 . The viscosity value, which needs to be known in order to model the process, is, therefore, possible to be determined on the basis of the performed viscosity tests in DMTA analysis with assumed parameters.
In order to determine the longitudinal viscosity η e for compressive/tensile applications, Trouton's definition of viscosity can be used, according to which [46,71]: On this basis, using the dependencies (15) and (22) (Fig. 18) it will be possible to adopt the appropriate viscosity value to model belt plasticizing.
Some of the most interesting results were achieved during the attempts at plasticizing the belt material on the hot plate (Tables 3 and 4). In II phase of this operation, regardless of the plasticizing velocity v pl and the band diameter d, a proportional, approximately straight-line and repeated increase in the material plasticizing force F pl can be distinguished with the specimen displacement s relative to the hot plate.
In the further course of belt plasticizing, it is also possible to distinguish the third phase (III), in which the value of the plasticizing force F pl changes in most cases in a nonlinear way (Tables 3 and 4). The characteristic of its course, depending on the plasticizing time t or the specimen displacement s, depends on the plasticizing velocity v pl and the belt diameter d.
Analyzing the course of phase III of plasticizing a belt with a diameter d = 12 mm (Table 3), it can be observed that for velocities v pl equal to 2 mm/min and 4 mm/min, the value of the force F pl decreases non-linearly with an increase in plasticizing time t. In the case of velocities v pl equal to 12 mm/min and 16 mm/min, the value of the force F pl increases non-linearly with an increase in plasticizing time t. An interesting observation is also the course of the plasticizing force F pl at velocity v pl of 8 mm/min, when the force F pl assumes values that are almost constant, independent of plasticizing time t.
Similar dependencies can be observed when the belt with a diameter d = 18 mm is plasticized (Table 4), the difference being that for plasticizing velocity v pl = 8 mm/min, the value of the force F pl increases non-linearly, in contrast to the belt with a smaller diameter d. At other values of plasticizing velocity v pl , similar phenomena were observed as in the case of the belt with a diameter d = 12 mm.
Analyzing the course of phases II and III, it is possible to approximate the obtained characteristics, using the least squares method, with a correlation coefficient R 2 > 0.996, resulting in a diagram of the correlation of the force plasticizing the belt material and time t or displacement s, based on a mathematical description (Fig. 20).
The interpretation of the phenomenon of proportional increase in the value of the plasticizing force F pl in the phase II of the plasticizing operation can be the fact that the process of heating the belt material is not instantaneous. In contrast to the often accepted assumptions about the heating of the material by the hot plate, which consist in modeling this operation as a one-sided heating of a semi-infinite body with the application of a heat source to its only symmetry surface [9,13,18,19], the heat exchange between the specimen and the hot plate under real conditions is accompanied by a certain thermal resistance [39,42,45]. In addition, each material is characterized by a certain thermal inertia resulting from a finite value of the thermal diffusion coefficient a [72].
Therefore, assuming that the plasticizing velocity v pl is constant and non-zero, at this phase, there is a phenomenon of partial compression of the still unplasticized solid, so that an increase in the pressing force F pl with the displacement s is the expected effect. According to the results of the compression tests (Fig. 10), the application of the compressive force F pl to the material in the range of small deformations is characterized by a linear range, which is reflected in phase II of plasticizing the material (Tables 3 and 4). It is obvious, however, that the directional coefficient of this straight a II (Fig. 20) changes its value with the change of velocity v pl and diameter d, whereas the correlation between these values is positive.
In the case of phase III of this operation, the observed change in the characteristics of the course of the plasticizing curve (Tables 3, 4 and Fig. 20) depending on the constant and non-zero velocity v pl and diameter d, is mainly due to the limited material heating velocity. This is related to a finite and approximately constant value of the thermal diffusion coefficient a, approximately independent of plasticizing velocity v pl .
Therefore, in phase III there are 3 types of phenomena, the effects of which can be observed on the basis of the analysis of the course of the plasticizing force F pl (Fig. 20) and the shape of the specimens after plasticizing (Fig. 21 (Fig. 20), as a result of compression of the material subject to plastic flow. Analyzing the shape of the specimen after the plasticizing operation (Fig. 21), no clear boundary can be observed between the core (2) and the flash (1). This is due to excessive deformation of the material, its melting and flashing from the working area. It can therefore be concluded that such plasticizing is excessive and causes unnecessary loss of energy and material. In industrial practice, however, this phenomenon can cause additional design and technological problems associated with the outflow of molten polymer from the working area of the machine, • in the case of a high value of the plasticizing velocity v pl3 (for the belt with the diameter d = 12 mm it is 12 mm/ min and 16 mm/min), the material is under-heated, since the velocity of pressing the specimen to the plate v pl is significantly greater than the velocity of heat dissipation in the polymer structure. As a result, it is insufficiently plasticized, which in principle means the material is being compressed in the solid state. As a result, the plasticizing force increases significantly in the further course of this operation (Fig. 20). The observable effect is the shiny structure of the core (2), from which the plasticized material is completely extruded into the flash area (1) (Fig. 21). This phenomenon is disadvantageous because the core of the material does not have sufficient temperature and degree of plasticity to form a permanent connection. In addition, the value of the technological forces of the process is adversely increased, which also worsens its energy efficiency, • in the case of an intermediate value of the plasticizing velocity v pl2 (for a belt with a diameter d = 12 mm, this is about 8 mm/min), the material is properly plasticized. By analyzing the correlation between the plasticizing force F pl and time t (Fig. 20), an almost constant course of its value can be observed, which indicates that there is a balance between heat conduction and deformation of the material. The observable effect of this phenomenon is uniform melting of both the flash area (1) and the core (2) (Fig. 21). Such a course of the plasticizing operation can be considered effective.
The next very interesting issue is the correlation between plasticizing velocity v pl and plasticizing time t pl . Due to controlling the v pl velocity in the automatic welding machine, and also during the plasticizing tests conducted on the MTS durometer, the plasticizing velocity fulfills the following condition: Therefore, the negative correlation between displacement s and time t is obtained. Obviously, the velocity in the real application is not constant due to the speeding up and braking of the mechanism. As it can be noticed in the process cyclogram (Fig. 3), the shape of the velocity diagram is not rectangular, but trapezoidal. The same situation is present during the test in MTS Insight durometer. However, this phenomenon during process analysis can be omitted in both stages for the following reasons: • during speeding up, because of the initial clearance between the end of the belt and the hot plate; increasing the velocity is realized during gap deleting, • during braking, because of the fact that the required plasticizing distance is exceeded in this time; so, this phenomenon does not have any influence on necessary plasticizing distance.
On the basis of the obtained results (Tables 3 and 4), the analysis of the plasticizing time t pl can be conducted. Consideration of the overall plasticizing time t pl can be pointless, because of the fact that the value of this parameter depends on displacement s which was set. Its value in the practical approach depends on the required plasticized distance h, whose value can depend on, e.g., sizing of the inequalities on the surface of the belt which should be melted. On the other hand, time t pl depends also on plasticizing velocity v pl . It is obvious that increasing the plasticizing velocity v pl gives the effect of decreasing the plasticizing time t pl , considering the constant value of the required displacement s. The essence of the work conducted by the authors concerns formulation of the mathematical model of plasticizing, to find the most effective value of the plasticizing velocity v pl which is necessary to obtain the required plasticization of the belt with possibly lowest energy consumption. Additionally, performing an analysis of the time t pl to obtain a hypothetical correlation between maximal plasticizing force F plMAX and plasticizing time t pl may not give expected results, because of a different plasticizing course in the III phase, for different plasticizing velocity v pl . In the authors' opinion, better results can be obtained by designation of a new parameter-time to plasticizing phase change t II-III . This parameter can give information about the moment of changing the plasticizing phase, from the second to the third. Analysis of this parameter can give an insight into the performance of the plasticizing operation. Its correlation with plasticizing velocity v pl and sample diameter d is given in the Fig. 22.
Analysis of the presented dependencies points out that phase change time t II-III decreases with increasing plasticizing velocity v pl , which is an expected effect. Both obtained characteristics can be approximated using the least squares method with a correlation coefficient R 2 > 0.9996 to the Fourier's model, which can be given by the equation: where a 0 , a 1 , a 2 , b 1 and b 2 are constant factors, determined empirically. The variation of the phase change time t II-III becomes slighter with increasing plasticizing velocity v pl . This is due to the fact that the increase of the plasticizing velocity v pl , causes blurring of the differences between II and III phases. This fact can be observed in the gained plasticizing curves (Tables 3 and 4). The mathematical proof of this phenomenon can be obtained by the analysis of the derivative of the gained function (24), which tends to zero while the value of the independent variable (v pl ) goes to 16 mm/min.
The physical interpretation of this phenomenon may be the fact that increasing the plasticizing velocity gives an effect of insufficient heating, which causes compression of material at almost solid state, rather than pure plasticizing. In this case, particular plasticizing results (for different diameters d) become similar to each other and similar to the results from compression tests of a solid-state material. Due to this, for higher plasticizing velocity v pl , the situation of less influence of sample diameter d on plasticization effects becomes observable. Of course, obtained values of plasticizing force F pl are different, but the shape of the plasticizing graph (rather a quasi-compression process) is similar.
Using the extrapolation method, the prediction of the analyzed function can be obtained (Fig. 22). According to the predicted models (such as before and after curve from tests), it can be observed that time t II-III grows up for smaller and also for higher plasticizing velocities v pl . For lower plasticizing velocity v pl , this effect is expected, because with decreasing velocity v pl , the dynamics of the process is lower so the heating operation lasts longer. For higher plasticizing velocity (in the predicted model for velocity higher than 16 mm/ min in the Fig. 22), the effect of the increase of the time t II-III is unexpected. This fact may result from the changing nature of the process in this range of velocities, from plasticizing to pure compression. However, this fact should be checked in an experimental way in the future.
As can be seen in the studied plasticizing velocity v pl range, the correlations between obtained time t II-III for both diameters d are very similar. What is more, the difference between them is almost constant, which can be proved by obtaining the relative difference between the values of the phase change time at δt II-III for all of plasticizing velocities v pl and two diameters d. Its value can be determined from the dependency: where t II-III (18) is the time of phase change for 18 mm of diameter and t II-III (12) is the time of phase change for 12 mm of diameter. Obtained differences can be presented on a correlation diagram (Fig. 23).
As can be seen, in the obtained dependency, the relative difference between time to phase change for both diameters is almost constant and reaches around 40%. For belts with d = 18 mm the time t II-III is bigger. The (25) (18) × 100% , Fig. 22 Correlation between phase change time ( II → III ) and plasticizing velocity v pl plasticizing velocity, t II-III time of plasticizing phase change physical interpretation of this phenomenon concerns the bigger volume of heated material (of about 56%), which results from higher diameter d. Consequently, this sample demands more energy to increase the temperature. Assuming the constant power of the hot plate, which results in constant heat flux, according to thermodynamics' principles the heating time should be higher. This phenomenon is responsible for increasing the phase change time t II-III with increasing plasticizing velocity v pl [72].
On the other hand, it should be noticed that the outer surface area is mainly responsible for energy losses, which is connected with convection and radiation [18,35,45,72]. Increasing the sample diameter d, results in decreasing the ratio between outer surface area and volume, which can be described by the flowing formula: where d is the diameter of the belt, A is the outer surface area, V is the volume, h is the the unitary height of the sample. As a result, the influence of surface area A is smaller for belts with higher diameter d, so the relative energy losses for bigger diameters d can be smaller. This phenomenon may cause a slight decrease of the time t II-III for bigger diameters d. However, this phenomenon is not balanced with the volume V increase, so the effect of the time t II-III growth is observable.
Despite these facts, the most important issue is that the difference between phase change time t II-III is almost constant for these two diameters d. This observation proves that the heat exchange conditions for both types of samples were almost the same, so the results of the test could be compared.

Conclusions
The conducted research made it possible to identify the basic thermomechanical properties of drive belts undergoing the welding process. Their values will be used to formulate a mathematical model of the belt plasticizing operation, as the most important phase of the technological process of hot plate welding. The determined characteristics: • compressive stresses σ c and proportionality modulus at compression E c from deformation velocity ε c and temperature T, • dynamic longitudinal viscosity η e from the temperature value at the specified deformation velocity, • plasticizing force F pl as a function of time t or displacement s at a constant plasticizing velocity v pl , will allow to formulate a mechanical model of the material during the plasticizing process, using the appropriate combination of basic models: Newton's, and Hooke's, especially the Kelvin-Voigt type with controlled parameters.
To carry out these activities in an appropriate manner, the following activities are planned in further research work: • detailed testing of the belt material plasticizing will be performed, and the refinement will consist in narrowing the velocity range v pl to values that allow obtaining a defined course of plasticizing in phase III of the operation (Fig. 20), for both 12 mm and 18 mm diameter belts. In addition, the results will be expanded with additional values of the hot plate temperature T p , • the temperature distribution T(x) during belt heating will be determined using mathematical models of heat conduction in a rod with a moving heat source [72]. The adopted model will be verified empirically during detailed plasticizing studies, • additional tests of the belt material compression will be performed, at compression velocities v c as close as possible to the belt material plasticizing velocities v pl , The obtained results will make it possible to determine a mathematical model of the plasticizing operation, based on a mechanical model, in which the technological parameters of this operation will be taken into account.