^ Since 1955 _ Strojniški vestnik Journal of Mechanical Engineering Strojniški vestnik - Journal of Mechanical Engineering (SV-JME) Aim and Scope The international journal publishes original and (mini)review articles covering the concepts of materials science, mechanics, kinematics, thermodynamics, energy and environment, mechatronics and robotics, fluid mechanics, tribology, cybernetics, industrial engineering and structural analysis. The journal follows new trends and progress proven practice in the mechanical engineering and also in the closely related sciences as are electrical, civil and process engineering, medicine, microbiology, ecology, agriculture, transport systems, aviation, and others, thus creating a unique forum for interdisciplinary or multidisciplinary dialogue. The international conferences selected papers are welcome for publishing as a special issue of SV-JME with invited co-editor(s). Editor in Chief Vincenc Butala University of Ljubljana, Faculty of Mechanical Engineering, Slovenia Technical Editor Pika Škraba University of Ljubljana, Faculty of Mechanical Engineering, Slovenia Founding Editor Bojan Kraut University of Ljubljana, Faculty of Mechanical Engineering, Slovenia Editorial Office University of Ljubljana, Faculty of Mechanical Engineering SV-JME, Aškerčeva 6, SI-1000 Ljubljana, Slovenia Phone: 386 (0)1 4771 137 Fax: 386 (0)1 2518 567 info@sv-jme.eu, http://www.sv-jme.eu Print: Grafex, d.o.o., printed in 380 copies Founders and Publishers University of Ljubljana, Faculty of Mechanical Engineering, Slovenia University of Maribor, Faculty of Mechanical Engineering, Slovenia Association of Mechanical Engineers of Slovenia Chamber of Commerce and Industry of Slovenia, Metal Processing Industry Association President of Publishing Council Branko Širok University of Ljubljana, Faculty of Mechanical Engineering, Slovenia Vice-President of Publishing Council Jože Balič University of Maribor, Faculty of Mechanical Engineering, Slovenia Cover: Experimental identification of mechanical characteristics of sheet metal is typically performed in the in-plane direction, demanding special preparation of specimens. An out-of-plane approach through sheet-metal blanking has been proposed. A cutting module (top) is used on a purposely-build laboratory blanking apparatus to determine material characteristics in the out-of-plane direction through analysis of full and partial (bottom) blanking experiments. Image Courtesy: LADISK, Faculty of Mechanical Engineering, University of Ljubljana ISSN 0039-2480 International Editorial Board Kamil Arslan, Karabuk University, Turkey Josep M. Bergada, Politechnical University of Catalonia, Spain Anton Bergant, Litostroj Power, Slovenia Miha Boltežar, UL, Faculty of Mechanical Engineering, Slovenia Franci Čuš, UM, Faculty of Mechanical Engineering, Slovenia Anselmo Eduardo Diniz, State University of Campinas, Brazil Igor Emri, UL, Faculty of Mechanical Engineering, Slovenia Imre Felde, Obuda University, Faculty of Informatics, Hungary Janez Grum, UL, Faculty of Mechanical Engineering, Slovenia Imre Horvath, Delft University of Technology, The Netherlands Aleš Hribernik, UM, Faculty of Mechanical Engineering, Slovenia Soichi Ibaraki, Kyoto University, Department of Micro Eng., Japan Julius Kaplunov, Brunel University, West London, UK Iyas Khader, Fraunhofer Institute for Mechanics of Materials, Germany Jernej Klemenc, UL, Faculty of Mechanical Engineering, Slovenia Milan Kljajin, J.J. Strossmayer University of Osijek, Croatia Janez Kušar, UL, Faculty of Mechanical Engineering, Slovenia Gorazd Lojen, UM, Faculty of Mechanical Engineering, Slovenia Thomas Lübben, University of Bremen, Germany Janez Možina, UL, Faculty of Mechanical Engineering, Slovenia George K. Nikas, KADMOS Engineering, UK José L. Ocana, Technical University of Madrid, Spain Miroslav Plančak, University of Novi Sad, Serbia Vladimir Popović, University of Belgrade, Faculty of Mech. Eng., Serbia Franci Pušavec, UL, Faculty of Mechanical Engineering, Slovenia Bernd Sauer, University of Kaiserlautern, Germany Rudolph J. Scavuzzo, University of Akron, USA Arkady Voloshin, Lehigh University, Bethlehem, USA General information Strojniški vestnik - Journal of Mechanical Engineering is published in 11 issues per year (July and August is a double issue). Institutional prices include print & online access: institutional subscription price and foreign subscription €100,00 (the price of a single issue is €10,00); general public subscription and student subscription €50,00 (the price of a single issue is €5,00). Prices are exclusive of tax. Delivery is included in the price. The recipient is responsible for paying any import duties or taxes. Legal title passes to the customer on dispatch by our distributor. Single issues from current and recent volumes are available at the current single-issue price. To order the journal, please complete the form on our website. For submissions, subscriptions and all other information please visit: http://en.sv-jme.eu/. You can advertise on the inner and outer side of the back cover of the journal. The authors of the published papers are invited to send photos or pictures with short explanation for cover content. We would like to thank the reviewers who have taken part in the peerreview process. © 2015 Strojniški vestnik - Journal of Mechanical Engineering. All rights reserved. SV-JME is indexed / abstracted in: SCI-Expanded, Compendex, Inspec, ProQuest-CSA, SCOPUS, TEMA. The list of the remaining bases, in which SV-JME is indexed, is available on the website. The journal is subsidized by Slovenian Research Agency. Strojniški vestnik - Journal of Mechanical Engineering is available on http://www.sv-jme.eu, where you access also to papers' supplements, such as simulations, etc. Contents Strojniški vestnik - Journal of Mechanical Engineering volume 61, (2015), number 4 Ljubljana, April 2015 ISSN 0039-2480 Published monthly Papers Špela Bolka, Janko Slavič, Miha Boltežar: Identification of Out-of-Plane Material Characteristics through Sheet-Metal Blanking 217 Dario Croccolo, Massimiliano De Agostinis, Stefano Fini, Alfredo Liverani, Nicolò Marinelli, Eugenio Nisini, Giorgio Olmi: Mechanical Characteristics of Two Environmentally Friendly Resins Reinforced with Flax Fibers 227 Yabing Cheng, Yang Wang, Lei Li, Shuaibing Yin, Lichi An, Xiaopeng Wang: Design Method of Dual Phase Hy-Vo Silent Chain Transmission System 237 Franci Čuš, Uroš Župerl: Surface Roughness Control Simulation of Turning Processes 245 Janez Lupše, Leopold Škerget, Jure Ravnik: BEM-Based Algorithm for URANS Simulations of Flow over a Square Cylinder 254 Wang Jun, Sun Binyu, Huang Qingxue, Li Hongzhou, Han Heyong: Research on the Position-Pressure Master-Slave Control for a Rolling Shear Hydraulic Servo System 265 Florencio Sanchez-Silva, Ignacio Carvajal-Mariscal, Ariel E. Moreno-Cordobés, Pedro Quinto Diez, Miguel Toledo Velazquez: Study of an Annular Two-Phase Thermosyphon Used as an Isothermal Source in Thermometry 273 Strojniški vestnik - Journal of Mechanical Engineering 61(2015)4, 217-226 © 2015 Journal of Mechanical Engineering. All rights reserved. D0l:10.5545/sv-jme.2014.2302 Original Scientific Paper Identification of Out-of-Plane Material Characteristics through Sheet-Metal Blanking Špela Bolka1 - Janko Slavič2 - Miha Boltežar2* 1 Hidria Rotomatika, Spodnja Idrija, Slovenia 2 University of Ljubljana, Faculty of Mechanical Engineering, Slovenia The mechanical characteristics of sheet metals are typically identified in the in-plane direction, although the sheet-metal forming processes (e.g., blanking, deep-drawing) are normally applied in the out-of-plane direction. As the mechanical characteristics are not necessarily constant, their direct experimental evaluation through the forming process would enable material monitoring and process optimization, and, additionally material characterization in the out-of-plane direction. Full, partial (to a certain depth) and sequential (consecutive partial steps to full penetration) blanking experiments are performed on a laboratory blanking apparatus to correlate the out-of-plane material characteristics with the in-plane ones. The well-established in-plane approach to damage is introduced for the out-of-plane direction to determine the isotropic Lemaitre damage variable. Furthermore, yield and ultimate shear stresses are determined and correlated to their respective in-plane counterparts, offering a new insight in the sheet-metal blanking process. Keywords: cutting and forming, constitutive behavior, elastic-plastic material, mechanical testing, ductile damage, out-of-plane experiment Highlights • Experimental out-of-plane sheet-metal material characterization. • Use of in-plane approach for out-of-plane characterization. • Direct material characterization and correlation to in-plane parameters. • Experimental Lemaitre damage characterization in the out-of-plane direction. • Full, partial and sequential blanking experiments on laboratory blanking apparatus. 0 INTRODUCTION In sheet-metal forming, for instance blanking, shearing or deep drawing, identification of material mechanical properties (e.g. yield stress, ultimate stress, fracture point or material anisotropy) is crucial for manufacturing of high quality parts. The material parameters are typically determined experimentally, through in-plane tensile or shear tests, which require special preparation of specimens and cannot be applied continuously. As material characteristics may vary, a direct characterization process would enable fast correlation of current material characteristics to those of the standardized in-plane experiments and, thus, opened possibilities for process optimization. The identification of yield and ultimate stresses in the in-plane direction is well documented and standardized. On the other hand, the evaluation of fracture point in ductile metals can be approached through several mechanisms [1] and [2], leading to different experimental evaluation procedures. Among the mechanisms, roughly grouped into micromechanical models (e.g. [3] to [6]), uncoupled models (e.g. [1] and [2]) and coupled phenomenological models (e. g. [1], [7] and [8]), the Lemaitre coupled phenomenological model, following the continuum damage mechanics (CDM) may be considered as most appropriate for direct and continuous material characterization. The model considers material damage through degradation of macroscopic material properties, for instance Young's modulus of elasticity or microhardness [9] and [10]. Damage is described as gradual material weakening due to initiation, growth and propagation of microscopic cracks and voids, resulting in macroscopic crack [1] and [11]. The Lemaitre damage model, based on experimental in-plane identification of parameters or on an inverse numerical analysis, has been successfully implemented in both in-plane (e.g. [7] and [12]) and out-of-plane (e.g. [13] and [14]) numerical applications. Nevertheless, extensions of the model have been proposed, for example including damage anisotropy (e.g. [15] and [16]) or Lode angle (e.g. [1] and [17]) for low stress-triaxiality applications. Experimental identification of Lemaitre damage variable is typically performed through observation of variation in modulus of elasticity [11] and [12] or microhardness [15] and [18] by the uniaxial tensile experiment. The out-of-plane damage value is then estimated through the in-plane measurements [10]. In the out-of-plane direction, several attempts of experimental characterization have been reported by the use of shear-punch testing, mainly as a validation of a numerical analysis. In most researches, yield *Corr. Author's Address: University of Ljubljana, Faculty of Mechanical Engineering, 6 Aškerčeva cesta, Ljubljana, Slovenia, miha.boltezar@fs.uni-lj.si 217 and ultimate shear stress are identified and correlated to tensile stresses [19] and [20] and Abendroth and Kuna [21] used the shear punching process for inverse identification of the micromechanics-based GTN damage model parameters. But direct identification of the fracture and damage parameters has not yet been fully concluded [20]. Therefore, sheet-metal blanking would offer similar options to the shear-punch testing, resulting in total separation of the sheet metal, and thus enabling fracture point determination. Furthermore, the identification and optimization of the blanking process, although one of the most frequently used industrial processes for sheet laminations, is mostly focused on the optimization of the tools [22] to [24] or electromagnetic properties of the finished part [25]. Therefore would a direct material characterization procedure enable further understanding and optimization of the process. In this research, sheet-metal blanking process was used to identify the out-of-plane material characteristics. Full, partial (to a certain depth) and sequential (consecutive partial steps up to full penetration) quasi-static and dynamic blanking experiments were performed in order to obtain the yield and ultimate stress and the material-damage progression directly from an out-of-plane experiment. A novel blanking apparatus (presented in detail in [26]) offers an out-of-plane experimental analysis of the sheet metal. The identified material parameters were compared to the material parameters obtained from classic uniaxial tensile tests. The article is organized as follows. In Section 1 the theoretical background of the Lemaitre damage model and the small-strain plasticity are presented. Section 2 details the experimental approach. The experimental results are presented and discussed in Section 3. The conclusions are drawn in Section 4. 1 THEORETICAL BACKGROUND In sheet-metal blanking process, a non-linear material response is accompanied with different phenomena of material degradation. The Lemaitre phenomenological damage model has an advantage of observing the degradation from an averaged, macroscopic point of view, therefore characterization of different damage phenomena is not required. In this section, the constitutive equations of plasticity and the original formulation of the Lemaitre damage model are briefly summarized. The small-strain rate-independent plasticity with isotropic material hardening and isotropic damage is considered. The Lemaitre damage model defines the isotropic scalar damage variable D (0 < D < 1) as the ratio of the damaged area S and the virginal (undamaged) area S0 at any given plane in the volume of the material [27]: D = 1 - S. (1) Using the hypothesis of strain equivalence [28], the stress tensor ст of the undamaged material is replaced by the effective stress tensor CT in the damaged material: CJ = 1 - D (2) In the framework of small-strain rate-independent plasticity, the total strain after the yield can be decomposed into the elastic £e and plastic £p (see [29] for details): £ = £e + £p . (3) The Helmholtz free energy у is defined as the free potential, depending on the internal variables, i.e., the elastic strain £e, the internal variables associated with isotropic hardening к and the isotropic damage variable D. The free energy can be decomposed into the elastic-damage fed and the plastic free energy [28]: у(,к,D) = yei (,D) + ^p = -ee : (1 - D)D : ee -W (к)], (4) where p represents the material density and у1(к) is the free energy associated with isotropic hardening. D is the standard elasticity tensor. The thermodynamic forces associated with the internal variables (see [29]) are determined as: erđ (,D) a = p-^-= (1 - D )D : e\ R = p Bee _ d^p (к) дк =R (к). D) 1 Y = -'- = --Ee : D : se dD 2 (5a) (5b) (5c) where ст is the stress tensor, R is the thermodynamic force associated with isotropic hardening and Y is the thermodynamic force associated with damage. 0 For the Lemaitre damage model, the thermodynamic damage force, or the energy density release rate, is determined as [26]: -Y = - о e 2E (1 - D ) _ (1 + v) + 3(1 - 2v) ( \2 P oea V eq у (6) where: = V 2|s|1 (7) is the von Mises equivalent stress with s as the stress deviator tensor. p = 1tr(a) represents the hydrostatic stress, E is the Young's modulus of elasticity and v is the Poisson's ratio. The evolution of the internal variables is derived by assuming the existence of a scalar potential of the dissipation W that can be decomposed into the plastic Wp and damage Wd components (see [10] for details): (1 - D)( +1) r -Y Г1 — I , (8) where Ф is the yield function and r and s are the material parameters of the damage model. For the von Mises yielding, the yield function is given as: Ф (Pq, D ) "eq 1—D (9) where uy represents the current yield stress: °y = ^y0 + R (к), (10) with ay0 as the initial yield stress and R(k) is the isotropic hardening function, with the hardening variable к equal to the accumulated plastic strain ep (see [27] for details): (11) The evolution of the internal variables is determined by: -p ■ dw -ЖТ Ep =Y— = yN da . .dw . .p К = -Y-= Y = £p, dR D _ - 1 f-Y ~Y dY ~J 1 - D \ r (12a) (12b) (12c) where j represents the plastic multiplier and N is the plastic flow vector. Assuming an associative plastic flow it follows: -fi] N 4 2 Ы. (13) The evolution of the internal variables is subjected to the Kuhn-Tucker loading/unloading conditions: Y > 0, Ф< 0, у Ф = 0. (14) 2 EXPERIMENTAL PROCEDURE Two sets of experiments were performed - the classic uniaxial tensile experiments and the blanking experiments. 2.1 Material A low-alloyed steel sheet metal M400-50A HP of nominal thickness 0.5 mm was used for all the experiments in this investigation. All the specimens were obtained from the same part of the sheet-metal roll. The mechanical characteristics will be presented later on, as obtained from uniaxial tensile experiments. The uniaxial tensile test were performed in the steel-sheet rolling direction. 60 mm i < steel sheet rolling direction Fig. 1. Specimen used for blanking experiments with the indicated circular blanking shape Standard specimens (ISO 6892-1) for the tensile testing were obtained by blanking. The edges of the samples were additionally smoothed to reduce the edge effects in the tensile tests. For the blanking experiments, specimens, an example is shown in Fig. 1, were obtained by water-jet cutting. The central hole with a diameter of 1.2 mm, used for the punch-displacement measurement (explained later), was obtained by drilling. r i 2.2 In-Plane Experimental Setup (Tensile Experiments) Full and stepwise uniaxial tensile experiments were performed using an Instron universal hydraulic machine. Full tensile tests were performed according to the ISO 6892-1 standard. Stepwise uniaxial tensile tests were based on the same standardized procedure. For every 2 % of deformation the sample was unloaded and loaded again up to full rupture. Strain, elongation, force and stress were measured at the hydraulic head with the built-in sensors. 2.3 Out-of-Plane Experimental Setup (Blanking Experiments) Blanking experiments were performed using the custom-built laboratory blanking apparatus [24]. The apparatus enables quasi-static (blanking speeds of about 0.1 mm/s) and dynamic blanking experiments (blanking speeds depending on the thickness of the sheet metal and the circumference of the cutting tools; typically about 200 mm/s or more). The quasi-static blanking is performed by turning a fine-threaded screw in the fixation block (Fig. 2) that propels the piston and the punch forward in the cutting module, shown in Fig 4. For dynamic blanking, the blanking apparatus is elevated at one end for an angle a (Fig. 3) and the impactor is released from rest from an arbitrary position on the slope. Typically, the blanking apparatus is raised for about 5°. The initial potential energy of the impactor is transformed to the initial kinetic energy of the piston and the punch at the base of the slope. The position of the impactor on the slope can be determined up to 1 mm accurately, resulting in ±5 mm/s control of the blanking speed. Fig. 2. Custom blanking apparatus in configuration for quasi-static experiments; g represents acceleration due to gravity The blanking is performed in the cutting module (Figs. 4 and 5). The sheet metal is held in place by 8 strong magnets that are inserted into the die. The punch displacement is measured on the punch tip, a 1-mm- diameter extension from the tip, using a Keyence LG-82 displacement laser. The punch force is measured by a Kistler 9061A ring-type force transducer, mounted behind the punch. The punch acceleration is measured using a Brül&Kjsr, type 8309 accelerometer. Fig. 3. Custom blanking apparatus: configuration for dynamic experiments; g represents acceleration due to gravity Fig. 4. The cutting module Fig. 5. The cutting module in an opened position For this research, a centric circular punch-die set with a nominal diameter of 10 mm and a clearance of 25 ^m per radius was used. The set was not lubricated and all the experiments were performed at room temperature (about 23 °C). The full experiments were performed quasi-statically with a blanking speed of about 0.1 mm/s and dynamically with blanking speeds of approximately (200, 250, 300, 350) mm/s. The partial and sequential quasi-static experiments were performed with a blanking speed of about 0.1 mm/s, whereas the dynamic blanking experiments were performed with a blanking speed just high enough for the punch to reach the designated depth. The partial experiments were performed up to depths of (100, 200, 300, 350) ^m and the sequential experiments were performed in steps of 50 ^m. The repeatability of the laboratory blanking apparatus is above 95 % and is discussed in detail in [24]. In all the out-of-plane experiments, at least five trials (to obtain minimal statistical requirements) were performed and averaged. Data was smoothed using a Hann window over 11 measurement points. 3 RESULTS 3.1 In-Plane Experimental Results 3.1.1 Yield Stress and Material Hardening Law The material hardening data was obtained from standard uniaxial tensile tests (results shown in Fig. 5). As the material exhibits a Lüder plateau, a combined isotropic hardening material model was defined. After the initial yield stress &y0 the material exhibits a linear hardening behavior with a hardening coefficient K: where E and E0 are the current and initial Young's modulus of elasticity. Under the assumption of 5 = 1, the damage parameter r was obtained from the sequential tensile tests [12]. For the uniaxial tensile test, Eq. (6) is simplified to: -7 = —-72 . (18) 2E (1 - D)2 Using Eq. (18) in Eq. (12c), damage propagation in relation to the accumulated plastic deformation can be written as: dD а d£p 2Er (1 - D )2 (19) CTy =CTyo + K£p, (15) Fig. 6. Averaged full uniaxial tensile experiments with modelled hardening curves where e p is the accumulated plastic strain. After the limit accumulated plastic strain e p0 is reached, the material exhibits an exponential Voce/saturation-type hardening law [28]: ay =ayi )( 1 - e" (16) where сту1 = oy0 + Kep0, is the saturation stress. The parameters of the combined material model for the rolling direction of steel are given in Table 1. The fitting of the combined model to the experimental data is shown in Fig. 5. 3.1.2 Damage Parameters As shown in [12], it follows from Eqs. (1), (2) and (5a): d = i - E, (17) It follows that the damage parameter r is equal to: dD 2 ЛёР ' 2E (1 - D )2 d£ (20) For statistical reasons as mentioned above, five sequential tests were performed, as shown in Fig. 6. As recommended by [10], the elastic modulus was determined on the unloading part of each step. Fig. 7 shows the degradation of the elastic modulus and the accumulation of damage in relation to the accumulated plastic strain. For each unloading point the minimum and maximum values of the parameters are shown by the error bar. Using Eq. (20), it follows r = 0.953 MPa. 3.2 Out-of-Plane Experimental Results Full, partial and sequential experiments were performed quasi-statically and dynamically. The 2 а r = mechanics of the blanking experiments is different to the mechanics of the in-plane experiments. In contrast to the clamped specimen in the tensile test, in the blanking experiment the cutting tools and the specimen initially come into contact and the cutting tools affect the material's response in the initial elastic zone. Therefore, the resulting blanking force vs. punch penetration (P - d) diagrams do not show all the distinct characteristics of the total stress vs. total strain diagrams resulting from the tensile experiment. The elastic response of the material was therefore identified during the unloading phase when the cutting tools are retracting from the material. 0.05 0.10 0.15 0.20 0.25 0.30 e [mm/mm] Fig. 7. Sequential uniaxial tensile experiments, compared to averaged full uniaxial tensile experiment Fig. 8. Identification of damage for uniaxial tensile experiments Although the inertial forces are accounted for (see [24] for details), complete elimination is not possible at higher blanking speeds, as the dynamics of the process becomes more pronounced (Fig. 8). However, during the unloading part, a similar slope of the P - d diagram is expected. 3.2.1 Yield and Ultimate Stress The shear stress can be estimated as (see [19]): P t = - 2nW ' (21) where ravg = (rdie+rpunch) / 2 is the average radius of the cutting tools and t is the thickness of the specimen. The correlation between the ultimate and the yield stress in the tensile and shear directions can be obtained from the P - d diagram for the full blanking experiments. Averaged diagrams (at least five measurements for each measurement type) for different blanking speeds are shown in Fig. 8 and the effect of the blanking speed on the shear ultimate tu and the yield ту stress during blanking in comparison to the tensile ultimate CTu and yield oy stress is shown in Fig. 9. Yielding was determined by the use of partial blanking experiments (discussed later) and a linear approximation. It was determined from partial blanking experiments that the plastic deformation of the specimens occurred at an approximate punch penetration of 30 дт. A linear approximation was used on the P - d data to acquire the blanking force at the yield. The quasi-static ultimate shear stress is lower than the ultimate tensile stress (ratio about 1.06), whereas the dynamic ultimate shear stress is slightly higher than the ultimate tensile stress (ratio of approximately 0.99). In contrast, the yield shear stress is lower than the tensile yield stress (1.45), but it increases with the blanking speed (from about 1.30 at a blanking speed of 200 mm/s, and to 1.15 at a blanking speed of 350 mm/s). Increase in process speed may lead to thermal softening of the material due to the heat generated inside the sheet-metal. The temperature increase during blanking was discussed in the literature from experimental and numerical point of view (see [24]), but for the material, presented in this study, no significant effects of generated heat were observed at different blanking speeds. Therefore the thermal effects will be disregarded in this study. Table 1. Material hardening parameters in the rolling direction, obtained from tensile tests Parameter tolling [GPa] CTy0 [MPa] sp° [mm/mm] K [MPa] [MPa] m Value 175.6 273 0.0168 321.698 396.623 20.704 3.2.2 Damage Identification For damage identification, partial (Fig. 10) and sequential (Fig. 11) quasi-static and dynamic blanking experiments were performed. As can be seen from Fig. 10, the partial dynamic experiments conform well with the full dynamic experiments. By lowering the initial kinetic energy of the punch for the stepwise experiments, the stepwise steps show a more quasi-static-like curve, but still follow closely the full experiments (Fig. 11). 7 65- CL, 3 21-1 °0 100 200 300 400 500 600 d [/mi] Fig. 9. Averaged blanking diagrams at different blanking speeds Fig. 10. Correlation of uniaxial tensile and blanking yield and ultimate stresses The material degradation was determined, similar to the uniaxial tensile tests, by observing the gradient G=dP / dd during the unloading phase of the blanking for the sequential steps (Fig. 12). In contrast to the tensile experiments, where the unloading slope of the curve is approximately linear, the blanking partial and sequential unloading curves exhibit non-linearities that are attributed to the contact conditions between the cutting tools and the specimen. Fig. 13 shows the identified unloading gradients for partial and sequential experiments at different punch penetration depths. The gradients for the unloading data increase approximately to the maximum blanking force and then start to decrease rapidly. d [/im] Fig. 11. Averaged partial experiments d [firn] Fig. 12. Averaged sequential dynamic blanking experiments For an approximately equal punch-penetration depth, similar gradient values are obtained from the quasi-static and the dynamic partial and sequential experiments. Only the small, 50 pm steps reveal the peak gradient at the maximum force depth. The P - d diagram up to the maximum force is usually assumed to show the hardening of the sheet-metal, followed by damage initiation and shearing after the force peak [29]. The shearing zone ends with an abrupt fracture of the sheet metal. By investigating metallographic samples, the fracture point was determined to be at approximately 350 дт. This can also be seen from the sequential experiments - even for sequential steps with a nominal depth of 50 pm, the final fracture occurred at a depth of approximately 350 дm. A comparison of the damage accumulation during the uniaxial tensile experiment and of the damage-like - 0.1 mm/s --- 200 m m/s 250 mm/s — 300 mm/s - 350 mm/s variable during the blanking experiment is shown in Fig 14. Fig. 13. Identification of damage from sequential blanking experiments Fig. 14. Damage at different punch penetrations for sequential and partial experiments As the strains during the blanking cannot be directly characterized, neither as pure tensile nor as pure shear strains [19] and [29], relative plastic displacement is used to compare the values of the damage variable during tensile and blanking experiments instead of the accumulated plastic strain. The relative plastic displacement for the tensile experiment drel,t is determined as: d - rel,t _ , , d (22) f,t where dyt is the displacement after the material yield and df,t is the displacement at fracture for the tensile experiment. The relative plastic displacement for the blanking experiment is determined as: dP where dp,b is the displacement after the peak force and t is the sheet-metal thickness. The damage coming from the blanking is characterized as: Я = 1 - G, b Go (24) drel b — P,b t ' (23) G0 represents the unloading gradient at the highest peak. A further analysis was performed, the same as with the data from the tensile experiments. In the tensile experiments the Young's modulus was reduced from an initial 177 MPa to a final 94 MPa, or, about 46 %. Similarly, in the blanking experiments the gradient was reduced from the steepest 376 kN/mm to 160 kN/mm, which is about 57 %. This is consistent with the results obtained by [17]. However, the observation of the relative plastic displacement shows that, while during the uniaxial tensile experiment, damage accumulation starts almost immediately with the plastification, the accumulation of damage in the out-of-plane direction starts after the maximum force is reached, resulting in a smaller relative plastic displacement at fracture. 4 CONCLUSIONS In this research, an experimental procedure to directly identify material parameters in the out-of-plane direction is proposed. The in-plane procedures from uniaxial tensile testing are implemented in the out-of-plane blanking experiments to determine shear yield and ultimate stresses and material damage accumulation up to fracture point. The out-of-plane values are compared to the values obtained from the uniaxial tensile experiment. Full, partial and sequential blanking experiments at various blanking speeds were performed to obtain the ultimate and the yield stress and the damage accumulation in relation to the punch depth in the out-of-plane direction. The damage was identified from the blanking force vs. punch penetration unloading gradient of the partial and sequential steps. The damage initiation in the out-of-plane direction can only be observed with sufficiently small sequential steps, as can be seen from the results above. Furthermore, out-of-plane damage dependency on the step size and/or process speed needs to be assessed for the tested material using partial and sequential experiments with different step sizes. Similarly, the yield and ultimate shear stresses should be tested at different process speeds in order to perform full correlation with the in-plane data. Furthermore, as the yield stress evaluation in the out-of-plane direction is influenced by the material-punch contact conditions, additional care should be taken in order to ensure equal contact conditions (e.g. lubrication, tool wear, etc.), otherwise the parameter might monitor the variation of the contact conditions instead. Fig. 15. Correlation of tensile and blanking damage propagation up to fracture With the integration of the proposed procedure into an industrial blanking press, an on-line monitoring of the material quality could be obtained from which process monitoring and optimization could be achieved, for example by correlating the blanking speed to the material fracture point. Furthermore, with appropriate extensions to the Lemaitre damage model formulation, damage parameters for a fully anisotropic Lemaitre damage model could be determined. 5 ACKNOWLEDGEMENTS This operation was partially financed by the European Union, European Social Fund (Grant nr. MR-10/103). 6 REFERENCES [1] Cao, T.-S., Gachet, J.-M., Montmitonnet, P., Bouchard, P.-O. (2014). A Lode-dependent enhanced Lemaitre model for ductile fracture prediction at low stress triaxiality. Engineering Fracture Mechanics, vols. 124-125, p. 80-96, D0l:10.1016/j. engfracmech.2014.03.021. [2] Li, H., Fu, M., Lu, J., Yang, H. (2011). Ductile fracture: Experiments and computations. International Journal of Plasticity, vol. 27, no. 2, p. 147-180, DOI:10.1016/j. ijplas.2010.04.001. [3] Gurson, A.L. (1977). Continuum theory of ductile rupture by void nucleation and growth. Part I: yield criteria and flow rules for porous ductile media. Journal of Engineering Materials and Technology, vol. 99, no. 1, p. 2-15, DOI:10.1115/1.3443401. [4] Tvergaard, V., Needleman, A. (1984). Analysis of the cup-cone fracture in a round tensile bar. Acta Metallurgica, vol. 32, no. 1, p. 157-169, DOI:10.1016/0001-6160(84)90213-X. [5] Nahshon, K., Hutchinson, J. (2008). Modification of the Gurson model for shear failure. European Journal of Mechanics - A/Solids, vol. 27, no. 1, p. 1-17, DOI:10.1016/j. euromechsol.2007.08.002. [6] Starman, B., Vrh, M., Halilovič, M., Štok, B. (2014). Advanced modelling of sheet metal forming considering anisotropy and Young's modulus evolution. Strojniški vestnik - Journal of Mechanical Engineering, vol. 60, no. 2, p. 84-92, DOI:10.5545/sv-jme.2013.1349. [7] Saanouni, K. (2008). On numerical prediction of ductile fracture in metal forming. Engineering Fracture Mechanics, vol. 75, no. 11, p. 3545-3559, DOI:10.1016/j. engfracmech.2007.02.023. [8] Lemaitre, J. (1985). A continuous damage mechanics model for ductile fracture. Journal of Engineering Materials and Technology vol. 107, p. 83-89, DOI:10.1115/1.3225775. [9] Ghozzi, Y., Labergere, C., Saanouni, K. Parrico, A. (2014). Modelling and numerical simulation of thick sheet double slitting process using continuum damage mechanics. International Journal of Damage Mechanics, vol. 23, no. 8, pp. 1150-1167, DOI:10.1177/1056789514520799. [10] Lemaitre, J., and Desmorat, R. (2005). Engineering Damage Mechanics, Springer, Berlin, Heidelberg. [11] Celentano, D.J., Chaboche, J.L. (2007). Experimental and numerical characterization of damage evolution in steels. International Journal of Plasticity, vol. 23, no. 10-11, p. 17391762, DOI:10.1016/j.ijplas.2007.03.008. [12] Mashayekhi, M., Ziad-Rad, S., Parvizian, J., Niklewicz, J., Hadavinia, H. (2007). Ductile crack growth based on damage criterion: Experimental and numerical studies. Mechanics of Materials, vol. 39, no. 7, p. 623-636, DOI:10.1016/j. mechmat.2006.10.004. [13] Hambli, R. (2001). Comparison between Lemaitre and Gurson damage models in crack growth simulation during blanking process. International Journal of Mechanical Sciences, vol. 43, no. 12, p. 2769-2790, DOI:10.1016/S0020-7403(01)00070-4. [14] Teixeira, P., Santos, A.D., Pires, F.A., de Sa, J.M.A.C. (2006). Finite element prediction of ductile fracture in sheet metal forming processes. Journal of Materials Processing Technology vol. 177, no. 1-3, p. 278-281, DOI:10.1016/j. jmatprotec.2006.04.059. [15] Ganjiani, M. (2013). Identification of damage parameters and plastic properties of an anisotropic damage model by microhardness measurement. International Journal of Damage Mechanics, vol. 22, no. 8, p. 1089-1108, DOI:10.1177/1056789513482598. [16] Desmorat, R., Otin, S. (2008). Cross-identification isotropic/ anisotropic damage and application to anisothermal structural failure. Engineering Fracture Mechanics, vol. 75, no. 11, p. 3446-3463, DOI:10.1016/j.engfracmech.2007.05.011. [17] Malcher, L., Mamiya, E. (2014). An improved damage evolution law based on continuum damage mechanics and its dependence on both stress triaxiality and the third invariant. International Journal of Plasticity, vol. 56, p. 232-261, DOI:10.1016/j.ijplas.2014.01.002. [18] Mkaddem, A., Gassara, F., Hambli, R. (2006). A new procedure using the microhardness technique for sheet metal damage characterization. Journal of Materials Processing Technology, vol. 178, no. 1-3, p. 111-118, D0l:10.1016/j. jmatprotec.2006.02.018. [19] Guduru, R.K., Darling, K.A., Kishore, R., Scattergood, R.O., Koch, C.C., Murty, K.L. (2005). Evaluation of mechanical properties using shear-punch testing. Materials Science and Engineering: A, vol. 395, no. 1-2, p. 307-314, DOI:10.1016/j. msea.2004.12.048. [20] Sellamuthu, P., Collins, P.K., Hodgson, P.D., Stanford, N. (2013). Correlation of tensile test properties with those predicted by the shear punch test. Materials & Design, vol. 47, p. 258-266, DOI:10.1016/j.matdes.2012.11.057. [21] Abendroth, M., Kuna, M. (2006). Identification of ductile damage and fracture parameters from the small punch test using neural networks. Engineering Fracture Mechanics, vol. 73, no. 6, p. 710-725, DOI:10.1016/j. engfracmech.2005.10.007. [22] Mole, N., Cafuta, G., Štok, B. (2013). A method for optimal blank shape determination in sheet metal forming based on numerical simulations. Strojniški vestnik - Journal of Mechanical Engineering, vol. 59, no. 4, p. 237-250, DOI:10.5545/sv-jme.2012.989. [23] Vrh, M., Halilovič, M., Štok, B. (2008). Impact of Young's modulus degradation on springback calculation in steel sheet drawing. Strojniški vestnik - Journal of Mechanical Engineering, vol. 54, no. 4, p. 288-296. [24] Gantar, G., Sterzing, A. (2008). Robust design of forming processes, Strojniški vestnik - Journal of Mechanical Engineering, vol. 54, no. 4, p. 249-257. [25] Vegelj, D., Zajec, B., Gregorčič, P., Možina, J. (2014). Adaptive pulsed-laser welding of electrical laminations. Strojniški vestnik - Journal of Mechanical Engineering, vol. 60, no. 2, p. 106-114, DOI:10.5545/sv-jme.2013.1407. [26] Slavič, J., Bolka, Š., Bratuš, V., Boltežar, M. (2014). A novel laboratory blanking apparatus for the experimental identification of blanking parameters. Journal of Materials Processing Technology, vol. 214, no. 2, p. 507-513, DOI:10.1016/j.jmatprotec.2013.10.006. [27] Kachanov, L. (1958). On creep rupture time. Proc Acad Sci USSR DIv Engn Sci, p. 26-31, [28] Lemaitre, J. (1992). A Course on Damage Mechanics. Springer, Berlin, Heidelberg, DOI:10.1007/978-3-662-02761-5. [29] de Souza Neto, E.A., Perić, D., Owen, D.R.J. (2001). Computational Methods for Plasticity: Theory and Applications. John Wiley & Sons, Chichester, DOI:10.1002/9780470694626. [30] Voce, E. (1948). A Practical Strain-Hardening Function. Acta Metallurgia, vol. 51, p. 219-226. [31] Klingenberg, W., Singh, U.P. (2005). Comparison of two analytical models of blanking and proposal of a new model. International Journal of Machine Tools and Manufacture, vol. 45, no. 4-5, p. 519-527, DOI:10.1016/j. ijmachtools.2004.09.002. Strojniški vestnik - Journal of Mechanical Engineering 61(2015)4, 227-236 © 2015 Journal of Mechanical Engineering. All rights reserved. D0l:10.5545/sv-jme.2014.2248 Original Scientific Paper Mechanical Characteristics of Two Environmentally Friendly Resins Reinforced with Flax Fibers Dario Croccolo* - Massimiliano De Agostinis - Stefano Fini - Alfredo Liverani -Nicolò Marinelli - Eugenio Nisini - Giorgio Olmi University of Bologna, Department of Industrial Engineering, Italy In this paper, a comparison between the mechanical characteristics of two different natural fiber composites has been carried out. The two composites were made using two different environmentally friendly resins, but with the same flax fiber reinforcement. Thirty-two specimens were tested for tension, bending and shear. The results obtained in the tensile and bending tests are acceptable according to the Standards, whereas those of the shear tests were discarded, since the samples exhibited unacceptable failure modes. In light of the results obtained, both in terms of strength and stiffness, the isophthalic resin (181EN2X) outperformed the vinyl ester resin (VEef220ST). In the case of the isophthalic resin, microscopic observations showed a good adhesion between the matrix and the fiber, with a small amount of air inclusions. The present study demonstrates that the mechanical characteristics of flax reinforced and environmentally friendly resins are similar to those achievable with the same fibers using conventional resins. Keywords: natural fibre composites, mechanical properties, flax, environmentally friendly resin Highlights • Two different natural fiber composites have been compared. • The two resins have wide applications in the marine industry. • The two composites have been obtained with two different environmentally friendly resins. • The mechanical properties of the manufactured composites have been determined. • The mechanical properties of the Resin 181EN2X are significantly improved by the presence of the fiber reinforcement. 0 INTRODUCTION Over the past twenty years, polymer matrix composites reinforced with vegetable fibers have attracted great interest [1] and [2]. The increasing attention on these materials arises from both the excellent mechanical properties of vegetable fibers, such as high tensile strength for low density, low environmental impact, renewability, recyclability and low cost. Despite the good qualities of some types of vegetable fibers, such as flax, hemp and jute, vegetable fiber reinforced plastics are still mainly used for non-structural applications only. This is due to the low strength properties that can be achieved for the composite, despite its generally good impact resistance and low density. The mechanical response of a composite is strictly related to its structure [3] to [7]. The relatively poor mechanical properties of the vegetal fiber reinforced composites are usually due to the following: i) damage to the plant fibers during the extraction process, ii) variability of fiber characteristics depending on the source. Other drawbacks are the difficulty in finding fabrics specially made for composites, resulting in non-optimal orientation and volume fractioning of the fibers that can be obtained, and a poor affinity between the fiber and the matrix. Moreover, due to the short length of technical plant fibers, the reinforcement needs to be in the form of staple fiber yarns, which have a twisted structure. Although twisting facilitates yarn processing, it has several detrimental effects on the composites produced [8]. Natural fibers also have limited thermal stability and harsh environments may seriously affect their mechanical properties. There is, finally, a major drawback related to the application of natural fibers for the reinforcement of polymeric matrices: natural fibers exhibit high moisture absorption [9] due to the presence of hydroxyl and other polar groups in their constituents. In the present study, the mechanical characteristics of composites obtained using environmentally friendly polymer matrices, reinforced with flax fibers derived from the textile industry, were evaluated. The aim of the research is to assess the potential of low cost vegetable fiber reinforced plastics. An set of experiments was carried out in order to show whether the combination of flax fiber reinforcement and environmentally friendly resins could lead to high-performance materials, with comparable structural properties to other composites made using natural fibers and non-environmentally friendly polymer matrices. *Corr. Author's Address: University of Bologna, Department of Industrial Engineering, Viale del Risorgimento 2, 40136 Bologna, Italy, dario.croccolo@unibo.it 227 1 STATE OF THE ART: VEGETABLE FIBERS AND RESINS 1.1 Vegetable Fibers There are six basic types of natural fibers. They are classified as follows: bast fibers (jute, flax, hemp, ramie and kenaf), leaf fibers (abaca, sisal and pineapple), seed fibers (coir, cotton and kapok), core fibers (kenaf, hemp and jute), grass and reed fibers (wheat, corn and rice) and all other types (wood and roots) [10]. The best fibers are bast fibers and they are extracted from the stems of plants such as flax, jute and hemp. These fibers are composed of clusters made up of long elementary fibers. The elementary fibers have varying lengths depending on the kind of plant, ranging from 20 mm to 70 mm for flax against 2 mm to 3 mm for softwood. They have extremely thin cell walls (5 дт up to 15 дт) and a diameter between 15 ^ and 35 дm [11]. These clusters of elementary fibers, whose cross section consists of 12 to 48 cells, are called technical fibers. Technical fibers are the basic element for both the textile industry and the preparation of composites. For example, flax technical fibers have lengths in the range 0.3 m to 0.6 m and diameters falling between 50 дт and 500 дm. As shown in Fig. 1, inside each technical fiber, the elementary fibers are joined together by a matrix (middle lamella) primarily consisting of pectin and hemicellulose (structures composed of low molecular weight polysaccharides). Typically, elementary fibers display a 5- to 7-edge polyhedral section, which improves the packaging capability with respect to a circular section. They consist mainly of cellulose, lignin and non-structural compounds such as waxes, salts and organic nitrogen compounds [11]. The structural characteristics of elementary fibers are mainly due to cellulose, which has a high molecular weight and a Young's modulus of E = 135 GPa [13]. pectin and hemicellulose interphase (middle lamella) Elementary fibres MM/' Fig. 1. Cross section of a technical fiber (bundle of elementary fibers) Inside the cell wall, cellulose can be mainly found in the form of bundles of cellulose polymers, joined together by means of polysaccharides (hemicellulose) and lignin. The elementary fibers have an organized structure of multiple layers, with a very thin outer skin (the primary cell wall, for flax it is about 2 ^m thick) and an intermediate zone (secondary cell wall), which is the main bearing element split into four parts. These are named S1 to S3, along with a hollow central part called the lumen (Fig. 2). The so-called S2 part is the most important in terms of volume amount, as it occupies about 80 % of the fiber volume. The secondary cell wall consists of cellulose microfibers spiral-wound with a certain tilt angle with respect to the elementary fiber axis. The secondary cell wall is included in a matrix of hemicellulose and lignin. In the S1 and S3 areas, microfibers are wound with both right hand and left hand helixes, hence forming a lattice structure. Conversely, in the S2 area, they are wound in one direction only and with a reduced helix angle in comparison with other areas. Referring to the S2 area, fibers derived from wood exhibit values for the tilt angle between 10° and 30°, whereas bast fibers have values lower than 10°, sometimes dropping down to as low as 2°. Lumen Fig. 2. Structure of a bast elementary fiber It has been found that the mechanical properties of the fibers are related to the tilt angle in the S2 area. Small values of the tilt angle are generally associated with high values of tensile strength and of the elastic modulus of the elementary fiber. Therefore, the architecture of the S2 area is responsible for the tensile strength and the stiffness of the elementary fiber, the S1 area enhances the stability of the fiber against the compressive loads, whereas the S3 area reacts to the internal pressure of the cell [14]. The density of the cell wall is about 1500 kg-m-3 and appears to be almost independent of the plant species the fiber is extracted from. Conversely, the density of different plant fibers is mainly influenced by the extension of the hollow central area (lumen) and, consequently, by the plant age. For instance, for fibers extracted from a mature plant, the lumen area is less than 10 % of the whole section of the fiber itself. Therefore, the average density of elementary fibers is quite low, on the order of 1350 kg-m-3. As for the mechanical properties of the technical fibers, significant discrepancies can be found between the values reported by different authors [11] and [15] to [17]. This occurrence is due to one or more of the following reasons [18]: (a) cultivation of the fiber, growth conditions of the plant, position of the fiber in the plant, maturity of the fiber and collection technique; (b) machining and extraction technique or mechanical modification of the fiber [8]; (c) testing methodologies and measurements performed according to different standards, and different environmental conditions during the tests. Table 1 summarizes the main physical and mechanical properties found in the literature [18] to [20] for some natural fibers. The large scattering in the tensile strength values may depend on the clamping length. In fact, for clamping lengths greater than 25 mm, values of tensile strength of approximately UTS = 500 MPa can be found, whereas for lower clamping lengths (lower than 3 mm) the tensile strength may increase up to UTS = 800 MPa to 1000 MPa. Such behavior is due to the transition from a cohesive failure of the interface made of pectin adhesive, which joins the elementary fibers together, to the rupture of the elementary fibers themselves. Moreover, when the clamping length increases, more defects are found in the specimen [16]. 1.2 Resins The resins that are used in fibre-reinforced composites are composed of long chain-like molecules consisting of many simple repeating units. Resins can be classified under two types, 'thermoplastic' and 'thermosetting', according to the effect of heat on their properties. Thermoplastics are softened by heating and harden again after cooling. Widely used thermoplastics include nylon, polypropylene and ABS, which can be reinforced, usually by short, chopped fibres, such as glass. Thermosetting materials, or 'thermosets', are formed from an in situ chemical reaction, where the resin and hardener or the resin and catalyst are mixed together and undergo a nonreversible chemical reaction to form a hard, infusible product. Once cured, thermosets do not become liquid again if heated, although their mechanical properties experience a significant change above a certain temperature. This temperature is known as the glass transition temperature (Tg) and depends on the particular resin system used, on its degree of cure and on whether it was mixed correctly. Although there are many different types of thermosetting resins being used in the composite industry, the structural parts are mainly made of polyester, vinylester and epoxy. Most polyester resins are of the 'unsaturated' type. Unsaturated polyester resin is a thermoset, capable of being cured from a liquid or solid state under suitable conditions. Orthophthalic and isophthalic are the main types of polyester resins. The orthophthalic polyester resin is the most conventional and economical. The isophthalic resin is now becoming the preferred material in marine industries, due to its higher water resistance. The addition of styrene enables the resin to cure from a liquid to a solid by 'cross-linking' the molecular chains of the polyester. These resins can therefore be moulded without the use of pressure and are called 'contact' or 'low pressure' resins. Vinylester resins have a similar molecular structure, but different locations of their reactive sites, posed only at the ends of the molecular chains. An important issue is related to the molecular chain being able to absorb shock loads over its whole length. As a consequence, vinylester resins are made tougher and more resilient than polyesters. The vinylester molecule also features fewer ester groups. These ester groups commonly have a detrimental effect, as they suffer water degradation by hydrolysis. Therefore, a further important feature of vinylesters is their better resistance to water and many other chemicals than their polyester counterparts. Therefore, they are widely used in applications involving liquids, such as pipelines and chemical storage tanks. Epoxy resins are the resins exhibiting the best mechanical properties and a good resistance to environmental degradation, with many applications in aircraft construction. When used as laminating resins, their enhanced adhesive properties and resistance to water degradation, make them suitable to many applications in the marine industry. 2 MATERIALS AND METHODS As mentioned in the introduction, the present study focuses on flax reinforced composites. In fact, flax is characterized by good mechanical and physical properties, such as tensile strength, elastic modulus, and density (Table 1), moreover its technical fibers are available in a wide variety of tissues. We chose a flax warp with a weight of 478 g-m-2, whose properties are reported in Table 2 and provided by the Ta-Bru company. Since lamination was performed by Table 1. Physical and mechanical properties of some natural fibers _Diameter pum]_Length [mm]_Density [kg-m-3] sisal_20_1A8_1330 jute_20_2J0.5_1450 flax 13 to 22 20 to 70 1500 hand, reduced environmental impact resins suitable for hand lamination were selected from among the class of thermosetting polymers. Thermosetting polymers have, in fact, good mechanical strength, chemical resistance, thermal stability and durability in comparison with thermoplastic resins. The latter, in spite of their better flexibility in terms of producing even complex shapes, have to reach process temperatures higher than the degradation threshold for flax (about 230 °C). In the set of thermosetting polymers available on the market and responding to the needs of low environmental impact, the authors have chosen two resins produced by the Polynt company. The first (Resin 181EN2X) is an isophthalic with 25% polyethylene terephthalate (PET), obtained from recycled material. The second (Resin VEef220ST) is a vinyl ester obtained from a bisphenolic epoxy resin with a low styrene content. Unlike the epoxy resins these two types of resins do not need thermal treatment after cure. In fact, thermal treatment may damage the flax fibres. The mechanical characteristics of the two resins are reported in Table 3 [21]. Considering the data reported in Table 3, it is possible to observe some differences between the isophthalic and the vinyl ester resins. In particular, the vinyl ester resin exhibits a lower elastic modulus and a higher ultimate strength than the isophthalic one. These different properties are in agreement with the related remarks in the Subsection 1.2, where it has been pointed out that a vinyl ester resin is more resilient than an isophthalic one. In order to investigate the mechanical properties of the manufactured composites, the tensile, the shear and the bending tests have been performed, following the most relevant standards [22] to [25]. The following instrumentation has been used: Instron 8033 standing press with a 25 kN load cell, MTS 3.1 station manager software, adjustable grips for tensile testing and an ad-hoc support for three-point loading, used for both bending and shear tests, according to [24] and[25]. 3 EXPERIMENTAL 3.1 Specimen Preparation The specimens were prepared by hand, following the procedure described below: (a) spreading of a release S2 spiral angle [°] Tensile strength [MPa] Elastic modulus [GPa] 22 500 to 700 9 to 20 8_393 to 773_27_ 5 to 10 345 to 1100 28 to 100 agent on a glass plane, used as a base to facilitate the material removal at the end of the operations; (b) application of the first fabric layer, manually and uniformly wetted with the resin; (c) application of the second fabric layer, followed by manual application of the resin using a brush. The last operation was repeated until the desired thickness was reached. After wetting the last layer, a layer of peel ply was applied on the top of the stack, in order to obtain a good surface finish. Table 2. Fabric properties Units Flax fabric 8615 Fiber content [%] 100% flax Specific weight [g/m2] 478 Yarns in warp [threads/cm] 10 Warp count [-] 22 Yarns in weft [threads/cm] 5 Weft count [-] 10 Table 3. Mechanical properties of the resins 181EN2X VEef220ST Tensile strength [MPa] 60 65 Young's modulus [MPa] 4000 3500 Elongation at break [%] 2 2 Flexural strength [MPa] 98 105 Flexural modulus [MPa] 4100 3600 After that, the peel ply was covered with a micro perforated film and an aerator layer. Then, the vacuum sack and the hoses for air extraction were applied. With a suitable choice of the porosity of the micro-perforated film, the excess resin distributed during the process was finally extracted, thus achieving the desired weight ratio between the matrix and the reinforcement. In the present experimentation, the weight was equally divided between the fibers and the resin, corresponding to a 0.5 matrix-reinforcement ratio. For tensile and bending tests, specimens made up of four layers, corresponding to an overall thickness of 3.5 mm ± 0.20 mm have been manufactured, whereas nine-layer, 8 mm ± 0.20 mm thick specimens, were prepared for the shear tests. In order to obtain the shapes prescribed by the standards, the plates were machined at low speed without cutting oil [23], type 1 specimens with the dimensions shown in Fig. 3 were realized for the tensile tests. [24], method A (three-point bending) specimens with the dimensions shown in Fig. 4 were manufactured for the bending tests. In light of the difficulties encountered in manufacturing thin specimens, the specimens for the shear test were manufactured as shown in Fig. 5, following the procedure suggested by [24] for nonstandard specimens. The whole set of specimens is characterized by a principal fiber orientation parallel to the greater dimension (length). The only fibers with perpendicular orientation with respect to the aforementioned direction are those needed for laying out the warp: they are on the order of 15% of the fabric by weight. Fig. 3. Specimen for tensile test Fig. 4. Specimen for bending test Fig. 5. Specimen for interlaminar shear strength test (short-beam method) 3.2 Testing The tensile tests were performed according to [22] to [23]. Seven specimens for each resin were tested, at a room temperature (RT) of RT = 23°C, and a relative humidity (ф) ф = 50 %. The standing press was equipped with a 25 kN load cell. The nominal strain was determined according to [22], method A. The same standard was used to evaluate the Young's modulus. Since results may be significantly affected by misalignments [26], the specimens were clamped by means of self-aligning grips, with a clamping length of 115 mm ± 1 mm. The velocity of the crosshead was set at 1.8 mm/min throughout the test. The bending tests were run according to [24] standard, with a three-point bending loading configuration. The nominal span length of the specimen was L = 60 mm. Six specimens for each resin were tested, at a room temperature of RT = 23 °C and a relative humidity ф = 50 %. The standing press was equipped with a 25 kN load cell. The flexural modulus Ef was determined according to [24], method A. The crosshead velocity was set at 1mm/ min throughout the test. The shear tests were run according to [25] standard. Six specimens were tested at a room temperature of RT = 23 °C and a relative humidity ф = 50 %. The standing press was equipped with a 25 kN load cell. The crosshead velocity was set at 1 mm/min. At the end of each test, conducted under the aforementioned loading modes, the fracture surfaces were carefully examined in order to verify that the incurred failures conformed to the standard prescriptions. 3.3 Results 3.3.1 Tensile Tests The results obtained for Resin 181EN2X and Resin VEef220ST are collected in Tables 4 and 5 and Figs. 7 and 8 respectively. Some data have been discarded (Rejected: Y) and were not considered for further analyses, as the cracking did not conform to [25] requirements. 3.3.2 Bending Test Tables 6 and 7 and Figs. 9 and 10 show the bending test results for Resin 181EN2X and Resin VEef220ST, respectively. Table 4. Tensile test results for Resin 181EN2X Specimen ID ob [MPa] Et [MPa] Rejected Gbm [MPa] Std. dev. [MPa] Em [MPa] Std. dev [MPa] 1 76.91 5076 N 2 77.43 4662 N 3 79.94 5169 Y 4 77.85 4898 N 75.24 3.80 4860 207 5 65.18 4679 Y 6 68.64 4629 N 5036 N Table 5. Tensile test results for Resin VEef220ST Specimen ID ob [MPa]_Et [MPa]_Rejected_obm [MPa] Std. dev. [MPa] Etm [MPa] Std. dev [MPa] 1_6611_4698_Y _2_72.14_4649_N 3_70.48_4749_N 4_73.39_4654_N 72.27 1.03 4822 158 5_71.95_4988_N 6_73.01_4985_N 7 72.65 4907 N Table 6. Bending test results for Resin 181EN2X Specimen ID Fmax [N] ofM [MPa] Ef [MPa]_Rejected ofMm [MPa] Std. dev. [MPa] Efm [MPa] Std. dev. [MPa] 1_0.232 113.6_4166_N 2_0.235 115.1_4659_N 3 0.201 98.45 3720 N -;---тг::-;-г;- 106.52 9.42 4437 414.2 4_0.187 91.59_4614_N 5_0.222 108.7_4815_N 6 0.228 111.7 4653 N Table 7. Bending test results for Resin VEef220ST Specimen ID Fmax [N] ofM [MPa] Ef [MPa]_Rejected ofMm [MPa] Std. dev. [MPa] Efm [MPa] Std. dev.[MPa] 1_0.167 81.93_4368_N 2_0.165 80.75_4503_N 3 0.169 82.98 4301 N n„ „c c ,0 1I(, --83.05 5.73 4513 214 4_0.192 94.14_4906_N 5_0.158 77.52_4435_N 6 0.165 80.98 4570 N until an interlaminar failure mode is obtained. A comparison between the characteristics of the resins with and without reinforcement is reported in Table 8. The data referring to the resin without reinforcement are derived from [21] and are the same as those listed above in Table 3. In the case of resin 181 EN2X, the reinforcement increases the ultimate tensile strength UTS by about 25.4 %, the Young's modulus Et by about 21.5 %, the flexural strength f by about 8.7 % and the flexural modulus Ef by about 8.2 % with respect to the resin itself. In the case of resin VEef220ST, the beneficial effect of the reinforcement is a bit less evident in tension strength. The results are slightly controversial in the case of the flexural response, since reinforcement leads to an increase in the elastic modulus, but produces a reduction in strength. Nevertheless, the reinforcement has a significant effect on improving the strength of the composites especially in the case of resin 181 EN2X, which behaves better than resin VEef220ST. The results obtained for the first one are consistent with those retrieved in the literature, even if some authors obtained better properties for similar composites [17], [18], [27] and [28]. The data for Table 8 are compared Fig. 6. Specimen fractured in shear test 3.4 Discussion The data reported in Tables 4 and 5 for tensile tests and in Tables 6 and 7 for bending tests prove that the results are quite consistent with each other and repeatable. Moreover, the failure modes obtained comply with the relevant standards. Therefore, it may be concluded that both the matrix and the reinforcement contribute to the overall strength of the composite. On the other hand, results obtained in the shear test are not satisfactory due to an unacceptable failure mode of the whole set of specimens (Fig. 6). This outcome could be due to the characteristics of the laminated composite. The issue may be fixed by choosing a span length different from that suggested by the aforementioned standard, e [-] Fig. 7. Tensile test for resin 181EN2X О 0.005 0.01 0.015 0.02 0.025 S H Fig. 8. Tensile test for resin VEef220ST to each other in the histograms of Figs. 11 and 12. Confidence intervals (±2 Std. Dev.) are added to account for the scatter in the results reported in Tables 4 to 7. The values of strengths and of elastic moduli for composites, having the same feature of being reinforced by flax fiber, but manufactured by conventional epoxy resins, are added in the bar graphs of Figs. 11 and 12 (data from [28]), as reference values. It can be remarked that, thanks to the enhancement provided by the reinforcement, the strength and the stiffness of the environmentally friendly composites made using flax fiber are comparable to those of the reinforced composites made using conventional resins. Fig. 9. Bending test for resin 181EN2X 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 E H Fig. 10. Bending test for resin VEef220ST In order to assess the matrix-reinforcement adhesion and the presence of defects, some specimens manufactured with the resin 181EN2X have been sectioned and then observed with an optical microscope. The assessment showed that the warp yarns have diameters of 500 цш to 600 цш (Fig. 13), whereas the weft yarns have diameters between 200 цш and 250 цш (Fig. 14). Air inclusions, shown in Fig. 15, are fewer and smaller with respect to the available data in the literature for composites manufactured for boat hulls. The analysis of the fracture surfaces of specimens tested in tension (Fig. 16) shows a homogeneous and regular surface. Such an occurrence demonstrates that both the matrix and the reinforcement participate together in the overall Table 8. Comparison of the characteristics of the resins with and without reinforcement Resin 181EN2X [MPa] Resin 181EN2X + 50% flax [MPa] Variation [%] Resin VEef220ST [MPa] Resin VEef220ST + 50% flax [MPa] Variation [%] UTS 60 75.24 + 25.4 65 72.27 + 11.2 Et 4000 4860 + 21.5 3500 4822 + 37.7 98 106.52 + 8.7 105 83.05 -20.9 Ef 4100 4437 + 8.2 3600 4513 + 25.4 S 80 s on 60 ti ^ 40 20 I I I I I I to è! ■ Resin 181EN2X ■ Resin 181EN2X+ 5C% flax ■ Resin VEef220ST • Resin VEef220ST + 5(Ж, flax ■ Reinfcrced composite with conventional resin ■ Resin 181EN2X ■ Resin 181EN2X+ 50% flax ■ Resin VEef220ST ■ Resin VEef220ST + 50% flax ■ Reinforced composite with conventional resin Fig. 11. Comparison between a) the ultimate strengths and b) the Young's moduli of unreinforced and reinforced composites with environmentally friendly and conventional epoxy resins [28] 5000 " 4000 3000 ^ 2000 1000 b) 0 ■ Resm 181EN2X ■ Resin 181EN2X+ 50% flax ■ Resin VEef220ST • Resin VEef220ST + 50% flax ■ Reinforced composite with conventional resin ■ Resin 181EN2X ■ Resin 181EN2X+ 50% flax ■ Resin VEef220ST ■ Resin VEe£220ST + 50% flax ■ Reinforced composite with conventional resin Fig. 12. Comparison between a) the flexural strengths and b) the flexural Young's moduli of unreinforced and reinforced composites with environmentally friendly and conventional epoxy resins [28] Fig. 14. Cross section of a weft yarn Fig 15. Air inclusion in a sectioned specimen Fig. 13. Cross section of a warp yarn Fig. 16. Fracture surface of a specimen broken in tension tensile strength of the composite. In other words, a good adhesion has been achieved between the two components. From the manufacturing standpoint, a drawback of the fibers used in the present study is that the fibers absorb a high amount of resin, making it difficult to achieve a high reinforcement to the matrix ratio. 4 CONCLUSIONS The present study compares the mechanical characteristics of two different composites made of natural fiber and obtained from two different environmentally friendly resins combined with the same flax fiber reinforcement. Thirty-two specimens were manufactured by hand. They were then tested for tension, bending and shear according to the relevant international standards. The results obtained in the case of tensile and bending tests are acceptable according to the standards, whereas those coming from the shear tests have been discarded, due to unacceptable failure modes observed. Therefore, the comparison was carried out based on solely the tensile and bending test results. In light of this, the isophthalic resin (181EN2X) showed a more uniform improvement in its properties with respect to the vinylester resin (VEef220ST) in terms of strength and stiffness. In the case of the isophthalic resin, microscopic observations showed a good adhesion between the matrix and the fiber, with a small amount of air inclusions. The outcome is that the mechanical properties of the isophthalic resin are remarkably enhanced by the presence of the fiber reinforcement. In other words, this kind of resin is able to take advantage of the presence of flax fibers. Conversely, the addition of the same fibers may even have a detrimental effect with regard to the vinyl ester resin. The results of the present study show that the mechanical characteristics of composites created using environmentally friendly resins can be compared with those achievable with similar fibers but created using conventional resins. 5 NOMENCLATURE stb nominal strain at break [-]; ab stress at break [MPa]; Et slope of the stress/strain curve in the strain interval between sj = 0.05 % and s2 = 0.25 % [MPa]; abm average value of ab [MPa]; Etm average value of Et [MPa]; f flexural strength of the test specimen at maximum load [MPa]; Ef flexural modulus of elasticity [MPa]; ajMm average value of a^ [MPa]; Efm average value of Ef [MPa]; Fmax failure or maximum load [N]; UTS ultimate tensile strength of the elementary fiber [MPa]; RT room temperature [°C]; Ф relative humidity [-]; L span length [mm]; b specimen cross section width [mm]; h specimen cross section height [mm]. 6 REFERENCES [1] Hornsby, P.R., Hinrichsen, E., Tarverdi, K. (1997). Preparation and properties of polypropylene composites reinforced with wheat and flax straw fibers - Part I Fibre characterization. Journal of Materials Science, vol. 32, no. 2, p. 443-449, D0l:10.1023/A:1018521920738. [2] Hornsby, P.R., Hinrichsen, E., Tarverdi, K. (1997). Preparation and properties of polypropylene composites reinforced with wheat and flax straw fibers - Part II Analysis of composite microstructure and mechanical properties. Journal of Materials Science, vol. 32, no. 4, vol. 1009-1015, D0I:10.1023/A:1018578322498. [3] Croccolo, D., De Agostinis, M., Olmi, G. (2013). Experimental characterization and analytical modelling of the mechanical behaviour of fused deposition processed parts made of ABS-M30. Computational Materials Science, vol. 79, p. 506518, D0I:10.1016/j.commatsci.2013.06.041. [4] Croccolo, D., De Agostinis, M. (2013). Analytical solution of stress and strain distributions in press fitted orthotropic cylinders. International Journal of Mechanical Sciences, vol. 71, p. 21-29, D0I:10.1016/j.ijmecsci.2013.03.002. [5] Kopfer, H., Friedrich, C., De Agostinis, M., Croccolo, D. (2012). Friction characteristics in light weight design focusing bolted joints. Proceedings of ASME International Mechanical Engineering Congress and Exposition, Houston, vol. 3, p. 839846, D0I:10.1115/IMECE2012-85940. [6] Troiani, E., Donati, L., Molinari, G., Di Sante, R. (2014). Influence of plying strategies and trigger type on crashworthiness properties of carbon fiber laminates cured through autoclave processing. Strojniški vestnik - Journal of Mechanical Engineering, vol. 60, no. 6, p. 375-381, D0I:10.5545/sv-jme.2013.1506. [7] Movaghghar, A., Lvov, G.I. (2012) Theoretical and experimental study of fatigue strength of plain woven glass/ epoxy composite. Strojniški vestnik - Journal of Mechanical Engineering, vol. 58, no. 3, p. 175-182, DOI:10.5545/sv-jme.2011.135. [8] Shah, D.U., Schubel, P.J., Clifford, M.J., (2013). Modelling the effect of yarn twist on the tensile strength of unidirectional plant fibre yarn composites. Journal of Composite Materials, vol. 47, no. 4, p. 425-436, D0I:10.1177/0021998312440737. [9] Mukhopadhyay, S., Fangueiro, R. (2009). Physical modification of natural fibers and thermoplastic films for composites - A review. Journal of Thermoplastic Composite Materials, vol. 22, no. 2, p. 135-162, D0I:10.1177/0892705708091860. [10] Faruk, O., Bledzki, A.K., Fink, H., Sain, M. (2012). Biocomposites reinforced with natural fibers: 2000-2010. Progress in Polymer Science, vol. 37, no. 11, p. 1552-1596, D0I:10.1016/j.progpolymsci.2012.04.003. [11] Hughes, M. (2012). Defects in natural fibres: their origin, characteristics and implications for natural fibre-reinforced composites. Journal of Materials Science, vol. 47, no. 2, p. 599-609, DOI:10.1007/S10853-011-6025-3. [12] Baley, C. (2002). Analysis of the flax fibres tensile behavior and analysis of the tensile stiffness increase. Composites -Part A: Applied Science and Manufacturing, vol. 33, no. 7, p. 939-948, D0I:10.1016/S1359-835X(02)00040-4. [13] Sakurada, I., Nukushina, Y., Taisuke, I. (1962). Experimental determination of the elastic modulus of crystalline regions in oriented polymers. Journal of Polymer Science, vol. 57, no. 165, p. 651-660, D0I:10.1002/pol.1962.1205716551. [14] Booker, R.E. Sell, J. (1998). The nanostructure of the cell wall of softwoods and its functions in a living tree. Holz als Roh - und Werkstoff (European Journal of Wood and Wood Products), vol. 56, no. 1, p. 1-8, D0I:10.1007/s001070050255. [15] Symington, M.C., Banks, W.M., West, O.D., Pethrick, R.A. (2009). Tensil testing of cellulose based natural fibres for structural composite applications. Journal of Composite Materials, vol. 43, no. 9, p. 1083-1108, D0I:10.1177/0021998308097740. [16] Bos, H.L., Van Den Oever, M.J.A., Peters, O.C.J.J. (2002). Tensile and compressive properties of flax fibres for natural fibre reinforced composites. Journal of Materials Science, vol. 37, no. 8, p. 1683-1692, D0I:10.1023/A:1014925621252. [17] Zhu, J., Zhu, H., Njuguna, J., Abhyankar, H. (2013). Recent development of flax fibres and their reinforced composites based on different polymeric matrices. Materials, vol. 6, no. 11, p. 5171-5198, D0I:10.3390/ma6115171. [18] Yan, L.A., Chouw, N.A., Jayaraman, K.B. (2014). Flax fibre and its composites - A review. Composites Part B: Engineering, vol. 56, p. 296-317, D0I:10.1016/j.compositesb.2013.08.014. [19] Oksman, K., Mathew, A.P., Längström, R., Nyström, B., Joseph, K. (2009). The influence of fibre microstructure on fibre breakage and mechanical properties of natural fibre reinforced polypropylene. Composites Science and Technology, vol. 69, no. 11-12, p. 1847- 1853, D0I:10.1016/j. compscitech.2009.03.020. [20] Ku, H., Wang, H., Pattarachaiyakoop, N. Trada, M. (2011). A review on the tensile properties of natural fiber reinforced polymer composites. Composites Part B: Engineering, vol. 42, no. 4, p. 856-873, D0I:10.1016/j.compositesb.2011.01.010. [21] Polynt (2014), from http://www.polynt.it/ENG/homepage. aspx, accessed on 2014-06-10. [22] ISO 2012.527-1 (2012). Plastics-Determination of tensile properties-General principles, II ed. International Organization for Standardization, Geneva. [23] ISO 1997.527-4 (1997). Plastics-Determination of tensile properties-Test conditions for isotropic and orthotropic fobre-reinforced plastic composites, I ed. International Organization for Standardization, Geneva. [24] ISO, 1998.14125 (1998). Fibre-Renforced plastic composites-Determination of flexural properties, I ed. International Organization for Standardization, Geneva. [25] ISO 2003.14130 (2003). Technical corrigendum 1. Fibre-Renforced plastic composites-Determination of apparent interlaminar shear strength by short-beam method. International Organization for Standardization, Geneva. [26] Olmi, G. (2011). A new loading-constraining device for mechanical testing with misalignment auto-compensation. Experimental Techniques, vol. 35 no. 6, p. 61-70, D0I:10.1111/j.1747-1567.2010.00678.x. [27] Petrucci, R., Santulli, C., Puglia, D., Sarasini, F., Torre, L., Kenny, J. M. (2013). Mechanical characterisation of hybrid composite laminates based on basalt fibres in combination with flax, hemp and glass fibres manufactured by vacuum infusion. Materials and Design, vol. 49, p. 728-735, D0I:10.1016/j. matdes.2013.02.014. [28] Fiore, V., Valenza, A., Di Bella, G., (2011). Mechanical behavior of carbon/flax hybrid composites for structural applications. Journal of Composite Materials, vol. 46, no. 17, p. 2089-2096, D0I:10.1177/0021998311429884. Strojniški vestnik - Journal of Mechanical Engineering 61(2015)4, 237-244 © 2015 Journal of Mechanical Engineering. All rights reserved. D0l:10.5545/sv-jme.2014.2318 Original Scientific Paper Design Method of Dual Phase Hy-Vo Silent Chain Transmission System Yabing Cheng* - Yang Wang - Lei Li - Shuaibing Yin - Lichi An - Xiaopeng Wang Jilin University, School of Mechanical Science and Engineering, China Based on the theory that the Hy-Vo silent chain meshes with the sprocket, and the related condition that the Hy-Vo silent chain sprocket meshes with a tool, multi-variations of the Hy-Vo silent chain and sprocket are designed. The meshing design system of the dual phase Hy-Vo silent chain plate-sprocket-tool is built. The design method of the dual phase Hy-Vo silent chain transmission system is presented. Combined with a concrete example, the dual phase Hy-Vo silent chain transmission system model is established, and the dual phase transmission system is decomposed into two single-phase transmission systems. The research shows that the instantaneous transmission ratio of a dual phase transmission system may be different from the instantaneous transmission ratio of a decomposed single-phase transmission system. If the two pitches of the dual phase sprocket are equal and the two groups of pressure angles are the same, the largest fluctuation quantity of the dual phase transmission system is less in comparison with the single-phase transmission system, and the polygon effect is smaller. The design method will provide a corresponding theoretical guidance for the design of other multi-phase chain systems. Key words: dual phase transmission, Hy-Vo silent chain, multi-variation, fluctuation quantity, polygon effect Highlights • Designed the multi-variation Hy-Vo silent chain and sprocket. • Established the design system of the dual phase Hy-Vo silent chain plate-sprocket-tool. • Studied the polygon effect through mathematical analysis and simulation analysis. • Reduced the largest fluctuation quantity and the polygon effect of drive system. • Provided the theoretical guidance for the design of multi-phase chain system. 0 INTRODUCTION The Hy-Vo silent chain is a series of high-end products in the field of mechanical transmission. In recent years, with the continuous development of highspeed chain transmission technology and the growth of different personalized needs regarding automotive engines, the Hy-Vo silent chain and sprocket have undergone a continuous variation development trend; this trend represents significant change and innovation for Hy-Vo silent chain transmission technology. Research on Hy-Vo silent chain transmission systems has a long history. Meng et al. proposed the proper meshing conditions for the Hy-Vo silent chain and the sprocket [1], and the meshing system of the new silent chain plate-sprocket-hob was built; meanwhile, the design system of Hy-Vo silent chain transmission system was established [2] to [3], and the wear failure mechanism was explored through the road experiments [4] to [5]. Wang et al. designed a new sprocket tooth profile to reduce the polygonal action and meshing impact. If this method were adopted, the meshing impact and friction of the chain would be decreased [6]. Liu et al. changed the plate with a nano-structured metal mesh-polyurethane composite material, and proposed a new type of double pitch silent chain and a new silent chain assembly method [7] to [8]. Sun et al. obtained the ANSI sprocket tooth profile and involute sprocket tooth profile, and studied the effects of the dynamic tension of the chain in the silent chain drive [9]. Bucknor et al. developed a quasi-static model for a silent chain drive and modelled the kinematic effects of rocker-pin clearances [10]. Troedsson and Vedmar [11] and Stephenson et al. [12] researched the oscillations and forces of the chain drive; the model has been established to work at moderate and high speed conditions. Schwab et al. presented [13] a procedure to estimate the maximum contact force of the impact model, and Flores et al. [14] presented a general methodology to constrain rigid multibody systems. Based on these findings, Pereira et al. [15] built chain drive automatic multibody models from a minimal set of data, and aimed to overcome the difficulty of building manually complex models of chain drives. At present, the multi-variation design on Hy-Vo silent chain and sprocket includes meshing mechanism variation, shape variation and parameter variation. The change process for meshing mechanism variation is outer meshing, inner-outer compound meshing and outer meshing with inner-outer compound meshing. "Shape variation" refers to major parts structure variations of Hy-Vo silent chains, such as plate hole shape variation, shape variation about pin *Corr. Author's Address: Jilin University, School of Mechanical Science and Engineering, Renmin street 5988 , Changchun, China. chengyb@jlu.edu.cn 237 cross section, guide plate shape variation. "Parameter variation" refers to the main parameter variations of the Hy-Vo silent chain and sprocket, including the positioning offset angle variation, the distance variation of benchmark circle centre, apothem variation, big side surface curvature radius variation of rocker-pin, the variation of plate tooth profile angle, the variation of sprocket pressure angle, and the variation of tool normal tooth angle. It should be noted that regardless of what the abovementioned variation is, they are all variations in the single-phase transmission system. In this paper, the variation of the transmission system is studied based on the single-phase Hy-Vo silent chain transmission system, the dual phase Hy-Vo silent chain transmission system is designed, and the multivariation design method of the dual phase Hy-Vo silent chain plate and sprocket are proposed, and the fluctuation quantity of the dual phase Hy-Vo silent chain transmission system is researched. 1 THE MESHING DESIGN OF THE DUAL PHASE Hy-Vo SILENT CHAIN TRANSMISSION SYSTEM 1.1 The Parameter Design of the Dual Phase Hy-Vo Silent Chain Plate Supposing that A is the hole pitch, P is the basic pitch. r is the big side surface curvature radius of the rockerpin. Sm is the distance between the circle centre O of the chain plate hole and the big side surface of the rocker-pin. у is the positioning offset angle when the rocker-pin is locked in the chain plate. a is the tooth profile angle of the Hy-Vo silent chain plate. f is the benchmark apothem between the circle centre O and chain plate outside linear profile. Fig. 1 shows the schematic diagram of the dual phase Hy-Vo silent a) b) Fig. 1. The schematic diagram of the dual phase Hy-Vo silent chain in straightening c) chain in straightening. The coordinate system xoy is established on the circle centre O. O1 is the extension intersection of the adjacent benchmark circle centre lines. O2 is the contact point of two rocker-pin at the action of tension. Defining the distance of O1 - O1 as the pitch of the dual phase Hy-Vo silent chain, in straightening: P' = P = A + 2[(r + Sm )cosY-r], (1) where f1 is the equivalent apothem between the point O1 and chain plate outside linear profile: fi = f-[(r + sm )cosy- r ] cos«. (2) The inside profile overhang about the internal meshing plate is S, and the tooth profile curvature is R. A straight line goes through the inside tooth profile curvature centre and the circle centre O, and another straight line is not only perpendicular to chain plate outside linear profile but also passes through the inside tooth profile curvature centre, and the angle between the two straight lines is X. The coordinates of the inside tooth profile curvature centre is as follows: (R - f -5)sin (X + 60°) cos X (R - f-5)cos(X- y =---—— 60° (3) cos X where f = f-(p' - ^)cosa. The load is transferred through the rolling role of the two rocker-pins in the dual phase Hy-Vo silent chain transmission system. When the adjacent links rotate mutually, the contact point in the big side surface of the rocker-pin will gradually move up, as shown in Fig. 2. The hy is the distance between the point O1 and the point O2. When the relative angle of the adjacent links is 26: hy =(Sm + Г )sin (у-в) + + tan0[(Sm + Г )cos (у-в)-Г ]. (4) When hy> 0, the point O1 is below the point O2. When hy=0, the point O1 coincides with the point O2. When hy < 0, the point O1 is above the point O2. At this moment, the working chain plate pitch is: P" = A + 2 r + Sm cos(в-у)-r cos0 (5) Equivalent apothemf is changed to: (r + Sm )cos (в-у)- r Л = f - cose cosa. (6) 1.2 The Design System of the Dual Phase Hy-Vo Silent Chain Plate-Sprocket-Tool The dual phase sprocket structure (with 35 teeth, for example) is shown in Fig. 3. Because the dual phase sprocket teeth have a phase difference, each phase tooth may have a difference in pitch, module, pressure angle parameters, and the tools processing each phase tooth are different. Supposing the dual phase sprocket pitches are P11 and P12, the modules are mjj and mj2, and the pressure angles of sprockets are a11 and a12. Fig. 2. The rotation schematic of dual phase Hy-Vo silent chain Fig. 3. The schematic diagram of the dual phase sprocket The dual phase sprocket is projected onto xoy, as shown in Fig. 4. The solid line represents phase teeth whose sprocket tooth pitch is P11, and the dotted line represents phase teeth that sprocket tooth pitch is P12, and the phase difference of the sprocket tooth is o. In the process of the chain plates are meshing with the dual phase sprocket, the working chain plate is regarded as a tool of the involute tooth. When the Hy-Vo silent chain plates are surrounding the dual phase transmission sprocket in the proper meshing place, the twining round is formed, and the twining round diameter is dc . The tooth slot symmetric centreline is n -n. The tooth slot symmetric centreline is intersected with the twining round one point O1. The first tool normal tooth angle is a21, the first tool normal pitch is P21, and the first tool module is m21. According to conjugate meshing condition P11 cos a11 = P21 cos a21 and making a21 = a, the first design system of the dual phase Hy-Vo silent chain plate-sprocket-tool is established. The dual phase sprocket centre is regarded as coordinate origin, the symmetrical centreline of the dual phase sprocket tooth is regarded as the vertical axis, the horizontal line that is through the centre of sprocket and vertical to the symmetrical centreline of the dual phase sprocket tooth is regarded as the horizontal axis, and the coordinate systems is established. Supposing the Hy-Vo silent chain initial pitch is P01 (P01 = P21) and the initial apothem is f01, the meshing between the tool (P21, a21) and the sprocket (P11, a11) could be regarded as the meshing between the Hy-Vo silent chain (P01, af01) and the sprocket (P11, a11). If the dual phase sprocket tooth number is Z, and the relative rotating angle of the adjacent working plates is 26 = 2n/z. The initial apothem f01 is: foi = f -I cotП-Тз)(Р"-P21 )/4. (7) The modification coefficient x corresponds to the tool (P21, a2i) is: ncota21 nfo 01 2 2tan(n /z) P21sma21 (8) When the outside linear profile of working link is fully meshed and located, providing AC is the involute tooth profile tangent line, the sprocket base circle radius is rb, and the involute tooth profile is: [X = rb (sin(0 + 0)-0cos(0 + 0)) |y = rb (cos (0 + 0) + 0sin (0 + 0)) Г< (9) The coordinate of the point C is: P r xC = — + f2 cosa P n yC = — / tan--f2 sina 2 z (10) the coordinate of the tangent point A is: fxA = rb sina - rb0A cosa [yA = rb cosa + rb0A sina (11) in which, 0 A = yC srna- xC cosa Assuming the diameter of the column is dR, and when Z is an even number, the measuring distance of sprocket column is: m21z cosa21 , MR = —-— + dR. cosa (12) B1 When Z is an odd number, the measuring distance of sprocket column is: Mr = m21z cosa cos aB1 21 cos- 90° (13) aB1 is the involute pressure angle when the involute tooth profile is through the centre of sprocket column, 2 x tana2 mvaBl = -mva 21 m2z cosa 21 п 2z ' Assuming the first tool normal pitch P21 is equal to Pjj and the first tool normal tooth angle a21 is equal to a1b the modification coefficient x1 corresponds to the tool (P21 = Pn, a21 = Яц) can be obtained by a transformation formula which is: 2x1 tana21 = mvaB1 - mva21 - 21 m21z cosa 21 -п -<14> surface curvature radius of the rocker-pin r is 7.4 mm; the distance between the circle centre O of chain plate hole and the big side surface of the rocker-pin Sm is 0.7 mm; the positioning offset angle that rockerpin is locked in the chain plate у is 3.667°; the tooth profile angle of Hy-Vo silent chain plate a is 30°; the benchmark apothem between the circle centre O, and the chain plate outside linear profile f is 5 mm. According to the Eq. (1), P' = A + 2[(r + Sm )cos у - r] = 9.5268 mm. Making 5 = 0.10 ~0.30 mm, R = 93 mm, X = 1.6°, according to the Eq. (3), it could obtain: У1 = (R - f -S)sin (Я + 60°) cos X (R - f-5)cos(X + 60°) cos X = 78.287 = 42.33 Making P11 = P12 = P = 9.523 mm, an = a12 = 31.5° a= n / z. When Z=35, 20 = 2n / z. According to the Eq. (5), P" = A + 2 r + Sm cos(9-y)-r cos в = 9.5603 mm. According to the Eq. (7), f01 = f2 -1 cot j(P" - P21 ) / 4 = 3.966 mm. According to the Eq. (8), the modification coefficient x corresponds to the hob (P21, a21) is: яcota21 01 2 2tan(n /z) P21sina21 = -1.344, so the modification coefficient x corresponds to the hob (P22, a22). Making dR = 5.225 mm, according to the Eq. (13), Similarly, the second design system of dual phase Hy-Vo silent chain plate-sprocket -tool is established by using another dual phase sprocket pitch P12, module m12, pressure angle of sprocket a12. 2 DESIGN EXAMPLE The basic pitch of dual phase Hy-Vo silent chain P is 9.525 mm; the hole pitch A is 8.16 mm; the big side Mr = m21z cosa21 cos aB1 cos- 90° + dR = 103.239 mm, according to Eq. (14), the modification coefficient xj corresponds to the tool (P2i = Pn, a21 = an) can be got, x1 = -1.531, so the modification coefficient x1 corresponds to the tool (P22 = P12, a22 = a12). When Z=35, x=-1.319, MR = 109.456 mm, x1 = -1.523. r, b x = z 3 THE FLUCTUATION ANALYSIS OF THE DUAL PHASE TRANSMISSION SYSTEM 3.1 Fluctuation Quantity Mathematical Analysis Based on the above design example, the dual phase Hy-Vo silent chain transmission system is built, as shown in Fig. 5. The transmission system is decomposed into two single-phase Hy-Vo silent chain transmission systems. Fig. 6 shows a single-phase Hy-Vo silent chain transmission system that the top of sprocket tooth is up, while Fig. 7 shows a single-phase Hy-Vo silent chain transmission system so that the slot of the sprocket tooth is up. Fig. 5. The dual phase Hy-Vo silent chain transmission system Fig. 6. Single-phase Hy-Vo silent chain transmission system that the top of the sprocket tooth is up The instantaneous transmission ratio in Fig. 6 is i1 = m1/m 2, and the instantaneous transmission ratio in Fig. 7 is i2 = m3 /©4. When two single-phase Hy-Vo silent chain transmission systems in Figs. 6 and 7 are engaging independently and simultaneously, at any time, i1 is not necessarily equal to i2, but in the dual phase transmission system, i1 = i2. The instantaneous transmission ratio of dual phase transmission system may be different from the instantaneous transmission ratio of decomposed independently operating singlephase transmission systems, and the different instantaneous transmission ratio has a different effect on the tight-side centreline fluctuation quantity and the polygon effect; therefore, the tight-side centreline fluctuation quantity in the dual phase transmission system will be different from the fluctuation quantity in the independent operation single-phase transmission system . The Fourier series of period T=2n about two hanging tight side fluctuations is established. Because а = n / z, the Fourier series phase difference is n. It is assumed that fluctuation functions are: fl =Yj]^an C0s (kt ) + bn sin (kt )] 1 да f2 =^,\_an cos (kt + n) + bn sin (kt-n)] . (15) Fig. 7. Single-phase Hy-Vo silent chain transmission system that the slot of the sprocket tooth is up Therefore, for the round pin silent chain system, the outer meshing Hy-Vo silent chain system and the inner-outer compound meshing Hy-Vo silent chain system, f1 +f2=0, namely, the fluctuation quantity of two tight sides can be superimposed on a straight line in the dual phase Hy-Vo silent chain transmission system. For the new type Hy-Vo silent chain systems that include outer meshing Hy-Vo silent chain and the inner-outer compound meshing Hy-Vo silent chain, f1 +f ~ 0, namely, the fluctuation quantity of two tight side can be superimposed on a straight line approximately. 3.2 Fluctuation Quantity Simulation Analysis The silent chain polygon effect has an important role in the stability of chain transmission process and the reduction of vibration noise. The tight-side centreline fluctuation quantity is the important indicator of the silent chain polygon effect. Based on the above design example, the new type Hy-Vo silent chain system model of the single-phase transmission is established, as shown in Fig. 8; so does the new type Hy-Vo silent chain system model of the dual phase transmission, as shown in Fig. 9. The tight-side centreline fluctuation quantity is extracted, and the silent chain polygon effect is studied. Fig. 8. The new type Hy-Vo silent chain system model of the singlephase transmission The largest fluctuation quantity in a meshing period is extracted, Fig. 10 shows the largest fluctuation quantity about the single-phase transmission system. The horizontal axis is the time t, and the vertical axis Y is the coordinate whose link plate centroid is in the vertical direction. The parameter Y is illustrated in Figs. 6 and 7, and the maximum fluctuation quantity of the measuring link plate is 0.40 mm. Fig. 11 shows the largest fluctuation quantity about the dual phase transmission system and the maximum fluctuation quantity of the measuring link plate is 0.26 mm. In comparison with the new type Hy-Vo silent chain system of the single-phase transmission, the largest fluctuation quantity about the new type Hy-Vo silent chain system of the dual phase transmission is reduced by 35 percent, and the polygon effect of dual phase transmission system is smaller. Fig. 9. The new type Hy-Vo silent chain system model of the dual phase transmission 4 CONCLUSIONS (1) Based on the theory that the Hy-Vo silent chain plate profile is regarded as a tool profile of involute tooth, The meshing design system of Fig. 10. The largest fluctuation quantity about the single-phase transmission system 0.0134 0.0136 0.0138 0.0140 f[sj Fig. 11. The largest fluctuation quantity about the dual phase transmission system Hy-Vo silent chain-sprocket-tool is established. It is not only suitable for the design of the singlephase Hy-Vo silent chain transmission system, but also for the design of the dual phase Hy-Vo silent chain transmission system, the design method of the dual phase Hy-Vo silent chain plate-sprocket-tool is a new kind of design method. (2) If Pn ФP12 or a11 Фa12, the two design system of the dual phase Hy-Vo silent chain plate- sprocket-tool will be different. If P12 Ф P22, a12 Ф a22, the design system of the dual phase Hy-Vo silent chain plate-sprocket-tool can be built through the relevant parameters transformation. (3) When P11 = P12 or a11 = a12, the instantaneous transmission ratio of dual phase transmission system may be different from the instantaneous transmission ratio of decomposed independent operation single-phase transmission system. In comparison with the single-phase transmission system, the polygon effect of dual phase transmission system is smaller. 5 ACKNOWLEDGEMENTS The research presented in this paper was supported by National Natural Science Foundation of China (No. 51305154). 6 REFERENCES [1] Meng, F., Li, C., Cheng, Y. (2007). Proper conditions of meshing for Hy-Vo silent chain and sprocket. Chinese Journal of Mechanical Engineering, vol. 20, no. 4, p. 57-59. (in Chinese) [2] Liu, X. Ding, S. Sun, W. Liu, J. Lu, J. Xu, M. (2011). Study on join fastness of silent chain. Proceedings of International Conference on Management Science and Intelligent Control, p. 275-278. (in Chinese) [3] Meng, F., Feng, Z., Chu, Y. (2006). Meshing theory and design method of new silent chain and sprocket. Chinese Journal of Mechanical Engineering, vol. 19, no. 3, p. 425-427, D0l:10.3901/CJME.2006.03.425. (in Chinese) [4] Meng, F., Liu, X., Xu, S. (2008). Meshing analysis and design of rotundity-datum-aperture Hy-Vo silent chain. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 222, no. 7, p. 12971303, D0I:10.1243/09544062JMES865. [5] Meng, F., Wang, S., Lan, H., Xu, R., Xu, H. (2005). Wear failure mechanism and multi-impact property of automotive engine chain. Chinese Journal of Mechanical Engineering, vol. 18, no. 3, p. 401-404, DOI:10.3901/CJME.2005.03.401. (in Chinese) [6] Cheng, Y., Meng, F., Xu, H., Wu, J. (2007). Experimental study on high speed characteristics of automotive engine Oil-pump chain. Chinese Journal of Mechanical Engineering, vol. 20, no. 1, p. 36-39, D0I:10.3901/CJME.2007.01.036. (in Chinese) [7] Wang, Y., Ji, D., Zhan, K., (2013). Modified sprocket tooth profile of roller chain drives. Mechanism and Machine Theory, vol. 70, p. 380-393, D0I:10.1016/j.mechmachtheory.2013.08.006. [8] Liu, X., Wang, W., Sun, W., Wu, T., Liu, J.., Liu, J.. (2012). Design and experimental analyse of low noise double-pitch silent chain for conveyor. Procedia Engineering, vol. 29, p. 21462150, D0I:10.1016/j.proeng.2012.01.277. [9] Sun, W., Liu, X., Liu, J. (2012). Experimental investigation of sprocket tooth form effect on dynamic tension of silent chain. Research Journal of Applied Sciences, Engineering and Technology, vol. 4, no. 7, p. 846-850. [10] Bucknor, N.K., Freudenstein, F. (1994). Kinematic and static force analysis of rocker-pin jointed silent chains with involute sprockets. ASME Journal of Mechanical Design, vol. 116, no. 3, p. 842-848, D0I:10.1115/1.2919458. [11] Troedsson, I., Vedmar, L. (2001). A dynamic analysis of the oscillations in a chain drive. ASME Journal of Mechanical Design, vol. 123, no. 3, p. 395-401, D0I:10.1115/1.1374196. [12] Stephenson, R., Glennie, D., Fawcett, J.N., Hale, J.M. (2000). A method of measuring the dynamic loads in high-speed timing chains. Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, vol. 214, no. 2, p. 217-226, D0I:10.1177/095440700021400211. [13] Schwab, A.L., Meijaard, J.P., Meijers, P. (2002). A comparison of revolute joint clearance models in the dynamic analysis of rigid and elastic mechanical systems. Mechanism and Machine Theory, vol. 37, no. 9, p. 895-913, D0I:10.1016/ S0094-114X(02)00033-2. [14] Flores, P., Ambròsio, J., Claro, J.C.P., Lankarani, H.M., Koshy, C.S. (2009). Lubricated revolute joints in rigid multibody systems. Nonlinear Dynamics, vol. 56, no. 3, p. 277-295, D0I:10.1007/s11071-008-9399-2. [15] Pereira, C.M., Ambròsio, J.A., Ramalho, A.L. (2010). A methodology for the generation of planar models for multibody chain drives. Multibody System Dynamics, vol. 24, no. 3, p. 303-324, D0I:10.1007/s11044-010-9207-x. Strojniški vestnik - Journal of Mechanical Engineering 61(2015)4, 245-253 © 2015 Journal of Mechanical Engineering. All rights reserved. D0l:10.5545/sv-jme.2014.2345 Original Scientific Paper Surface Roughness Control Simulation of Turning Processes Franci Čuš - Uroš Župerl* University of Maribor, Faculty of Mechanical Engineering, Slovenia The aim of this paper is to present a surface roughness control in turning with an associated simulation block diagram. The objective of the new model based controller is to assure the desired surface roughness by adjusting the machining parameters and maintaining a constant cutting force. It modifies the feed rate on-line to keep the surface roughness constant and to make machining more efficient. The control model was developed based on simplified models of the turning process and the feed drive servo-system. The experiments were performed to find the correlation between surface roughness and cutting forces in turning and to provide functional correlation with the controllable factors. Simulation setup and results are presented to demonstrate the efficiency of the proposed control model. In terms of surface roughness fluctuations and cutting efficiency, the suggested control model is much better than a conventional CNC controller alone. Integrating the developed control model with the CNC machine controller significantly improves the quality of machined components. Keywords: machining, turning, surface roughness, model based control, simulation Highlights • A surface roughness control simulation of turning is accomplished. • A transfer function that describes the controlled plant model is carried out. • A correlation is found between surface roughness and the cutting force. • Simulation results confirm the efficiency of the control simulation model. • The system is insensitive to disturbances due to tool wear and variations of cutting depth. • Maintaining constant cutting force leads to a better (constant) quality of surface. 0 INTRODUCTION The machining quality on computerized numerically controlled (CNC) machine tools is sensitive to the machining parameters. With modern machine tools, an operator still manually adjusts controlling parameters, such as feed rate and cutting speed. The adjusted values mainly depend on operator experience and knowledge. Standard machining catalogues and commercial cutting condition prediction software [1] generally offer only recommended values, which are far from optimum. Usually, technicians set the highest feed rate possible. The limitations are the available machine power, tool strength and surface roughness tolerance. The increase of the cutting force due to the higher feed rates increase the tool and work piece deflection, machine chatter; consequently, work piece precision is reduced. Therefore, maintaining the cutting force on the tool tip at the appropriate value is a method to ensure dimensional accuracy. To solve this problem, it is necessary to modify the controlling parameters on-line [2]. Adaptive control (AC) accomplishes this. The AC of turning operations is a logical continuance of the evolution of CNC systems. The goal of adaptive control is to control the turning process via the on-line adjusting of controlling parameters [3], which are subject to machining constraints. An AC system is introduced in the cutting process by Stute and Goetz [4]. Such systems to maximize the feed rate for a specified cutting force in turning operations have been proposed in the past [5]. Many AC systems that employed on-line adjustment of the feed rate to ensure a stable and efficient turning have been proposed [6] and [7]. Thus far, no control system has been proposed to adjust cutting condition to maintain surface roughness constant in prescribed tolerances. Model reference adaptive control (MRAC) systems have been also developed for turning [8]. MRAC systems, developed from AC theory, are widely used due to their robustness and disturbance rejection capability. The basic adaptive model reference adaptive controller approach was originally investigated by Landers and Ulsoy [9]. These controllers were simulated and evaluated and physically implemented by [10]. Numerous forms of MRAC system have been developed, especially in milling [11] and [12]. Non-deterministic MRAC controllers are not always suitable for real time control [13] because the controller requires an explicit model and significant computation. The neural network-based MRAC systems [14] and [15] allow much faster development of the models and thus make these controllers practical. The main drawback of integrating MRAC systems in turning is that many algorithms generate significant oscillatory behaviour at low depths of cutting [10]. The model-based control system developed in our research overcomes this *Corr. Author's Address: University of Maribor, Faculty of Mechanical Engineering, Smetanova 17, 2000 Maribor, Slovenia, uros.zuperl@um.si 245 problem by employing off-line numerical simulations. Via simulations, the performance and stability of the control system are determined. The optimum control parameters for turning have been investigated. Via simulations, we also attempted to determine if the proposed control model has been justified for all test conditions. The proposed control model is conceived in such a way that there is no need to perform complicated modelling since the simulations depend only on easily measurable cutting force signals. Recent studies show that force signals contained the most useful information for determining the output process variables, especially surface roughness [16]. However, in many cases, the use of force sensors is not practical for retrofit applications and the spindle power signal is often used as an alternative. By analysing the cutting force characteristics, it is possible to assess the changes in the quality of surface finish [17]. Therefore, the experiments were performed to find the correlation between surface roughness and the resultant cutting force. Correlations found for various machining conditions are incorporated into simulation models. The resulting control model maintains an almost constant surface roughness throughout the cutting process regardless of the work piece material used, cutting tool and machining parameters, especially depth of cutting. In this paper, an AC system is described for turning on a lathe with a constant cutting force constraint. The detailed structure of an adaptive surface control system for turning is presented. The model based controller modifies the feed rate on-line in order to keep the surface roughness constant and to maintain the cutting efficiency. Due to the very encouraging simulation results, the control model has been implemented successfully on a commercial lathe. The experimental results from an actual machining operation will be presented in the future work. 1 STRUCTURE OF ADAPTIVE CONTROL SYSTEM The idea of this research is to implement on-line AC in conjunction with an off-line cutting condition prediction algorithm. A model-based control system is a controller that modifies its behaviour in response to change in the dynamics of the process. The structure of the proposed control system of turning is presented in Fig. 1. The system consists of a cutting conditions prediction model, a reference, a comparison block and control block, a cutting process, a CNC servo system, a feed rate override command element and a sensor that measures the output process variable. The defined constraints of all output process variables are fed into the comparison block. Then they are compared with current feed values and send to the control block. The control block algorithm adapts cutting parameters in such a way that it improves the cutting tool life and surface roughness. The objective of the proposed control is, therefore, to adjust the feed rate in order to maintain the turning resultant cutting force and surface roughness as constant. If the resultant cutting force is maintained as constant during the process of turning, then the surface finish also remains stable. In order to control the process, force information must be available to the control algorithm at every 10 ms. LabVIEW data acquisition software is used to Digital adaptation of cutting parameters Fig. 1. Structure of surface roughness control in turning provide this information. The cutting force resultant is obtained using a Kistler force dynamometer. These measured cutting force signals are used in the control model to regulate the surface roughness. The command feed rate calculated from the control model algorithm is converted to a percentage of the programmed feed rate. The control model adjusts the feed rate by assigning an override percentage to the CNC controller of the turning machine, based on a measured cutting force. The actual feed rate is the product of the feed rate override percentage f% and the recommended feed rate f. The output block determines output process variables (WB tool wear [18], T tool life, CP operation cost, F resultant cutting force). A simulation block diagram for testing the stability and fine tuning of the control parameters is developed. After simulations, the control model can be installed on real CNC turning machine tool. 1.1 Sequence of Steps for Adaptive Control of Turning Process The sequence of steps for surface roughness control of the turning process is presented below. 1. The desired value Raref is initiated; 2. Based on desired Raref, the recommended cutting conditions are determined off-line by an artificial neural network model [19]; 3. Based on desired Raref, the reference force Fref is predicted according to the reference block B1; 4. When the force Fref is known, the command value fc is determined; 5. The other recommended cutting conditions are sent to the CNC of the machine tool; 6. The measured resultant cutting forces are compared with reference value and sent to the control (block for modification of cutting parameters); 7. Actual Ra is determined based on the measured cutting force F; 8. The control block adjusts the feed rates and sends it back to the machine; 9. Steps 6 to 9 are repeated until termination of machining. The algorithm of the proposed control structure is presented in Fig. 2. The objective of control algorithm is to adjust the turning process parameters to achieve the required values of surface finish. Fig. 2. Flow chart of adaptive surface roughness control 2 CONTROL MODEL SIMULATION BLOCK DIAGRAM The block diagram of the proposed control model is presented in Fig. 3. It enables a closed loop control of cutting force and thus produces a desired surface roughness in turning. The simulator input is the desired surface roughness, and the output is the actual surface quality. The control model is developed in Matlab Simulink. It comprises of simulator of cutting, a simulator of the feed servo-drive, a reference block, an output variables block, a block for modification of cutting parameters and the models determining the mutual relations between influencing cutting Fig. 3. Simulation block diagram for the surface roughness control model quantities. These elements describe the dynamics of the control model. The objective of the block for the modification of cutting parameters is to minimize control errors and to provide the appropriate feed rate for a specific machining situation. Digital adjustment of feed rate is accomplished in such a way that the output process variable is first compared with the reference value. If the value of actual Ra exceeds the desired surface roughness Raref then the control algorithm immediately reduces the feed rate in order to reduce the produced surface roughness, while the cutting speed remains constant. If the value of actual Ra does not reach the desired surface roughness, then the control algorithm immediately increases the feed rate. The plant to be controlled includes the cutting process and the CNC servo system. The machine tool used in the study was a typical CNC lathe controlled by a CNC whose model allows feed rate override commands ranging from 0 % to 150%. The cutting process and the CNC servo system model can be represented as: Hp (s) = F ( s) = kc f(s) i 1 H s (s) = 1+1„ 1 f(s) =_ fc (s) 1 + ts ■ s ' (1) (2) time constant of the cutting process, and ts is the time constant of the CNC servo system. The feed servo-drive model simulates the machine response to changes in the desired feed rate. The objective of the feed servo-drive system is to minimize the error and provide the appropriate and actual feed rate for a specific machining situation. The step response of the CNC system used in our research is shown in Fig. 4. The pre-programmed feed rate (0.2) is overridden from 0 % to 150 %. The time to reach 0.632 of the step input is found graphically. The time constant ts is measured to be 0.13 s. Fig. 4. Step response of the x-axis feed drive system where n is the spindle speed in [rpm], kc is the specific cutting force [N/mm2], a is the depth of cut, tc is the The whole turning plant system can be represented with the transfer function: H (s) = i--^- (3) n 1 + (0.13 + tc) • s + 0.13 • tc • s2 The uncontrollable measurement noise is not included in the above transfer function due to the difficulties of describing its random nature mathematically. The transfer functions of the machine dynamics are based on the technical specifications of the GF NDM 16 CNC lathe. The predicting models of the output variables developed in this work are based on experimental results performed in research [20]. The experimental datum contains the value of the predicted (modelled) quantity and the relevant influencing parameters (cutting parameters). Statistical and mathematical modelling is used to provide functional relationships with the most important parameters and variables, i.e. feed rate, spindle speed, surface roughness and the resultant cutting force. The results of modelling are the simulation models, which can be easily included into the Simulink simulation package. These models are depicted in Fig. 3. Prior machining the required quality of the surface finish Raref is defined. The model defining the correlation between the desired surface roughness and the appurtenant resultant cutting force is labelled as M1. In order to develop the M1 model, the experiments were carried out for all combinations of the chosen parameters that significantly affect the turning process. These input parameters are cutting conditions, depth of cut, feed rate, and cutting speed. Other parameters, such as cutting tool geometry, type of raw material, etc., are kept constant for the scope of this study. A total of 1500 combinations of input parameters (depth of cut, feedrate, and spindle speed) were selected to measure cutting forces and the corresponding roughness of the machined surface. The measurements of the tool wear during machining were not carried out, despite the fact that it has a significant effect on the rise of cutting force and on the achieved surface roughness. The aim of this study is to take into account the tool condition in the prediction of cutting forces without expensive in-process tool wear measurements. Recent studies show that cutting force signals and the surface roughness values contain the most useful information for determining the tool condition (tool wear) [5]. Therefore, the influence of the tool wear in the prediction of cutting forces is captured indirectly through increased surface roughness due to changes of the cutting edge geometry. Based on the values of these 1500 experimental data, the relation between the cutting force and the surface roughness was determined according to Eq. (4). This equation was obtained by using single regression formula for steel machining at low depths of cut. Fref = -279.26 • Raref5 + 2946.9 • Raref 4 --12032 • Raref3 + 23636 • Raref 2 -- 22062 • Raref + 8250.8. (4) Eq. (4) determines the reference cutting force necessary to reach and maintain the desired surface roughness regardless of the cutting tool condition. It indirectly takes into account the condition of the tool in the calculation of cutting forces. The output variables block (B4) tests if the simulated Ra corresponds with the desired Ra. The produced surface roughness was calculated from the above mentioned experimental data according to Eq. (5): Ra = 0.2401- e°-0034'f. (5) Similar to Eq. (4), Eq. (5) is obtained by using an exponential regression formula. Both equations are based on the same experimental data. The outputs of block B4 were also compared to the predictions of the artificial neural network surface roughness model developed by [21]. The predictive capability of the block B4 and neural network approach is compared using statistics, which showed that neural network predictions for surface roughness were only 2 % closer to the experimental measurements, in comparison to block B4. A much simpler and faster B4 is used in the simulation model, despite the lower accuracy of predictions. When the simulation model of the feed drive system (shown in Fig. 3) and the simulator of cutting process are replaced by the real machine tool, a harmonized system for dynamic adjustment of feed rate is obtained. The system with an adjustment of feed rate override assures the required surface roughness. The following algorithm is used to on-line adjust the feed rate. fc 1.1 - f(t -1) for^f > 1 F 0.9 - f(t -1) for Ff < 1 (6) where, fc is a feed rate command [mm/rev], and f(t-1) previous actual feed rate command. This algorithm determines the optimal feed rate required to obtain the controllable surface roughness. 3 SIMULATION RESULTS AND DISCUSSION Via simulations, the efficiency and stability of the control model with different requirements of the surface quality are tested. The criterion for the efficiency of the control model is the difference between the desired and simulated Ra. The starting cutting conditions and the desired Raref are the input data. The following paragraphs present the simulation example (Sim. No. 1 and Sim. No. 3; Table 1) in which the model-based control was applied to demonstrate its performance. The simulation outputs are shown in Fig. 5. Recommended cutting conditions for each cut are first determined off-line, and then simulation tests are made without and with control model action. The recommended cutting parameters are determined with the prediction model [19]. The spindle speed is set at 610 rpm, and the resultant cutting force undergoes random changes up to 10 %. The simulation time is 9) 5 h) 10 ■) Fig. 5. Simulated control model response with in without feed rate adjustment Table 1. The simulation requirements and result Desired surface Results without control model action Results with control model action Sim. No. roughness f [mm/rev] f [mm/rev] Actual surface Raref [mm] F [N] F[N] roughness Ra [mm] 1 1.0 460.63 0.17 448.06 0.15 0.95 2 1.3 465.54 0.18 522.86 0.17 1.28 3 1.6 556.23 0.2 582.12 0.2 1.59 4 1.9 628.78 0.23 615.5 0.22 1.87 5 2.2 641.5 0.24 638.28 0.23 2.1 6 2.5 661.83 0.25 683.73 0.25 2.45 7 2.8 735.77 0.28 727.0 0.27 2.77 8 3.1 762.35 0.29 750.0 0.28 2.98 15 s. In Table 1, the requirements and the results of simulations are indicated. The undesirable aspects of the experimental results of the measured cutting forces were represented by the model of machine noise, which is shown on Fig. 3. This model has been generated to simulate the uncontrollable measurement noise resulting especially from the non-homogeneity of raw material, tool wear and tool breakage. The machine noise model is generated by using a random of variation, which relates to the settings of cutting force data acquisition software. It adds noise in the simulation model in order to test the reference of the output with the effect of the noise. When the noise is eliminated from the model, an ideal simulation result was obtained. The signals of resultant cutting force and the quality of surface are mutually related and have identical trends. Fig. 5d and e show the surface roughness of the workpiece without and with control model action. The Ra in Fig. 5d is not constant when machining without control action, because the feed rate is not adjusted when the output exceeds the reference value. Fig. 5d is the response of the surface roughness where there is no feed rate adjustment. The Ra is constant when using dynamic adjusting of feed rates (Fig. 5d; Section B to D). The control model is initiated by the selection of the reference value Raf afterwards, reference cutting force Fref is predicted according to reference block B1. To use the proposed control system and to adjust the feed rate, the desired surface roughness is 1.6 mm, the recommended feed is 0.2 mm/rev and its allowable adjusting rate is [0 % to 150 %]. The data of simulation example is marked as No. 3 and is presented in Table 1. The depth of cutting was set to 1 mm. When the reference force is known, the command value fc is determined according to the B3 block. The initial value fc is 0.2 mm/rev, thus allowing the feed to be overridden from 0 mm/rev to 0.3 mm/rev. The control bock is changing the initial value fc, until the optimum controlled feed rate of 0.2 mm/rev has been reached. The CNC servo drive model generates the actual feed rate f from the command signal fc. Dynamic adjustment of feeding is a prerequisite for maintaining constant surface roughness of 1.6 mm. The actual Ra is determined by the output variables block B4. The result is a roughness of 1.59 mm, which is acceptable if compared with the desired value of 1.6 mm (Fig. 5e). In the first simulation using constant feed rates the Ra reaches its proper value only in the last second of Sections A and B (Fig. 5d). However, by applying control action, the average Ra achieved is much more close to the reference Ra. It reaches and maintains the value of 1.59 mm (Fig. 5e). The resultant cutting force for the control model is maintained at about 570 N (Fig. 5b). The control system responds to the fall of the desired Ra by immediate reduction of feeding (Fig. 5f; Section C to D). As a result, the cutting force drops to value of 448 N (Fig. 5c). It can be seen that by continuously adjusting the feed rate (section B to D), the cutting force is maintained at a constant so that the desired surface roughness can be attained. The results confirm that the proposed control model is efficient in obtaining the desired surface roughness and that it is efficient in fine machining. Constant cutting forces lead to better surface quality. Although the system remains stable in all performed simulation tests, one obvious problem appears if the measured resultant cutting force approaches zero. In this situation, the feed rate command values would increase to infinity. A very small value of cutting force would lead to incorrect increase in the feed rate, causing excessive forces when the tool enters the workpiece. This problem has been efficiently eliminated by restricting the controller action to sample periods when the measured forces are above a certain user-specified threshold for a specified time. In our study, the threshold to start the control action is set to be 100 This solution ensures that the control algorithm is activated only when the tool is engaged in the workpiece. When the forces are bellow a user-defined threshold, the tool is assumed to be out of the work piece, and the control actions are paused. The developed solution also includes the restrictions of the control actions for feedrates limited by chip thickness. After simulations, the simulation model of CNC servo system (shown in Fig. 5) and the simulator of cutting process were replaced by the real CNC machine tool GF NDM 16. The performance of the suggested system was then verified by a serial of experimental works on the lathe for Ck 45 steel workpieces with variation of cutting conditions. The machine tool was connected with a personal computer. The connection was designed to allow GF NDM16 to be under digital control implemented in the PC. The model control was developed in Matlab. The experimental testing results are not described in this paper. 4 CONCLUSION In this research, the model-based controller is proposed to control the turning process with constant surface roughness. The control model modifies the feed rate in real time in order to keep the surface roughness constant and to make machining more efficient. The surface roughness control model was developed based on the simplified models of the feed drive servo system and the turning process. The correlation between the process output variables (cutting forces, surface roughness) and cutting parameters is determined by using statistics. A high degree of correlation is found between surface roughness and the resultant cutting forces for various machining conditions. The effects of the developed control model on Ra control have been studied through simulations. They have shown that the proposed control model makes the system insensitive to disturbances due to tool wear, non-homogeneity of work piece material and variations of depth of cutting. The performance of the developed control model is much better in comparison to classic CNC turning in the aspects of surface roughness fluctuation. The control system was also realized experimentally due to the excellent results of the simulations. 5 REFERENCES [1] Qolak, O. (2014). Optimization of machining performance in high-pressure assisted turning of Ti6Al4V alloy. Strojniški vestnik - Journal of Mechanical Engineering, vol. 60, no. 10, p. 675-681, D0l:10.5545/sv-jme.2013.1079. [2] Grigoriev, S.N., Starkov, V.K., Gorin, N.A., Krajnik, P., Kopač, J. (2014). Creep-feed grinding: an overview of kinematics, parameters and effects on process efficiency. Strojniški vestnik - Journal of Mechanical Engineering, vol. 60, no. 4, p. 213-220, D0l:10.5545/sv-jme.2013.1547. [3] Koren, Y. (1983). Computer Control of Manufacturing Systems. McGraw-Hill, New York. [4] Stute, G., Goetz, F.R. (1995). Adaptive control system for variable gain in ACC systems. Proceedings of the 16th International Machine Tool Design and Research Conference, Manchester, p. 117-121. [5] Pitstra, W.C., Pieper, J.K. (2000). Controller designs for constant cutting force turning machine control. ISA Transactions, vol. 39, no. 2, p. 191-203, D0I:10.1016/S0019-0578(00)00008-2. [6] Hwang, C.L. (1993). Adaptive turning force control with optimal robustness and constrained feed rate. International Journal of Machine Tools and Manufacture, vol. 33, no. 3, p. 483-493, D0l:10.1016/0890-6955(93)90054-X. [7] Ruey, J.L. (2012). Intelligent control of a constant turning force system with fixed metal removal rate. Applied Soft Computing, vol. 12, no. 9, p. 3099-3111, D0l:10.1016/j.asoc.2012.05.012. [8] Daeshmend, L.K., Pak, H.A. (1986). Model reference adaptive control of feed force in turning. ASME Journal of Dynamic Systems, Measurement and Control, vol. 108, no. 3, p. 215222, D0l:10.1115/1.3143770. [9] Landers, R.G., Ulsoy, G.A. (1998). Supervisory machining control: design approach and experiments. CIRP Annals - Manufacturing Technology, vol. 47, no. 1, p. 301-306, D0l:10.1016/S0007-8506(07)62838-8. [10] Landers, R.G., Ulsoy, A.G., Ma, Y.-H. (2004). A comparison of model-based machining force control approaches. International Journal of Machine Tools and Manufacture, vol. 44, no.7-8, p. 733-748, D0l:10.1016/j. ijmachtools.2004.02.005. [11] Tomizuka, M., Oh, J.H., Dornfeld, D.A. (1983). Model reference adaptive control of the milling process. ASME Control of Manufacturing Processors and Robotic Systems, vol. 11, p. 55-64. [12] Lauderbaugh, L.K., Ulsoy A.G. (1989) Model reference adaptive force control in milling. ASME Journal of Engineering for Industry, vol. 111, no. 1, p. 13-21, D0l:10.1115/1.3188726. [13] Yamazaki, K., Kojima, N., Sakamoto, C., Saito, T. (1991). Real-time model reference adaptive control of 3-D sculptured surface machining. CIRP Annals - Manufacturing Technology, vol. 40, no. 1, p. 479-482, D0I:10.1016/S0007-8506(07)62034-4. [14] Psaltis, D, Sideris, A, Yamamura, A.A. (1988). A multilayered neural network control. IEEE Control Systems Magazine, vol. 8, no.2, p.17-21, DOI:10.1109/37.1868. [15] Tarng, Y.S., Hwang, S.T., Wang, Y.S., (1994). A neural network controller for constant turning force. International Journal of Machine Tools and Manufacture, vol. 34, no. 4, p. 453-460, D0l:10.1016/0890-6955(94)90077-9. [16] Zuperl, U., Cus, F. (2004). Tool cutting force modeling in ballend milling using multilevel perceptron. International Journal of Materials Processing Technology, vol. 153-154, p. 268-275, D0l:10.1016/j.jmatprotec.2004.04.309. [17] Korkut, I., Donertas, M.A. (2007). The influence of feed rate and cutting speed on the cutting forces, surface roughness and tool-chip contact length during face milling. Materials & Design, vol. 28, no. 1, p. 308-312, D0l:10.1016/j. matdes.2005.06.002. [18] Brezak, D., Majetić, D., Udiljak, T., Kasać, J. (2014). Tool wear estimation using an analytic fuzzy classifier and support vector machines. Journal of Intelligent Manufacturing, vol. 23, no. 3, p. 797-809, D0l:10.1007/s10845-010-0436-x. [19] Cus, F., Zuperl, U. (2006). Approach to optimization of cutting conditions by using artificial neural networks. Journal of Materials Processing Technology, vol. 173, no. 3, p. 281-290, D0l:10.1016/j.jmatprotec.2005.04.123. [20] Zuperl, U., Cus, F. (2003). Optimization of cutting conditions during cutting by using neural networks. Robotics and Computer-Integrated Manufacturing, vol. 19, no. 1-2, p. 189199, D0l:10.1016/S0736-5845(02)00079-0. [21] Asiltürk, I., Qunka§, M. (2011). Modeling and prediction of surface roughness in turning operations using artificial neural network and multiple regression method. Expert Systems with Applications, vol. 38, no. 5, p. 5826-5832, D0l:10.1016/j. eswa.2010.11.041. Strojniški vestnik - Journal of Mechanical Engineering 61(2015)4, 254-264 © 2015 Journal of Mechanical Engineering. All rights reserved. D0l:10.5545/sv-jme.2014.2185 Original Scientific Paper BEM-Based Algorithm for URANS Simulations of Flow over a Square Cylinder Janez Lupše* - Leopold Škerget - Jure Ravnik University of Maribor, Faculty of Mechanical Engineering, Slovenia The boundary element method (BEM) is extended for the discretisation of the incompressible Navier-Stokes equations in the velocity-vorticity form for the URANS formulation. BEM uses fundamental solutions as weighting functions, thus the approach enables inclusions of some of the governing physics in the method itself. For testing purposes the flow around the square cylinder in the two-dimensional channel with a blockage ratio of 1/8 was chosen as the appropriate test case. The flow was calculated for the Reynolds numbers ranging from 100 to 1000. At Reynolds number 1000, the governing equations were solved directly in a quasi DNS manner and in the RANS form using the Spalart-Allmaras turbulence model within an URANS simulation. Keywords: boundary element method, velocity-vorticity form, incompressible flow, URANS, square cylinder Highlights • Extension of the boundary element based algorithm is proposed for simulation of turbulent flows. • Velocity vorticty formulation of Navier - Stokes equations is employed. • Simulation of flow over a square cylinder inside a channel at blockage of 1/8. • Simulations were performed within the unsteady laminar regime of the flow and within the low turbulent regime of the flow. • Comparison between the unsteady RANS model and the quasi-DNS simulation. 0 INTRODUCTION The boundary element method has been seen as a very successful method for solving various problems. In general, this method has an advantage over other numerical methods when the fundamental solution to the governing equations exists, thus effectively reducing the dimension of the problem [1]. Using the description of the governing physics regarding the fundamental solution, a stable and highly accurate numerical scheme is obtained, even on coarse numerical meshes. Unlike most of the other discretisation methods, using fundamental solutions this method relies on capturing the governing physics instead of the brute force. Furthermore, the derivatives of the transported variables, normal to the boundary of the domain or sub-domain, are a part of the solution of the system. This enables accurate representation of these values and the corresponding integral parameters such as the Nusselt number for example. Unfortunately, the governing equations for the fluid flow in general are of convective-diffusive type with important and non-trivial source terms. Even when the convective-diffusive fundamental solution is used, discretisation of the whole solution domain is needed. In addition, when a fundamental solution that includes more of the governing physics of the flow is used (with the exception of the Laplace fundamental solution), the integrals that arise from the discretisation of the governing equations need to be re-evaluated every time the corresponding field of variables (e.g., velocity) is changed. However, the inclusion of the flow physics into the discretisation usually leads to coarser numerical meshes. Whilst there are various approaches for the discretisation of domain terms (dual reciprocity method [2], the green element method [3], etc.) arising regarding the discretisations of the Navier-Stokes equations the so-called boundary-domain integral method (BDIM) [4] is used in our case. This approach gives discretisation of the solution domain similar to the finite element method, with each cell representing a separate sub-domain. The appropriate compatibility and equilibrium conditions have to be applied in order to interconnect the sub-domains. This results in an over-determined linear system of equations with a non-symmetric sparse system matrix that is solved by the LSQR solver [5] using diagonal-preconditioning. In Lupše et al. [6] the method was extended for solving the Reynolds averaged Navier-Stokes equations (RANS). This article is aimed at testing the developed algorithm for the intrinsically unsteady flows that can also be described using the unsteady RANS (URANS) approach. In regard to this purpose the flow over a square cylinder inside a channel was chosen as a relevant test case. The confined flow around the square cylinder is a challenging test for numerical methods as a multitude of interesting flow phenomena are present within relatively simple geometry. This type of flow 254 *Corr. Author's Address: University of Maribor, Faculty of Mechanical Engineering, Smetanova 17, 2000 Maribor, Slovenia, janez.lupse@yahoo.ca is important during many engineering applications ranging from the cooling of electronics to building aerodynamics. Although free-flow around the circular cylinders has been the focus of extensive research [7], less effort has been put into researching the confined flow around bluff bodies, especially around the square cylinders. Additional parameters need to be defined for the confined flow, namely those inflow velocity profiles and blockage ratios of the channel that do not appear in free-flow cases. Davis et al. [8] performed experiments and numerical simulations of the flows for two different blockage ratios (B = 1/4 and B = 1/6) and a wide variety of Reynolds numbers ranging from 100 and up to 1800. Later, [9] performed numerical simulations for blockage ratio of 1/8 and Reynolds numbers of the flow up to 900. A common argument about simulations of [9] is that the meshes used were too coarse around the vicinity of the cylindrical obstacle, hence the discrepancies of the obtained results with other authors. Breuer et al. [10] tested two different numerical methods on this type of flow at the blockage ratio of 1/8 for Reynolds numbers up to 300. Their results are widely considered as the benchmark for the laminar part of the flow. The results of [10] were later confirmed by the so-called multi-particle collision dynamics method tested at Reynolds numbers lower than 130 [11]. Galletti et al. [12] performed a test of a proper orthogonal decomposition-based model on the confined flow over a square cylinder and obtained comparable results to [10] having a blockage ratio of 1/8. Later the flow was also simulated for various laminar Reynolds numbers and extended those simulations to different blockage ratios and non-isothermal flows [13] and [14]. The flow chosen for the test case featured many interesting phenomena at different flow regimes. At the Reynolds number equal to 1 creeping flow is observed. The flow separates from the trailing edge of the cylinder with any increase in the Reynolds number. Although two recirculation bubbles are present and increase in length as the Reynolds number of the flow increases, the flow is steady up to a critical point. Authors in literature have reported different values for this critical number but agree that it lies somewhere between Re = 54 [15] and Re = 70 [16]. After this point the flow becomes unsteady but periodic. Separation occurs at the trailing edge at first and with increasing Reynolds number of the flow moves to the leading edge of the cylinder. Flow separation causes a distinct pattern of flow, known as the Von Karman vortex street. Spatial effects occur at about the Reynolds number 300 [10]. As the Reynolds number of the flow is defined by the cylinder diameter and not the channel width it is necessary to have in mind that with a used blockage ratio of 1/8 at about this Reynolds number, transition to turbulence in the channel flow is slowly starting to occur. Thus in reality, the flow beyond the Reynolds number 300 is be-coming increasingly turbulent. Consequentially, for the inlet profile the turbulent velocity profile would be more appropriate passed these Reynolds numbers of the flow. The data available in literature for the Reynolds numbers up to 1000 [8] is obtained by imposing the standard parabolic velocity profile. As such, the parabolic velocity profile in this work, used at the inlet boundary condition, was used purely for the testing purposes of the developed numerical algorithm and code in order to enable easier and more accurate comparison with the available data of [9]. 1 GOVERNING EQUATIONS Whilst the more common approach to describing the fluid flow is the use of the Navier-Stokes equations in their primitive forms, the velocity-vorticity formulation is used in our work. In this formulation the flow equations are split into kinematic and kinetic parts. In general, the approach gives us the advantage of eliminating the pressure gradient from the governing equations. However, the boundary conditions for the vorticity field in Eqs. (1) and (2) are generally unknown intuitively, unlike the velocity boundary condition in the momentum transport equation. The flow kinetics is described by the vorticity transport equation obtained from the momentum transport equation by applying the curl operator [17]. For the incompressible planar flow the following equation is obtained: da TT da — + U-= v0- dt dxj dxj dxj da 1 dpg, --e„ Po " dxj (1) where ю is used for the vorticity, Uj for the components of the velocity vector and Xj the spatial coordinates; v0 the kinematic viscosity, p0 the fluid density and g, the components of the gravity acceleration vector. The vorticity vector is treated as a scalar within the planar flow since only the component of the vorticity normal to the plane of the flow has a non-zero value. In the RANS form of Eq. (1), an additional source term is obtained on the right hand side and the RANS form of the vorticity transport equation is then [6]: да TT da — + U — = v0 d2 a 1 dpgi 1 d/r dt dx; dxj dx; p0 dx; p0 dx; (2) where the additional source term fm is defined by: e dvtw dvt 2 dvt dUi Ji = -ßij-+ -® + 2--. (3 ) dxj dxj dxj dxj and constants of the used model see the original reference [20]. 2 BOUNDARY-DOMAIN INTEGRAL EQUATIONS The effective viscosity is now comprised of the molecular part v0 and the modelled part vt . The overbar indicating the averaging process is omitted for the sake of brevity. The kinematics part is formed by the Poisson equation that links the velocity and vorticity fields. It is derived at from the continuity equation and the vorticity definition [1]. For the planar flow the following form of the kinematics equation is obtained: d2Ui dm n -—+eij — = 0. dx ; dx ; dxk (4) j j 4 Prior to discretisation using the boundary element method, the governing equations need to be rewritten in integral form. If we use Q for the solution domain and Г for its boundary, then the integral form of the kinematics Eq. (4) can be written as c № (O+Ju du* d r-ßu*d г ey fan,u d Г- eУ f ю^— d Q 1J 1 У J Ят Q dxj (8) There exist very efficient algorithms for solving the velocity-vorticity forms of the equations [18] but they are only applicable for special cases using simplified geometry. Care has to be taken however since Eq. (4) admits solutions in which neither solenoidality of the velocity field nor vorticity definition hold true [19]. In order to obtain a general solution, an additional equation has to be included [1]: dU: dUj -+ e:jn : ю = -e:j- dn 1 1 j dt (5) where the normal and tangential derivatives in respect to domain boundary are defined as: щ dn DX : (6) where we use % for the location of the source point, c for the geometric coefficient and u* to denote the Laplace fundamental solution. * 1 1 u = — ln —. 2n r (9) In Eq. (9) r is used for the distance between the integration and source points. The detailed derivation of the integral form of Eq. (8) can be found in [21]. Eq. (8) is then used for calculation of the velocity field inside the solution domain. In order to calculate the boundary conditions for the vorticity transport equation in a general case Eq. (4) has to be combined with Eq. (5). A new integral form of the kinematics Eq. (4) is obtained; c m («и и dn d г and n and t are the unit normal and unit tangent vectors. Finally, in order to close the system of equations, the Spalart-Allmaras model [20] is used as a closure model. The main equation of this closure model without the trip term is: D = Си [1 - Ft 2 ] Д-f- Dt a dx : t ~\ dv ( +V t . j / +-C a b2 dX: Cw1Fw ' C b1 21 Ft2 V d (7) Eq. (7) is the transport equation for the so-called modified turbulent viscosity from which modelled viscosity vt, as used in Eq. (2) is obtained by simple algebraic expression. For the details on the functions dU-dr-eJ«^ (10) a dxj with the tangential derivative defined as in Eq. (6). The vorticity transport equation is a non-stationary convective-diffusive equation with nonlinear source terms. As such, several fundamental solutions can be used as weighting functions. If the parabolic diffusion fundamental solution is used the following integral form of the vorticity transport equation is obtained [6]; du dea c tk )+ la-dT - |— u dT = V W V k! T dn T d" -\ak-iul-id n '0 n = n 2 J wUjnjWdГ - J mU du d Q vr Q dxj У voPo J PginjUd r-Jpg, du \ vr ( dX; J У J frnjud r-J fi du d Q * \ d Q vr Q dXj У (11) The index k is used for denoting the current time step, thus tk denotes the current time. Instead of the parabolic diffusion fundamental solution the convection-diffusion fundamental solution can also be used. The following form of the equation is obtained [6]; с (š)ffl(š, tk ) + \wQfd Г-jfVd Г = Г Г = - — \(®k-1 + ß2®k-2 )u"dfi-V0 fi -V\ю( -UJ°hUdГ- Vo Г 1 М( - и o )) d о- ofi VoPo Jpginju*d Г-jpgl du \ Чг дХ; J J frnju*d Г-J fi du d fi * Л dfi Чг _ dx; fi j у (12) where we use ßj = -2/At, ß2 = -1/(2At) and Q* = du*/dn + (Uj0nj0/v0)u*, for the time discretisation. As the transport equation of the turbulence model used is of the same type as the vorticity transport equation, a similar integral equation is derived at: [c фу 3.0.C0;2-6. [3] Taigbenu, A.E. (1995). The green element method. International Journal for Numerical Methods in Engineering, vol. 38, no. 13, p. 2241-2263, DOI:10.1002/nme.1620381307. [4] Škerget, L., Rek, Z. (1995). Boundary-domain integral method using a velocity-vorticity formulation. Engineering Analysis with Boundary Elements, vol. 15, no. 4, p. 359-370, D0I:10.1016/0955-7997(95)00036-N. [5] Paige, C.C., Saunders, M.A. (1982). LSQR: An algorithm for sparse linear equations and sparse least squares. ACM Transactions on Mathematical Software (TOMS), vol. 8, no. 1, p. 43-71, DOI:10.1145/355984.355989. [6] Lupše, J., Škerget, L., Ravnik, J. (2014). Velocity-vorticity RANS turbulence modeling by boundary element method. Engineering Analysis with Boundary Elements, vol. 39, p. 4452, DOI:10.1016/j.enganabound.2013.11.004. [7] Ternik, P., Rudolf, R. (2014). Laminar forced convection heat transfer characteristics from a heated cylinder in water based nanofluids. International Journal of Simulation Modelling, vol. 13, no. 3, p. 312-322, DOI:10.2507/IJSIMM13(3)5.271. [8] Davis, F.W., Moore, F., Purtell, P. (1984). A numerical-experimental study of confined flow around rectangular cylinders. Physics of Fluids, vol. 27, no. 1, DOI:10.1063/1.864486. [9] Mukhopadhyay, A., Biswas, G., Sundararajan, T. (1992). Numerical investigation of confined wakes behind a square cylinder in a channel. International Journal for Numerical Methods in Fluids, vol. 14, no. 12, p. 1473-1484, DOI:10.1002/fld.1650141208. [10] Breuer, M., Bernsdorf, J., Zeiser, T., Durst, F. (2000). Accurate computations of the laminar flow past a square cylinder based on two different methods: lattice-boltzmann and finite volume. International Journal of Heat and Fluid Flow, vol. 21, no. 2, p. 186-196, DOI:10.1016/S0142-727X(99)00081-8. [11] Lamura, A., Gompper, G. (2002). Numerical study of the flow around a cylinder using multiparticle collision dynamics. The European Physical Journal E: SoftMatter and Biological Physics, vol. 9, no. 5, p. 477-485, DOI:10.1140/epje/i2002-10107-0. [12] Galletti, B., Bruneau, C.H., Zannetti, L., Iollo, A. (2004). Low-order modelling of laminar flow regimes past a confined square cylinder. Journal of Fluid Mechanics, vol. 503, p. 161170, DOI:10.1017/S0022112004007906. [13] Turki, S., Abbassi, H., Ben Nasrallah, S. (2003). Two-dimensional laminar fluid flow and heat transfer in a channel with a built-in heated square cylinder. International Journal of Thermal Sciences, vol. 42, no. 12, p. 1105-1113, DOI:10.1016/S1290-0729(03)00091-7. [14] Turki, S., Abbassi, H., Nasrallah, S.B. (2003). Effect of the blockage ratio on the flow in a channel with a built-in square cylinder. Computational Mechanics, vol. 33, no. 1, p. 22-29, DOI:10.1007/s00466-003-0496-2. [15] Kelkar, K.M., Patankar, S.V. (1992). Numerical prediction of vortex shedding behind square cylinders. International Journal of Numerical Methods for Fluid, vol. 14, no. 3, p. 327-341, DOI:10.1002/fld.1650140306. [16] Okajima, A. (1982). Strouhal numbers of rectangular cylinders. Journal of Fluid Mechanics, vol. 123, p. 379-398, DOI:10.1017/S0022112082003115. [17] Škerget, L., Samec, N. (2005). BEM for the two-dimensional plane compressible fluid dynamics. Engineering Analysis with Boundary Elements, vol. 29, no. 1, p. 41-57, DOI:10.1016/j. enganabound.2004.08.007. [18] Kim, J., Moin, P., Moser, R. (1987). Turbulence statistics in fully developed channel flow at low Reynolds number. Journal of Fluid Mechanics, vol. 177, no. 1, p. 133-166, DOI:10.1017/ S0022112087000892. [19] Wu, J.C., Thompson, J.F. (1973). Numerical solutions of time-dependent incompressible Navier-Stokes equations using an integro-differential formulation. Computers & Fluids, vol. 1, no. 2, p. 197-215, DOI:10.1016/0045-7930(73)90018-2. [20] Spalart, P.R., Allmaras, S.R. (1994). A one-equation turbulence model for aerodynamic flows. La Recherche Aérospatiale, vol. 1, no. 1, p. 5-21, DOI:10.2514/6.1992-439. [21] Škerget, L., Ravnik, J. (2009). BEM simulation of compressible fluid flow in an enclosure induced by thermoacoustic waves. Engineering Analysis with Boundary Elements, vol. 33, no. 4, p. 561-571, DOI:10.1016/j.enganabound.2008.08.003. [22] Ravnik, J., Škerget, L., Hriberšek, M. (2004). The wavelet transform for BEM computational fluid dynamics. Engineering Analysis with Boundary Elements, vol. 28, no. 11, p. 13031314, DOI:10.1016/j.enganabound.2004.05.002. [23] Orlanski, I. (1976). A simple boundary condition for unbounded hyperbolic flows. Journal of Computational Physics, vol. 21, no. 3, p. 251-269, DOI:10.1016/0021-9991(76)90023-1. Strojniški vestnik - Journal of Mechanical Engineering 61(2015)4, 265-272 © 2015 Journal of Mechanical Engineering. All rights reserved. D0l:10.5545/sv-jme.2014.2291 Original Scientific Paper Research on the Position-Pressure Master-Slave Control for a Rolling Shear Hydraulic Servo System Wang Jun1 - Sun Binyu1,* - Huang Qingxue1 - Li Hongzhou2 - Han Heyong1 1 Taiyuan University of Science and Technology, China 2 Hebei WenFeng Iron and Steel Co., China The position-pressure master-slave control of hydraulic servo system is a dual closed-loop system, which takes position closed loop as the main control system, pressure closed loop as the supplementary, bringing these two seemingly contradictory theories together while simultaneously control the system. When the position error value is bigger than the threshold value, the system will transfer the pressure signal to the signal which will be added to the position signal, and regulate the system. This paper first builds a mathematical model of the hydraulic servo system position-pressure master-slave control, then simulates the system with AMEsim and Matlab softwares, which verifies the feasibility of this function. Finally, experiments are conducted using the new function, and the results show that position-pressure master-slave control can improve the precision and stability of hydraulic servo system; this new function can also be applied to other high speed, high precision and heavy equipment. Keywords: position control, pressure control, master-slave control-hydraulic servo system, simulation, dual closed-loop system Highlights • Proposed a new position-pressure master-slave control function for double-sided rolling shear. • Calculating a position-pressure conversion model. • Building a mathematical model of the hydraulic system of rolling shear. • Improving the control precision and stability of hydraulic system. 0 INTRODUCTION Hydraulic double-sided rolling shear, the most advanced shearing machine, cuts off steel plates into a specified size [1], which is an important part of rolling production line and it will determine whether the quality of plates is good or not. Asymmetrical, valve controlling an asymmetrical cylinder (shown in Fig. 1) is the key to the hydraulic system of double-sided rolling shear, which uses a position-closed loop to make the cylinders move in order to ensure the rolling shear cutting off the steel plates according to the routine action. Hydraulic rolling shear cuts off the steel plates with two cylinders moving harmonically along the given curve; the process is shown in Fig. 2. First, the left cylinder moves (shown as Fig. 2a). Then, when it moves to the designated place, the right cylinder begins to move (shown in Fig. 2b). Next, when the blade is moving to the place as shown in Fig. 2c, the rolling shear completes the action. Finally, the two cylinders retract to the initial position; a full action is completed, as shown in Fig. 2d. This equipment is simple, powerful, and highly reliable. Usually, hydraulic rolling shear adopts position control to make the cylinder move along the specified curve, while (when shearing thick steel plates) the pressure of the two cavities of the cylinder is always P1 il---1 T II-- Fig. 1. The hydraulic scheme of hydraulic rolling shear too low to push the upper blade moving along the desired curve to cut off the steel plates. As a result, the quality of the sheared plates will be poor (shown as Fig. 3). Therefore, how to improve the effect of controlling of rolling shear is a widely studied issue *Corr. Author's Address: Taiyuan University of Science and Technology, No. 66 Waliu Street, Taiyuan, China, sunby_sby@126.com 265 C D Fig. 2. The shearing process of hydraulic rolling shear [2]. Some researchers want to improve the structure of the rolling shear. Huang et al. design a new mechanism [3], which is powerful but heavy. Optimizing the mechanism of rolling shear cannot solve the problem of powerless; some researchers propose new control algorithms to improve the performance of the position control: Quan and Xu [4] put forward a new control method combining pressure and position in series and parallel, while Madani and Moallem [5] propose position/force control, but both of them are unstable when switching the pressure either to position control or position to pressure control. Bessa et al. [6] put forward a new type of control using sliding mode control, while Kim et al. [7] use QFT to perfect the position control. These are too complex to apply to the production in factories. This paper puts forward a new control method: pressure-position master-slave control for double-sided rolling shear, which takes the position closed loop as the main control system, and the pressure closed loop as supplementary. This new function can make the cylinder move along the desired curves and will not be affected by the external loads, so the shape of the sheared plates will be improved. It has both theoretical significance and practical significance for improving the control precision and performance of heavy and high-speed equipment. Fig. 3. Defect of steel plate 1 BUILD THE MATH MODEL The structure diagram of position-pressure master-slave control is shown in Fig. 4. The pressure signal from cylinder is detected by the pressure sensors and sent to the pressure-position conversion model; the master-slave controller will Fig. 4. Structure diagram of position-pressure master-slave control then judge when the pressure signal will be added to the position signal, if the position signal error is bigger than the threshold, the pressure closed loop is connected, the converted pressure signal is a compensation for the position signal to adjust the input signal, rapidly enlarging the flow of servo valve , and accelerating the speed of cylinder, so the position error of cylinder will be small, the pressure and position of cylinder will be in control. Meanwhile, the position control will not be affected by the pressure control. Thus, this new function can simultaneously control the position and the pressure of the cylinder. The schematic of the asymmetrical servo valve controlling asymmetrical cylinder is shown in Fig. 4. This paper researches the mathematical model when the blade begins to contact the steel plates, FL is a constant, the equilibrium equation for the servo cylinder and the loads at two adjacent moments, which are the beginning and ending points of the shortest time that the pressure sensor can obtain: AiPt-i = mtat-i + Bpxt-i + Kxt-i + fl , (!) AlPt = mtat + Bpxt + Kxt + fl . (2) Eq. (3) is subtracted from Eq. (2) to Eq. (1): - - Ai M t ( At - At-0 Bp (Vt - Vt-1) Pt - P,-i K Pt - Pt-i Pt - P-i (3) When shearing the plates, Mt (At _ At _ ) and Bp (Vt _ Vt_j ) will be ignored because the accelerated speed and speed are changing so small. The simplified Eq. (3) is shown as Eq. (4). k = Xt - X<- Pt - Pt- A к. (4) Thus, the computation formula of position-pressure conversion coefficient is shown as Eq. (5): Kt = AX Xt - Xt-i AP Pt - Pt-1 (5) where Kt is the position-pressure conversion gain, AP is the difference in pressure between two moments, AX is the difference in position between two moments, Xt is the position signal of cylinder at NO.t moment, Xt-1 is the position signal of cylinder at NO.t-1 moment, Pt is the pressure signal of cylinder at NO.t moment, and Pt-1 is the pressure signal of cylinder at NO.t-1 moment. When the upper blade begins to contact the steel plate, the equilibrium equation is shown as Eq. (6): where p is the pressure of the blind end of the cylinder, p2 is the pressure of the rod end of the cylinder, A1 is the area of the piston, A2 is the area of the piston rod, and FL is the load of servo cylinder. The positive direction is where the rod stretches out, the pressure of system load is PL, flow is qL: Pl = Fl , L 4 and F = p Ax - p A2. So, the Eq. (7) can be transferred to be: (7) P = Fl = P A - P2 A A where X = — A A = P - P — = P-IP, A 2 (8) Fig. 5. The schematic of servo asymmetrical controlling asymmetrical cylinder The flow equation of the valve is: Il = Kqxv - KcPL . (9) Because this system uses asymmetrical valve to control asymmetrical cylinder, the flow gain Kq and flow-pressure coefficient Kc are shown as following: dxv K =dq± =_ C dpL л/2р(1 + Я)(Л - pi )' (10) (11) p A - Pi A = fl , (6) Assuming that the connect pipelines between servo valve and cylinder are short, symmetric and thick; the pressure loss due to frictional resistance and dynamic effect of pipe are ignored; the oil bulk modulus and the oil temperature are constants; the internal and external leakage are laminar flow. According to Fig. 5, the flow into the cylinder is q1 out of cylinder is q2, so, the flow equations are: The cylinder flow continuity equation can be obtained by linearizing Eq. (19) : , dxp ^ V dpL qL = A~t + CpPL +—<--Г\--M dt 2(1 -Я2dt (20) * = f + ^ ( -A) + ^ +ld~ì> (12) dV2 t W n ^ V2 dP2 q2 ~ + CP ( -p2 ) + C^ — , (13) where V1 is the effective volume of the blind end of the cylinder, V2 is the effective volume of the rod end of the cylinder, Cip is the internal leakage coefficient, Cep is the external leakage coefficient, and ße is the oil bulk modulus of elasticity. The volumes of the cavities of cylinder in the initial state are V10 and V20. When the piston of the cylinder is moving xp, the equations are as follows: V = V10 + Axp (14) where: Ctp is the leakage coefficient of cylinder and Cf„ = -1+4 C -C„, * 1 -42 p 1 -42 Equilibrium equation between cylinder and load [8] to [10] is Eq. (21): d2xp dxp APx - A P2 = + + Kxp + Fl , (21) dt dt andA2 = AAb the left of the Eq. (21) can be simplified as: A1P1 -A2P2 = Ax(px - Xf2) = A1pL. Then, we arrive at the Eq. (22) [11] to [13]: and V = V20 + A2 xp, dV1 = A dt dt dx„ dx— (15) (16) ^ = - A2= - A, ——. (17) dt dt dt The flow of cylinder is shown as following: qL = q1 +lq2. (T8) Put Eqs. (10) to (15) into (18), and simplify it: qL = qi +^2 = Cip (Pi- P2)(1+Я) + +C (Pi-Xp2) + (1+ Я2)ax ^ + YV0 dpL - ep\r1 ГЦ K J 1 dt ße dt -ьХж^г + (dpL + Я2 % (19) ße dt ße dt dt' d xp dxp p i j> p APl = + BP4t + Kxp + Fl , (22) dt dt where mt is the total mass of piston and loads refferred to piston, Bp is the viscosity damping coefficient, K is the spring stiffness of the load, and FL is the external load. Eqs. (9), (20) and (22) are Laplace transformed [14]: Ql = KqXv - KCPL , Ql = aisxp + C1pPl + Vt 2 (l-A2 ) (23) sPl , (24) APl = mts2 Xp + BpSXp + KXp + FL . (25) According to the equations above, the transfer function of the servo cylinder can be obtained, Eq. (25a): V. - K A A XP = 1 + Vt 2 (l-A2 )ßeKce Fr (25a) mtVt 2 (l-A2 ) A2 Kcemt aA VtBp 2 (l-A2 ) A,2 1 + A,2 KVt 2 (l-A2 ) A,2 5 +- KK A,2 3 2 + s The gain of position sensor is [15]: UA Kd = Y (26) where Ud is the output of the position closed loop. The transfer function of slave pressure control is: U p - PL ' KP ' Kt U e = U g - Up - Ud, (27) (28) where UP is the output of slave pressure control, KP is the gain of the pressure sensor, and Kt is the position-pressure conversion coefficient. 2 SIMULATIONS This paper uses AMEsim software to simulate the hydraulic system of rolling shear, and SIMULINK software to simulate its control system. The hydraulic system model of rolling shear in AMEsim is shown in Fig. 6, and the control system in SIMULINK is shown in Fig. 7. The system parameters in the simulation model are shown in Table 1. First, only the position closed loop is used in the simulation system. The position and pressure curves are shown in Fig. 8, which indicates that when the Fig. 6. The hydraulic system model in AMEsim Fig. 7. The control system model in SIMULINK external load is changing fast, the pressure will lag behind, and the cylinder cannot move along with the given curve, so that the rolling shear is not able to cut off the steel plates. Table. 1. Parameters of the system Working pressure 28 MPa Oil elastic modulus 14000 Parameters of 510/380 Cylinder stroke 477 mm cylinder mm Oil density 900 kg/m3 Dynamic mass 17000 kg Frequency respon- 35 Hz Flow P-A of servo 1000 L/min se of servo valve valve Pressure drop of 0.5 MPa Flow B-T of servo 500 L/min servo valve valve Damping ratio of 0.8 Flow A-T of servo 1000 L/min servo valve valve Rated current of 40 mA Flow P-B of servo 500 L/min servo valve valve Gain 1 10 Gain 2 0.5 Gain 3 3 Gain 4 0.1 Gain 5 0.01 Gain 4 Flow gain of the 1.7 m2/s Flow-pressure gain 2.1x10-12 servo valve of the servo valve m5/Nxs Fig. 9 shows the position and pressure curves when using position-pressure master-slave control, which indicates that the new function makes the system react fast and the position curve can follow the track of the given curve perfectly. The results indicate that the cylinder can move along with the given curve using position-pressure master-slave control, and the pressure of the cylinder is stable, which is not affected by the external load. Thus, the pressure-position master-slave control can improve the precision and the stability of the hydraulic servo system. 3 EXPERIMENTS The experiments are conducted on experimental double-sided rolling shear. The sheared steel plates are 30 mm in thick, 2100 mm in width. First, the experiments are conducted with the position closed loop only, and then using the position-pressure master-slave control instead. Finally, comparing the curves obtained by these two functions, the merits and demerits are analyzed. The servo valve used is 4WRTE32 produced by Rexroth; the pressure sensor used is HDA3840-A-400-000 produced by HYDAC. Table 2 is the experimental data using the closed loop position only, and the Fig. 10 is the position a) b) Fig. 8. a) The position curves and b) the pressure curves using position closed loop only a) Time [s] kj Fig. 9. a) The position curves and b) the pressure curves using position-pressure master-slave control curves and pressure curves fitted by software Origin according to the experimental data. Table 2. Experiment data Given displacement [mm] Actual displacement [mm] The rod The blind Time [s] end of the cylinder [mm] end of the cylinder [mm] 0 0 0 4.935 8.498 1 71.314 69.091 9.262 8.276 2 164.937 157.841 8.741 9.960 3 264.411 207.361 7.010 10.001 4 320 310.491 6.114 11.170 5 237.714 263.194 6.946 10.539 6 144.457 169.937 8.741 10.112 7 61.074 76.314 9.262 8.276 8 0 8.113 5.351 7.302 0 2 4 E 8 Time [s] Fig. 10. The position and pressure curves of rolling shear with position closed loop only Table 2 and Fig. 10 show that the actual position curve is defective; at about 2.3 s, the actual position curve is seriously deflected from the given position curve, the biggest error of position is 32 mm, and the pressure fluctuation is about 2.1 MPa. Table 3. Experiment data Time [s] Given displacement [mm] Actual displacement [mm] The rod end of the cylinder [mm] The blind end of the cylinder [mm] 0 0 0.264 4.9347 10.498 1 71.314 69.21 9.262 8.2764 2 164.937 163.18 9.7412 9.9609 3 264.411 262.832 6.2698 11.7554 4 320 320.424 5.9774 12.8276 5 237.714 241.722 6.9641 11.845 6 144.457 141.964 10.8895 6.982 7 61.074 59.064 10.4502 8.594 8 0 0.125 4.502 6.215 Then, the pressure closed loop is added to the system, using the position-pressure master-slave control. Table 3 shows data from laboratory, and the Fig. 11 shows position and pressure curves according to the data with the new method proposed in this paper. Time [s] Fig. 11. Position and pressure curves Table 3 shows the data of given displacement, actual displacement, pressure of the blind end of the cylinder and the pressure of the rod end of the cylinder in a whole process of rolling shear. Fig. 11 and Table 3 indicate that the biggest error between actual displacement and given displacement is 3 mm, the pressure fluctuation is less than what in which uses the closed loop position only. Compared to the results in Fig. 10, the position curve in Fig. 11 is perfect, and the following performance is better, so the rolling shear can cut off the steel plate to a level of quality that meets the standards of industrial production. Simply put, when double-sided rolling shear using position-pressure master-slave control to cut off the steel plates, the position curve of the cylinder can follow the calculated curves well, the pressure is stable, the quality of sheared plates and the production efficiency are both improved, which can meet the demands of clients. Position-pressure master-slave control can improve the precision of the position and pressure of cylinider, making the rolling shear cut thicker plates. 4 CONCLUSIONS Based on a practical issue, for solving the problems when the rolling shear cuts the steel plates, this paper puts forward a new position-pressure master-slave control function, which can improve the precision of position and pressure; putting these two seemingly contradictory theories together and improving the quality of steel plates. This paper builds a mathematical model of the hydraulic system of rolling shear, and then simulates the hydraulic system by using AMEsim and SIMULINK; the results indicate that the new theory can be applied to practice. Experiments are conducted; the pressure curves and position curves obtained via simulation and experiment show that the position-pressure master-slave control is feasible and verify the theory put forward by this paper has practical value. The position-pressure master-slave control proposed in this paper improves the anti-jamming capability, precision and stability of the system; thickness of sheared steel plates is increased by 35 % and the quality also be improved. It is easy to use in the plants and with other hydraulic servo systems with high speed, massive flows and heavy loads. 5 ACKNOWLEDGEMENT This paper was supported by Technology Innovation Project for Postgraduate of Taiyuan University of Science and Technology (20134023); Ph.D. Programs Foundation of Ministry of Education of China (20111103); The Initial Funding of Doctor Research of Taiyuan University of Science and Technology (20122047); Provincial Fund for Young Scholars (2013021019-2), Provincial Fund for Technological Innovation (2013101013). 6 REFERENCES [1] Han, H. Y., Huang, Q. X., Wang, J., Wang, J. (2013). The analysis of new hydraulic rolling shear servo system dynamic characteristics. Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture, vol. 227, no. 3, p. 453-459, D0l:10.1177/0954405412472667. [2] Yang, G. J., Yu, S. H., Zhang, H. C., Chen, G. D. (2009). Construction on the optimization designing system of shearing mechanism of rolling shear. Journal of Machine Design, vol. 29, no. 3, p. 63-66. (in Chinese) [3] Huang, Q. X., Ma, L. F., Li, J. B. (2008). Principle of asymmetric crank mechanism of new-type rolling shear. Chinese Journal of Mechanical Engineering, vol. 44, no. 5, p. 119-123, D0I:10.3901/JME.2008.05.119. (in Chinese) [4] Quan, L., Xu, X.Q. (2008). Simulation and test of electro-hydraulic servo position and pressure hybrid control principle. Chinese Journal of Mechanical Engineering, vol. 44, no. 9, p. 100-105 D0I:10.3901/JME.2008.09.100. (in Chinese) [5] Madani, M., Moallem, M. (2011). Hybrid position/force control of a flexible parallel manipulator. Journal of the Franklin Institute, vol. 348, no. 6, p. 999-1012, D0I:10.1016/j. jfranklin.2011.03.005. [6] Bessa, W.M., Dutra, M.S., Kreuzer, E. (2010). Sliding mode control with adaptive fuzzy dead-zone compensation of an electro-hydraulic servo-system. Journal of Intelligent and Robotic Systems, vol. 58, no. 1, p. 3-16, D0I:10.1007/s10846-009-9342-x. [7] Kim, J.W., Xuan, D.J., Kim, Y.-B. (2008). Design of a force control system for a dynamic road simulator using QFT. International Journal of Automotive Technology, vol. 9, no. 1, p. 37-43, D0I:10.1007/s12239-008-0005-x. [8] Tai, N.T., Kha, N.B., Ahn, K.K. (2010). Predictive position and force control for shape memory alloy cylinders. Journal of Mechanical Science and Technology, vol. 24, no. 8, p. 17171728, D0I:10.1007/s12206-010-0504-3. [9] Zhao, W. Z., Sun, P. K., Liu, S., Lin, Y. (2012). Multi-objective optimization of active steering system with force and displacement coupled control. Journal of Central South University, vol. 19, no. 4, p. 974-981, D0I:10.1007/s11771- 012-1099-x. [10] Murakawa, M., Lu, Y. (1997). Precision cutting of sheets by means of a new shear based on rolling motion. Journal of Materials Processing Technology, vol. 66, no. 4, p. 232-239, D0I:10.1016/S0924-0136(96)02529-0. [11] Lei, J. B., Wang, X. Y., Pi, Y. J (2013). Sliding mode control in position control for asymmetrical hydraulic cylinder with chambers connected. Journal of Shanghai Jiaotong University (Science), vol. 18, no. 4, p. 454-459, D0I:10.1007/s12204- 013-1419-9. [12] Huang, F. (2011). The Acceleration and Deceleration Algorithm and Simulation of Position and Pressure Control Conversion Process of Electro-Hydraulic Servo System. Wuhan University of science and technology, HuBei. [13] Han, H. Y., Wang, J., Wang, J., Huang, Q.X. (2012). Analysis of unsymmetrical valve controlling unsymmetrical cylinder stability in hydraulic leveler. Nonlinear Dynamics, vol. 70, no. 2, p. 1199-1203, D0I:10.1007/s11071-012-0523-y. [14] Yu, J., Li, Y. (2011). Simulation on electro-hydraulic proportional position control system under pressure boundary conditions. Advanced Materials Research, vol. 186, p. 21-25, D0I:10.4028/www.scientific.net/AMR.186.21. [15] Yi, J.G. (2015). Modelling and analysis of step response test for hydraulic automatic gauge control. Strojniški vestnik -Journal of Mechanical Engineering, vol. 61, no. 2, p. 115-122, D0I:10.5545/sv-jme.2014.2046. Strojniški vestnik - Journal of Mechanical Engineering 61(2015)4, 273-282 © 2015 Journal of Mechanical Engineering. All rights reserved. D0l:10.5545/sv-jme.2014.2123 Original Scientific Paper Study of an Annular Two-Phase Thermosyphon Used as an Isothermal Source in Thermometry Florencio Sanchez-Silva* - Ignacio Carvajal-Mariscal - Ariel E. Moreno-Cordobés -Pedro Quinto Diez - Miguel Toledo Velazquez National Polytechnic Institute, School of Mechanical and Electrical Engineering, Mexico The feasibility of using an annular two-phase thermosyphon (TPT) in thermometry was studied experimentally; the main results are presented in this paper The TPT was conceived as an isothermal source for calibration purposes, using the comparison method. For this purpose, an annular TPT was designed, built and characterized, using water as a working fluid. In order to determine experimentally how isothermal the TPT annular surface can be, eight thermocouples were installed axially on its surface. Tests were performed for two different cooling water flows 16.67 cm3/s and 4.17 cm3/s. For each cooling water flow, four heat power levels were supplied in the evaporation zone: (38, 154, 350, 1,385) W. The experimental results have shown that the thermal device works isothermal for power levels higher than 350 W, where the so-called bar effect is not present and the temperature measurement error is reduced. Key words: thermometry, two-phase thermosyphons, isothermal, annular, calibration Highlights • Using an annular two-phase thermosyphon as isothermal source in thermometry. • Immersion depth of the temperature sensor is 21 cm. • Higher values of heat supplied promoted isothermicity. • Device efficiency was higher than 94%. • The system accomplished all the requirements for calibration of temperature instruments using the comparative method. 0 INTRODUCTION The definition of a temperature scale, independent of the substance thermometric qualities used, was proposed by Kelvin in 1898 and was named the thermodynamic temperature scale, which coincides with a thermometer scale filled with a perfect gas. However, achieving this scale is practically impossible; therefore, another temperature scale was developed, the nearest possible to the thermodynamic one and fully characterized to remain convenient and maintain its precision. Taking into account all these specifications, the international temperature scale (ITS) was created; it is updated every 20 years. This scale is based on a primary set of fixed points of temperature, for instance: the triple point of water, the melting point of gallium, the freezing point of zinc among many others. One of the calibration methods is based on the use of cells containing the fixed points cited in the ITS. The other method of calibration is the comparison method, in which both the thermometer to be calibrated and the reference thermometer (calibrated one) are placed in an isothermal bath. The American Society for Testing and Materials (ASTM) established the calibration procedure to be followed when using this method [1]. For this purpose, a highly stable isothermal reservoir is required in order to avoid the so-called bar effect; this thermal phenomenon produces heat gains or losses due to the axial thermal gradient that appears when one end of the bar is in contact with the thermal source and transfers heat as a fin [2]. This situation leads to significant errors in temperature measurements, according to European Association of National Metrology Institutes EURAMET [3] and [4]. Reiss was a pioneer: he suggested using heat pipes to compensate the error induced by the bar effect in the calibration of thermometers [5]. Sostmann [6] and Bienert [7] proposed the use of heat pipes in thermometry; they pointed out that heat pipes are normally designed to transport the maximum heat flow with the minimum possible temperature difference. However, in thermometry applications, the heat transport capacity is less relevant than the device's isothermal behavior. Tamba et al. evaluated a water heat pipe, controlling the operation pressure, when it was used as a comparative furnace in order to calibrate a standard platinum resistance thermometer (SPRT) in a range of 338.15 K to 430.15 K, in a direct comparison against the reference SPRT, the combined uncertainty was estimated as 3.1 mK. These results show that *Corr. Author's Address: National Polytechnic Institute, School of Mechanical and Electrical Engineering Campus Zacatenco, Mexico, fsnchz@yahoo.com.mx 273 his system had a good performance for calibration purposes [8]. Noorma et al. used a water heat pipe as a black body, i.e. as a reference spectral radiator in order to be used as heat source between 323.15 K and 523.15 K; this was adopted by the National Institute of Standards and Technology NIST [9]. In contrast, isothermal furnaces have also been developed on heat pipe bases, but they require 30 minutes to achieve a stable operation temperature; Tasneem has proposed a cylindrical cavity in which the temperature is stabilized. The calibration is done using the cavity as a black body; the cylindrical cavity is filled with graphite in order to make the temperature uniform [10]. In this article, the experimental results of an annular TPT thermal characterization are presented; the goal is to determine its isothermal behavior and the time it takes to achieve this state of operation. The latter is the TPT's most significant characteristic for thermometry applications. 1 THEORY OF THERMOSYPHONS 1.1 Working Principles of Thermosyphons Closed TPT, also known as wickless heat pipe, is a simple but effective two-phase heat transfer device. It performs better in a vertical position, with the evaporator or zone where the heat is supplied at the bottom; in this position, it takes advantage of the gravity acceleration for bringing down the condensed vapor as a liquid film [11]. The working fluid is stored at the TPT bottom as a liquid pool. The container is an evacuated closed pipe filled with a certain amount of a suitable working fluid. The TPT is divided into three main sections; the evaporator section, the adiabatic section, and the condenser section. If heat is added to the evaporator section, the working fluid inside the pipe vaporizes and carries heat from the heat source to the condenser section, where heat is removed by a cooling system. The condensed working fluid returns to the evaporator section, taking advantage of the gravity force, and closing the cycle. Therefore, the condenser section of the TPT must be located above the evaporator section. 1.2 Fundamentals for Calculation In order to fully characterize the TPT, an experimental and analytical study must be undertaken. The condensation can be modeled as taking place inside a vertical pipe but with the interface liquid-vapor strain stress taken into consideration. The convection coefficient in the condenser is normally compared to the one given by the Nusselt theory or against the extended Nusselt theory [11]. In 1916, Nusselt was the first to consider the condensation phenomenon in a vertical plate without stress in the interface. Faghri et al. extended the Nusselt theory to a laminar and turbulent flow inside a thermosyphon, flow reduction due to the condensation rate was evaluated [12]. Spendel [17] also studied the condensation problem in the TPT; the momentum equation in the steam region was numerically solved for the laminar case and film condensing analysis in the condensing region was to coupled it. In contrast, in the specialized literature there are two general approaches that are extensions of Nusselt's film condensation theory and are used to include the variable strain stress in the TPT condenser interface. The first method empirically stablishes the friction factor between the two phases in order to satisfy the mass balance to take into account the axial steam flow. The second numerically solves the steam flow in 2D and couples it with the Nusselt analysis for the film condensation. When TPT is at steady state, the total pressure drop in a closed cycle is the addition of the both phases pressure drop in the different TPT regions. The total balance must be equal to zero [13]. Pg +(Pvc - Pve ) + {Plc - Pvc ) + + (Ple - Plc ) + (Pve - Ple ) = 0. (1) In the above equation, (PVC - PVE ) = APr is the pressure drop in the vapour phase flow, between the condenser and the evaporator. (PLE - PLC ) = APl is the pressure drop in the liquid phase flow between the evaporator and the condenser. (PV - PL ) - (PV - PL )C = 0 represents the pressure difference in the interfaces liquid-vapour in the evaporator region and condenser region respectivelly, so that we can write: Pg + APV +APl = 0. (2) This means that the total pressure drop in the three TPT regions must be equal to zero. Thereofre, the necessary condition for the TPT to function is: Pg >APr +APl . (3) This means that gravity is the system driving force, which must always be equal or superior to the total pressure drop in the system. The pressure drop in the liquid phase flow of the TPT working fluid can be determined using Darcy's equation: AP =-ELmL 1 P L K Al' (4) Because TPT has not a wick on its inner wall, permeability K is a parameter depending on the system characteristics and the working fluid. Moreover, in the TPT, there is a countercurrent annular liquid-vapor flow, so the wavy interface is an important factor for the permeability magnitude. In contrast, f is the TPT effective length and is given by: 1 - l-E +1 + C -eff - + -a + . (5) In order to determine the pressure drop in the vapour phase, the following equation can be applied: 8Mv mV leff (6) щ =- Pv П 4 in which L is the phase change enthalpy given by the expression: L = hv (Ts)-hL (Ts). (10) In order to maintain the TPT steady state operation, the liquid mass flow that evaporates in the evaporator must be equal to the vapour mass flow that condenses in the condenser, that is mL = m V. Next, the liquid mass flow is obtained by adding Eqs. (4), (6) and (8): Pp gl KAL nPv rV leff (ßLnPvrV + 8VVPlKAL ) (11) and the heat flow transported by a TPT is easily computed: Q = (hV - hL )Plgl KALnPvrY leff (ßLnPvrV + 8VVPlKAL ) (12) Gravity effects are given by: Pg = ±Pl gl sind, (7) where в is the TPT inclination angle with respect to the horizontal. In our case, the TPT position is vertical with the evaporator at the bottom so, в = 90° and the TPT gets the maximum driving force: Pg =±PLg1 ■ (8) In contrast, it is well known that the transported heat flow by a TPT depends on the quantity of liquid working fluid that changes phase in the evaporator: Q = m L, (9) 2 EXPERIMENTAL SECTION 2.1 TPT Design According to Faghri [11] and Sanchez et al. [14], the starting point in the TPT design for this kind of application is to know the temperature range in which the device is expected to work. With this information, the working fluid and the TPT body material are selected; both must be chemically compatible; moreover, in thermometry applications, the needed device dimensions are generally known, so it is desirable to use commercial pipes to build the annular TPT. Fig. 1. a) A conventional TPT; b) A two-phase closed annular thermosyphon Based on this information, the heat flux supplied to the evaporator zone can be determined as well as the magnitude of cooling water flow to obtain the stable temperature operation in the system. Fig. 1 shows both a conventional TPT and the annular TPT for thermometry applications. Both the conventional and the annular TPT are pipes sealed on both sides and, using a pump, a vacuum was created before the introduction of the working fluid; the difference is that the annular TPT has an annular cross section in the condenser zone, and the annulus is used as a volume where the thermocouples will be installed for calibration. The annular TPT technical specifications are summarized in Table 1. Table 1. Annular TPT technical specifications Total length 500 mm Evaporator zone length 200 mm Adiabatic zone length 150 mm Condenser zone length 150 mm External diameter of the TPT internal pipe 28.6 mm Internal diameter of the TPT internal pipe 26.2 mm External diameter of the TPT external pipe 54 mm Internal diameter of the TPT external pipe 50.8 mm External diameter of the condenser sleeve 79.4 mm Internal diameter of the condenser sleeve 74 mm Thermal insulation material Mineral wool TPT body material Copper Working fluid Water Volume of the working fluid in the TPT 130 cm3 For the 323.15 K to 523.15 K temperature range, the working fluid must be water. Copper type L was selected as the TPT body material; it withstands a pressure of 35 bar, even though the water saturation pressure at 373.15 K, which is the expected operation temperature, is 101.325 Pa. Both water and copper are chemically compatible so they can work together. The filling volume of working fluid was taken as 15 % of the TPT free volume according to Faghri [11], Sanchez [15] and Carvajal et al. [16]. The annular TPT dimensions are shown in Table 1. 2.2 Experimental System For the experimental system, a 373 W centrifugal pump was selected to guarantee the cooling water circulation in the condenser. The experimental rig has two reservoirs to store cooling water for condensation with a volume of 1 m3 each; both reservoirs are interconnected with a pipe located at the bottom of the tanks. For the heat supply system, a band electrical resistance was used; the resistance has an electrical power of 1,550 W and a resistance of 10.4 Q. Fig. 2 shows a technical scheme of the experimental set up with all the elements. In Fig. 3 the annular TPT, the band type electrical resistance's installation, and the cooling system at the TPT top part are shown (all elements are without thermal insulation). In order to change the power supply, a 120 VCA potentiometer was used. The potentiometer provides Fig. 2. Technical scheme of the experimental set up 0 V to 120 V, and it had connected a voltmeter (±0.05 V accuracy) in order to measure the voltage difference at the exit. A 0 cm3/s to 16.67 cm3/s range rotameter (±0.5 cm3/s accuracy) was used to measure the cooling water flow. The thermocouples used were of the K type (±2.2 °C accuracy in the temperature range studied). They were used to measure the temperature on the surface of the internal pipe wall that is the annular zone, and they were also used to monitor the temperature of the cooling water at the inlet and outlet of the condenser. Fig. 3. The annular TPT experimental system showing the band electrical resistance on the evaporator A digital thermometer (±0.5 °C accuracy in the temperature range studied) with an Omega model HH-22 microprocessor was used for thermocouples type J/K. The thermocouples used in the experiments were located on the internal wall surface of the annular region, all the long its length and with a 3 cm separation between them, as shown in Fig. 4. The annular TPT was thermally insulated with mineral wool and was also covered with a sheet of aluminum. 3 EXPERIMENTAL RESULTS AND DISCUSSIONS Experimental runs were undertaken in order to determine the time in which the annular TPT reaches a working steady state under different operating conditions. Therefore, thermocouples t1, t5 and t8 (see Fig. 4) were monitored during 30 minutes, i.e. the period in which the heating and cooling conditions were maintained at constant levels. Experiments were run for two different cooling water flows (16.67 cm3/s and 4.17 cm3/s), and for each of them the power supplied was varied (38, 154, 350, 1,385) W using the potentiometer to change the applied voltage according to Eq. (13) Ф = —. R (13) 3.1 Experimental Results The experimental results with cooling water flows equal to 16.67 cm3/s and 4.17 cm3/s are shown in Fig. 5a and b. In these figures, it is possible to observe the temperature variation of the thermocouples t1, t5 and t8, located at the bottom, center and at the top of the annular region, respectively (Fig. 4), for three heat flow levels supplied to the evaporator: 154 W, 615 W and 1,385 W. For 16.67 cm3/s, the cooling water flow temperature is slightly bigger for the thermocouple t1, but the difference against t5 and t8 is not more than 5K for the same power level supplied (Fig. 5a). In contrast, for the cooling water flow of 4.17 cm3/s (Fig. 5b), the temperature difference between t1, t5 and t8 is uneven; this means that the small cooling water flow is not enough to make the TPT temperature uniform. Practically, the TPT operation temperature for the small cooling water flow is always higher, and its behavior is uneven, especially in the case of 615 W, in which t1, t5 and t8 show values of (370, 351, 370) K, respectively. In contrast, in Fig. 6a and b, the percentage of the variation in temperatures t1, t5 and t8 in respect with lectures obtained with the same thermocouples, 30 minutes after the beginning of the experiment, is plotted against time. Therefore, the time needed for steady state operation was obtained using the relation: %AT30_i = ( - Ti ) 100 %. (14) 30 Fig. 4. Location of the thermocouples in the TPT annular region Fig. 6a and b show the TPT temperature behavior for the cooling water flows of 16.67 cm3/s, and 4.17 cm3/s, respectively. For larger cooling water flows, there is a smooth tendency to the steady state, which is attained in (18, 28, 18) min for all the power supplied magnitudes, for t1, t5 and t8, respectively (Fig. 6a). In the case of the small cooling water flows, the behavior is more uneven; t1, t5 and t8 attain the steady state in (20, 24, 18) min, respectively (Fig. 6b). For the cooling flow of 4.17 cm3/s, there is a difference of 5 % in the thermocouple t1 for 615 W that vanishes after 28 min of operation; in contrast, the thermocouple t8 shows that the steady state takes place after 18 min. As expected, the steady state is reached faster with larger cooling water flows. To investigate in detail the isothermal behavior, the range of power supply values was broadened including 615 W and 960 W. In Fig. 7, t1 to t8 indicate the thermocouples installed inside the annular region (see Fig. 4). The wall's isothermal behavior inside the annular region is evident, and this condition is easily attained for power supplied in magnitudes greater than 350 W in the evaporator. In this work, isothermicity is defined as such state when the temperature difference is not bigger than 3 % along the TPT surface. It is important to note that the operation temperature values are always higher for the lower cooling water flow. For the power supplied in magnitudes of 38 W and 154 W, the uniformity of the temperature profile is never achieved; therefore, there is an axial temperature gradient, and the annular TPT efficiency is reduced. Practically, for the two cooling water flows used in the experiments (16.67 cm3/s and 4.17 cm3/s) the steady state was obtained 18 min after the experiment a) b) Fig. 5. Temperature behavior measured by thermocouples t1, t5 and t8, for three heating conditions in the evaporator (154, 615,1,385 W) versus time; a) for a cooling water flow of 16.67 cm3/s and b) for a cooling water flow of 4.17 cm3/s Thermocouple 1 I ......154 W--615 W —1385 W 0 2 4 6 S 10 12 14 16 18 20 22 24 26 28 30 t [min] Thermocouple 5_ • 154 W--bl5W -^—1385 W 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 t [min] Thermocouple 1 ■ 154 W--615 W —1385 W ?65 \ [ i 45 о35 f-S25 15 5 -5 3 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 t [min] Thermocouple 5 • 154 W--615 W —1385 75 65 ] &45 О S 35 25 s-?15 5 -5 0 2 4 6 S 10 12 14 16 18 20 22 24 26 28_ 30 t [min] Thermocouple 8 ■ 154 W---615 W 1385 W 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 _ t [min] a) 75 R i 45 35 e-?25 15 5 -5 Thermocouple 8 ■ 154 W--615 W -—1385 W £ 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 t [min] b) Fig. 6. Percent variation of temperatures t1, t5 and t8 for three heat flow levels; a) for a cooling water flow of 16.67 cm3/s and b) for 4.17 cm3/s T [K] —-3BW---154 W 350 W -- 615W -^960 W ....... 1385 W | 380 370 360 350 340 330 320 310 300 290 t1 t2 t3 t4 t5 t6 t7 T [K] I —*—3SW---154 W —*-350 W -- 615W -*-960W | 390 370 350 330 310 290 fS t1 t2 t3 t4 t5 t6 t7 t8 a) b) Fig. 7. Temperature values measured for the different thermocouples t1 to t8 for 38,154, 350, 615, 960 and 1385 W of heat flow levels supplied, with cooling water flow of a) 16.67 cm3/s and b) 4.17 cm3/s, 30 minutes after the experiment start started (heat flow supply and cooling water flow maintained constant). For this period, there was a variation in temperature values of 2 % in all analyzed cases. 3.2 Results Discussion One of the objectives of this research was to evaluate the annular TPT as an isothermal source to calibrate temperature instruments, using the comparison method. In the experimental test, it was proved that TPT operation is isothermal along the 21 cm characterized annular zone length; consequently, the bar effect that produces a temperature measurement error, due to the axial temperature gradient, was avoided. Therefore, the device fulfills the requirements to be used as an isothermal reservoir for calibration purposes using the comparative method in the range from 313.15 K to 373.15 K. Of course, when the required range of application is larger, the TPT working fluid should be changed, as well the cooling fluid and its flow. In this research, the annular TPT immersion depth was 21 cm, which can be easily extended by changing the device characteristics. For the small heat flow magnitudes (38 W and 154 W), the heat flow supplied to evaporator is not sufficient to guarantee that the axial pressure gradient is near zero and consequently the axial isothermal behavior is not reached. Nevertheless, for heat flow supplied magnitudes higher than 350 W this effect disappears, as shown in Fig. 7. However, it is important to note that the annular TPT region operates at higher temperatures when the water cooling flow is reduced. The other goal of this research is to propose the use of the TPT as a secondary reference for temperature instrument calibration. However, this application implies a guarantee of the TPT stability and repeatability for different TPT working conditions; therefore, there should be the same values of temperature in different schedules using the same operation conditions. Figs. 8a and b show the temperature profile on the annular region wall of the TPT, for three distinct experimental tests carried out in different days, but keeping the same operation condition (960 W and 1,385 W, with 16.67 cm3/s of cooling water flow), after 30 minutes of the experiment start. In the above-mentioned graphs, it obvious that the temperature values are not strictly the same, and the differences were approximately of 4 K (i.e. 2.9 %): therefore, strictly speaking, the results are not reproducible. This difference could be attributable to the cooling water temperature, which was not constant, because it depends on environmental conditions. In contrast, the other involved parameters (e.g. the electrical power supply and the cooling water flow) need to be controlled with a higher degree of precision. In contrast, considering that the measured values obtained in different instants are not the same, it was not possible to obtain a calibration curve for the annular TPT in order to use it as a secondary reference for the calibration of temperature instruments. In order to obtain a reliable TPT calibration curve and propose its use as a secondary reference for metrology purposes, it is necessary to keep constant, precise and stable the following parameters: • Inlet cooling water temperature; • Cooling water flow; • Heat flow supplied magnitude in the evaporator. Furthermore, even though the TPT is thermally insulated, it is also important to avoid any stream of air in the room where the experiments take place to avoid changing the surrounding conditions. -27/07/00 1/08/00 26/07/00 T [K] T [K] 378 377 376 375 374 373 372 371 370 369 368 ■ 27/07/00 1/08/00 - 26/07/00 tl t2 t3 t4 t5 t6 t7 t8 a) b) Fig. 8. Repeatability of the axial temperature profile. Replay of the same experimental condition; a) 960 W and b) 1,385 W with 16.67 cm3/s of cooling water flow, carried out in different days and on a different schedule In order to evaluate the radial heat leakages, the thermal TPT efficiency was computed; this parameter was expected to be high because a TPT is similar to an ideal fin. The abovementioned is shown in Fig. 9, where the TPT efficiency is plotted at different operation power magnitudes, and was always higher than 94 %. П = Qk Ф (15) Therefore, the device fulfills the requirements of an isothermal reservoir for calibration purposes using the comparison method. The range of the device's stable operation can be broadened by making some changes as previously mentioned. In contrast, if the variations of the temperature are reduced by controlling the main parameters affecting the TPT stable operation and their precision, it would be possible to use the annular TPT as a secondary reference for calibration purposes. The mean uncertainty of measured temperature in the range studied was ± 6 %. i 0.95 0.9 0.85 0.8 0.75 0.7 0.65 0.6 0.55 0.5 • ••4.17 :m3/s 7 cm/s • [W] Fig. 9. TPT Efficiency for different operation conditions Table 2. Uncertainty of voltage V and power supplied Ф Voltage [V] uv [V] иф [%] 40 ±0.00124 ±0.514 60 ±0.000814 ±0.508 80 ±0.000598 ±0.496 100 ±0.000482 ±0.490 120 ±0.000410 ±0.488 Table 3. Uncertainty of heat flow transported, Q and efficiency ц Voltage [V] u0 [%] Un [%] 40 ±4.20 ±4.23 60 ±2.88 ±2.92 80 ±2.61 ±2.66 100 ±2.55 ±2.60 120 ±2.53 ±2.58 The maximum uncertainties obtained for main variables calculated in the present study is uncertainty of electrical resistance, R: uR = ±0.00481 a 4 CONCLUSION An annular two-phase thermosyphon was designed, built and characterized, in order to evaluate its use in thermometry. Experiments were carried out for 16.67 cm3/s and 4.17 cm3/s cooling water flows with heat flow supplied magnitudes in the evaporator of (38, 154, 350, 615, 960, 1,385) W for each cooling water flow. Except in the cases of 38 W and 154 W, the isothermal TPT characteristic was observed throughout the annular 21 cm length. The system accomplishes all the requirements for an isothermal source for calibration of temperature instruments using the comparative method. The bar effect does not appear or is reduced. Therefore, temperature measurements do not need to be corrected to take in account this thermal effect. According to isothermal behavior shown by TPT, it can be proposed as an isothermal source for the calibration of temperature instruments using the comparative method. Another advantage is the immersion depth of the temperature sensor, which is 21 cm, in this case. This value could be increased by changing the geometry of the annular TPT. Regarding the repeatability of experiments, the cooling water temperature was not constant, because it depends on the surrounding conditions and other parameters involved, which were not fully controlled in the experiments (voltage, electric supply to the resistance terminals and cooling water flow). It was, thus, not possible to obtain the device calibration curve. Therefore, for experiments with the same operation conditions (heat flow supplied to evaporator and cooling water flow in the condenser), but in different schedules and days, there was not a good reproduction of the annular region temperature where the steady state was attained. Under these conditions, it is not possible to use the annular TPT as a secondary reference for calibration purposes. 5 ACKNOWLEDGMENT The authors would like to thank the National Polytechnic Institute of Mexico and the National Science and Technology Agency of Mexico (CONACYT) for funding this research. П 6 NOMENCLATURE A area [m2] K permeability [m2] L change phase specific enthalpy [J/kg] P pressure [Pa] Q heat flow transported [W] R electrical resistance [Q] T temperature [K] V electrical potential [V] g gravity acceleration [m/s2] h specific enthalpy [J/kg] l length [m] m mass flow rate [kg/s] r radius [m] t time [s] u uncertainty [-] Greek A difference [-] Ф heat power supplied [W] n efficiency [-] dynamic viscosity [Paxs] P density [kg/m3] q inclination angle [°] Sub-index a adiabatic eff effective g gravity i instant s saturation C condenser E evaporator L liquid LC liquid in the condenser LE liquid in the evaporator Q heat flow transported R resistance S saturation V vapor, voltage VC vapor in the condenser VE vapor in the evaporator 30 after 30 minutes of exposition 30-i difference between 30 minutes and the time n efficiency 7 REFERENCES [1] ASTM E 220-13, (2011). Standard method for calibration of thermocouples by comparison techniques. American Society for Testing and Materials, Philadelphia, D0l:10.1520/E0220. [2] Mills, A.F. (1999). Heat Transfer, 2nd ed., Prentice Hall, Homewood, Boston. [3] EURAMET (2011). Calibration of Thermocouples, EURAMET cg-8, Version 2.1., from: http://www.euramet.org/fileadmin/ docs/Publications/calguides/EURAMET_cg-8__v_2.1_ Calibration_of_Thermocouples.pdf, accessed on 2015-03-10. [4] Bernier, M.A., Baliga, B.R. (1993). Effects of heat losses (or gains) from insulated portions of closed loop thermosyphons with vertical heat transfer sections. Journal of Heat Transfer, vol. 115, no. 2, p. 480-482, D0I:10.1115/1.2910705. [5] Reiss, F.E. (1974). Heat pipe as a means of avoiding emergent stem correction for liquid in glass Thermometers. Review of Scientific Instruments, vol. 45, no. 9, p. 1157-1158, D0I:10.1063/1.1686831. [6] Sostmann, H.E. (1991). The heat pipe and its use in thermometer calibration. ISOTECH Journal of Thermometry, vol. 2, no. 1, p. 31. [7] Bienert, W.W. (1991). Isothermal heat pipes and pressure controlled furnaces. ISOTECH Journal of Thermometry, vol. 2, no. 1, p. 32-49. [8] Tamba, J., Kishimoto, I., Arai, M., (2002). Evaluation of a Pressure-Controlled Water Heat Pipe for Accurate Comparison of Platinum Resistance Thermometers. Proceedings of the 41st Society of Instrument and Control Engineers, Osaka, vol. 4, p. 2435-2438, D0I:10.1109/SICE.2002.1195792. [9] Noorma, M., Mekhontsev, S., Khromchenko, V., Litorja, M., Cagran, C., Zeng, J., Hanssen, L. (2006). Water heat pipe blackbody as a reference spectral radiance source between 50°C and 250°C. Proceedings of the International Conference on Thermal Sensing and Imaging Diagnostic Applications, vol. 6205, D0I:10.1117/12.667257. [10] Tasneem, K.A. (2009). Calibrating a thermocouple and investigating heat transfer. Technical Report: Heat transfer and thermocouples. Lahore University of Management Sciences, Lahore, from http://lums.edu.pk/docs/Calibration. pdf, accessed on 2015-03-10. [11] Faghri, A., (1995). Heat Pipe Science and Technology. Taylor & Francis Group, London. [12] Faghri, A., Chen, M.M., Morgan, M. (1989). Heat transfer characteristics in two-phase closed conventional and concentric annular thermosyphons. Journal of Heat Transfer, vol. 111, no. 3, p. 611-618, D0I:10.1115/1.3250726. [13] Chi, S.W. (1976). Heat Pipe Theory and Practice, McGraw-Hill, New York. [14] Sanchez-Silva, F., Carvajal Mariscal, I., Tolentino Eslava, G., Abugaber Francis, J. (2001). Thermal and mechanical design of a two phase flow thermosyphon heat exchanger. Información Tecnológica, vol. 12, no. 3, p. 63-70. (in Spanish) [15] Sanchez-Silva, F. (1980). Contribution to the study of two-phase water thermosyphons, Ph.D. Thesis, Université de Poitiers, Poitiers. (in French) [16] Carvajal-Mariscal, I., Sanchez-Silva, F., Polupan G., (2012). Development of High Efficiency Two-Phase Thermosyphons for Heat Recovery. Mitrovic, J. (ed.). Heat Exchangers - Basics Design Applications, p. 97-116, D0I:10.5772/32693. [17] Spendel, T. (1984). Laminar film condensation heat transfer in closed two-phase thermosyphons. Proceedings of 5th International Heat Pipe Conference, p. 208-213. Vsebina Strojniški vestnik - Journal of Mechanical Engineering letnik 61, (2015), številka 4 Ljubljana, april 2015 ISSN 0039-2480 Izhaja mesečno Razširjeni povzetki Špela Bolka, Janko Slavič, Miha Boltežar: Identifikacija materialnih lastnosti pločevine v prečni smeri glede na ravnino valjanja s procesom prebijanja SI 43 Dario Croccolo, Massimiliano De Agostinis, Stefano Fini, Alfredo Liverani, Nicolò Marinelli, Eugenio Nisini, Giorgio Olmi: Mehanske lastnosti dveh okolju prijaznih smol, ojačenih z lanenimi vlakni SI 44 Yabing Cheng, Yang Wang, Lei Li, Shuaibing Yin, Lichi An, Xiaopeng Wang: Postopek snovanja dvofaznega tihega verižnega prenosnika Hy-Vo SI 45 Franci Čuš, Uroš Župerl: Simulator sistema vodenja hrapavosti površine pri procesu struženja SI 46 Janez Lupše, Leopold Škerget, Jure Ravnik: Algoritem na osnovi metode robnih elementov za URANS simulacije toka preko kvadratne ovire SI 47 Wang Jun, Sun Binyu, Huang Qingxue, Li Hongzhou, Han Heyong: Raziskava glavnega/podrejenega krmiljenja položaja/tlaka za servohidravlični sistem potujočih škarij SI 48 Florencio Sanchez-Silva, Ignacio Carvajal-Mariscal, Ariel E. Moreno-Cordobés, Pedro Quinto Diez, Miguel Toledo Velazquez: Študija obročastega dvofaznega termosifona kot izotermnega vira pri merjenju temperature SI 49 Osebne vesti Doktorske disertacije, magistrska dela, diplomske naloge SI 50 Identifikacija materialnih lastnosti pločevine v prečni smeri glede na ravnino valjanja s procesom prebijanja Špela Bolka1 - Janko Slavič2 - Miha Boltežar2* 1 Hidria Rotomatika, Slovenia 2 Univerza v Ljubljani, Fakulteta za strojništvo, Slovenija Karakterizacija materialnih lastnosti pločevine se tipično izvaja v ravnini valjanja pločevine s pomočjo standardiziranih testov, denimo enoosnega nateznega ali strižnega preizkusa. Pločevino za različne aplikacije, na primer prebijanje ali globoki vlek, se praviloma dobavlja v večjih kolutih. Materialne lastnosti se spreminjajo že v posameznem kolutu, zagotovo pa med posameznimi koluti, kar otežuje tako proizvodnjo kvalitetnih izdelkov, kot tudi eksperimentalno karakterizacijo materiala. Poleg tega se preoblikovanje pločevine pretežno izvaja v prečni smeri glede na ravnino valjanja pločevine, zato bi bilo smotrno spreminjanje materialnih lastnosti spremljati v tej smeri, vrednosti identificiranih parametrov pa primerjati z vrednostmi, dobljenimi na podlagi standardnih enoosnih nateznih testov v ravnini valjanja pločevine. V prispevku obravnavamo pristop k eksperimentalni karakterizaciji mehanskih lastnosti pločevine v prečni smeri glede na ravnino valjanja pločevine. Za eksperimentalno identifikacijo je uporabljena na novo razvita laboratorijska naprava za prebijanje pločevine, ki omogoča polne, delne (do določene globine) ter zaporedne (zaporedni delni preboji do polnega pretrga materiala) preboje pločevine pri različnih hitrostih procesa. Naprava omogoča veliko ponovljivost in natančnost meritev. Mogoča je uporaba različnih geometrij rezilnih elementov -v obravnavanem delu so bili uporabljeni rezilni elementi krožne oblike nazivnega premera 10 mm ter zračnosti 0.025 mm na polmer. Polni preboji so bili izvedeni kvazistatično s hitrostjo približno 0.1 mm/s ter dinamično s hitrostmi (200, 250, 300, 350) mm/s. Enako so bili kvazistatično in dinamično izvedeni tudi delni preboji do globin (100, 200, 300, 350) ^m ter zaporedni preboji s korakom 50 ^m oziroma 100 ^m. Pri tem so bili delni in zaporedni kvazistatični preboji izvedeni z enako hitrostjo kot polni, pri dinamičnih eksperimentih pa je bila začetna hitrost prebijanja ravno dovolj velika, da je nož prebil material do želene globine. Za vse eksperimente je bila uporabljena končno-žarjena elektro pločevina debeline 0.5 mm, eksperimenti pa so bili izvedeni brez dodatnega mazanja in pri sobnih temperaturah. Eksperimenti polnih prebojev so služili kot osnova za identifikacijo strižne meje tečenja in strižne trdnosti materiala, ki jihje bilo nato mogoče primerjati z natezno mejo tečenja in natezno trdnostjo materiala v smeri valjanja pločevine. Z uporabo delnih in zaporednih prebojev pa je bilo mogoče iz testa prebijanja identificirati akumulacijo poškodovanosti materiala. Tudi te rezultate je bilo moč primerjati z rezultati nateznih testov - s pomočjo koračnih enoosnih nateznih testov je mogoče, z opazovanjem spreminjanja naklona krivulje razbremenjevanja, določiti napredovanje parametra poškodovanosti, popisanega z Lemaitreovim izotropnim zakonom poškodovanosti. Enak pristop je bil uveden tudi v prečni smeri glede na ravnino valjanja pločevine. Predstavljena metoda tako omogoča nov način karakterizacije in neposrednega spremljanja materialnih lastnosti iz samega procesa prebijanja, kar nakazuje možnost širše industrijske uporabe za sprotno identifikacijo materialnih lastnosti. Ključne besede: rezanje in preoblikovanje, konstitutivne relacije, elasto-plastični material, mehanski testi, žilavi lom, eksperiment v smeri, prečni na ravnino valjanja *Naslov avtorja za dopisovanje: Univerza v Ljubljani, Fakulteta za strojništvo, Aškerčeva 6, 1000 Ljubljana, Slovenija, miha.boltezar@fs.uni-lj.si SI 43 Mehanske lastnosti dveh okolju prijaznih smol, ojačenih z lanenimi vlakni Dario Croccolo* - Massimiliano De Agostinis - Stefano Fini - Alfredo Liverani -Nicolò Marinelli - Eugenio Nisini - Giorgio Olmi Univerza v Bologni, Oddelek za industrijski inženiring (DIN), Italija Članek podaja primerjavo mehanskih lastnosti dveh različnih kompozitnih materialov z naravnimi vlakni. Kompozita sta izdelana iz dveh različnih okolju prijaznih smol in ojačena z enakimi lanenimi vlakni. Prva (181 EN2X) je smola na bazi izoftalne kisline s 25 % polietilen tereftalata, pridobljenega iz recikliranega materiala. Druga (VEef220ST) je vinil ester, pridobljen iz bisfenolne epoksi smole z majhno vsebnostjo stirena. Ti dve smoli sta se uveljavili v pomorstvu zaradi visoke obstojnosti proti vplivom okolja, še posebej proti korozivnemu učinku slane vode. Mehanske lastnosti kompozitov so bile preiskane z nateznimi, strižnimi in upogibnimi preizkusi po veljavnih standardih (ISO 527-1 in -4, ISO 14125, ISO 14130). Eksperimenti so bili opravljeni na stiskalnici Instron 8033, opremljeni s 25 kN obremenitveno celico, nastavljivimi prijemali za natezni preizkus in provizorično podporo za tritočkovno obremenitev pri upogibnih in strižnih preizkusih, skladno z ISO 14125 in ISO 14130. Preizkušenih je bilo skupno 32 vzorcev, pripravljenih po omenjenih standardih. Pri nateznih in upogibnih preizkusih so bili pridobljeni sprejemljivi rezultati, kar pa ne velja za strižne preizkuse, kjer so vzorci odpovedali na nesprejemljiv način - verjetno zaradi lastnosti laminiranega kompozita. Rezultati kažejo, da je izoftalna smola boljša od vinil estrske smole (VEef220ST) tako po trdnosti, kakor tudi po togosti. Ojačitev poveča natezno trdnost smole 181 EN2X za približno 25,4 %, modul elastičnosti za približno 21,5 %, upogibno trdnost za približno 8,7 % in upogibni modul elastičnosti za približno 8,2 %. Pri vinil estrski smoli je ugoden učinek ojačitve nekoliko manj izražen, predvsem pri natezni trdnosti. V diskusiji je predstavljena primerjava s podobnimi študijami v dostopni literaturi, ki obravnavajo kompozite, ojačene z rastlinskimi vlakni. Omeniti je treba, da je učinek ojačitve z lanenimi vlakni pri okolju prijaznih smolah primerljiv z učinkom ojačitve pri kompozitih, izdelanih s konvencionalnimi epoksi smolami. Za oceno adhezije med ojačitvijo in matriksom ter prisotnosti napak je bilo nekaj vzorcev iz izoftalne smole razrezanih in pregledanih pod optičnim mikroskopom. Analiza je pokazala, da je osnovna nit premera 500 ^m do 600 ^m, medtem ko je votkovna nit premera 200 ^m do 250 ^m. Zračnih vključkov je manj in so manjši v primerjavi s podatki iz literature za kompozite, iz katerih se izdelujejo trupi plovil. Analiza površine zloma pri vzorcih nateznega preizkusa je pokazala homogeno in pravilno površino, kar pomeni, da matriks in ojačitev skupaj prispevata k celotni natezni trdnosti kompozita. Takšen rezultat je dokaz za dobro adhezijo med komponentama. Najpomembnejši sklep te študije je, da se mehanske lastnosti izoftalne smole občutno izboljšajo s prisotnostjo vlaknene ojačitve. Pomembna lastnost te smole je, da je zelo občutljiva na prisotnost lanenih vlaken. Nasprotno pa dodajanje teh vlaken ni nujno koristno tudi pri matriksu iz vinil estrske smole. Ključne besede: kompoziti z naravnimi vlakni, mehanske lastnosti, lan, okolju prijazna smola, napetosti, deformacije SI 44 *Naslov avtorja za dopisovanje: Univerza v Bologni, Oddelek za industrijski inženiring (DIN), Viale del Risorgimento 2, 40136 Bologna, Italija, dario.croccolo@unibo.it Postopek snovanja dvofaznega tihega verižnega prenosnika Hy-Vo Yabing Cheng* - Yang Wang - Lei Li - Shuaibing Yin - Lichi An - Xiaopeng Wang Univerza Jilin, Fakulteta za tehniške vede in strojništvo, Kitajska Dvofazni tihi verižni prenosni sistem Hy-Vo je zasnovan za zmanjšanje fluktuacij in poligonskega efekta pri pogonskih sistemih. Članek podaja postopek snovanja dvofaznega tihega prenosnega sistema s ploščicami, verižnikom in orodjem kot teoretično osnovo za projektiranje drugih večfaznih verižnih prenosnikov. Predstavljen je postopek snovanja dvofaznega tihega verižnega prenosnika Hy-Vo. Dvofazni verižnik ima za razliko od enofaznih prenosnikov dva razdelka in dve skupini kotov pritiska. Zaradi fazne razlike med zobmi dvofaznega verižnika imajo lahko zobje obeh faz različne razdelke in kote pritiska. Postavljen je model dvofaznega tihega verižnega prenosnika Hy-Vo in dvofazni prenosnik je razstavljen v dva enofazna prenosna sistema. Poligonski efekt dvofaznega prenosnega sistema, ki ga opisuje veličina fluktuacije, je analiziran matematično in s simulacijo. Članek preučuje metodo snovanja dvofaznega tihega verižnega prenosnika Hy-Vo ter poligonski efekt pri dvofaznem prenosniku. Na podlagi teorije ubiranja tihe verige Hy-Vo z verižnikom in povezanega pogoja ubiranja tihega verižnika Hy-Vo z orodjem sta zasnovana multivariabilna tiha veriga in verižnik Hy-Vo ter je postavljen sistem ubiranja dvofazne tihe verige Hy-Vo. Razlika med fluktuacijami pri dvofaznih in enofaznih prenosnih sistemih je bila analizirana s Fourierjevo vrsto, za oceno poligonskega efekta pa je bila uporabljena simulacija. Sistem za snovanje ubiranja med tiho verigo, verižnikom in orodjem Hy-Vo je primeren tako za snovanje enofaznih tihih verižnih prenosnih sistemov Hy-Vo, kakor tudi za snovanje dvofaznih tihih prenosnih sistemov Hy-Vo. Metoda snovanja ploščic tihe verige, verižnika in orodja dvofaznega sistema Hy-Vo je v tem primeru novost. Če sta razdelka dvofaznega verižnika enaka in sta enaki tudi skupini kotov pritiska, je lahko trenutno prenosno razmerje dvofaznega prenosnega sistema drugačno od trenutnega prenosnega razmerja razstavljenega neodvisnega enofaznega prenosnega sistema. Največja fluktuacija pri dvofaznem prenosnem sistemu je manjša kot pri enofaznem sistemu, manjši pa je tudi poligonski efekt. Članek preučuje metodo snovanja dvofaznega tihega verižnega prenosnika Hy-Vo ter poligonski efekt pri dvofaznem prenosniku. V prihodnjih raziskavah bi bilo mogoče preučiti tudi hrup in napetosti dvofaznega tihega verižnega prenosnika Hy-Vo, kakor tudi sistem večfazne verige. Predstavljen je nov postopek snovanja dvofaznega tihega verižnega prenosnika Hy-Vo. Uporaba tihe verige in dvofaznega verižnika Hy-Vo v prenosnih sistemih, npr. za krmiljenje odmične gredi, oljne črpalke in izravnalne gredi, lahko zmanjša hrup ter pripomore k bolj gladkemu delovanju prenosnega sistema. Ključne besede: dvofazni prenosnik, tiha veriga Hy-Vo, multivariabilnost, snovanje sistema ubiranja, veličina fluktuacij, poligonski efekt *Naslov avtorja za dopisovanje: Univerza Jilin, Fakulteta za tehniške vede in strojništvo, Renmin street 5988, Changchun, Kitajska, chengyb@jlu.edu.cn SI 45 Simulator sistema vodenja hrapavosti površine pri procesu struženja Franci Čuš - Uroš Župerl* Univerza v Mariboru, Fakulteta za strojništvo, Slovenija Cilj članka je predstaviti modelno podprt sistem vodenja struženja s pripadajočim simulacijskim blokovnim diagramom. Namen na modelu temelječega sistema vodenja je zagotavljati želeno hrapavost obdelane površine s sprotnim prilagajanjem rezalnih parametrov in ohranjanjem konstante rezalne sile. Sistem z dinamično adaptacijo podajanja obvladuje hrapavost površine ter rezalne sile na struženem nožu. Namen izdelave predlaganega sistema vodenja je poiskati omejitve takšnega principa vodenja. Modelno podprt sistem vodenja je izdelan na osnovi poenostavljenega modela procesa struženja in podajalnega servo-sistema. Celotni sistem vodenja se sesoji iz modela napovedovanja rezalnih pogojev, referenčnega bloka, bloka vodenja, procesa odrezavanja, CNC servosistema, elementa za korekcijo podajanja in senzorja, ki meri izhodno veličino procesa. Glavna predpostavka raziskave je bila, da signali izmerjenih rezalnih sil vsebujejo največ uporabnih informacij o kakovosti obdelane površine. Zato so izvedeni eksperimenti, s katerimi so določene funkcijske korelacije med hrapavostjo obdelane površine, rezalnimi silami in parametri vodenja. Korelacija med procesnima izhodnima veličinama in rezalnimi parametri je določena z regresijo. Pri vnaprej definirani globini rezanja je eksperimentalno raziskan vpliv obdelovalnega materiala in rezalnih parametrov (podajanje, globina rezanja) na omenjeno korelacijo. Ugotovljena je visoka stopnja ujemanja med hrapavostjo površine in rezultirajočo rezalno silo pri različnih rezalnih pogojih. Te povezave se nato uporabljene pri izdelavi modelno podprtega sistema vodenja v katerem se z reguliranjem rezalnih sil izboljšuje zahtevana kakovost obdelane površine. Z namenom demonstracije učinkovitosti predlaganega sistema vodenja je predstavljena zgradba simulatorja in simulacijski rezultati. Model vodenja je izdelan v simulacijskem paketu Matlab-Simulink. Izpeljana je prenosna funkcija celotnega procesa struženja. Rezultati simulacij pokažejo, da je z vidika fluktuacij hrapavosti površine in učinkovitosti struženja predlagan sistem vodenja precej bolj učinkovit kot sistem z konvencionalnim CNC krmiljem. S simulacijami je potrjena adekvatnost in stabilnost sistema vodenja. Rezultati zagotovijo načine za povečanje učinkovitosti procesa z izboljšanjem kakovosti površine, minimiziranjem posledic spremenljivosti procesa in zmanjšanjem stroškov napak pri operacijah končne obdelave. Dokazano je, da lahko z obvladovanjem rezalnih sil uspešno nadzorujemo hrapavost površine, ki je bistven kazalec kakovosti procesa. Z ohranjanjem konstantne rezalne sile se zagotovi konstantna kakovost obdelane površine. Po uspešno izvedenih simulacijah je bil simulacijski model CNC servo-sistema in simulator odrezavanja nadomeščen z obdelovalnim strojem GF NDM 16. Stružnica je vodena s CNC krmiljem, katero omogoča korekcije podajalne hitrosti v območju od 0 % do 150 %. Povezana je s osebnim računalnikom na katerem je implementiran sitem vodenja. Zmogljivosti predlaganega sistema so testirane s serijo eksperimentov na Ck45 jeklenih testnih preizkušancih pri spremenljivih rezalnih pogojih. Z integracijo izdelanega sistema vodenja s CNC krmiljem stružnice se poveča kakovost obdelanih komponent. Obdelovalni sistem postane neobčutljiv na motnje zaradi obrabe orodja, nehomogenosti obdelovalnega materiala ter variacij globine struženja. Sistem je zasnovan za preoces struženja, vendar ga je možno prilagoditi za vse postopke obdelave z odrezavanjem. Ključne besede: odrezavanje, hrapavost površine, modelno podprto vodenje, simulator SI 46 *Naslov avtorja za dopisovanje: Univerza v Mariboru, Fakulteta za strojništvo, Smetanova ulica 17, 2000 Maribor, Slovenija, uros.zuperl@um.si Algoritem na osnovi metode robnih elementov za URANS simulacije toka preko kvadratne ovire Janez Lupše* - Leopold Škerget - Jure Ravnik Univerza v Mariboru, Fakulteta za strojništvo, Slovenija Danes so numerične metode za preračun toka tekočin v industriji v vsakodnevni uporabi. Vendar pa je večina industrijsko pomembnih tokovnih problemov računsko zahtevnih, kar je povezano z velikimi računskimi zmogljivostmi in s tem veliko porabo časa in energije. Večina današnjih metod v uporabi deluje na principu »surove sile«, kjer z zelo gostimi numeričnimi mrežami opisujemo fizikalno dinamiko pojavov. Glavna prednost metode robnih elementov (MRE), v primerjavi z večino ostalih numeričnih metod, je vključitev fizikalnih mehanizmov toka že med samim postopkom diskretizacije, kar vodi v zmanjšanje potreb po gosti numerični mreži, hkrati pa poveča stabilnost numeričnega algoritma. V članku predstavljamo numerični algoritem, osnovan na metodi robnih elementov, za računanje nestacionarnih tokov v direktnem in URANS načinu. Algoritem je osnovan na hitrostno-vrtinčnem zapisu Navier-Stokesovih enačb. Tak zapis uporabljajo tudi drugi najbolj učinkoviti algoritmi za izvajanje direktnih numeričnih simulacij. Pomemben del algoritma je določitev robnih vrednosti vrtinčnosti, kjer se uporaba MRE izkaže. Izračun območnih vrednosti mora biti opravljen globalno na celotni domeni reševanja, kar pripelje do reševanja polnih sistemov enačb in s tem velike uporabe računalniškega spomina ter posledično praktičnih omejitev pri velikosti numeričnih mrež. Delno se da problem omiliti z raznoraznimi metodami stiskanja polne matrike, saj se izkaže, da imajo členi, ki ležijo daleč stran od diagonale, zanemarljiv vpliv na natančnost rešitve, ali pa z uporabo drugačne oblike vodilnih enačb. Zaradi lastnosti Greenovih funkcij, uporabljenih kot utežne funkcije v obravnavani metodi, je možno uporabiti relativno redke numerične mreže v primerjavi z drugimi metodami. Hkrati so del rešitve tudi normalni odvodi funkcij polja, kar zmanjša napake pri računanju le-teh in je pomembno npr. pri natančni določitvi Nusseltovega števila. Algoritem je nadgradnja našega algoritma, razvitega za RANS modele z integracijo do zida. Ker so ti modeli nezadostni v primeru, ko je tokovni problem že v sami osnovi nestacionaren, smo v predloženem članku algoritem razširili še za URANS simulacije. Preverjanje algoritma je izvedeno na primeru toka okoli kvadratne ovire v ravnem kanalu, pri čemer je bilo razmerje med širino ovire in širino kanala 1:8. Tok postane nestacionaren že pri Reynoldsovem številu med 54 in 70, kar omogoča preverjanje algoritma že v laminarnem režimu. V laminarnem režimu smo izvedli numerične simulacije pri Reynoldsovih številih 100, 150, 200 in 250 ter v nizko turbulentnem režimu pri Reynoldsovem številu 1000. Simulacije laminarnega toka pokažejo dobro ujemanje rezultatov z rezultati izračunov iz literature za komponento hitrosti v smeri toka in tudi za Strouhalovo število ter koeficient upora. V nizko turbulentnem režimu obravnavamo tok s pomočjo kvazi-direktne numerične simulacije in s Spalart-Allmarasovim modelom turbulence v URANS načinu, kjer primerjamo rezultate za oba integralna parametra z rezultati iz literature. Zaradi dušenja toka z modelirano viskoznostjo v URANS simulaciji ne vidimo manjših struktur v toku, kar se pozna tudi pri obeh izračunanih parametrih, tok pa postane navidezno podoben laminarnemu. V nasprotju pa pri kvazidirektni simulaciji tega dušenja nimamo in so simulirane strukture omejene z velikostjo mreže, kar se kaže v veliko bolj bogatem spektru struktur v toku. Ključne besede: metoda robnih elementov, numerični algoritem, hitrostno-vrtinčna oblika, turbulentni tok nestisljive tekočine, URANS, kvadratna ovira *Naslov avtorja za dopisovanje: Univerza v Mariboru, Fakulteta za strojništvo, Smetanova ulica 17, 2000 Maribor, Slovenija, janez.lupse@yahoo.ca SI 47 Raziskava glavnega/podrejenega krmiljenja položaja/tlaka za servohidravlični sistem potujočih škarij Wang Jun1 - Sun Binyu1* - Huang Qingxue1 - Li Hongzhou2 - Han Heyong1 1 Tniverza za znanost in tehnologijo v Taiyuanu, Kitajska 2 Hebei WenFeng železo in jeklo, Kitajska Na univerzi za znanost in tehnologijo v Taiyuanu so bile razvite potujoče hidravlične škarje, ki se uporabljajo za rezanje jeklenih plošč debeline do 35 mm. Če debelina jeklenih plošč presega 35 mm, imajo potujoče škarje težave z rezanjem in podatki kažejo, da ni mogoče dovolj hitro vzpostaviti tlaka. Razlog je v tem, da škarje uporabljajo le zaprtozančno krmiljenje položaja in tlak ni pod nadzorom. Za razrešitev problema predlagamo novo funkcijo za glavno/podrejeno krmiljenje položaja/tlaka, ki lahko izboljša točnost položaja in tlaka. Zahteve krmiljenja položaja in tlaka so si navidezno v nasprotju, zato uspešno krmiljenje obeh parametrov predstavlja poseben izziv. Nova krmilna funkcija je zasnovana kot dvojen zaprtozančni sistem, kjer ima glavno vlogo krmiljenje položaja in sekundarno vlogo krmiljenje tlaka. Tlačni signal iz valja zaznajo tlačni senzorji in ga pošljejo v model za pretvorbo tlaka in položaja. Krmilnik nato presodi, ali se tlačni signal doda položajnemu signalu. Če je napaka položajnega signala večja od praga, se priključi zaprta zanka za krmiljenje tlaka in tlačni signal se pretvori v kompenzacijski signal za prilagoditev vhodnega položajnega signala. Pretok servoventila se hitro poveča in valj postane hiter, tako da je napaka položaja valja majhna, tlak in položaj valja pa sta pod nadzorom. Krmiljenje položaja ne vpliva na krmiljenje tlaka. Z novim pristopom je mogoče istočasno krmiliti položaj in tlak valja. Glavno/podrejeno krmiljenje položaja/tlaka, ki je opisano v članku, izboljšuje varnost pred zatikanjem, natančnost in stabilnost sistema. Debelina rezanih jeklenih plošč je lahko večja za 35 %, izboljša pa se tudi kakovost. Krmilni sistem je mogoče enostavno vgraditi v omenjene in druge hidravlične servosisteme z visokimi hitrostmi, velikimi pretoki in velikimi bremeni. Članek podaja matematični model glavnega/podrejenega krmiljenja položaja/tlaka za servohidravlični sistem. Sistem je bil simuliran s programsko opremo AMEsim in Matlab, pri čemer je bila potrjena uporabnost pristopa. Končno so bili opravljeni tudi poizkusi na eksperimentalni opremi z novo funkcijo. Rezultati eksperimentov kažejo, da lahko glavno/podrejeno krmiljenje položaja/tlaka izboljša natančnost in stabilnost servohidravličnega sistema, novo funkcijo pa je mogoče uporabiti tudi za drugo visokohitrostno, visokonatančno in težko opremo. Omejitev nove funkcije je v funkciji pretvorbe. Zdaj se uporablja sprotni nadzor položaja in tlaka valja ter izračunavanje količnika pretvorbe, ki je počasno in zahteva visokokakovosten nadzorni sistem. Učinkovitost bi se bistveno izboljšala, če bi bilo pretvorbo mogoče zapisati z izrazom. Članek predstavlja nov pristop k združevanju položaja in tlaka, za katera sicer velja, da ju ni mogoče krmiliti istočasno. Za izboljšanje natančnosti in stabilnosti je uvedena funkcija za pretvorbo tlačnega signala, ki se lahko doda položajnemu signalu. Nov krmilni sistem razširja priložnosti za odkrivanje novih načinov krmiljenja servohidravličnih sistemov. Uporaben je v vseh tehničnih sistemih, ki istočasno zahtevajo visoko natančnost položaja in tlaka, zlasti pri opremi z visoko hitrostjo, velikimi pretoki in velikimi bremeni. Ključne besede: potujoče škarje, hidravlični sistem, krmilni sistem, glavno/podrejeno krmiljenje položaja/ tlaka, pretvorba, natančnost, stabilnost SI 48 *Naslov avtorja za dopisovanje: Taiyuan University of Science and Technology, No. 66 Waliu Street, Taiyuan, Kitajska, sunby_sby@126.com Študija obročastega dvofaznega termosifona kot izotermnega vira pri merjenju temperature Florencio Sanchez-Silva* - Ignacio Carvajal-Mariscal - Ariel E. Moreno-Cordobés -Pedro Quinto Diez - Miguel Toledo Velazquez Nacionalni politehnični institut, Šola za strojništvo in elektrotehniko, Mehika Opravljena je bila študija uporabnosti obročastega dvofaznega termosifona (TPT) pri merjenju temperature in glavni rezultati toplotne zmogljivosti so predstavljeni v tem članku. Obročasti TPT je zasnovan kot izotermen vir za umerjanje merilnikov temperature po metodi primerjave. V ta namen je bil zasnovan, zgrajen in eksperimentalno okarakteriziran obročast TPT z vodo kot delovnim fluidom. Za eksperimentalno opredelitev delovanja obročaste površine TPT je bilo aksialno po notranji obročasti površini v območju kondenzatorja nameščenih osem termoparov. Toplotno moč je zagotavljal električni upornik v območju uparjalnika in toploto je zbiral pretočni sistem hladilne vode okrog kondenzatorja. Toplota je bila merjena po načelih kalorimetrije. Preizkusi so bili opravljeni pri dveh različnih pretokih hladilne vode: 16,67 cm3/s in 4,17 cm3/s. Pri vsakem pretoku hladilne vode so bile vzpostavljene štiri različne ravni toplotne moči v območju uparjanja: (38, 154, 350, 1385) W. Rezultati eksperimentov so pokazali, da je delovanje toplotne naprave izotermno pri ravneh moči nad 350 W, kjer se aksialni temperaturni gradient zmanjša na zelo nizko raven in torej ni več vedenja kot pri palici, ki je na enem koncu v stiku s toplotnim virom in oddaja toploto kot rebro. Manjša je tudi napaka meritev temperature. Izotermni pogoji se po zaslugi visoke prevodnosti naprave TPT vzpostavijo že v par minutah in tako je mogoče v kratkem času opraviti več umerjanj. Pri manjših toplotnih močeh (38 W in 154 W) toplotni tok v uparjalnik ni dovolj velik, da bi se aksialni gradient približal vrednosti nič in zato ni doseženo aksialno izotermno vedenje. Ta težava izgine pri toplotnih močeh nad 350 W. Poudariti je treba, da deluje obročasti del TPT pri višjih temperaturah, ko se zmanjša pretok hladilne vode. Drugi cilj raziskave je bil preučitev možnosti uporabe TPT kot sekundarnega vira za umerjanje merilnikov temperature. Za takšno uporabo je nujna stabilnost in ponovljivost TPT pri različnih delovnih pogojih. Vrednosti temperature morajo biti torej enake pri različnih režimih ob enakih delovnih pogojih. Temperaturni profil stene obročastega dela TPT je bil pridobljen v treh ločenih eksperimentalnih preizkusih, opravljenih na različne dni in v enakih delovnih pogojih. Pridobljene vrednosti temperature niso strogo enake - razlika je bila približno 4 K ali 2,9 % in rezultati torej niso ponovljivi. Potopno globino temperaturnih tipal, ki je bila 21 cm, je mogoče povečati s spremembo parametrov geometrije obročastega termosifona. Uporabno temperaturno območje bi bilo mogoče razširiti na poljubno območje s pomočjo primernega delovnega fluida in primerne zmogljivosti ogrevalnega/hladilnega sistema. Na tem področju preostaja še obravnava stabilnosti in ponovljivosti z nadzorom nad vsemi glavnimi parametri eksperimenta. Z vzpostavitvijo takšne ravni nadzora bi bilo obročasti TPT mogoče uporabiti kot sekundarno referenco pri umerjanju temperaturnih instrumentov. Ključne besede: merjenje temperature, dvofazni termosifoni, izotermni, obročasti, umerjanje, eksperiment *Naslov avtorja za dopisovanje: Nacionalni politehnični institut, Šola za strojništvo in elektrotehniko, Campus Zacatenco, Mehika, fsnchz@yahoo.com.mx SI 49 Doktorske disertacije, magistrska dela, diplomske naloge DOKTORSKE DISERTACIJE Na Fakulteti za strojništvo Univerze v Ljubljani sta obranila svojo doktorsko disertacijo: • dne 6. marca 2015 Vanja ČOK z naslovom: »Implementacija Kansei inženiringa v razvojno-konstrukcijski proces« (mentor: prof. dr. Jožef Duhovnik, somentorica: izr. prof. dr. Metoda Dodič Fikfak); Obstoječa metodologija kansei inženiringa tipa 1 upošteva čustvene odzive določene skupine kupcev. Cilj doktorske disertacije je dopolnitev metodologije kansei inženiringa tipa 1, tako da bo metodologijo mogoče uporabiti za analizo dveh ali več skupin kupcev. Razvoj izdelkov za globalni trg zahteva mehanizme, ki omogočajo preverbo odziva lokacijsko različnih kupcev na izdelek. Zanimalo nas je kako z uporabo dopolnjene metodologije kansei inženiringa tipa 1 odkrijemo razlike med dvema skupinama kupcev. S pilotno študijo smo analizirali obstoječi postopek kansei inženiringa tipa 1. V pilotni študiji smo ugotavljali tudi kje se kansei inženiring tipa 1 vključuje v razvojno-konstrukcijski proces. Razvili smo metodološki model KCUL, ki upošteva čustvene odzive več skupin kupcev do vizualnih elementov izdelka v zgodnji fazi koncipiranja, natančneje v nivojih abstrakcije. Model KCUL smo razvili, preverili in potrdili na skupini študentov iz Srednje Evrope in Južne Azije. V doktorskem delu je razvita metoda implementacije nadgrajenega kansei inženiringa tipa 1 v razvojno-konstrukcijski proces; • dne 30. marca 2015 Matej SITAR z naslovom: »Določitev geometrijskega stanja konstrukcijskega elementa iz funkcionalnega gradiva po razbremenitvi in nelinearnim elasto-plastičnem utrjevanju« (mentor: prof. dr. Franc Kosel); Obravnavan je problem določanja geometrijskih in mehanskih stanj enoosnih funkcionalno grajenih nosilcev med in po procesu upogibnega obremenjevanja in razbremenjevanja v elasto-plastičnem območju. Pri tem so nosilci nesimetričnega prečnega prereza izdelani iz gradiv z nelinearno utrjevalnim reološkim modelom. Upoštevana sta tako izotropni kot tudi kinematični model utrjevanja. Stanje specifičnih deformacij je definirano kot funkcija zasuka in premika nevtralne osi ter ukrivljenosti vzdolžne osi nosilca. Upoštevana je celotna zgodovina obremenjevanja, skupaj z efektom lokalnega obremenjevanja določenih vlaken nosilca pri monotono padajoči obremenitvi in obratno. Na osnovi numeričnih primerov je predstavljen vpliv velikosti maksimalne obremenitve na končen in (elastično) povračljiv zasuk nevtralne osi, njen premik ter končno in (elastično) povračljivo ukrivljenost osi nosilca ob različnih oblikah prečnega prereza in gradivih. Ustreznost predstavljenega teoretičnega modela je preverjena z eksperimenti. Ob upoštevanju evolucije elastičnega modula se izkaže, da relativna razlika med numerično in eksperimentalno določenimi rezultati znaša le 0,273 %. * Na Fakulteti za strojništvo Univerze v Mariboru je obranil svojo doktorsko disertacijo: • dne 19. marca 2015 Milan KAMBIČ z naslovom: »Raziskava primernosti ionskih tekočin za uporabo v hidravličnih sistemih« (mentor: izr. prof. dr. Darko Lovrec); Vsaka vrsta hidravlične tekočine ima svoje prednosti za določen primer uporabe. Voda, kot prva in tudi v zadnjih letih ponovno aktualna hidravlična tekočina je negorljiva, praktično nestisljiva, neškodljiva za zdravje človeka ter okolje, cenena in razpoložljiva skoraj povsod. Mineralno olje je možno uporabiti skoraj povsod in v raznovrstni opremi. Težko vnetljive hidravlične tekočine zaradi večje varnosti pred požarom uporabljamo predvsem v metalurških obratih in premogovnikih, biološko hitreje razgradljive hidravlične tekočine pa zaradi manjše nevarnosti za okolje predvsem v mobilni hidravliki. Poleg prednosti pa ima vsaka vrsta hidravlične tekočine tudi slabosti, tako da je pri izbiri najbolj primerne vedno potreben nek kompromis. Idealne rešitve ni, radi pa bi se ji vsaj približali. Cilj doktorske disertacije je bil, da v praktično neomejenem številu različnih ionskih tekočin, skušamo poiskati nekaj primernih za uporabo v hidravličnih sistemih in na ta način predlagati bodočo alternativo današnjim, običajno uporabljanim hidravličnim tekočinam. V ta namen so obravnavane osnovne lastnosti ionskih tekočin in njihove prednosti. Primernost smo ugotavljali tako na osnovi meritev fizikalno kemijskih lastnosti, kot tudi na osnovi praktičnih testov, aktualnih za hidravlične tekočine. V disertaciji so uporabljene laboratorijske in praktične testne metode tudi podrobneje opisane. Veliko pozornosti je posvečeno sposobnosti korozijske zaščite, saj večina testiranih vzorcev ni nudila zadostne zaščite pred korozijo. Prav tako je podrobneje obravnavana stisljivost, saj ima neposreden vpliv na porabo energije in odzivnost hidravličnih sistemov. Pri iskanju ionske tekočine, ki bi bila primerna za uporabo v hidravličnih sistemih se je izkazalo, da kljub nekaterim boljšim rezultatom v primerjavi s konvencionalnimi hidravličnimi tekočinami, (za sedaj) ni veliko takih, ki bi vse potrebne lastnosti imele vsaj na trenutnem nivoju konvencionalnih hidravličnih tekočin. * MAGISTRSKA DELA Na Fakulteti za strojništvo Univerze v Ljubljani so pridobili naziv magister inženir strojništva: dne 2. marca 2015: Samo ČARMAN z naslovom: »Presoja meroslovnih značilnic merilnika padavin« (mentor: izr. prof. dr. Ivan Bajsić); Simon GODEC z naslovom: »Analiza vpliva tehnoloških parametrov visokotlačnega struženja pri obdelavi zlitine Inconel 718« (mentor: doc. dr. Davorin Kramar, somentor: izr. prof. dr. Peter Krajnik); Andrej KASTELIC z naslovom: »Razvoj vsebnika ultra tanke toplotne cevi« (mentor: prof. dr. Iztok Golobič, somentor: doc. dr. Matija Jezeršek); dne 27. marca 2015: Klemen KOTNIK z naslovom: »Modeliranje torzijskih vibracij motorja z notranjim izgorevanjem ter njihov prenos prek pogonskega sistema« (mentor: prof. dr. Miha Boltežar, somentor: prof. dr. Jörg Wallaschek); Marko BRUS z naslovom: »Nadgradnja pravil načrtovanja izdelave oblikovnih plošč za brizganje komponente avtomobilskega žarometa« (mentor: izr. prof. dr. Joško Valentinčič, somentor: doc. dr. Andrej Lebar); Nejc JEGLIČ z naslovom: »Tehnologija izdelave orodja za brizganje plastike in avtomatizacija proizvodnega procesa« (mentor: izr. prof. dr. Joško Valentinčič, somentor: izr. prof. dr. Niko Herakovič). * Na Fakulteti za strojništvo Univerze v Ljubljani je uspešno zagovarjal svoje magistrsko delo magistrskega študijskega programa II. stopnje: dne 2. marca 2015: Pablo GUTIERREZ RUIZ z naslovom: »Povezava med tlačnimi nihanji in kavitacijsko erozijo / Relationship between pressure oscialtions and cavitation erosion« (mentor: izr. prof. dr. Matevž Dular). * Na Fakulteti za strojništvo Univerze v Mariboru sta pridobila naziv magister gospodarski inženir: dne 25. marca 2015: Nejc CVÖRNJEK z naslovom: »Uporaba algoritmov po vzorih iz narave pri napovedovanju cene delnic in optimizaciji portfelja« (mentorja: prof. dr. Miran Brezočnik, prof. dr. Timotej Jagrič); Sandi KAČIČNIK z naslovom: »Dolgoročno planiranje produktnih in storitvenih platform na področju pomivalnih strojev« (mentorja: izr. prof. dr. Iztok Palčič, doc. dr. Igor Vrečko). DIPLOMSKE NALOGE Na Fakulteti za strojništvo Univerze v Ljubljani so pridobili naziv univerzitetni diplomirani inženir strojništva: dne 2. marca 2015: Nejc DRAKSLER z naslovom: »Izdelava vodilnika aksialnega ventilatorja sadjarskega pršilnika« (mentor: izr. prof. dr. Marko Hočevar, somentor: prof. dr. Branko Širok); Gal HLEBANJA z naslovom: »Izdelava in tribološke lastnosti polimernih kompozitov z dodatki grafita in polidimetilsiloksana« (mentor: prof. dr. Mitjan Kalin, somentor: doc. dr. Miroslav Huskić); Alja PUHEK z naslovom: »Vpliv sekundarnega toka na tokovne razmere ob površini krila« (mentor: prof. dr. Branko Širok, somentor: doc. dr. Tom Bajcar); dne 27. marca 2015: Jani BOLE z naslovom: »Modifikacija navoja prednapetega visokotrdnostnega vijaka« (mentor: doc. dr. Boris Jerman, somentor: prof. dr. Marko Nagode); Miha ZUPANC z naslovom: »Analiza mehanskih lastnosti lepljenega spoja med NdFeB magnetom in jeklenim rotorjem« (mentor: doc. dr. Viktor Šajn, somentor: prof. dr. Franc Kosel); Jera AVŠIČ z naslovom: »Primerjava med točkovnim uporovnim in točkovnim laserskim varjenjem« (mentor: prof. dr. Janez Tušek); Jernej OBLAK z naslovom: »Optimiranje nastavitev 3D skenerja za dosego najboljših rezultatov meritev« (mentor: doc. dr. Andrej Lebar, somentor: doc. dr. Davorin Kramar); dne 30. marca 2015: Klemen SPRUK z naslovom: »Uporaba mikro-parne turbine v sistemu SPTE z izkoriščanjem lesne biomase« (mentor: prof. dr. Mihael Sekavčnik). * Na Fakulteti za strojništvo Univerze v Ljubljani je uspešno zagovarjal svoje diplomsko delo univerzitetnega študija : dne 26. marca 2015: Albert PEREZ ESTEVE z naslovom: »Razvoj dogodkovno gnanega simulacijskega modela proizvodne delavnice / Application of a discrete event simulation in a manufacturing workshop« (mentor: prof. dr. Alojzij Sluga). * Na Fakulteti za strojništvo Univerze v Ljubljani sta pridobila naziv diplomirani inženir strojništva: dne 5. marca 2015: Alenka FILIPAŠIĆ z naslovom: »Izbor turbinskega agregata HE Završnica« (mentor: prof. dr. Branko Širok); Uroš KERMC z naslovom: »Toplozračni solarni ogrevalni sistem z latentnim hranilnikom toplote« (mentor: prof. dr. Sašo Medved, somentor: doc. dr. Ciril Arkar). * Na Fakulteti za strojništvo Univerze v Mariboru je pridobil naziv diplomirani inženir strojništva (UN): dne 25. marca 2015: Nejc GRAJŽL z naslovom: »Izdelava hladilnih celic« (mentor: izr. prof. dr. Iztok Palčič). * Na Fakulteti za strojništvo Univerze v Mariboru je pridobil naziv diplomirani gospodarski inženir (UN): dne 25. marca 2015: Samo FIJAVŽ z naslovom: »Priprava projekta selitve podjetja iz najetih v lastne prostore na primeru podjetja MF-CT d.o.o.« (mentorja: izr. prof. dr. Iztok Palčič, doc. dr. Igor Vrečko). * Na Fakulteti za strojništvo Univerze v Mariboru sta pridobila naziv diplomirani inženir strojništva: dne 26. marca 2015: Mitja BLAGOTINŠEK z naslovom: »Primerjava plinskega ogrevanja, ogrevanja s spremenljivim volumnom hladila in njuna kombinacija« (mentor: izr. prof. dr. Jure Marn); Tadej FIGEK z naslovom: »Tehnološke izboljšave pri izdelavi zavorne zagozde za avtomatski menjalnik« (mentor: prof. dr. Miran Brezočnik, somentor: prof. dr. Jože Balič). * Na Fakulteti za strojništvo Univerze v Ljubljani so pridobili naziv diplomirani inženir strojništva (VS): dne 5. marca 2015: Luka KOVAČIČ z naslovom: »Tehnično-ekonomska zasnova male fotonapetostne elektrarne« (mentor: prof. dr. Mihael Sekavčnik); Rok PANDEL z naslovom: »Razvlaknjenje fruktoze na enokolesni centrifugi« (mentor: prof. dr. Branko Širok); Anže SLIVNIK z naslovom: »Analiza aeroprofila krila s spremenljivo obliko« (mentor: izr. prof. dr. Tadej Kosel). * Na Fakulteti za strojništvo Univerze v Mariboru je pridobil naziv diplomirani inženir strojništva (VS): dne 27. marca 2015: Tomaž GRACEJ z naslovom: »Rekonstrukcija dodajalnika okroglih surovcev na 28MN stiskalnici za namene zmanjšanja mrtvega časa med cikli« (mentor: izr. prof. dr. Darko Lovrec, somentor: dr. Vito Tič). Information for Authors All manuscripts must be in English. Pages should be numbered sequentially. The manuscript should be composed in accordance with the Article Template given above. The maximum length of contributions is 10 pages. Longer contributions will only be accepted if authors provide juastification in a cover letter. For full instructions see the Information for Authors section on the journal's website: http://en.sv-jme.eu . SUBMISSION: Submission to SV-JME is made with the implicit understanding that neither the manuscript nor the essence of its content has been published previously either in whole or in part and that it is not being considered for publication elsewhere. All the listed authors should have agreed on the content and the corresponding (submitting) author is responsible for having ensured that this agreement has been reached. The acceptance of an article is based entirely on its scientific merit, as judged by peer review. Scientific articles comprising simulations only will not be accepted for publication; simulations must be accompanied by experimental results carried out to confirm or deny the accuracy of the simulation. Every manuscript submitted to the SV-JME undergoes a peer-review process. The authors are kindly invited to submit the paper through our web site: http://ojs.sv-jme.eu. The Author is able to track the submission through the editorial process - as well as participate in the copyediting and proofreading of submissions accepted for publication - by logging in, and using the username and password provided. SUBMISSION CONTENT: The typical submission material consists of: - A manuscript (A PDF file, with title, all authors with affiliations, abstract, keywords, highlights, inserted figures and tables and references), - Supplementary files: • a manuscript in a WORD file format • a cover letter (please see instructions for composing the cover letter) • a ZIP file containing figures in high resolution in one of the graphical formats (please see instructions for preparing the figure files) • possible appendicies (optional), cover materials, video materials, etc. Incomplete or improperly prepared submissions will be rejected with explanatory comments provided. In this case we will kindly ask the authors to carefully read the Information for Authors and to resubmit their manuscripts taking into consideration our comments. COVER LETTER INSTRUCTIONS: Please add a cover letter stating the following information about the submitted paper: 1. Paper title, list of authors and their affiliations. 2. Type of paper: original scientific paper (1.01), review scientific paper (1.02) or short scientific paper (1.03). 3. A declaration that neither the manuscript nor the essence of its content has been published in whole or in part previously and that it is not being considered for publication elsewhere. 4. State the value of the paper or its practical, theoretical and scientific implications. What is new in the paper with respect to the state-of-the-art in the published papers? Do not repeat the content of your abstract for this purpose. 5. We kindly ask you to suggest at least two reviewers for your paper and give us their names, their full affiliation and contact information, and their scientific research interest. The suggested reviewers should have at least two relevant references (with an impact factor) to the scientific field concerned; they should not be from the same country as the authors and should have no close connection with the authors. FORMAT OF THE MANUSCRIPT: The manuscript should be composed in accordance with the Article Template. The manuscript should be written in the following format: - A Title that adequately describes the content of the manuscript. - A list of Authors and their affiliations. - An Abstract that should not exceed 250 words. The Abstract should state the principal objectives and the scope of the investigation, as well as the methodology employed. It should summarize the results and state the principal conclusions. - 4 to 6 significant key words should follow the abstract to aid indexing. - 4 to 6 highlights; a short collection of bullet points that convey the core findings and provide readers with a quick textual overview of the article. These four to six bullet points should describe the essence of the research (e.g. results or conclusions) and highlight what is distinctive about it. - An Introduction that should provide a review of recent literature and sufficient background information to allow the results of the article to be understood and evaluated. - A Methods section detailing the theoretical or experimental methods used. - An Experimental section that should provide details of the experimental set-up and the methods used to obtain the results. - A Results section that should clearly and concisely present the data, using figures and tables where appropriate. - A Discussion section that should describe the relationships and generalizations shown by the results and discuss the significance of the results, making comparisons with previously published work. (It may be appropriate to combine the Results and Discussion sections into a single section to improve clarity.) - A Conclusions section that should present one or more conclusions drawn from the results and subsequent discussion and should not duplicate the Abstract. - Acknowledgement (optional) of collaboration or preparation assistance may be included. Please note the source of funding for the research. - Nomenclature (optional). Papers with many symbols should have a nomenclature that defines all symbols with units, inserted above the references. If one is used, it must contain all the symbols used in the manuscript and the definitions should not be repeated in the text. In all cases, identify the symbols used if they are not widely recognized in the profession. Define acronyms in the text, not in the nomenclature. - References must be cited consecutively in the text using square brackets [1] and collected together in a reference list at the end of the manuscript. - Appendix(-icies) if any. SPECIAL NOTES Units: The SI system of units for nomenclature, symbols and abbreviations should be followed closely. Symbols for physical quantities in the text should be written in italics (e.g. v, T, n, etc.). Symbols for units that consist of letters should be in plain text (e.g. ms-1, K, min, mm, etc.). Please also see: http://physics.nist.gov/cuu/pdf/sp811.pdf . Abbreviations should be spelt out in full on first appearance followed by the abbreviation in parentheses, e.g. variable time geometry (VTG). The meaning of symbols and units belonging to symbols should be explained in each case or cited in a nomenclature section at the end of the manuscript before the References. Figures (figures, graphs, illustrations digital images, photographs) must be cited in consecutive numerical order in the text and referred to in both the text and the captions as Fig. 1, Fig. 2, etc. Figures should be prepared without borders and on white grounding and should be sent separately in their original formats. If a figure is composed of several parts, please mark each part with a), b), c), etc. and provide an explanation for each part in Figure caption. The caption should be self-explanatory. Letters and numbers should be readable (Arial or Times New Roman, min 6 pt with equal sizes and fonts in all figures). Graphics (submitted as supplementary files) may be exported in resolution good enough for printing (min. 300 dpi) in any common format, e.g. TIFF, BMP or JPG, PDF and should be named Fig1.jpg, Fig2.tif, etc. However, graphs and line drawings should be prepared as vector images, e.g. CDR, AI. Multi-curve graphs should have individual curves marked with a symbol or otherwise provide distinguishing differences using, for example, different thicknesses or dashing. Tables should carry separate titles and must be numbered in consecutive numerical order in the text and referred to in both the text and the captions as Table 1, Table 2, etc. In addition to the physical quantities, such as t (in italics), the units [s] (normal text) should be added in square brackets. Tables should not duplicate data found elsewhere in the manuscript. Tables should be prepared using a table editor and not inserted as a graphic. REFERENCES: A reference list must be included using the following information as a guide. Only cited text references are to be included. Each reference is to be referred to in the text by a number enclosed in a square bracket (i.e. [3] or [2] to [4] for more references; do not combine more than 3 references, explain each). No reference to the author is necessary. References must be numbered and ordered according to where they are first mentioned in the paper, not alphabetically. All references must be complete and accurate. Please add DOI code when available. Examples follow. Journal Papers: Surname 1, Initials, Surname 2, Initials (year). Title. Journal, volume, number, pages, DOI code. [1] Hackenschmidt, R., Alber-Laukant, B., Rieg, F. (2010). Simulating nonlinear materials under centrifugal forces by using intelligent cross-linked simulations. Strojniški vestnik - Journal of Mechanical Engineering, vol. 57, no. 7-8, p. 531-538, D0I:10.5545/sv-jme.2011.013. Journal titles should not be abbreviated. Note that journal title is set in italics. Books: Surname 1, Initials, Surname 2, Initials (year). Title. Publisher, place of publication. [2] Groover, M.P. (2007). Fundamentals of Modern Manufacturing. John Wiley & Sons, Hoboken. Note that the title of the book is italicized. Chapters in Books: Surname 1, Initials, Surname 2, Initials (year). Chapter title. Editor(s) of book, book title. Publisher, place of publication, pages. [3] Carbone, G., Ceccarelli, M. (2005). Legged robotic systems. Kordić, V., Lazinica, A., Merdan, M. (Eds.), Cutting Edge Robotics. Pro literatur Verlag, Mammendorf, p. 553576. Proceedings Papers: Surname 1, Initials, Surname 2, Initials (year). Paper title. Proceedings title, pages. [4] Štefanić, N., Martinčević-Mikić, S., Tošanović, N. (2009). Applied lean system in process industry. MOTSP Conference Proceedings, p. 422-427. Standards: Standard-Code (year). Title. Organisation. Place. [5] ISO/DIS 16000-6.2:2002. Indoor Air — Part 6: Determination of Volatile Organic Compounds in Indoor and Chamber Air by Active Sampling on TENAX TA Sorbent, Thermal Desorption and Gas Chromatography using MSD/FID. International Organization for Standardization. Geneva. WWW pages: Surname, Initials or Company name. Title, from http://address, date of access. [6] Rockwell Automation. Arena, from http://www.arenasimulation.com, accessed on 200909-07. EXTENDED ABSTRACT: When the paper is accepted for publishing, the authors will be requested to send an extended abstract (approx. one A4 page or 3500 to 4000 characters). The instruction for composing the extended abstract are published on-line: http://www.sv-jme.eu/information-for-authors/ . COPYRIGHT: Authors submitting a manuscript do so on the understanding that the work has not been published before, is not being considered for publication elsewhere and has been read and approved by all authors. The submission of the manuscript by the authors means that the authors automatically agree to transfer copyright to SV-JME when the manuscript is accepted for publication. All accepted manuscripts must be accompanied by a Copyright Transfer Agreement, which should be sent to the editor. The work should be original work by the authors and not be published elsewhere in any language without the written consent of the publisher. The proof will be sent to the author showing the final layout of the article. Proof correction must be minimal and executed quickly. Thus it is essential that manuscripts are accurate when submitted. Authors can track the status of their accepted articles on http://en.sv-jme.eu/. PUBLICATION FEE: Authors will be asked to pay a publication fee for each article prior to the article appearing in the journal. However, this fee only needs to be paid after the article has been accepted for publishing. The fee is 240.00 EUR (for articles with maximum of 6 pages), 300.00 EUR (for articles with maximum of 10 pages), plus 30.00 EUR for each additional page. The additional cost for a color page is 90.00 EUR. These fees do not include tax. Strojniški vestnik - Journal of Mechanical Engineering Aškerčeva 6, 1000 Ljubljana, Slovenia, e-mail: info@sv-jme.eu http://www.sv-jme.eu Contents Papers Spela Bolka, Janko Slavič, Miha Boltežar: Identification of Out-of-Plane Material Characteristics through Sheet-Metal Blanking Dario Croccolo, Massimiliano De Agostinis, Stefano Fini, Alfredo Liverani, Nicolò Marinelli, Eugenio Nisini, Giorgio Olmi: Mechanical Characteristics of Two Environmentally Friendly Resins Reinforced with Flax Fibers Yabing Cheng, Yang Wang, Lei Li, Shuaibing Yin, Lichi An, Xiaopeng Wang: Design Method of Dual Phase Hy-Vo Silent Chain Transmission System Franci Čuš, Uroš Župerl: Surface Roughness Control Simulation of Turning Processes Janez Lupše, Leopold Škerget, Jure Ravnik: BEM-Based Algorithm for URANS Simulations of Flow over a Square Cylinder Wang Jun, Sun Binyu, Huang Qingxue, Li Hongzhou, Han Heyong: Research on the Position-Pressure Master-Slave Control for a Rolling Shear Hydraulic Servo System Florencio Sanchez-Silva, Ignacio Carvajal-Mariscal, Ariel E. Moreno-Cordobés, Pedro Quinto Diez, Miguel Toledo Velazquez: Study of an Annular Two-Phase Thermosyphon Used as an Isothermal Source in Thermometry 273 9770039248001