Microsoft Word - ETASR-khromov.doc ETASR - Engineering, Technology & Applied Science Research Vol. 2, No. 6, 2012, 320-324 320 www.etasr.com Khromov and Kawalla: Simulation of a Steel Wire Straightening Taking into Account … Simulation of a Steel Wire Straightening Taking into Account Nonlinear Hardening of Material Ilya Khromov Dept. of Technical Mechanics and Science of Machines Sevastopol National Technical University Sevastopol, Ukraine i.v.khromov@sevntu.com.ua Rudolf Kawalla Institute of Metal Forming TU Bergakademie Freiberg Freiberg, Germany Rudolf.Kawalla@imf.tu-freiberg.de Abstract—Straightening is used to ensure the straightness of a rod workpiece as well as to redistribute or reduce the residual stresses in a finished product. The accuracy of modeling the operation of this kind plays an important role when manufacturing high quality ropes and cables of a small diameter. The physical model of alternating wire bending in a straightener has been considered in the paper. Computer methodology has been developed to calculate deflected mode and internal bending moment in a cross section of a wire at the different stages of deformation, which allowed taking into account nonlinear hardening of the material. Keywords- elastic-plastic deformation; deflected mode; nonlinear hardening; straightener I. INTRODUCTION Straightness is one of the most vital quality factors when manufacturing metal ropes and cables of a small diameter. It is evaluated by the residual curvature diameter of the loose length of a finished product. The value of the parameter mentioned depends significantly on the initial curvature of a wire used in cable/rope twisting. With the view of curvature reduction and technological stress redistribution different types of machining are applied in the process of production, straightening relates to the operations of this kind. Despite the fact that straightening seems to be a rather simple process in terms of mechanical implementation, it is difficult enough to find thorough analysis of the process simulation related to steel wire rope production. The earliest works on solving technological problems and designing the forming process of twisted products were focused on calculations in the absence of computers, so they were based on a very rough simplification. For example, Dorodnykh [1] and Tulenkov [2] used a non-hardening model of material when modeling straightening process. Guericke and et al. [3, 4] proposed the use of a linear hardening model to describe the behavior of a wire being deformed by rollers of a straightener, which refined significantly the numerical calculations. Lee et al. [5] applied FEM to simulate the straightening of metal cord using the software package Abacus. The same approach was employed by K. K. Adewole and S. J. Bull [6] when studying alternative bending of a flat wire. Very often there is no possibility to apply existent program packages for technological problem solving and a research team has to develop their own specialized procedures for finding the optimal deformation parameters when designing technological processes for steel rope manufacturing. However, in such cases, technical articles are usually limited to a brief description of the models used and focus more on explaining the interface features of the developed software [7]. In this paper, the authors tried to describe in more detail the main stages of developed computer methods for calculating the deflected mode of a steel wire while straightening and taking into account its real stress-strain diagram. II. OPERATION PURPOSE AND PHYSICAL MODEL OF THE PROCESS Straightening is a multiple alternating elastic-plastic bending of a rod workpiece on the rolls of a straightener, as shown in Figure 1. This kind of machining is used to ensure the straightness of a rod workpiece as well as to redistribute or reduce the residual stresses in a finished product. Theoretical analysis and calculation of setup modes for a straightener (i.e. center-to-center distances horizontally hA and vertically vA ) are based on the physical model of the process, its graphical scheme shown in Figure 2 as a multistep diagram )(kM x . While being straightened, the rod is subjected to pure bending deformation. However, alternating deformation behavior leads to a considerable complication of the mathematical model of this process. Intermediate positions of random cross section of a rod workpiece moving along the straightener are marked with points 1,2,…,7 in Figure 1a. Respective deflected modes of the cross section are shown in Figure 2b as points 1-7. Separate stages of elastic-plastic deformation (lines 02, 23, 34, 45) appear to be described by different nonlinear integral functions as ),( 1, kkM iix  , where i is the number of the stage of bending deformation. At the two final stages (lines 56 and 67 in Figure 2b) the diagram has linear character ),( 55, kkM x . In this case, the physical model allows to set down an implicit equation for defining unknown values of the workpiece Research has been carried out under financial support by DAAD. ETASR - Engineering, Technology & Applied Science Research Vol. 2, No. 6, 2012, 320-324 321 www.etasr.com Khromov and Kawalla: Simulation of a Steel Wire Straightening Taking into Account … curvature at the four intermediate stages of bending 5432 ,,, kkkk . It might be given on condition that there is no residual curvature of a finished product (straightness condition) 0),(),,,,( 55,5432057  kkMkkkkkMM x . (1) Missed relations between parameters 5432 ,,, kkkk can be ascertained from kinematic and technological problem specification. Due to the implicit form of the equation (1) mathematical modeling of the process comes to step-by-step numerical computation and construction of the curve under Figure 2b for the variable set of required parameters 5432 ,,, kkkk . The algorithm of iterative procedure should ensure the fulfillment of (1) as the final result. Thus the main problem of the straightening process design is to guarantee that there is no internal bending moment and residual curvature in the finished product through several stages of elastic-plastic alternating bending of a rod workpiece. a) b) Fig. 1. a) Kinematic scheme of the process and b) appearance of the device for straightening a long-length product Fig. 2. a) Loading condition and b) deformation diagram of a wire in the course of straightening III. EXPERIMENTAL STRESS-STRAIN DIAGRAM OF A STEEL WIRE, APPROXIMATION FUNCTION Stress-strain and shear diagrams have fundamental significance not only for obtaining initial information relating to the material behavior under simple types of loading, they also form the basis for constructing mathematical models for deformable body of any configuration under different types of complex loading. A real stress-strain diagram of a wire of 1,1 mm in diameter made of steel C60 with ultimate strength of 1750u  MPa has been applied for further simulation. It was obtained in the set of tests conducted at the Institute of Metal Forming, TU Bergakademie Freiberg. Mechanical properties for steel wire samples are shown in a table 1. TABLE I. MECH ANICAL PROPERTIES OF WIRE M AD E OF STEEL C60 Young’s modulus, MPa Tensile strength, MPa Elongation, % Yield strength, MPa 190000 1750 1,5 1020 According to the analysis made in [8], it is reasonable to describe the experimental stress-strain diagram of a steel wire using the following combined analytical function             )]1(1[,,)( )1( yy y   abEif , (2) where y , y - yield stress and strain, a, b - invariables determined by three base points of the experimental diagram (for this curve 312,0a , 821,0b ). In this case approximation function almost completely matches the experimental points of the )( diagram (Figure 3). Fig. 3. Experimental stress-strain diagram and approximation function ETASR - Engineering, Technology & Applied Science Research Vol. 2, No. 6, 2012, 320-324 322 www.etasr.com Khromov and Kawalla: Simulation of a Steel Wire Straightening Taking into Account … IV. METHODS OF COMPUTER SIMULATION OF A WIRE STRAIGHTENING PROCESS In accordance with the classification provided in [9] modeling of the straightening process comes to solution of the direct problem of technological mechanics of a rod. In this case, having known the equations for the trajectory of deformation of the longitudinal axis of a rod it is necessary to determine the stress state and internal forces under the arbitrary value of a certain process parameter (e.g. curvature of a rod) The final stress state and rod form i.e. the technical parameters of the finished product are assumed to depend on the history of its deformation in all operations of technological process. In this connection, the method of mathematical simulation should include stepwise calculation of deflected mode changes of the rod material taking into consideration its non-linear properties. When modeling plane sections hypothesis is applied [10], as well as the following hypotheses and theories are accounted:  discrete model is used when cross section of a rod is divided into finite number of elementary areas. Stresses at every point of the area are supposed to be constant. Therefore, one part of the cross section might be at the elastic stage while the other is at the plastic one;  material behavior under pure tension is described by piecewise function (2);  at the plastic stage, the theory of plastic flow is used for calculating stresses. Due to this theory, the final stresses depend on the trajectory of deformation while equations describing plastic deformation are differential relations that are nonintegrable analytically [11];  as criterion for plasticity, the condition of tangential stress intensity (von Mises condition) is applied. Every stage of deformation is divided into a number of steps with curvature increment k . The algorithm supports multiple step-by-step calculation of the stresses on the basis of the unified procedure elaborated previously by Khromov [12]. Having calculated the stress value the internal bending moment can be defined according to the integral formula   F x dFyM  . (3) The main difference of the proposed procedure for stress calculation at the plastic deformation stage is the fact that the nonlinear hardening of the material has been taken into account when solving plastic flow equations. In this case hardening modulus and hardening coefficient depend on accumulated intensity of plastic deformation. Considering (2) at the optional loading stage, they are determined as follows:                        1 y y )ln(, ,       aabEEif d d E , (4)                          1 y y )ln(, 1 , )(      aabif E E . (5) The flowchart of the general algorithm is shown in Figure 4. It includes several one-type cycles executed in consecutive order, their number being equal to the number of alternating bending stages. Data files of internal bending moment, curvature and normal stresses are formed in the course of calculation. After program finishing, the files can be read and processed using built-in functions of standard mathematical program packages. Fig. 4. Flowchart of the general algorithm for solving the problem of alternating bending. ETASR - Engineering, Technology & Applied Science Research Vol. 2, No. 6, 2012, 320-324 323 www.etasr.com Khromov and Kawalla: Simulation of a Steel Wire Straightening Taking into Account … V. RESULTS AND DISSCUSION Figures 5 and 6 represent an example of simulating the fourfold straightening process for a steel wire of 1,1 mm in diameter. Calculations and data file forming have been carried out for three values of initial curvature of a wire ( 5,20 k m-1, 50 k m-1, 100 k m-1) using the procedure developed in Microsoft Visual Basic. The curvature value at the end of a single stage of straightening is unknown. It should be determined by means of the condition that internal bending moment xM and residual wire curvature at the end of the process would be equal to zero under optional value of the initial curvature 0k . In this example, the deformation value for all three initial curvatures remained the same at the separate step of bending. The data files with the calculated arrays of normal stress, bending moment and curvature were processed using the mathematical package MathCad. As it is shown in Figure 5, in this numerical calculation the function )(kM x comes to the coordinate origin with the accuracy sufficient for practical calculations when initial curvature 5,20 k m-1, 50 k m-1. For 100 k m -1 the final value of the moment is not zero. In this case, it is necessary either to increase the value of the bending deformation at the first stage, or add one more step of bending to ensure the output condition 0xM . Fig. 5. Dependence )(kM x calculated for three values of initial curvature Figure 6 represents diagrams for ),( yx stress distribution along the cross section of a wire for three values of initial curvature. The left part of it shows initial stresses when curvilinear wire is straightened at the entrance to the straightener. Stress diagram is linear, outer fiber stresses reach maximum values. After straightening normal stresses are redistributed, stress diagram becomes complex. However, initial maximum stresses are reduced 2-3 times. Even for the case of 100 k m -1 residual stresses do not exceed 300 MPa and this value could be decreased more by adding one more stage of bending. a) b) Fig. 6. Graphical stress analysis of a wire when a) entering and when b) leaving the straightener VI. CONCLUSION The technique developed provides a method for choosing the necessary number of bending cycles at the stage of technological procedure designing based on real properties of a material as well as for visual evaluation of distribution of stresses in cross section of a wire at different stages of deformation. Further the different elements described can be used for modeling one of the perspective operations of machining which is a rotating straightening. ACKNOWLEDGMENT Authors appreciate Dipl.-Ing. Marcel Graf for technical support while preparing and making experiments. ETASR - Engineering, Technology & Applied Science Research Vol. 2, No. 6, 2012, 320-324 324 www.etasr.com Khromov and Kawalla: Simulation of a Steel Wire Straightening Taking into Account … REFERENCES [1] B. Dorodnykh, “Method for calculation of set-up parameters of deforming device when manufacturing two lay ropes”, Mining Journal, No. 9, pp. 60-66, 1959 [in Russian] [2] F. Tulenkov, “Change of stress condition of a wire in the process of straightening”, Steel Ropes, No. 2, pp. 367-372, 1965 [in Russian] [3] W. Guericke, M. Paech, E. Albert, “Simulation of the wire straightening process”, Wire Industry, No. 8, pp. 613-620, 1996 [4] W. Guericke, “Material model describing cyclic elastic-plastic deformation of roller levelling and straightening processes”, Steel Research Int., Vol. 80, No. 4, pp. 281-287, 2009 [5] J. S. Lee, H. Huh, J. Bae, J. Lee, D. Kim, “Design optimization of roller straightening process for steel cord using response surface methodology”, Materials Processing And Design; Modeling, Simulation and Applications; NUMIFORM '07; Proceedings of the 9th International Conference on Numerical Methods in Industrial Forming Processes. AIP Conference Proceedings, Volume 908, pp. 581-588, 2007 [6] K. K. Adewole, S. J. Bull, “Simulation of the wire reverse bending test”, Journal of Civil Engineering and Construction Technology, Vol. 3, No. 1, pp. 116-126, 2012 [7] M. Paech, “The positioning of straightening rolls”, Wire, Vol. 52, No. 2, pp. 82-83, 2002 [8] I. V. Khromov, “Methods for selecting parameters of exponential for stress-strain diagram of a steel rope wire”, SevNTU Journal. Mechanics, Power Engineering, Ecology, No. 87, pp. 12  15, 2008 [in Russian] [9] I. V. Khromov, “Classification of processes of plastic working of rod workpieces with permanent cross-sectional shape”, SevNTU Journal. Mechanics, Power Engineering, Ecology, No. 106, pp. 31-34, 2010 [in Russian] [10] V. D. da Silva, Strength of materials, 2006 [11] A. S. Kobayashi, S. N. Atluri, “Analytical mechanics of solids”, Springer handbook of experimental solid mechanics, pp. 3-16, 2008 [12] I. V. Khromov, “Analysis of the stressed state of a circular thin-walled cylinder under complex loading and with nonlinear hardening of the material”, Strength of Material, Vol. 41, Issue 3, pp. 278-284, 2009