Microsoft Word - numero 25 art 7 J. Tong et alii, Frattura ed Integrità Strutturale, 25 (2013) 44-49; DOI: 10.3221/IGF-ESIS.25.07 44 Special Issue: Characterization of Crack Tip Stress Field Near tip strain evolution under cyclic loading J. Tong, Y.-W. Lu, B. Lin University of Portsmouth, UK Y. H. Tai Rolls-Royce plc, UK J.R. Yates University of Manchester, UK ABSTRACT. The concept of ratchetting strain as a crack driving force in controlling crack growth has previously been explored at Portsmouth using numerical approaches for nickel-based superalloys. In this paper, we report the first experimental observations of the near-tip strain evolution as captured by the Digital Image Correlation (DIC) technique on a compact tension specimen of stainless steel 316L. The evolution of the near-tip strains with loading cycles was studied whilst the crack tip was maintained stationary. The strains were monitored over the selected distances from the crack tip for a given number of cycles under an incremental loading regime. The results show that strain ratchetting does occur with load cycling, and is particularly evident close to the crack tip and under higher loads. A finite element model has been developed to simulate the experiments and the simulation results are compared with the DIC measurements. KEYWORDS. DIC; FE; Ratchetting; Crack tip mechanics; Fatigue crack growth. INTRODUCTION he concept of ratchetting strain as a crack driving force in controlling crack growth has been investigated utilising a variety of constitutive models, including elastic-plastic, visco-plastic and crystal-plastic formulations, for nickel- based superalloys [1-3]. Crack tip deformation fields were examined for both stationary and growing cracks at room and 650C using the finite element method. Distinctive strain ratchetting behaviour near the crack tip was identified in all cases, leading to progressive accumulation of tensile strains normal to the crack growth plane. It was hypothesised that this tensile strain may be responsible for material separation leading to crack growth. Most recently the concept has been applied successfully [4] to fatigue crack growth of a nickel alloy in vacuum at a range of temperatures, where the influence of oxidation is removed. Although this latest work appears to be very encouraging, no direct experimental evidence is yet available to support this line of reasoning. In this paper, we present the first series of experimental results on the measurement of the near-tip strains as a function of the number of loading cycles in a compact tension specimen of stainless steel 316L using the DIC method. Finite element analysis has also been carried out on the specimen and the near-tip strain results from the simulation are compared with those obtained experimentally. T http://www.gruppofrattura.it http://dx.medra.org/10.3221/IGF-ESIS.25.07&auth=true J. Tong et alii, Frattura ed Integrità Strutturale, 25 (2013) 44-49; DOI: 10.3221/IGF-ESIS.25.07 45 EXPERIMENTAL METHODS he material used for the investigation was stainless steel 316L. A compact tension (CT) specimen was employed with a width of 56 mm, a thickness of 7 mm and a machined notch size of 24 mm. An average grain size of the material along the rolling direction was found to be approximately 28.5 µm. The specimen was pre-cracked first in a servo-hydraulic test machine under cyclic tension, allowing a crack growth of ~2.05 mm to obtain a final crack length of ~26.05 mm. The cyclic testing was then conducted at constant load amplitude with the range of load increased after a given number (20) of cycles, as shown in Fig. 1. The evolution of the strain fields near the crack tip was monitored during the cyclic loading using the digital image correlation (DIC) technique. The DIC system employed is a stereomicroscope system, Vic-3D Micro™ by Correlated Solutions [5]. Speckle patterns were painted on the specimen surface to facilitate the image analysis, and a resolution of 1224 x 1024 pixels was achieved through the use of the cameras coupled with the microscope. Four series of tests were conducted where the near-tip strains were monitored and recorded using the DIC method. A loading frequency of 0.1 Hz was used for the testing to allow sufficient time to collect the data. The images were captured at a framing rate of 5 per second for each cycle and about 50 images were recorded for each cycle. The field of view of the images was approximately 2.15 x 1.80 mm, giving a pixel size of 1.76 µm. Figure 1: The cyclic loading scheme used for the experiment. Digital image correlation has been used extensively to determine fracture parameters using the displacement data extracted from digital images [6, 7]. The basic principles of DIC are to take a set of sequential images for a deformed object, with the first image taken before deformation serves as a reference and the subsequent images acquired at different deformation stages from the same region correlated with the reference image. The algorithm is based on the mathematical correlation of the intensity changes of the sequentially recorded digital images, and implemented through a procedure of finding the best correlation between the two images. Digital images are usually divided into smaller interrogation windows, or sub-sets, a matching process is performed on each of these sub-sets. A full field map of displacements/strains of each subset may be obtained when the correlation process is completed successfully. In the current work, LaVision DaVis 8.1.1 was employed to carry out the image correlation [8]. The size of the subsets was chosen as 28 x 28 pixels with a step size of 3 pixels, sufficient to provide a high spatial resolution and good image correlation quality as well as acceptable computational cost. To study the near-tip strain evolution with load cycles, regions of interest near the crack tip were selected and strain distributions studied. Specifically, two points ahead of the crack and on the crack plane, R2 and R4, were monitored, where the distance to the crack tip, R2 = 28.5 µm and R4 = 57 µm, as illustrated in Fig. 2. These values were selected to be either the same as the average grain size of the material (R2) or about double the average grain size (R4). The strain data at each location were obtained from an average value of multiple points within a 25 µm x 25 µm square. During the loading cycles, the crack length was also monitored and corrected post testing when micro crack growth was detected, such that the values of R2 and R4 stay the same throughout the tests. Micro-crack growth was indeed detected during some of the loading sequences (Test series 3-1 to 3-3, omitted in Fig. 1). These were excluded in the subsequent ratchetting analyses. 0 1 2 3 4 5 6 60 80 100 120 140 160 L oa d [ K N ] Number of Cylces Test 3-4 Test 3-5 Test 3-6 Test 3-7 T http://www.gruppofrattura.it http://dx.medra.org/10.3221/IGF-ESIS.25.07&auth=true J. Tong et alii, Frattura ed Integrità Strutturale, 25 (2013) 44-49; DOI: 10.3221/IGF-ESIS.25.07 46 Figure 2: The random speckle pattern on the specimen and a typical displacement map extracted from the DIC measurement. The locations of R2 and R4 are illustrated with respect to the crack tip. FINITE ELEMENT ANALYSES inite element analyses were carried out on the CT specimen using ABAQUS [9] under plane stress loading conditions to obtain the near-tip strain distribution on the specimen surface, and the results were compared with those measured by the DIC. Material model The material model by Lemaitre and Chaboche [10] was adopted where both isotropic and nonlinear kinematic hardening rules were used to describe the monotonic and the cyclic behaviour of SS316L [11, 12]. The equivalent von Mises stress is defined as:      3 : 2 dev devf S S       (1) where dev is the deviatoric part of the back stress and S is the deviatoric stress tensor. The size of the yield surface is defined using a simple exponential law:  0 0 1 plbQ e      (2) where 0 is the size of the yield surface at zero plastic strain, Q and b are isotropic hardening material parameters. The equivalent plastic strain is given by:      2 2 21 2 2 3 3 1 1 1 1 2 pl v                (3) The evolution of the kinematic hardening component is defined as:   0 1 1pl plC C C                (4) where  is the back stress, and are material parameters, and is the rate of change of C with respect to temperature and field variables. A total of five materials parameters are required to run the FE analyses: Initial yield stress 0; kinematic hardening parameters C and  and isotropic hardening parameters Q∞ and b. The values of these parameters for SS316L were obtained using the experimental data for a strain range 4.0%  [13], and are summarised in Tab. 1. F http://www.gruppofrattura.it http://dx.medra.org/10.3221/IGF-ESIS.25.07&auth=true J. Tong et alii, Frattura ed Integrità Strutturale, 25 (2013) 44-49; DOI: 10.3221/IGF-ESIS.25.07 47 E [GPa] 0 [MPa] Kinematic hardening Isotropic hardening C [MPa]  Q∞ [MPa] b 193 100 60000 280 200 6 Table 1: The material parameters for SS316L [13]. The finite element model The compact tension specimen as used in the experiment was modelled using the finite element method (ABAQUS). Due to symmetry of the geometry and the loading, only half of the specimen was meshed with 4-noded quadrilateral plane stress elements, as shown in Fig. 3. Finer elements were generated around the crack-tip area ( ≈ 1.8 µm) and a rigid line was attached along the symmetry line to prevent the potential penetration of the crack flanks due to crack closure under cyclic loading. A mesh convergence study was carried out, and a mesh size of 3.6 µm was found to be adequate such that the effect of mesh size on the stress-strain responses became negligible. The loading pin (marked in red in Fig. 3) was modelled as elastic and no slip was allowed between the loading pin and the specimen [14]. Cyclic tensile loading was applied at the centre of the pin according to the loading scheme (Fig. 1). A sequential crack tip node-release technique was applied to model the micro-crack growth, where nodes were released incrementally and set to the micro-crack growth length detected, otherwise the crack was assumed stationery. The strain values were calculated at the integration points and the average strain values were obtained for R2 and R4 over a square of 25 µm x 25 µm, similar to those obtained by the DIC method. Figure 3: The FE mesh for the half-model of a CT-specimen, where a finer mesh ( ≈ 1.8 µm) was used in the crack tip region. RESULTS AND DISCUSSION he strain evolutions with cycles are shown in Fig. 4 at R2 and R4 for the four test series (as shown in Fig. 1) from both the DIC and the FE analyses. It is evident from both analyses that tensile strains increase with the number of cycles in all load cases, although the FE analyses appear to predict higher responses than those from the DIC. A closer agreement between the two is achieved at R2, as opposed to R4, indicating that grain size might serve as a suitable candidate as a “critical distance”. Micro-crack growths were detected at some of the early test series, such that ratchetting behaviour could not be monitored on these occasions. These may be due to the pre-cracking load history or the notch effect. Further work is being carried out on controlled crack growths so that the relation between the ratchetting strain and the micro-crack growth may be explored. T Rigid line http://www.gruppofrattura.it http://dx.medra.org/10.3221/IGF-ESIS.25.07&auth=true J. Tong et alii, Frattura ed Integrità Strutturale, 25 (2013) 44-49; DOI: 10.3221/IGF-ESIS.25.07 48 (a) (b) Figure 4: Strain responses at (a) R2 and (b) R4 with the number of cycles for the four test series (Fig. 1), comparison of the results from the DIC and the FE analysis. CONCLUSIONS e report the first experimental evidence of the near-tip strain evolution as captured by the Digital Image Correlation on a compact tension specimen of stainless steel 316L. The evolution of the near-tip strains with loading cycles was monitored at selected load levels whilst the crack tip remains stationery. The results clearly show that strain ratchetting occurs with cycle, and is particularly evident close to the crack tip and under higher loads. Finite element analyses have also been carried out, and the results compare favourably with those measured at a distance about the grain size. REFERENCES [1] Tong, J., Zhao, L.G., Lin, B., Ratchetting strain as a driving force for fatigue crack growth, Int. J. Fatigue, 46 (2013) 49–57. [2] Zhao, L.G., Tong, J., A viscoplastic study of crack-tip deformation and crack growth in a nickel-based superalloy at elevated temperature, J. Mech. Physics Solids, 56 (2008) 3363-3378. [3] Zhao, L.G., Tong, J., Byrne, J., The evolution of the crack tip stress-strain fields and plasticity induced crack closure revisited, Fatigue Fracture Engng. Mater. Struc., 27 (2003) 19-29. [4] Tong, J., Cornet, C., Lin, B., Zhao, L.G., Unpublished results, (2012). [5] http://www.correlatedsolutions.com/index.php/products/vic-3d-2012/microscopy, (2012) [6] McNeill, S.R., Peter, W.H., Sutton, M.A., Estimation of stress intensity factor by digital image correlation, Engineering Fracture Mechanics, 28(1) (1987) 101–112. [7] Durig, B., Zhang, F., McNeill, S.R., Chao, Y.J., Peters III, W.H., A study of mixed mode fracture by photoelasticity and digital image analysis, Optics and Lasers in Engineering, 14(3) (1991) 203–215. [8] LaVision GmbH, StrainMaster Mannual DaVis 8.1, Gottingen, Germany, (2012). [9] ABAQUS 6.8, Hibbitt Karlsson and Sorensen Inc, Providence, RI, (2009). [10] Lemaitre, J., Chaboche, J.L., Mechanics of Solid Materials, Cambridge University Press, (1990). W http://www.gruppofrattura.it http://dx.medra.org/10.3221/IGF-ESIS.25.07&auth=true J. Tong et alii, Frattura ed Integrità Strutturale, 25 (2013) 44-49; DOI: 10.3221/IGF-ESIS.25.07 49 [11] Chaboche, J.L., Rousselier, G., On the plastic and viscoplastic constitutive equations–Part II: Application of internal variable concepts to the 316 stainless steel, Transactions of the ASME, Journal of Pressure Vessel Technology, 105 (1983) 159–64. [12] Chaboche, J.L., Constitutive equations for cyclic plasticity and cyclic viscoplasticity, International Journal of plasticity, 5 (1989) 247–302. [13] Van Eeten, P., Nilsson, F., Constant and variable amplitude cyclic plasticity in 316L stainless steel, Journal of Testing and Evaluation, 34 (2006) 298–311. [14] Kuntiyawichai, K., Burdekin, F.M., Engineering assessment of cracked structures subjected to dynamic loads using fracture mechanics assessment, Engineering Fracture Mechanics, 70 (2003) 1991–2014. http://www.gruppofrattura.it http://dx.medra.org/10.3221/IGF-ESIS.25.07&auth=true