Microsoft Word - numero_30_art_22 E. Sgambiterra et alii, Frattura ed Integrità Strutturale, 30 (2014) 167-173; DOI: 10.3221/IGF-ESIS.30.22 167 Focussed on: Fracture and Structural Integrity related Issues Stress induced martensite at the crack tip in NiTi alloys during fatigue loading E. Sgambitterra, L. Bruno, C. Maletta Dept. of Mechanical, Energy and Management Engineering, University of Calabria, 87036 Rende (CS), Italy. esgambitterra@unical.it ABSTRACT. Crack tip stress-induced phase transformation mechanisms in nickel-titanium alloys (NiTi) were analyzed by Digital Image Correlation (DIC), under fatigue loads. In particular, Single Edge Crack (SEC) specimens, obtained from a commercial pseudoelastic NiTi sheet, and an ad-hoc experimental setup were used, for direct measurements of the near crack tip displacement field by the DIC technique. Furthermore, a fitting procedure was developed to calculate the mode I Stress Intensity Factor (SIF), starting from the measured displacement field. Finally, cyclic tensile tests were performed at different operating temperature, in the range 298-338 K, and the evolution of the SIF was studied, which revealed a marked temperature dependence. KEYWORDS. Shape Memory alloys; Fracture; Digital Image Correlation; Stress Intensity Factor. INTRODUCTION n last decades, shape memory alloys (SMAs), and in particular the Nickel Titanium-based ones (NiTi), received big attention from scientific and engineering communities thanks to their unique characteristics, namely shape memory effect (SME) and superelastic effect (SE) [1]. In particular, these properties allow large recoverable strains or large induced internal forces due to a reversible solid-state phase transformation between austenite and martensite; this transformation can be activated by a temperature variation (TIM, thermally induced martensite) or by the application of external forces (SIM, stress-induced martensite). Due to these interesting features, as well as to their good mechanical performances and biocompatibility, NiTi alloys have seen growing use in many branches of engineering and medicine [2- 3]. As a direct consequence of this interest, many studies were carried out to better investigate on the thermo-mechanical behavior of SMAs, in terms of both SME and SE [1]. Their mechanical properties were investigated using tensile (e.g. [4–5]) and fatigue testing (e.g. [6–11]), however, many aspects are still unknown, especially because the hysteretic stress and/or thermally induced phase transformations significantly affect the damage mechanisms occurring under fatigue loadings, i.e. the crack formation and propagation mechanisms. From a materials science point of view, an accurate knowledge of these topics is essential for predicting the functional and structural life of damaged structures as well as their failure modes, with the aim to improve the overall performances of NiTi-based components or structures. The fracture behavior of austenitic NiTi alloys strongly depends on the SIM transformation which occurs in the crack tip region, as a consequence of the high values of local stresses. As a direct consequence of the marked non-linear and hysteretic behavior, classic elastic and/or elastic-plastic theories cannot be directly applied to SMAs. The constitutive laws implemented in commercial finite element software are based on phenomenological approaches and on theory of plasticity, i.e. they use plasticity-like concepts to describe the effects of phase transformation mechanisms on I E. Sgambiterra et alii, Frattura ed Integrità Strutturale, 30 (2014) 167-173; DOI: 10.3221/IGF-ESIS.30.22 168 the macroscopic response of NiTi alloys [12]. Even though these numerical methods represent useful design tools to simulate the macroscopic response of simple or complex SMA based systems, special care should be taken when they are used to study the local effects in the proximity of high stress concentration regions. In particular, it was demonstrated that the high values of local stresses arising in the crack tip region of NiTi alloys cause stress-induced phase transition mechanisms [13-30], which significantly affect the crack tip stress distribution and, consequently, the crack evolution under both static and fatigue loading conditions. Despite the increasing number of research activities on fracture and fatigue of NiTi alloys in recent years, much effort should be devoted for an effective understanding of the role of the phase transformations in the crack formation and propagation mechanisms and in the stress state generated at the crack tip. Within this context, the development and application of full field techniques to analyze the local transformation mechanisms near geometrical discontinuities and, in particular, in the crack tip region, represent a highly challenging scientific goal. For this purpose, synchrotron X-ray micro-diffraction (XRD) [13-15], infrared thermographic (IR) [16] and Digital Image Correlation (DIC) [18, 19] techniques were recently applied, to better understand the mechanisms of phase transformation in the notch and/or crack tip proximity. In particular, a pseudoelastic NiTi alloy for medical applications were analyzed in [13] by using miniature compact-tension (CT) specimens, which were directly obtained from thin-walled tubes, similar to those used for manufacturing self-expanding stents. XRD microdiffraction investigations of fatigue pre- cracked specimens revealed that the crack tip local strain are due to both B2 to B19’ transformation and to the subsequent loading of the martensitic phase. Strain and texture evolution near the crack tip of a martensitic NiTi alloy were analyzed in [14], by synchrotron X-ray experiments, after fatigue crack propagation in a compact-tension (CT) specimen; it was found that texture evolution is mainly due to detwinning, the main deformation mechanism in martensitic NiTi alloys. Both martensitic and austenitic alloys were analyzed in [15], by using miniaturized CT specimens [19] after fatigue crack propagation, which revealed the presence of detwinned martensite at the crack tip of both martensitic and austenitic specimens. An austenitic NiTi alloy was analyzed in [17] by DIC analysis of thin edge cracked specimen, which allowed direct measurement of the crack tip strain field related to stress-induced transformation. These experimental investigations, together with other several recent numerical [20-22] and analytical [23-30] studies, provide very useful information about the occurrence of crack tip transition mechanisms in NiTi alloys subjected to static and/or monotonic loads, and the effect of the operating temperature on the stress intensity factor (SIF) and on the size of the transformation region was also numerically investigated in [30]. However, the role of these mechanisms on the fracture properties on NiTi alloys is not yet completely defined, i.e. no reliable methodologies were developed for an effective design against fracture. In addition, the evolution of microstructural changes at the crack tip occurring during fatigue loading, i.e. the thermal and mechanical hysteretic behavior, was not yet investigated. This topic is of major concern as the hysteretic nature of phase transitions is expected to play a significant role on the crack growth mechanisms. In this work the stress intensity factor evolution in a NiTi pseudoelastic alloy was investigated by means of a full field experimental techniques In particular, DIC method was used to analyze the displacement field in the crack tip region by which, with a proper fitting procedure based on the William’s series expansion [31], the mode I stress-intensity factor (SIF) was evaluated. Different cyclic tests, a different operating temperatures, in the range 298-338 K, were performed and the effect of the temperature on the SIF was investigated. MATERIAL AND METHODS commercial nickel-titanium sheet (thickness t=0.5 mm) with pseudoelastic properties at room temperature (Austenite finish temperature Af=286.7 K) was used in this investigation. Fig. 1a shows the isothermal (T=298 K) stress strain () response of the material obtained from a complete loading-unloading cycle up to a maximum deformation of about 6.2 %, corresponding to a complete stress-induced martensite transformation. The figure also reports the values of the main thermo-mechanical parameters of the alloy: transformation stresses (sAM, fAM, sMA and fMA), tranformation strain (L), Young’s moduli of austenite and martensite (EA and EM) and Clausius-Clapeyron constants (CA and CM). Edge-notched tension specimens were cut from a commercial plate by electrical discharge machining (EDM) and a 0.1 mm EDM wire was used to machine the notch. The specimen surface has been properly treated in order to provide a suitable speckle pattern for using DIC, therefore, the specimens were not painted. A E. Sgambiterra et alii, Frattura ed Integrità Strutturale, 30 (2014) 167-173; DOI: 10.3221/IGF-ESIS.30.22 169 Specimens were fatigued by using a servo hydraulic testing machine (Instron 8500), at rate of 5 Hz, to initiate and grow a crack from the EDM notch tip at constant load amplitude and a load ratio, R=min/max, of roughly zero, in order to get a length ratio, a/W, close to 0.3 (see Fig. 2). Pre-cracked specimens were subsequently subjected to complete tensile loading-unloading cycles by using an electro-mechanical testing machine (MTS Criterion s42) equipped with 5 kN load cell; different tests were carried out at different operating temperature, in the range 298-338 K and with a maximum load P=300 N, corresponding to a maximum stress max=Pmax/Wt=60 MPa. Figure 1: Loading-unloading isothermal stress-strain cycle (298 K) Figure 2: Single Edge crack specimen A digital camera (Sony XCD- X910 model) with a resolution of 1280 by 960 pixels of 4.65 m was used to capture images throughout measurement tests. The focus of the images was performed using a Linos Photonics microscope objective with a 4x magnification and a numerical aperture of 0.1, which ensures, in conditions of correct illumination, a resolution of approximately 2.5 m, and therefore lower than the pixel size of the camera. DIC was performed on images from each test carried at different operating temperature using a special tool available in the Matlab® software platform. The first image in the measurement cycle (at minimum load) was used as the reference image, and terms up to first order displacement gradients were used in all correlations. DIC was used to obtain full-field displacements for each image throughout the cycle and the correlations were performed using a subset size with a radius of 42 pixels and a spacing of 2 pixels between subset centers. RESULTS AND DISCUSSIONS yclic tensile tests were performed on a single edge crack specimen, at different operating thermal conditions in order to understand the effect of the temperature on the stress induced transformation mechanisms at the crack tip. Three different values of temperature, in the range of 298-338 K, were adopted during the tests. To determine the effective stress intensity factor the specimen actually experiences, which could be different than the theoretical one determined by the Linear Elastic Fracture Mechanic (LEFM) theory due to the marked non-linear behavior of NiTi alloys, a least squares regression was performed on the DIC measured displacements. The displacement field, at the crack tip, proposed by Williams [31],was used to fit the experimental data obtained from the tests. As the displacement field is not unique unless rigid body motion is also specified, two in-plane rigid motion terms, rotation, A, and rigid translation perpendicular to the crack line, B, have to be specified. Furthermore, due to the low magnitude of the pictures, if the T-stress term is included, in order to improve the fitting of the analytical and experimental results, the asymptotic crack tip displacement equations become:      1 1sin 1 cos sin cos 2 2 2 2 2 1 IK r vv k Tr Ar B v                                 (1) where KI is the effective stress intensity factor, r is the distance from the crack tip, C E. Sgambiterra et alii, Frattura ed Integrità Strutturale, 30 (2014) 167-173; DOI: 10.3221/IGF-ESIS.30.22 170  is the angle from the crack line ahead of the tip,  is the Poisson’s ratio,  is the shear modulus given by:  2 1 E v    (2) and k is given by: 3 1 v k v    (3) The vertical displacement filed, recorded at the maximum applied load (P=300 N) for a specimen with a/W=0.32, is given in Fig. 3. The experimentally obtained and regressed displacements were plotted together to demonstrate the accuracy of the regression technique. In particular, blue contours represent the experimentally found displacements and the red contours represent the regressed displacement contours. As observed, the experimental and regressed displacement contours show good agreement. For a given value of the Young’s modulus, which was assumed to be constant and equal to that of austenitic untransformed structure (EA=68GPa), the fitting procedure allows a direct estimation of the mode I stress-intensity factor (KI). However, it is important to underline that, due to the transformation mechanisms occurring in the crack tip region, the effective elastic properties changes, as consequence of the generation of the new phase and its variants reorientations, therefore further estimations should be carried out in this way to well characterize the real stress state involving the crack tip in a pseudoelastic NiTi alloy. Figure 3: Comparison between the experimental and regresses vertical displacements (v) near the crack tip. Fig. 4, shows the evolution of the calculated KI as a function of the applied load during a complete loading-unloading cycle for all the investigated operating temperatures. As shown, all the curves exhibit two different slopes: a lower slope for load values lower than 50 N and an higher one for load values between 50 N and 300 N. This behavior can be attributed to the unique non-linear stress-strain behavior of NiTi alloys as well as to the crack opening and closure mechanisms. Furthermore, it is possible to observe that the higher the operating temperature the lower the recorded stress intensity factor. In Fig. 5, the stress intensity factor, recorded at the maximum applied load (P=300 N), is plotted in function of the temperature. Results revealed that KI tends to decrease by increasing the temperature until the material exhibits pseudoelastic properties (T<320 K). For higher values of temperature (T>320 K), transformation mechanisms tends to E. Sgambiterra et alii, Frattura ed Integrità Strutturale, 30 (2014) 167-173; DOI: 10.3221/IGF-ESIS.30.22 171 vanish because the energy required to create slip plane is lower than the one needed to induce phase transformation, therefore the stress intensity factor (KI) tends to stabilize. Figure 4: Evolution of the calculated Stress Intensity Factor (SIF) in a cyclic test for different operating temperatures. Figure 5: Evolution of the Stress Intensity Factor (SIF), recorded at the maximum applied load (P=300 N), as a function of the operating temperature. CONCLUSIONS Igital image correlation technique was adopted in this investigation to study the phase transition mechanisms at the crack tip in NiTi alloys. To this aim Single Edge Crack (SEC) specimens, obtained from a commercial NiTi sheets with pseudoelastic behavior, were analyzed. DIC method was used to analyze the displacement field in the crack tip region and a fitting procedures, based on the William’s series expansion, were properly implemented to estimate the mode I Stress-Intensity Factor (SIF). The effect of the operating temperature on the SIF was investigated and results showed that the stress intensity factor decrease by increasing the temperature and it tends to stabilize when the pseudoelastic properties of the alloy vanish. D E. Sgambiterra et alii, Frattura ed Integrità Strutturale, 30 (2014) 167-173; DOI: 10.3221/IGF-ESIS.30.22 172 REFERENCES [1] Otsuka, K., Ren, X., Physical metallurgy of Ti-Ni based shape memory alloys, Progress in Material Science, 50 (2005) 511. [2] Otsuka, K., Wayman, C.M., editors. Shape memory materials, Cambridge: Cambridge University Press, (1998). [3] Duerig, T.W., Melton, K.N., Stokel, D., Wayman, C.M., editors, Engineering aspects of shape memory alloys, London: Butterworth, (1990). [4] Schmahl, W.W., Khalil-Allafi, J., Hasse, B.,Wagner, M., Heckmann, A., Somsen, C., Investigation of the phase evolution in a super-elastic NiTi shape memory alloy (50.7 at.%Ni) under extensional load with synchrotron radiation, Materials Science and Engineering A, (2004), 378:81–5. [5] Chen, J.H., Sun, W., Wang, G.Z., Investigation on the fracture behavior of shape memory alloy NiTi, Metallurgical and Materials Transactions A, (2005) 36:941–55. [6] Eggeler, G., Hornbogen, E., Yawny, A., Heckmann, A., Wagner, M., Structural and functional fatigue of NiTi shape memory alloys, Materials Science and Engineering A, (2004) 378:24–33. [7] Maletta, C., Sgambitterra, E., Furgiuele, F., Casati, R., Tuissi, A., Fatigue of pseudoelastic NiTi within the stress- induced transformation regime: A modified Coffin-Manson approach, Smart Materials and Structures (21), (2012) 1– 9. Art. no. 112001. [8] Maletta, C., Sgambitterra, E., Furgiuele, F., Casati, R., Tuissi, A., Fatigue properties of a pseudoelastic NiTi alloy: Strain ratcheting and hysteresis under cyclic tensile loading, International Journal of Fatigue (66) 2014, 78-85. [9] Melton, K.N., Mercier, O., Fatigue of NiTi thermoelastic martensites, Acta Metallurgica, (1978), 27:137–44. [10] Sawaguchi, T., Kaustrater, G., Yawny, A., Wagner, M., Eggeler, G., Crack Initiation and Propagation in 50.9 At. pct Ni-Ti Pseudoelastic Shape-Memory Wires in Bending-Rotation Fatigue, Metallurgical and Materials Transactions A (2003), 34:2847–60. [11] Stankiewicz, J.M., Robertson, S., Ritchie, R.O., Fatigue-crack growth properties of thin-walled superelastic austenitic Nitinol tube for endovascular stents, Journal of Biomedical Materials Research Part A (2006), 81:685–91. [12] Paiva, A., Savi, M.A., An overview of constitutive models for shape memory alloys, Mathematical Problems in Engineering, art. no. 56876 (2006). [13] Robertson, S.W., Mehta, A., Pelton, A.R., Ritchie, R.O., Evolution of crack-tip transformation zones in superelastic nitinol subjected to in situ fatigue: a fracture mechanics and synchrotron X-ray microdiffraction analysis, Acta Materialia, 55(18) (2007) 6198–6207. [14] Daymond, M.R., Young, M.L., Almer, J.D., Dunand, D.C., Strain and Texture Increment during Mechanical Loading of a Crack Tip in Martensitic Shape-memory NiTi, Acta Materialia, 55 (11) (2007) 3929-3942. [15] Gollerthan, S., Young, M.L., Baruj, A., Frenzel, J., Schmahl, W., Eggeler, G., Fracture Mechanics and Microstructure in NiTi Shape Memory Alloys Acta Materialia, 57 (4) (2009) 1015-1025. [16] Gollerthan, S., Young, M.L., Neuking, K., Ramamurty, U., Eggeler, G., Direct physical evidence for the back- transformation of stress-induced martensite in the vicinity of cracks in pseudoelastic NiTi shape memory alloys Acta Materialia, 57 (2009) 5892–5897. [17] Daly, S., Miller, A., Ravichandran, G., Bhattacharya, K., An experimental investigation of crack initiation in thin sheets of Nitinol, Acta Materialia, 55 (2007) 6322-6330. [18] Gollerthan, S., Herberg, D., Baruj, A., Eggeler, G., Compact tension testing of martensitic/pseudoplastic NiTi shape memory alloys Materials Science and Engineering A, 481–482 (2008) 156–159. [19] Bewerse, C., Gall, K.R., McFarland, G.J., Zhu, P., Brinson, L.C., Local and global strains and strain ratios in shape memory alloys using digital image correlation, Materials Science and Engineering A, 568 (2013) 134–142. [20] Wang, X.M., Wang, Y.F., Baruj, A., Eggeler, G., Yue, Z.F., , On the Formation of Martensite in Front of Cracks in Pseudoelastic Shape Memory Alloys Materials Science and Engineering A, 394(1–2) (2005) 393–398. [21] Maletta, C., Falvo, A., Furgiuele, F., Leonardi, A., Stress-Induced Martensitic Transformation in the Crack Tip Region of a NiTi Alloy, Journal of Materials Engineering and Performances, 18(5–6) (2009) 679–685. [22] Freed, Y., Banks-Sills, L., Crack Growth Resistance of Shape Memory Alloys by Means of a Cohesive Zone Model, Journal of the Mechanics and Physics and Solids, 55 (2001) 2157–2180. [23] Yi, S., Gao, S., Fracture Toughening Mechanism of Shape Memory Alloys Due to Martensite Transformation, International Journal of Solids and Structures, 37(38) (2000) 5315–5327. [24] Xiong, F., Liu, Y., Effect of Stress-Induced Martensitic Transformation on the Crack Tip Stress-Intensity Factor in Ni-Mn-Ga Shape Memory Alloy, Acta Materialia, 55(16) (2007) 5621–5629. E. Sgambiterra et alii, Frattura ed Integrità Strutturale, 30 (2014) 167-173; DOI: 10.3221/IGF-ESIS.30.22 173 [25] Maletta, C., Furgiuele, F., Analytical modeling of stress-induced martensitic transformation in the crack tip region of nickel–titanium alloys, Acta Materialia, 58 (2010) 92-101. [26] Maletta, C., Furgiuele, F., Fracture control parameters for NiTi based shape memory alloys, International Journal of Solids and Structures, 48 (11) (2011) 1658-1664. [27] Maletta, C., Young, M.L., Stress-induced martensite in front of crack tips in NiTi shape memory alloys: modeling versus experiments, Journal of Materials Engineering and Performances, 20 (2011) 597–604. [28] Baxevanis, T., Lagoudas, D., A mode I fracture analysis of a center-cracked in NiTi shape memory alloy panel under plane stress International Journal of Fracture, 175 (2012) 151–166. [29] Maletta, C., A novel fracture mechanics approach for shape memory alloys with trilinear stress–strain behaviour, International Journal of Fracture, 177 (2012) 39–51. [30] Maletta, C., Sgambitterra, E., Furgiuele, F., Crack tip stress distribution and stress intensity factor in shape memory alloys, Fatigue and Fracture of Engineering Materials and Strictires, 36 (9) (2013) 903-912. [31] Eftis, J., Subramonian, N., Liebowitz, H., Crack Border Stress and Displacement Equations Revisited, Engineering Fracture Mechanics, 9 (1977) 189-210.