See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/278816725

Flatness Defect in Thin Strip Cold Rolling and the Friction Impact on it

Article · June 2012

CITATIONS READS

4 1,176

3 authors:

Rebecca Nakhoul Pierre Montmitonnet

The MathWorks, Inc MINES ParisTech
18 PUBLICATIONS   52 CITATIONS    162 PUBLICATIONS   1,709 CITATIONS   

SEE PROFILE SEE PROFILE

Sami Abdelkhalek
Military Academy, Tunisia
20 PUBLICATIONS   197 CITATIONS   

SEE PROFILE

Some of the authors of this publication are also working on these related projects:

Modelling Sintering process of printed ink View project

Roll Stack deformation models View project

All content following this page was uploaded by Sami Abdelkhalek on 18 August 2015.

The user has requested enhancement of the downloaded file.

 Author manuscript, published in "40th NAMRC - North American Manufacturing Research Conference, Notre-Dame, Indiana :
États-Unis (2012)"

Proceedings of NAMRI/SME, Vol. 40, 2012

Flatness Defects in Thin Strip Cold Rolling

and the Friction Impact on it

Rebecca Nakhoul, Pierre Montmitonnet

CEMEF, UMR CNRS 7635
MINES ParisTech
Sophia-Antipolis, France

Sami Abdelkhalek
IORC Research Centre
ArcelorMittal
Maizières-les-Metz, France
hal-00724786, version 1 - 23 Aug 2012

ABSTRACT
Flatness defects are one of the major problems in strip rolling. They are manifested by a wavy shape on the edge, in the
centre or in between. Waves are most of the time transverse, but all directions can be observed. These defects come from the
heterogeneity of the stress field and the resulting buckling of the compressive areas out of the roll bite. This paper is based on
the approach proposed by [1-3] and [4], and programmed previously [5-7] in the FEM software LAM3/TEC3 [8]. In the
present paper, the latter is enhanced and applied to the impact of friction and strip tension on flatness of a rolled thin strip.
The study shows e.g. that the optimal setting of Work Roll Bending force (WRB) should be changed when friction varies.

KEYWORDS
Rolling, Thin Strips, Finite Element Method, Friction, Flatness Defect.

INTRODUCTION
The origin is the heterogeneous distribution of residual
In thin sheet metal forming processes, buckling may stress: buckling occurs whenever compression exceeds a
occur and result in major defects. This is especially the case certain critical value over a significant area. Of course, this
in rolling of thin strips or sheets, e.g. tinplate in the steel out-of-plane displacement can happen only out of the roll
industry or foil in the aluminium industry. The quality of bite. Therefore, the measurement of this stress profile by
the product is affected by waviness of diverse directions shape-meter rolls (Figure 2) is a central tool in strip shape
and amplitudes, known as flatness defects (Figure 1) [9]. control. Only the distribution of the longitudinal stress in
the transverse direction, σ xx ( y ) , can be measured.

Figure 1. Schematic view of flatness defects during strip

Figure 2. An example of stress profile (“latent flatness”)
rolling.
where T is the front tension stress applied on the strip.
 Proceedings of NAMRI/SME, Vol. 40, 2012

On the rolling mill, a high tension stress is applied on penetrate the roll surface, or on the contrary lose contact
the strip (1/10th to 1/4th of the yield stress typically). It often artificially.
prevents the stress profile from transforming into a wavy Another important point is the thermal – mechanical
shape on line, as it brings the stress in the positive range coupling. Due to the high Peclet number (advection
everywhere. In this case, the defect is “latent”: it is not dominates conduction heat flow), a Streamline Upwind
eliminated, and may show up as soon as the tension force is method is used [15].
cancelled, or upon cutting blanks. It then becomes a
“manifested” defect. The term “latent defect” is therefore
often used in place of “stress profile”, whereas “manifested Mesh generation, contact initialization, boundary conditions
defect” points to the shape of the sheet.
Inversely, for thinner sheets, part of the defect may Velocity and state variables computations (Newton-Raphson)
show even under tension, with generally very local defects.
This paper addresses this situation in priority. A simple Strip temperature computation (SUPG)
approach of buckling has been proposed [1-3] and first
implemented in a Finite Difference Method (FDM) context Roll temperature modelling;
[4], then by two of the present authors [5-7] in a Finite Roll and stand elastic displacement
Element Model (FEM) called Lam3/Tec3. It has allowed a
fully coupled model to be built. The stress field resulting Updating of roll surface
from rolling is computed, with strip plastic deformation and Updating of streamlines and mesh
hal-00724786, version 1 - 23 Aug 2012

roll elastic deformation taken into account, and the out-of- Updating of contact variables
bite stress map is examined for buckling in the same model.
The stress redistribution due to buckling is a result of this
analysis and the effect on the whole system can be studied. Convergence tests
This is different from most of the published, decoupled (loads, temperature, geometry)
approaches [10-12] which address the problem in two steps:
• Evaluation, measurement or computation of the
(post-bite) residual stress resulting from the plastic
deformation of the strip; End
• Semi-analytical or shell FEM modelling of the effect
of this post-bite stress field for a buckling / post- Figure 3. General algorithm of the Lam3/Tec3 FEM strip
buckling analysis of the structure. rolling model.
In the present paper, a summary of the existing model
Lam3/Tec3 is presented first. The buckling model is
described, and an enhanced algorithm correcting a defect of
the method is applied and proven efficient. The coupled
model is used to study the effect of friction on the stress
profile in a cold rolled thin strip, and to show how flatness
actuators (WRB force, strip tension) can be set up for on-
line control of (always possible) variations of friction.

THE FEM MODEL LAM3/TEC3

The basis of the present approach is an implicit FEM-

based rolling model called Lam3/Tec3, developed in the
90’s and early 2000’s. It couples strip and roll stack
deformation models, as described in [8, 13]. Its general
flow chart is given in Figure 3.
For the strip deformation (Lam3), the most salient
feature is a steady state formulation based on streamline
integration to correct the shape for spread, anticipation of
deformation at bite entry etc... This can be considered as a
variant of Eulerian – Lagrangian formulation. A great care
is devoted to the determination of the contact onset and
exit, a difficulty in streamline techniques [14], since, due to Figure 4. The structured mesh used in Lam3. Note axes x =
the space integration, strip surface streamlines may Rolling Direction [RD], y = Transverse Direction [TD], z =
Normal Direction [ND].
 Proceedings of NAMRI/SME, Vol. 40, 2012

r r
As many aspects of the formulation rely on streamlines, n1 .σ .n1 = 0
a structured mesh has been preferred. It is based on 8-node, r r
tri-linear hexahedra, with reduced integration of the n 2 .σ .n 2 > 0 (1)
r r
pressure in the Principle of Virtual Power [16]. Figure 4 n1 .σ .n 2 = 0
illustrates the structured mesh, formed by “extrusion in the
rolling direction” of a rectangular grid of the upstream r r
where n1 and n 2 are the directions of the principal
plane. Structuring the mesh allows a very efficient local
refinement of the mesh, in particular at bite entry and exit. Cauchy stress tensor in the buckled structure (hence the
As residual stresses are essential here, an elastic- third equation). This means that when a tension is applied
viscoplastic constitutive model is used. It is based on in a direction, the membrane is stiff; if the stress becomes
Prandtl - Reuss additive decomposition of strain rate. negative, it gets slack and in fact, the corresponding stress
Jaumann objective derivative is used to write the elastic is put to 0.
model in rate form, and associated von Mises behaviour is The essence of the model consists in determining an
assumed for plasticity. The incremental consistency is extra deformation which elastically brings the stress in the
based on the standard radial return technique. buckled direction back to 0. It may be interpreted as the
In the principle of virtual work, the updated stress is shortening of a material line due to buckling of the
obtained by streamline integration, where the time needed structure. This is more or less analogous to elastic-plastic
for matter to move from an integration point to the next in decomposition, but is activated only out of the roll bite, i.e.
hal-00724786, version 1 - 23 Aug 2012

the streamline is a substitute for time step [8] – since time where buckling may manifest.
does not exist properly speaking in a steady state
formulation. As the pseudo-time step is therefore point- Practical application
dependent due to the adapted mesh, the formulation has
been termed “Generalized Heterogeneous Time Stepping” In the context of small incremental deformation, the
(GHTS). strain tensor is the sum of two components:
The roll stack deformation model is another essential
feature. Like most of the previous ones [17], the single roll
bending and flattening model is based on Timoshenko beam ∆ε = ∆ε el + ∆ε bu (2)
theory, Boussinesq solution for a half-space under general
where ∆ε is the elastic and ∆ε is the “buckling
el bu
loading, combined after the results in [18]. Based on
extensive FEM simulations, corrections have been brought strain” increment. Plane stress is assumed (out of bite). If
for end effect and the barrel / axle transition. Hertzian buckling occurs in direction 1 (respectively 2), the
contact mechanics is assumed for work roll (WR) / back-up following conditions hold:
roll (BUR) contact. The Influence Function Method (IFM)
r r r r
is used to discretize the system, with particular refinement n1 .σ .n1 = σ c n1 .σ .n1 > 0
near the edge of the strip – WR contact. A global non-linear r r r r
system is formed with all displacements of all contact lines, n 2 .σ .n 2 > 0 respectively n 2 .σ .n 2 = σ c (3)
r r r r
with external forces (rolling load, WRB or BURB) in the n1 .σ .n 2 = 0 n1 .σ .n 2 = 0
right-hand side. This non-linear system is solved by
Newton-Raphson method. Details can be found in [13]. The extra deformation representing buckling is
computed in the principal axes then transported to the
A SIMPLE BUCKLING MODEL
reference frame. Let λi , i = I,II be the deformation
Principle representing buckling in the principal directions. It is
deduced from σi, i = I, II as follows:
Counhaye [4] has proposed a method to deal with sheet
buckling in a FDM rolling model, which seems quite σi −σc
similar to the one introduced in a more general context by λi = i = I , II (4)
Roddeman et al. [1]. The same has been implemented in E
Lam3/Tec3 by Abdelkhalek [5,6].
In [1-3], it is proposed for the membrane theory, and Moving back to the reference frame, the buckling strain
forbids the appearance of a negative stress: every time a is added to the global strain increment (u and v are the two
negative stress is about to appear, the structure buckles; this in-plane incremental displacements, θ is the angle between
principal and reference frames, ν is Poisson's ratio and E is
means that σ c = 0 ( σ c is the critical buckling stress). The
Young's modulus):
following critical conditions are therefore introduced:
 Proceedings of NAMRI/SME, Vol. 40, 2012

∂u Recalculating σ (G ) in the same way shows that if θ = 0 :

∆ε 11 = + λ I cos 2 θ + λ II sin 2 θ
∂x
∂v σ (G ) = E∆ε + σ c (10)
∆ε 22 = + λ II cos 2 θ + λ I sin 2 θ
∂y (5)
In this case, the stress is brought back on the buckling
1  ∂u ∂v 
∆ε 12 =  +  + λ II − λ I cos θ sin θ criterion as soon as the strain becomes small. This may not
2  ∂y ∂x  be exactly at exit of the bite, because of a transition zone
ν where the stress reorganizes as the velocity field,
∆ε 33 = − (ε 11 + ε 22 ) heterogeneous at bite exit, evolves to a rigid body
1 −ν
movement [7].
This strain increment replaces the standard one fed into
the module solving the constitutive differential equations.
Table 1. Simulated rolling operation.
An improved algorithm
Width of the metal sheet 855 mm
The results presented in [5,6] show good agreement Entry Thickness 0.355 mm
with available experimental results, i.e. σxx(y) at the shape- Exit Thickness 0.252 mm
meter roll position, about 1 m after the stand (see Figure 8 Back and front tension 170 MPa – 100 MPa
hal-00724786, version 1 - 23 Aug 2012

below). Yet, it has been noticed that the non-buckling Rolling velocity 20.5 m/s
criterion is not respected until ~500 mm after roll bite exit. Work roll crown 0.0322% of radius
The numerical result is therefore locally in contradiction Sheet crown 4.81% of thickness
with the nature of the model [6]. In the following, a simple Bending force 482 KN
1D analysis of the origin is presented. Let G be a Gauss Critical stress σc -10 MPa
integration point. The algorithm implemented by Friction Law Coulomb µ = 0.03
Abdelkhalek [5] is such that: Poisson's ratio ν = 0.3
Young's modulus E = 210 GPa
σ (G ) = σ (G −1) + ∆σ (G ) = σ (G −1) + E.∆ε el (6)

∆εel is the elastic strain increment as a material point

moves from G-1 to G. From equation (2),
∆ε el = ∆ε − ∆ε bu . It was assumed that buckling takes
place at the end of the increment:

σ (G ) − σ c
∆ε bu = (7)
E

Reporting equation (7) into equation (6):

σ (G −1) + E∆ε + σ c
σ (G ) = ≠ σc (8)
2

It shows that this implementation does not lead

immediately to the criterion, but decreases the stress
gradually – which is indeed observed in the results. A
change in the stress integration method along streamlines is
proposed to improve this. Considering that ∆ε bu occurs
during the increment, one may assume:

∆ε bu =
[θσ (G )
]
+ (1 − θ )σ ( G −1) − σ c
(9)
E
Figure 5. The mill structure (4-Hi) and dimensions.
 Proceedings of NAMRI/SME, Vol. 40, 2012

The algorithm has been modified accordingly. To show

the effect of this modification, a rolling operation described
in Table 1 and Figure 5 is simulated (the same is used
throughout the paper). The critical buckling stress is set
to σ c = −10 MPa .
Results are presented as a comparison between the
curves obtained according to equations (6) and (10),
focusing on the post-bite area where the inconsistency
showed. In Figure 6, only one streamline is shown, but an
identical behaviour is found on all streamlines. The new
algorithm now practically respects the criterion everywhere.
These changes do not affect the global shape of the stress
maps, as shown by the stress profile downstream (Figure 7).
Figure 8 shows that this stress profile is in good agreement
with measurements reported in [4].

Figure 8. Comparison between the experimental and

hal-00724786, version 1 - 23 Aug 2012

numerical results (at shape-meter roll position).

APPLICATION: FRICTION AND FLATNESS

As an example of the applicability of this model, the

effect of several parameters on the flatness of a thin strip is
studied hereafter. The goal is not to provide a complete
abacus to run a specific rolling mill, but to highlight
interactions between several actuators.
Friction on the one hand, strip tensions on the other
hand, have been proved to impact strip profile after cold
rolling [19-21]. In [22], the effect of a parabolic variation of
Figure 6. Stress along streamline 19 (edge) after bite exit, friction in the transverse direction has been evaluated. The
before and after the change of integration scheme. Stress following study examines the effects of friction, WRB force
scale is zoomed around σ c = −10 MPa . and strip tension on the shape of the strip. It is shown in
particular that for each friction coefficient, there is a best-
adapted WRB force in the sense of a minimally
heterogeneous stress profile.

Influence of friction on flatness

In this section, a quick study on how friction influences

flatness of a thin strip is presented. The rolling operation
detailed in Table 1 has been modelled using coefficients of
friction between 0.01 and 0.035. Figures 9 and 10 present
respectively the strip thickness profile and the stress profile
in a cross section about 1 m after roll bite exit.
Figure 9 shows that as friction increases, the thickness
of the strip on the edges decreases. Note that a perfect strip
gauge control has been assumed, so that the central
thickness is maintained at its nominal value. The reason for
the growing profile defect is the higher rolling load
Figure 7. Comparison between the two algorithms on a resulting in more WR bending, hence a thicker centre and a
cross section chosen far from the bite (at shape-meter roll larger edge-drop defect. This occurs in spite of the roll
position). crown, which is normally planned to compensate for
bending, but is not adaptable when the roll load varies.
 Proceedings of NAMRI/SME, Vol. 40, 2012

Figure 9. The impact of friction on strip thickness profile.

Figure 11. Roll active generator shape z(y), showing
flattening and bending.
hal-00724786, version 1 - 23 Aug 2012

Figure 10. The impact of friction on stress distribution

(WRB force = 482 kN).
Figure 12. Roll load transverse distribution F(y).

Proceeding with the stress distribution downstream of

the roll bite, Figure 10 shows that changing friction may Finally, the lowest friction µ = 0.01 leads to a dramatic
transform the on-line residual stress distribution in the strip. change in the stress profile: even under strip tension, the
The profiles must be compared with the average tension centre is practically slack, whereas the edges are extremely
100 MPa: wherever the stress is smaller, there is a suspicion tensile: the WRB force 482 kN, added to the roll crown, is
of compression after tension release, therefore of waviness. too large for the lower roll load, the gap is thinner in the
For the central value µ = 0.025, the profile is slightly centre, the elongation of the strip is too large there (this is
concave in the centre, where σxx < 100 MPa. The strip also where the strip is thickest at entry) and a wavy centre
might remain flat on the mill, but a wavy centre is not will obviously result. It can be said that the bending force
excluded after tension release – wavy edges may also occur used is poorly adapted to such a low friction.
simultaneously as the stress there is close to σxx = 0.
For the highest friction, µ = 0.035, the tension is slightly Best-adapted WR bending force for a given friction level
above 100 MPa in the centre, rather constant, and drops
only in the last 100 mm. The curve is concave except near This analysis shows that friction affects the strip
the edge where waviness of limited dimension might occur. thickness profile, the stress distribution, and will therefore
Note that as shown by the roll profiles (Figure 11) and the impact the defect type and amplitude. One of the flatness
roll load transverse profile (Figure 12), “roll kiss” has actuators which can be used to improve flatness is the WRB
occurred (mutual contact of the two WR on either side of force. Its effect, at constant µ, is pictured in Figure 13.
the strip). Increasing it to large values clearly exaggerates counter-
 Proceedings of NAMRI/SME, Vol. 40, 2012

bending, with a too large reduction in the centre resulting in when the latter varies, e.g. during mill acceleration or
a concave profile, leading to a wavy centre. deceleration.
hal-00724786, version 1 - 23 Aug 2012

Figure 15. Evolution of the slip rate as friction changes.

Figure 13. The σxx(y) stress profile for varying bending
force and fixed friction coefficient (µ = 0.025).
The relationship has been established for µ = 0.02 to
0.035. µ < 0.02 is not desirable, due to a risk of skidding
shown by the negative forward slip below µ = 0.015 in
Figure 15. “High” friction, µ > 0.03, leads to 'Roll Kiss'.

The effect of strip tension on the stress profile

Another possible actuator is the strip tension itself. It is

examined in the following if, for a given rolling operation,
it can be used to control – and minimize - the stress
variations σxx(y). Average tension stresses of 50 and 150
MPa have been tested. Curves of Figures 16 and 17 present
respectively the stress profile (average tension T has been
subtracted) and the strip thickness profile in a cross section
∼1 m after roll bite exit.
Figure 14. The relationship between friction and optimal Decreasing the tension (50 MPa) does not change the
bending force. stress profile too much, but increasing it above 100 MPa
clearly affects the stress distribution (Figure 16). As tension
increases, the stress on the edge is kept at the buckling
Thanks to a set of numerical experiments, the relation threshold (-10 MPa), but in relative terms, it is more
between friction, bending force and stress profile σxx(y) has compressive after tension subtraction. Increasing the
been established for the rolling operation investigated. For tension to 150 MPa increases σxx(y) everywhere, but the
each value of the friction coefficient, the “optimal bending profile is almost three times more concave, suggesting a
force” is defined here as the one giving the most flat stress wavy centre upon cancellation of the tension force.
profile (as judged by eye). For instance, from Figure 13, it
can be concluded that 350 kN is the optimal bending force Figure 17 shows that as the tension decreases, the
for µ = 0.025. Inversely, in Figure 10, F = 482 kN gives a thickness of the strip on the edges decreases (the edge-drop
flat profile with µ = 0.03, so that the optimal bending force defect is larger). This is because lower tension results in
for µ = 0.03 is 482 kN. The results obtained for the whole increased roll load, giving more roll flattening, also shown
range of parameters are given in Figure 14. It shows how by the large peak at the edge of the roll load profile (Figure
the bending force could be preset as a function of friction 18). The total rolling force increases linearly from 8380 to
9180 and to 10100 kN as tension decreases from 150 to 100
and to 50 MPa.
 Proceedings of NAMRI/SME, Vol. 40, 2012

Figure 16. The effect of increasing strip tension on the Figure 18. The effect of increasing strip tension on the roll
stress profile (µ = 0.025, WRB force = 482 kN). Strip load transverse distribution F(y) (µ = 0.025, WRB force =
tension stress T has been subtracted for easier comparison. 482 kN).
hal-00724786, version 1 - 23 Aug 2012

CONCLUSION

The problem of flatness defects has been addressed

using a completely coupled model combining the FEM for
strip elastic-viscoplastic deformation, a powerful semi-
analytical model of roll elastic deformation, and a simple
model of buckling based on elastic / buckling deformation
decomposition. The main result addressing flatness is the
computed stress distribution downstream of the roll bite,
together with the buckling indicator λ. The stress map is
shown consistent with the buckling stress criterion. The
stress profile at the position of shape-meter rolls is in good
agreement with measurements reported in [4] (see Figure
8). A study of the effect of friction concludes that, in the
thin sheet (tinplate) case, friction is an essential parameter
for flatness. An adaptive set up of the bending force is
shown to be able to compensate for unavoidable friction
variations: to each coefficient of friction corresponds a
bending force ensuring the best possible flatness. For the
time being, the expected optimal flatness is judged only by
the homogeneity of the stress profile, but this criterion
could be refined in the future.
Front tension has been studied as well. The effect of an
increased tension is positive on the roll load and roll
deformation, but raises a flatness issue. Probably, this could
be again partly compensated by changing simultaneously
the WRB force. Playing with two actuators together
however raises the question of the local / global character of
their effects.
In this study, the effect of ∆εth, the thermal dilatation
Figure 17. The effect of increasing strip tension on the strain, has not been investigated. In a thin strip cold rolling
thickness profile of the strip (µ = 0.025, WRB force = 482 problem, the temperature increase in the roll bite may reach
kN). 100 K, and, due to the differential reduction, may not be
homogeneous. Cooling may also be heterogeneous, with
 Proceedings of NAMRI/SME, Vol. 40, 2012

edges being cooled faster. This may bring in thermo-elastic [10] Bush A, Nicholls R, Tunstall J (2001) Stress levels for
stresses, which will be the subject of future investigations. elastic buckling of rolled strip and plate. Ironm. Steelm.
28, 481-484.

ACKNOWLEDGEMENT [11] Rammerstorfer FG, Fisher FD, Friedl N (2001)

Buckling of free infinite strips under residual stress and
The authors wish to thank the French National global tension. ASME Trans., J. Appl. Mech. 68, 399-
Research Agency (ANR) for its financial support and the 404.
partners of the Platform Project for the authorization to
publish this work. [12] Fisher FD, Rammerstorfer FG, Friedl N, Wisser W
(2000) Buckling phenomena related to rolling and
levelling of sheet metal. Int. J. Mech. Sci. 42, 1887-
REFERENCES 1910.

[1] Roddeman DG, Drukker J, Oomens CWJ, Janssen JD [13] Hacquin A, Montmitonnet P, Guillerault JP (1998) A
(1987) The Wrinkling of Thin Membranes: Part I- 3D semi-analytical model of rolling stand deformation
Theory. ASME Trans., J. Appl. Mech. 54, 884-887. with finite element validation. Eur. J. Mech. A (Solids)
17, 1, 79-106.
[2] Roddeman DG, Drukker J, Oomens CWJ, Janssen JD
hal-00724786, version 1 - 23 Aug 2012

(1987) The Wrinkling of Thin Membranes: Part II- [14] Yamada K, Ogawa S, Ataka M (1992) 3D analysis of
Numerical Analysis. ASME Trans., J. Appl. Mech. 54, flat rolling using rigid-plastic FEM coupled with roll
888-892. deformation analysis. Proc. Numiform 92, Chenot JL,
Wood RD & Zienkiewicz OC, eds, Balkema,
[3] Roddeman DG (1991) Finite Element Analysis of Rotterdam, 1992.
Wrinkling Membranes. Comm. Appl. Numer. Meth. 7,
299-307. [15] Pardo E, Weckman DC (1990) A fixed grid FE
technique for modelling phase change in steady state
[4] Counhaye C (2000) Modélisation et contrôle industriel conduction – advection problems. Int. J. Num. Meth.
de la géométrie des aciers laminés à froid. PhD Thesis, Engg 29, 969-984.
Liège University (in French).
[16] Gratacos P, Montmitonnet P, Chenot JL (1992) An
[5] Abdelkhalek S (2010) Un exemple de flambage sous integration scheme for Prandtl-Reuss elastoplastic
contraintes internes: Étude des défauts de planéité en constitutive equations. Int. J. Num. Meth. Eng. 33
laminage à froid des tôles minces. PhD Thesis, MINES (1992) 943-961.
ParisTech (in French).
[17] Shohet KN, Townsend NA (1968) Roll bending
[6] Abdelkhalek S, Montmitonnet P, Legrand N, Buessler P methods of crown control in four-high plate mills. J.
(2011) Coupled approach for flatness predictions in thin Iron Steel Inst. 1088-1098.
strip cold rolling. Int. J. Mech. Sci. 53, 651-675.
[18] Berger B, Pawelski O, Funke P (1976) Die elastische
[7] Abdelkhalek S, Zahrouni H, Montmitonnet P, Legrand Verformung der Walzen von Vierwalzengerüsten. Arch.
N, Potier-Ferry M (2011) Manifested flatness Eisenhüttenwes. 47, 6, 351-356 (in German).
predictions in thin strip cold rolling using a general
rolling FEM model. Steel Res. Int. Sp. Issue ICTP 2011, [19] Jiang ZY, Tieu AK, Zhang XM, Lu C, Sun WH
111-116. (2003) Finite element simulation of cold rolling of thin
strip. J. Mat. Proc. Tech. 140, 542–547.
[8] Hacquin A, Montmitonnet P, Guilleraut JP (1998) A
steady state thermo-elastoviscoplastic finite element [20] Jiang ZY, Zhu HT, Wei DB, Tieu AK (2007)
model of rolling with coupled thermo-elastic roll Asymmetric Cold Rolling of thin Strip with Roll Edge
deformation. J. Mat. Proc. Tech. 60, 109-116. Kiss. Proc. ICIEA 2007 (Harbin, PR China, May 23-25,
2007), 2781-2786, published by IEEE.
[9] Roberts WL (1983) Hot Rolling of steels.
Manufacturing Engineering and Materials Processing [21] Wang XC, Yang Q, Du XZ, Jiang ZY (2010)
Series, Vol. 10. Marcel Dekker, New York. Allowable variation of cold-rolled strip transverse
profiles in high tension. Int. J. Minerals, Metallurgy
and Materials 17, 5, 608-616.
 Proceedings of NAMRI/SME, Vol. 40, 2012

[22] Jiang ZY, Xiong SW, Tieu AK, Wang QJ (2008)

Modelling of the effect of friction on cold strip rolling.
J. Mat. Proc. Tech. 201, 85-90.
hal-00724786, version 1 - 23 Aug 2012

View publication stats