Squeezing Mechanical Analysis and Model Establishment of the Viscoelastic Rubber-Strip-Feeding Process of the Cold-Feed Rubber Extruder

In the process of rubber extrusion, the feed structure directly affects the extrusion quality, extrusion uniformity, screw lateral force, and feed power consumption. Until now, the feed structure was mainly based on empirical designs, and there was no theoretical model for the optimal design of a feed structure. This paper focused on the squeezing mechanical analysis and model establishment of the feeding process in which viscoelastic rubber strips are passed through feed-wedge clearance in cold-feed extruders. The screw flight rotation squeezing process was simplified into a disc rotation squeezing process; the instantaneous squeezing velocity h˙(t) in the disc rotation squeezing model was derived according to feed wedge clearance geometry and the disc rotating speed. By transforming rotation squeezing into differential slab squeezing, mathematical expressions of the velocity distribution, pressure distribution, total squeezing force, and power consumption in the feeding process were derived in a rectangular coordinate system under isothermal and quasi-steady assumptions and certain boundary conditions by using balance equations and a Newtonian viscous constitutive relation. Theoretical calculations and experimental values showed the same trend. Through comparison, it was found that the power consumption (P3) caused by sliding friction is about 200–900 W according to theoretical calculations, while the experimental test results show it to be about 300–700 W. Additionally, the difference between theoretical pressure value and the experimental pressure value can be controlled within 5–15%. This could reflect the main factors that affect the feeding process, so could be used for analyses of actual feeding problems, and to contribute to rough quantitative descriptions of the feeding process, finite element simulation, and the optimization of the feeding structure.


Introduction
The extrusion process is an important step in rubber processing, and the vast majority of rubber compounds are extruded at least once during the molding process [1]. The feeding process has a significant impact on extrusion quality [2,3], extrusion stability, the side force applied to screw, and extruder power consumption. In modern research and applications, in order to realize a uniform feed rate in rubber extruders, rubber strip feeding has been used in either rarely used hot-feed extruders or widely used cold-feed extruders. Because rubber has relatively strong adhesiveness at room temperature in a highly elastic state, the cost of the granulation process will increase if a rubber extruder is fed granulated material. A large number of separating agents can be added to prevent the adhesion phenomenon of rubber granulates, resulting in changes in the formulation properties [4]. For a cold-feed extruder fed with a rubber strip form, regardless of whether a free-feeding structure or the widely used forced-feeding structure with a feed roll is used, the feed wedge clearance is always designed at the feed port. Figure 1a shows the common feeding model (feed roll structure) of a cold-feed rubber extruder. The feed wedge clearance is formed between the screw flight crest and the undercut groove on the inner wall of the feed barrel. The feed wedge clearance is a key structure when feeding the rubber strip into the screw channel: the rubber strip is gradually squeezed (compressed) and thinned when the rotating screw flight drags the rubber strip through the feed wedge clearance. At the same time, the screw flight is more deeply embedded into the rubber strip (Figure 1b). When the rubber strip approaches the minimum wedge gap (the design gap between the outer diameter of the screw and the inner wall of the feed barrel), it is longitudinally broken off [5], thus entering the channels at both sides of the screw flight and completing the feeding process of the rubber strip. In a combined structure of a feed roll/screw system and the feed wedge clearance, the feed roll only increases the pushing force [6] when feeding the rubber strip or creating an additional pressure build-up [7]. Therefore, when using rubber extruders, the rubber strip is dragged through the feed wedge clearance and broken off longitudinally at the end of the feed wedge clearance; this is the key to successfully carrying out rubber strip feeding. In the vicinity of the minimum wedge gap, the thin strip easily causes longitudinal break-off under the axial thrust action of the screw flight. Since rubber is essentially a viscoelastic liquid at room temperature [8], a special kind of squeezing (compression) flow occurs in the feed wedge clearance. The feed wedge clearance geometry is a nonlinear wedge region with a narrow width (the width of the screw flight). The squeezing force is generated by dynamic viscous pressurization in the nonlinear wedge region through the rotation drag motion of the screw flight, and rubber flows in the nonlinear wedge region occur only in the directions towards both sides of the screw flight, because the minimum wedge gap is very small, and can be neglected for rubber compounds with very high viscosity.
Polymers 2022, 14, x FOR PEER REVIEW 2 the feed wedge clearance is always designed at the feed port. Figure 1a shows the com feeding model (feed roll structure) of a cold-feed rubber extruder. The feed wedge c ance is formed between the screw flight crest and the undercut groove on the inner of the feed barrel. The feed wedge clearance is a key structure when feeding the ru strip into the screw channel: the rubber strip is gradually squeezed (compressed thinned when the rotating screw flight drags the rubber strip through the feed w clearance. At the same time, the screw flight is more deeply embedded into the ru strip ( Figure 1b). When the rubber strip approaches the minimum wedge gap (the d gap between the outer diameter of the screw and the inner wall of the feed barrel) longitudinally broken off [5], thus entering the channels at both sides of the screw and completing the feeding process of the rubber strip. In a combined structure of a roll/screw system and the feed wedge clearance, the feed roll only increases the pus force [6] when feeding the rubber strip or creating an additional pressure build-u Therefore, when using rubber extruders, the rubber strip is dragged through the wedge clearance and broken off longitudinally at the end of the feed wedge clearance is the key to successfully carrying out rubber strip feeding. In the vicinity of the mini wedge gap, the thin strip easily causes longitudinal break-off under the axial thrust a of the screw flight. Since rubber is essentially a viscoelastic liquid at room temperatur a special kind of squeezing (compression) flow occurs in the feed wedge clearance feed wedge clearance geometry is a nonlinear wedge region with a narrow width width of the screw flight). The squeezing force is generated by dynamic viscous pres zation in the nonlinear wedge region through the rotation drag motion of the screw f and rubber flows in the nonlinear wedge region occur only in the directions towards sides of the screw flight, because the minimum wedge gap is very small, and can b glected for rubber compounds with very high viscosity. The history of rubber extruders dates back more than 140 years, to the inventi hot-feed extruders in 1879. The continuous development of feeding technology is m reflected in the improvements to the feeding structure [9][10][11][12][13] and the control of fee uniformity [4,[14][15][16]. However, there has been little research into the feeding theo rubber extruders, especially cold-feed extruders. The existing research on the theo the rubber-strip feeding process has either analyzed the mechanical conditions of ru strips being drawn into the feed wedge clearance [6], or used a disc rotation sque model device to carry out experimental research on the strip feeding process [5]. Sinc study by Jianbin Li did not consider the contribution of the flight flanks to the addit drag effect [6], it reached the incorrect conclusion that the feeding mechanical cond of the friction coefficient μs between the flight crest and the rubber strip was greater the friction coefficient μb between the barrel inner wall and rubber strip. Other ex mental results in the studies of Yanchang Liu [5] included single-peak pressure dist tion and the obvious power consumption when the rubber strip passed through the m The history of rubber extruders dates back more than 140 years, to the invention of hot-feed extruders in 1879. The continuous development of feeding technology is mainly reflected in the improvements to the feeding structure [9][10][11][12][13] and the control of feeding uniformity [4,[14][15][16]. However, there has been little research into the feeding theory of rubber extruders, especially cold-feed extruders. The existing research on the theory of the rubber-strip feeding process has either analyzed the mechanical conditions of rubber strips being drawn into the feed wedge clearance [6], or used a disc rotation squeezing model device to carry out experimental research on the strip feeding process [5]. Since the study by Jianbin Li did not consider the contribution of the flight flanks to the additional drag effect [6], it reached the incorrect conclusion that the feeding mechanical condition of the friction coefficient µ s between the flight crest and the rubber strip was greater than the friction coefficient µ b between the barrel inner wall and rubber strip. Other experimental results in the studies of Yanchang Liu [5] included single-peak pressure distribution and the obvious power consumption when the rubber strip passed through the model clearance. To date, little theory of the rubber strip feeding process exists, which makes the feeding theory for rubber extruders lag far behind engineering practice. The side force applied to the screw and the feed power consumption have not been calculated theoretically, and the feed structure and technology have not been improved or innovated by theoretical guidance. In order to establish a feeding process theory of cold-feed rubber extruders, the squeezing flow dynamics during this process must be analyzed. However, most studies of squeezing flow deal with the pure squeezing flow of viscoelastic materials [17][18][19][20][21][22][23][24][25][26][27][28], or viscoelastic materials [29][30][31][32][33][34][35][36][37][38][39][40], or non-Newtonian liquids with a squeezing-extrusion combined flow [41] between two parallel plates; very little research has considered the squeezing flow between two non-parallel plates. The squeezing and sliding flow of Oldroyd-B fluid was examined by N. Phan-Thien [42] in a wedge geometry of semi-wedge angle a(t) with a wedge apex, using a numerical method. In this wedge, the flat plane boundaries of the wedge were closing at a rate . a(t) and were sliding along the direction perpendicular to the two-dimensional wedge with a constant velocity. However, there was no consideration of the squeezing flow of non-Newtonian liquids in a nonlinear wedge region, in which the lubrication approximation was not valid due to large wedge angle or slip boundary condition. The squeezing flow in feed wedge clearance belongs to the latter category, because the rubber strip slides along the cylinder wall in the feeding process. Therefore, in this study, we attempted to analyze the dynamics of squeezing flow caused by the drag action of the screw flights in the nonlinear wedge region under sliding boundary conditions. In order to obtain an approximate analytical solution of the motion equation which reflects the basic characteristics of the squeezing flow, under some assumptions (Newton material, isothermal quasi-steady flow), the physical model was simplified (disc rotation squeezing model) and the motion transformation (differential slab squeezing flow) was applied. To verify the accuracy of the theoretical model, a comparison between the theoretical calculations and the experimental data was carried out. The theoretical model that was established could provide guidance and suggestions for the design and optimization of actual feed structures.

Physical Model and Squeezing Velocity
In this section, a physical model of feeding process is established and the calculation formula of squeezing velocity in the feed wedge clearance is obtained by analyzing the relationship between the rotation angle and the wedge gap.

Physical Model
At the feed port, the basic screw configuration of the cold feed extruder is a single screw section with double-flighted design. When rubber strips pass through the feed wedge clearance, the screw flight crest will always contact and squeeze the rubber strip. If expanded along the outer diameter of the screw, the contact traces between the screw flight crest and the rubber strip are two inclined narrow strips (the black areas in Figure 2). In Figure 2, β is the helical angle and w is the axial width of the flight. To describe the behavior of the feeding process, the geometry of the helix feed wedge clearance is simplified, and then motion transformation is applied, as discussed in Section 3.1. As the rubber strip passes through the feed wedge clearance, the effective flight compression path is perpendicular to the screw axis. Therefore, each feed's screw flight can be simplified to a disc (helical angle β = 0 • ), and the diameter and width of the disc are screw diameter D and screw flight axial width w, respectively. An effective feed wedge clearance is formed between the disc and the undercut groove on the inner wall of the feed barrel. In this way, a disc rotation squeezing model is built with a disc and a barrel ( Figure 3). Because the two screw flights have the same geometry, one was chosen for mathematical description.
For all general specifications of cold-feed extruders, a standardized design of the feed wedge clearance has been achieved. In Figure 3   For all general specifications of cold-feed extruders, a standardized design of the feed wedge clearance has been achieved. In Figure 3, a = 0.03 − 0.04D, b = 0.06 − 0.08D, R = 0.5D, R1 ≈ 0.5D, is the start position or maximum value of the wedge gap (on the extended line of O and O1, that is = ), ℎ is the end position or minimum value of the wedge gap (the typical value is 0.0045R), and is the center angle of the circular arc of the outer diameter of the disc in the whole range of the feed wedge clearance. For example, for a cold-feed extruder with D = 65 mm, a = 2 mm, b = 5 mm, R = R1 = 32.5 mm, H = 5.4 mm, hmin = 0.146 mm, and θ0 = 85° ≈ 1.483 rad.

Squeezing Velocity
An analysis of squeezing flow needs to determine either the squeezing force or the squeezing velocity. In the disc rotation squeezing model, the squeezing force is unknown, and the squeezing velocity can be derived according to the geometry of the feed wedge clearance and the disc speed. Figure 4 shows the relationship between the gap ℎ( ) and the angle .   For all general specifications of cold-feed extruders, a standardized design o wedge clearance has been achieved. In Figure 3, a = 0.03 − 0.04D, b = 0.06 − 0.08D, R1 ≈ 0.5D, is the start position or maximum value of the wedge gap (on the e line of O and O1, that is = ), ℎ is the end position or minimum valu wedge gap (the typical value is 0.0045R), and is the center angle of the circu the outer diameter of the disc in the whole range of the feed wedge clearance. For for a cold-feed extruder with D = 65 mm, a = 2 mm, b = 5 mm, R = R1 = 32.5 mm mm, hmin = 0.146 mm, and θ0 = 85° ≈ 1.483 rad.

Squeezing Velocity
An analysis of squeezing flow needs to determine either the squeezing for squeezing velocity. In the disc rotation squeezing model, the squeezing force is u and the squeezing velocity can be derived according to the geometry of the fee clearance and the disc speed. Figure 4 shows the relationship between the gap ℎ the angle .

Squeezing Velocity
An analysis of squeezing flow needs to determine either the squeezing force or the squeezing velocity. In the disc rotation squeezing model, the squeezing force is unknown, and the squeezing velocity can be derived according to the geometry of the feed wedge clearance and the disc speed. Figure 4 shows the relationship between the gap h(θ) and the angle θ.
If the initial thickness of the fed rubber strip is equal to the maximum value H of the wedge gap, the position of the line segment ss can be set as the initial point at which the rubber strip enters into the wedge clearance. When the disc rotates through an angle θ (rad) from time t = 0 to t(s), the strip is squeezed (compressed) from thickness H to h(θ) (Figure 4). The mathematical relationship between h(θ) and θ is where c = OO 1 = √ a 2 + b 2 . By expanding the term with the root sign on the right side of Equation (1) into a power series, By ignoring the terms greater than the quadratic term of the power series on the right side of Equation (2) and substituting it into Equation (1) and considering R ≈ R 1 , one obtains If the rotating speed of the disc is Nrpm, then θ = πNt/30. Substituting the expression of θ into Equation (3) gives By differentiating both sides of Equation (4), one can obtain the instantaneous squeezing velocity where "−" indicates squeezing (compression). From Equation (5), .
h is directly related to the geometric parameters of the feed wedge gap clearance (c and R 1 ) and the operating parameters (N and t).  If the initial thickness of the fed rubber strip is equal to the maximum value H of the wedge gap, the position of the line segment ss can be set as the initial point at which the rubber strip enters into the wedge clearance. When the disc rotates through an angle θ (rad) from time t = 0 to t(s), the strip is squeezed (compressed) from thickness H to h(θ) ( Figure 4). The mathematical relationship between h(θ) and θ is By expanding the term with the root sign on the right side of Equation (1) into a power series, By ignoring the terms greater than the quadratic term of the power series on the right side of Equation (2) and substituting it into Equation (1) and considering R ≈ R1, one obtains If the rotating speed of the disc is Nrpm, then θ = πNt/30. Substituting the expression of θ into Equation (3) gives To study the change in instantaneous squeezing velocity, according to Equation (5), MATLAB software was used to calculate and draw . h − t curves at different disc speeds (N = 30, 45 and 60 rpm) for the disc rotation squeezing model of D = 65 mm (c = √ 2 2 + 5 2 ≈ 5.39 mm, R 1 = 32.5 mm). Figure 5 shows that the . h − t curve is approximately composed of two parts: a constant acceleration squeeze in the early stage and a constant velocity squeeze in the later stage. The former accounts for most of the whole rotational squeezing process, and the latter accounts for a small part. Setting t 0 and t t to represent the time used during the whole rotational squeezing process and the transition time from the constant acceleration squeeze to the constant velocity squeeze, respectively, t 0 can be calculated according to the formula t 0 = 30θ 0 /πN, and t t can be measured from the asymptotic transition part of the . h − t curve in Figure 5. When N = 30, 45, and 60 rpm, t 0 is 0.47, 0.31, and 0.24 s, t t is 0.39, 0.27, and 0.21 s, and t t / t 0 is 0.83, 0.87, and 0.875, respectively. Figure 5 also shows that the squeezing velocity increases faster with increasing disc speed. Similar results can be obtained by using the same method to calculate other disc rotation squeezing models with different diameters (D = 90, 120, 150, 200, and 250 mm). acceleration squeeze to the constant velocity squeeze, respectively, t0 can be calculated according to the formula t0 = 30θ0/πN, and tt can be measured from the asymptotic transition part of the ℎ − curve in Figure 5. When N = 30, 45, and 60 rpm, t0 is 0.47, 0.31, and 0.24 s, tt is 0.39, 0.27, and 0.21 s, and tt/ t0 is 0.83, 0.87, and 0.875, respectively. Figure 5 also shows that the squeezing velocity increases faster with increasing disc speed. Similar results can be obtained by using the same method to calculate other disc rotation squeezing models with different diameters (D = 90, 120, 150, 200, and 250 mm). The above analysis shows that, for a given disk rotating squeezing model and disc rotating speed, the instantaneous squeezing velocity ℎ ( ) changes rapidly for most of the time taken for the whole squeezing process. This presents difficulties when analyzing the squeezing flow. To simplify the theoretical analysis, the average squeezing velocity (−V) was introduced, where '−' also represents the squeezing direction, = |− |. According to N, H, and θ0, it is easy to calculate V.
Equation (6) shows that, for a given disc rotation squeezing model, the average squeezing velocity V is only proportional to N. For example, for the disc rotation squeezing model with D = 65 mm, the average squeezing velocities calculated by Equation (6) at N = 30, 45, and 60 rpm are −11.4, −17.1, and −22.8 mms −1 , respectively. For the disc rotation squeezing models for other general specifications of cold-feed extruders, when the rotational speed N = 30 rpm, the average velocities calculated by Equations (5) and (6) in the whole feed wedge gap are V(5) and V(6), respectively. The range of [(V(5) − V(6))/V(5)] × 100% is approximately −0.51~1.30%, which strongly indicates that the simplified Equation (5) has a very high level of accuracy. The above analysis shows that, for a given disk rotating squeezing model and disc rotating speed, the instantaneous squeezing velocity . h(t) changes rapidly for most of the time taken for the whole squeezing process. This presents difficulties when analyzing the squeezing flow. To simplify the theoretical analysis, the average squeezing velocity (−V) was introduced, where '−' also represents the squeezing direction, V = |−V|. According to N, H, and θ 0 , it is easy to calculate V.
Equation (6) shows that, for a given disc rotation squeezing model, the average squeezing velocity V is only proportional to N. For example, for the disc rotation squeezing model with D = 65 mm, the average squeezing velocities calculated by Equation (6) at N = 30, 45, and 60 rpm are −11.4, −17.1, and −22.8 mms −1 , respectively. For the disc rotation squeezing models for other general specifications of cold-feed extruders, when the rotational speed N = 30 rpm, the average velocities calculated by Equations (5) and (6) in the whole feed wedge gap are V (5) and V (6) , respectively. The range of [(V (5) − V (6) )/V (5) ] × 100% is approximately −0.51~1.30%, which strongly indicates that the simplified Equation (5) has a very high level of accuracy.

Mathematical Model
The disc rotation squeezing model established in Section 2 was further simplified into a differential slab squeezing model. According to the balance equations, mathematical formulas including velocity, pressure distribution, total squeezing force, and power consumption in the wedge clearance were obtained.

Kinematic Exchange
During the actual feeding process, the rubber strip that is adhered to the screw flight rotates with the rotational screw flight and slides along the barrel surface. Therefore, in the disc rotation squeezing model, the disc drags the rubber strip to rotate through an angle θ (rad) within time t, and the rubber strip is squeezed (compressed) from the initial thickness H to the thickness h(t). This process can be regarded as one in which a "differential slab" of arc length ∆l (Figure 4) squeezes the strip in parallel from the initial thickness H to the thickness h(t) within time t at the squeezing velocity generated by the disc rotation. The squeezing flow obtained from the kinematic exchange is called the differential slab squeezing flow, which is shown in Figure 6. Since screw flight is always embedded into the rubber strip during the feeding process, this flow belongs to a constant area squeezing flow [43].
rotates with the rotational screw flight and slides along the barrel surface. Therefore, in the disc rotation squeezing model, the disc drags the rubber strip to rotate through an angle θ (rad) within time t, and the rubber strip is squeezed (compressed) from the initial thickness H to the thickness h(t). This process can be regarded as one in which a "differential slab" of arc length ∆ (Figure 4) squeezes the strip in parallel from the initial thickness H to the thickness h(t) within time t at the squeezing velocity generated by the disc rotation. The squeezing flow obtained from the kinematic exchange is called the differential slab squeezing flow, which is shown in Figure 6. Since screw flight is always embedded into the rubber strip during the feeding process, this flow belongs to a constant area squeezing flow [43].

Velocity and Pressure Distributions
The theoretical analysis used the rectangular coordinate system shown in Figure 7. First, the following assumptions about the differential slab squeezing flow were made to solve the balance equations for the analytic solutions of the velocity and pressure distributions:

Velocity and Pressure Distributions
The theoretical analysis used the rectangular coordinate system shown in Figure 7. First, the following assumptions about the differential slab squeezing flow were made to solve the balance equations for the analytic solutions of the velocity and pressure distributions: the disc rotation squeezing model, the disc drags the rubber strip to rotate through an angle θ (rad) within time t, and the rubber strip is squeezed (compressed) from the initial thickness H to the thickness h(t). This process can be regarded as one in which a "differential slab" of arc length ∆ (Figure 4) squeezes the strip in parallel from the initial thickness H to the thickness h(t) within time t at the squeezing velocity generated by the disc rotation. The squeezing flow obtained from the kinematic exchange is called the differential slab squeezing flow, which is shown in Figure 6. Since screw flight is always embedded into the rubber strip during the feeding process, this flow belongs to a constant area squeezing flow [43].

Velocity and Pressure Distributions
The theoretical analysis used the rectangular coordinate system shown in Figure 7. First, the following assumptions about the differential slab squeezing flow were made to solve the balance equations for the analytic solutions of the velocity and pressure distributions:  (i) The flow of the rubber compound is isothermal, has a quasi-steady state [42,44], and is laminar. In the quasi-steady state, locally and instantaneously, the squeezing flow may be regarded as a steady flow between fixed parallel plates [45].
(ii) The inertia, gravity, and normal stress terms in the momentum equation can be ignored.
(iii) The rubber compound is a noncompressible Newtonian liquid.
To simplify the analysis of the parallel plate squeezing flow, the component of momentum in the direction parallel to the plates (here, it is the x-component of momentum) is usually assumed to be the most important, and the component of momentum in the direction normal to the plates (here, it is the z-component of momentum) is not used [45]; that is, the component of velocity v z normal to the plates is neglected [46].
Because the minimum value h min of the feed wedge gap is very small (h min ≈ 0) and there is an additional pressure build-up created by the feed roll/screw system in the upstream position of the maximum wedge gap, the component of velocity v y in the direction of the differential dimension (y-direction) can be assumed to be zero in the differential slab squeezing flow. In other words, the rubber compound can only be squeezed out along the directions on both sides of the disc (or flight) (x-direction) during the feeding process. This is an important boundary condition for the disc rotation squeezing model or the differential slab squeezing flow.
Based on the above assumptions and analyses, and considering the isotropic pressure and the component of velocity to be constant at the differential dimension dy, the equations of momentum can be reduced to a one-dimensional form as follows: where v x is the x-component of velocity and p is isotropic pressure. Because the rubber strip slides along the barrel surface during the feeding process, it is assumed that there is no slippage at the moving slab (screw flight crest) and at the stationary slab of the x-direction, that is, the boundary condition z = 0, v x = 0 and z = h, v x = 0. Integrating Equation (7) twice with respect to z and applying the boundary conditions can derive Equation (8): This v x velocity distribution must satisfy the following overall continuity relationship: By substituting Equation (9) into Equation (8) and considering dy as a constant value: Usually, the screw channel at the feed port is not fully filled with rubber [5]. Therefore, it can be assumed that the pressure on both sides of the disc (screw flight) is atmospheric pressure, that is, x = ±w/2, p = 0. Under this boundary condition, integrating Equation (10) can obtain: Equation (11) shows that the pressure distribution p in the model wedge p is related to V, w, h, x, and η, and the effects of h, x, and w on p are particularly significant. For a given disc rotation squeezing model, when the disc rotating speed N is constant, p is only related to x at any squeezing thickness h, showing the parabolic distribution shown in Figure 8. Figure 8 shows that, for an arbitrary h, the maximum pressure p max can be obtained at x = 0. Equation (11) shows that the pressure distribution p in the model wedge p is related to V, w, h, x, and η, and the effects of h, x, and w on p are particularly significant. For a given disc rotation squeezing model, when the disc rotating speed N is constant, p is only related to x at any squeezing thickness h, showing the parabolic distribution shown in Figure 8. Figure 8 shows that, for an arbitrary h, the maximum pressure pmax can be obtained at x = 0.
According to Equation (11), it is easy to calculate the average pressure ̅ along w in Figure 8: Equation (12) shows that, for a given disc rotation squeezing model and rubber compound, when the disc rotation speed is constant, ̅ is inversely proportional to the cubic power h, proportional to the quadratic power w, and proportional to η. As h decreases, ̅ increases rapidly, and when ℎ → 0, p → ∞ . However, or ̅ is always a finite value because ℎ has a minimum design value hmin (at θ = θ0). In addition, Liu Y. C. and Yu F [3] used a disc rotation squeezing model device with D = 65 mm and found that the rubber strip would be longitudinally broken off at a position of approximately h = 1 mm before the design value of the minimum wedge gap hmin(=0.16 mm). hb and θb are set as the wedge gap value and the corresponding center angle when the rubber strip is longitudinally broken off, respectively. For the disc rotation squeezing model with D = 65 mm and hb = 1 mm, the θb = 73° can be obtained by calculating or drawing, but the regularity of hb and θb needs further experimental research. 6 ( ) From Equation (13), vx depends on V, h, x, and z. When V, h, and x are fixed, the relationship between vx and z shows a parabolic curve, and shear flow occurs along the zdirection. When V, h, and z are fixed, vx is proportional to x and produces elongational flow in the x direction. When V, x, and z are fixed, vx is inversely proportional to the cubic power of h; that is, when h decreases, vx increases rapidly. According to Equation (11), it is easy to calculate the average pressure p along w in Figure 8: Equation (12) shows that, for a given disc rotation squeezing model and rubber compound, when the disc rotation speed is constant, p is inversely proportional to the cubic power h, proportional to the quadratic power w, and proportional to η. As h decreases, p increases rapidly, and when h → 0 , p → ∞ . However, p or p is always a finite value because h has a minimum design value h min (at θ = θ 0 ). In addition, Liu Y. C. and Yu F [3] used a disc rotation squeezing model device with D = 65 mm and found that the rubber strip would be longitudinally broken off at a position of approximately h = 1 mm before the design value of the minimum wedge gap h min (=0.16 mm). h b and θ b are set as the wedge gap value and the corresponding center angle when the rubber strip is longitudinally broken off, respectively. For the disc rotation squeezing model with D = 65 mm and h b = 1 mm, the θ b = 73 • can be obtained by calculating or drawing, but the regularity of h b and θ b needs further experimental research.
A combination of Equations (8) and (10) obtains: From Equation (13), v x depends on V, h, x, and z. When V, h, and x are fixed, the relationship between v x and z shows a parabolic curve, and shear flow occurs along the z-direction. When V, h, and z are fixed, v x is proportional to x and produces elongational flow in the x direction. When V, x, and z are fixed, v x is inversely proportional to the cubic power of h; that is, when h decreases, v x increases rapidly.

Total Squeezing Force
The total squeezing force refers to the resultant force of the pressures exerted on the rubber strip in the radial direction of the disc (screw) within the whole feed wedge clearance. In the design of a rubber extruder, for example, the total squeezing pressure is used to calculate the elastic flex deformation of the screw caused by the side force, because excessive flex deformation will cause scraping between the screw and the barrel.
The total squeezing force F can be deduced from the disc rotation squeezing model. The method for calculating F is as follows: first, the coordinate system y Oz shown in Figure 9 was established, and the force dF acting on the area element wdl(= wRdθ) at position θ in the radial direction of the disc was decomposed into dF y and dF z components, as shown in the enlarged view. Second, the magnitudes of the F y and F z components of the total squeezing force F were calculated by integral, that is, F y = ∑ dF y and F z = ∑ |dF z |. Finally, the magnitude and tangential direction angle of the total squeezing force F were calculated, that is, F = F y 2 + F z 2 and tan ψ = F z /F y . ϕ + arctanψ is the angle between the direction of the total squeezing force and the horizontal line.
As mentioned in Equation (7), the pressure in the model wedge is isotropic. Therefore, the pressure in the radial direction in Figure 9 is equal to the pressures calculated by Equations (11) and (12). To simplify the calculation, the average pressure p was used. Hence: Putting Equation (12) into Equations (14) and (15) and integrating them obtains: where ponents, as shown in the enlarged view. Second, the magnitudes of the and components of the total squeezing force were calculated by integral, that is, = ∑ and = ∑| |. Finally, the magnitude and tangential direction angle of the total squeezing force were calculated, that is, = + and tan = ∕ . + arctan is the angle between the direction of the total squeezing force and the horizontal line. As mentioned in Equation (7), the pressure in the model wedge is isotropic. Therefore, the pressure in the radial direction in Figure 9 is equal to the pressures calculated by Equations (11) and (12). To simplify the calculation, the average pressure ̅ was used. Hence: Putting Equation (12) into Equations (14) and (15) and integrating them obtains: Substituting Equations (16) and (17) into the above expressions of F and tan , respectively, can obtain: Substituting Equations (16) and (17) into the above expressions of F and tan ψ, respectively, can obtain: For each given model wedge, ϕ is known; therefore, according to Equation (19), the angle ϕ + arctanψ between the direction of the total squeezing force and the horizontal line can be obtained.

Power Consumption
When the rubber strip passes through the model wedge, the power consumption comes from the following three parts: (i) The increase in kinetic energy when the rubber compound is squeezed out from both sides of the disc (screw flight); (ii) Viscous dissipation in squeezing flow; (iii) Friction loss of the rubber strip sliding along the barrel.
(1) Power consumption P 1 caused by an increase in kinetic energy As h decreases, v x shows a rapid increasing trend (Equation (13)), which increases the kinetic energy of the rubber in the extrusion direction (x direction); the kinetic energy reaches its maximum at the edges at both sides of the disc (screw flight). In Figure 7, the volumetric flow rate in the direction of +x or −x caused by the moving slab squeezing downwards is Vxdy (V is the average squeezing velocity, which is shown in Equation (6)). Therefore, the average extrusion velocity v x at the x position is: If the flow between the approaching slabs is regarded as potential flow, the average velocity v x in the x direction is a function of x, not h [16]. The average velocity of the rubber extruded from the edges of the slabs (x = ±w/2) is v w/2 = Vw/2h (only the magnitude of the velocity considered). Therefore, the rate of work performed on a moving slab with dy length elements is equal to the rate of change of the kinetic energy of rubber dP 1 . where dm/dt is the mass flow rate of the extruded rubber, which is equal to 2hdy(Vw/2h)ρ, and ρ is the density of rubber, so Equation (21) becomes: According to Figure 2, considering dy = Rdθ and integrating Equation (22) along the whole wedge gap yields P 1 : Considering the expression of h(θ) (Equation (3)) will make the integral operation of Equation (23) very complex. Therefore, the simple arithmetic mean H/2 of h(θ) was used to replace h(θ) to obtain the simple approximate solution of P 1 : (2) Power consumption P 2 caused by viscous dissipation. The velocity gradient of v x in the z direction (Equation (13)) will cause viscous heat generation. Therefore, the shear rate . γ can be obtained by differentiating Equation (13) with respect to z: From Equation (25) γ max = 3Vw/h 2 . To simplify the calculation, the corresponding shear rate at h = H/2 can be taken as the average shear rate . γ in the whole squeezing flow process: that is, . γ = 12Vw/H 2 . Therefore, the rate of viscous dissipation of rubber compound per unit volume is: If the volume of the whole model wedge is V w and the very small volume remaining after the longitudinal break-off of rubber strip is ignored, the power consumption P 2 caused by viscous dissipation is: Through the geometric relationship in Figure 3, V w can be obtained: (3) Power consumption P 3 caused by sliding friction. During the sliding process of rubber strips dragged by the disc (screw flight) along the barrel, the shape of the rubber strip is bent, and the local pressure is continuously changed. As shown in Figure 10, for a rubber element at position θ, setting ∠OnO 1 = β, according to the cosine theorem: (3) Power consumption caused by sliding friction. During the sliding process of rubber strips dragged by the disc (screw flight) along the barrel, the shape of the rubber strip is bent, and the local pressure is continuously changed. As shown in Figure 10, for a rubber element at position , setting ∠ = , according to the cosine theorem: Figure 10. Velocity vector analysis of the rubber strip during squeezing; (a) enlarged drawing of velocity vector analysis. Equation (29) shows that, in the feeding process, the value of β will change with h or θ. Since the center and radius of the cylinder arcê s are O 1 and R 1 , respectively, the tangential velocity v bt of the rubber element sliding along the cylinder at position θ is: Obviously, v bt changes with h or θ.
Assuming that the total squeezing force F (Equation (18)) is evenly distributed on arc e s , and the sliding friction coefficient f between the rubber strip and the barrel is a constant, the friction force dF f generated by the barrel on the rubber element at position θ is: where α = ∠eO 1 s for a given model wedge gap and α is a known constant. From Equations (30) and (31), the friction power dP 3 generated by the barrel to the rubber element at position θ can be obtained: The power consumption P 3 caused by sliding friction in the whole feeding wedge is: To simplify the integral operation of Equation (33), the average value β of β is used, that is, β = ∠OeO 1 /2 . For a given model wedge gap, ∠OeO 1 is known, and β and cos β are known constants. For example, for the disk rotating squeezing model with D = 65 mm, ∠OeO 1 ≈ 9.5 • , β = 4.75 • , and cos β ≈ 0.9966. So Equation (34) becomes where By combining Equations (24), (27) and (34), the power consumption P in the process of the rubber strip passing through the model wedge gap is obtained: That is: It should be noted that the above mathematical expressions of total squeezing force and power consumption only consider a single disc. However, for an actual cold feed extruder, there are usually at least double-threaded flights at the feed port, and the calculation of the total squeezing pressure and the power consumption needs to multiply the above expressions by the number of flights.
The theoretical analysis model of the feeding process of cold-feed rubber extruders has important guiding significance for analyzing feeding problems and optimizing and innovating feeding structures: (1) the side force applied to the screw and the feeding power consumption can be roughly calculated without a feed roll, which provides basic calculation parameters for rubber extruder design and overcomes the pure experience of feed structure design. (2) The width of the screw flight at the feed screw segment should be as narrow as possible, because the total squeezing force and power consumption are proportional to the third power of the width of the screw flight. However, a screw width that is too narrow can reduce the drag-in force of the screw flight crest so that the rubber strip cannot be dragged through the feed wedge clearance. For a screw flight with a narrower width, shallow grooves can be opened on the top or the side of the screw flight to increase friction drag action, or else the screw flight can be interrupted to increase "penetration" drag action. The latter is often used in the feed screw segment of modern cold-feed extruders. (3) Under the condition of the same feed wedge clearance length, the maximum value (H) in the feed wedge clearance should be as large as possible, because the total squeezing force and the power consumption are proportional to the reciprocal of the power of H. In other words, a thicker rubber strip should be used for feeding. Therefore, the improvement direction of the structure of the feed wedge clearance was to increase the eccentricity values a and b and select the appropriate eccentric arc radius R. This design can shorten the feeding time and improve the feeding efficiency. (4) Because the power consumption mainly depends on the frictional resistance of the barrel to the rubber strip (see Section 5), the inner wall of the barrel should have a lower friction coefficient at the feed port; that is, the inner wall of the barrel should be smoother than it has been in previous designs. (5) If the feed wedge clearance, feed screw flight, and feed barrel are fully optimized, the effect of feed roll can be reduced or even rendered unnecessary, and the feeding power consumption can be further decreased.

Experiment and Materials
The pressure and power data obtained by using the disc rotation squeezing model device ( Figure 11) were chosen. Test device: this model device comprises disc 2 driven by motor 1, semicylinder 4 with an eccentric undercut groove on the inner surface, pressure sensor 5, and a control system (which is not shown in the figure). The top surface of the disc, and the inner wall of the groove under the semicylindrical feed wedge clearance 6 and 7, using feeding rubber strips. The disc and drive device are supported by bracket 3. The DJYZ-10 cylindrical pressure sensors were installed in the upper (U), middle (M), and lower (L) positions of the semicylinder, and were used to measure the pressures of the rubber strip through the feed wedge clearance. The pressure sensor had a measurement range of 0-500 kg and a nonlinearity of 0.3-0.5%. The pressure and power data were automatically calculated, recorded, and stored by upper computer software. The disc speed was changed by a variable frequency speed-regulating motor. The temperature of the laboratory was set at 23 ± 2 • C. The geometric parameters of the model device are shown in Table 1.  Test materials: truck radial tire compound, ruck tire inner liner compound (TTI) ((51 ML1 + 4 (100 °C); the Mooney viscosity of the rubber obtained by testing at 100 ° C is 51), truck tire sidewall compound (TTS) (55 ML1 + 4 (100 °C)), and truck tire tread compound (TTT) (65 ML1 + 4 (100 °C)), provided by Shandong Anchi Tire Co., Ltd. The size of the feeding rubber strip was approximately 5.4 (thick) × 65 (width) × 100 (length) mm.
Test method: rotation squeezing tests were carried out at room temperature using the above three kinds of rubber compounds to measure the pressure distribution and power consumption of the feed wedge clearance under different disc speeds (30,45, and 60 rpm).

Measurements of Viscosity and the Sliding Friction Coefficient
The shear viscosity and sliding friction coefficient of the above three compounds at room temperature are two physical parameters required by the theoretical model. These two parameters must be measured.
(1) Viscosity Due to the high viscosity of the rubber compound, it was necessary to use a parallel plate plastometer to measure the shear viscosity. The model of parallel plate plastometer used is the MZ-4014. The measurement was made at 23 ± 2 °C. The measurement principle is that, under a given temperature, for a squeezing flow with constant volume, the shear viscosity can be calculated with the following formula using applied load , load time t, measured sample heights of ℎ and ℎ before and after the load , and measured sample volume [47]: The sample was a cylindrical sample with a diameter of 16 mm and a thickness of 3 mm. The applied load was 49 N, and the applied load time was 60 s. For the above three rubber compounds, the shear viscosity values obtained by calculation are shown in Table  2.   Test materials: truck radial tire compound, ruck tire inner liner compound (TTI) ((51 ML1 + 4 (100 • C); the Mooney viscosity of the rubber obtained by testing at 100 • C is 51), truck tire sidewall compound (TTS) (55 ML1 + 4 (100 • C)), and truck tire tread compound (TTT) (65 ML1 + 4 (100 • C)), provided by Shandong Anchi Tire Co., Ltd. The size of the feeding rubber strip was approximately 5.4 (thick) × 65 (width) × 100 (length) mm.
Test method: rotation squeezing tests were carried out at room temperature using the above three kinds of rubber compounds to measure the pressure distribution and power consumption of the feed wedge clearance under different disc speeds (30,45, and 60 rpm).

Measurements of Viscosity and the Sliding Friction Coefficient
The shear viscosity and sliding friction coefficient of the above three compounds at room temperature are two physical parameters required by the theoretical model. These two parameters must be measured.
(1) Viscosity Due to the high viscosity of the rubber compound, it was necessary to use a parallel plate plastometer to measure the shear viscosity. The model of parallel plate plastometer used is the MZ-4014. The measurement was made at 23 ± 2 • C. The measurement principle is that, under a given temperature, for a squeezing flow with constant volume, the shear viscosity η can be calculated with the following formula using applied load F N , load time t, measured sample heights of h 0 and h before and after the load F N , and measured sample volume V r [47]: The sample was a cylindrical sample with a diameter of 16 mm and a thickness of 3 mm. The applied load was 49 N, and the applied load time was 60 s. For the above three rubber compounds, the shear viscosity values obtained by calculation are shown in Table 2. Since the temperature changes from 23 • C to 50 • C during the experiment in previous studies, the viscosity corresponding to the temperature change can calculate according to the theoretical formula.
(2) Sliding friction coefficient The sliding friction coefficient of the rubber compounds at room temperature was measured by Anton Paar TRB 3 , and the measurement was made at 23 ± 2 • C. To obtain the approximate value of the sliding friction coefficient between the rubber strip and the inner wall of the barrel, the measured linear speed was made close to the linear speed of the disc rotation squeezing model, and the unit area pressure applied to the test metal element close to the average pressure in the disc rotation squeezing model. The test metal element material is Q235-A. The values of the sliding friction coefficient obtained by measurement are shown in Table 3. Table 3. Sliding friction coefficients of different rubber compounds at room temperature.

Results and Discussion
(1) Pressure distribution In Section 2.2, according to Equation (6), the average squeezing velocity V of the model device in Table 1 at N = 30, 45, and 60 rpm was −11.4, −17.1, and −22.8 mm·s −1 , respectively. In Figure 11, the h 1 (= H), h 2 , and h 3 of the design wedge gap at the different positions of the pressure sensor installed in the upper (U), middle (M), and lower (L) sections were 5.4, 4.7, and 2.6 mm, respectively.
The w value in Table 1, the η value of different rubber compounds in Table 2, the V value at different disc speeds, and the h value at different positions were substituted into Equation (12) ( p = ηVw 2 /h 3 ) to obtain the pressure distributions of different rubber compounds at different disc speeds. Figure 12 shows the pressure distribution of TTI, TTS, and TTT at N = 30 rpm. Figure 13 shows the pressure changes of TTI, TTS, and TTT in the position of h 2 = 4.7 mm at different disc speeds. Figure 12 shows that, under a given disc speed, for different kinds of rubber strips, both theoretical prediction and experimental data had approximate pressure distribution curves, and the increase in pressure with the decrease in h in the region of large h was much slower than that in the region of small h, because p was proportional to the reciprocal of the third power of h (Equation (12)). However, except for the start position of squeezing h 1 (= H), the difference between the theoretical predicted pressure and the experimental value is 0.5-2 Mpa. This is the result of two main factors: average squeezing velocity V and viscosity η. In the initial stage of squeezing, V is greater than the instantaneous squeezing velocity . h ( Figure 5), and the theoretical pressure obtained using Equation (12) is higher than the test data. In the middle and end stages of squeezing, V is less than . h ( Figure 5), and there are also obvious shear thinning and viscous heat generation effects. Although the pressure in the last two stages was calculated by using V, which was smaller than . h, the original viscosity data were larger, and the viscosity decrease caused by shear thinning and viscous heat generation was not considered; therefore, the theoretical prediction pressure was higher than that attained in the experimental data.
(1) Pressure distribution In Section 2.2, according to Equation (6), the average squeezing velocity of the model device in Table 1 at N = 30, 45, and 60 rpm was −11.4, −17.1, and −22.8 mm·s −1 , respectively. In Figure 11, the ℎ (= ), ℎ , and ℎ of the design wedge gap at the different positions of the pressure sensor installed in the upper (U), middle (M), and lower (L) sections were 5.4, 4.7, and 2.6 mm, respectively.
The value in Table 1, the value of different rubber compounds in Table 2, the value at different disc speeds, and the ℎ value at different positions were substituted into Equation (12) ( ̅ = ∕ ℎ ) to obtain the pressure distributions of different rubber compounds at different disc speeds. Figure 12 shows the pressure distribution of TTI, TTS, and TTT at N = 30 rpm. Figure 13 shows the pressure changes of TTI, TTS, and TTT in the position of ℎ = 4.7mm at different disc speeds.   Figure 12 shows that, under a given disc speed, for different kinds of rubber strips, both theoretical prediction and experimental data had approximate pressure distribution curves, and the increase in pressure with the decrease in ℎ in the region of large ℎ was much slower than that in the region of small ℎ, because ̅ was proportional to the reciprocal of the third power of ℎ (Equation (12)). However, except for the start position of squeezing ℎ (= ), the difference between the theoretical predicted pressure and the experimental value is 0.5-2 Mpa. This is the result of two main factors: average squeezing velocity and viscosity . In the initial stage of squeezing, is greater than the instantaneous squeezing velocity ℎ (Figure 5), and the theoretical pressure obtained using Equation (12) is higher than the test data. In the middle and end stages of squeezing, is less than ℎ ( Figure 5), and there are also obvious shear thinning and viscous heat generation effects. Although the pressure in the last two stages was calculated by using , which was smaller than ℎ , the original viscosity data were larger, and the viscosity decrease caused by shear thinning and viscous heat generation was not considered; therefore, the theoretical prediction pressure was higher than that attained in the experimental data. Figure 13 shows that both pressures obtained by theoretical calculation and experimental measurement increase almost linearly with increasing disc speed; the reason for this phenomenon is that ̅ is proportional to the disc speed (Equations (6) and (13)). Additionally, the difference between theoretical pressure and the experimental value is about 0.05-0.52 Mpa. However, Figure 13 also shows that the theoretical value was greater than the experimental data, and the difference increased significantly with increasing . This is because the position of ℎ = 4.7 mm occurs in the middle stage of squeezing. With the increase in disc speed , the effects of shear thinning and viscous heat generation caused by the increase in shear rate (Equation (27)) increase. It has been reported in the literature [3] that the temperature of the squeezed TTI rubber compound is approximately 27 °C higher than room temperature when the disc rotation squeezing model test for TTI Figure 13. Comparison of calculated and experimental pressure changes of h2 = 4.7 mm for TTI, TTS, and TTT at different disc speeds. Figure 13 shows that both pressures obtained by theoretical calculation and experimental measurement increase almost linearly with increasing disc speed; the reason for this phenomenon is that p is proportional to the disc speed (Equations (6) and (13)). Additionally, the difference between theoretical pressure and the experimental value is about 0.05-0.52 Mpa. However, Figure 13 also shows that the theoretical value was greater than the experimental data, and the difference increased significantly with increasing N. This is because the position of h 2 = 4.7 mm occurs in the middle stage of squeezing. With the increase in disc speed N, the effects of shear thinning and viscous heat generation caused by the increase in shear rate (Equation (27)) increase. It has been reported in the literature [3] that the temperature of the squeezed TTI rubber compound is approximately 27 • C higher than room temperature when the disc rotation squeezing model test for TTI strips was carried out under N = 60 rpm without cooling the disc and the semicylinder.
(2) Power consumption The density of the rubber compound is approximately 1.5 × 103 kg·m −3 . The relevant parameter values were substituted into Equation (37) to obtain the power consumptions of P1, P2, and P3 for the TTI, TTS, and TTT at different disc speeds. The range of P1 is 1.8 × 10 −6~1 .5 × 10 −5 W, P2 is 2.5~13 W, and P3 is 200~900 W. Compared with P3, P1 and P2 were small and negligible. In other words, during the feeding process, the power consumption of the strip through the feed wedge clearance mainly comes from the sliding friction of the strip along the barrel. This conclusion is of great significance for the feed structure design of a cold-feed extruder.  Figure 14 shows that the power consumption obtained by theoretical calculation and by experimental measurement increased linearly with increasing disc speed, and the former increased faster than the latter. The difference between the theoretical predicted pressure and the experimental value is about 40-200 W. Moreover, from the data comparison, it was found that the theoretical value is almost greater than the experimental value. This is mainly because the large viscosity data were used for the calculation of the total squeezing force without considering the significant decrease in viscosity caused by shear thinning and viscous heat generation, especially in the middle and end stages of the feeding process at a high disc speeds. In other words, during the test, with the increase in disc speed, the decrease in rubber viscosity would make the squeezing force small, eventually resulting in a decrease in power consumption. However, in the theoretical calculation, a constant and high value of rubber viscosity was used, so the theoretical value of the power consumption was greater than its experimental value. However, in the actual feeding process of the cold feed extruder, due to the cooling of the screw and the feed barrel, the temperature rise in the rubber compound in the feed wedge clearance would be greatly limited. The effect of the decrease in viscosity caused by the temperature rise would be smaller, and the viscosity decrease mainly arises due to the contribution of shear thinning. This would increase the measured values of the pressure and power consumption, and lessen the difference between the theoretical prediction and the experimental data.
Therefore, for the feeding process of a cold feeding rubber extruder, the piecewise change in the squeezing velocity, shear thinning, viscous heat generation, and heat transfer must be considered if we are to develop a more accurate mathematical model. Since the Deborah number of the rubber compound passing through the feed wedge clearance at a high screw speed is probably large, feeding analysis theory should also consider the elastic response of the rubber compound [23,44,48].

Conclusions
This paper offered a theoretical analysis of the feeding process of rubber strips passing through a feed wedge clearance in a cold-feed extruder. First, by simplifying the screw flight rotation squeezing process into a disc rotation squeezing process, the instantaneous squeezing velocity ℎ ( ) in the disc rotation squeezing model was derived according to the feed wedge clearance geometry and the disc rotating speed. The ℎ − curve comprised two general parts: a constant acceleration squeeze in the early stage and a constant velocity squeeze in the later stage. To simplify the theoretical analysis, the squeezing process of the whole feed wedge clearance was regarded as constant velocity squeezing, and the average squeezing velocity was used to replace the instantaneous squeezing velocity ℎ . Second, by transforming rotation squeezing into differential slab squeezing, mathematical expressions of the velocity distribution, pressure distribution, total squeezing force, and power consumption in the process of the rubber strip passing through the model wedge were derived; this analysis utilized the rectangular coordinate system under  Figure 14 shows that the power consumption obtained by theoretical calculation and by experimental measurement increased linearly with increasing disc speed, and the former increased faster than the latter. The difference between the theoretical predicted pressure and the experimental value is about 40-200 W. Moreover, from the data comparison, it was found that the theoretical value is almost greater than the experimental value. This is mainly because the large viscosity data were used for the calculation of the total squeezing force without considering the significant decrease in viscosity caused by shear thinning and viscous heat generation, especially in the middle and end stages of the feeding process at a high disc speeds. In other words, during the test, with the increase in disc speed, the decrease in rubber viscosity would make the squeezing force small, eventually resulting in a decrease in power consumption. However, in the theoretical calculation, a constant and high value of rubber viscosity was used, so the theoretical value of the power consumption was greater than its experimental value. However, in the actual feeding process of the cold feed extruder, due to the cooling of the screw and the feed barrel, the temperature rise in the rubber compound in the feed wedge clearance would be greatly limited. The effect of the decrease in viscosity caused by the temperature rise would be smaller, and the viscosity decrease mainly arises due to the contribution of shear thinning. This would increase the measured values of the pressure and power consumption, and lessen the difference between the theoretical prediction and the experimental data.
Therefore, for the feeding process of a cold feeding rubber extruder, the piecewise change in the squeezing velocity, shear thinning, viscous heat generation, and heat transfer must be considered if we are to develop a more accurate mathematical model. Since the Deborah number of the rubber compound passing through the feed wedge clearance at a high screw speed is probably large, feeding analysis theory should also consider the elastic response of the rubber compound [23,44,48].

Conclusions
This paper offered a theoretical analysis of the feeding process of rubber strips passing through a feed wedge clearance in a cold-feed extruder. First, by simplifying the screw flight rotation squeezing process into a disc rotation squeezing process, the instantaneous squeezing velocity . h(t) in the disc rotation squeezing model was derived according to the feed wedge clearance geometry and the disc rotating speed. The . h − t curve comprised two general parts: a constant acceleration squeeze in the early stage and a constant velocity squeeze in the later stage. To simplify the theoretical analysis, the squeezing process of the whole feed wedge clearance was regarded as constant velocity squeezing, and the average squeezing velocity V was used to replace the instantaneous squeezing velocity . h. Second, by transforming rotation squeezing into differential slab squeezing, mathematical expressions of the velocity distribution, pressure distribution, total squeezing force, and power consumption in the process of the rubber strip passing through the model wedge were derived; this analysis utilized the rectangular coordinate system under isothermal and quasi-steady assumptions and certain boundary conditions by using balance equations and the Newtonian viscous constitutive relation. Third, the shear viscosity and sliding friction coefficient of three kinds of rubber compounds were measured at room temperature and under conditions close to the actual feeding process by using a parallel plate plastometer and a CSM Anton Paar TRB 3 . The measured shear viscosity and sliding friction coefficient were substituted into the theoretical formulas to compare the theoretical prediction with previous experimental data. The comparison of pressure distribution and power consumption showed that the established theoretical model of the feeding process of cold-feed rubber extruders was not very accurate, on a quantitative level; however, it might reflect the main factors that affect the feeding process, and can reveal the main trends in feeding behavior. It was also noted that the piecewise change of the squeezing velocity, shear thinning, viscous heat generation, and heat transfer must be considered in a more accurate mathematical model, and the elastic response of the rubber compound should also be considered at a high screw speeds.