TEMPERATURE AND HEAT EFFECTS ON POLYETHYLENE BEHAVIOUR IN THE PRESENCE OF IMPERFECTIONS

Similar documents
MAE 171A MECHANICAL ENGINEERING LABORATORY Materials Testing Laboratory Week 1 - LINEAR ELASTIC FRACTURE MECHANICS

UNIT 1 (Class4) Introduction and Concepts

Infrared Imaging for Material Characterization in Fracture Mechanics Experiments

EXPERIMENTAL AND NUMERICAL INVESTIGATIONS OF CONCRETE BEHAVIOUR AT MESO-LEVEL DURING QUASI-STATIC SPLITTING TENSION

FACTORS INFLUENCING CRB TEST RESULTS

Tensile Testing of Plastics

INGE Engineering Materials. Chapter 7: Part 2. Mechanical Failure. INGE Engineering Materials. Failure Analysis

Infrared Imaging for Material Characterization in Fracture Mechanic Experiments

TI Typical HDPE raw material and pipe test methods EDITION 0607 PAGE 1/10

4 Image Analysis of plastic deformation in the fracture of paper

True Stress and True Strain

Chapter Outline: Failure

THE USE OF FIELD INDENTATION MICROPROBE IN MEASURING MECHANICAL PROPERTIES OF WELDS

ENGR 151: Materials of Engineering LECTURE #12-13: DISLOCATIONS AND STRENGTHENING MECHANISMS

Mechanical Behaviour of Materials Chapter 10 Fracture morpholgy

1) Fracture, ductile and brittle fracture 2) Fracture mechanics

Description of concrete fracture at meso-level using Finite Element Method based on X-ray micro-ct images

Infrared Inspection of Steel Structures Under Load by Barbara Szymanik and Grzegorz Psuj

AN INVESTIGATION OF THE EFFECT OF THE MACHINING PROCESSES OF CIRCULAR HOLES ON THE MECHANICAL BEHAVIOUR OF FLAX FIBER REINFORCED POLYMER

EFFECT OF LOCAL PLASTIC STRETCH OM TOTAL FATIGUE LIFE EVALUATION

FSW TIG Welding of Cu 99 Copper

ME 207 Material Science I

Polymers Performance Overview. Bruce Duncan. Wednesday 18 th April 2007

Characteristics. PE-GFR pressure pipes

Fourier analysis applied to Infrared Thermography of Fiber Composites used for the Strengthening of Structural Elements

Fourier analysis applied to Infrared Thermography of Fiber Composites used for the Strengthening of Structural Elements

MSE200 Lecture 9 (CH ) Mechanical Properties II Instructor: Yuntian Zhu

DETERMINING THE RESIDUAL QUALITY OF PE PIPES USING THE STRAIN HARDENING TEST

Abstract. 1. Introduction. 2. Literature Review

Investigations of fracture process in concrete using X-ray micro-ct

INFLUENCE OF LONG-TERM ANNEALING ON RESIDUAL STRESS DISTRIBUTION AND QUASI-BRITTLE FAILURE PROPERTIES OF TALCUM REINFORCED PIPE GRADE POLYPROYPLENE

PASSIVE INFRARED THERMOGRAPHY FOR DAMAGE MONITORING DURING STRUCTURAL TESTING OF CFRP PARTS

Simulation of the Cutting Process Basic Instability Using Molecular Dynamics Technique

Influence of Ductility Ingredients of Structural Adhesives on Fracture Energy under Static Mixed-Mode Loading 1. Introduction

FATIGUE BEHAVIOUR OF CLINCHED JOINTS

Numerical Modelling and Analysis of Wall loss defects in Thin cylindrical GRP pipes using Flash Thermography

SAFETY ASSESMENT OF PRESSURE VESSELS

Introduction to Engineering Materials ENGR2000 Chapter 8: Failure. Dr. Coates

Thermography- A Tool for Understanding Dynamics of Destructive Testing in Composites

DIRECT MEASUREMENT ON FRACTURE TOUGHNESS OF CARBON FIBER ABSTRACT

Mechanical behavior of crystalline materials- Comprehensive Behaviour

Analysis of Perforated Steel Channel Members When Subjected to Compression M. Jayakumar BhaskarM.Tech 1 P. Sai Pravardhan ReddyM.

PLASTIC PIPE TERMS & DEFINITIONS

Chapter 14 Fracture Mechanics

MECHANICAL PROPERTIES OF MATERIALS

FRICTION STIR PROCESSING OF THE CAST ALUMINUM ALLOY EN AW 5083 (AlMg4,5Mn0,7)

Stress-Strain Behavior

CHAPTER 4 COMPRESSIVE STRENGTH OF HYBRID FIBRE REINFORCED CONCRETE

High Density Polyethylene (HDPE) As a Marine Material. Written by: Tobias Clarke

INVESTIGATIONS OF THE MIXED MODE CRACK GROWTH BEHAVIOR OF AN ALUMINUM ALLOY

Taylor & Francis Taylor & Francis Croup Boca Raton London New York

Recent Advances in NDE Technologies for Turbines and Generators

Engineering. Composiies. CRC Press. Taylor & Francis Croup. Boca Raton London New York. CRC Press is an imprint of the

Pulsed LASER-TIG hybrid welding of coated unalloyed steel thin sheets

The Future of Mix Performance Testing: Fatigue/Thermal Cracking. Presented by Tom Brovolod

Welcome to ENR116 Engineering Materials. This lecture summary is part of module 2, Material Properties.

Study the pattern of J-profile along the crack front through 3D finite element analysis

Tensile Strength and Pseudo-elasticity of YAG Laser Spot Melted Ti-Ni Shape Memory Alloy Wires

FEA simulation and geometric calibration of Arcan fixture for butterfly specimen of RP material

MANUFACTURING TECHNOLOGY

Ductile to Brittle Transition Behavior of Super Duplex Stainless Steels

MECHANICAL BEHAVIOR OF GAS METAL ARC WELDS AND LITHIUM- ION BATTERY MODULE SPECIMENS. Catherine M. Amodeo

COMPARISON BETWEEN DUCTILE TEARING ANALYSIS AND LINEAR ELASTIC FRACTURE MECHANICS ANALYSIS

ACOUSTIC EMISSION DURING STRUCTURE CHANGES IN SEMI-CRYSTALLINE POLYMERS

NDE Inspections of the NRU Calandria Vessel Part 1: NRU and ET Inspections

FAILURE ANALYSIS OF THE TITANIUM ALLOY HORN USED IN ULTRASONIC PROCESSING OF POLYMERIC MATERIALS IN THE AUTOMOTIVE INDUSTRY

MACROSTRUCTURE, MICROSTRUCTURE AND MICROHARDNESS ANALYSIS

Equation 1 Proposed Method for Evaluating Complex Molded Polyethylene Geometries

Introduction to Joining Processes

Module 8: Composite Testing Lecture 36: Quality Assessment and Physical Properties. Introduction. The Lecture Contains

Research in Experimental Solid Mechanics. Krishna Jonnalagadda Mechanical Engineering Indian Institute of Technology Bombay

Structural Vessel Repairs Using Automated Weld Overlays

Impact Tester. Ball Drop Easy operating falling dart instrument according ASTM D 1709 and ISO Instruments for mechanical testing.

INFLUENCE OF MATERIAL VELOCITY ON HEAT GENERATION DURING LINEAR WELDING STAGE OF FRICTION STIR WELDING

Material HDT Continuous Use Temp. Radal (Polyphenylsulfone) 400 F 300 F

How do we find ultimate properties?

Effect of Heat Cycle in Multi-layer Welding on Charpy Absorbed Energy in Low and High Toughness Steel

DNVGL-CP-0347 Edition May 2016

Introduction. 1. Testing procedure. 1.1 Specimen shape

Dynamic Shear Rheometer: DSR

3. Residual Stresses

Development of self-reference lock-in thermography and its application to remote nondestructive inspection of fatigue cracks in steel bridges

Fracture. Brittle vs. Ductile Fracture Ductile materials more plastic deformation and energy absorption (toughness) before fracture.

Characterization of Chip & Cut Behavior with an Instrumented Laboratory Device

Fundamental Course in Mechanical Processing of Materials. Exercises

DELAMINATION CRACK GROWTH OF UNIDIRECTIONAL CFRP IN THERMO-MECHANICAL FATIGUE

Study of Thickness Effect on Fracture Toughness of High Grade Pipeline Steel

FRACTURE CHARACTERIZATION OF BIODEGRADABLE PLLA POLYMER BLENDS

STUDY OF SENT SPECIMENS WITH A TILTED NOTCH TO EVALUATE DUCTILE TEARING IN SPIRAL WELDED PIPELINE APPLICATIONS

HARMONISATION OF POLYETHYLENE PIPE BUTTFUSION PROCEDURES AND TEST METHODS

Flexural Behavior of Casing Joint of Square Steel Tube

CH 6: Fatigue Failure Resulting from Variable Loading

STP1220-EB/Aug Subject Index

Fracture Behaviour of Three-layer PE-HD Pipe

RESEARCH ON WELDING OF PLASTIC MATERIALS

Prediction of Residual Stresses after Local Post-Weld Heat Treatment

Effects of γ-irradiation on the Mechanical and Relaxation Behaviour of High and Low Density Polyethylene.

MACHINES DESIGN SSC-JE STAFF SELECTION COMMISSION MECHANICAL ENGINEERING STUDY MATERIAL MACHINES DESIGN

Types of Fatigue. Crack Initiation and Propagation. Variability in Fatigue Data. Outline

INFLUENCE OF MEMBER COMPONENTS ON THE STRCUTURAL PERFORMANCE OF BEAM-TO-COLUMN JOINTS OF PITCHED ROOF PORTAL FRAMES WITH CLASS 3 AND 4 SECTIONS

Transcription:

THERMAL SCIENCE: Year 2016, Vol. 20, No. 5, pp. 1703-1712 1703 TEMPERATURE AND HEAT EFFECTS ON POLYETHYLENE BEHAVIOUR IN THE PRESENCE OF IMPERFECTIONS by Alin Constantin MURARIU a * and Jasmina V. LOZANOVIĆ ŠAJIĆ b National R&D Institute for Welding and Materials Testing ISIM Timisoara, Timisoara, Romania b Innovation Center of the Faculty of Machanical Engineering, University of Belgrade, Belgrade, Serbia Original scientific paper DOI:10.2298/TSCI151110220M This paper highlights the changes of polyethylene behaviour during various loading rates as well as the influence of test temperature on the material characteristics. Passive infrared thermography method and a high speed infrared camera were used to observe the temperature changes of the sample surface during the tests. The experimental program was carried out on samples taken from PE80 polyethylene gas pipes with simulated imperfections with bilateral V-notch, U-notch, and central hole. Samples have been tensile tested and the results are correlated with the temperature distribution of the samples surface. Key words: tensile tests, polyethylene pipes, infrared thermography, imperfections Introduction High density polyethylene (HDPE) is a product of particular importance, currently used to manufacture pipelines for gas and water supply. Polyethylene has a molecular structure in which monomers are linked together by attractive intermolecular forces, thus resulting polymers. Unlike metals, polymers do not contain a crystallographic plan and grain boundaries; they rather consist of long molecular chains. The factors that govern the toughness and ductility of polymers include the strain rate, temperature and molecular structure. Tests carried out by conventional destructive methods provide information on the structure and mechanical characteristics of the virgin material or estimate them after a certain operation interval in different loading conditions and environments [1-4]. In order to achieve a product, in many cases the material choice is influenced by its long term properties [5-7]. As a result of this trend, in order to study the new materials and for their full characterization, there have been developed new techniques based on the concepts and theories of fracture mechanics [8, 9]. In this paper passive infrared thermography is used to monitor mechanical tensile tests performed under different conditions (loading rate and temperature). This method allows us to evaluate the fracture behaviour of the material, without or in the presence of stress concentrators and localization of failure initiation. * Corresponding author; e-mail: alin@isim.ro

1704 THERMAL SCIENCE: Year 2016, Vol. 20, No. 5, pp. 1703-1712 Materials and methods The experimental program used samples extracted from PE 80 SRD 11 GAS ϕ160 15.5 mm polyethylene pipes. Strip samples with width b = 20.0 mm and thickness a = 15.5 mm were cut from sample pipes with length L = 200 mm. The experimental sets consisted of tree types of test samples: (1) Sample set for evaluation of the influence of test temperature on the mechanical characteristics. The testing conditions were: loading rate: v = 50 mm/min; test temperatures: T = 1 ºC; 5 ºC; 10 ºC; 15 ºC; 20 ºC; 23 ºC; 25 ºC; 30 ºC; 40 ºC; 50 ºC, test environment: water, and number of sets: 10 sets of 5 samples (1 set for each test temperature). (2) Sample set for assessment of the influence of imperfection type and dimension on the failure dynamics and fracture character. The simulated imperfections were obtained by milling or drilling: bilateral V-notched samples with notch depth h = 1.0; 1.5; 2.0; 2.5; 3.0 mm, bilateral U-notched samples with notch depth h = 1.0; 1.5; 2.0; 2.5; 3.0 mm, and sample with central hole with diameter d = 2.0; 3.0; 4.0; 5.0; 6.0 mm. The testing conditions were: loading rate: v = 50 mm/min.; test temperature: T = 23 ± 2 ºC; test environment: air, and number of sets: 15 sets of 5 samples (1 set for type/dimension of simulated imperfection). Notes: In case of V-notched specimens, the notches were obtained using a milling tool of 2 mm width and a radius of 0.25 mm at the notch top. In case of U-notched specimens, a cylindrical milling tool of 2 mm width was used. In case of specimens with central hole, the holes were drilled in mid width. There are many imperfection simulation techniques. For this experimental research the mechanical processing technique has been applied. It should be stated that in order to avoid the difficulties that may arise in processing of these simulated imperfections, optimization of the splintering process is necessary to minimize the sample heating. In the case of accentuated material warming, the splintering becomes cumbersome; the superheated material adheres to the splintering tool resulting in alteration of shape and dimensional accuracy of the simulated imperfections. Tensile tests (TT) carried out in the experimental program were monitored by infrared thermography (IRT) by following the deformation and degradation process of the bilateral notched specimens up to fracture. The test method considered that during the mechanical loading process followed by plastic deformation, the majority of the materials heat up. Following this observation the thermography method was used, being based on the temperature distribution analysis of the specimen surface (analysis of the cold and hot points). The heat energy emission during the sample s loading time period is highlighted, through visualization and quantification by non-contact measuring of the temperature zone of interest (zone of stress raisers created by the presence of plane imperfections). This hybrid TT-IRT method was used for a better understanding of the fracture process in samples with simulated imperfections. Using this method the thermal images of the tested specimens are recorded on-line and analyzed. Thermal image distributions show continuous changes during the test, due to the heat released normally during the static tensile loading. The computerized infrared thermography is used to scan the thermal values of the areas with simulated imperfections and based on the values obtained the severity and progress phase are evaluated.

Murariu, A. C., et al.: Temperature and Heat Effects on Polyethylene Behaviour... THERMAL SCIENCE: Year 2016, Vol. 20, No. 5, pp. 1703-1712 1705 The experimental set-up (fig. 1) consisted of: universal testing equipment ZD 10/90 type, air conditioned room, infrared camera FLIR System A40, and a personal computer with appropriate software for data acquisition and processing: field of view minimal focal distance: 24 x 18 / 0.3 m, instantaneous field of view (IFOV): 1.3 mrad, thermal sensitivity at 50/60 Hz: 0.08-30 C, focus: automatic type of infrared detector: Focal Plane Array (FPA), atmospheric window: 7.5-13 µm, temperature ranges: 40-+55 C; 0-120 C; 0-1500 C, automatic emissivity correction: variable form 0.1-1.0, and port (IEEE-1394): 16 bit images. (a) (b) Figure 1. Experimental set-up; (a) ZD10/90 Universal Testing Machine and IR camera FLIR System Type A40, (b) EDZ 40 Universal Testing Machine and IR camera FLIR System Type A40 Experimental results The experimental results obtained by tensile testing, performed on sample sets designed to evaluate the influence of test temperature on the tensile strength, are presented in fig. 2. Influence of the test temperature on the tensile strength Rm [N/mm2] 30 28 26 y = -0.3056x + 27.067 R2 = 0,9961 24 22 20 18 16 14 12 10 0 5 10 15 20 25 30 35 40 45 50 Test temperature, T [ C] Testthe temperature T [ºC] Figure 2. Variation of tensile strength for PE 80 with test temperature

1706 THERMAL SCIENCE: Year 2016, Vol. 20, No. 5, pp. 1703-1712 Evolution of the samples temperature during the tensile testing, performed on strip samples extracted from the base material of the pipe, is presented in fig. 3. Figure 4 presents IR image of the tensile testing for loading rates 25, 75, and 150 mm/min. (a) Loading rate: v =1 mm/min. (b) Loading rate: v =5 mm/min. (c) Loading rate: v =10 mm/min. Evolution of the base material sample s temperature during the tensil (a (d) Loading rate: v =25 mm/min (e) Loading rate: v =50 mm/min. (f) Loading rate: v =100 mm/min. (g) Loading rate: v =150 mm/min Evolution of the base material sample s temperature during the tensile tests Figure 3. Evolution of the base material sample s temperature during the tensile tests

THERMAL SCIENCE: Year 2016, Vol. 20, No. 5, pp. 1703-1712 1707 Evolution of the samples temperature during the tensile testing, performed on the bilateral V-notched sample set, is presented in fig. 5. Figure 6 presents IR image of the tensile testing for notch depths 1, 2, and 3 mm. Evolution of the samples temperature during the tensile testing, performed on the bilateral U-notched sample set, is presented in fig. 7. Figure 8 presents IR image of the tensile testing for notch depths 1, 2 and 3 mm. (a) (b) (c) Figure 4. The IR image for different loading rates; (a) 25 mm/min, (b) 50 mm/min, (c) 150 mm/min 100.0 C 23.0 C (a) notch depth h = 1.0 ) (b) notch depth h = 1.5 (c) notch depth h = 2.0 (d) notch depth h = 2.5 (e) notch depth h = 3.0 notched sample s temperature during the tensile tests Figure 5. Evolution of the V-notched sample s temperature during the tensile tests

1708 THERMAL SCIENCE: Year 2016, Vol. 20, No. 5, pp. 1703-1712 100.0 C Figure 6. The IR image for different notch depths; (a) 1 mm, (b) 2 mm, (c) 3 mm (а) (b) (c) 23.0 C (a) notch depth h = 1.0 (b) notch depth h = 1.5 (c) notch depth h = 2.0 (d) notch depth h = 2.5 (e) notch depth h = 3.0 Figure 7. Evolution of the U-notched sample s temperature during the tensile tests Evolution of the samples temperature during the tensile testing, performed on the sample set with the central hole, is presented in fig. 9, whereas IR images are given in fig. 10.

THERMAL SCIENCE: Year 2016, Vol. 20, No. 5, pp. 1703-1712 1709 100.0 C Figure 8. The IR image for different notch depths; (a) 1 mm, (b) 2 mm, (c) 3 mm (a) (b) (c) 23.0 C (a) hole diameter d = 2.0 (b) hole diameter d = 3.0 (c) hole diameter d = 4.0 (d) hole diameter d = 5.0 (e) hole diameter d = 6.0 Figure 9. Evolution of the sample s temperature with the central hole, during the tensile tests Discussion As it can be seen in fig. 2, the test temperature plays a decisive role, since the tensile strength of PE80 thermoplastic material is linearly dependent with it, in the temperature range in which experimental tests were conducted (range included in the field of viscoelastic behaviour).

1710 THERMAL SCIENCE: Year 2016, Vol. 20, No. 5, pp. 1703-1712 100.0 C Figure 10. The IR image for different hole diameters; (a) 2 mm, (b) 4 mm, (c) 6 mm (a) (b) (c) 23.0 C The tensile strength determined at a temperature of 10 C is about two times higher than the one determined at a temperature of 50 C, while in the temperature range of 23 ± 2 ºC it varies by approximately ±5%, thus affecting the measurement uncertainty. This demonstrates both the need for specimen conditioning before the test and a rigorous control of the test temperature. At higher loading rate the material does not have enough time to relax and so forced displacement of molecules appears that produces friction and a higher load is required for material deformation. Thus, by friction some of the energy is converted into heat which dissipates into the material. This aspect could be evidenced by the maximum temperature obtained during the test, as seen in fig. 11. This temperature is much lower at a testing rate of 1 mm/min. (44.8 C) compared to the one measured at a testing rate of 50 mm/min. (81.7 C) and the one measured at a testing rate of 150 mm/min. (93 C). Complex behaviour in the range 50-150 mm/min might be attributed to different, competing mechanisms of internal heating. Temperature, T [ C] 100,0 90,0 80,0 70,0 Influence of the loading rate on the maximal temperature of the sample surface y = 7E-05x 3-0.0176x 2 + 1.4722x + 42.989 R 2 = 0.9959 60,0 50,0 40,0 0 25 50 75 100 125 150 Loading rate, v [mm min 1 ] Loading rate v [mm/min] Figure 11. Evolution of the maximum temperature vs. loading rate recorded during the tensile test performed on samples extracted from the base material, without imperfections To highlight how the imperfections effects on the tensile testing behaviour of polyethylene, by processing the data given within figs. 5, 7, and 9, fig. 12 presents the evolution of the maximum recorded temperature vs. the percentage of cross section occupied by the imperfections during the tensile test, performed on samples extracted from the base material with

THERMAL SCIENCE: Year 2016, Vol. 20, No. 5, pp. 1703-1712 1711 different type of simulated imperfections. Once again, complex behaviour might be attributed to competing mechanisms of internal heating. Temperature, T [ C] 100,0 90,0 80,0 70,0 Influence of the imperfection type and dimensions on the maximal temperature of the sample surface V-notches U-notches Central Hole y = 0.0072x -0.0072x 3 + 0.4631x 2-8.6957x + 110.84 R 2 = 0.9929 y = 0.003x 3-0.1674x 2 + 2.3421x + 63.72 R 2 = 1 y = 0.0149x 3 + 0.8654x 2 15.414x + 142.62 R 2 = 0.9504 60,0 50,0 40,0 5 10 15 20 25 30 35 The percentage of cross-section occupied by the imperfections [%] Figure 12. Evolution of the maximum temperature vs. the percentage of cross-section occupied by the imperfections during the tensile test performed on samples with different types of imperfections Conclusions It was found that recorded temperatures during the tensile tests are correlated with the local stress condition, induced by loading of the specimen with imperfections (geometrical discontinuities obtained by mechanical machining). The maximum temperature recorded during the tensile test varies with the loading rate, the curve shape being altered by the imperfections depending on their type and the percentage of the specimen s cross-section occupied by them. Further investigation is needed to understand mechanisms of complex temperature behaviour in presence of imperfections and with different loading rates. Acknowledgement The author greatly appreciates the support of the National Research and Development Institute for Welding and Materials Testing - ISIM Timisoara, Romania. The author would like to thank Professor Emeritus V. I. Safta, member of Romanian Academy of Technical Sciences, for co-ordination of the experimental research. References [1] Kratochvilla, T. R. U., et al., Experiences with Modified Test Condition for Notch Pipe Testing. Polym. Test., 27 (2008), 2, pp. 158-160 [2] Prabhakaran, R., et al., Notch Sensitivity of Polymers, J. Appl. Polym. Sci., 22 (2003), 10, pp. 3011-3020 [3] Peres, F. M., et al., Essential Work of Fracture Testing Method Applied to Medium Density Polyethylene, Procedia Materials Science, 3 (2014), July, pp. 756-763 [4] Linul, E., et al., Determination of Mixed Mode Fracture Toughness of Pur Foams, Structural Integrity and Life, 14 (2014), 2, pp. 87-92 [5] Pinter, G., et al., Cyclic Crack Growth Tests with CRB Specimens for the Evaluation of the Long-Term Performance of PE Pipe Grades, J. Polym. Test. 26 (2007), 2, pp. 180-188

1712 THERMAL SCIENCE: Year 2016, Vol. 20, No. 5, pp. 1703-1712 [6] Frank, A., et al., Prediction of the Remaining Lifetime of Polyethylene Pipes After up to 30 Years in Use, J. Polymer Testing, 28 (2009), 7, pp. 737-745 [7] Murariu, A. C., et al., Long-Term Behaviour of Polyethylene PE 80 Pressurized Pipes, in Presence of Longitudinal Simulated Imperfections, Materiale Plastice, 47 (2010), Mar., pp. 263-266 [8] Anderson, T. L., Fracture Mechanics. Fundamentals and Applications, 3 rd ed., CRC Press, Boca Raton, Fla., USA, 2005 [9] Uchida, M., Tada, N., Micro-, Meso- to Macroscopic Modeling of Deformation Behavior of Semicrystalline Polymer, International Journal of Plasticity, 49 (2013), Oct., pp. 164-184 Paper submitted: November 10, 2015 Paper revised: December 14, 2015 Paper accepted: December 23, 2015