- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Thermo-mechanical modelling of hot flat-rolling of...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Thermo-mechanical modelling of hot flat-rolling of components with curved profiles 1993
pdf
Page Metadata
Item Metadata
Title | Thermo-mechanical modelling of hot flat-rolling of components with curved profiles |
Creator |
Xiao, Ming |
Date Created | 2008-08-27 |
Date Issued | 2008-08-27 |
Date | 1993 |
Description | A complete mathematical model for implementing the VGR process (a flat rolling process for producing components with variable thickness) under hot working conditions has been established. The work is based on an idea that a combination of separate submodels, such as, deformation, flow stress, roll force, temperature, etc., should be used, so that new development in each of these areas can be easily incorporated. A deformation submodel, based on the upper bound theorem, is given for the analysis of three-dimensional deformation of the work piece. A simple velocity field is proposed. To preserve the theoretical consistency, an equivalent coefficient of friction is adopted for the roll force calculation. Furthermore, the basic assumption, i.e., the rigid perfectly-plastic material assumption, is modified by introducing a concept of an isotropic rate-dependent material. Satisfactory results were obtained in spread, torque and force prediction. To obtain the mean temperature of the work piece, a temperature submodel is formulated based on one-dimensional transient flow in the roll-bite and two-dimensional flow outside the roll-bite. The model is capable of predicting the through-thickness temperature distribution in the roll-bite, estimating the mean temperature of the deforming body, and roughly calculating the mean temperature distribution along the work piece. To characterize the high temperature behaviour of steels, the well-known unified creep relationship is chosen as the flow stress submodel. Reasonably accurate prediction of the flow stress is achieved by using some experimental data reported in the literature. |
Extent | 4019281 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | Eng |
Collection |
Retrospective Theses and Dissertations, 1919-2007 |
Series | UBC Retrospective Theses Digitization Project |
Date Available | 2008-08-27 |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0081027 |
Degree |
Master of Applied Science - MASc |
Program |
Mechanical Engineering |
Affiliation |
Applied Science, Faculty of |
Degree Grantor | University of British Columbia |
Graduation Date | 1993-11 |
Campus |
UBCV |
Scholarly Level | Graduate |
URI | http://hdl.handle.net/2429/1546 |
Aggregated Source Repository | DSpace |
Digital Resource Original Record | https://open.library.ubc.ca/collections/831/items/1.0081027/source |
Download
- Media
- ubc_1993_fall_xiao_ming.pdf
- ubc_1993_fall_xiao_ming.pdf [ 3.83MB ]
- Metadata
- JSON: 1.0081027.json
- JSON-LD: 1.0081027+ld.json
- RDF/XML (Pretty): 1.0081027.xml
- RDF/JSON: 1.0081027+rdf.json
- Turtle: 1.0081027+rdf-turtle.txt
- N-Triples: 1.0081027+rdf-ntriples.txt
- Citation
- 1.0081027.ris
Full Text
We accept this thesis as conforming to the required standard THERMO-MECHANICAL MODELLING OF HOT FLAT-ROLLING OF COMPONENTS WITH CURVED PROFILES By Ming Xiao B. ASc. (Electronic-Mechanical Engineering) Guilin Electronic Industry College, 1984 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF APPLIED SCIENCE in THE FACULTY OF GRADUATE STUDIES MECHANICAL ENGINEERING THE UNIVERSITY OF BRITISH COLUMBIA June 1993 © Ming Xiao, 1993 In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for extensive copying of this thesis for scholarly purposes may be granted by the head of my department or by his or her representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. Mechanical Engineering The University of British Columbia 2075 Wesbrook Place Vancouver, Canada V6T 1Z1 Date: / 9 9 3 ABSTRACT A complete mathematical model for implementing the VGR process(a flat rolling process for producing components with variable thickness) under hot working conditions has been established. The work is based on an idea that a combination of separate submodels, such as, deformation, flow stress, roll force, temperature, etc., should be used, so that new development in each of these areas can be easily incorporated. A deformation submodel, based on the upper bound theorem, is given for the analysis of three-dimensional deformation of the workpiece. A simple velocity field is proposed. To preserve the theoretical consistency, an equivalent coefficient of friction is adopted for the roll force calculation. Furthermore, the basic assumption, i.e., the rigid perfectly-plastic material assumption, is modified by introducing a concept of an isotropic rate-dependent material. Satisfactory results were obtained in spread, torque and force prediction. To obtain the mean temperature of the workpiece, a temperature submodel is formu- lated based on one-dimensional transient flow in the roll-bite and two-dimensional flow outside the roll-bite. The model is capable of predicting the through-thickness temper- ature distribution in the roll-bite, estimating the mean temperature of the deforming body, and roughly calculating the mean temperature distribution along the workpiece. To characterize the high temperature behaviour of steels, the well-known unified creep relationship is chosen as the flow stress submodel. Reasonably accurate prediction of the flow stress is achieved by using some experimental data reported in the literature. 1 1 Table of Contents ABSTRACT^ ii List of Figures^ vii Notation^ xi ACKNOWLEDGEMENT^ xiii 1 INTRODUCTION 1 1.1 Deformation Analysis of Rolling Process — An Overview ^ 1 1.2 Manufacturing of Single Leaf Springs ^ 2 1.3 Aim and Approach ^ 4 2 BASIC CONSIDERATIONS 6 2.1 Available Analytical Methods ^ 6 2.2 Assumptions ^ 10 2.3 Coordinate Systems ^ 12 3 CONSTRUCTION OF VELOCITY FIELDS 14 3.1 The Dual-Functional Method ^ 14 3.2 A Brief Derivation of the Solution to the Velocity Fields ^ 15 3.3 The Derivation of the Velocity Fields ^ 16 3.4 Discussion of the General Solution 20 3.5 Further Considerations ^ 22 Hi 4^3.5.1^Theoretical Aspect ^ ^ 3.5.2^Technical Aspect FORMULATION THROUGH UPPER BOUND APPROACH 22 23 26 4.1 The Strain Rates ^ 27 4.2 The Energy Rates 28 4.2.1^The energy dissipation rate due to deformation(Ed) ^ 28 4.2.2^The energy dissipation rate due to friction between roll surfaces and workpiece(Ef) ^ 29 4.2.3^The energy dissipation rate due to velocity discontinuities (E,) 30 4.3 The Minimization Procedure ^ 31 4.3.1^The derivatives of the geometry ^ 32 4.3.2^The neutral point ^ 32 4.3.3^Solutions ^ 33 5 FLOW STRESS MODEL 34 5.1 Dynamic Flow Stress ^ 34 5.2 Incorporation of Rate Dependency ^ 36 5.3 Mean Values of Strain and Strain Rate 38 5.3.1^Mean Strain ^ 38 5.3.2^Mean Strain Rate 40 5.4 Temperature Prediction ^ 41 5.4.1^Basic Considerations 41 5.4.2^Heat Conduction Model ^ 43 5.4.3^Assumptions ^ 43 5.4.4^Initial and Boundary Conditions ^ 46 5.4.5^Solution to the Heat Conduction Equation ^ 47 iv 5.4.6 Heat Transfer Coefficient at Workpiece-Roll Interface ^ 49 5.4.7 Convection and Radiation ^ 52 5.4.8 Mean Temperature of the Deforming Body in the Roll-Bite^53 5.5 Flow Stress Models ^ 54 5.5.1 Model for Steels ^54 5.5.2 Using Lead as a Modelling Material ^ 55 6 MODEL VERIFICATION^ 57 6.1 Flow Stress Model Verification 57 6.1.1 Empirical Constants in Flow Stress Model ^ 57 6.1.2 Determination of Flow Stress at any Desired Conditions ^ 60 6.2 Temperature Model Verification ^ 61 6.2.1 Verification Using Data from Conventional Rolling ^ 61 6.2.2 Simulation of First Pass in VGR Process ^ 62 6.3 Deformation Model Verification ^ 63 6.3.1 The Side Spread Estimation 63 6.3.2 Torque Estimation ^ 64 6.3.3 Roll Separating Force Estimation ^ 64 6.4 Concluding Remarks ^ 66 6.4.1 Summary and Conclusions ^ 66 6.4.2 Suggestions for Further Work 68 Bibliography^ 77 Appendices^ 82 A ANALYTICAL SOLUTION TO THE UPPER BOUND EQUATIONS 82 A.1 The neutral point ^ 82 A.1.1 The derivative of Ed against w1 ^ 82 A.1.2 The derivative of Ef against w1 ^84 A.1.3 The derivative of t, against w1 ^85 A.1.4 The derivatives of td, E1 and ti, against xn ^ 86 B TRUE STRESS & STRAIN CORRELATION FOR C-Mn STEELS 88 vi List of Figures 1.1 The schematic illustration of the profiles of single leaf springs. a. Tapered profile. b. Parabolic profile. c. Arbitrary profile ^3 2.1 The schematic illustration of ring-rolling configuration ^7 2.2 (a) The illustration of fishtail-shaped sides in ring-rolling. (b) The illus- tration of the slip-line field for flat-end tool forging. ^8 2.3 Schematic illustration of the part of a bar in rolling ^11 2.4 The coordinates system ^13 3.1 Schemetic illustration of the top view of a rectangular bar between the rolls 18 3.2 Side views(a, c) and plan views(b, d) of 1S alloy specimens deformed at 40% reduction; height: 26 mm, width: 26 mm ^ 19 3.3 The schematic illustration of the simplified side profile of the section be- tween the rolls ^20 3.4 The schematic illustration of a parabolic profile in rolling configuration ̂ 24 4.1 The schematic illustration of the velocity vector space in conventional rolling process ^30 5.1 Schematic flow curve for metals at elevated temperature. From point 0 to 1, the plastic strain rates increase rapidly. Note that the scale for the strain is arbitrarily chosen from typical C-Mn steels[27] to illustrate a possible range in applications ^35 5.2 Schematic setting of plane compression test ^38 vi i 5.3 Schematic relation between mean yield stress and mean equivalent strain in plane compression test ^ 39 5.4 Thermal conductivity as a function of temperature for 3 carbon content 0.08 C, 0.4 C, 0.8 C in austenite phase. (From Devadas, 1989) ^ 45 5.5 Specific heat as a function of temperature for 3 carbon content 0.08 C, 0.4 C, 0.8 C in austenite phase. (From Devadas, 1989) ^ 46 5.6 Thermal conductivity as a function of temperature for C-Mn steels .^47 5.7 Specific heat as a function of temperature for C-Mn steels ^ 48 5.8 Example for heat transfer coefficient at workpiece-roll interface (lubricant: water) ^ 50 5.9 Flow stress model: the correlation between yield-stress and reduction ratio of pure lead. Data source: [56] ^56 6.1 Relationship between ln(sinh(ao-)) and ln e for a C-Mn steel (0.03 C 0.62 Mn) at a strain of 0.35. Experimental data: o — (i = 2s-1), * — (i =- 20s-1) and + — (e = 140s-1). Data source: [53] ^58 6.2 Relationship between ln(sinh(ao-)) and 1/T for a C-Mn steel (0.03 C 0.62 Mn) at a strain of 0.35. Experimental data: o — (i = 2s-1), * — (i = 20s-1) and + — (i = 140s-1). Data source: [53] ^59 6.3 Flow stress model verification: C-Mn steel (0.03 C 0.62 Mn). Data source: [53]. +, *, x — experimental data. Solid or dotted lines — model predictions. 69 6.4 Temperature distribution in roll-bite. Data source:[37]. Top: 3D-mesh plot. Middle: grid used. Lower left: front view of the 3D-mesh. Lower right: surface temperature comparison. Reduction ratio: 35%. Entry temperature: 1012 °C. ^70 6.5 Sample calculation: mean temperature distribution of the workpiece for three different rolling conditions. Workpiece dimension(before rolling): 21.6 x 21.6 x 800mm. Reduction ratio: 8% - 50.8%. Furnace temperature: 1200 °C. Top(case A), time for passl: 6.015 sec. Middle(case B), time for passl: 3.008 sec. Bottom(case C), time for passl: 3.008 sec. 71 6.6 Sample calculation: effects of changes in geometry and roll velocity on mean temperature of the deforming body in the roll-bite during the first pass. Furnace temperature: 1200 °C. Workpiece dimension(before rolling): 21.6 x 21.6 x 800mm. Reduction ratio: 8% - 50.8% 72 6.7 Spread prediction - comparison between simulation results and experi- ments[551. Roll diameter: 101.6 mm(4 in). Specimen dimension: 9.525 x 9.525 mm(0.375 x 0.375 in) 73 6.8 Torque prediction: comparison between simulation results and experi- ments[57]. Roll diameter: 127 mm(5 in). Specimen dimension: 12.7 x 19.05 mm(0.5 x 0.75 in) 74 6.9 Schematic illustration of the roll force vectors layout in flat rolling process 75 6.10 Roll force prediction: comparison between simulation results and experi- ments[57]. Roll diameter: 127 mm(5 in). Specimen dimension: 12.7 x19.05 mm(0.5 x 0.75 in) 76 B.1 True- stress and strain are experimental data B.2 True- stress and strain are experimental data B.3 True- stress and strain are experimental data ^ relationships for 0.01 C 0.19 Mn steels. 89 90 91 relationships for 0.03 C 0.62 Mn steels. relationships for 0.19 C 0.64 Mn steels. ix B.4 True- stress and strain relationships for 0.38 C 0.64 Mn steels. * 0 X2^2^2 + are experimental data ^ 92 x Notation 0) 1 ^subscripts for entrance and exit, respectively i) f ^subscripts indicating the initial and final conditions surface area of the workpiece, m2 Stefan-Boltzmann constant, 5.67x 10- 8 wm-2K-4 control parameter in the roll gap function specific heat of the workpiece, Jkg-1K-1 Ed, E1, E3 energy dissipation rates due to deformation, friction and shear, respectively roll separating force FT^frictional force at the interface between rolls and workpiece FR^overall radial force on the roll G(C)^the roll gap function convection heat transfer coefficient h(x)^the thickness function H(r)^heat transfer coefficient J.^upper-bound on power the shear yield stress k,^thermal conductivity of the workpiece, kWm-2K-1 1^the length of the projected arc in contact friction shear factor heat generation the prescribed flow factor xi Q.^activation energy, kJ/mol the reduction ratio: r = (to — ti )/t0 Rg^the universal gas constant, 8.31 J/mol K roll radius the thickness function the rolling torque constant temperature, with different subscripts (used in Chapter 5) T,.^roll surface temperature, °C T,^surface temperature of the workpiece, °C T(x, 7-)^temperature field in workpiece, °C ^ volume of the workpiece VR^roll peripheral velocity Vx, Vy, V,^velocity components in the x, y, z directions w(x)^the width function Xn^neutral point AV^velocity discontinuity emmisivity of the workpiece, Jkg-1K' strain rate tensors mean strain mean strain rate /leg^equivalent coefficient of friction P.^density of the workpiece, (for steel, 7600 kgm-3) the flow stress of the material To^the friction shear stress time variable xi i ACKNOWLEDGEMENT I would like to express my sincere thanks to Professor F. Sassani for taking the author as a graduate student and providing the inspiration for this work. His valuable discussions and helpful suggestions throughout the project are particularly apperciated. I would like to thank Professor I. Yellowly for his suggestion of using the upper bound approach in this work at first place, and his stimulating discussions at the times when I was in need of help. Financial asistance from NSERC is gratefully acknowledged. I would also like to thank computer system managers, G. Rohling and A. Steeves, for their help in using the computers in the Department. I wish to express my appreciations for my fellow students (also my friends) in the Department for their constant interactions and encouragement, in particular, Ramin Ardekani and Fariba Aghdasi. I take this opportunity to express my gratitude for the support and encouragement shown by my parents, to which I am in debt for a lifetime. TO MY PARENTS xiv Chapter 1 INTRODUCTION 1.1 Deformation Analysis of Rolling Process — An Overview The history of the rolling process can be traced to more than four hundred years back, whereas, the theoretical analysis, mainly within the plane-strain assumption, of the pro- cess began much later, probably in the 1920's[1]. Apparently because the plane-strain analysis was much simpler and was able to meet the requirements in steel production, the most widely known and used theoretical models to date were derived by assuming plane-strain deformation, such as Orowan's model [1] and Von Karman's equation [1). Usually, all the earlier passes in the hot rolling of slabs were scheduled by means of empirical methods, and only the few final passes were of greater concern, where the parameters of the products were obtained. In such situations, the plane-strain assumption worked very well, since in most cases, the width-to-thickness ratio of the material in rolling is much larger than 10 and thus the spreads in width can be ignored. It was observed in the operation of the strip mills that the spreads in width seldom exceed 1 to 2 percent [2]. With increasingly stringent demands for accuracy, innumerable efforts have been put into the development of theoretical models for flat rolling during the last seven decades. Almost all the successful comprehensive(namely, various factors were included, and usu- ally also combined hot and cold rolling) models for the strips and plates rolling were 1 Chapter 1. INTRODUCTION^ 2 evolved from plane-strain assumption, not only because it was most reasonable, but also because the resulting models were readily adaptable for engineering use and rather con- venient for the programming of control computers. It would not be an exaggeration that if one would state that the plane-strain analysis in the rolling process has culminated with the help of the rapid advances in the computer technology in the recent years. Or in other words, if one were seeking a solution to a plane-strain problem in the traditional rolling process, where the roll-gap is constant, most likely he would find out, amongst the vast literature on the related subjects, that the work has already been done for some time. Although the research and development of the rolling process have reached such a encouraging stage, however, all the theoretical works, including quite a few successful attempts made in the last three decades, where three dimensional deformation analyses were worked out, say, [3-6], were originated from constant roll-gap geometries, for a uniform thickness of the products has been one of the primary goals in flat rolling. Having considered this situation, one may realize that a detailed review on the literature would not be helpful, and therefore is omitted here. 1.2 Manufacturing of Single Leaf Springs At the beginning of 70's, the automobile industry in Britain began to make tapered leaf springs using the flat rolling process, where the roll-gap was designed as a variable so that it could be changed continuously while rolling and produced a tapered leaf, which had thickness variations along its longitude, figure 1.1a. Ever since, the world-wide demand for the tapered leaf springs has been increasing. The distinct advantage of such a process is its flexibility, in producing an arbitrary set of parameters of the tapered shape, that Chapter 1. INTRODUCTION ^ 3 can not be achieved by any other tooling structure. It has been well-known for a long time that the most efficient type of leaf springs is those with a parabolic profile [7], figure 1.1b, but the preparation for such shaped materi- als is much too difficult. The reason that the industry wound up in the manufacturing of the tapered leaf springs is that the tapered leaf is a compromise between the difficulties involved in the manufacturing and its efficiency. C Figure 1.1: The schematic illustration of the profiles of single leaf springs. a. Tapered profile. b. Parabolic profile. c. Arbitrary profile. Although, in the last two decades, a number of mills were set up around the world specifically for making of the tapered leaf springs, no theoretical work on the subject has appeared in the literature. The only research work on this subject that one may find out in the literature, [8], [9] and [10], was based on the empirical formula for side spreads prediction, which was originated from factor-analysis in 1950's, proposed by Sparling[11] in 1961 and improved by El-Kalay and Sparling[12] in 1968. It is well-known that empirical formulae are often Chapter 1. INTRODUCTION^ 4 simple and handy to use. Chitkara and Johnson[13] found out that the formula gave reasonable results within the range of experimental conditions under which they were devised, "however, their accuracy under different conditions is open to question." The present author considers that it could be a good idea to take advantage of the formula, if the purpose of the research is only to give fair estimates in the deformation, that is, the spreads and the elongation. As a matter of fact, good results in their work were reported[8][9][10]. However, the formula may not be adequate for machine set-up and process design, where a consistent model is needed in order to calculate deformation along with the force and torque variations during production. Moreover, temperature effects are dominant in hot rolling, to which much attention should be given. It would be most desirable to have a model capable of incorporating the most recent developments, say in tribology and material science, so that a thorough knowledge of the deformation and force of the process may be obtained, thereby confidence may be gained for the design purpose. After a careful study of the important previous works on the flat rolling process, the author believes that any reasonable theoretical analysis of the variable gap rolling (VGR) process' has to start with the ever-changing geometry of the workpiece, and one has to give up the idea of making a direct use of any available theoretical models in flat rolling analysis. 1.3 Aim and Approach As mentioned in Section 1.1, only in cases where the spread in width is negligible, the problems can be considered in plane-strain. 'From now on, for the sake of convenience, the rolling process with a variable roll-gap geometry will be referred to as the VGR process. Chapter 1. INTRODUCTION^ 5 It is conceivable that two schemes may be used in the design of the production rolling mills for leaf springs. One is that the rolling is to start with a rectangular bar, of which the width-to-thickness ratio may be 1 to 3, as those required by some already existing production mills for tapered leaf springs. In this scheme, the leaf spring is produced one at a time, and the main problems are that the width of leaf must be constant and also the thickness variation is usually from 30% to 70%[14]. Apparently, in such cases the deformation analysis has to be considered in three dimensions, except when considering edge-rolling for the purpose of width correction(see Section 2.1). The other scheme is to start with a wide slab, the width-to-thickness ratio may be larger than 10. In this case, the slab is to be rolled to the desired dimensions and more than 10 leaves can be obtained at a time by cutting the rolled blank. Accordingly, the spread in width is at least no longer a serious problem, and the process may be considered in plane-strain. However, the machine set-up in this case may be much more costly. Besides, at present, cutting process is usually detrimental to the final properties of the products in terms of their strength and fatigue resistance, and therefore is not preferable in leaf spring production. The present work is aimed to propose a generic theoretical model, three dimensional strains will be considered. The model should be capable of giving a good description of the variation in deformation, force and torque, under conditions of hot working and varying geometry. Thus, the model may also be used in producing such components shown in figure 1.1c. Chapter 2 BASIC CONSIDERATIONS 2.1 Available Analytical Methods In the past half century, a number of approximate analytical methods have been devel- oped and applied to the rolling process, amongst which the most well known are the following. 1. The slab analysis This method is very restricted by its nature and is usually used for a rough pre- diction of the required load, with a reasonable stress distribution; but it is not a suitable method for a rigorous deformation analysis. 2. The slip line field Since its derivation is based on the plane-strain assumption for perfectly plastic ma- terials, generally speaking, it is not reasonable to apply this method to three dimen- sional problems. In addition, the construction of slip-line fields is still quite limited in predicting results that offer good correlations with experimental evidence[16]. One of the applications of the method that should be mentioned here is the one used by Hawkyard et al [9] in ring-rolling. See figure 2.1. Despite its formidable- seeming complexity, that is, ring-rolling has the similarities with asymmetric rolling processes, it could still be tackled since the work-roll feed rate is generally small and is a known constant. It is worth noting that the reduction is usually very 6 Chapter 2. BASIC CONSIDERATIONS^ 7 Roll velocity Feeding velocity Figure 2.1: The schematic illustration of ring-rolling configuration small for each turn and the ratio of 1/h is usually large, figure 2.1. These lead to such an effect that the plastic deformation between the rolls is far from penetrating the ring thickness, and instead of the usual bulges in rolling, fishtail-shaped sides result, figure 2.2a. Based on the above observations, they defined the plastic deformation between the rolls by the kind of slip-line field applicable to flat-end rigid tools, which was developed by Hill [17], when forging a block whose depth is larger than that of the tool width, figure 2.2b, was considered. Apparently, such a method is not applicable to the deformation in the present process, but it may be a good choice for the edge-rolling, which is usually used to obtain a uniform width. Chapter 2. BASIC CONSIDERATIONS^ 8 (a) (b) Figure 2.2: (a) The illustration of fishtail-shaped sides in ring-rolling. (b) The illustration of the slip-line field for flat-end tool forging. 3. The visioplasticity This approach is experimental and could be useful if theoretical or numerical rela- tions cannot be worked. 4. Hill's general method[3] This method is derived from the virtual work-rate principle, and will be cited in the related part of the next chapter. 5. The finite -element method(FEM) The applications of this method to the analysis of three dimensional rolling process were started in early 80's [18]. One of the main advantages of the method is that it can provide detailed information, point to point, of the deforming body between the rolls, which, however, requires a step-by-step computational procedure and Chapter 2. BASIC CONSIDERATIONS^ 9 therefore, requires a much higher data storage and much longer computing time. Fortunately, the conventional rolling process is a steady-state process, and only a few samples of the deformation analysis are needed to obtain a good prediction of the whole process. The formidable difficulty in the analysis of the VGR process is that the pattern of the deformation between the rolls changes continuously throughout the process, and inevitably, a calculation for every point along the product has to be done. For instance, to roll a one-meter long component, the step to generate a set of reasonable interpolation data for the digital controller of the machine is 0.0005 meter, then, at least 2000 steps is needed (for one pass), which means that if finite- element method were to be used, the demands for data storage and computing time would increase a thousand times! Incidentally, Grober [19] reported in 1986 that a simulation of transient hot rolling of steel was performed, where a CPU time of 7 hours on a microVAX II computer was consumed for the simulating of a total length of 51.25 mm, and only a two-dimensional finite-element model was used. This, most likely, were the very reason that no significant theoretical work on the subject, ie., the VGR process, by using FEM (nor by any other method) was found in the English publications(a rather great amount of efforts has been put into the search). Conscious of the foreseeable great demands for the computation, the method is not to be used in the present problem. 6. The upper bound approach (UBA)1 The upper bound approach is valuable to mechanical or production engineers, since it offers rather good results in the roll mill operations with lower costs, and from the viewpoint of machine design, the data produced by the method, though somewhat 1This name was first used by W. Johnson and was considered being well-expressed[55]. Chapter 2. BASIC CONSIDERATIONS^ 10 approximate in nature, are certainly on the safe side. It can also be found in the literature that many researchers used this approach in conjunction with other methods, such as least squares, finite-element, viscoplasticity and the weighted residuals, etc., with which the difficulties to obtain proper velocity fields may be reduced. Perhaps the earliest application of the UBA to the three dimensional analysis of the conventional rolling process was made by Oh and Kobayashi in 1975[5], where quite accurate solutions to the average spread and the torque requirements were given. By extending the work of Oh and Kobayashi[5], Lahoti, et al, developed a production software package for the rolling of plates and airfoils in 1978 [20], and Sevenler, et al, managed to analyze shape rolling with flat rolls [21]. With the method, bulge predictions were also attempted by Kennedy[6] in 1986, though fairly noticeable difference between the theoretical predictions and the experimental results, cited in the work, exists. The major difference between the model by Oh and Kobayashi[5] and the one by Kennedy[6] is that two totally different solutions to the velocity fields were used. Considering all the facts involved, it is reasonable to use the upper bound approach in the analysis of the present problem. 2.2 Assumptions The present analysis is based on the mathematical consideration proposed by C. S. Yih [22] in 1957. To be practical, the following assumptions have to be made. 1. The material is isotropic and incompressible. 26 î I I- / / / / /, 2 t1 2w , Chapter 2. BASIC CONSIDERATIONS^ 11 Z^2t Figure 2.3: Schematic illustration of the part of a bar in rolling 2. The material assumes perfect rigid-plastic properties, that is, the flow stress of the material a is constant2. 3. The deformations meet Mises yield criterion and assume unrestricted plastic flow, therefore, the Levy-Mises equations can be used. The inertial forces are negligible, since in this type of rolling lower speeds are used. 4. The friction shear stress To is constant in the roll-workpiece interface and is given by a constant friction shear factor m defined by the relation, To = m • a/./a. 5. The bar to be rolled can be considered as a set of flow solenoids with different rectangular cross sections, figure 2.3 shows one such element. For the sake of simplicity, any line in z-direction is assumed to remain straight. Since one can never have the exact solution to the flow patterns in question, good approximation has to be dependent upon two properties — theoretical tenability and mathematical feasibility. In order to avoid overwhelming mathematical difficulty, it would be 2This assumption is only based on the strict restraint in the upper-bound theorem. In fact, in real applications, this restriction has to be relaxed and the concept of instantaneous yield conditions has to be adopted. Detailed discussion may be found in Sections 5.2 and 5.3. Chapter 2. BASIC CONSIDERATIONS^ 12 better to give up the bulge prediction at this stage. 3 6. The profile of the product can be expressed in a continuous function or may be defined in several continuous intervals with different functions, and the slopes of the function, have to be very small with no sharp change on the boundaries of the intervals, if any. 2.3 Coordinate Systems It can be seen that the cartesian coordinate system,(x, y, z), is convenient in developing the velocity fields, figure 2.4. The x, y and z directions are defined along the length, width and thickness of the product, respectively. For convenience in design and control of the machine, the coordinates (X ,Y, Z), figure 2.4, are used. And C is a constant dependent upon the dimensions of the products and the configuration of the rolling mill. Henceforth, the roll gap function can be simply expressed as G(x + C) in the analysis coordinate system (x, y, z). 3 Actually the author considers that it should be a good choice to make use of the available methods for bulge prediction with this model, amongst which there are several good methods. Besides, the author believes that accurate bulge prediction that is really helpful in production is more an art than a science, by its nature. x ( X ) Figure 2.4: The coordinates system Chapter 2. BASIC CONSIDERATIONS^ 13 Chapter 3 CONSTRUCTION OF VELOCITY FIELDS 3.1 The Dual-Functional Method As mentioned Chapter 2, the analysis will start with Yih's mathematical considera- tions[16], of which the related results are to be re-introduced here, however, detailed discussions should be referred to the original paper. Yih derived that for a continuous flow, the velocity vector can be expressed as the following: V = (grad x (grad g)^ (3.1) where f and g are two artificially chosen functionals that may represent stream surfaces in which streamlines are embedded, and along any streamline the two functionals assume the following form, f(x,y,z) = a^g(x, y, z) = b^ (3.2) where a and b are some constants. In addition, since the tool surface is the actual flow surface, f can be so chosen as to make a solid boundary an [surface, and g the free sides of the being rolled material and can be solved by the equation, V • (grad g) = 0^ (3.3) which is one of the properties of continuous flow.' lYih did not indicate this last approach explicitly, which, however, can be easily derived from the original idea.(See Section 3.2) 14 Chapter 3. CONSTRUCTION OF VELOCITY FIELDS^ 15 3.2 A Brief Derivation of the Solution to the Velocity Fields Incidentally, a much shorter and simpler mathematical derivation of the aforementioned idea was discovered during this work, which is the following. Assuming an incompressible material(nondiffusible), the following holds, ^ay. av^avz^ + Y^+ ^ = 0Oxay az (3.4) which can be written in the following form, div V = 0^ (3.5) Note that we have the following mathematical relationship, div (A x B) = B - (curl A) — A - (curl B)^(3.6) curl (grad 0) = 0^ (3.7) where A, B and 0 are any differentiable functionals. Suppose that there are two contin- uous functionals, f (x, y, z) and g(x,y, z), the gradients of which exist. Let A = grad f ^ B = grad g^ (3.8) Substituting Eq.3.8 into 3.6, and using Eq.3.7, the following results, div [(grad f) x (grad g)] = (grad g). [curl (grad f)] — (grad f) x [curl (grad g)] = 0^(3.9) Comparing Eq.3.9 with 3.5, we have V = (grad f) x (grad g)^ (3.10) Eq.3.10 suggests that V could be constructed by two functionals, and the only re- straint to these two functionals is that their gradients exist. In other words, f and g Chapter 3. CONSTRUCTION OF VELOCITY FIELDS^ 16 need to be constants. Henceforth, f and g can be artificially chosen. Naturally, f can be chosen as the tool surface, meanwhile, adopting the prescribed flow factor (see the deriva- tions below); and g can be so chosen as to represent the free surfaces of the material. The latter property leads to the following, V • n = 0 (3.11) where n is the normal vector on free surface of the material between the rolls. Since n and grad g have the same direction, Eq.3.11 can be replaced by V (grad g) = 0 (3.12) which is the same as Eq.3.3 and can be used to determine g. 3.3^The Derivation of the Velocity Fields Using the notation of gradient, V a f By^az ag^Bg By^az Eq.3.1 1+ may be written aLL az^ax ag^ag az^as j+ as (2L Ox^By ag^ag Os^By (3.13) where i, j, k are unit vectors along x, y and z directions, respectively. Thus, the velocity fields can be expressed as the following. (3.14) With the idea described in previous Section in mind, it can be easily understood that the assumption 5 in Section 2.2 is theoretically tenable. The remaining problem will be Chapter 3. CONSTRUCTION OF VELOCITY FIELDS^ 17 to choose the proper functionals, f and g, which should assume the simplest forms but satisfy Eq.3.3. As indicated earlier, the functional f can be chosen as the tool surface. From figure 2.4, under assumption 5 the roll surface can be taken as part of a cylinder, and expressed as z = t(x) = G(x +C)-F R — VR2 — ( 1 — x)2^ (3.15) Obviously, t(x) is a non-zero function in this case and Eq.3.15 therefore can be written as t(x) Since we are talking about flow, it will be seen that to relate Eq.3.15 to the roll peripheral velocity at this stage is a reasonable choice (also see next section), thus, Eq.3.16 becomes Qz t(x) = Q^ (3.17) where Q is an arbitrary constant related to the flow and its physical meaning will turn out later. Now define Q z f(x,y,z) =^ (3.18) t(x) which then meets the restraint that f should be a constant. Before using Eq.3.3 to determine functional g, some explanation relating to the as- sumption 5 should be given. As can be seen in figure 3.1, g(x,y, z) is assumed as the flow surfaces constituted by the free sides of the part of a bar being rolled and are parallel to z-direction. Note that the curvature of the sides is usually very large (the dotted line in figure 3.1) in reality and very close to a straight line (the solid line in figure 3.1), which may be confirmed by many photographs in the literature reflecting the deformation of the workpiece between the rolls that showed the projected profiles in the x-y plane remain 2=11 (3.16) Chapter 3. CONSTRUCTION OF VELOCITY FIELDS^ 18 g(x,y,z) Figure 3.1: Schemetic illustration of the top view of a rectangular bar between the rolls almost a straight line, and which, in fact, is not unexpected as a result of the requirement for keeping the deformation as uniformly as possible in the roll-pass design, and figure 3.2b and 3.2d, from [17], may serve as a good example amongst many. However, it seemed that researchers preferred to attempt an exact estimation of the curvature in the calculations with a complicated model. It is imaginable that the con- tribution of such estimation to the overall value of the energy dissipation would be very small, which leads to the idea that assuming a straight line as the side profile in the present problem is a good approximation. The main purpose of such consideration, of course, is to reduce as much as possible the difficulty in the mathematical handling and the computational demands. Now, the equation for a straight taper profile, figure 3.3, can be written as the fol- lowing, , ‘ wi — woy . wkx ) = 1 x + wo Suppose that g(x, y, z) will assume the form, (3.19) 4Y^g(x,y, z) =^= q w(x) (3.20) 171/ = tW2 v. = _-Qqtw Qqyw' (3.22 ) Chapter 3. CONSTRUCTION OF VELOCITY FIELDS^ 19 Skta *ad pkka viatqa at 1$ mitoy ay.amainens do-far:nod at 40% reductlea; St ^mm OV/h,.: :^ .^. Figure 3.2: Side views(a, c) and plan views(b, d) of 1S alloy specimens deformed at 40% reduction; height: 26 mm, width: 26 mm where q is some constant, therefore Eq.3.20 meets the constancy restraint of stream surface functional. Substituting Eq.3.14, Eq.3.18 and Eq.3.20 into Eq.3.3, the following results, Q Yw'e^Yw' 0(^ ^tw w2^w2^ (3.21) which indicates that q can be any constant. Substituting Eq.3.18 and Eq.3.20 into Eq.3.14, the following are obtained, Qqzti = wt2 Through observation, it can be decided that q = — 1 is the best choice, simply be- cause Vx can never be negative in the current coordinate. Hence, the velocity fields are Vs =-- Vz =-- Q tw QYWI VY = tW2 Q zt' wt2 (3.23) Chapter 3. CONSTRUCTION OF VELOCITY FIELDS^ 20 Figure 3.3: The schematic illustration of the simplified side profile of the section between the rolls determined as 3.4 Discussion of the General Solution At this point, it would be convenient to make a few comments about history. Hill[3] first proposed a similar velocity field solution for the flat rolling and bar drawing in 1963, in which researchers have shown great interests ever since. He must have found a particular way to reach his solution, but, unfortunately, he didn't indicate it2. Lahoti, et al, con- sidered Hill's work as a "unfinished" one, and continued the work in 1973[4]. However, 'For the sake of convenience, Hill's solution is quoted in the following. "The simplest class of approximating velocity fields with these properties (i.e., the similar assumptions as those of 1, 2 and 5 in Section 2.2, and the same coordinate system, (x,y,z), as in Section 2.3.) has Chapter 3. CONSTRUCTION OF VELOCITY FIELDS^ 21 as Hill himself didn't explain how to relate the solution to the rolling velocity(VR), nei- ther did Lahoti, et al. This, however, becomes significant when a solution to the rolling process, as general as it is in the present work, is sought. It can be noted from Eq.3.23 that for flat rolling, the following holds, Q vo = vrl..0 =^ (3.25) and this may be rearranged into, Q , Vo toW0^ (3.26) which indicates an interesting fact that the constant in the solution is the discharge. Accordingly, in the case of (a) a bar pushed in between two idle rolls, Q , Votowo^ (3.27) where Vo is a known constant. (b) a bar pulled through two idle rolls(Steckel rolling, see reference[1]), Q = vit(ow(i) ^ (3.28) where V1 is a known constant. (c) the conventional rolling process, Q = 0(xn, VR)^ (3.29) where 0 is some functional. Therefore, in the case of (a) and (b) , there is no need to worry about the estimation of the neutral points. This conclusion is different from the components(here, h =2.), 1^YID'^zh' (3.24)tly = h—, Vy = hw2) vz = Wh2 within the region bounded by the planes of entry and exit, the faces of the dies, and the lateral surfaces y = ±w(x), Such fields are clearly solenoidal and preserve the configuration; The immaterial constant of proportionality(the emphasis is added here) has been chosen so that the total rate of flow across each quarter-section is unity. " towo Chapter 3. CONSTRUCTION OF VELOCITY FIELDS^ 22 one by Oh and Kobayashi[5]. There are some velocity field solutions for the conventional rolling, say, [6], in which the relation, Q = Votowo, was used. This can be considered obscure, because in the conventional rolling, 170 is an unknown constant before starting the process, besides, I/0 may be temperature dependent. Thus, it is reasonable and instructive to adopt the prescribed factor Q in the solution, especially when considering the variation of the rolling geometry in the VGR process, where Q may no longer be a constant. In other words, Q may be an instantaneous constant (See next Section). 3.5 Further Considerations 3.5.1 Theoretical Aspect In the foregoing sections, the solution to the velocity fields has been treated in a general manner. It should have been noted that the derivation of Eq.3.23 is started from the original geometry without any restriction on the motion of the rolls. This is one of the nice things about the solution, as it would be very awkward and unnecessary to construct the component for the up-or-down velocity of rolls at this stage, which could be non-linear. The implication is that when Eq.3.15 is used, which includes a roll-gap function, the variations in the velocity components will naturally be taken into account. It should also be noted here that in the VGR process, owing to the ever-varying geometry, the velocity vectors in the deformation space become time-dependent, thus, it is a non-steady state process, which is the essential difference from the traditional rolling process, where the geometry is fixed and therefore it is a steady state process. In the present case, however, during the whole process, the geometrical resemblance is basically maintained, which may be the theoretical validity that one set of formulation can be used Chapter 3. CONSTRUCTION OF VELOCITY FIELDS^ 23 in the modelling of the present process. Accordingly, it is both reasonable and necessary to divide the process into a series of small pseudo-steady state deformation steps, that is, by focusing attention on each sampling instant, the geometry may be frozen, then a set of instantaneous solutions can be obtained. And by updating the geometry, a set of approximating numerical solutions may be finally produced, which may serve as a description of the dynamic characters of the process under study. This idea should also be considered as the basis for the discretization procedure in the design of the machine control system. 3.5.2 Technical Aspect In the VGR process, when the conventional rolling scheme is used, the discharge Q is a variable if the roll peripheral velocity is to be designed as a constant, which, apparently, is preferable; when Steckel rolling scheme is used, if Q were to be designed as a constant, the pulling mechanism would be very complicated. In any case, since Q is preferred to be a variable of the geometry, it is clear that the geometrical parameters have to be determined in advance in order to update Q, which will be a crucial issue in the control scheme. The roll-gap function G From machine design standpoint, it is important to consider the feasibility of a simple control scheme at this stage. In figure 3.4, the two coordinate systems, (x, y, z) and (X ,Y, Z), have the following Chapter 3. CONSTRUCTION OF VELOCITY FIELDS^ 24 Figure 3.4: The schematic illustration of a parabolic profile in rolling configuration simple relationship, X = C (3.30) Z^z Here, in order that a typical set of data may be produced, a parabolic type of product has been chosen, although the roll-gap is generically defined. In coordinates (X,Y,Z), the profile assumes the following form, (X — a)2 = —2p(Z — h) (3.31) which can be reduced to a2 (X — a)2 = b h(Z — h) where, (3.32) a2 = —2p(b — h)^ (3.33) Chapter 3. CONSTRUCTION OF VELOCITY FIELDS i.e., Let then Eq.3.32 becomes Z = e(X — a)2 + h Using Eq.3.30, Eq.3.36 has the following form in coordinates (x, y, z), z = e(C — a)2 + h For notational consistency, Eq.3.37 should be written as G(x + C) = e(C — a)2 + h a2 P = 2(b — h ) b — he = a2 25 (3.34) (3.35) (3.36) (3.37) (3.38) It can be seen that Eq.3.38 ensures a simple but reliable control scheme. However, detailed discussion is beyond the scope of present work and is left for later; for now, Eq.3.38 is adequate for the present mathematical analysis. The length of plastic zone Note that in present problem, the length of the projected arc in contact is also a variable — a function of the roll-gap and to, i.e., /(x). From figure 2.3, 1(x) = VR2 — KR— to) — GI2^ (3.39) With the above geometrical relationships determined, a new velocity field solution is established, which will be the basis for the next chapter. Chapter 4 FORMULATION THROUGH UPPER BOUND APPROACH The upper-bound theorem states that among all kinematically admissible strain rate fields, the actual one minimizes the following expression, using the same notation in [24], p.419, J* = a,*343dV + ElD k^idyl; _^TitcisT^(4.1)Sr where, ê j is the strain-rate field derivable from v:;^is a kinematically admissible veloc- ity field; lAV*1 is the amount of velocity discontinuity along the surfaces of discontinuity Si); and o-:3 satisfies the yield criterion and is associated with i:i. In present case, only the conventional rolling scheme is considered, therefore, the back and front tensions are zero, i.e., T, = 0. SD is the surface where velocity discontinuities occur, and consists of three parts in present problem: the friction surfaces, S f , between the rolls and the material; the entrance and exit planes, Sdo Sdi, where shearing occurs; accordingly, SD = Sf Sdo Sdi) J* =u,dV +^TolAV sjelS f +^klAVdIdSd3 V Jsf Or using the conventional notations, Eq.4.2 may be written as, = Ed +Ef +E8 where, Ed = J èdV 26 Eq.4.1 becomes (4.2) (4.3) (4.4) Chapter 4. FORMULATION THROUGH UPPER BOUND APPROACH^27 =mkIVsIdSf SI =kIAV didSd 19.10+Sdl Definitions of various terms in Eqn.4.4 through 4.6 is given in the following sections. 4.1 The Strain Rates The six strain tensors are defined as the following, v.^a x1 /ay,/ ay., E. = ^ax isv — 2 ax + ay ^avy^.^1 (avy ay.,Ey= ^ evz = --‘ az + ay 'ay ^ ay.^..^1 , ay. ay.,ez =^Ezx — k^ + —)^az 2 az^ax (4.7) Substituting Eq.3.23 into Eq.4.7 and re-arranging it, the following results, i t^w' ex = —Q —(--- + —)tw t = Q wtw2 (4.8) (4.5) (4.6) Q ywi^211/ 21w2 t^) w Q z^2t' 212 [tEIZX iyz = 2 [(u1)2 1_ (t')2 Hh w)^tw= (t' 2w)M =^t ^w) IT t"^t' 12t'^w'\ T^-t ) Substitution of Eq.4.10 into Eq.4.9 results the following, Ed = 40.fift,weQ Nip 4. 1NI y 2 + N2z2dydzdx o o o^3 tw After integration against y and z, Eq.4.12 becomes Ed = 2 .1 —Qa Z(x)dx 3^o where, (4.11) (4.12) (4.13) Chapter 4. FORMULATION THROUGH UPPER BOUND APPROACH^28 4.2 The Energy Rates 4.2.1 The energy dissipation rate due to deformation(Ed) Using the Levy-Mises equations, Eq.4.4 may be written as, •^2Ed = —0 \11 • •J 2 EljEljdV (4.9) Since, 1 . 2—E .^.E2 E2 _4_ 2 2 y^zj^; zs -r- c yz 1 Q2 2^2 2 2 t2w2 " M^N2z2) where, (4.10) Z(x) = ^ 12 + m2w2 Ni^+ m2w2 + N212 V12 M2y2 N2z2 ^ inNt^v12 m2w2 212MNtw [1,1(02 m2w2, Mw, Nt) —^, 0 , NO] -F Nm 2 tw (02 + m2w2, Mw, Nt) — (/, 0, NO] (4.14) Chapter 4. FORMULATION THROUGH UPPER BOUND APPROACH^29 where, foc Pi(a, b, c)^ln(b + Va2 + x2)dx = cln(b + Va2+ c2) -c-f-bln (Va2 + C2 + a ^C) (4.15) + 2 va 2 - b2 tan' iVa2 - b2 Wa2 + c2 - al L a + b^c 10. P2(a, b, c)^x2 ln(b + Va2 + x2)dx 1 ^1= d {c3 ln(b + Va2 + c2) - c3 -I- c(a2 - b2) 1 1^(Va2-F-bcVa2 + c2 - -b(3a2 - 2b2)1n ' ^+ c2 _4_ c) 2^2 a (4.16) - 2(a2 - b2) tan-1^a (Va2 - 112 Wa2 + c2 - all -I - b^c 4.2.2 The energy dissipation rate due to friction between roll surfaces and workpiece(Éf) From figure 4.1, the velocity discontinuity between the surfaces of roll and material can be expressed as the following, IIEf = 4mkQJ fw o Jo \ ,\2 [(YW ) + (VR 1 wt /1 + t'2) 21 (1 + t12)dxdy AVs =- Vy + VR - (V. + V.) (4.17) and the magnitude, 1,6Ars I = VVy2 + ( VR — 1/1/2 + 1/2 )2 QN1 (yw i 2)^(VR V1 + t' 2 \ 2 ^tI.V2 )^I\ Q^wt ) Note that on the surface, z = t(s). Then Eq.4.5 becomes, (4.18) (4.19) - (V1 + 111 ) --i - - - V, + Vy + V,^AV x V, + Vz ,/^ I 1 1 Chapter 4. FORMULATION THROUGH UPPER BOUND APPROACH^30 Figure 4.1: The schematic illustration of the velocity vector space in conventional rolling process and after integration against y and z, Eq.4.19 becomes Ef -= 2mkQ 1 1 F(x)0 + ts2dx o where, w (wi ) 2 ( VR Wt^VRV1 + 0 ) 2 V]. + ti2 ) 2 .F(x) = — 2 1 \ tto) + Q + ws^Q ( wt^) wt^) (4.20) in [{ Wi^ ()2 (V^R ^N/1 + ti2)2 . VR^■,/1. + t'2 Wi +^iW + Q^wt^• Q^wt^ (4.21) 4.2.3 The energy dissipation rate due to velocity discontinuities (E',) From figure 4.1, at the entrance plane, AVdo = V — V, = V1, + V,^ (4.22) Chapter 4. FORMULATION THROUGH UPPER BOUND APPROACH^31 and the magnitude, 1AVdol = VVy2^ (4.23) And at the exit plane, since the width of the deforming body is almost the same as the width of the rigid body, we may assume Vy = 0; thus, from figure 4.1, the following holds, AVdi = V — Vz = V, = 0^ (4.24) Substitution of Eq.4.23 and Eq.4.24 into Eq.4.6 gives, note that, Vo = Q/woto, = 4k j V(Vy2 1/2).=ods, = 4kVo ft.() ft. Jo^\ (ytv, ) 2 (Zii ) 2 dzdy w x=o^t (4.25) The integration of Eq.4.25 results in the following, kQS(x,wi)^ (4.26) where, t12^+ 012 + W12 S(X, WI)^[0'2 +^ln 4^F I^I^2w'2^1 I tlw7P2(t t^) z=o (4.27) 4.3 The Minimization Procedure Since the unknown parameters in the velocity fields solution are w1 and xn, the problem is reduced to solving the following non-linear equations, awl^'^ax„ ^ (4.28) Chapter 4. FORMULATION THROUGH UPPER BOUND APPROACH^32 By substituting Eq.4.3 into Eq.4.28, the equations to be minimized are obtained, atd atf at,^, +^ + ,^ =o (4.29) owl mill owl atd DE1 at, ^ + ^ + ^ = 0 (4.30) axn axn axn 4.3.1 The derivatives of the geometry From Eq.3.19, the following is obtained, wi = (w1 — w0)-11, along with aw x ^, awl^1 and Ow' _ 1 awl — / through which the following results readily, (aaww:).=0 = 1 and Eq.4.31 also gives the following, 1,^1^\ wo --= kwi — wo )—1 7 (4.31) (4.32) (4.33) (4.34) (4.35) 4.3.2 The neutral point Apparently, in present case, the neutral point is ever-changing. However, the following still holds instantaneously for each frozen moment, figure 4.1, YR = Yx(Xn) + Yz(Xn)^ (4.36) V1 + ti(x,) w(x)t(x) cl ( x n ) = vR (4.39) Chapter 4. FORMULATION THROUGH UPPER BOUND APPROACH^33 Using Eq.3.23 and note that, z = t(x„), we have the following, 4? /VR =-- - V 1 + ii Wi X=Xn (4.37) or it may be written as, Or, VR V1 + e(x) Q^w(x)t(x) (4.38) 4.3.3 Solutions The analytical solution to Eqn.4.28 is too complicated and may be found in Appendix A. The numerical solution is obtained by coding in MATLAB. A CPU time of 15 seconds for each cycle of calculation is achievable on Sparc-I SUN workstation. The results are given in Chapter 6. Chapter 5 FLOW STRESS MODEL In the foregoing chapters, a fundamental geometric description of the process has been given. To achieve a reasonable control of the roll-gap in the VGR process, a detailed and accurate roll force change prediction is required. It should also be noted that, in order that the newest developments in material science may be adopted in the calcula- tion, the assumption that the materials are rigid-perfectly plastic is in need of further modification. In other words, it is necessary to establish a correlation describing the changes in flow stress as the workpiece moves through the rolls and between the passes, so the changes in rolling forces may be predicted. An adequate mathematical model of flow stress incorporating the thermal characteristics of materials would be of great importance in implementing the VGR process and in controlling the final properties of the rolled products. 5.1 Dynamic Flow Stress Unlike cold working, the hot deformation of steels at elevated temperatures (above the recrystalization temperature) may be divided into three transient stages[2511261[27] — 1. Microstrain deformation. This is the interval during which the plastic strain rate increases from zero to the approximate forming rate, and the stress in the workpiece rises rapidly, Fig. 5.1. The loading slopes are rate controlling in this interval. 34 True stress Work hardening] Dynamic recovery High strain rate Dynamic recrystallization Low strain rate Chapter 5. FLOW STRESS MODEL^ 35 _ 0 0.2 0.4 0.6 True strain Figure 5.1: Schematic flow curve for metals at elevated temperature. From point 0 to 1, the plastic strain rates increase rapidly. Note that the scale for the strain is arbitrarily chosen from typical C-Mn steels[27] to illustrate a possible range in applications 2. Plastic yielding. The region of the yield stress marks the beginning of work hard- ening. The slope is temperature and rate sensitive, but in general, is gradually decreasing. 3. Dynamic restoration. This stage is most complicated due to the concurrent opera- tions of the strain hardening and softening. Generally speaking, if the strain rate is low and the hardening effect is only offset by the dynamic recovery process alone, a steady state will quickly be reached, so a simple flow curve results. The alternative restoration process is dynamic recrystallization, which modifies the flow curve in the following ways. At high strain rates, the flow stress rises to a maximum at peak strain(point 2 in figure 5.1), and then it diminishes to a value intermediate between the yield stress and the peak stress(See figures in Appendix Chapter 5. FLOW STRESS MODEL^ 36 B). At low strain rates, the dynamic recrystallization process causes a softening followed by a oscillatory type of flow curve before a steady state is reached, figure 5.1. In [27], I. Tamura, et al, stated the following. For most commercial steels, dynamic recovery behaviour is observed in iron and carbon steel deformed in the ferritic region, or in ferritic alloy steels such as silicon iron. Dynamic recrystallization takes place in the hot-deformation of -y including HSLA steels and austenitic alloy steels. This may serve as a general guideline, although care should be taken for each different material. The reduction ratio in the VGR process can be anything between 8% and 50%, which corresponds to the true strain range from 0.051 to 0.693. From figure 5.1, it is clear that the flow stress will vary significantly during each pass. In other words, in reality, the materials are rate dependent within the strain range of hot rolling. Accordingly, instead of using a flow stress value corresponding to the steady state, which itself may result in conflicts, it is necessary to have a flow stress model to count in such variations. 5.2 Incorporation of Rate Dependency It has long been concluded that the flow stress model that gives a complete description of the hot deformation of steels in the roll-bite assumes the form[28]: = f(ccrmposition of the alloying elements,T,e,e)^(5.1) where T is temperature. As is stated in Section 3.5, the VGR process is a nonsteady state forming process, and in the last Section we concluded that it is not reasonable to use a constant yield stress Chapter 5. FLOW STRESS MODEL^ 37 throughout the whole process. A yield function similar to the above form would be very important in obtaining a solution close to reality. Therefore, modifications are needed since, theoretically, the upper bound approach is only valid in rigid perfectly plastic materials. It is clear that, unless the strain is beyond the steady state point, which is not the case in rolling process, the state of rigid-perfectly plastic deformation will never be reached. In order that rate dependent' effects can be incorporated, the concept of instantaneous yield conditions has to be adopted and further assumptions must be made. First of all, it is interesting to note that the available flow stress data are usually obtained from three sources: (1) uniaxial compression tests; (2) thin-walled tortion tests; (3) plain strain compression tests. A close examine of those tests tells one that it is impossible to get rid of the rate effects. Therefore, the results from these tests are usually similar to those in figure 5.1. Figures B.1—B.4 in Appendix B give some results from experiments that clearly show irregularities in flow stress due to the restoration process. This indicates that even though a perfectly plastic material is assumed, flow stress value at the steady state can not be used. Instead, one has to use the data at certain strain points. The implication is that the solution calculated from the available experimental yield stress data can only be an approximation. Nevertheless, this is the usual practice in modelling of the conventional rolling process. Then, the usual assumption that the material is perfectly plastic may be modified, that is, the appropriate assumption in the real situation may be that the material has a isotropic rate dependency'. Thus, the material is still a Mises material, and the yield locii become a set of uniform expansions of the previous one, i.e., a set of circular cylinders. Furthermore, a similar assumption adopted by Hi11[17], when an extrusion process was considered, should also be made. That is, the mean strain of the deforming body remains 1The definition for rate dependency may be given differently by different authors. In this work, the term, "rate", should include the temperature decreasing rate, the rate of restoration and the strain rate. 2The condition of isothermal is implied. Normal platen pressure Specimen strip Platen Chapter 5. FLOW STRESS MODEL^ 38 Figure 5.2: Schematic setting of plane compression test unchanged regardless of the rate dependent characteristics of the workpiece material. This implies that the mean strain of such a rate dependent material is the same as that of a perfectly plastic material, since isotropy and volume constancy are assumed. Based upon this assumption, the concept of mean equivalent strain, or effective strain, E, and of mean strain rate can be introduced. 5.3 Mean Values of Strain and Strain Rate In Eqn.5.1, it should be clear that for a practical purpose, a, e, é, and T can only be given in mean values. 5.3.1 Mean Strain To explain the concept of mean equivalent strain, it is helpful to look at the approach applied in the plane compression test. In figure 5.2, the lateral strain may be taken as zero, i.e., ez = 0, and volume constancy is assumed. The thickness strain ez is Chapter 5. FLOW STRESS MODEL^ 39 Figure 5.3: Schematic relation between mean yield stress and mean equivalent strain in plane compression test numerically equal to longitudinal strain ey. These results may then be used in the Levy- Mises relationship, 1 ^ E = DI2[(e s — Eyr + (Ey - Ezy + (Ez - E.)21 2 vaez which indicates that the mean equivalent strain is proportional to the strain in the direction of reduction. So the mean flow stress am in such tests is usually given as = F1(e) = F(e)^ (5.3) where F is determined by the test. The results would be similar to those in figure 5.3. The above approach reveals an important idea. For an isotropic rate dependent material, the mean yield stress may be obtained by establishing first an appropriate procedure to calculate mean equivalent strain. This idea was successfully used by many researchers in the past, such as, Kudo[29][30], when upper bound approach was applied. Johnson[31] and Dodeja and Johnson[32] used it in an analysis of forging and extrusion process in an upper bound solution, where the E was calculated through the estimation of the total power consumed while assuming the workpiece is perfectly plastic. = (5.2) Chapter 5. FLOW STRESS MODEL^ 40 For the present case, since no bulge is assumed, the following method described in [33][34] can be used in estimating the mean strain, E^j1 VA.^A. 1 0 ^2 2^r2 dx^ (5.4) where the components, ex, e ^ez are taken as the principle strains. The assumptions in Section 2.2 imply proportional straining, i.e., EX + Ey Ez = 0^ (5.5) which leads to — 1 I ‘IE2 E E^e2 dx1 0^y^Y z z where the instantaneous strains are usually defined as, ez^ln t(x) = ln w(x) to^wo Note that sometimes the stress-strain tests were carried out in uniaxial compression. Then, the following conversion should be done before using the results from Eqn.5.4, eu ^2 (5.8) It is noticed in the literature that the stress testing data for rolling process were tradi- tionally arranged in a correlation between stress and reduction ratio, which is convenient in application. In such cases, a mean reduction ratio may be calculated first as follows, r„, 1 — ee^Or rmu 1 — elu ^ (5.9) (5.6) (5.7) 5.3.2 Mean Strain Rate In the present work, a similar hypothesis in calculating the mean strain described by Chakrabarty[35] will be used, which is the following. Chapter 5. FLOW STRESS MODEL^ 41 Considering the second invariant of the plastic strain increment tensor, an equivalent plastic strain increment is defined as 2 = 3 –i (5.10) Therefore, =^dV—1 I In V V (5.11) By substituting Eqn.5.11 into Eqn.5.3, a rate dependent hypothesis, which is exactly the same as the strain-hardening hypothesis described in [35], may be introduced as the following, crm=F(ëm) (5.12) where F should also be considered as a function of the other two variables(the mean strain and the mean temperature). Accordingly, the strain rate increments may be calculated by firstly assuming that the material is a perfectly plastic, so Eqn.4.3 can be used to obtain a reasonable mean strain. 5.4 Temperature Prediction Although it is known that temperature has a dominant effect on the magnitude of the flow stress, not many models incorporating temperature changes have been formulated for the purpose of predicting the flow stress[36]. A large amount of research work on the thermal effects has been done in the past two decades. An excellent review may be found in reference [37]. 5.4.1 Basic Considerations Since the VCR process is designed for the rolling of rather small workpieces, the tem- perature of such small workpiece will drop very quickly. Accordingly, drastic changes in flow stress during rolling will be expected. Chapter 5. FLOW STRESS MODEL^ 42 The important variables that have great influence on temperature changes are: • Processing time per pass. • Interpass time. • Temperature of the rolls (determined by cooling methods). • Conditions of the roll-workpiece interface. It is clear that the VGR process planning will largely depend on the above variables. A valid correlation is needed. The present work is based on the idea that a combination of separate submodels, such as, deformation, flow stress and force, etc., should be used, so that new developments can be easily adopted. Fundamental analysis is preferred, since all the empirical formulae were meant for large mill operations and are difficult to make necessary changes. All the available fundamental analyses, e.g., [37][38][39], are similar in that a heat transmission from the furnace to the finish stands was considered. In the present case, a relatively small machine is considered, and the temperature will drop due to contacting the rolls and the run-out table, convection to the spray and the air, and radiation. It is conceivable that in such a heat transfer chain, the temperature distribution of the rolls is affected by too many factors, such as the cooling method, the size and structure, the surface condition and the number of rolls, etc. It is almost impossible to give a accurate prediction due to the variation of the interface area, even if one has a detailed specifications of the process. Therefore, the present author believes that, when the VCR process is considered, the temperature of the rolls has to be obtained from some sensor reading. Chapter 5. FLOW STRESS MODEL^ 43 5.4.2 Heat Conduction Model The heat conduction in a. isotropic solid, the rectangular bar being rolled, may be mod- elled by the following equation: a2T .92T 52T). 8^ay2^2 aT k (— +^+^+ q = Pweiwax2 az (5.13) The coordinates are the same as before. The present work will only consider one dimensional transient heat transfer, and necessary simplifications to Eqn.5.13 have to be made. 5.4.3 Assumptions 1. In order to simplify the calculation in the temperature model, it is assumed that the workpiece is always a rectangular bar when convection and radiation are considered, that is to say that the influence on heat transfer due to the changes in the geometry is ignored. Mean values of the dimensions will be used in calculation. 2. Heat conduction along the length of the product is negligible in comparison with the heat transfer by the movement of the product. 3. The temperature gradients in the roll-gap along the width direction is negligible since no contact is made between the sides and the rolls. In the roll-bite, only one-dimensional transfer will be considered. It is assumed that the small amount' of heat generated by the deformation, and of heat generated by friction between the rolls and the workpiece will somewhat offset the heat loss from the two open sides in the roll-bite. 31t can be easily found out that the heat generated by deformation in the roll-bite is very small. Some researchers, say, [41], calculated the amount of temperature raised due to deformation in plate and strip rolling and the result is only about 20 °C. The contribution of such a variation in temperature to the change in flow stress is negligible. Chapter 5. FLOW STRESS MODEL^ 44 Then, it may be assumed that 4 = 0. Thus, Eqn.5.13 is reduced to: 02T^OTk.^ pwCw az2^Or (5.14) 4. The temperature gradient at the centerplane through the thickness is assumed as zero. Previous work[40] showed that no temperature difference exists between the top and bottom surfaces of the workpiece. This indicates that no heat flows across the centerplane. That is, aT(0,T) Oz — (5.15) 5. The effect of the scale layer on the surfaces of the workpiece can be ignored. This is a usual practice in the research works, due to the lacking of experimental data. 6. It is assumed that the stages of cooling under consideration are: radiation and convection before entry, cooling by contacts to the rolls, radiation and convection after exit. The effect of contacting to the entry and exit tables is ignored, because no suitable data can be applied. Besides, when the workpiece is out of the roll- bite, the body temperature gradient is ignored, only surface temperature will be considered. 7. Thermal conductivities, k., and specific heat, G, of steels with various carbon con- tents as function of temperature was established by Devadas[371 by fitting available data from BISRA Tables[42]. The results are shown in figures 5.4 and 5.5. These two figures indicate that, with the temperature increasing to about 1200 °C, the effect of varying the carbon contents on k, and CI, decreases. In cases of 0.4 C and 0.8 C, the differences are negligible. Carbon Content • 0.08%C + 0.4%C * 0.8%C 0.1 Chapter 5. FLOW STRESS MODEL^ 45 31 30 29 2: 26 25 24 E5 PE 22 23 750 850 950 1050 1150 TEMPERATURE (°C) Figure 5.4: Thermal conductivity as a function of temperature for 3 carbon content 0.08 C, 0.4 C, 0.8 C in austenite phase. (From Devadas, 1989) Figure 5.6 and 5.7 show the results from fitting some data in [43] for two different C-Mn steels. When choosing the appropriate functions for thermal conductivity and specific heat, reference may be made to these results: = —1.11024(10-4)T2 + 0.3670T + 399.5600 (for 0.19C 1.39Mn steel) = 0.2511T + 411.4217^(for 0.37C 1.44Mn steel) (5.16) (5.17) = —0.0001T2 + 0.4103T + 390.9195 (for 0.23C 1.51Mn steel) (5.18) k, --- —0.0246T + 49.202^(for all C > 0.15 C — Mn steels)^(5.19) Chapter 5. FLOW STRESS MODEL^ 46 690 680 - 670 - 660 650 - 640 630 - 620 610 600 - 590 - 580 ^ 600 553.867+0.164T-5X 10E-5T21 1652.886-0.1927T+17.3X10E-5T21 1431.25+0.247T+2.2X10E-5T2 800 ^ 1000 Carbon Content • 0.08% C + 0.4%C * 0.8%C 1200 TEMPERATURE CO Figure 5.5: Specific heat as a function of temperature for 3 carbon content 0.08 C, 0.4 C, 0.8 C in austenite phase. (From Devadas, 1989) 5.4.4 Initial and Boundary Conditions The initial temperature of the workpiece is taken as that when it is moved out of the furnace. i.e., T(z, 0) = To^ (5.20) In addition to Eqn.5.15, Eqn.5.14 is subjected to the following conditions, which can be derived from heat balance. Owing to symmetric cooling, at the interfaces between roll and workpiece, only con- sider convection, we have the following, aT aZ^ = 11(T)[T(i, T) Tr]^(5.21) z=t(x) where 1/(r) is the heat transfer coefficient, to which further definition is given in Section 5.4.6. Experimental results: o --- 0.23 C 1.51Mn steel + --0.37 C 1.44 Mn steel Chapter 5. FLOW STRESS MODEL^ 47 Temperature -- Celsius degree Figure 5.6: Thermal conductivity as a function of temperature for C-Mn steels 5.4.5 Solution to the Heat Conduction Equation In terms of computing speed, it is preferred to have an analytical solution to Eqn.5.15. The following is a brief account of the method used for finding the solution. A detailed description may be found in [51]. 1. Assuming the solution has two parts — one is steady and the other is transient, T(z,r) = U(z)+ v(z,r)^ (5.22) then, using the boundary conditions, the steady part, U(z), may be determined as U() Tr^ (5.23) 2. Setting v(z,r) = T(z,r)— U(z), using Eqn.5.14 and the boundary conditions, the following results, av(0, r )^0 az Chapter 5. FLOW STRESS MODEL^ 48 Temperature -- Celsius degree Figure 5.7: Specific heat as a function of temperature for C-Mn steels av(t, r) ^ -I- hkv(t, r) = 0az v(z,0) = To — U(z) where, hk = H(r) K = k, ks^p.,,,Cu, 3. Eqn.5.14 with homogeneous conditions of Eqn.5.24 has a solution as (5.24) (5. 25) v(z,r) = f(z)e —AKI- ^ (5 .26) where, K = ka pwC., and A is an eigenvalue and f(z) is an eigenfunction of Sturn-Liouville problem: (5.27 ) 02f az2 + Af = 0, with 494°)^ _ 0 az — - (5.28) v(z, 7-) = E BnCOS(ZIZI)e- n=1^ AnKr 00 (5.31) Chapter 5. FLOW STRESS MODEL^ 49 The general solution to Eqn.5.29 is f(z) = Asin(z\a) Bcos(zVX). Using the boundary condition, one may obtain Nasin(t VX) = hkcos(tV)'), then we have f(z) Bcos(zA^ (5.29) VAT, = hkcot(tX) (n = 1, 2, 3 ...)^(5.30) From Eqn.5.26, Eqn.5.29, Eqn.5.30, 4. When T = 0, from Eqn.5.25, To — U(z) E BnCOS(ZVAn)^ (5.32) n=1 Using orthogonality, Jot [To — U(z)]cos(z,VAn)dz = Bn cos2(zOn)dz^(5.33) Bn may be decided as, ^ ^ 2(T0 — Tr)sin(t0,-,) n + So, finally the desired solution is v(z,T) = E 2(To — Tr) i sin(t t, ) cos(z\IZ.,)e-AnKr ^ n t^+ s n(2VAT) (5.34) (5.35) 5.4.6 Heat Transfer Coefficient at Workpiece-Roll Interface Recent experiments[46] showed that the heat transfer coefficient at the workpiece-roll interface is a variable. In [46], when summarizing the experiments, it is stated that, The heat transfer coefficient is seen to vary from 17 to 57 kW/m2K during the first 0.015 second in the roll gap and then to remain constant thereafter. Chapter 5. FLOW STRESS MODEL^ 50 65^ 60 55 50 45 '4- 4° 35 30 25 20 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02 Time in seconds Figure 5.8: Example for heat transfer coefficient at workpiece-roll interface (lubricant: water) The progressive increase in the heat-transfer coefficient in the first 0.015 sec- ond is thought to be due to the decrease in contact resistance with the increase in specific rolling pressure, which reaches a maximum at the neutral point. This may be the reason that the coefficient is expressed as a function of time in a more recent work[471 as the following, H(T) = 19.561 + 42.47(102)-r — 1.09(105)7-2 (kW/m2K) (5.36) where T is in seconds. Figure 5.8 shows the plot. It is clear that the transient period is very short. Therefore, in present study the coefficient may be taken as constant corresponding to the mean steady state value. It should be noted that lubricants have significant effects on the magnitude of the coefficient. Using the results from compression tests in studying the contacts between high temperature (780 °C) specimens and low temperature specimens, Murata and his 15 0 Chapter 5. FLOW STRESS MODEL^ 51 Table 5.1: Heat-transfer coefficients at roll-gap interface for different lubricants Lubrication condition Heat-transfer coefficient, No scale kWm-2K-1 Scale No lubricant 29.1 - 34.9 7 - 10.6 Water 23.3 - 81.4 10.6 Hot-rolling oil 200 - 460 5.8 Hot-rolling oil+20% CaCO3 69.8 - 175 12.8 - 23.3 Hot-rolling oil+40% CaCO3 12.79 - 17.4 ... KPO3 5.8 ... coworkers[49] calculated the heat transfer coefficients at the interface under different lubrication conditions, and suggested some useful values for hot rolling processes. The results are reproduced in Table 5.1. Devadas et al[36} used 37 kWm-2K-1in their work. They stated that according to Stevens et al's estimation[50], the heat-transfer coefficient at the roll-gap interface is 37.6 kWm-2K-1during the first 30ms and 18 kWm-2K-1 thereafter. They examined the coefficient in the range 30-50 kWm-2K-1, and concluded the following. • A 25% variation in the coefficient causes a 2% change in surface temperature. • 40% CaCO3 in a hot-rolling oil reduces the coefficient to 13 kWm-2K-1from an average value of 50 kWm-2K-1for water lubrication. The difference of the surface temperature at the exit for the two cases is 175 °C. These results should be taken into account in the VCR process design and the produc- tion. The present study uses water as coolant and an average value of 50 kWm-2K-1as the coefficient, so the roll surface temperature may be taken as around 100 °C. After integration and some tedious manipulation, the following results, 11-iT = {th {e2hcs,T' (1^hc 713 EB^ 1 Chapter 5. FLOW STRESS MODEL^ 52 5.4.7 Convection and Radiation When the workpiece is moved out of furnace, temperature starts dropping because of convection and radiation. The method for estimating heat loss by convection and ra- diation used in this work is similar to that suggested in [44][45], which consists of two steps. 1. Establishing a model for one-dimensional flow According to assumption 6 in Section 5.4.3, when taking its two pairs of the two parallel surfaces as infinitely large, an instantaneous heat balance of the bar in open air holds for each pair: Rate of heat loss from surfaces = —rate of heat accumulating in the bar Thus, we have the following equation, 2ha.A.(T. — Ta) 26.13A3(T: — 11) = — pCialiddTs^(5.37) where Ta is the ambient temperature which, in comparison with the workpiece tempera- ture, may be taken as zero. Assuming Sb = pC — V . = pCL A where L will be substituted by either t or w. Eqn.5.37 becomes, 2(haT. EBT) = S- dr by rearranging, Tf^dT. —= 7, dr IT, 2h,T, E.1371,4 (5.38) (5.39) (5.40) (5.41) Chapter 5. FLOW STRESS MODEL^ 53 where, AT is the time delay for temperature dropping to Ti; h, is convection heat transfer coefficient and is chosen to use forced convective heat transfer correlations for flow over a plate. Since the rolling velocity is not high and the cooling can be treated as laminar flow over isothermal plates, the average heat transfer coefficient is given by[48]: h, = 0.664R.0.5 po.33kc r Lc (5.42) where Re and Pr are the Reynolds and Prandtl numbers that are defined as usual; L, is the characteristic length, determined as the length of the plate in contact with the air flow; Ice is the thermal conductivity of the air. 2. Establishing a model for two-dimensional flow Using the idea introduced in [45], due to J. S. Langston, the following may be written, (Tf)= (Tf)^Tf \^total^width^) thickness (5.43) Substituting to and wo or ti and w1 into Eqn.5.38, respectively, corresponding values of Sb may be obtained. From values of Sb, one will be able to compute (Tf)wzdth and (Tf)thidknese before the entry and after the exit. Thereafter, Eqn.5.43 may be solved under different conditions. 5.4.8 Mean Temperature of the Deforming Body in the Roll-Bite With all the above relationships established, it is very convenient to compute the desired mean temperature for each pass. Eqn.5.35 gives the temperature distribution along the z-direction at certain instant. By chosing a finite number of points equally spaced along z-axis, i.e., ;, and time-axis — same as the x-axis —, i.e., T3 , the mean temperature may be computed in the following Chapter 5. FLOW STRESS MODEL^ 54 manner, 1 1 ni T„, =^E E T(zi,r;)^ (5.44) ni ni J=1 i=1 which is similar to finite element analysis. In this study, it is considered that accurate enough solution should be obtained by using i = 10 and j = 20. 5.5 Flow Stress Models 5.5.1 Model for Steels The current analysis will consider low to medium operation speed with a strain rate range of 1-50s-', while the temperature may vary in the commercial hot rolling range, i.e., 800-1250 °C. The flow stress model used should apply to this conditions. It can be easily found out in the literature that numerous researchers utilized the unified creep relationship that works well in hot rolling operations, which is the following, or rearrange it, = A sinh(aa) exp^Qa R T (5.45) Z = i exp Q.R ) = A sinh(ao-)^ (5.46) gT where A, n and a are constants independent of temperature, Z is the well-known Zener- Holloman temperature compensated strain rate parameter. A reasonable value of the activation energy, Qa, for a wide range of C-Mn steels(0.05- 0.68 C, 0.44-1.64 Mn), has been reported[52] to be approximately 312 kJ/mol, although it was observed that some variations did exist, such as 270 and 286 kJ/mol cited in [37]. Devadas[37] showed that both A and Q. are variables of strain, though the variations are small. Chapter 5. FLOW STRESS MODEL^ 55 Eqn.5.46 may be linearized as, Qa lni = lnA nln[sinh(aa)] RgT (5.47) The empirical constants may be determined by holding either the strain rate or the temperature as a constant and using the experimental data(further explanation on the procedures to determine the constants is given in Section 6.1.1). For C-Mn steels, the data published in [53] are used in present study. In order to incorporate the data in the computer algorithm, a fairly large amount of data is read from [53], and then fitted by least square approximation into some 80 different polynomials, each of which is within 5th to 7th order and is in a matrix form. Figures B.1—B.4 in Appendix B show the fitting results. Such a wide range of data is needed to suit the variation in the C and Mn composition, and the variation of the temperature, strain and strain rates. Such a data-base is also convenient for interpolation or extrapolation in determining Eqn.5.45. 5.5.2 Using Lead as a Modelling Material It is undisputed that an accurate flow stress model is absolutely necessary for a desirable solution to the hot rolling problem. Unfortunately however, since it is generally not feasible to obtain such information from production mill trials, a complete set of data on the deformations of various metals under the processing conditions does not exist. "It was observed that lead strain hardens in a similar manner to steel at elevated temperatures, and that the static softening behavior can be modelled using equations similar to these for C-Mn and stainless steels." [54], Hence the laboratory techniques of simulation and modelling have been developed by using commercially pure lead [54][55]. It should be noted here that lead can be only used to simulate the dimensional changes during hot rolling, but not the temperature effects on the rolling loads. However, using Chapter 5. FLOW STRESS MODEL^ 56 i 30 28 26 z 24 22 . cg-0 Ti 20 ''-, 18 16 14 01 ^ 0.15^0.2^0.25^0.3^0.35^0.4 ^ 0.45 ^ 0.5 ^ 0.55 Reduction ratio --- (Ho - h) / Ho Figure 5.9: Flow stress model: the correlation between yield-stress and reduction ratio of pure lead. Data source: [56] lead as a modelling material may be a good approach to reduce the experimental errors to minimum caused by the thermal effects and the difficulties in handling when operating at high temperature. The data used in [55] came from [56], which was the only source that can be found at this time providing a wide range of compression testing data. By using least square data fitting technique, the following stress model is obtained, ^a = 742.5r4 — 1028.9r3 -I- 485.2r2 — 55.8r + 17.2^(5.48) and figure 5.9 shows the fitting results. In Eqn.5.48, instead of using strain as variable, the reduction ratio, r, is used, which was the testing results before conversion. Chapter 6 MODEL VERIFICATION In the above chapters, a complete model for implementing the VGR process has been fully described, which consists of the deformation model, the flow stress model and the temperature model. It is important that these models be verified in experiments before the machine design process may be carried out. Due to the lacking of the experimental equipments, a complete set of testing data, regarding a certain workpiece composition, temperature changes, strain-stress correlations, roll force and torque, etc., is not available. Nevertheless, the model may be verified separately, by using some appropriate data reported in the literatures. 6.1 Flow Stress Model Verification 6.1.1 Empirical Constants in Flow Stress Model The linearized flow stress model, Eqn.5.47, is used to determine the empirical constants, A, n and a. Due to the availability of experimental data, C-Mn steel is used for a sample calculation. A value of 312 kJ/mol for the activation energy, Qa, is chosen, as is explained in Section 5.5.1. By holding T as a constant, n and A may be determined graphically by the relation- ship between ln(sinh(acr)) and In i. The results are shown in figure 6.1. Or similarly, by holding i as a constant, figure 6.2, n and A may be determined by the relationship between ln(sinh(cw)) and 1/T. It is observed, after many testing runs in the simula- 57 Chapter 6. MODEL VERIFICATION^ 58 55 50 45 40 35 30 25 20 15 10 50. L. 0 - _____ 5^1^1 5^5^3 5^4^45 ,C cc = 0.15 in t C cc = 0.01 - -0.5 - • ^ _________0 -^ -4- 1.5 ^ 1.5^4^4.5^5 In e Figure 6.1: Relationship between ln(sinh(acr)) and in i for a C-Mn steel (0.03 C 0.62 Mn) at a strain of 0.35. Experimental data: o — (i = 2s-1), * — (i = 20s-') and + — (i = 140s-1). Data source: [53]. tion with different values of a, that the magnitude of a has no effect on the final stress calculation. In other words, the value for a may be arbitrarily chosen. Note that, since the experimental data shows significant irregularity in the flow stress of this type of steel, the slopes in different regions of the relationship, ln(sinh(acr)) — in or ln(sinh(aa)) — 1/T, are not the same. So n and A are no longer constants. This results in a complication in the modelling. It is also confirmed that, when holding a as a constant, both n and A are non-linear Chapter 6. MODEL VERIFICATION^ 59 55 50 45 40 35 30 ez• 25 20 15 10 5 10000/T --- 1/K 2.5 2 1.5 t5 :44^1 .-E 0.5 0 —0.5 —1 6.5 .s^g.s^4.5^10 1 0000rr --- 1/K 10.5 Figure 6.2: Relationship between ln(sinh(aa)) and 1/T for a C-Mn steel (0.03 C 0.62 Mn) at a strain of 0.35. Experimental data: o — = 2s-1), * (i = 20s-1) and + — (i = 140s-1). Data source: [53]. functions of strain and temperature, which can not be written in a simple analytical expression. In order to calculate the instantaneous values of n and A, the two non-linear functions were established by a computer program based on a simple surface fitting algorithm[59]. A brief description of the algorithm follows. 1. Assuming the surface is a two dimensional polynomial, [n lnA] f(e,T) (6.1) = ao + ale a2T a3e2 a4ET a5T2 Chapter 6. MODEL VERIFICATION^ 60 „,2rp^,,3 E3T2 4-(2,6c.^--r- u7c.. -r as_^, ..• (6.2) where a, (i = 0,1,2, ...) are constants to be determined, and by holding temperature as a constant, an one dimensional polynomial may be obtained, [n lnA] = bo bie b2E2 b3E3^ (6.3) where b, (i = 0,1,2, ...) are constants to be determined. 2. Use Eqn.6.3 to generate enough points at certain strains, including the desired strains points. 3. Repeat the first two steps at different typical temperatures. 4. By holding strain at the desired value to interpolate n and lnA along temperature direction, n and lnA as functions of temperature at the desired strain is then established. 6.1.2 Determination of Flow Stress at any Desired Conditions With the establishment of functions for n and A, Eqn.5.46 may be rearranged as the following, —nh1 i-1 a^ {i exp (Q„k1T-1) a = s ] A (6.4) Now, by using Eqn.6.1 to calculate the values of n and A at the required strain and temperature, the parameters in Eqn.6.4 may all be determined. So a may be determined for that point, (e,i,T). Note that a value for a has to be chosen when determining Eqn.6.1. Figure 6.3 shows a comparison between the model prediction and experimental data for a certain steel. For the working range of the VGR process, i.e., a strain rate Chapter 6. MODEL VERIFICATION^ 61 range of 1-50s-', a temperature range of 800-1250 °Cand a strain range of 0.051— 0.693(reduction ratio 8% - 50%), the errors in the model prediction are within 5%. It is noticed that at 700 °C, the errors are larger. This may be a confirmation to a remark by Hosford and Caddell[60]: such correlations as Eqn.5.46 "may break down if applied over too large a range of temperatures, strains, or strain rates". 6.2 Temperature Model Verification 6.2.1 Verification Using Data from Conventional Rolling The testing data in [37] for workpiece surface temperature and the same set of data of the rolling conditions have been used in order to solve for the temperature model, Eqn.5.22, 5.41 and 5.43. It is found out that the converging speed of Eqn.5.35 is largely dependent upon the size of time steps and the magnitude of the rolling speed. The smaller the size of the time step and the faster the rolls rotate, the slower Eqn.5.35 converges, and vice versa. The calculation of the first node converges very slowly. As the calculation moves on to the second and the rest nodes, the equation converges faster and faster. Figure 6.4 shows the simulation results. For the grid shown and for a converging accuracy of 0.1 °C, the number of terms taken in calculating Eqn.5.35 for the first node is 97, and for the last node is 9. This results in a relatively long time in computation. The mesh plot shows the roll chilling effects clearly. The temperature gradients only exist in the 4-5 layers to the surface. Such a shallow gradient is due to the high rolling speed(note that only about 33 milliseconds of in-bite time). In the surface temperature prediction, large errors, compared to the experimental data, exist in the first few nodes, which most likely are caused by taking a contact of 50(kWm-2K-1) for the heat transfer coefficient at the interface between the workpiece and rolls, H(r). As is indicated in figure 5.8, H(T) is very small at the beginning, which Chapter 6. MODEL VERIFICATION^ 62 would result in a very small temperature gradient at the first few nodes, as confirmed by the experimental data. For the rest of the nodes, good agreement is obtained. It should be indicated that such a fine grid as was used in simulation is not necessary for real applications. Besides, for the VGR process, the rolling speed will be much slower. These would contribute to a significant reduction in computing time. 6.2.2 Simulation of First Pass in VGR Process In order to further verify the capability of the model, simulations for the mean temper- ature distribution along a whole workpiece during and after the first pass of the VCR process have been carried out. The simulation follows the sequence in a real production. That is, first, allowing 5 seconds for transferring the workpiece from furnace(1200 °C) to the mill, the actual entrance temperature may be calculated by the radiation and convection model. Then rolling starts and the temperature at each node in the roll-bite is calculated by the conduction model, while continuing the estimation of temperature drop in the workpiece outside the entrance. At the exit, the mean temperature is estimated by averaging the temperature at the nodes in vertical direction. And then, right out of the exit, starts the radiation and convection calculation again to count in the temperature drop. Finally, allows 5 seconds handling time for the workpiece to be transferred to the entrance for second pass. Figures 6.5 and 6.6 are the results from three different runs, showing the effects of varying rolling conditions on a same workpiece. The same heat transfer coefficient(50 kWm-2K-1) at the interface between workpiece and rolls used in the last Section was chosen. Some enlightening observations for this particular workpiece and working tem- perature range may be summarized in the following. Chapter 6. MODEL VERIFICATION^ 63 1. The rate of temperature drop due to radiation and convection is about 2 - 3 °C/sec. 2. The mean temperature drop is mainly due to conduction in the roll-bite. This again confirms that the roll chilling effects is dominant. The higher reduction ratio, the larger the amount of mean temperature drops. Since higher reduction ratio implies longer time contacting to the rolls. 3. Increasing the rolling speed reduces the contacting time between the workpiece and the rolls, which results in smaller amount of temperature drop. 4. There are two different ways of increasing the rolling speed. One is to increase roll diameter(case B in figure 6.5). The other is to increase roll velocity(case C in figure 6.5). It is clear from figure 6.5 that the second method(case C) is preferred, because in that way the contact length 1 is not increased. 5. Figure 6.6 shows the mean temperature in the roll-bite along the whole pass. Again, it shows that case C(the top curve) is preferred, as smaller changes in mean tem- perature is resulted, which, in turn, would result in smaller variations in flow stress. The calculations were performed by Eqn.5.22 and 5.44. The above observations are important and may serve as guidelines in the VGR process implementation. 6.3 Deformation Model Verification 6.3.1 The Side Spread Estimation In Chapter 4, a deformation model based on upper bound theorem has been worked out. As has been explained in Section 5.2, when looking merely at the deformation, the material may be assumed rigid perfectly-plastic. This means that the flow stress becomes Chapter 6. MODEL VERIFICATION^ 64 a proportional constant and has no influence on the amount of side spread. Figure 6.7 shows the results from the simulation run and a comparison is made with the data from [55]. The agreement is excellent. 6.3.2 Torque Estimation The total net power consumption is given by the minima of Eq.4.3, and torque may be calculated from the following equation, M = —.I* • VR min (6.5) Figure 6.8 gives the results from the simulation run and a comparison is made with the data from [57]. The agreement is very good, as is expected that the results from Eq.4.3 should be the upper bound. 6.3.3 Roll Separating Force Estimation Traditionally, the upper-bound approach can not be directly used in determining the rolling force. The reasons are apparent — the frictional force at the interfaces changes direction due to the changes in relative speeds between the rolls and the workpiece, and there is very little knowledge on this particular subject available. It is unfortunate that, one may find out in the literature, many researchers had to use the slab method to calculate the roll force, when everything else had been worked out by the upper bound approach. Since slab method can only be used under plain strain condition, some researchers, such as [20][6], after devoted painstaking efforts in three dimensional solution through upper bound approach, chose to further simplify the problem into a two dimensional one so the slab method can be used for the solution to the roll force. This resulted in inconsistency in the theoretical analysis and may be considered a shortcoming of the UBA approach. Chapter 6. MODEL VERIFICATION^ 65 In order to overcome such a difficult situation, a concept called equivalent coefficient of friction(peq), similar to the one proposed by D. Y. Yang and J. S. Ryoo[58] when ring rolling process was investigated, will be adopted. This coefficient is defined as the ratio of the frictional force(FT), which is the difference between the foreward friction and the backward friction, to the overall radial force(FR) upon the rolls, figure 6.9. That is, FT peg R (6.6) where FT is also the tangential force on the rolls. Accordingly, (6.7) Therefore, if /leg were determined by experiments, the roll separating force may be de- termined by using Eqn.6.5, 6.6 and 6.7 as, F = V.F1, +^ ( 6.8) Since the actual friction is affected by various factors, such as surface conditions, roll size, temperature, workpiece dimension, etc., it is virtualy impossible to obtain the ex- act distribution of friction. By considering the overall friction effect on the roll force and torque, the extremely complicated situation is simplified. The value of peg may be determined from a set of measurements for the roll force and torque, as suggested in [58]. Note that geometrically(figure 6.9) peg^tanO ^X d^ (6.9) VR2 — then, there must be 0 < xd < 1. As the normal pressure is not uniformly distributed along the contact length(/), the value of xd has to be determined experimentally. This may suggest another approach in computation. Chapter 6. MODEL VERIFICATION^ 66 In the present study, the experimental results for roll force and torque reported in [57] are used as a demonstration, and the values of peg are computed for the same set of data as was used in computing the torque(Section 6.3.2). Figure 6.10 shows the simulation output. For the same reason as described in Section 5.5.2, lead was chosen as the experimental material, and the flow stress model, Eqn.5.48, is used. For reduction ratios ranging from 8% — 50%, the best results are obtained when the values of peg were chosen as a variable within 0.108 — 0.135, which gives the values of the ratio, xd//, within 0.8506 — 0.4285. 6.4 Concluding Remarks 6.4.1 Summary and Conclusions A complete mathematical model for implementing the VCR process under hot work- ing conditions has been established through a fundamental analysis. The model takes into account the three major aspects of the process that the machine designers and the production engineers have to deal with. 1. The deformation model is capable of producing the upper bounds of the rolling loads, namely, the roll separating force and the torque, that will be needed in machine design task. Using the model, the instantaneous rolling loads, which are important information for the roll-gap control, can also be obtained with the help of accurate flow stress data. The force calculation is based on a new concept called the equivalent coefficient of friction, peg. The adoption of this concept leads to two significant advantages. Firstly, the tests for the actual friction coefficient at the interface between the rolls and the workpiece, as in the traditional rolling process, may be avoided. It is known Chapter 6. MODEL VERIFICATION^ 67 that there is no reliable method for such testing. Besides, the available methods usually involve many facilities that are costly, time consuming and inconvenient in practice. The test(See Section 6.3.3) to the /leg, however, will be very simple and reliable. Secondly, this concept leads to a theoretical consistency in the application of the upper bound approach. Further simplifications are no longer needed, as were done by some researchers(See Section 6.3.3). Another new feature related to the deformation model is that a simpler velocity field solution has been proposed. This largely reduced the work in derivation and computation. By such a velocity field, the side spread predictions were shown to be successful. 2. A temperature model has been formulated based on one-dimensional flow assump- tions in the roll-bite and two-dimensional flow outside the roll-bite. The model is capable of predicting the through-thickness temperature distribution in the roll- bite, estimating the mean temperature of the deforming body that is one of the important inputs to the flow stress model, and roughly calculating the the mean temperature distribution along the workpiece. Based on the mean temperature of the workpiece, roll pass scheduling may be performed, along with the side spread estimation. Some preliminary estimation of the final properties of the rolled products may also be possible. 3. The modelling of flow stress for steels is very involved. An accurate flow stress model, which is of great importance in implementing the VCR process, has been chosen. It has been shown that the unified creep relationship, Eqn.5.45, works very well, as is widely reported in the literature. Chapter 6. MODEL VERIFICATION^ 68 It is observed that, for the certain type of C-Mn steel used in this work, the param- eters, n and A are non-linear functions of strain and temperature. Such functions may be determined by numerical methods. Also, the value of 312 kJ/mol for the activation energy, Qa, has been selected, as was reported in the literature. It worked well in this study. 6.4.2 Suggestions for Further Work 1. Further study of the equivalent coefficient of friction by experiments is needed. Without experiments, it is difficult to decided the actual range of this variable. 2. A temperature model for the rolls is needed in order to estimate accurately the temperature changes of the rolls. Although many reasons that in this study the temperature of rolls was assumed as a constant is given in Section 5.4.1, attempts can still be made. 3. The final microstructure of the rolled products will be another issue in the VGR process planning. A mathematical submodel of this purpose may be incorporated in the present model. Strain rate: 20/s 400 350 c" 3001:* 250 g 200 150 100 50 ea 300 250 col 200 150 E-• 100 50 0 Chapter 6. MODEL VERIFICATION 69 200 180 160 140 120 100 80 60 40 20 0 Strain rate: 2/s 14.1^.2^0.3^6.4 0 0 0 0 0 0 0 ()opt* 6.s 6.6^0.7 True strain( ot = 0.025) Figure 6.3: Flow stress model verification: C-Mn steel (0.03 C 0.62 Mn). Data source: [53].^*, x — experimental data. Solid or dotted lines — model predictions. Chapter 6. MODEL VERIFICATION Figure 6.4: Temperature distribution in roll-bite. Data source:[37]. Top: 3D-mesh plot. Middle: grid used. Lower left: front view of the 3D-mesh. Lower right: surface temper- ature comparison. Reduction ratio: 35%. Entry temperature: 1012 °C. 70 1100 1080 1090 104-0 al 1 020 0.1^0.2^0.3^04^0.5^0.0 Position along the workplace. rn 960o aria 0.7 0 Start 1000 980 Roll diameter 0.0635 truster (5 ni) Flail velocity: 20 r-r-np (Q.133 m/s) emperatureSt the. skit of the first pass -- temperature. St the en-hence to the saccinci pass . 1080cr. :g 1000 10-40 1 020 1000 0 Encl 0.1^0.0^0.3^0.4^0.5^0.0 Position along the vvorkpiece 0.7 0 Start 1120 1100 it-1)ctimmeter: O.127 rneter (1-0 velocity: 20 rma (0.206 m/$) Flail Fl all tamperispAre. St trvia exit of the first pass tamper-Attire aat^entriarsCa 1.0 the. sec.:v.-1cl • pass 1110 C..1 1 100 : 1 090 .4M ES' 1 oso 1070 1060 1050 1040 0 S rid 0.1^0.2^0.3^0.4^0.5^0.6 1=.c.sition elc.ing the. workplace 0.7 0 Start 1130 1100 exit of ths first passtemperature. St ths Flail c9airneter: p.osas rrie.t.ar (s it-I) Flail -velocity: 4.0 rrrip (o.kss m/s) emparahrre St the entreince to 'Wis siscc.nd pease Chapter 6. MODEL VERIFICATION ^ 71 Figure 6.5: Sample calculation: mean temperature distribution of the workpiece for three different rolling conditions. Workpiece dimension(before rolling): 21.6 x 21.6 x 800mm. Reduction ratio: 8% — 50.8%. Furnace temperature: 1200 °C. Top(case A), time for passl: 6.015 sec. Middle(case B), time for passl: 3.008 sec. Bottom(case C), time for passl: 3.008 sec. Roll -diameter: .0. Roll velocity: 20 Roll- diameter-:- Roll velocity: 20. RolL diameter:. a. Roll velocity: 40 0635 meter (5 in) ma-1p (0.133 m/s) .127.rneter (10 in) rmp (0.266 m/s) .0635 -meter -(5• iri) rmp (0.266 m/s) 1150 C-) ,th,1140 a> a> 15 1120 .1= a> a 00 cT, 1090 a> 1080 1070 1060 0 End 0.1^0.2^0.3^0.4^0.5^0.6 Position along the workpiece m 1130 111 0 0.7 ^ 08 Start Chapter 6. MODEL VERIFICATION^ 72 Figure 6.6: Sample calculation: effects of changes in geometry and roll velocity on mean temperature of the deforming body in the roll-bite during the first pass. Furnace temper- ature: 1200 °C. Workpiece dimension(before rolling): 21.6 x 21.6 x 800mm. Reduction ratio: 8% — 50.8%. Chapter 6. MODEL VERIFICATION 73 1.45 1.4 1.35 1.3 g 1.25 1.2 1.15 1.1 1.05 1 0.05^0.1^0.15^0.2^0.25^0.3^0.35 Reduction ratio Figure 6.7: Spread prediction — comparison between simulation results and exper- iments[55]. Roll diameter: 101.6 mm(4 in). Specimen dimension: 9.525 x 9.525 mm(0.375 x 0.375 in). 0.4 ^ 0.45 ^ 05 Chapter 6. MODEL VERIFICATION ^ 74 0.1^0.15^0.2^0.25^0.3^0.35 ^ 0.4 ^ 0.45 ^ 0.5 Reduction ratio Figure 6.8: Torque prediction: comparison between simulation results and experi- ments [57]. Roll diameter: 127 mm(5 in). Specimen dimension: 12.7 x 19.05 mm(0.5 x 0.75 in). Chapter 6. MODEL VERIFICATION^ 75 Figure 6.9: Schematic illustration of the roll force vectors layout in flat rolling process Chapter 6. MODEL VERIFICATION^ 76 Figure 6.10: Roll force prediction: comparison between simulation results and experi- ments[57]. Roll diameter: 127 mm(5 in). Specimen dimension: 12.7 x 19.05 mm(0.5 x 0.75 in). Bibliography [1] Roberts, W. L., Cold Rolling Of Steel, Marcel Dekker, INC.1978. [2] Theocaris, P.S., "A Study of the Contact Zone and Friction Coefficient in Hot- Rolling", Metal Forming and Impact Mechanics, Pergamon Press, Ed. by S.R. Reid, 1985 [3] Hill, R., "A General Method of Analysis for Metalworking Processes", Mech. Phys. Solids, Vol. 11, 1963, pp. 305. [4] Lahoti, G. D., Kobayashi, S., "On Hill's General Method of Analysis of Metalworking Processes", International Journal of Mechanical Sciences, Vol. 16, 1974, pp. 521-540. [5] Oh, S. I., Kobayashi, S., "An Approximate Method For A Three-Dimensional Anal- ysis of Rolling", International Journal of Mechanical Sciences, Vol. 17, 1975, pp.293- 305. [6] Kennedy, K. F., "A Method for Metal Deformation and Stress Analysis in Rolling", Doctoral Dissertation, The Ohio State University, 1986. [7] Design And Application Of Leaf Springs — SAE J788a Handbook Supplement, 1978. [8] Sepehri, N., "Computer-Aided Rolling of Parts With Variable Rectangular Cross- Section", MASc. Thesis, The University of British Columbia, 1986. [9] Sassani, F., Sepehri, N., "Prediction of Spread in Hot Flat Rolling under Variable Geometry Conditions", J. Materials Shaping Texhnology, Vol. 5, No. 2, pp.117-123, 1987. [10] Sassani, F., Sepehri, N., "Computer-Aided Process Planning for Rolling of Parts Having Smoothly Varying Rectangular Cross-Section", Int. J. Mach. Tools Manu- fact., Vol. 29, No. 2, pp.257-266, 1989. [11] Sparling, L. G. M., "Formula for Spread in Hot Flat Rolling", Proc. Inst. Mech. Engrs., Vol. 175. pp. 604-640, 1961. [12] El-Kalay, A. K. E. A., Sparling, L. G. M., "Factors Affection Friction and their Effects upon Load, Torque, and Spread in Hot Flat Rolling", J. Iron Steel Inst., Feb. 1968, No. 206, pp. 1110-1117. 77 Bibliography^ 78 [13] Chitkara, N. R. and Johnson, W., "Some Experimental Results Concerning Spread in the Rolling of Lead", Journal of Basic Engineering, Vol. 88, pp. 489-499, 1966. [14] Takahashi, J., Sato, T., Sakai,Y., and Ayada, M., "Production of Taper Leaf Spring By S.P.M.",UDC 621.771.07, 1980. [15] Hawkyard, J. B., Johnson,W., Kirkland, J., Appleton, E., 1973 "Analyses for roll force and toque in ring-rolling, with some supporting experiments", Int. J. Mech. Sci. 15, pp.873-893. [16] Kobayashi, S., Oh, S., Altan, T., Metal Forming and the Finite-Element Method, Oxford University Press, 1989. [17] Hill, R., The Mathematical Theory of Plasticity, Clarendon Press, Oxford, 1950. [18] Li, G.-J., Kobayashi, S., "Spread Analysis in Rolling by the Rigid-Plastic Finite El- ement Method", Numerical Methods in Industrial Forming Processes, p.777, Piner- idge Press, Swansea UK, 1982. [19] Grober, H., "Finite Element Simulatoin of Hot Flat Rolling of Steel", NUMIFORM 86 — Numerical Methods in Industrial Forming Processes, ed. by Mattiasson, K., et al, p.225, 1986. [20] Lahoti, C.D., Akgerman, N., Altan, T., "Computer-aided analysis and design of the shape rolling process for producing turbine engine air-foils", Final report, 1978 December, Contract Nas 3-20380. [21] Sevenler, S., Raghupathi, P. S., Altan, T., "Spread and bulge in Bar and Rod Rolling Using Flat Rolls", Iron and Steel Engineer, March 1986, pp. 57-62. [22] Yih, C. S., "Stream Functions in Three-Dimensional Flow", La Houille Blanche, Vol. 12, 1957, pp. 445-450. [23] Sheppard, T., Wright, D. S.,"Parameters affecting lateral deformation in slabbing mills", Metals Technology, February 1981, pp. 46-57. [24] W. Johnson, P.B. Mellor, Engineering Plasticity, Ellis Horwood Limited,1983 [25]M. J. Luton, J. J. Jonas, Proc. Int. Conf. Strength Metals Alloys, 2nd, 1970, pp. 1100-1105. [26]H. J. McQueen, J. J. Jonas, "Recovery and Recrystallization during High Temper- ature Deformation", Plastic Deformation of Materials Vol. 6, Ed. R. J. Arsenault, 1975, pp. 393-493. Bibliography^ 79 [27] I. Tamura, C. Quchi, T. Tanaka, H. Sekine, Thermomechanical Processing of High Strength Low Alloy Steel, 1988. [28] W. J. McG. Tegart, A. Gittins, "The Hot Deformation of Austenite", The Hot Deformation of Austenite, Ed. J. B. Ballance, 1977, pp.1-46. [29] H. Kudo, "Upper Bound Approach To Metal Forming Processes — To Date And In The Future", Metal Forming and Impact Mechanics, Ed. S.R. Reid, Pergamon Press, 1985, pp. 19-45. [30] H. Kudo, "Some Analytical And Experimental Students of Axi-symmetric cold forg- ing and Extrusion—I", International Journal of Mechanical Sciences, Vol. 2, 1960, pp. 102-127. [31] W. Johnson, J. Inst. Met., Vol. 85, 1956, pp. 403. [32] L. C. Dodeja, W. Johnson, J. Mech. Phys. Solids, Vol. 5, 1957, pp. 281. [33] Avitzur, B., Handbook of Metal-Forming Processes, John Wiley & Sons, 1983. [34] Thomsen, E. G., Yang, C. T., Kobayashi, S., Mechanics of Plastic Deformation in Metal Processing, the Macmillan Ltd., 1965. [35] J. Chakrabarty, Theory of Plastisity, McGraw-Hill, 1987. [36] C. Devadas, I. V. Samaradekera, "Heat transfer during hot rolling of steel strip", Ironmaking and Steelmaking, Vol.13, No. 6, 1986, pp.311-321. [37] C. Devadas, "The Prediction of the Evolution of Microstructure During Hot Rolling of Steel Strips", Ph.D. thesis, the University of British Columbia, June 1989. [38] D. Partington, L. Talbot, "Computer-aided draftiong control in plate rolling", Hot Working and Forming Processes, Ed. C. M. Sellers, et al, 1979, pp716-180. [39] G. F. Bryant, M. 0. Heselton, "Roll-Gap Temperature Models for Hot Mills", Metals Texhmology, Vol. 9, Dec. 1982, pp.469-477. [40] F. Hollander, "A model to calculate the complete temperature dsitribution in steel during hot rolling", The Iron and Steel Institute, London, 1969, pp.46-78. [41] Roberts, W. L., Hot Rolling Of Steel, Marcel Dekker, INC.1983. [42] B.I.S.R.A., Ed., Physical Constants of some Commercial Steels at Elevated Temper- atures, Butterworths, London, 1953, pp.1-38. Bibliography^ 80 [43] J. Woolman, R. A. Mottram, The Mechanical and Physical Properties of the British Standard EN Steels, the Macmillan Co., 1964. [44] R. I. L. Guthrie, Engineering in Proess Metallurgy, Oxford Science Publications, 1989. [45] J. P. Holman, Heat Transfer, 6th. edn., McGraw-Hill, 1986. [46] C. Devadas, I. V. Samaradekera, and E. B. Hawbolt, "The Thermal and Metal- lurgical State of Steel Strip during Hot Rolling: Part I. Characterization of Heat Transfer", Metallurgical Transactions, Vol. 22A, Feb. 1991, pp.307-319. [47] Ashok Kumar, I. V. Samaradekera, and E. B. Hawbolt, "Roll-bite deformation furing the hot rolling of steel strip", Journal of Materials Processing Technology, Vol.30, 1992, pp.91-114. [48] F. Kreith, W. Z. Black, Basic Heat Transfer, Harper and Row, New Youk, 1980. [49] K. Murata, H. Morise, M. Mitsutsuka, H. Haito, T. Komatsu, and S. Shida, Trans. Iron Steel Inst. Jpn., Vol. 24, (9), 1984, pp.B309. [50] P. G. Stevens, K. P. Ivens, and P. Harper, Journal of Iron and Steel Institute, 1971, Vol. 209, pp.1-11. [51] M. A. Pinsky, Introduction to Partial Differential Equations with Application, McGraw-Hill Book, 1984. [52] C. M. Sellars, "The Physical Metallurgy of Hot Working", Hot Working and Forming Processes, Ed. C. M. Sellers, et al, 1979, pp.3-15. [53] M. J. Stewart, "Hot Deformation of C-Mn Steels from 1100 to 2200 °F (600 to 1200 °C) with Constant True Strain Rates from 0.5 to 140 s'", The Hot Deformation of Austenite, Ed. J. B. Ballance, 1977, pp.47-67. [54] S. F. Wong, P. D. Hodgson, P. F. Thomson, "Lead As A Model Material For The Hot Rolling Of Steel", Mathematical Modelling of Hot Rolling of Steel, Ed. S. Yue, 1990, pp.281-289. [55] Chitkara, N. R. and Johnson, W., "Some Experimental Results Concerning Spread in the Rolling of Lead", Journal of Basic Engineering, Vol. 88, pp. 489-499, 1966. [56] Loizou, N. and Sims, R. B., "The Yield Stress of Pure Lead in Compression", Journal of the Mechanics and Physics of Solids, Vol. 1, pp. 234 to 243, 1953. [57] Sims, R. B., "The Calculation of Roll Force and Torque in Hot Rolling Mills", Proc. Inst. Mech. Engng, Vol. 168, pp. 191. 1954. Bibliography^ 81 [58] Yang, D. Y., Ryoo, J. S., "An Investigation into the Relationship between Torque and Load in Ring Rolling", Transactions of ASME, Journal of Engineering for In- dustry, Vol. 109, pp. 190-196, Aug. 1987. [59] Gerald, C. F., Wheatley, P. 0., Applied Numerical Analysis, 4th edn. 1989. [60] Hosford, W. F., Caddell, R. M., METAL FORMING — Mechanics and Metallurgy. Prentice-Hall, Inc., 1983. S=--Xn Appendix A ANALYTICAL SOLUTION TO THE UPPER BOUND EQUATIONS A.1 The neutral point Differentiating Eqn.4.39 gives, aQ(s,i) = VRt(xo) Ow(x) awi^t' (x„) awi Or, 5Q (x ) ^ aw(x,i) awl ^w(x) awl and 5Q(x) n (1 Si + 1 aW^t' at' ax„^t\t axn w ax„ 1 + axn) A.1.1 The derivative of Ed against w1 (A.4) where, for the sake of simplicity, define Eq.4.14 as Z(x)^Z3(x) Z4(x) Z5(x), and after some manipulations, one may have the following, Z3(x ) = Z4(x) = (/2^ln + Z2a^3^f3 32 {z2 Ay(2.12 + B2)] (A.5) (A.6) Z5(x) = 12 I ( 2/ —1) tan_1 D T 2) tgia^3^a 3'y (A.7) From Eq.4.13 and Eq.A.2, the following results, aEd 1^aw(xn) E 2 (az3 az, az5 clsawl d w(x)^awl Q o,., owl owl owl 82 Appendix A. ANALYTICAL SOLUTION TO THE UPPER BOUND EQUATIONS83 where, let -y = Nt^a = Mw^ (A.8) 0 = V/2 +72 (A.9) B = VI2 + a2^ (A.10) (Z2 – B)I D = (A.11) (B + a)-y Z2 --= VB2 + 72^ (A.12) A = in '..r + Z2 B (A.13) tg= tan- ^ (A.14) 7 and differetiation of Eq.A.8 — A.14 against w1 results in the following, ON — ton^ = tn^ (A.15)uwi^uwi 0p 1 (aI2^ON),=^+ 27t (A.16) owl^20 awi^awl Oa am A, Ow^ „ (A.17) awl^w awl + 1" awl aB2^ar^aa,^. ,, + 2a ^ (A.18) owl^owl^owl. 0z21^ (aB2^ON),.,^= + 2-yt (A.19) awl^2Z2 atpi^0w1 aA, ^1 (3y az2) 1 aB awl^7 + Z2 alVi + aWi ) Bawl^ (A.20) atg 1 ^f 03 1) a-y^7 [-yt ON + 2 (i.^1)0121} Owl^72 + (0 –1)2 1^awl^1. (A.21) 0 awl^/3 I) awl ap aB) a' I. ,^+(z2 B) , wawl^[(B + 1 42)7]2 {(B + a)7 [I (0,Z2 owl owl ol –(Z2 – B)I [ry ( aawB 4_ Oa \ + (B + a) 0-y 11 (A.22) 1^awl)^awl.) j Now, with the help of the Eq.A.8 — A.21, the partial derivative of Eq.A.5, A.6 and Q^wt (p4 = NIA + (Pi W (P2 -Wt VR V 1 + t '2 (P3 = Appendix A. ANALYTICAL SOLUTION TO THE UPPER BOUND EQUATIONS84 A.7 may be expressed as the following, az3 1 { Z3 Oa + (312 +,-), al ) ln a + Z2= awl^a^awl \awl 3 awo^0 + (12 + -12.) r 1 ^( Oa + 3z2) 1 aol 1 (A.23) 3 ) [a + Z2 W1 OW1 ) 13 awl jj 0Z4^2 1 3z2 +1 [(212 + B2) (aA^ A a-y ),^+ A (2°,,I2 + 3,B2)] }(A.24)= awl^3 awl 'I^awl / owl owl owl az5 25 ar + 12 ( 1 { 1 (^aa^1= c, . a arawi 1 atvi) [2 tan-1 D + (3-7 — 2) tgi 0wi^/2 awl +I Ill _ 2) Otg^tg a-y 1 1 ^1 ^(21 1) aD) (A.25) R3-y^) Owl 323J f 1 + D2 a 3) awi Thus, Eq.A.4 is completely determined. A.1.2 The derivative of Ef against wi. From Eq.4.20, one may have the following, anEf — M k (naQ 11 F(x)Vi + t̂'2dx + Q II V1 + t12 !F dX) 01111^OW1 0^ 0^OW1 where, let, (A.26) (A.27) (A.28) (A.29) then, Eq.4.21 becomes, F(x) =w (v4+ C-a-1 + ln CP2 + (P4 cp2^l(P3i Partial differetiating of Eq.A.27 — A.30 against w1 results the following, a(p2 = 1 ( 312^,Ow awl^tw2 ell awl w atvi ) (A.30) (A.31) Appendix A. ANALYTICAL SOLUTION TO THE UPPER BOUND EQUATIONS85 VR a Q^aw Q2 awl + tw2 awl 1^5992^aso3 (104 402 ^ +403 +493atoi) ( 404+^+ln Co2 + (iO4 Ow^1 ^(9:03)21 (102^l^ (P2(P31^awl + (104 (2CO3^1 3403 +(+ ^1 ^5404 k 402^403) °W1^402 + (P4) OW1 1 note that for Eq.A.30, when co3 <0, 14031 = —w3, then, 595,3 atvi 3994 awl SF awl (A.32) (A.33) 5402 awl (A.34) 514031 _ _aw3 (A.35) awl — awl therefore, we have, 1 314031 _ 1 3403 (A.36) 1403I awl^(,03 awl A.1.3 The derivative of È, against From Eq.4.26 and Eq.4.27, one may have the following, ^aQ^as kQ (A.37) awl^Q awi^awl where, let, S = 1/1D'2 t'2^ (A.38) then, Os^w' Ow'„^ (A.39) awl^s owl and the following may be obtained with the help of the above two equations, as^ , 1= 1 ^w 2 t 7 -= e ( [2 ln(ti + s) — w^s)sawl^3 Appendix A. ANALYTICAL SOLUTION TO THE UPPER BOUND EQUATIONS86 w + s+(ti )2 ( ^1 ^+ 1) + 2 in s —̂ in ktvi)^ivi(tvi+s)^s)^3^w'^It'l 2 ( 1 \ 1^&Lilo' 3 .s(s — ts) w1) fs=o awl (A.40) A.1.4 The derivatives of Ed, Ef and E, against xn 1. From Eq.4.13, one may have the following, atd^[ 5Q axn^3 axn o = —a-- Z(x)dx„ by substituting Eq.A.3 into the above equation, one may obtain, aEd^(1 at^1 aw^t' at' ) axn^t axn + w axn 1 + t'2 axn Ed 2. From Eq.4.20, one may have the following, atf^c2 F(x)i+ti2dx +Q vq+ti2aF dx) axn mn' aOxn 0V ^ 0^axn and note that from Eq.A.28 and Eq.A.29, VR aQ Q2 ax„ VR (p3 aQ Q2 (,04axn X=Xn ap3 awl acp4 awl (A.41) (A.42) (A.43) (A.44) (A.45) Now, with the help of Eq.A.44 and Eq.A.45, partial differentiation of Eq.4.21 results the following, SF VR CO3 ( _4_^1 2W3 1 1 aQ axn 111(22 L(P4 (P2 + CO4 + (P2 (p3.1 axn thus, by substituting Eq.A.46 and Eq.A.3 into Eq.A.43, one may obtain, (A.46) atf^(1 Si^1 aw^t' at')^{" i — VR \ + t2 • Ox,. aXn + W aXn^ii2 aXn ) 0 X=Xn [(p3 (1 + ^1 ^)^4)3^1 I cis} CO4^CO2 + (004^(102^(t03 (A.47) Appendix A. ANALYTICAL SOLUTION TO THE UPPER BOUND EQUATIONS 87 3. Partial differetiating Eq.4.26 and then using Eq.A.3, one may have the following, at,^ac2= kS ax„^ax„ ^ (1 at^1 aw^t'^at'= t ax,, +w ax„ 1 + t'2 ax„ .... t„^(A.48) Appendix B TRUE STRESS & STRAIN CORRELATION FOR C-Mn STEELS 88 30 20 700 C 18 25 16 20 14Strain rare: 20/s 15 10 12 _.1000 C 8^10 ^•^ 800C . —7-77 .... •^ '900.C.^ 8 1200 C 6 5 4 2 0.1 0.2 0.3^0.4 03 0.6 oo07 0.1 0.2 0.3^0.4 0.60.5 07 True strain True strain 30 22 20 - 18 - 16 r 14 - 121- 10- ......... : ... . r . 900C _ 1000 C ................... • 800 C 15 1200 C 5 0 2 Oo0.20.1 0.3^0.4 0.603 07 0.1 0.60.30.2 07 True strain True strain Strain rate: Strain rate: 2/s -- --^...... . . •^....^.... .• .... . • 0.3 0.4 0 .50.1 0.2 25 •10 0.1 0.2 0.3^0.4 'True strain ---------- • 900C • ^•^ 800 C• ^ _^- - -a ---^- Haw - 0.1 05 0.60.2 070.3^0.4 40 35 30 20 15 10 700 C. ------------ ----- • ---- Strain rate: 20/s 1200 C True strain 55 700 C Strain rate: 140/s 900C • •__^-^ ---------- ------ •• ••^ 1 860 C............. 15 v." 35 5 0 05 070.6 1100^_ ......^ ......... • ---------- --------- 1200 c oo 07 •• 4 0. 1 0.2 700 C • • •^• / 57.. if if .5. Strain rate: 0.5/s 0.3^0.4 True strain 2 0 0.5 • 800 C 0.6 0 7 18 16 14 12 10 6 4 900 C• 1000 C ..,-^, o'w^l'.. ...^x,.. a. .....___,.....--.4. ........... - -^...... 1000 C 800 C ..... • ..... • ..... , ....^ _,. ...... .;% .....4.-..4-..-*-.- 4-------------.------.-1-1.°-. 0 C l ?' f:!1 4. i ...''. 4.--.--&-----•____, ...... ..... . .... .^•^• 1?..00 C True strain 30 25 20 15 10 5 30if • 800 C• Strain rate: 20/s 0.6 070.1 ^ 0.2 ^ 0.3^0.4 True strain 45 40 35 25 20 15 10 5 0 45' • 0 60 50 q=- 30 coco 4 20 0 10 oo 700 C 0.1 ^ 0.2 ^ 0.3^0.4 ^ 05 ^ 0.6 ^ 07 True strain • C 900C 700 C Strain rate: 1401s • ................ ^........ ^• 1000 Ca • 15 e" 7?(.: 10 -Is# ........ oo ........ -&1000 C 1100 C .... ...... . • .. .......^1200^_ 5 • - Strain rate: Vs 0.1 0.2 0.3^0.4 0.5 0.70.6 30 700 C ....*, .......... .100 ._ .... ... 900 C 25 20 True strain _________^.•^ 900 C• CO3 45 ^ c ao 2535 Strain rate: 20/s 30 . .^1000C^1 20 • . 900C^i I1100 C - a 10(10.0 a• . 4 10 15 1200C 10 1 05^0.6^07 0.1^0.2^0.3^0.4 0.5^0.6^07 True strain Trite strain 60 oo 0.1^0.2^0.3^0.4 35 Strain rate: 140/s SOO C 50 30 ao ...... - - - ^.1000 30 ..^• • .......^. ......^ ....... • . .^900 C^_ .^ 1000C V. 20 .1100C 1200 C 10 ........... . ^)200 C 00 0.1^0.2^0.3^0.4 True strain 0.5 0.6^07 0.3^0.4 0.6 07 True strain .......................... ..... ..25
Cite
Citation Scheme:
Usage Statistics
Country | Views | Downloads |
---|---|---|
United States | 8 | 0 |
Czech Republic | 6 | 0 |
France | 3 | 0 |
Syria | 1 | 0 |
China | 1 | 11 |
Japan | 1 | 0 |
India | 1 | 0 |
City | Views | Downloads |
---|---|---|
Ceska | 6 | 0 |
Unknown | 4 | 3 |
Charleston | 4 | 0 |
Mountain View | 3 | 0 |
Noida | 1 | 0 |
Beijing | 1 | 0 |
Tokyo | 1 | 0 |
Ashburn | 1 | 0 |
{[{ mDataHeader[type] }]} | {[{ month[type] }]} | {[{ tData[type] }]} |
Share
Share to: