Research Article  Open Access
Rudolf Pernis, Tibor Kvackaj, "New Numerical Solution of von Karman Equation of Lengthwise Rolling", Journal of Applied Mathematics, vol. 2015, Article ID 843720, 20 pages, 2015. https://doi.org/10.1155/2015/843720
New Numerical Solution of von Karman Equation of Lengthwise Rolling
Abstract
The calculation of average material contact pressure to rolls base on mathematical theory of rolling process given by Karman equation was solved by many authors. The solutions reported by authors are used simplifications for solution of Karman equation. The simplifications are based on two cases for approximation of the circular arch: (a) by polygonal curve and (b) by parabola. The contribution of the present paper for solution of twodimensional differential equation of rolling is based on description of the circular arch by equation of a circle. The new term relative stress as nondimensional variable was defined. The result from derived mathematical models can be calculated following variables: normal contact stress distribution, front and back tensions, angle of neutral point, coefficient of the arm of rolling force, rolling force, and rolling torque during rolling process. Laboratory cold rolled experiment of CuZn30 brass material was performed. Work hardening during brass processing was calculated. Comparison of theoretical values of normal contact stress with values of normal contact stress obtained from cold rolling experiment was performed. The calculations were not concluded with roll flattening.
1. Introduction
It is likely that the author Karman [1] in the year 1925 was the first to submit equation of power balance in the roll gap and derived the twodimensional differential equation which described lengthwise rolling process. The following simplifications were received:(i)Rolled material has rectangular cross section with initial thickness and is deformed by cylindrical rolls to final thickness without any action of front and back tensions.(ii)Spread is negligible.(iii)Friction coefficient between rolls and rolled material is constant.(iv)Rolls are rigid without any elastic deformation and bending.(v)Rolled material is homogenous and without any elastic deformation.(vi)Tresca condition of plasticity: .(vii)Flow stress of material in rolling gap is constant.(viii)Rolling speed is constant.
Geometrical description of roll gap is given in Figure 1 and essential formulas are presented in the appendix. From graphical scheme of roll gap the result is that all forces have to be at balance which was the base for derivation twodimensional differential equation describing lengthwise rolling process. The procedure of derivation of the twodimensional differential equation can be found in the classical literature of rolling, for example, Avitzur [2], Hensel and Spittel [3], and Mielnik [4]. More recent literatures are the publications Hajduk and Konvičný [5] and Kollerová et al. [6]. The stress state in rolling gap describes differential equation:where plus sign (+) is for forward slip zone, minus sign (−) is for backward slip zone, is normal contact stress on rolls, is the shear stress between rolls and rolling material, and , are the coordinates of the cylinder touching the rolled material.
The Tresca condition of plasticity is used for solution of (1):where is the maximal principal stress (vertical direction), is the minimal principal stress (horizontal direction), is the flow stress [7–9], and is the shear stress that is proportional to normal contact stress according to formula:where is the friction coefficient between rolls and rolled material.
Substituting (2) into (3) the following formula will be obtained with only one stress variable :Flow stress and friction coefficient are as a constant in (4). If differential equation (4) is divided by flow stress it is possible to obtain relative normal contact stress :where is a function.
Substituting (5) into (4) the following form is obtained:The application equation (5) can be excluded from (4) material constant (flow stress ). Thus (6) become independent from the rolled material. This differential equation is a linear differential equation of first order with the right side and has a general transcription: . The coordinates and are defining the geometry of the cylinder which is connected to the rolling material. The variable is a function of the coordinate. If equation of a circle is substituting into (6), so differential equation is obtained which has no analytical solution till now. The author Orowan [10] presented the solution of differential equation (1) based on analyticalgraphical principle. He defined procedure for calculation rolling force, rolling torque, and electric input for hot and cold rolling processes. In the literature [11] showed solution of differential equation (1) by simulation on analog computer. The authors of [12] solved this case with the use of the circular arc but the result is purely numerical solution. The author of [13] presented simplified analytical solution with acceptance of material work hardening during cold rolling process. The development in computational software allowed to solve this task by using finite element method as was given in the literature [14].
2. Experimental Material and Methods
The investigated material was annealed brass CuZn30 with chemical composition corresponding to standard DIN 17660, W.No. 2.0265. The sample sizes before rolling were = 3,4 × 30 × 150 mm. The samples were laboratory rolled using duo mill (the diameter of rolls = 210 mm) with constant circumferential speed of rolls = 0,66 m/s at ambient temperatures. The temperature of a sample surface immediately after rolling was measured using a thermovision camera and was from interval [°C]. The cold rolling reductions of thickness were made in an interval as one pass deformation. Deformation forces during samples processing by cold rolling were measured by tensometric elements with signal registration in Spider apparatus. The static tensile test was performed in accordance with STN EN 10 0021 standard.
The differential equation of lengthwise rolling process given by formula (1) has no analytical solution when contact arc is described with the equation of the circle. For solution of this problem a new combined analytical and numerical approach is presented.
3. Analysis and Results
3.1. Analysis of the Present Solution of Differential Equation Describing Stress Condition in Rolling Gap
The present solutions of differential equation of rolling were based on the same simplifications because the form of equation is complicated. More than twenty solutions of differential equation are described in the literature. The present solutions of differential equation (1) can be divided into following groups:(i)Geometry of cylinder: the circular arc of the cylinder is being approximated by simplifying curves as straight line, polygonal line, and parabola or in some cases using equation of a circle.(ii)Friction conditions between rolls and rolled material:(a)The shear stress between rolls and rolling material is constant: .(b)The shear stress between rolls and rolling material is changing: .(c)Friction coefficient is .(iii)Mathematical methods of solution: analytical, analyticalgraphical, analyticalnumerical, numerical, and semiempirical.(iv)Work hardening of rolled material: most of the solution is based on assumption that material during processing is without work hardening. In a few isolated cases solution with work hardening of rolled material was used.
In the following part most frequently ways for calculation of function which are fall to first group will be analyzed.
3.1.1. The Solution according to Korolev
The author Korolev [15] presented solution based on approximation of circular arch by parallel step function which is shown in Figure 2. The circular arch (BA) is approximated by polygonal line. Neutral point N divides space of rolling gap on parts: backward slip zone (BN) and forward slip zone (NA). Position of neutral point is determined by coordinate : The distribution functions describing relative normal contact stress in rolling gap were divided separately for backward and forward slip zone.
The distribution function describing relative normal contact stress in backward slip zone of rolling gap is as follows:whereThe distribution function describing relative normal contact stress in forward slip zone of rolling gap is as follows: whereThe author of [15] on the basis of the previous distribution function derived an equation for calculation function as follows:where is constant:
3.1.2. The Solution according to Tselikov
The authors Tselikov et al. [16] presented solution based on approximation of circular arch by straight line function which is shown in Figure 3. The circular arch (BA) is approximated by straight line. Neutral point N divides space of rolling gap on parts: backward slip zone (BN) and forward slip zone (NA). Position of neutral point is determined by coordinate and thickness in neutral section was derived as follows: where is constant:The distribution functions describing relative normal contact stress in forward slip zone and in backward slip zone of rolling gap are as follows: The authors of [16] on the basis of the previous distribution functions derived an equation for calculation function as follows:From (17) the result is that, for condition , function is discontinuous.
3.1.3. The Solution according to Bland and Ford
The authors Bland and Ford [17] presented solution based on approximation of circular arch by parabolic function which is shown in Figure 4. The circular arch (BA) is approximated by part of parabola. Neutral point N divides space of rolling gap on parts: backward slip zone (BN) and forward slip zone (NA). The authors defined neutral angle according to the following formula:The distribution functions describing relative normal contact stress in forward slip zone and in backward slip zone of rolling gap are as follows: whereThe authors of [17] on the basis of the previous distribution functions derived an equation for calculation function as follows:where is constant and , are limits of integration, :
3.1.4. The Solution according to Sims
The author Sims [18] presented solution of differential equation without consideration of friction coefficient. Roll radius is considered with elastic flattening . The ratio was derived from the following formula: whereThe distribution functions describing relative normal contact stress in forward slip zone and in backward slip zone of rolling gap were calculated according to the following formulae: The author of [18] on the basis of the previous distribution functions derived an equation for calculation function as follows:
3.1.5. The Modification Smiryagin Solution by Pernis
The authors of [19] for solution of (1) presented modification of Tselikov solution using a parabolic approximation function of the circular arc as is shown in Figure 4. The distribution functions describing relative normal contact stress in forward slip zone and in backward slip zone of rolling gap were calculated according to the following formulae: where , are constants:At this point the Tselikov solution of differential equation (1) was completed. The function for calculation of neutral angle and average relative normal contact stress has not been determined. Additional calculations were carried out by the author of [20] which are derived by additional equations describing neutral angle and average relative normal contact stress. Location of the neutral point was determined through using the value (dimensionless number) and is expressed by an explicit equation:The neutral point can be obtained from the following formula: The authors of [20] on the basis of the previous distribution functions derived an equation for calculation function as follows:
3.2. New Approach for Solution of von Karman Differential Equation
The following approach represents analytic and numerical solution of (1) where for description of the circular arc equation of circle is applied as is shown in Figure 5. Equation (6) can be classified as the firstorder linear differential equation with right side:Many authors [21–23] have dealt with the solution of this type of differential equation. Firstly (35) will be solved as equation without right side:which can be divided into the two following differential equations: The following solution of (37) is valid for backward slip zone: The calculation of integral of the left side of the equation can be carried out as follows:If the integral is the exponent and also is a known function so it is possible to writeEquation (40) is representing general form of (37) and is function of variable :Differentiating the equation is obtained fromEquations (42) and (43) are substituting into (35):from which function can be determined as follows: where is an integration constant that must be determined from the boundary conditions for each segment separately (backward and forward slip zone).
Substituting (45) into (42) will obtain the general form of (35) for backward slip zone:Determination of ratio and function can be done using the equation of the circle:Separating the variable from (47) and its differentiation is obtained:The next solution is apparent from the description: where the constant values are included in one constant : The simplification of (50) can be obtained with transformation from Cartesian coordinates to polar coordinates with use of the following substitution with respect to Figure 6:Substituting (52) into (50) will obtain the formwhere a new function is expressed as follows:The function can be expressed from (41) according to the solution defined in work [24] as follows:Transformation of function from Cartesian coordinates to function in polar coordinates will be obtained: whereusing the following substitutions:which are inserted into (57) and the integral is obtained: whereBy decomposition of the integral of (59) to partial fractions the following formula is obtained:and calculation of subintegrals in (61) enters to the form Substituting (60) into (62) will obtain the formand also substituting (58) into (63) integral will be in the form Integral represents the analytical form of the function :Transformation Cartesian coordinate to polar coordinate in (46) and substituting function are obtained from the form Limits of integration in (66) are . Solution of (38) describing the forward slip process can be obtained by the analogical way as reported for solving (37). The final solution of (38) of relative normal contact stress in forward slip zone will be as follows:Limits of integration are . From (66) and (67) the result is that the relative normal contact stress is proportional to the product of two functions: .
Graphical visualization of function depending on angle coordinate and ratio is given in Figure 7. Functions throughout the observed interval are monotonically increasing. With increasing of ratio also functions are rising. The function during rolling process represents wrapping angle of roll by rolled material ( is not gripping angle). Graphical visualization of function depending on angle coordinate is given in Figure 8. Local maximum of functions increases with increasing ratio . Location point expressing the maximum value of the function is determined from the condition :where is coordinate in maximum and is constant of differential equation given by (51).
Maximal value of function can be determined from (54) for coordinate as follows: The broken line in Figure 8 shows the positions of the maximal values. The integration constants in (66) and (67) must be determined from Tresca condition of plasticity. Stress state according to Figure 1 can be described for point A (output of material from rolling gap) and for point B (input of material to rolling gap) using the following formulae:Equations (70) are normalized according to flow stress:The determination of integration constants is performed with assumption without work hardening of material during rolling process, without forward and backward stretching, and without roll flattening. In the input and output plane of rolling material the stress is . The integration constant for the backward slip zone () is found from (66) as follows:where is a known function given by (65). The integration constant for the forward slip zone () is found from (67) as follows:The total form of equation for the calculation of the relative normal contact stress for the backward slip zone can be written as follows: and the equation for the calculation of the relative normal contact stress for the forward slip zone can be written as follows: The analytical solution of (74) and (75) is problematic because the integrals are very complicated and has no primitive function. Therefore, it is necessary to solve integrals and using numerical methods:In (6) the first integral is calculated analytically from (36) and second integral is calculated numerically from (45). This combined method for solving differential equations can be called as analyticalnumerical method for solving differential equation of the rolling process described in (1).
3.3. Calculation of Neutral Section
The two auxiliary constants and were used for reducing of calculated time of function :Equation (65) is considerably simplified after the introduction of constants:The numerical integration was used for calculation of definite integrals described by (74) and (75).
Graphical visualization of the distribution of the relative normal contact stress in rolling gap depends on relative coordinate and relative deformation is given in Figure 9. From graphical dependence the result is that neutral point N is significantly dependent on the relative deformation. Material output from the rolling gap is represented by a point A and the angle coordinate . Material input to the rolling gap is represented by a point B and the angle coordinate . The influence of friction coefficient on position of neutral point and distribution of the relative normal contact stress in rolling gap is given in Figure 10 with constant value of relative deformation (). If the friction coefficient is increased, the peak value of relative normal contact stress is raising but weak changes were observed in shifting of neutral point depending on relative coordinate . The influence of roll radius on position of neutral point and distribution of the relative normal contact stress in rolling gap is given in Figure 11 with constant value of relative deformation and friction coefficient. Significant influence of the roll radius to the peak value of the relative normal contact stress was observed but weak effect was observed in shifting of neutral point depending on relative coordinate .
Determination of the neutral point N described by coordinate is based on equality of the relative normal contact stress in backward and forward slip zone as follows: Substituting (74) and (75) into (79) will obtain the formulawhere function is given by (65) and function is given by (66). The coordinate of neutral angle is integration limit. The analytical determination of the angle of neutral point resulting from (80) is practically impossible and therefore must be used numerical solution. The numerical solution is based on the step method where interval is scanned with the step . The general solution is represented with the change of variables as follows: and . The conception of numerical solution is shown in Figure 12.
The curves of the relative normal contact stress in the surrounding of point N are substituted by equation of straight lines in the form . If then coordinates of points will be and . The slope of the straight lines ( and ) in interval can be described as follows: Equations of the straight lines in interval can be described as follows:By comparing the equations of the straight lines are obtained refined local coordinate of the neutral point :The improved value of the angle of neutral point N () will be obtained from formulaThe calculation of the neutral angle according to the described algorithm was performed for the following rolling conditions: roll radius mm, thickness of material before rolling mm, relative deformation , and friction coefficient = 0,3. Result convergence of calculation of neutral angle is given in Table 1. The precision of calculation is evaluated by difference of values . Graphical visualization of the angle of neutral point in dependence on relative deformation and ratio is given in Figure 13. Work hardening of rolled material is not considered. The neutral angle increases with rising of relative deformation and decreasing of the ratio . Approximate formula for ratio was presented by the authors of [16].

The validity of formula is shown in Figure 14 where ratio depends on relative deformation and ratio . The local maximum is significantly dependent on the ratio .
The material thickness in neutral point can be determined by the neutral angle as follows:Dividing (85) by output thickness the ratio can be obtained in the form Graphical visualization ratio in dependence on relative deformation and ratio for friction coefficient = 0,1 is shown in Figure 15. The ratio increases with rising of relative deformation and decreasing of the ratio .
If the coefficient of friction so the status can be named as the ideal rolling. This condition can be used for simplified solution of differential equation (4) when it is possible to continue in solving (80) as follows: The value of variable of neutral angle is one of integration limit of function equation (54). The calculation of integral using substitution and their differential will be as follows:Integration in limits for forward slip zone obtains the form and for backward slip zone in limits The neutral angle of ideal rolling process can be calculated from comparison of (89) and (90) as follows: Substituting (51) to (91) and using simplification the next formula will be obtained: Graphical visualization equation (92) is given in Figure 16. The neutral angle is increased with rising of relative deformation and decreasing of ratio . From comparison of dependences shown in Figure 13 () Figure 16 () is the result that neutral angle is slightly higher if friction coefficient . Equation (92) can be considered as simplest form for direct calculation of neutral angle with satisfying results.
3.4. Calculation of Average Normal Contact Stress
The distribution of relative normal contact stress to rolls in plastic deformation zone is shown in Figures 9, 10, and 11 by the curves ANB. The rolling force can be calculated from the average relative normal contact stress along gripping angle (surface of rectangular OEFGO) as is given in Figure 17.
It is supposed the egality of surfaces OANBGO ≡ OEFGO. If it is valid that then the following can be written: The calculation of must be done independently for backward and forward slip zone according to the following formula: Substituting (74) and (75) to (94) is obtained in the following form:The limits of integration for a forward slip zone are defined by interval and calculation of is performed according to (75). Identically it is valid for backward slip zone where and calculation of is performed according to (74). The definite integrals in (95) can be performed by some numerical mathematical method. Graphical visualization of (1) is given in Figure 18. Function throughout the observed interval is monotonically increasing. The relative normal contact stress is increased with rising of relative deformation and decreasing of ratio . Dependence of relative normal contact stress on ratio and relative deformation for three groups of friction coefficient are given in Figure 19. From dependence the result is that the strongest effect on relative normal contact stress has friction coefficient.
For numerical integration of (95) Simpson’s quadrature rule can be used. The author of [24] presented Weddle formula which is more effective than Simpson’s quadrature rule for numerical calculation of general function . The advantage of this numerical method is faster convergence of calculation. The continuous function is numerically integrated in interval where which is divided on final number of subinterval with step as follows:where is the whole positive number:where , , , , , and is the whole positive number (). The calculation of relative normal contact stress in backward slip zoneequation (66) according to Weddle formula is given in Table 2. The precision of calculation is evaluated according to difference of values: . The satisfied precision of calculation is obtained after where .

Boundary conditions of rolling process were given as follows: relative deformation , ratio , friction coefficient = 0,2, and angle coordinate = 0,05 rad. From comparison of numerical calculations the result is that numerical integration according to Weddle formula is four times faster than numerical integration realized according to Simpson’s quadrature rule.
The previous calculations were carried out provided without work hardening of material during rolling process.
In the next part the rolling process with acceptance of isotropic work hardening of material will be described. Alexander [25] mentions Swift described work hardening curves according to formulawhere is the uniaxial yield stress from static tensile test depending on true strain , is the uniaxial yield stress from static tensile test resulting from annealing state of material, and , is the material constants.
The average value of flow stress can be calculated as follows: where is the true strain before actual rolling and is the final true strain after rolling.
For calculation of flow stress the author of [25] used the formula in which effective stress and effective strain are given:Substituting (100) into (99) the following is obtained:and after integration average value of flow stress has a form Let us assume that material is processing by rolling from annealed state for which is valid and transform upper limit of integral , (102) is simplified to form The measurements values of yield stress depending on true strain of laboratory cold rolling experiment are given in Figure 20. For approximation of measurement values of stressstrain curve (98) was used. Final regression equation of stressstrain curve for brass CuZn30 has a formGraphical visualization of (104) is given in Figure 20.
Substituting (104) into (103) the following formula is obtained: Measured and calculated data from experimental cold rolling schedule are given in Table 3. The calculation of deformation resistance was carried out according to formula where is the rolling force and is the contact surface of rolled material with roll.

The average relative normal contact stress calculated from measurement data is given as follows:Average value of flow stress is calculated according to (105). Theoretical value of function is calculated according to (95) with friction coefficient = 0,2. Numeric comparison of measured and theoretically calculated values of function is given in Table 3 and graphical comparison is shown in Figure 21.
The differences between measured and calculated values of the average relative normal contact stress are from interval [%]. These deviations are resulting from the calculation which does not take into account the flattening rolls.
3.5. Calculation of Front and Back Tensions
Back tension is actuating in direction of axis and marked by label . Front tension is actuating in direction of axis and marked by label . According to Figure 1 point B is valid: and tension stress . The tension stresses and can be described by using flow stress as follows:where and are coefficient for back and front tensions, (for without the back tension and for maximal back tension have to be valid ), and (for without the front tension and for maximal front tension have to be valid ).
Using the Tresca condition of plasticity formulae of plasticity of material for backward and forward slip zones will be obtained as follows:Editing of (109) is received: The integration constants and which are describing front and back tensions can be determined from limit conditions given by (111) and (113). The integration constant for forward slip zone and for resulting from (67) will be as follows:The integration constant for backward slip zone and for resulting from (66) will be as follows:The formula for calculation of relative normal contact stress at forward slip zone with consideration of front tension will be as follows:And for calculation of relative normal contact stress at backward slip zone with consideration of back tension will be as follows:The function was described by (65) and function by (64). The distribution of the relative normal contact stress in rolling gap depends on relative coordinate and back tension is given in Figure 22. From graphical dependence the result is that neutral point N is significantly dependent on the back tension. If coefficient of front tension is then front tension stress will be . The back tension stress described by flow stress will be given by coefficient . The value of back tension is characterizing by point B. For example, for point B_{4} being = 0,4 it means that back tension will be as follows: .
If the back tension increases then average relative normal contact stress is decreased. The front and back tensions are not active in neutral point N_{0}. The neutral point is shifting to exit side of rolling gap with increase of back tension.
The distribution of the relative normal contact stress in rolling gap depends on relative coordinate and front tension is given in Figure 23. From graphical dependence the result is that neutral point N is significantly dependent on the front tension. If coefficient of back tension is then back tension stress will be . The front tension stress described by flow stress will be given by coefficient . The value of front tension is characterized by point A. For example, for point A_{4} being = 0,4 it means that back tension will be as follows: .
The common action of back and front tension on the distribution of the relative normal contact stress in rolling gap depends on relative coordinate that is given in Figure 24. The indexing of point is as follows: is the tension stress on enter side of material to rolling gap, is the tension stress on exit side of material from rolling gap. This means that N_{0,0} is representing intersection point of curves without action of front and back tension: ; for N_{6,6} it is for maximal values of front and back tension: .
3.6. Calculation of Friction Coefficient
The simple determination of friction coefficient can be used with graphical dependence given in Figure 25 for constant ratio . More precisely determination of friction coefficient depends on deformation resistance of rolled material. The deformation resistance of rolled material is characterized by threeaxial stress state and depends on friction coefficient between rolls and rolled material. In the case of action only oneaxial stress state is function of friction eliminated and stress is defined as flow stress . If ratio and deformation are constant then the relative normal contact stress in rolling gap depends only on friction coefficient. Mathematical expression describing relationship between and can be given as follows:From (118) the following can be determined:Substituting (119) into (95) and subtraction of (119) from (95) is obtained:Graphical dependence of is given in Figure 26. Intersection point of curve describing function with zero value of axis is the determining value of friction coefficient . The value of friction coefficient valid for given conditions of rolling process. For numerical calculation of friction coefficient Regula Falsi method was used. Maximal value of friction coefficient results from(i)Tresca condition of plasticity: ,(ii)HMH condition of plasticity: .
The first approximation for calculation of friction coefficient using Regula Falsi method will be as follows:The calculation of friction coefficient using Regula Falsi method for rolling conditions given in Table 3; line number 4 (, ratio: ,