Ratcheting and Rolling Contact Fatigue Crack Initiation Life of Rails under Service Loading. Wenyi YAN Monash University, Australia

Transcription

1 Ratcheting and Rolling Contact Fatigue Crack Initiation Life of Rails under Service Loading Wenyi YAN Monash University, Australia Chung Lun PUN Peter Mutton Qianhua Kan Guozheng Kang

2 Contents Introduction Research objectives A comprehensive approach Experimental study Development of ratcheting model for rail steels Determination of contact pressure and traction distribution Prediction of crack initiation life Summary 2

3 Heavy haul railway in Australia Iron ore and coal High axle load Iron ore (35-40 t) Coal (30 t) passenger trains (18-20 t) Rail Steel High strength, e.g. eutectoid grade Example Rio Tinto in Pilbara Track Services 14 mines 1,400 km A single train 2.4 km long 234 ore cars (31,000 t) 3

4 Rail degradation Rolling contact fatigue Stress level > plastic shakedown limit Kapoor and Johnson (1994) ratcheting Ratcheting strain > ductility limit Std C Std HH HE HT Initiation of rail degradation. i.e. fatigue cracks Wessels L, Oswald S, Welsby D and Mutton P (2015), Managing the transition from wear to rolling contact fatigue in a heavy haul environment, Proc. International Heavy Haul Conference (IHHA 2015), Perth 4

5 Research objectives 1. Evaluate the ratcheting behaviour of three high strength rail steels currently used in heavy haul railways in Australia 2. Develop a ratcheting-based approach to predict crack initiation life of those rail steels 5

6 A Comprehensive approach 1.Experimental program 2.Material model development 4. Ratcheting and crack initiation life 3. Pressure & traction distributions 6

7 Hardness HV (kg/mm 2 ) Hardness HV (kg/mm 2 ) 1. Materials Three high strength pearlitic rail steels Rail C (%) Si (%) Mn (%) P (%) S (%) Cr (%) LAHT HE HE h LAHT HE1 HE LAHT 180 HE1 HE Depth below gauge corner h (mm) Angular displacement (degree) 7

8 Experimental Program Monotonic Tensile Tests Uniaxial Strain Cycling Uniaxial Stress Cycling MTS kN machine Bi-axial Compression- Torsion Stress Cycling 8

9 Axial Stress (MPa) Monotonic Tensile Test Axial Strain (%) LAHT HE1 HE2 D ln L L ln o 1 1 R Elastic Modulus (GPa) Nominal yield strength (MPa) Ultimate tensile strength (MPa) Reduction of area R (%) Ductility D (%) LAHT HE HE

10 Uniaxial symmetrical strain cycling, a = 0.8% Stress amlitude a Mean Stress m LAHT HE1 HE2 (MPa) σ a = 1 2 (σ max σ min ) σ m = 1 2 (σ max + σ min ) N (cycles) All three materials features cyclic softening a with N LAHT & HE2 m HE1 m 10

11 Equivalent shear strain 3-1/2 (%) Equivalent shear strain 3-1/2 (%) Biaxial stress cycling Rectangular loading path Elliptical loading path 0.50 HE2 Loading cycle, N 0.50 HE2 Loading cycle, N Axial Strain (%) Axial Strain (%) 11

12 Biaxial stress cycling 0.0 LAHT 0.0 HE r (%) r (%) Loading Path: -2.0 Linear Oblique Rectangular Butterfly Elliptical N (cycles) Loading Path: -2.8 Linear Oblique Rectangular Butterfly Elliptical N (cycles) HE2 r (%) Loading Path: -2.8 Linear Oblique Rectangular Butterfly Elliptical N ε r N dε r dn Loading paths influence ratcheting Elliptical loading paths lowest axial ratcheting strain N (cycles) ε r = ε max + ε min 2 12

13 2. Ratcheting model for rail steels Abdel-Karim and Ohno model Modification of non-proportional parameter Φ proposed by Tanaka, 1994 by fourth order tensor C. tr( C : C) n : C : C : n 2 tr( C : C) Pun CL, Kan Q, Mutton PJ, Kang G and Yan W (2014). Ratcheting Behaviour of High Strength Rail Steels under Bi-axial Compression-Torsion Loadings: Experiment and Simulation. International Journal of Fatigue 66,

14 Verification of the ratcheting model 0.0 LAHT 0.0 HE r (%) r (%) Exp. Simu. Linear Oblique Rectangular Butterfly Elliptical N (cycles) Exp. Simu Linear Oblique -3.6 Rectangular Butterfly Elliptical N (cycles) HE2 r (%) Exp. Simu Linear Oblique -3.6 Rectangular Butterfly Elliptical N (cycles) Quantification of ratcheting with acceptable accuracy Overprediction / underprediction < 10% 14

15 3. Determination of pressure & traction distributions Non-Hertzian pressure distribution More realistic Obtained from a static 3D FE simulation 15

16 Strip Theory Divide the contact area into thin strips Depends on traction coefficient ξ = P t fl ξ = 0 Free rolling 0 < ξ < 1 Partial slip ξ = 1 Full slip 16

17 Define the stick zone size Relationship between load and size of stick zone, a 0 (Haines and Ollerton (1963) Proc. IMechE) 3 ξ = 1 2 where 2K K K K2 1 K sin 1 2K K 2 K - a 0 a0 a o semi width of the stick zone in x direction a o longest semi width of the contact area in x direction (measured from FEA) Define minor width of stick zone b 0 b 0 = b 0 2 a o a 0 a o a 0 2 where b 0 is the longest semi width of the contact area in y direction (measured from FEA) 17

18 Define stick zone for each strip Stick leading edge Slip trailing edge Define size of stick zone for each strip a i a i a i = a 0 a 0 where a i semi width of the contact area of each strip along x direction (measured directly from FEA) a i semi width of the stick zone of each strip along x direction 18

19 Longitudinal Tangential Traction Carter s theory: For stick zone, τ x x i, y i = fp x i, y i a i fp x a i, y i i For slip zone, τ x x i, y i = fp x i, y i 19

20 4. Numerical Model for Cyclic Rolling contact Fine mesh contact zone 60mm long 30mm wide 18mm depth degrees of freedom Cyclic rolling contact Translate from left to right within the fine mesh contact zone 20

21 Ratcheting strain ε r in rail head Both normal and shear plastic strains accumulate Ratcheting strain p ε r = ε eff max Effective Plastic strain p ε eff = 2 ε 3 ij P P ε ij Ratcheting strain rate dε r dn ε r N ε r N 1 21

22 Ratcheting behaviour dε r dn (a) d r / dn LAHT f = 0.2 f = f = 0.4 f = 0.5 f = 0.6 L = 35 tonnes, = (b) d r / dn 10-2 HE f = 0.2 f = 0.3 f = 0.4 f = 0.5 f = 0.6 L = 35 tonnes, = 0.5 (c) d r / dn HE N (cycles) f = 0.2 f = 0.3 f = 0.4 f = 0.5 f = 0.6 L = 35 tonnes, = N (cycles) N (cycles) Non-zero ratcheting strain rate non-zero net plastic deformation Ratcheting occurs 22

23 Estimation of crack initiation life N i Based on dε r /dn & ductility of the material Stabilized max. dε r /dn of the entire model Criterion for stabilized max. dε r /dn dε r/dn i dε r /dn i 1 row dε r /dn i 1 Ductility ln ( L L 0 ) ln L L 0 = ln 1 1 R Crack initiation life N i N i = ln L L0 dε r /dn < 0.5 % for 5 cycles in a, R = reduction of area Materials ln L L 0 (%) LAHT HE HE

24 N i (10 6 cycles) Crack initiation life influence of LAHT HE1 HE2 L = 35 tonnes, f = 0.4 Practical Crack Initiation Life for HE2 L = 35 tonnes, f = 0.4 with traction ξ N i When 0.5, the LAHT has the longest N i among all three rail steels When > 0.5, the HE2 has the longest N i 24

25 N i (10 6 cycles) Crack initiation life influence of friction LAHT HE1 HE2 L = 35 tonnes, = 0.5 Practical Crack Initiation Life for HE2 L = 35 tonnes, f = 0.4 with traction 2 1 f N i Rapid reduction of N i for LAHT when f > 0.4 f HE1 and HE2 more consistent decreasing rate 25

26 N i (10 6 cycles) Crack initiation life influence of axle load σ y (MPa) σ ys (MPa) σ ys/ σ y 5 4 LAHT HE1 HE2 f = 0.4, = 0.5 Practical Crack Initiation Life for HE2 L = 35 tonnes, f = 0.4 with traction LAHT HE HE L N i L (tonnes) When L 35, the LAHT has the longest N i among all three rail steels When L > 35, the HE2 has the longest N i The HE1 steel has almost constant N i. 26

27 Summary 1.Experimental program 2.Material model development 4. Ratcheting and crack initiation life 3. Pressure & traction distributions 27

28 Acknowledgments Australian Research Council Linkage Project (LP ) Rio Tinto Iron Ore National Computational Infrastructure (NCI), Canberra, Australia Pun CL, Kan Q, Mutton PJ, Kang G and Yan W (2014). Ratcheting Behaviour of High Strength Rail Steels under Bi-axial Compression-Torsion Loadings: Experiment and Simulation. International Journal of Fatigue 66, Pun CL, Kan Q, Mutton PJ, Kang G and Yan W (2015). An Efficient Computational Approach to Evaluate the Ratcheting Performance of Rail Steels under Cyclic Rolling Contact in Service. International Journal of Mechanical Sciences. Available online. 28